PML性能的测试,安照Taflove书上的方法做,似乎很不理想

关于PML性能的测试,安照Taflove书上的方法做,似乎很不理想啊.
测试区域100x100x50,参考区域150x150x150,时间步100 PML厚度6
源
Is=round(nx/2);Js=round(ny/2);Ks=round(nz/2);
 if iteration <= 40
    stimulus = (1/32)*(10-15*cos(pi*iteration/20) + ...
                6*cos(2*pi*iteration/20)-cos(3*pi*iteration/20));
 else 
    stimulus=0;
 end  
 Ez(Is,Js,Ks)=stimulus;
全局误差
 e_global=sum(sum(sum(abs(Ez_T-Ez_B).^2)))
ex_global= 0.0316
ey_global=0.007
Ez_global=0.0319
大家实现的效果如何?>?
附图(mesh的那个)是x=50面上的ex,ey,ez,

申明：网友回复良莠不齐，仅供参考。如需专业解答，请学习本站推出的微波射频专业培训课程。