[求助]平面波程序调试遇到问题了～～急，急，急需各位高手师兄指点

  我根据allen taflove的书写的三维平面波的程序，可是出来的图结果明显不对啊，我检查了程序，感觉连接边界设置的没问题啊，是不是网格不对应啊大家看看下面的图给点意见，急啊！！！！！！！！！！！11 Q32GI,M%B  
[attachment=20319] ka<rlh<h  
  [attachment=20320] Bw4PxJs-  
  [attachment=20321] _))--+cL  
  [attachment=20 .. H1e^/JD)  
`_E@cZ4  

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

[attachment=20324] b&j}f  
根据上图的公式我编写的程序为： W1T% Q88  
ex(ie_jmin:ie_jmax-1,je_jmin,ke_jmin:ke_jmax)=ex(ie_jmin:ie_jmax-1,je_jmin,ke_jmin:ke_jmax)-cb(1)*Hinczd_Left(ie_jmin:ie_jmax-1,ke_jmin:ke_jmax); 2`;&Uwt  
是不是我的网格不对应啊

不知道你图什么意思，每组的三个图上的标注都是一模一样的。

哦，我这三个图是这么来的 0'Si ^>bW  
  timestep=int2str(n); ~_-]> SI  
    figure(1) Bb:C^CHIQm  
    tview(:,:)=squeeze(ez(is,jh_bc:upml+je,kh_bc:upml+ke)); F~8'3!<9  
    uview(:,:)=squeeze(ey(is,jh_bc:upml+je,kh_bc:upml+ke)); Z;%uDlcXI  
    wview(:,:)=squeeze(ex(is,jh_bc:upml+je,kh_bc:upml+ke)); -?-XO<I  
     +4.s4&f)  
    subplot(3,1,1), imagesc(tview'); $'%GB $.  
    shading flat; 4WJY+)  
    caxis([-1.1786e+003  1.1786e+003]); >UMxlvTg&  
    colorbar; "bIb?e2h9G  
    axis image; axis xy; Bz<hP*.O  
    title([' time step = ',timestep]); +?Ii=*7n  
    xlabel('i coordinate'); QJ3#~GYNr  
    ylabel('j coordinate'); Rh ^(91d  
     Oj;*Gi9E  
    subplot(3,1,2), imagesc(uview'); +bS\iw+  
    shading flat; $uZmIu9Bi+  
    caxis([-1.0501e+004  1.0501e+004]); MzD1sWmK  
    colorbar; .P$m?p#  
    axis image; axis xy; '< U&8?S  
    title([' time step = ',timestep]); )"S%'myj  
    xlabel('i coordinate'); 3l~+VBR_  
    ylabel('j coordinate'); 7ocUFY0"  
     N#K)Z5J)b  
    subplot(3,1,3), imagesc(wview'); 4LKpEl.=  
    shading flat; u<Kowt<ci  
    caxis([-9.7198e+003  9.7198e+003]); (T_-`N|  
    colorbar; >UvP/rp  
    axis image; axis xy; Qds<j{2  
    title([' time step = ',timestep]); bGCC?}\  
    xlabel('i coordinate'); 3'Y-~^ml|  
    ylabel('j coordinate'); H-PVV&r  
    xgcxA:  
分别表示固定x时yz面上ez,ey,ex的变化，只是程序刚开始调试我就没有清楚的标注

麻烦师兄给看看吧，要中期检查了，程序还有问题，再调不通，就得延期毕业了～～

你总场里面放了金属没？

没有啊，总场里面是空的啊，我想先把平面波的程序调通了，再加金属，金属那块调的差不多了

里面是真空，那你平面波加的肯定不对了。 5OzEY7K)  
你看出来的图形几乎还是点源一样扩散出去的。 ;j7G$s9  
连接边界是加在20，80的地方么？

我的连接边界是这么设置的 n\3#69VY  
%X方向连接边界的设置 J=t}9.H~=  
%ie_j=60;%X方向的网格数 -0x Q'1I  
ie_jmin=30; -Ep#q&\  
ie_jmax=90; (:HT|gKoE  
%Y方向连接边界的位置 0%dOi ko  
%je_j=60;%Y方向的网格数 `;J`O02  
je_jmin=30; p@nj6N.--  
je_jmax=90  ; 8B:y46
  
%Z方向连接边界的位置 &9fQW?Czs  
%ke_j=60;%Z方向的网格数 T7^;!;i`X  
ke_jmin=30; @uldD"MJ<]  
ke_jmax=90; N+@ Ff3M  
因为图只是显示真空区域的范围，PML的厚度没有显示，PML厚度为10个网格，所以图是从20－80，而我设置的是30－90 1l"2 ~k  
yCvtglAJ4  
后面运行的时候确实发现不对，我是对照allen taflove的书上写的程序，目前怀疑网格设置有问题，不过不知道具体错在什么地方，该怎么改，很迷茫

这个和吸收边界也没关系的，按道理来说，模型里面都是真空，那么电磁场应该不会传播到连接边界以外的地方。 HuevDy4  
-'&/7e6>y  
对这种模型是有解析解的，每一步的电磁场的值都可以预先知道的，应该不难调试啊，就是三维的运行时间长点而已。 \vA*dQ-  
N6Mr#A-{  
比如平面波(XZ)向Y方向传播，正入射，那么在XZ平面的值应该都是一样的，对吧？ I^
y<W%Et  
yKc-:IBb{u  
你调试一下，为什么值不一样？怎么计算出来不一样的？ 什么地方导致不一样的？ 特别是角点的加加减减符号，是不是有一个弄反了？