偏微分方程边值反问题的数值方法研究

易苗,刘扬

(1.武汉大学数学与统计学院,湖北武汉430072)

(2.武汉理工大学理学院,湖北武汉430070)

偏微分方程边值反问题的数值方法研究

易苗1,刘扬2

(1.武汉大学数学与统计学院,湖北武汉430072)

(2.武汉理工大学理学院,湖北武汉430070)

本文研究了奇异积分方程在反边值问题中的应用问题.利用圆周上的自然积分方程及其反演公式,把Laplace方程的边值反问题转化为一对超奇异积分方程和弱奇异积分方程的组合,通过选取三角插值近似奇异积分的计算并构造相应的配置格式,并使用Tikhonov正则化方法求解所得到的线性方程组.数值实验表明了该方法的有效性.

边值反问题;奇异积分方程;三角插值;Tikhonov正则化

1 引言

考虑二维平面圆形区域Ω上的Laplace方程的边值反问题

Laplace方程的边值反问题也称为Cauchy问题,这类问题在工程技术如：地质勘探、卫星测量、空间遥感、目标识别以及医学成像等领域有着深刻的应用背景.根据已有的关于解的信息反演初值或边界条件通常是很困难的,目前理论上还没有统一的方法,因此偏微分方程的反问题研究具有很大的挑战.

问题1.1已经被证明了是严重不适定问题,即初值的任意很小的改变都会引起解的巨大改变.Belgacem[1]对Cauchy问题的严重不适定性进行了理论分析;Calder´on和Engl[2,3]研究了Cauchy问题弱解的存在性与唯一性;文献[4]证明Cauchy问题是条件稳定的,即在一个附加的有界条件下,问题的解连续依赖初始数据.关于Laplace方程边值反问题的数值方法目前已经有了一些成果[5-7],其中Li[7]利用圆周上Hilbert变换的特殊性质构造了Laplace方程边值反问题的积分方程解法.本文主要研究奇异积分在偏微分方程边值反问题中的应用,构造求解这类问题的一种新的数值方法,并用数值实验验证方法的有效性.

2 边值反问题的数值方法

设u0=u|Γ0和g0=un|Γ0是部分边界Γ0上的已知函数Z,u1=u|Γ1和g1=un|Γ1是部分边界Γ1上的待定函数,u和un分别满足周期性条件=0和相容性条件.根据自然边界归化原理[8]可知,单位圆内区域Ω上的调和方程满足自然积分方程

其反演公式

是一个弱奇异积分方程.由方程(2.1)和(2.2)可得

2.1 奇异积分的数值方法

考虑超奇异积分

其中φ(θ)是以2π为周期的函数.对于任意的周期函数φ(θ),其n阶三角多项式插值可定义为

用φn(θ)分别替换(2.5)和(2.6)式中的φ(θ),得到超奇异和弱奇异积分的计算公式分别为

其中

当φ(θ)足够光滑时,公式(2.8)和(2.9)是非常精确的,有关它们的详细讨论可以参见文献[9].

2.2 积分方程的配置法

并假设它满足相容性条件利用三角插值积分公式(2.8)逼近超奇异积分方程(2.12),取配置点βi(0≤i≤2n)并排列所得到的方程,由此可以得到如下线性方程组

或写成矩阵的形式

并且φi(i=0,1,···,2n)表示函数φ在配置点βi处的近似值.矩阵A一般是奇异的,为了求解线性方程组(2.13)或(2.14),通常需要选取某种正则化方法去求解.

对于弱奇异积分方程

利用三角插值积分公式(2.8)逼近超奇异积分方程(2.15),取一系列配置点βi(0≤i≤2n)并排列所得到的方程,由此可以得到如下线性方程组

或写成矩阵的形式

其中

并且φi(i=0,1,···,2n)表示函数φ在配置点βi处的近似值.矩阵B一般也是奇异的,为了求解它也需要正则化之后才能求解.

2.3 正则化

利用配置法求解(2.3)或(2.4)会得到两个线性方程组,其矩阵形式为

其中u是方程(2.3)或(2.4)中未知量近似值的列向量,F是由方程(2.3)或(2.4)右端表达式计算的边界Γ0上的已知值.方程(2.18)通常是一个严重的病态方程组,假如边界Γ0上没有足够多的信息或方程的右端项含有一定的噪声,无法用方程(2.18)求解另一部分边界Γ1上所需要的信息.

为了求解方程(2.18),采用Tiknonov正则化方法.令

其中α＞0是正则化参数,L是离散的一阶线性微分算子.上式的解可以写成矩阵的形式如下

求解这个方程即可得到边值反问题的近似解.这里关于最优正则化参数α的选取可以参考文献[10].

3 数值实验

下面给出两个数值算例来验证上节中提出的边值反问题的数值计算方法.

数值算例1考虑边值反问题(1.1),其中u0=3cos2x+4sin2x,un=6cos2x+8sin2x,分别计算了Γ0=[0,π],Γ0=[0,π/2]和Γ0=[0,π/4]三种情况,其结果见图1-3.

