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

我是初学者，想用二维FDTD做小孔衍射的模拟，设置的小孔如下： t9`.bx8  
-----------                 ---------- ? (Oy\  
   10            20             10 AT
3cc  
6<SAa#@ey  
但是画出来的磁场的图（matlab）却在我设置的孔下面多出来一个单元格，不知道为什么！ J1|\Q:-7p  
我的程序就是在经典的Hagness的上面改的，只是把小孔重新设置了一下，就把小孔当做是金属吸收介质设置的： 7kLz[N6Ll  
for 21：25 \ a<h/4#|  
    for j=1:bankuan ,kGc]{'W  
         caex(ii,j)=ca(2); <z&/L/bl"  
         cbex(ii,j)=cb(2); @})|Z}~  
         caey(ii,j)=ca(2); 6I4\q.^qw  
         cbey(ii,j)=cb(2); J?1 uKR  
    end <?6|.\&  
    for j=(je-bankuan+1):je L|+~"'l  
         caex(ii,j)=ca(2); GW@;}m(  
         cbex(ii,j)=cb(2); 6 7.+ .2  
         caey(ii,j)=ca(2); wE>\7a*P%  
         cbey(ii,j)=cb(2); dr}`H,X"3  
    end [A~xy'T  
end< .. O, wJR  
~"H,/m%2o  

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

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

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

这个是什么啊。？ p?%y82E