拋物方程的Landweber迭代正則化方法的后驗(yàn)誤差估計(jì)

摘要：考慮二維拋物方程Cauchy問(wèn)題的反問(wèn)題，該問(wèn)題是嚴(yán)重不適定的.首先，用Landweber迭代正則化方法得到該問(wèn)題的一個(gè)正則近似解，用Fourier變換求出該問(wèn)題的精確解；其次，在后驗(yàn)正則化參數(shù)的選取規(guī)則下，給出精確解和正則解之間的Holder型誤差估計(jì)，并使用更強(qiáng)的先驗(yàn)條件給出端點(diǎn)x=1處的誤差估計(jì)；最后，給出數(shù)值實(shí)例說(shuō)明該方法的有效性.結(jié)果表明，該方法比已有方法收斂速度更快。

關(guān)鍵詞：Cauchy問(wèn)題；不適定問(wèn)題；Landweber迭代正則化；誤差估計(jì)；后驗(yàn)估計(jì)中圖分類號(hào)：O241文獻(xiàn)標(biāo)志碼：A文章編號(hào)：1671-5489（2024）05-1113-09

Posterior Error Estimation of Landweber Iterative Regularization Method for Parabolic Equations

SHEN Yu，XIONGXiangtuan

（College of Mathematics and Statistics，Northwest Normal University，Lanzhou 730070，China）

Abstract：We considered the inverse problem of Cauchy problem of two dimensional parabolic equations，which was seriously ill-posed.Firstly，a regular approximate solution of the problem was obtained by using Landweber iterative regularization method，and Fourier transform was used to obtain the exact solution of the problem.Secondly，the H?lder type error estimation between the exact solution and the regular solution was given under the selection rules of the posterior regularization parameters，and stronger prior conditions were used to give the error estimation at the end point x=1.Finally，numerical examples were given to demonstrate the effectiveness of the proposed method.The results show that the proposed method has a faster convergence rate than existing methods.

Keywords：Cauchyproblem;ill-posedproblem;Landweber iterative regularization;errorestimation;posterior estimation

0引言

假設(shè)（y，t）∈L2（RXR）為給定函數(shù)，且當(dāng)lt;0時(shí)，（y，）=0.定義g（y，）的范數(shù)為收稿日期：2023-10-25.網(wǎng)絡(luò)首發(fā)日期：2024-07-05.

考慮二維拋物型方程的Cauchy問(wèn)題：

其中函數(shù)u（x，，和a（）都屬于L2（RXR），精確數(shù)據(jù)函數(shù)（）及其噪聲數(shù)據(jù)函數(shù)48（，）滿足：

這里δgt;0表示輸入數(shù)據(jù)的噪聲水平.記f（y，t）：=u（1，y，t）且有一個(gè)先驗(yàn)界

其中Egt;0是一個(gè)常數(shù).需從給定Cauchy數(shù)據(jù)[u，u]在x=0處的值反演u（，，）（0≤x≤1）.

