[十五度]

目前正在学习偏移处理算法，并自己用程序实现。最简单的15°时域差分法、相移法、SSF法。都没有问题。

但是在频率空间域有限差分法上面却一直过不去。这也直接导致了FFD也没法继续。

一个最大的问题就在于求解复系数的三对角方程组上。

翻遍很多资料，对于求解这个方程组，其说法都是大同小异的。

不外乎就是追赶法。。。当然，也都是以实系数为例的，并假设满足求解条件。

要理解实现这个追赶法，并不困难。但现在的关键就在于边界条件不知道应该怎么去调节。

我到目前为止没有找到任何一本书提到过这个问题。

也就是说怎么样才能够避免实际中不满足求解三对角线方程组的假设条件？

克莱伯特的书上提到过所谓的elift和eright。这个东西在他的书上，只是提到必须要给，但是却没说怎么给。着实令我十分费解。

因为在实际做得过程中，很可能求解三对角线方程组的条件是不满足的。

到底应该怎么去调节这个边界条件呢？用什么理论方法？

我看了Su的代码（sumigfd.c），里面对于边界的处理十分玄乎，完全搞不清楚为什么要这样。

还请各位牛人帮帮小弟，到底应该如何理解、解决这个问题呢？

[十五度]

没人回答呀，
幸好自己找到了，原来用的是：Clayton and Engquist 1980年提出的频率空间域波动方程边界吸收条件。
真该死克莱伯特的书不讲清楚。

[ytlee]

引用框（十五度 ＠ 03-09-2011 - 08:07）

没人回答呀，
幸好自己找到了，原来用的是：Clayton and Engquist 1980年提出的频率空间域波动方程边界吸收条件。
真该死克莱伯特的书不讲清楚。

Good work!!!!!

Way to go. Keep it up, you're on the right path becoming a great scientist.http://en.wikipedia.org/wiki/Scientist