[求助]请高手帮忙调试下这个程序

调了好久，也问了很多人，都说是角点反射的问题，这几天书都看了好几遍了，也翻了很多资料，改了很多地方，可是总是不对，求FDTD高手帮忙调试一下。最好告诉小弟问题出在什么地方，谢 .. YKS'#F2  
y:``|*+  

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

俺是没有时间给你调试了 &4iIzw`  
谁给你解决这个问题，奖励一个威望

重赏，重赏啊。 852$Ui|I  
可别浪费机会了。 `Q}.9s_ri  
  

高手在哪里？

Sullivan的书上的程序对角区没有特殊处理，我粗略看了一下，要重新推导他的公式才可以，原来也调过他的程序，后来改自己写了，自己写的改起来更容易些。建议楼主可以试试自己写新的。

我推导他的公式感觉也很有道理啊，而且二维的程序也没有问题啊，为什么三维的就是不正确？

不好意思，昨天太匆忙了，后来仔细看了一下，我觉得像Sullivan那种按六个面赋PML值的话，角区棱区也是被处理到了啊，应该反射不大。，不知道为什么反射大，难道他公式里面用的近似有问题？x(n)的经验值取0.3333，换一个会不会好些？你试过没有？

谢谢楼上提醒，程序我重新改过了。总该觉三维PML似乎没有二维PML好些。 svaclkT=  
下面是我最新做得图，大家看看还有没有问题？