二維拋物方程Cauchy問(wèn)題的逆問(wèn)題，即在二維拋物方程中用x=0處的值去反演x∈（0，1]處的函數(shù)值，應(yīng)用廣泛.在實(shí)際應(yīng)用中，給定源參數(shù)在固定位置（x，y）的測(cè)量數(shù)據(jù)的傳熱過(guò)程可用問(wèn)題（1）描述，物理學(xué)中的許多分支，如流體動(dòng)力學(xué)、等離子體動(dòng)力學(xué)、光學(xué)、場(chǎng)論、凝聚態(tài)物理等領(lǐng)域也常涉及這類方程.而這類問(wèn)題是嚴(yán)重不適定問(wèn)題，其數(shù)據(jù)的一個(gè)小擾動(dòng)可能在解中引起極大誤差，因此需要一種合適的正則化方法解決此問(wèn)題.目前，關(guān)于拋物方程Cauchy問(wèn)題正則化的理論研究已有很多結(jié)果.對(duì)于一維情形，Knabner等[2]利用Fourier變換對(duì)一維方程的精確解給出了最優(yōu)穩(wěn)定性估計(jì)；文獻(xiàn)[3]用Meyer小波正則化方法討論了該問(wèn)題，并給出了先驗(yàn)的正則化參數(shù)選取規(guī)則及誤差估計(jì)；文獻(xiàn)[4]將Meyer小波正則化方法推廣到N維情形，并給出了先驗(yàn)的正則化參數(shù)選取規(guī)則及誤差估計(jì).對(duì)于二維情形，文獻(xiàn)[5-6]提出了兩種正則化策略：基于方程的修正核方法和基于高頻分量的截?cái)喾椒ǎ晃墨I(xiàn)[7]提出了一種改進(jìn)的擬邊值正則化方法處理不適定問(wèn)題，通過(guò)選擇適當(dāng)?shù)恼齽t化參數(shù)和引入一些技術(shù)不等式，得到了近似解與其精確解之間精確的誤差估計(jì)；文獻(xiàn)[8]從不同角度提出了兩種正則化方法，Tikhonov方法和Fourier截?cái)喾椒?，給出并證明了精確解與其正則化近似之間的收斂估計(jì).Landweber迭代正則化方法也可用于分析和解決逆問(wèn)題的不適定性[912]，本文考慮用Landweber迭代處理二維拋物方程Cauchy問(wèn)題，給出后驗(yàn)的正則化參數(shù)選取規(guī)則及誤差估計(jì)，并通過(guò)使用更強(qiáng)的先驗(yàn)條件得到端點(diǎn)x=1處的誤差估計(jì).

1預(yù)備知識(shí)及不適定性分析

本文考慮問(wèn)題（1）的不適定性并求解該問(wèn)題.對(duì)輸入數(shù)據(jù)（y，t），定義它的Fourier變換為

通過(guò)對(duì)問(wèn)題（1）的變量y和t進(jìn)行二維Fourier變換，在頻率空間中得到如下二階常微分方程的Cauchy問(wèn)題：

從而可得問(wèn)題（5）的解：

其中0=0（n，E）是i+2的主平方根.易證0的實(shí)部和虛部由下式給出：

且可得｜（）｜≥、（）≥0.

由Fourier逆變換，可得問(wèn)題（5）的精確解：

特別地，由式（6）可得

下面假設(shè)問(wèn)題（5）的解滿足如下先驗(yàn)界：

將定義為精確解的核由于ok）在gt;0是無(wú)界的因此數(shù)據(jù)中的小錯(cuò)誤可能會(huì)爆炸并完全破壞解決方案.此外，高頻分量中的誤差被該因子exp2++放大，因此，不能用經(jīng)典的數(shù)值方法計(jì)算這個(gè)問(wèn)題.本文采用Landweber迭代正則化方法，即通過(guò)修改核降低問(wèn)題（1）的難度.

引理1[5]若a≥3≥0，x≥0，o=sign（），∈R，則有

2 Landweber代正則化方法

問(wèn)題（5）是不適定的.如果想要恢復(fù)其解的穩(wěn)定性，需要使用正則化方法.本文采用Landweber迭代正則化方法得到式（5）的正則化解.

根據(jù)式（6），可得

定義算子K：u→，因?yàn)槭浅朔ㄋ阕?，所以式?）可改寫(xiě)為如下算子方程：

則有KK*=“K”2，采用Landweber正則化方法求解Ku=的正則化解.將算子方程Ku（x，n，E）=（）替換為算子方程

得到迭代格式如下：

其中a是松因子且足TK

設(shè)算子Rm：L2（Ω）→L2（Ω）為

則用可測(cè)數(shù)據(jù)（，）進(jìn)行Landweber迭代求得解為

從而可得

再利用Fourier逆變換，可得問(wèn)題（1）的Landweber迭代正則化解為

3后驗(yàn)參數(shù)選取下的誤差估計(jì)

下面將在后驗(yàn)正則化參數(shù)的選取規(guī)則下，給出問(wèn)題（1）的誤差估計(jì)，設(shè)gt;1為固定常數(shù)，并在m=m（δ）∈N。第一次出現(xiàn)時(shí)停止算法：

其中‖‖≥rò.

引理2令B（m）≤Ri（x，ミ）-（，），則下列結(jié)論成立：

1）lim B（m）=llp（n，ll

2）limβ（m）=0;

3）B（m）是一個(gè)連續(xù)函數(shù)；

4）對(duì)任意m∈（0，+∞），B（m）是嚴(yán)格單調(diào)遞減的函數(shù).

