[求助]二维FDTD问题，小孔衍射设置

我是初学者，想用二维FDTD做小孔衍射的模拟，设置的小孔如下： ~rO&Y{aG#  
-----------                 ---------- p*jU)@a0  
   10            20             10 :_i1gY)  
2m*ugBO;  
但是画出来的磁场的图（matlab）却在我设置的孔下面多出来一个单元格，不知道为什么！ )|v du  
我的程序就是在经典的Hagness的上面改的，只是把小孔重新设置了一下，就把小孔当做是金属吸收介质设置的： !bLCha\  
for 21：25 <<S4l~"o  
    for j=1:bankuan vcy1itY  
         caex(ii,j)=ca(2); Tx&H1  
         cbex(ii,j)=cb(2); ESoqmCJjb:  
         caey(ii,j)=ca(2); !-T#dU  
         cbey(ii,j)=cb(2); bE% Hm!  
    end [V_
mF  
    for j=(je-bankuan+1):je fhZwYx&t  
         caex(ii,j)=ca(2); Bn8&~  
         cbex(ii,j)=cb(2); =>PX~/o  
         caey(ii,j)=ca(2); r)>'cjx/  
         cbey(ii,j)=cb(2); p&Nav,9x  
    end 6[XaIco=C  
end< .. {BM:c$3@j  
=.S2gO >  

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

哦，有几个地方忘记了说明，bankuan就是板宽，这里是10，上面孔的大小是30，写错了 Qw%0<~<  
这个论坛不知道怎么编辑自己的帖子，所以回一个在这里！！！

1：2表示的是1和2两个格点，但是你用pcolor的时候它只显示出一个格子。 Xw^:<Nx:  
你可以试一试d=[1 2 3;4 5 6;7 8 9];pcolor(d)；3*3的矩阵显示的只有四个格子 Yu
^}  
你的程序怎么样？全做好了吗？我现在也在做，能不能传一个给我看看！ 1^;&?E  
QQ号，24765428

这个是什么啊。？ RMrt4:-DI