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

模拟一维正弦平面波在自由空间的传播，只有Ei和Hi两个分量，如果入射波振幅为1，那么空间各点的振幅是不是应该为1呢，为什么我的结果是0.5左右，是不是我的入射波加错了？程序如下，matlab里运行，请大家指教，谢谢! O$ p  
%模拟一维平面波在自由空间的传播 cOhx  
%设置参数部分 T;?=,'u  
clear all; c~,OU7[  
c=3.0e8;                          %光速 %8U/!(.g  
mu0=4*pi*1.0e-7;                  %真空磁导常数 abs\Ku9  
e0=8.85e-12;                      %真空介电常数 '+%<\.$  
lambda=532e-9;                    %波长 g3fxf(iY(  
f=c/lambda;                       %频率 cD&QN9  
dw=10;                           %数值波长 'r/+za:2  
ds=lambda/dw;                     %网格空间步长 }M
W7,F
 
dt=ds/(2*c);                      %时间步长 ?o0ro?9j  
steps_t=200*dw;                    %模拟时间步的次数 8
:D|[u;iG  
source_location=2*dw;             %光源的坐标位置 y~16o   
t0=20;                            %高斯脉冲参数设置 {nr}C4]o  
spread=6;                         %高斯脉冲参数设置 U"SH fI:  
PML_layer=9;                      %PML层数 U;PGBoe  
n=4;                              %PML层电导率阶数 >VnBWa<j3  
grids=10*dw;                      %一维网格数目 |"gg2p  
Ei(grids+1)=0;                    %入射平面波电场分量 *!mT#Vm^  
Hi(grids)=0;                      %入射平面波磁场分量 .N:& {$o:  
sigma_m=1/(30*pi*ds);             %PML层电导率最大值 )iKV"jsC  
num=PML_layer*2:-1:1; cu~dbv6H  
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)); A(j9T,!  
sigma_x(PML_layer*2+1)=sigma_m*(ds/2)^(n+1)/((ds*PML_layer)^n*ds*(n+1)); )Id.yv}_  
rho_x=sigma_x*mu0/e0; Njp?/r  
E_M=length(Ei); yX`5x^wVw  
H_M=length(Hi); H{tG:KH  
for T=1:steps_t Y! 8 I  
    pulse=sin(2*pi*f*T*dt);
\Uz7ar#,  
    %pulse=exp(-0.5*((t0-T)/spread)^2); htgtgW9 ^P  
%%%%%%%%%%%%%%%%%%%%%%%%%%%  Ei  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% &>jSuvVT  
    for i=1:PML_layer+1                                                                                            %刷新上侧PML层Ei (vO\h8  
        if i==1 -a^%9 U  
            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); pgp@Zw)r)k  
        else ^/"[jq3F  
            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)); ^0x0 rY  
        end h vC gd^M  
    end obRYU
|T  
    for i=E_M-PML_layer:E_M                                                                                        %刷新下侧PML层Ei .()|0A B&g  
        if i==E_M `jI$>{oa  
            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)); &1Zq
C;  
        else s|cL mL[  
            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)); ff{L=uj  
        end VLL
CdZ%  
    end WUN|,P`b  
    for i=PML_layer+2:E_M-PML_layer-1                                                %刷新入射波电场分量 w#-J ?/m  
        Ei(i)=Ei(i)+0.5*(Hi(i-1)-Hi(i)); ;$il_xA)\>  
    end ~4T:v_Q7g  
    L lNd97Z  
%%%%%%%%%%%%%%%%%%%%%%%%%%%  Hi  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% d_ [l{  
    Hi(source_location)=pulse; ,\%qERk  
    for i=1:PML_layer                                                                                              %刷新左侧PML层Hi 0Z11V9Jk  
        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)); GC#s;X  
    end fV>12ici  
    for i=H_M-PML_layer+1:H_M                                                                             &nbs .. R&vV!d  
WVa-0;  

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