请大虾们帮我看看哪里出错了怎么画不出来图呀？MATLAB里面正负号格式是不是出错了？

Question

syms z ;%v=sym('v');
k=1.381e-23;
NA=6.022e23;%NA=6.022e23; %Mmol=28;mu=4.65e-26 %以氧气为例 %Mmol =mu/NA;
z = -0.01:0.001:0.01;
T=300;
Mmol=30;
m=(Mmol*10^(-3))/NA;
c1=m/(2*k*T);
boln=5.8*10^(-5)*1.6*10^(-19);
delB=400;
d=0.07;
D=0.42;
c2=boln*delB*d*D/m;
f='±'(c1*c2).^2*exp(('±'c1*c2)./z).*(abs(z)).^3
plot(z,f,'b:*')
f=('±'c1*c2).^2*exp(('±'c1*c2)./z).*(abs(z)).^3

Answer 1

syms z ;
k = 1.381e-23;
NA = 6.022e23;
z  =  -0.01:0.001:0.01;
T = 300;
Mmol = 30;
m = (Mmol*10^(-3))/NA;
c1 = m/(2*k*T);
boln = 5.8*10^(-5)*1.6*10^(-19);
delB = 400;
d = 0.07;
D = 0.42;
c2 = boln*delB*d*D/m;
f(1,:) = (c1*c2).^2*exp((c1*c2)./(z+eps)).*(abs(z)).^3;
f(2,:) = -(c1*c2).^2*exp((-c1*c2)./(z+eps)).*(abs(z)).^3;
plot(z,f,'b:*');

Answer 2

你那程序难以理解，一塌糊涂。取正负按照你的意思就类似按照y的取值的一个分段函数。y轴原点距离实际是x数据，x轴实际是y数据，在meshgrid时就应该改过来，这是原始物理模型的设置问题，我也懒得在前面改了，在最后用mesh(Y,X,Z)代替mesh(X,Y,Z)吧.

clear all
L=160;
x=linspace(0,150,50);
y=linspace(-80,80,50);

[X,Y]=meshgrid(x,y);
Z=((0.2*exp(-0.0221*X)+0.1).*(exp(-0.15*(L/2+Y))+1)).*(Y<0)+((0.2*exp(-0.0221*X)+0.1).*(exp(-0.15*(L/2-Y))+1)).*(Y>0);

mesh(Y,X,Z);
xlabel('采空区距x轴原点距离/m')
ylabel('采空区距y轴原点距离/m')
zlabel('孔隙率\phi / %')
axis ij本回答被网友采纳