matlab 错误使用 mesh (line 71)X、Y、Z 和 C 不能是复数。

Question

R1=0.5;
R2=0.6;
beta=0.21;
F=5;
a=R1*sin((pi/2)-beta);
sigmaMax2=(3*F)/(2*pi*a^2);
n1=0.5-0.24*exp((-15.08)*(1-(a/R2)));

theta1=0:pi/90:pi/2;
r4=R1*sin(theta1);
sigma3=subs(sigmaMax2*(1-((r4.^2)/(a.^2))).^n1);

figure(1)
plot(theta1,sigma3);
xlabel('theta1');
ylabel('sigma3');
figure(2)
axis square
[theta1,sigma3,z]=cylinder(sigma3,3000);
mesh(theta1,sigma3,z);
xlabel('theta');
ylabel('sigma3');
zlabel('z');
急求大神指教，MATLAB初学好多不懂,

Answer 1

修改源程序的第11行，求sigma3的模即可（abs函数）。

sigma3=abs(subs(sigmaMax2*(1-((r4.^2)/(a.^2))).^n1));

运行结果如下：