[求助]求助一维平面波的问题

模拟一维正弦平面波在自由空间的传播，只有Ei和Hi两个分量，如果入射波振幅为1，那么空间各点的振幅是不是应该为1呢，为什么我的结果是0.5左右，是不是我的入射波加错了？程序如下，matlab里运行，请大家指教，谢谢! G#>X~qk()  
%模拟一维平面波在自由空间的传播 uOx$@1v,  
%设置参数部分 -TWo-iu^  
clear all; 5+J/Qm8{bb  
c=3.0e8;                          %光速 glpdYg *  
mu0=4*pi*1.0e-7;                  %真空磁导常数 `{":*V   
e0=8.85e-12;                      %真空介电常数 'M{_S  
lambda=532e-9;                    %波长 Ws(>} qjy  
f=c/lambda;                       %频率 h&{pMmS3
,  
dw=10;                           %数值波长 1GxYuTZ{  
ds=lambda/dw;                     %网格空间步长 \I=:,cz*,  
dt=ds/(2*c);                      %时间步长 &ReIe>L  
steps_t=200*dw;                    %模拟时间步的次数 LWE !+(n  
source_location=2*dw;             %光源的坐标位置 -XBNtM_"  
t0=20;                            %高斯脉冲参数设置 zc"eSy< w$  
spread=6;                         %高斯脉冲参数设置 ;L7<mU  
PML_layer=9;                      %PML层数 W0|?R6|  
n=4;                              %PML层电导率阶数 ,LvJ'N  
grids=10*dw;                      %一维网格数目 E^!%m8--  
Ei(grids+1)=0;                    %入射平面波电场分量 o0q{:An_Z  
Hi(grids)=0;                      %入射平面波磁场分量 8b/yT4f  
sigma_m=1/(30*pi*ds);             %PML层电导率最大值 %rsW:nl  
num=PML_layer*2:-1:1; K67x.PZ  
sigma_x=sigma_m*((num*ds/2+ds/2).^(n+1)-(num*ds/2-ds/2).^(n+1))/((ds*PML_layer)^n*ds*(n+1)); Tm `CA0@  
sigma_x(PML_layer*2+1)=sigma_m*(ds/2)^(n+1)/((ds*PML_layer)^n*ds*(n+1)); 03xQ%"TU<  
rho_x=sigma_x*mu0/e0; Kh>^;`h  
E_M=length(Ei);  2r[,w]  
H_M=length(Hi); <8Ad\MU  
for T=1:steps_t bm^ou#]|  
    pulse=sin(2*pi*f*T*dt); 8LL);"$  
    %pulse=exp(-0.5*((t0-T)/spread)^2); ~(c<ioIf  
%%%%%%%%%%%%%%%%%%%%%%%%%%%  Ei  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @ 9uwcM1F  
    for i=1:PML_layer+1                                                                                            %刷新上侧PML层Ei P*}Oi7Z  
        if i==1 lL
?;?V~  
            Ei(i)=exp(-sigma_x(1)*dt/e0)*Ei(i)-(1-exp(-sigma_x(1)*dt/e0))*sqrt(e0/mu0)/(sigma_x(1)*ds)*Hi(i); t|//oEY  
        else &lD4-_2J  
            Ei(i)=exp(-sigma_x(2*i-1)*dt/e0)*Ei(i)-(1-exp(-sigma_x(2*i-1)*dt/e0))*sqrt(e0/mu0)/(sigma_x(2*i-1)*ds)*(Hi(i)-Hi(i-1)); pD{Li\LY  
        end w-1CA{"i7  
    end hOV+}P6  
    for i=E_M-PML_layer:E_M                                                                                        %刷新下侧PML层Ei z')'8155  
        if i==E_M 3k=q>~&@  
            Ei(i)=exp(-sigma_x(1)*dt/e0)*Ei(i)-(1-exp(-sigma_x(1)*dt/e0))*sqrt(e0/mu0)/(sigma_x(1)*ds)*(-Hi(i-1)); e&:fzO<~I  
        else k3Y>QN|q8  
            Ei(i)=exp(-sigma_x(2*(E_M-i)+1)*dt/e0)*Ei(i)-(1-exp(-sigma_x(2*(E_M-i)+1)*dt/e0))*sqrt(e0/mu0)/(sigma_x(2*(E_M-i)+1)*ds)*(Hi(i)-Hi(i-1)); 3wK)vW  
        end |Q{l]D  
    end c;Li~FLR  
    for i=PML_layer+2:E_M-PML_layer-1                                                %刷新入射波电场分量 v
e$P=ZuM  
        Ei(i)=Ei(i)+0.5*(Hi(i-1)-Hi(i)); ? in&/ZrB  
    end =I?p(MqW  
    6>l-jTM  
%%%%%%%%%%%%%%%%%%%%%%%%%%%  Hi  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ]fR 3f  
    Hi(source_location)=pulse; q+N}AKawB  
    for i=1:PML_layer                                                                                              %刷新左侧PML层Hi rh+OgKi  
        Hi(i)=exp(-rho_x(2*i)*dt/mu0)*Hi(i)-(1-exp(-rho_x(2*i)*dt/mu0))/(rho_x(2*i)*ds)*sqrt(mu0/e0)*(Ei(i+1)-Ei(i)); #xO`k1W.  
    end P~~RK&+i  
    for i=H_M-PML_layer+1:H_M                                                                             &nbs .. YpiSH(70`  
!nu#r$K(  

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