这是个椭球体的方程,直接用直角坐标表示的话,难免会出现开平方存在多值的问题,所以一般的做法是用球面坐标表示,然后再转换为直角坐标来绘图。
示例代码:
123456789101112131415161718% 网格数量n = 50;theta = (-n:2:n)/n*pi;phi = (-n:2:n)'/n*pi/2;cosphi = cos(phi); cosphi(1) = 0; cosphi(n+1) = 0;sintheta = sin(theta); sintheta(1) = 0; sintheta(n+1) = 0; % 变换到直角坐标x = sqrt(2)*cosphi*cos(theta);y = sqrt(4)*cosphi*sintheta;z = sqrt(6)*sin(phi)*ones(1,n+1); % 绘图mesh(x,z,y)axis equalxlabel('x')ylabel('y')zlabel('z')
事实上,完全可以更简单一些——直接使用ellipsoid函数:
1ellipsoid(0,0,0,sqrt(2),sqrt(4),sqrt(6),50)
关于函数的参数说明,请自行查阅文档。