[已解决]Mur边界条件角点处理求助

边界条件用的是Mur一阶,但角点处理效果不太理想,大家帮忙看看是不是公式有问题,万分感谢 %x;x_  
% 4.角点处理 6SJ"Tni8  
h2=(c0*dt-sqrt(2)*ds)/(c0*dt+sqrt(2)*ds); P=N$qz
$U  
Hz(1,1)=Hz(2,2)+h2*(Hz(2,2)-Hz(1,1));%左下角点 =+5,B\~q@C  
Hz(IE,1)=Hz(IE-1,2)+h2*(Hz(IE-1,2)-Hz(IE,1));%右下角点 Mj>}zbpk/  
Hz(1,JE)=Hz(2,JE-1)+h2*(Hz(2,JE-1)-Hz(1,JE));%左上角点 U8#xgz@  
Hz(IE,JE)=Hz(IE-1,JE .. P47V:E%  
A% Q!^d  

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

好想没有什么问题 qs\2Z@;  
可能是不是你重复计算了角点的吸收？

我查了下资料,原话如下,可能这是算法本身的不足. IcP)FB4  
"这就是Mur所建议的具有二阶近似的，适用于二维问题的近似吸收边界条件。他在FDTD中有广泛地应用。Mur吸收边界条件具有实施方便简单、吸收边界条件效果好的特点，然而在使用中注意到，一阶近Yee网格划分，在角区域存在较大误差，而二阶近似尽管就算精度较高，但编程复杂，且对三维情况还可能出现结果发散的现象"

应该没有问题吧？！我曾经编过类似的程序

我想问一下 上面得代码到底有没有错啊··为何我的运行出来会是这样的结果··是哪得错呢？角点？我也是按照上面的方式编的程序啊····还请大侠给点指点额···· !W?gR.0$=  
运行出得图形如下：  

图片:新建 文本文档 (2).jpg