方程的代码如下:r1=36.65;w=16;h=80;n=1;r2=5.15;theta0=atan((r1-r2)/h);phi1=0:pi/360:2*n*pi;x1=(r1-h*tan(theta0)/(2*pi*n)*phi1).*cos(-phi1);y1=(r1-h*tan(theta0)/(2*pi*n)*phi1).*sin(-phi1);z1=h/(2*pi*n)*phi1;plot3(x1,y1,z1,'r','linewidth',3)