第二类Wiener-Hopf积分方程的高精度自适应算法_论文.docVIP

学位论文原创性声明 本论文是我个人在导师指导下进行的工作研究及取得的研究成果. 论文中 除了特别加以标注和致谢的地方外, 不包含其他人或其它机构已经发表或撰写 过的研究成果. 对本文的研究做出贡献的个人和集体, 均已在论文中以明确方 式标明. 本人完全意识到本声明的法律责任由本人承担.  论文作者签名：  日期：  年  月  日  学位论文使用授权的说明 本人授权汕头大学保存本学位论文的电子和纸质文档, 允许论文被查阅和 借阅；学校可将本学位论文的全部或部分内容编入有关数据库进行检索, 可以 采用影印、缩印或其它复制手段保存和汇编论文；学校可以向国家有关部门 或机构送交论文并授权其保存、借阅或上网公布本学位论文的全部或部分内 容. 对于保密的论文, 按照保密的有关规定和程序处理. 本论文属于：保密（ ）, 在 年解密后适用本授权声明. 不保密（ ）. （请在以上括号内打“ ”）  论文作者签名：  指导教师签名： 日期： 年 月 日 日期： 年 月 日 汕 头 大 学 硕 士 学 位 论 文  摘  要  本文主要研究核函数 k(t ? s) 为半光滑函数的第二类 Wiener-Hopf 积分 方程  x(t) + ∫  ∞  k(t ? s)x(s)ds = y(t),  0 ≤ t ∞ 0 的数值解法. 针对该类方程, 已经有很多高精度数值解法被提出. 本文主要讨论应用 Nyst¨orm-Clenshaw-Curtis (NCC) 求积方法和 Clen- shaw-Curtis-Rational (CCR) 求积方法离散上述积分方程. 在简单介绍 NCC 求 积方法之后, 我们证明了其权系数矩阵的两个基本性质, 并提出 NCC 的一种 修正方法—核函数的一种新的延拓方法. 然后, 我们结合 CCR 求积方法和修 正的 NCC 求积方法, 提出一种混合的积分方法. 具体如下: 把积分方程写成 如下的积分方程组 ? ? ? ? ? ∫ R 0  k(t ? s)x(s)ds + k(t ? s)x(s)ds + ∫ ∞ R  k(t ? s)x(s)ds = y(t), k(t ? s)x(s)ds = y(t),  0 ≤ t ≤ R, R ≤ t ∞. 0 R  对以上方程组的各个积分应用合适的求积方法进行离散, 得到线性方程组并对 之求解, 从而得到精度比较高的数值解. 我们还考虑一些误差估计方法并由此推导 Wiener-Hopf 积分方程的自适 应求解方法. 最后, 我们给出数值例子检验本文提出的方法的有效性. 通过对 数值结果的分析比较, 可以看出本文提出的算法的优势. 关键词: 第二类 Wiener-Hopf 积分方程; 半光滑核函数; CCR 求积方法; 修正 NCC 求积方法; 自适应解法. I 汕 头 大 学 硕 士 学 位 论 文 Abstract This thesis considers numerical methods for the solution of Wiener-Hopf equations of the second kind  x(t) + ∫ ∞  k(t ? s)x(s)ds = y(t),  0 ≤ t ∞, 0 where k(·) are semi-smooth functions. Many numerical methods have been proposed and analyzed for such kind of integral equations. This thesis mainly discusses using the Nyst¨orm-Clenshaw-Curtis (NCC) quadrature and Clenshaw-Curtis-Rational (CCR) quadrature to discretize the above integral equation. After brie?y introducing the NCC quadrature, we prove two basic properties of the weighted matrices of the quadrature, and propose a modi?ed NCC based on a new extension method for kernel functions. Then, we propose a combined quadrature for the integral equation by applying the CCR quadrature and