尹忠旗
(四川師范大學數(shù)學與軟件科學學院,四川成都610066)
令X、Y、Z為可分的Hilbert空間,A:X→Z是一對一的有界線性算子.B:Y→Z是一個線性算子,它可能是無界的.考慮如下的算子方程
考慮算子A的值域不閉的情況,即這表明A的逆算子A-1:R(A)→X是不連續(xù)的.這是一個不適定的問題[1-3]:初始數(shù)據(jù)g的微小改變可能導致方程(1)的解f的巨大改變.于是,需要尋求一種正則化技術(shù)來確定方程(1)在噪聲水平為δ,精確數(shù)據(jù)g0∈Y,Bg0∈R(A)下的穩(wěn)定的近似解f0.所謂的噪聲水平δ是指‖gδ-g0‖Y≤δ.引入Tikhonov正則化[4-6]:
其中,實數(shù)α叫做正則化參數(shù).記(2)式的唯一極小值解為fα,δ,它是問題(2)的正則化解,滿足方程
感興趣的是與最優(yōu)收斂速度有關(guān)的正則化參數(shù)α的先驗選擇策略.文獻[7]討論了B=I,Y=Z,即Af=g;證明了對某些合適的函數(shù)Ψ有收斂速度‖fα,δ-f0‖X=O(Ψ0(δ)).文獻[8]考慮了算子方程Ax=y在不同光滑性條件下的正則化問題,其中方程的解屬于某一個自伴算子G:X→X的值域R(G)且有包含關(guān)系R(G)?φ(A*A),而φ是一個嚴格單調(diào)增加的連續(xù)函數(shù),它滿足初始條件φ(0)=0.關(guān)于先驗假設(shè)情況下的收斂問題,更多的文獻可以參考文獻[9-13]等相關(guān)專著和文章.本文的新穎之處在于所考慮的情形與值域包含有關(guān),即R(A)?R(B)而且B可能是無界的.容易知道,本文的結(jié)果已經(jīng)包含了文獻[7]的結(jié)果.在本文中,由于B可能無界,因此要求不同于以往的方法來處理相應(yīng)的收斂速度問題.還假設(shè)算子A*B可以分解為2個有界算子的復合.
先給出文中要使用的一些假設(shè),然后給出本文的主要結(jié)果.
令I(lǐng)為指標函數(shù)的集合,即
令G為X上的一對一的,正的緊自伴算子,它具有對應(yīng)于其特征值的規(guī)范的正交基{φn}n∈N,它的特征值是以單調(diào)增加的方式排列的
用{Xρ(G)}ρ∈I表示可變尺度的 Hilbert空間[14-15].Xρ(G)是如下集合的閉包:
在本節(jié)中,給出幾個預備引理.第1個引理刻畫了方程(3)在沒有擾動的情況下的正則解與在有δ擾動的情況下的正則解之間的誤差.第2個引理描述了方程(1)的解與其正則化解,即與方程(3)在沒有擾動下的解之間的誤差.它們的證明可以參考文獻[8].
引理 2.1令fα,δ是問題(2)的極小值解,fα是算子方程(A*A+αI)f=A*Bg0的解.那么有如下估計成立
[1]Kirsch A.An Introduction to the Mathematical Theory of Inverse Problems[M].New York:Springer-Verlag,2011.
[2]高常忠,宋惠元.一類確定線性偽拋物型方程邊界值的反問題[J].四川師范大學學報:自然科學版,2004,27(6):603-606.
[3]何歡,孫合明,左環(huán).對稱矩陣特征值反問題的最佳逼近解的一種數(shù)值解法[J].四川師范大學學報:自然科學版,2012,35(4):473-477.
[4]Tikhonov A N,Arsenin V Y.Solutions of Ill-Posed Problems[M].New York:John Wiley&Sons,1977.
[5]Engl H W,Hanke M,Neubauer A.Regularization of Inverse Problems[M].Dordrecht:Kluwer Academic Publishers,1996.
[6]Groetsch C W.The Theory of Tikhonov Regularization for Fredholm Integral Equations of the First Kind[M].Boston:Pitman,1984.
[7]Hofmann B,Yamamoto M.Convergence rates for Tikhonov regularization based on range inclusions[J].Inverse Problems,2005,21(5):805-820.
[8]B?ttcher A,Hofmann B,Tautenhahn U,et al.Convergence rates for Tikhonov regularization from different kinds of smoothness conditions[J].Appl Anal,2006,85(5):555-578.
[9]Baumeister J.Stable Solution of Inverse Problems[M].Braunschweig:Friedr Vieweg&Sohn,1987.
[10]Colton D,Kress R.Inverse Acoustic and Electromagnetic Scattering Theory[M].Berlin:Springer-Verlag,1992.
[11]Hofmann B.Regularization for Applied Inverse and Ill-Posed Problems[M].Leipzig:BSB B G Teubner Verlagsgesellschaft,1986.
[12]Hohage T.Regularization of exponentially ill-posed problems[J].Numer Funct Anal Optim,2000,21(3/4):439-64.
[13]Mair B A.Tikhonov regularization for finitely and infinitely smoothing operators[J].SIAM J Math Anal,1994,25(1):135-147.
[14]Hegland M.Variable Hilbert scales and their interpolation inequalities with applications to Tikhonov regularization[J].Appl Anal,1995,59(1/2/3/4):207-223.
[15]Neubauer A.Numerical realization of an optimal discrepancy principle for Tikhonov regularization in Hilbert scales computing[J].Computing,1987,39(1):43-55.