首先看公式
r 的取值为[0,1],且取值的间隔为0.01;
z 的取值为[0,1],且取值的间隔为0.01;
看matlab的m文件代码
r = 0:0.01:1;
z = 0:0.01:1;
meshgrid(r,z);
%以矩阵的形势进行运算出nn的值。
for row = 1:1:101
for col = 1:1:101
nn(row,col) = 0.255e-9*sin(r(row))*cos(z(col))+0.338e-9*sin(2*r(row))*cos(2*z(col))+0.367e-10*sin(3*r(row))*cos(3*z(col));
end
end
surf(r,z,nn);
xlabel('r');% 设置坐标轴显示
ylabel('z');
查看效果