本示例说明如何为表面设置动画。具体而言,此示例对球形谐波进行了动画处理。球形谐波是傅立叶级数的球形版本,可用于模拟地球的自由振动。
定义球面网格
在球形网格上定义一组点以计算谐波。
theta = 0:pi/40:pi; % polar angle
phi = 0:pi/20:2*pi; % azimuth angle
[phi,theta] = meshgrid(phi,theta); % define the grid
计算球谐函数
在半径等于5的球体表面上计算度数为6,阶数为1,振幅为0.5的球谐函数。然后,将值转换为笛卡尔坐标。
degree = 6;
order = 1;
amplitude = 0.5;
radius = 5;
Ymn = legendre(degree,cos(theta(:,1)));
Ymn = Ymn(order+1,:)';
yy = Ymn;
for kk = 2: size(theta,1)
yy = [yy Ymn];
end
yy = yy.*cos(order*p