李金,桑瑜,张晓蕾,苏晓宁,屈金铮
(华北理工大学理学院,河北唐山063210)
奇异积分在边界元方法特别是自然边界元方法[1]中尤为重要,超奇异积分方程算法的研究同样应用广泛[2-3]。 圆周上的超奇异积方程通常出现在圆域或椭圆域的边值问题中,求解超奇异积分方程首先要计算超奇异积分,对超奇异积分进行有效的计算是求解超奇异积分方程面临的一个重要问题[4]。
近年来,超奇异积分的近似计算得到了广泛的研究[5-7]。 LINZ[8]研究了计算区间上超奇异积分的牛顿科特斯公式,并给出了奇异点在子区间中点时区间上二阶超奇异积分的复化梯形公式及相应的误差估计。 邬吉明等[9]将超奇异积分作为以奇异点为变量的函数,并提出相应的计算方法。 余德浩[10]研究了圆周上超奇异积分的近似计算公式,并得到误差估计。 对于区间上和圆周上超奇异积分的超收敛现象也得到了快速的发展[11-12]。 对于超奇异积分方程,研究成果相对较少,FENG 等[13]研究了基于矩形公式的圆周上超奇异积分方程的数值方法。对于配点法求解三阶超奇异积分方程的文献还较少。 LI 等[14]研究了圆周上三阶超奇异积分的超收敛性,并基于超收敛点构造配置格式。
文章基于梯形公式近似计算圆周上三阶超奇异积分,在误差泛函特殊函数为0(零点即为超收敛点)的基础上,选取超收敛点为配点构造线性方程组求解超奇异积分方程;通过引入正则化因子将奇异线性方程组变为正定线性方程组,完成圆周上三阶超奇异积分方程的数值求解;……