證明：由

可知β（m）滿足以上4個(gè)條件.

引理3對(duì)任意x∈（0，1），正則化參數(shù)m滿足：

證明：根據(jù)式（16），有

從而有

因?yàn)椤?-a”K“2”lt;1，所以由式（22）可得‖KRm1-I‖≤1.根據(jù)式（19），可得

則

因此

此外，有

從而可得

由式（23），（24）可得

因此

定理1設(shè)（x，yt）是問(wèn)題（1）的精確解，un（y）是問(wèn)題（1）帶有噪聲數(shù)據(jù)的正則解.若噪聲假設(shè)式（2）和先驗(yàn)界式（9）成立，正則化參數(shù)由式（21）選擇，則對(duì)任意的0lt;xlt;1，如下誤差估計(jì)式成立：

證明：由三角不等式和Parseval等式得

記I1=‖i（x，ミ）ーi（x，ミ）‖，I2=in（x，ミ）i（x，）

首先估計(jì)I1.由式（21）及引理3可知

根據(jù)Bernoulli不等式可知

則‖um（x，n，）-um（x，n，）‖≤√amδ，從而可得

其次估計(jì)I2.根據(jù)Holder不等式可得

故

根據(jù)式（27），（28）可得

注1定理1只考慮了區(qū)間內(nèi)的誤差估計(jì)，而未考慮端點(diǎn)x=1處的誤差估計(jì)，因此只能說(shuō)明它是有界的，不能說(shuō)明它是收斂的，如果想獲得精確解和正則化解在x=1處的誤差估計(jì)，必須引入更強(qiáng)的先驗(yàn)假設(shè).為給出x=1處的誤差估計(jì)，先給出如下先驗(yàn)界：

其中（，）是P范數(shù)，且P=0是L2范數(shù).

引理4假設(shè)式（2）和先驗(yàn)條件（9）成立，對(duì)Pgt;0先驗(yàn)界（30）成立，則在x=1處的正則化參數(shù)m滿足：

證明令=oE則

從而

另一方面，有

其中/0=-（ー））

故

由式（32），（33）可得

從而可得

定理2設(shè)um3（x，y，t）為問(wèn)題（1）在x=1處的正則化解，若噪聲假設(shè)（2）和先驗(yàn)條件（30）對(duì)Pgt;0成立，將式（31）的解作為x=1處的正則化參數(shù)，則誤差估計(jì)為

證明：根據(jù)三角不等式和Parseval等式，由引理4可得l（1y.）-（1，y.t）=in（1）-i（1.，）≤

根據(jù)Holder不等式可得

因此

從而

其中C+-）\"+2證

注2對(duì)比文獻(xiàn)[8]中在端點(diǎn)/處的誤差估計(jì)結(jié)果

本文結(jié)果為

顯然本文結(jié)果的收斂速度比文獻(xiàn)[8]更快.

4數(shù)值實(shí)驗(yàn)

下面給出具體實(shí)例說(shuō)明本文方法的可行性和有效性.本文方法的數(shù)值過(guò)程基于離散快速Fourier變換和逆離散快速Fourier變換計(jì)算正則化解，并應(yīng)用引理3、引理4選擇正則化參數(shù).數(shù)據(jù)函數(shù)。不具有一般性，必須在有限域上采樣，假設(shè)為[-4，4]×[-4，4]，數(shù)據(jù)函數(shù)按101×101點(diǎn)等距采樣.隨機(jī)噪聲加到準(zhǔn)確的輸入數(shù)據(jù)φ上，即隨機(jī)數(shù)=+erand（size（φ）），其中函數(shù)rand（·）是由MATLAB在[0，1]上產(chǎn)生的隨機(jī)數(shù)組，將并得到噪聲數(shù)據(jù)，其中和滿足式（2）.

例1[8]令φ（y，t）=π2e-（y2+2/4），a=0.1，e=0.01，求方程（1）在x=0.5處的解.

