橢圓方程系數(shù)識(shí)別問題的正則化解

何 琴，王 謙

(蘭州交通大學(xué) 數(shù)理學(xué)院，蘭州 730070)

橢圓型方程[1-5]是偏微分方程領(lǐng)域內(nèi)的一個(gè)重要分支，是解決很多實(shí)際問題的有力工具.從整個(gè)區(qū)域上的不精確值識(shí)別橢圓方程Dirichlet問題中的系數(shù)的反問題，這引起了很多學(xué)者的極大關(guān)注[6-9].由于這些問題是不適定的，已經(jīng)有幾種穩(wěn)定的方法來解決他們，其中Tikhonov正則化方法是最著名的方法.雖然有很多論文專門研究這個(gè)問題，但是很少有人專門研究這些方法的收斂速度.文獻(xiàn)[10-14]研究了Bregman距離意義下線性不適定問題的凸變分正則化的收斂速度，獲得了全變分正則化的收斂速度.雖然有許多論文使用不適定問題的Tikhonov正則化，但很少有人專門研究收斂速度.Engl等[15]和Kanagaraj等[16]的論文專門討論了上述問題的Tikhonov正則化.這些作者使用輸出最小二乘法和非線性不適定問題的Tikhonov正則化，得到了某些源條件下的收斂速度.但是使用非凸泛函，很難找到全局極小值.在本文中，不使用輸出最小二乘法，而是將Tikhonov正則化方法應(yīng)用于求解原問題的新能量泛函，獲得該方法的收斂速度.

本文討論了具有齊次邊界Dirichlet橢圓方程的系數(shù)識(shí)別問題，陳述如下：

問題考慮如下橢圓型方程系數(shù)識(shí)別問題

(1)

u=0,在?Ω上.

(2)

其中:Ω是d,d≥1中具有Lipschitz邊界?Ω的開有界連通域，已知f∈L2(Ω)，c是未知系數(shù)，通過u在Ω上的觀測值來反演c.利用能量泛函和Tikhonov正則化方法，可獲得正則化解的收斂速度.構(gòu)造如下泛函

(3)

其中：ρ>0是正則化參數(shù)；c*是c的先驗(yàn)估計(jì).

1 問題設(shè)置

在本節(jié)中，主要證明了問題(1)～(2)存在唯一的解，給出解的估計(jì)，并提出了相應(yīng)的反問題.

(4)

A2={a∈(L∞(Ω))n|0<|a|≤a1},

(5)

(6)

其中:

(7)

(8)

2 Tikhonov正則化

(9)

(10)

由于問題是不適定的，使用Tikhonov正則化方法求解它.即求解最優(yōu)化問題

其中：ρ>0是正則化參數(shù)；c*是問題真實(shí)系數(shù)的先驗(yàn)估計(jì).

下面將介紹c*為最小范數(shù)解的概念，并概述了U(c)的一些性質(zhì).

引理1集合

在L2(Ω)-范數(shù)下是非空凸的，有界閉集.因此問題

有唯一解cz，它被稱為識(shí)別問題的c*是最小范數(shù)解.

對(duì)于此引理的證明，可參考文獻(xiàn)[17].

(11)

(12)

(13)

其中α>0.對(duì)任意的h∈L∞(Ω)，由Lax-Milgram引理，得到變分等式

(14)

‖η‖H1(Ω),

又因?yàn)?/p>

所以

‖h‖L∞(Ω)‖U(c)‖H1(Ω)‖η‖H1(Ω).

因此，

‖U(c)‖H1(Ω)‖η‖H1(Ω).

根據(jù)上一個(gè)不等式和式(6)，得下面不等式

(15)

由式(14)η的定義，得到

令v=U(c+h)-U(c)-η，由上式可得

所以

因此

‖η‖H1(Ω)‖U(c+h)-U(c)-η‖H1(Ω).

‖η‖H1(Ω)‖U(c+h)-U(c)-η‖H1(Ω).

根據(jù)上一個(gè)不等式和式(15)即得

引理證明完畢.

注1由于泛函Gzδ(c)的二階Fréchet導(dǎo)數(shù)的部分項(xiàng)

(c)h2·(U(c)-zδ),

對(duì)于此定理的證明，可參考文獻(xiàn)[17].

對(duì)于此定理的證明，可參考文獻(xiàn)[17].

3 收斂速度

因?yàn)長∞(Ω)=L1(Ω)*?L∞(Ω)*，對(duì)任意的c∈L∞(Ω)，有c∈L∞(Ω)*.于是對(duì)任意h∈L∞(Ω)有

(16)

和

‖c‖(L∞(Ω))*≤mes(Ω)‖c‖L∞(Ω).

是連續(xù)線性算子(見引理2).則U′(c)的對(duì)偶算子為

因此對(duì)任意的w*∈H-1(Ω)和h∈L∞(Ω)有

(17)

(18)

φρ=0,在?Ω上.

(19)

在H1(Ω)中的唯一解族.則存在一個(gè)常數(shù)M>0，使得對(duì)任意ρ∈(0,1)有

(20)

又因?yàn)?/p>

所以

使用分部積分公式和Cauchy-Schwarz不等式得

于是

(21)

根據(jù)假設(shè)條件，存在一個(gè)常數(shù)M1>0，使得對(duì)任意的ρ∈(0,1)有

由上式和不等式(21)，對(duì)任意的ρ∈(0,1)有

(22)

根據(jù)式(22)和Poincaré-Friedichs不等式，有式(20)成立.

引理證明完畢.

定理3假設(shè)存在函數(shù)w*∈H-1(Ω)使得

cz-c*=U′(cz)*w*.

(23)

則當(dāng)ρ→0,ρ～δ時(shí)，

于是

(24)

(25)

對(duì)不等式(24)右邊的第二項(xiàng)，根據(jù)等式(16)和(23)有

(26)

根據(jù)式(25)和式(17)可得

(27)

(28)

(29)

考慮下列橢圓方程Dirichlet問題

(30)

φρ=0,在?Ω上.

(31)

由分部積分得

(32)

(33)

由式(27)～(28)和(33)可得

(34)

根據(jù)不等式(29)～(34)得到

(35)

(36)

根據(jù)式(4)和(11)得

根據(jù)不等式(9)得到

(37)

由不等式(29)知，{φρ}ρ∈(0,1)在H1(Ω)-范數(shù)是有界的.因此由引理4知，存在一個(gè)只依賴于Ω的常數(shù)M>0，使得對(duì)任意ρ∈(0,1)有

(38)

根據(jù)式(38)和式(35)可得

(39)

聯(lián)立(35)～(39)推出

(40)

(41)

根據(jù)不等式(24)，(25)和(40)和(41)可得，當(dāng)ρ→0,ρ～δ時(shí)

即

定理證明完畢.