问一个二阶Mur边界条件的问题

12-08

最近用matlab编FDTD，用了《电磁波时域有限差分方法》中给出的
二阶Mur边界条件，但是运行150步之后就出问题了，在边界数值出
现强烈反射。
有人参考过这本书给出的方法吗？出现过这种问题吗？

我是完全按照书里面给的方法设置的，应该不会有问题吧。
另外，书里只推导了一个边界面，然后给出一个通用公式，说其他面也都适合，
我还没有检验，应该不是这个问题吧。
还有，关于棱边的处理，用什么方法比较好呢？