第一类强奇异积分方程的快速数值算法.pdfVIP

２０１５ 年６ 月 ： Ｊｕｎ ．２０１５ 第３２ 卷 第２ 期 Ｊｏｕｒｎａｌ ｏｆ Ｇｕａｎｇｘｉ Ｔｅａｃｈｅｒｓ Ｅｄｕｃａｔｉｏｎ Ｕｎｉｖｅｒｓｉｔｙ ：Ｎａｔｕｒａｌ Ｓｃｉｅｎｃｅ Ｅｄｉｔｉｏｎ Ｖｏｌ ．３２ Ｎｏ ．２ ：１００１‐８７４３（２０１５）０２‐００１８‐０７ * 李徘菱，吕 艾，隆广庆 （ ， ５３００２３） ：致力于讨论第一类强奇异积分方程的快速Ｆｏｕｒｉｅｒ配置法．首先介绍了配置法基本框架，其次，根据核 ２ 函数 ｂ（ｘ ， ）的性质提出了一种矩阵压缩策略，使得矩阵的计算量由原来的 （ｎ ）减到 （ｎｌｏｇ ｎ），并且压缩后的 ｙ 逼近方程几乎保持原逼近方程的性质．基于快速Ｆｏｕｒｉｅｒ变换给出一种求积公式，最后用数值例子检验算法的有效 性． ：强奇异；积分方程；配置法；压缩策略 ：Ｏ１７５ ．３ ：Ａ １ 引 言 在处理许多物理学问题时，我们发现虽然这些问题的背景各不相同，但却有着完全一样的数学形 式．特别地，物理学上的静磁场问题、无电荷的经典问题或热传导问题都可以归结为Ｌａｐｌａｃｅ方程的边界 值问题．在考虑其第一类边值的内部问题时，利用双层位势理论，可导出含有强奇异核第一类积分方程． ［１ ，２］ ［３］ 许多传统的数值方法都可以用来求解第一类强奇异积分方程，比如：Ｎｙｓｔｒöｍ方法 ，谱方法 Ｇａｌｅｒ‐ ［４］ ［５ ，６ ，７］ ｋｉｎ方法 以及配置法 等． Ａｔｋｉｎｓｏｎ在文献［８］选取 ｎ 次三角多项式作为逼近函数，利用梯形求积公式定义离散的正交投影， 构造离散的Ｇａｌｅｒｋｉｎ方法求解带有对数奇异核的第一类积分方程．事实上，由于Ｆｏｕｒｉｅｒ基底的特性， 离散的Ｇａｌｅｒｋｉｎ方法与配置法是等价的．１９８９ 年，Ｗ．ＭｃＬｅａｎ和Ｗ．Ｌ．Ｗｅｎｄｌａｎｄ在［９］中对周期的拟 微分方程的数值解进行研究，构造了ｎ 次三角多项式为试探函数的Ｇａｌｅｒｋｉｎ方法的理论框架．此外，对 ｒ － ｒ 该算法进行了收敛性分析，证明当方程的解属于 Ｃ 空间时，其收敛阶为 （ｎ ｌｏｇ ｎ）．同年，Ｗ． ＭｃＬｅａｎ在［１０］证明了以ｎ 次三角多项式为试探函数的配置法收敛阶几乎是最佳的，即方程的解属于 ｒ － ｒ Ｃ 空间时，其收敛阶为 （ｎ ｌｏｇ ｎ）．随后，又在［１１ ，１２］，利用快速 Ｆｏｕｒｉｅｒ变化获得全离散的配置法 以及收敛分析结果．１９９６ 年，Ｊ．Ｓａｒｓｎｅｎ在配置法的基础上，提出了用两层网格迭代法求解Ｓｙｍｍ方程， 该方法在第一次迭代时，已获得最佳的收敛阶，第二次迭代又对该结果进一步修正． 然而，无论是传统的数值方法还是经过修正、迭代处理的数值方法，离散方程获得的系数矩阵都是 满秩的，这就意味着矩阵生成的工作量非常庞大．为了获得节约计算时间的数值方法，一些学者开始研 究用快速算法求解拟微分方程，其思想都是利用适当的截断策略将稠密矩阵压缩为稀疏矩阵． １９９２ 年，Ｗ．Ｄａｈｍｅｎ在文献［１３］利用小波基底建立Ｐｅｔｒｏｖ‐Ｇａｌｅｒｋｉｎ方法的理论框架，并对该方法 做了稳定性与收敛性分析．随后，Ｗ．Ｄａｈｍｅｎ又在文献［１４］提出，基于小波基的特点，矩阵元素具有一 定的衰减规律，从而提出相应的截断策略，得到Ｐｅｔｒｏｖ‐Ｇａｌｅｒｋｉｎ的小波快速算法．１９９７ 年，Ｃ．Ａ．Ｍｉｃｃｈ‐ ｅｌｌｉ在求解含有弱奇异核的第二类积分方程组时，利用具有小支集的多尺度小波基底建立小波Ｇａｌｅｒｋｉｎ ：２０１５‐０４‐２３ * ：广西研究生教育创新计划资助项目（ＹＣＳＺ２０１４１８５）；国家自然科学基金项目（１１