Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Регулярные методы решения некорректно поставленных задач

..pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
11.88 Mб
Скачать

5. Остановимся на некоторых примерах, приводящих к решению основной задачи.

П р и м е р 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^15—) 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