[求助]电磁波在分层均匀磁化等离子层间的反射系数仿真

%epsilon1 为介电常数1； $Jn.rX0}$  
%ne       为电子数密度 OEzSItAI/[  
%qe       为电量 )k[XO  
%me       电子质量 EV* |\ te  
%Wp       等离子体频率 \>EUa}%xn  
%Wg       电子在外磁场中的自旋频率 mr_N
ArF  
%gamma    入射波在等离子体中的传播常数 qJK9
C`T%  
%W        入射波的频率 X|g5tnsj`  
%gamma1   入射波在第一种介质中的传播常数 vO"E4s  
%gamma2   入射波在第二种介质中的传播常数 J8@+)hn  
%R12      反射系数 =>Z4vWX*  
%k        有效碰撞频率 OI:=>Bk  
%B0       磁感应强度 q3-cWfU  
%d        等离子子层厚度 SX =^C  
%c        光在真空中的传播速度 *!Gb_!98  
clear; P@n rcgM.  
clc; -;HZ!Lf  
close all; N!v>2"x8q  
epsilon0=1/(36*pi)*10^(-9); rTJU)4I^h  
qe=1.6*10.^(-19); 'u%_Ab_H  
me=9.3*10^(-31); 'Cz]p~oF  
ne=linspace(10.^17,10.^18,1000); yDe*-N\'W  
B0=0.1; DIvxut  
mu0=4*pi*10^(-7); :x3DuQP  
k=linspace(10.^10,10.^11,1000); dvY3=~'  
x2l}$(7  
C0=3.0*10^8; kDE:KV<"c  
R(1001)=-1; JlE+CAny  
for W=1:70 -cZuP7oA  
d=0.01; /Jc^XWf  
Wp=sqrt(ne*qe^2/(me*epsilon0)); %Wp^2=ne*qe^2/(me*epsilon0) (5y*Btd=  
Wg=qe*B0/me; 8(A k  
epsilon1=1+(Wp.^2./W).*(W-i*k)./(Wg^2-(W-i*k).^2); 9WQ'"wyAQ  
epsilon2=(Wp.^2./W)*Wg./(Wg^2-(W-i*k).^2); %4#
ChlXB  
epsilon3=1-Wp.^2./(W.*(W-i*k)); nHfAx/9!  
gamma=sqrt(-W.^2*mu0*epsilon0.*(epsilon1-epsilon2.^2./epsilon1)); IsFL"Vx  
[5M!'  
gamma0=i*W*sqrt(mu0*epsilon0); {+kWK;1  
R(1)=(gamma0-gamma(1))/(gamma0+gamma(1)); /z)8k4  
Q(1)=R(1); 4IM&#_6  
for i=2:1000 K&=6DvfR  
R(i)=(gamma(i-1)-gamma(i))/(gamma(i-1)+gamma(i)); aK,\e/Oo  
Q(i)=R(i);  %\B?X;(  
end n4"xVDL  
for j=2:1000 #a8i($k{e  
    r(j)=Q(1002-j); 1OqVNp%K  
    r(j)=(r(j)+r(j-1)*exp(-2*gamm .. 4zBcq<R7  
,s@S`KS0  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

搂主，你的10.^11为什么要加点，还有M(:,10)貌似为一个矩阵的一列，但是前面定义的是M(j)=r(j)为向量

好久没有看等离子体的书了，几年前的积累基本快忘完了，看了一下你的代码，觉得两i 种原因，一是参数选择不合理，反射系数本身不会有什么大的差别， D,;6$Pvg^  
2， 仔细检查公式，epsilon1，epsilon2的数量级一致，导致gama的计算结果没有大的差别，最终的计算结果自然没有大的差别。