книги / Регулярные методы решения некорректно поставленных задач
..pdf5. Остановимся на некоторых примерах, приводящих к решению основной задачи.
П р и м е р 4. Пусть Н - L 2 [а, £] - пространство функций, сум-
ъ
мируемых по Лебегу с квадратом: И“ Н я “ / «* (х)dx, 2LD — мно-
a
жество функций из L 2[a, b], л-я производная которых суммируема
сквадратом [88]. Положим
d |
|
|
|
|
Аи= j к(х, £)м(£)с?£, |
Аи G F - L 2[C, d ] , |
|
(И ) |
|
С |
|
|
|
|
где к(х,%) |
— непрерывная функция своих |
аргументов. Тогда |
||
задача: найти |
u E D A u = f |
(fEF ) является некорректно |
постав |
|
ленной [90]. Можно положить Lu = d nu/dxn, |
и ED. Тогда |
основ |
ная задача заключается в вычислении л-й производной на решениях интегрального уравнения (11). Если для любого полинома Рп_ i(jc) степени не выше л - 1 из соотношения АРп_ х = 0 следует, что Рп__1(*) = 0»то выполнено условие дополнительности.
П р и м е р |
5. Пусть, как и выше, H = L2 [at b] , Lu = d nu/dxn, |
|
и ED, |
Аи = и. Основная задача заключается в вычислении элемента |
|
g - Lu, |
и ED |
(решением уравнения Аи = и = й является, очевидно, |
элемент и). Так как допускаются приближения к и из пространства L2 [а, Ъ], то эта задача о дифференцировании является некорректно поставленной.
П р и м е р 6. Снова пусть Н = L 2 [а, Ь]. Положим
Аи=(и(х1) , -----и{хт) ) т, uED,
где Х{ (/ = 1,. .. , т) |
— некоторая сетка узлов |
(которые могут и |
||
совпадать) |
на отрезке |
[а, Ъ]. Задача гладкой интерполяции, заклю |
||
чающаяся в определении элемента |
|
|||
uED: |
Аи =/, |
f E R m, |
(12) |
является, очевидно, некорректно поставленной, так как заведомо нарушается, например, условие единственности (а также условие существования при совпадающих узлах) интерполирующей функции.
Положим Lu = dnujdxn, |
и ED. Тогда |
основная задача заклю |
||
чается в |
выборе такой интерполирующей функции й = м ( х ) Е Д |
|||
для которой |
|
|
|
|
inf |
II L u - g II;■ |
= || Lu - g II , , |
g e L 2[a,b], |
(13) |
u ^ D : A u = f |
|
|
|
т.е. задача сводится к сплайн-интерполяции [70].
2t
Нетрудно дать обобщение задачи интерполирования и на много мерный случай. Общая теория сплайнов будет рассмотрена в гл. 4.
Другие применения основной задачи будут показаны по ходу изложения.
§ 2. Аппроксимация решения основной задачи
1. Для построения L -псевдорешений основной задачи необходи мо явное задание множества Uf, которое не всегда возможно осу ществить. Это затрудняет поиск //-псевдорешений. Естественно попытаться освободиться от указанного недостатка.
Определим обобщенный параметрический функционал Ти хонова'.
Фа [м] = I\ А и - f \\ F2 + а|| L u - g \ \ G,2 |
u&D, |
|
где а > 0 |
- параметр регуляризации, а / |
Е F и g € G - заданные |
элементы. |
Очевидно, Фа [м] = Фа [и\ f\A, g, L, D \. |
Рассматривается задача отыскания регуляризованных решений,
т.е. элементов йа таких, что |
|
inf Ф<*М =Фа [и«]- |
( 1) |
uGD |
|
Т е о р е м а 2. Пусть выполнены основные предположения (тео рема 1)- Тогда при любом а > 0 решение задачи ( 1) существует
иединственно.
До к а з а т е л ь с т в о . Нетрудно видеть, что справедливо тождество
1 |
1 |
\ |
и 1 +м2 |
|
= - |
Ф«[И|] + --Ф«[И2] - Фа |
|
|
|
для любых их, и2 £ А |
Если это решения задачи (1), то |
|||
II / |
«1 - и2 \ II2 |
L / и, |
- и 2\ \ \ г Л2 |
Г«1 + м21 |
1^1—5—) 1/ “ |Ч“Г-)1о'"“-ф* г Н <0-
ибо (их + и2)/2 G D. Но тогда \и х - и2 1= 0, поэтому, воспользо вавшись условием дополнительности, получаем их - м2Итак, решение ( 1), если оно существует, определяется однозначно.
Существование регуляризованного решения иа устанавливается аналогично теореме 1.
2. Пусть /= /; |
g - g и йа - |
решение задачи (1). |
Т е о р е м а 3. |
Если а 0, |
то I и& - и | 0. |
22
З а м е ч а н и е . Очевидно, lim || йа - и \\ц.= О в силу условия
а—*О
взаимной дополнительности.
Теорема 3 показывает, что регуляризованные решения иа при малых а аппроксимируют решение основной задачи. Так как при минимизации ( 1) не требуется знать множество Uf, то определение регуляризованных решений может быть эффективно выполнено.
Для доказательства теоремы 3 требуется ввести следующее опре деление. Операторы А и L называются совокупно слабо замкнуты ми на D, если из соотношений
Un e£), |
сл |
сл |
сл |
(2) |
ип — * /о , |
Лип— ►/„, |
Lun— >gо при |
||
где символ |
сл |
|
сходимость в соответствующем |
|
— ►означает слабую |
пространстве, следует: м0 6 D, Аи0 = / 0, Lu0 = £ 0- Множество D здесь необязательно выпуклое; линейность операторов А и L также необязательна.
Доказательство теоремы существенно опирается на следующую лемму.
Л е м м а 3. Если операторы А и L линейны, а множество D выпукло, то для совокупной слабой замкнутости А и L на D необ ходима и достаточна их совокупная замкнутость.
Д о к а з а т е л ь с т в о . Если операторы А иЬ совокупно слабо •замкнуты на D, то они, очевидно, будут и совокупно замкнуты, так как из сильной сходимости элементов всегда следует их слабая сходимость.
Пусть операторы А и L совокупно замкнуты. Покажем, что тог да они и совокупно слабо замкнуты на D. Действительно, пусть
выполнены соотношения (2). Обозначим и„m = ип+ m |
(m = 0, 1,...). |
|||||
Очевидно, |
|
|
|
|
|
|
сл |
|
» |
|
|
|
|
чпп>— "«о. |
т-*оо, |
|
|
|
||
для любого п = 1 , 2 , . . . . |
Согласно известной теореме Банаха - |
|||||
Сакса существуют выпуклые комбинации вида |
|
|
||||
1 |
k |
ип т . |
ms = ms(n)-><*>, |
. |
|
|
v„k = — |
2 |
|
||||
k |
s =l |
|
|
|
|
|
такие, что |
|
|
|
|
|
|
lim \\v„k - |
u0||W = 0, |
/7= 1, 2, . . . |
|
|
||
к —•°° |
|
|
|
|
|
|
Выбираем |
к = к(п) |
00 (п-*°°) таким образом, чтобы |
|
|||
vn =vnk{n)-»u0, |
п |
|
|
|
||
По условию множество D выпукло, поэтому элементы |
а также |
|||||
v„ содержатся в D по построению. |
|
|
23
Покажем, что выполняются соотношения
_ сл |
_ сл |
|
(3) |
Avn — +fo* |
L»n— |
«^°° - |
Пусть z —любой элемент из F. Используя линейность оператора А, имеем
{A vn —/о, z)^ - |
1 |
|
к |
(Аип +nis ~ /о» |
9 |
к к(п). |
|
|||||
у- |
|
X |
|
|||||||||
|
|
|
|
|
к |
s = i |
|
|
|
|
|
|
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
I (Avn - / 0,z)F |
| < |
sup |
| (А щ - f 0,z)F |-*0, |
n-* «>, |
|
|||||||
|
|
сл |
|
|
|
i > n |
|
|
|
|
|
|
|
- |
|
|
|
°°. Аналогично доказывается справедливость |
|||||||
т.е. Avn |
— ►/o, n |
|||||||||||
второго из соотношений (3). Утверждение (3) доказано. |
|
|||||||||||
Далее, положим f n = |
и построим элементы /„ аналогично |
|||||||||||
j7„. Пусть для простоты |
|
|
|
|
|
|||||||
_ |
_ . |
_ |
1 |
к |
f nт , |
|
|
|
|
|
||
fn ~fnk |
“ |
2 |
|
|
|
|
|
|||||
|
|
|
к |
5 = 1 |
|
5 |
|
|
|
|
|
|
где снова |
|
|
s -*°°, |
к - к{п) |
п -*°°. Имеем |
|
||||||
— |
|
л |
|
|
1 |
£ |
|
|
|
|
|
|
/ п =Ли„, |
|
и„ = — |
2 |
* |
|
|
|
|
||||
Тогда Avn -*f0, л |
A: |
|
i - i |
|
|
|
и0, и -^°°. |
|||||
|
|
по построению. Покажем, что |
||||||||||
Действительно, |
|
|
|
|
|
|
|
|
||||
II |
- |
м0 \\н |
< |
1 |
А |
|
Wo llw |
|
|
|
||
|
2 |
II Vnnts - |
|
|
|
|||||||
|
|
|
|
|
/: 5=1 |
|
|
|
|
|
< sup II Vi - Wo II//-^0, i> n
Аналогично предыдущему доказывается соотношение
. л «л
Lvn — *•£<), |
|
|
Таким образом, построены элементы |
такие, что |
|
v„-*u0, |
Av„ -> f 0, Lvn —+go, |
(4) |
Проводя очевидные дополнительные построения и не вводя новых обозначений, будем считать также, что
Lv„-+g0, п-+°°. |
(5) |
Используя условие совокупной замкнутости операторов А и L на D, из соотношений (4), (5) заключаем: и0 €Z), AUQ = / 0, L U0 = = g0. Это означает, что операторы А и L также и слабо совокупно замкнуты.
24
Д о к а з а т е л ь с т в о |
т е о р е м ы 3. |
Используя экстремаль |
|
ное свойство регуляризованных решений иа , получаем |
|||
Фсе[ма ]< Ф (У[м], |
|
|
(6) |
где ы —решение основной задачи. Так как |
йа € D, то из экстре- |
||
мального свойства решения основной задачи следует |
|||
\ \A u - f\\F <\\Aua - f \ \ F. |
|
(7) |
|
Из соотношений (6), (7) следует |
|
|
|
\\Аиа - / \ \ р < Ц д + а 1>\, |
II Lua - g ||G |
< |
|| Lu - g ||G = vL (8) |
для любого OL> 0.
Таким образом, семейства Аиа , Lua при 0 < а < а < «> ограни чены и в силу гильбертовости пространств F и G слабо компактны. Используя условие дополнительности, убеждаемся в том, что се
мейство |
иа также слабо компактно. |
|
|||
Пусть |
ап > 0 —любая последовательность, сходящаяся к нулю |
||||
при |
п -+00. Без ограничения общности можно считать, что последо |
||||
вательности йп = |
f n = ЛиП9 |
gn = Lun (я = 1 , 2 , . . . ) |
слабо |
||
сходятся: |
|
|
|
||
л |
СЛ |
А СЛ |
А |
СЛ |
(9) |
ип— >и0у Аип — ►/<), Ьип — >g0i л ^ ° ° . |
|||||
Используя совокупную слабую замкнутость операторов А и L |
|||||
на D (лемма 3), из этих соотношений получаем |
|
||||
м0 € А |
Аи0 =/<>> |
Lu0 =g0. |
(10) |
Так как норма в гильбертовом пространстве слабо полунепрерывна снизу, то из соотношений (8) —( 10) получаем
M w 0 - / I I F < |
]m _ \\A un -f\\F < |
Urn |
\\Аип - f \ \ F < nA, (11) |
||
|
П - * шсо |
|
n -► |
oo |
|
\\Lu0 - g \ \ G < |
lim |
II Lun - g ||G < |
Шп |
\\LuH- g \\o < VL .(12) |
|
|
П «> |
П |
00 |
||
Из ( 11) следует, что |
ы0 € Uf, а из (12) |
- |
что и0 является реше |
нием основной задачи и, следовательно, совпадает с м. Тогда, оче видно,
S/4'M - / I I f = |
lim |
lAun - f i F = цА , |
|
n~+~ |
(13) |
IILu - g llG = |
lim |
i Lun - g lG = vL. |
|
n-+°° |
|
Замечая, что в гильбертовом пространстве из слабой сходимости элементов и сходимости норм следует сильная сходимость, из (9) и (13) устанавливаем, что | ип - и I -*0 (и-*00). Из условия допол-
25
нительности следует также, что йп -+и в Я . Так как последователь ность {<*„} О выбиралась произвольной, то аналогичные соотно шения справедливы и для всего семейства регуляризованных ре
шений. |
1. Если ц А = 0, т.е. основная задача совместна, |
С л е д с т в и е |
|
то из (8) следует, что |
|
\Аиа - f \ \ F < |
\ / а vL -►0 при а - * 0. |
Если оператор А к тому же имеет ограниченный обратный, то спра ведлива оценка
||ы<» -ы ||я < y / a v ^ A ' 1И. |
|
|
||
С л е д с т в и е |
2. Если одновременно |
рА = vL |
= 0, го регуля- |
|
ризованные решения иа =и, сс> 0. |
|
|
||
§ 3. Вариационное неравенство Эйлера. Оценки точности |
||||
1. |
Скорость |
сходимости к нулю |
А |
А |
уклонения |
| иа - и | может |
быть как угодно медленной без дополнительных предположений. Поэтому целесообразно выделение таких случаев, когда можно гарантировать определенную скорость сходимости. Прежде чем переходить к решению поставленной задачи, выведем одно важ ное свойство регуляризованных решений (р.р.).
Т е о р е м а 4. Чтобы элемент ua GD был решением регуляризованной задачи, необходимо и достаточно, чтобы выполнялось вариационное неравенство Эйлера:
(Aua - f , |
A(v - |
ua))F + a(Lua - g , L ( v - u a))G > 0 |
Vu€£>. ( 1) |
|
Н е о б х о д и м о с т ь . |
Положим для краткости |
|
||
| и1а= ||.4и||£ + |
||1 и ||д, U&D, |
|
||
и пусть 0 < |
в < |
1. Используя экстремальное свойство иа, имеем |
||
Фа [ма ] < |
Ф*[(1 - в ) и а +ви] Vu€£>. |
(2) |
||
Так как |
|
|
|
|
Ф«[(1 - в ) и а + 0и] = |
|
|
||
= Ф<* [«а] |
+ 26{(Аиа |
A(v - u a))F +ot(Lua - g ,L ( v - ua))G} + |
||
+ 02| и - |
иа \а, |
|
|
|
то неравенство (2) равносильно следующему: |
|
|||
(Айх - }\ |
A(v - |
ua))F + ot(Lua - g, L(v - ua))c + |
|
|
+ в \ v - ua l l > 0 Vu€Z). |
|
|||
Переходя здесь к пределу при 0 -*0, получаем (1). |
|
26
Д о с т а т о ч н о с т ь . Пусть элемент za G D таков, что верно вариационное неравенство
(Aza - /, A(w> - za))F + a(Lza - g, L(vv - zc,))G > 0 Vw G Z>.
(3) Полагая в ( 1) v = za, а в (3) vv = ua, и складывая эти неравенства, получаем | иа — z&| < 0, т.е. Аиа = Лга , Lua = Zza . Используя условие дополнительности, получаем также, что йа = za.
2.Пусть Я - линейное множество. Полагая в (1) v = иа ± w,
видим, что элемент |
характеризуется соотношением |
|
(Аиа - /; Aw)F + |
- g, Lw)G = 0 V w e Z>, |
(4) |
которое называется тождеством Эйлера.
Рассмотрим случай, когда D = DA П DL = DAL . Предполагаем, что DAL = Я , где черта означает операцию замыкания. Это условие обеспечивает существование сопряженных операторов А* и Z*.
Пусть на D ограничен один из операторов А или L , например А.
Тогда очевидно, что функционал |
(Lua |
- g, |
Lw )G ограничен |
||
на D и, следовательно, |
Lua - g Е |
*, |
т.е. представим |
в виде |
|
(L*(Lua - g ) , w ) H Vw |
G D. Тогда тождество Эйлера (4) |
прини |
|||
мает вид |
|
|
|
|
|
(А*Аиа - A * f t w)H +a(L*(Lua - g ) , w)H = 0 |
VwGD. |
|
Так как по условию D = Я , то это соотношение равносильно сле дующему уравнению Эйлера:
Л * Л + aZ * (1 ^ - g) = А */. |
(5) |
В частности, если .4 = Я, т.е. рассматривается задача вычисления значений оператора L , уравнение Эйлера принимает вид
ua +OLL\Lua - g )= f . |
(6) |
Если элемент g достаточно гладкий, именно g € Z)^*, то урав нение (6) принимает вид
(£ + aZ*Z)wa = /+<*Z*g. |
(7) |
Аналогично в случае, когда оператор L ограничен, для определе ния р.р. иа уравнение Эйлера принимает вид
^ 4 ^ |
- / ) |
+ aZ*Zwa =aZ*g |
(8) |
и, если |
L = Я, т.е. рассматривается задача решения операторного |
||
уравнения A u - j \ уравнение (8) записывается как |
|
||
А *(Аиа - / ) |
+ аиа = ag. |
(9) |
27
Если |
то уравнение Эйлера принимает вид |
|
(аЕ + А*А) йа= А 7 + <*£. |
(10) |
Отметим, что уравнения (5) —(10) однозначно разрешимы при любом а > 0. Это следует из их эквивалентности соответствую щим тождествам Эйлера.
Следующее утверждение характеризует решение и основной задачи.
Т е о р е м а 5. Пусть и |
- решение асновной задачи; тогда |
||
справедливы неравенства Эйлера |
|||
( A u - f , |
A(V - U))f |
> 0 |
VyE£>, |
(L u -g, |
L(w - U ))c > 0 |
П 1) |
|
VwEUf. |
|||
Если D - линейное множество, то справедливы тождества Эйлера |
|||
(Аи - / , A V)F = 0 |
VuED, |
||
(Lu -g,L w )G = 0 |
V w e U f . |
||
Д о к а з а т е л ь с т в о . |
Достаточно перейти в соотношении (1) |
кпределу при а -^0 и воспользоваться теоремой 3.
3.В § 2 установлена скорость сходимости значений операто
ра А на р.р. иа к Аи при условии совместности основной задачи. Рассмотрим общий случай.
Пусть йр Е D - р.р. при значении параметра /3 > 0. В силу теоре мы 4 имеем
(Айр - / , A ( w - йр))р+ |
|
|
|
|||
+ P(Lup-g, L ( w -йр))о > 0 |
|
Vw eD . |
(12) |
|||
Полагая v = йр в |
(1), w = йа в (12) и складывая получившиеся |
|||||
неравенства, имеем |
|
|
|
|
|
|
|Ы(3- Mala < (a - |
0) (Lup - g, |
Lup - Lua)G. |
|
|||
Переходя здесь к пределу при /3 |
0 и применяя неравенство Ко |
|||||
ши - Буняковского, получим |
|
|
|
|||
\иа - и\а<: all Lu - #IIG II Lu - |
Lua \\G = avL II Lua - Lu\\G. |
(13) |
||||
Следствием (13) является |
|
|
|
|||
Т е о р е м а 6. Пусть й - решение основной задачи, иа - |
р.р. |
|||||
Тогда справедливы оценки |
|
|
|
|||
II Lu - Lua llG < |
vL , |
\\Au - |
Aua IIF < у/a. vL . |
(14) |
||
С л е д с т в и е . |
Если |
vL = 0, |
то ua =uf ot> 0. |
|
Этот факт можно доказать также и непосредственно. В самом деле, из экстремальных свойств решения основной задачи при
28
vL = О следуют неравенства |
|
|
|
|||
Ф* [иа] <11 A u - f |
I* + all L u - g\\2G = д2 < Jl Aua - f l l F>2 |
|||||
из которых вытекает |
|
|
|
|
||
HAua - f l l F = llA u - fl l F = pA , |
II Lua - g llG = 0, |
a > 0 , |
||||
и, следовательно, |
в |
силу теоремы |
1 о |
единственности решения |
||
основной задачи |
= |
что и утверждалось. |
|
|||
4. |
В настоящем пункте будем предполагать, что множество D — |
|||||
линейное и плотное в Я . Тогда справедливо тождество Эйлера (4). |
||||||
Полагая в нем иа |
= и + za , где и - решение основной задачи, и |
|||||
воспользовавшись теоремой 5, имеем |
|
|
||||
(AZOL,A W)F + a (Lza + Lu - g t Lw)G = 0 |
\fwED. |
(15) |
Следующие теоремы устанавливают скорость сходимости регуляризованных решений йа к решению основной задачи и, если по следнее обладает определенными свойствами гладкости.
Т е о р е м а 7. Пусть оператор А ограничен, Lu - |
g Е D * и |
существуетхотя бы один элементу Е F такой, что |
JL |
|
|
L*(Lu - g ) = A*y. |
(16) |
Тогда справедливы следующие оценки : |
|
IILua - I MIIG < \/a"ll>'llir, |
|
IUMa - ^ t/llF < |
(17) |
Д о к а з а т е л ь с т в о . Из тождества (15) при условии (16) получаем
(Aza + ayуAw)F + a (Lzat Lw)G = 0 Vw ED.
Отсюда согласно теореме 4 следует, что элемент za доставляет минимальное значение функционалу
II Az + a>Hl£ + all/,zllG, zED.
Но тогда, выбирая в качестве элемента сравнения z = 0, получаем
llAza +ayl\F + allZza llG < a2\\ylF.
Отсюда следуют соотношения (17). Теорема доказана. |
|
З а м е ч а н и е . Оценки (17) имеют место также |
при неогра |
ниченном операторе А , если только оператор L ограничен. |
|
Т е о р е м а 8. Пусть оператор А ограничен, Lu - |
g Е DL* и |
существуетхотя бы один элемент h E D L такой, что |
|
L*(Lu - g)-A*Ah. |
(18) |
Тогда справедливы следующие оценки:
\\Aua -Au\\F <
< alU/?llF + a 3/2 IIZ./HIG, |
\\Lua -Lu\\G < allZ7j!G. |
(19) |
Д о к а з а т е л ь с т в о . Из |
тождества (15) и условия |
(18) |
получаем |
|
|
(Aza + aAh, A W )F + a(Lzat Lw)G = 0 V w E D.
Согласно теореме 4 отсюда следует, что элемент za доставляет ми нимальное значение функционалу
II Az + aAhllF +allLzll^, z ED.
Выбирая в качестве элемента сравнения z = - a h , получаем
II Aza +aAh\\F +all Lza II2G < a 3II LhII2G.
Отсюда следуют требуемые оценки (19). Теорема доказана.
З а м е ч а н и е . Оценки (19) |
имеют место также и при неогра |
ниченном операторе А, если только оператор L ограничен. |
|
Пусть L = Е. Тогда условия (16) и (18) принимают ътусловий |
|
истокообразности решения основной задачи: |
|
u - g = A*y |
(20) |
или |
|
и - g = A*Ah. |
(21) |
Пусть А = Е. Тогда условия |
(16) и (18) принимают вид усло |
вий гладкости решения основной задачи: |
|
L u - g E D L *, |
(22) |
L * (L u - g )E D L . |
(23) |
Отметим, что условия (16), (18) заведомо выполнены, если вы
брать |
g = Lu. Тогда условия |
(16), |
(18) удовлетворяются при |
у = 0, |
h = 0. Легко видеть, что в этом случае |
||
LUu = Lu, Аиа ^А и, иа =и |
при |
a > 0. |
Достаточные условия, при которых имеют место соотношения (20), (21), приведены в главе 5.
§ 4. Устойчивость регуляризованных решений
1. Мы рассматриваем регуляризованное решение как функцию элементов /,# , операторов A, L , а также множества D.
Наименее трудно дать оценку устойчивости р.р. при возмущении элементов f u g . Будем считать, что вместо элементов / и g за-
30