根據(jù) Fourier變換由式 （8）可得當(dāng) φ（y，t）= π2e-（y 2+t 2/4）時(shí)在a=0.1，ε=0.01條件下的二維拋物方程 Cauchy問(wèn)題（1）的精確解u（0.5，·，·）， 結(jié)果如圖1所示.用噪聲數(shù)據(jù)函數(shù)φδ（·，·）且未使用正則化方法所求的解uδ （0.5，·，·）如圖2所示.為恢復(fù)解的穩(wěn)定性，用Landweber迭代正則化方法得到問(wèn)題 （1）在例 1 條件下的正則化解um，δ（0.5，·，·）如圖3所示.由圖1和圖2可見(jiàn)，問(wèn)題（1）是非常不適定的，小的誤差即會(huì)引起解的巨大振蕩.由圖3可見(jiàn)，問(wèn)題（1）經(jīng)過(guò)Landweber正則化后得到的正則解與精確解非常接近，說(shuō)明了該方法的有效性.此外，對(duì)比文獻(xiàn)[8]的實(shí)例結(jié)果表明，在選取參數(shù)一致的情況下本文實(shí)驗(yàn)結(jié)果優(yōu)于文獻(xiàn)[8]的結(jié)果.

綜上所述，本文解決了不適定二維拋物方程Cauchy問(wèn)題.先通過(guò)Fourier變換求出了問(wèn)題的精確解，再采用Landweber迭代正則化方法得到正則解，然后在后驗(yàn)正則化參數(shù)的選取規(guī)則下，得到了精確解與正則解之間的Holder型誤差估計(jì)，并且使用比已有結(jié)果更強(qiáng)的先驗(yàn)條件給出了端點(diǎn)x=1處的誤差估計(jì)，所得誤差估計(jì)的收斂速度比已有結(jié)果更快.數(shù)值實(shí)驗(yàn)結(jié)果表明了本文方法的有效性.

參考文獻(xiàn)

[1]DEHGHAN M.Fourth-Order Techniques for Identifying a Control Parameter in the Parabolic Equations[J].International Journal of Engineering Science，2002，40（4）：433-447.

[2]KNABNER P，VESSELLA S.The Optimal Stability Estimate for Some Ill-Posed Cauchy Problems for a Parabolic Equation[J].Mathematical Methods in the Applied Sciences，1988，10（5）：575-583.

[3]HAO D N，SCHNEIDERS A，REINHARDT H J.Regularization of a Non-characteristic Cauchy Problem for a Parabolic Equation[J].Inverse Problems，1995，11（6）：1247-1263.

[4]KNOSOWSKI Y，VON LIERES E，SCHNEIDER A.Regularization of a Non-characteristic Cauchy Problem for a Parabolic Equation in Multiple Dimensions[J].Inverse Problems.1999，15（3）：731-743.

[5]QIAN Z.FU C L.Regularization Strategies for a Two-Dimensional Inverse Heat Conduction Problem[J].Inverse Problems，2007，23（3）：1053-1068.

[6]QIAN Z.FU C L，XIONG X T.A Modified Method for Determining the Surface Heat Flux of IHCP[J].Inverse Problems in Science and Engineering，2007，15（3）：249-265.

[7]CHENG W.ZHAO Q.A Modified Quasi-boundary Value Method for a Two-Dimensional Inverse Heat Conduction Problem[J].Computersamp;.Mathematics with Applications，2020，79（2）：293-302.

[8]QIAN Z.Regularization Methods for a Cauchy Problem for a Parabolic Equation in Multiple Dimensions[J].Journal of Inverseamp;.Ill-Posed Problems，2009，17（9）：891-911.

[9] SCHERZER O.Convergence Criteria of Iterative Methods Based on Landweber Iteration for Solving Nonlinear Problems[J].Journal of Mathematical Analysis and Applications，1995，194（3）：911-933.

[10]XIONG X T，XUE X M，QIAN Z.A Modified Iterative Regularization Method for Ill-Posed Problems[J].Applied Numerical Mathematics，2017，122：108-128.

[11]JOSE J，RAJAN M P.A Simplified Landweber Iteration for Solving Nonlinear Ill-Posed Problems[J].International Journal of Applied and Computational Mathematics，2017，3：S1001-S1018.

[12]CHEN Y G，Y ANG F，DING Q.TheLandweber Iterative Regularization Method for Solving the Cauchy Problem of the Modified Helmholtz Equation[J/OL].Symmetry，（2022-06-11）[2024-03-19].https：/doi.org/10.3390/sym14061209.

（責(zé)任編輯：趙立芹）