MATLAB 半球面

这个题我先用的两个参数画的两个图形，但是我真不会用相同的参数把交线表示出来，问了问老师最后是暴力跑一遍点画的

 

r=linspace(0,1,100);
t=0:pi/50:2*pi;
[R T]=meshgrid(r,t);
x=R.*cos(T);
y=R.*sin(T);
z=sqrt((1-x.*x - y.*y)+eps);
mesh(x,y,z);
axis equal
hold on;
z1=R;
x1=cos(T).^2;
y1=sin(T).*cos(T);
mesh(x1,y1,z1);
shading flat
[x y]=meshgrid(-1:.002:1);
z=sqrt((1-x.*x - y.*y)+eps);
for i=1:1:1000
    for j=1:1:1000
        if((1-x(i,j).*x(i,j) - y(i,j).*y(i,j))>=0)
            if(abs(x(i,j)*x(i,j)+y(i,j)*y(i,j)-x(i,j))<=1e-3)
                plot3(x(i,j),y(i,j),sqrt((1-x(i,j).*x(i,j) - y(i,j).*y(i,j))),'b*');
            end
        end
        end
end

修改一下以上做法，实在是有点傻，把T带进去表示出来就完了，我怎么就想不起来呢，一直纠结那个R了

 

r=linspace(0,1,100);  
t=0:pi/50:2*pi;  
[R T]=meshgrid(r,t);  
x=R.*cos(T);  
y=R.*sin(T);  
z=sqrt((1-x.*x - y.*y)+eps);  
mesh(x,y,z);  
axis equal  
hold on;  
z1=R;  
x1=cos(T).^2;  
y1=sin(T).*cos(T);  
mesh(x1,y1,z1);  
shading flat  
x=cos(T).*cos(T);
y=cos(T).*sin(T);
z=(1-x).^0.5;
plot3(x,y,z,'bh','linewidth',5);