0

我通过使用参数方程进行了尝试，但无法生成它。如何改进此代码，以便为不同半径提供半球坐标？

 t=linspace(0,2*pi,50);
    i=1;
    for r = 1:1:50
         x=r*cos(t);
         y=r*sin(t);
         z= r*ones(size(x));
         plot3(x,y,z)
        hold on 
         x1(i,:)=x;
          y1(i,:)=y;
          z1(i,:)=z;
          i=i+1;
    end
    hold off
    j=1;
    
    for phi =linspace(0,pi/2,50)
        xc=x1(j,:).*sin(phi);
        yc=y1(j,:).*sin(phi);
        zc=z1(j,:).*cos(phi);
        zi(j,:)=zc;
        xi(j,:)=xc;
         yi(j,:)=yc;
        plot3(xc,yc,zc)
        hold on
         xlabel('x');
          ylabel('y');
           zlabel('z');
        j=j+1;
     
    end
4

0 回答 0