出现的这个地面如何去掉呢?
程序如下: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)

 Réponse acceptée

0 votes

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)

Plus de réponses (0)

Catégories

En savoir plus sur 计算几何学 dans Centre d'aide et File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!