[已解决] matlab使用mesh函数,出现了底面,如何去掉?

[复制链接]
123MATLAB呀 发表于 2021-11-24 11:28:47
出现的这个地面如何去掉呢?
程序如下:c=3*10^8;
lemdal=1550E-12;
ommiga=lemdal*c;
syms z;
z=90:1:150;
K=exp(-(z-125)/6);
Vin=100*K;
me=9.3E-31;
e=1.6E-19;
chi=45/180*pi;
q0=4.7E3;alpha=1.6E-7;
H=(z-105)./10;
Ne=sqrt(q0/alpha).*exp(1/2.*(1-H-1./cos(chi).*exp(-H)));
epsilong0=8.85E-12;
ommiga_p=sqrt((Ne.'*e^2)./(me*epsilong0));
elta=sqrt(1-ommiga_p.^2./(ommiga.^2+i*ommiga*Vin));
k=imag(elta);
Aref=lemdal*k;
Aabs=sqrt(-1/2.*(1-ommiga_p.^2./(ommiga.^2+Vin.^2))+1/2.*sqrt((1-ommiga_p.^2./(ommiga.^2+Vin.^2)).^2+(Vin./ommiga.*(ommiga_p.^2./(ommiga.^2+Vin.^2))).^2));  

fun=@(z)(Aref+Aabs);
Aatt=integral(fun,90,150,'ArrayValued',true);
vpa(Aatt);  
Ne1=Ne(1,:);
z1=z(1,:);
[z1,Ne1] = meshgrid(90:1:150, 2.5151e+04:2.9550e+04);  
mesh(Ne,z,Aatt)
brighten(0.5)


图片.PNG

最佳答案


WarnerChang 发表于 2021-11-24 11:56:13
c=3*10^8;
lemdal=1550E-12;
ommiga=lemdal*c;
syms z;
z=90:1:150;
K=exp(-(z-125)/6);
Vin=100*K;
me=9.3E-31;
e=1.6E-19;
chi=45/180*pi;
q0=4.7E3;alpha=1.6E-7;
H=(z-105)./10;
Ne=sqrt(q0/alpha).*exp(1/2.*(1-H-1./cos(chi).*exp(-H)));
epsilong0=8.85E-12;
ommiga_p=sqrt((Ne.'*e^2)./(me*epsilong0));
elta=sqrt(1-ommiga_p.^2./(ommiga.^2+i*ommiga*Vin));
k=imag(elta);
Aref=lemdal*k;
Aabs=sqrt(-1/2.*(1-ommiga_p.^2./(ommiga.^2+Vin.^2))+1/2.*sqrt((1-ommiga_p.^2./(ommiga.^2+Vin.^2)).^2+(Vin./ommiga.*(ommiga_p.^2./(ommiga.^2+Vin.^2))).^2));  

fun=@(z)(Aref+Aabs);
Aatt=integral(fun,90,150,'ArrayValued',true);
vpa(Aatt);  

[z1,Ne1] = meshgrid(z,Ne);  
mesh(z1,Ne1,Aatt)
brighten(0.5)
回复此楼

2 条回复


WarnerChang 发表于 2021-11-24 11:56:13
c=3*10^8;
lemdal=1550E-12;
ommiga=lemdal*c;
syms z;
z=90:1:150;
K=exp(-(z-125)/6);
Vin=100*K;
me=9.3E-31;
e=1.6E-19;
chi=45/180*pi;
q0=4.7E3;alpha=1.6E-7;
H=(z-105)./10;
Ne=sqrt(q0/alpha).*exp(1/2.*(1-H-1./cos(chi).*exp(-H)));
epsilong0=8.85E-12;
ommiga_p=sqrt((Ne.'*e^2)./(me*epsilong0));
elta=sqrt(1-ommiga_p.^2./(ommiga.^2+i*ommiga*Vin));
k=imag(elta);
Aref=lemdal*k;
Aabs=sqrt(-1/2.*(1-ommiga_p.^2./(ommiga.^2+Vin.^2))+1/2.*sqrt((1-ommiga_p.^2./(ommiga.^2+Vin.^2)).^2+(Vin./ommiga.*(ommiga_p.^2./(ommiga.^2+Vin.^2))).^2));  

fun=@(z)(Aref+Aabs);
Aatt=integral(fun,90,150,'ArrayValued',true);
vpa(Aatt);  

[z1,Ne1] = meshgrid(z,Ne);  
mesh(z1,Ne1,Aatt)
brighten(0.5)
回复此楼

123MATLAB呀 发表于 2021-11-24 17:29:15
WarnerChang 发表于 2021-11-24 11:56
c=3*10^8;
lemdal=1550E-12;
ommiga=lemdal*c;

感谢大佬!!!
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表