Question

在matlab中如何绘制部分球面网格图

球面方程为：x^2+y^2+z^2=300^2
x∈[-150~300]，y∈[0~300]，z∈[-150~150]
请给出具体程序，若正确追加分值，需要网格图，最好网格密度小一点
如何在这个球面基础上，绘制一些平面与该球面的交割线呢？这些平面经过Z轴，且每个平面夹角为14.3°，绘制15个这样的交割线，请高手指点

Answer 1

你的数据范围大，你要网格密度小，这得看你的机器的本事了。
n=50; %调大这个
[x,y,z]=sphere(n);
x=x*300;
y=y*300;
z=z*300;
z(x<-150 | x>300)=nan;
z(y<0 |y>300)=nan;
z(z<-105|z>150)=nan;
mesh(x,y,z)

追问

你这个程序出错了，怎么回事？

追答

你复制进去，
等一下再改n=50
我这里没有问题。

追问

??? Error using ==> sphere
Too many output arguments.
复制进命令框，回车之后出现以上错误信息，怎么回事？你说的等一下再改50是什么意思？

追答

你的MATLAB版本是多少？
我计算的结果如图。

n是将球多少等分，n增大，就提分辨率，但计算量也增大

追问

我发现在我的matlab7.0中一输入[x,y,z]=sphere(50)就会出错，而输入sphere(50)，就会出现半径为50的半球，不是一个完整的球，而且这个球的参数方程不是我现在设置的，而是我之前设置过的，我感觉自动默认了之前的参数设置，太奇怪了，请问怎么处理？

追答

你看一下
which sphere
返回的是什么，是不是你的同名函数文件。

Answer 2

x=-150:300;y=0:300;
[X,Y]=meshgrid(x,y);
z1=sqrt(90000-X.^2-Y.^2);
z2=-sqrt(90000-X.^2-Y.^2);
for i=1:length(x)*length(y)
if z1(i)>150
z1(i)=NaN;
end
if z2(i)<-150
z2(i)=NaN;
end
end
mesh(x,y,z1)
hold on
mesh(x,y,z2)
hold off

追问

如何在这个球面基础上，绘制一些平面与该球面的交割线呢？这些平面经过Z轴，且每个平面夹角为14.3°，绘制15个这样的交割线，请高手指点