图1：当Γ0=[0,π]时边值反问题的数值解与精确解：(a)u0(s)的图像,(b)un(x)的图像.

图2：当Γ0=[0,π/2]时边值反问题的数值解与精确解：(c)u0(s)的图像,(d)un(x)的图像.

图3：当Γ0=[0,π/4]时边值反问题的数值解与精确解：(e)u0(s)的图像,(f)un(x)的图像.

图4：当Γ0=[0,π/4]时边值反问题的数值解与精确解：(g)u0(s)的图像,(h)un(x)的图像.

图5：当Γ0=[0,π/4]时边值反问题的数值解与精确解：(i)u0(s)的图像,(j)un(x)的图像.

从图1可以看出,当Γ0=[0,π],即已知边界占整个边界一半时,由超奇异积分方程和弱奇异积分方程配置法复原得到的未知边界上的信息与其理论值高度一致,这说明本文所提出的算法是有效的.从图2和图3的结果可以看出,当已知边界小于一半时,由配置法复原得到的未知边界上的值与理论值有一定的出入,但近似计算结果仍然复原了原边界曲线上的波动特征.

在实际工程中,通常测量数据都会有一定的误差,即Γ0上的已知值会含有一定的噪声.由于边值反问题已经被证明是一个严重的病态问题,任何微小的扰动都会引起解的巨大改变.下面给出一个算例测试噪声对数值解得影响.

数值算例2考虑例1中的边值反问题.分别计算了Γ0=[0,π/2]和Γ0=[0,π/4]两种情况,同时对Γ0上的已知信息u0和g0加入2%的随机误差,其结果见图4-5.

从图4和图5的结果可以看出,由积分方程配置法复原的数值解与精确解匹配还是非常好.同时把图4和图5的结果与没有添加噪声的复原结果相比,其变化并不太大.这也从侧面印证了算法的有效性和可行性.

[1]Belgacem F B.Why is the Cauchy problem severely ill-posed[J].Inv.Prob.,2007,23(2)：823-826.

[2]Calder´on A P.Uniqueness in the Cauchy problem for partial dif f erential equations[J].Amer.J. Math.,1958,80(1)：16-36.

[3]Engl H W,Leitao A.A Mann iterative regularization method for elliptic Cauchy problems[J].Numer. Funct.Anal.Optim.,2001,22：861-864.

[4]Alessandrini G,Rondi L,Rosset E,Vessella S.The stability for the Cauchy problem for elliptic equations[J].Inv.Prob.,2009,25：123-134.

[5]Cao H,Pereverzev S V.The balancing principle for the regularization of elliptic Cauchy problems[J]. Inv.Prob.,2007,23(5)：1943.

[6]Chapko R,Johansson T.An iterative method based on boundary integrals for elliptic Cauchy problems in semi-inf i nite domains[J].J.Publ.Date.,2009,7(1)：1-12.

[7]Li H J,Feng X S,Xiang J,et al.New approach for solving the inverse boundary value problem of Laplace’s equation on a circle：Technique renovation of the Grad-Shafranov(GS)reconstruction[J]. J.Geophy.Res.：Space Phys.,2013,118(6)：2876-2881.

[8]Yu D.Natural boundary integral method and its applications[M].Holland：Springer Science& Business Media,2002.

[9]Belotserkovsky S M,Lifanov I K.Method of discrete vortices[M].Russian：CRC Press,1992.

[10]Ramm A G.Dynamical systems method for solving operator equations[J].Commun.Nonl.Sci. Numer.Simul.,2004,9(4)：383-402.

[11]孙萍,冯晓莉.一种求解修正的Helmholtz方程Cauchy问题的数值方法[J].数学杂志,2011,31(4)：756-762.

NUMERICAL METHODS FOR SOLVING INVERSE BOUNDARY VALUE PROBLEM OF PARTIAL DIFFERENTIAL EQUATION

YI Miao1,LIU Yang2
(1.School of Mathematics and Statistics,Wuhan University,Wuhan 430072,China) (2.School of Science,Wuhan University of Technology,Wuhan 430070,China)

In this paper,we study the application problem of singular integral equation in the inverse boundary value problem.Using the natural integral equation and its inversion formula on a circle,we transformed the Laplace equation inverse boundary value problem into a combination of a hypersingular integral equation and a weakly singular integral equation,then construct the corresponding collocation scheme based on the trigonometric interpolation,and use the Tikhonov regularization to solve the resulting linear equations.Numerical experiments show the ef f ectiveness of the method.

inverse boundary value problem;singular integral equation;trigonometric interpolation;Tikhonov regularization

O241.83

A

0255-7797(2017)05-1040-07

2016-01-09接收日期：2016-04-15

国家自然科学基金(11201358);中央高校基本科研业务费专项资金资助(2015IA007).

易苗(1990-),女,湖北宜昌,硕士,主要研究方向：偏微分方程数值解.

2010 MR Subject Classif i cation：65D30;65N21