例6-2.m

来自「这是一个MATLAB7.0基础与提高例题的所有源码」· M 代码 · 共 43 行

%输入参数，N为网格线数目
a=input('a=');
b=input('b=');
c=input('c=');
d=input('d=');
N=input('N=');                               

%建立x,y网格坐标
xgrid=linspace(-abs(a),abs(a),N);
ygrid=linspace(-abs(b),abs(b),N);
%确定N×N个点的x,y网格坐标
[x,y]=meshgrid(xgrid,ygrid);                    
z=c*sqrt(d-y.*y/b/b-x.*x/a/a);   
%u=1表示Z要取正负值              
u=1;    
%取Z的实部Z1                                           
z1=real(z);                                   

%for循环用于取消Z中含虚数的点	
for k=2:N-1
  for j=2:N-1
    if imag(z(k,j))~=0
       z1(k,j) = 0;
   end
   if all(imag(z([k-1,k+1],[j-1,j+1]))) ~=0
      z1(k,j) = NaN;
   end
  end
end
%画空间曲面
surf(x,y,z1),hold
%u=1时加画负半面，并加负坐标轴
if u==1                                        
   z2=-z1;
   surf(x,y,z2);
    axis([-abs(a),abs(a),-abs(b),abs(b),-abs(c),abs(c)]);
end
%加标注
xlabel('x'),ylabel('y'),zlabel('z')          		  

title(['[a,b,c,d,N]=[',int2str(a),',',int2str(b),',',int2str(c),',',int2str(d),',',int2str(N),']'])
hold off