Об оптимальности метода Тихонова нулевого порядка на некоторых классах корректности

УДК 517.948

ОБ ОПТИМАЛЬНОСТИ МЕТОДА ТИХОНОВА НУЛЕВОГО ПОРЯДКА НА НЕКОТОРЫХ КЛАССАХ КОРРЕКТНОСТИ

Н.М. Япарова

В работе проведено исследование величины погрешности метода регуляризации Тихонова нулевого порядка. Доказана оптимальность по порядку этого метода на некоторых классах корректности.

Ключевые слова: метод решения операторных уравнений, метод регуляризации, оптимальность по порядку, модуль непрерывности, погрешность метода.

При решении практических задач большое внимание уделяется выбору оптимального метода. В том случае, когда построение оптимального метода является трудновыполнимой задачей, вместо оптимального метода используют метод, имеющий тот же порядок погрешности, что и оптимальный, то есть применяют метод, оптимальный по порядку. Оптимальные по порядку методы имеют широкое применение при решении обратных задач математической физики неустойчивых относительно возмущения исходных данных.

Одним из распространенных методов решения подобного рода задач является метод регуляризации Тихонова нулевого порядка [1]. При численном решении обратной задачи физики твердого тела [2] было замечено, что этот метод позволяет получать более точные результаты, связанные с «тонкой» структурой решения, чем другие методы регуляризации. Это послужило поводом для исследования данного метода на оптимальность.

В работе установлено, что одним из факторов, влияющих на оптимальность метода, является структура множества, которому принадлежит искомое решение. Метод регуляризации Тихонова нулевого порядка был исследован на оптимальность по порядку на различных классах корректности и для этого метода доказана оптимальность по порядку на некоторых из рассматриваемых классов.

1. Постановка задачи

Пусть Н - гильбертово пространство, а А : Н  Н - инъективный линейный ограниченный оператор. Рассмотрим операторное уравнение первого рода

Аи = /, и,/ е Н.

Предполагается, что при / = /0 существует точное решение и0 е М , но вместо /0 дано некоторое приближение /5 е Н и 5 > 0 такие, что \/5 — /0) || < 5 .

Требуется по исходным данным задачи М, /5 и 5 определить приближенное решение и5 уравнения (1) и оценить его уклонение от точного решения.

Опишем множество М. Пусть оператор А1 определен соотношениями: А1 = А*А, спектр 5р(А1) = 0,||А||2 . Рассмотрим множество 8г = {у: V е Н, ||у|| < г} и линейный, ограниченный оператор В : Н  Н такой, что для оператора В1 = ББ* выполнено соотношение

В1/2 = с

где функция С(а) - непрерывна, строго возрастает на [0, ||А||J и С(0) = 0 . Тогда множество М есть множество Мг = ВБг.

Определение 1. Семейство непрерывных операторов {Т5 :0< 5 < 60} будем называть методом приближенного решения уравнения (1) на Мг, если "5 е (0,60 ], оператор Т5 : Н  Н и

Т5 /5  и0 при 5  0 равномерно на множестве Мг при условии, что \/5 — Аи01 < 5 .

Для дальнейшего изложения нам потребуется понятие и свойства модуля непрерывности а(т, Мг) обратного оператора в нуле [4]

Если а(т,г)  0 при т  0 , то Мг является классом корректности для уравнения (1).

Определение 2. Погрешностью метода {Т5 :0<5 <50} на классе корректности Мг назовем

функцию Л(Т5 ) такую, что "5 е (0,50 ]:

Л

и,/5

Определение 3. Метод {Тдр{: 0 < 5 < 50} оптимален на классе корректности Мг, если "5 е (0,50 ]

л

где Л5 (Т) определено формулой (3).

Теорема 1. Пусть Мг - класс корректности для уравнения (1), а Л Ор1 определена формулой

(4), тогда справедливы неравенства

2 а(5 , г) <Л°р1 < а(5 , г).

Доказательство приведено в [5].

Определение 4. Метод { Т5 :0 < 5 < 50} оптимален по порядку на множестве Мг, если Зк

такое, что " 5 е(0, 50 ] выполнено:

Л(Т) < ка ( 5 , г).

2. Метод регуляризации Тихонова нулевого порядка

Следуя [1], метод регуляризации Тихонова нулевого порядка заключается в сведении задачи приближенного решения уравнения (1) к вариационной задаче

{ Аи — / + а|и : и е Н}, а> 0.

и

т

(5), определяемое формулой

Из работы [7] следует, что для любой / е Н существует единственное решение иа задачи

u5 = Raf5 ,

где

Ra =

Для оператора Ra , определенного формулой (7), имеет место следующее соотношение [8]:

IRa II = ^= •

Va

Рассмотрим следующую величину:

D1

u0

В качестве параметра регуляризации a = a( 5 ) возьмем решение уравнения

Di

Пусть семейство операторов {Ra :0 < a < a0} определено формулой (7).Уравнение (10) будет

иметь единственное решение a(5) [7]. Оценим погрешность метода {R„(5 ):0< 5 <50} при вы-

боре параметра a , удовлетворяющего уравнению (10), на различных классах корректности.

3. Исследование точности метода регуляризации Тихонова нулевого порядка

Покажем, что этот метод не всегда будет оптимальным по порядку. Рассмотрим величину

Л0

а>0

где Л5 (Яа ) - функция, определенная формулой (3).

Из результатов, сформулированных в [9, с. 44-46], имеем следующее утверждение Замечание 1. Пусть а(5 ) - решение уравнения (10), тогда "5 е (0, 50 ] и 50 < гС(||а||)-||а|| , имеет место оценка

Л\

В качестве класса корректности рассмотрим множество Мг = ВБг, где оператор В определен соотношением (2), а

С

Соответствующий класс корректности Мг будем называть степенным классом.

Теорема 2. Если р е (0,2], то метод {Ла(- 5) :0 <5 < 50} оптимален по порядку на множестве

Мг .

Доказательство. Так как множества значений Я(В) и Я(В*) операторов В и В* всюду плотны в Н , то из леммы, сформулированной в [8], следует существование полярного разложения операторов В

В = QB11/2,

где Q - унитарный оператор.

Подставляя в Л2(а) представление (14), получаем

2

A,2 (a) = max i

ни r 0

a

2

a + s

G

где а0 = Q Ч”0, а {Еа : а е [0,||А||]} - разложение единицы Е , порожденное оператором А11/2 . Из (15) и теоремы, сформулированной в [1], следует, что

A1

I a + s2

Пусть p < 2, тогда из формул (16) и (13) следует, что

A1(a) = r max ----------------- G(s) = r max ------------------ sp.

r~ la + s2 se[0,|A|| ] a + s2

J(s) =--------2 sp .

Введем функцию

a a + s

Эта функция неотрицательная и непрерывная по s на отрезке [0,||A||] . Тогда существует значение s = 6(a) такое, что

у (s )= max у (ст) .

Заметим, что функция y(s) дифференцируема, и в точке локального экстремума должно быть выполнено соотношение

sp-1 (pa + ps2 -2s2) = 0.

Отсюда следует, что либо s = 0, но тогда y(s) = 0 , либо

s 2 = - pa

2 - p

Последнее соотношение будет иметь место при р < 2 и в точке

Ра

'2 - р

функция у (а) достигает максимума.

Таким образом, из (16) следует, что на степенном классе корректности

2 - Р

а

Из (8), (10) и (17) следует, что

а( 8 ) =

2 - Р

р

2

р+11 8 | р+1 г

На основании (12), (17) и (18) имеем, что для любого 8 є (0, 80 ]

• гР+1 8 Р+1

А 5

<

2 - Р

2

Заметим [9], что в случае функции С(а), определяемой формулой (13), имеет место соотношение

1р

Л°р1 = гр+1 5 р+1, где р е (0, да), а ЛУ определена формулой (4).

Из (19) и (20) следует оптимальность по порядку метода Ла(- д) для р е (0,2). Пусть р = 2 , тогда

2

ш\\А Лх(а) = 2

1411 + а

и при достаточно малых значениях 5 из (21) следует, что

*(5) ‘111' ( г 1'

1 1 2

Из (22) следует, что в этом случае

А*

Из (20) и (23) следует оптимальность по порядку этого метода.

Теорема 3. Если р є (2,да), то метод {^а(8 ): 0 < 8 < 80} неоптимален по порядку для Мг. Доказательство. Так как в этом случае

Из (24) следует, что при а < 11^|| а из (25) следует, что

А1(а) = 7^2-------

Ы + а

2

2

Из (26) имеем, что при достаточно малых значения 8

2

1 2

/ \ 1 и ц2“Р - -

А5(ящ5))>2ІИ 3 • г3 83.

Из (20) и (21) следует неоптимальность по порядку метода {л— 3) :0 < 8 < 80} на классе Мг. Тем самым, теорема доказана.

Теперь перейдем к исследованию поведения метода {л— 5) :0<8 <80} на других классах корректности. Пусть оператор В определен соотношением (2), а

и Мг = ВБг, то класс корректности Мг будет называться логарифмическим классом.

Оценим погрешность метода {R*(g ) :0 <5 < 50}, заданного формулами (7) и (10), на логарифмическом классе.

Из формулы (27) и свойств модуля непрерывности обратного оператора [4] следует сущест-

вование положительных чисел /і и ¡2 таких, что

l1 • r ln-q \ — |< w (т, r) < /2 • r ln-q I — I .

Т 0 ^ Т

Оценим А1(а).

Теорема 4. Существуют положительные числа І3 и І4 такие, что при достаточно малых значениях а

І3 • г 1п~ч ^—0 <А1(а) < І4 • г 1п-д ^^.

Доказательство. Из формул (2), (27), (7) и (28) следует существование положительных чисел І5 и І6 таких, что

l5r sup

a

ln q | — j<A1(a) < l6r sup _a + s' ~ "

Сначала получим оценку снизу. Из (16) имеем $/ > 0 :

a

a

a + s

ln - qi.

A1(a) > lr-

a + s

G (s).

А так как С (а) - возрастающая и при а  0 выполнено (27), то при а е [0,|| А|| ] имеет сто оценка

ме-

a

-ln-qI - I.

A1(a) > lr-

a + s V s

Рассмотрим два случая. Пусть a = max (1,|| A||) = 1. Тогда

ln - q fa 1>_L 1„ 4 f i a + s V s 0 2 q V a

во втором случае a = max (1,||A||) = ||A||. Тогда будет верным соотношение

a

- q I a

1, I 1

ln A + - ln I -2 l a

а + а ( а Из (30), (31), (32) имеем оценку снизу

Л1(а) ‘115 • г 1п~ч — . 2а

Теперь перейдем к оценке сверху. Введем функцию

>-^rn-q I і

2-q+1 I a

^Т 1п а [ — I, при а > 0,

У

0, при а = 0.

Функция у (а) непрерывна и неотрицательна по а на отрезке [0,|| А|| ] . Тогда существует значения а = а (а ) такое, что

у (у) = вир у (а).

Точка а = у*) либо совпадает с одним из концов отрезка [0,||А||] , либо совпадает с точкой локального максимума.

Предположим, что у (а) достигает наибольшего значения в точке локального максимума а . Оценим значение у (а) в точке а .

Так как у (а) дифференцируема при любом а е ^ 0,|| А|| ] , то в точке локального максимума а = а должно выполняться

2а

а + а -, а а1п —

Вторая производная от функции у (а ) в точке 7 имеет вид

у"(а ) =

4 у

(а + (7 21

-(а2 — а) .

Возможны два случая. В первом случае при а >л[* эта производная положительна и, следовательно, при а >л/* функция у (а) не может достигать локального максимума. Тогда точка локального максимума должна удовлетворять следующим условиям:

(7 <4ш 2(7 а

а + су -, а

су 1п —

а

I— а 11

Но при а < Vа имеем — > — >—= и, следовательно,

а а у/а

.а 1. 1

1п — >— 1п—. а 2 а

Отсюда следует, что при а <у* будет справедлива оценка

1п—а — < 2а 1п—а -,

7 а

а так как при любом значении а имеет место неравенство

а

-----^< 1,

а + а

то при а <у/* имеет место следующая оценка:

а

Л1 (а)< 13г вир

а + а

1п—а [ — |< 13г1п—а — < 2а13г 1п—а I.

а + а2

7 0 а + су а (а

Таким образом, в точке локального максимума будет выполнено требуемое неравенство. Исследуем поведение функции у (а), заданную формулой (33) на концах отрезка [0,|| А|| ] .

При а  0 имеем у (а)  0 . Исследуем случай, когда а = а = | |А||. Для этого рассмотрим величину

у И

а

а

а (а +ИІ2 )іп’-а

Покажем, что это - ограниченная величина. Предположим противное. Однако из того, что а  0 следует, что ЗК такое, что

|а +1ІЛІІ2) 1п9 А = К .

Отсюда при а  0 имеем

а

(а +|\Л\\2 ) 1п

0,

что противоречит (36). Таким образом, из (37) следует, что существует число /5 такое, что

у (о) =

а

_2 а + о

- -‘5' О

а

Из (29), (35), (38) следует утверждение леммы.

Теорема 5. Пусть а и г > 1, а метод ) :0 < 8 < 80} определен формулами (7) и (10). Тогда этот метод оптимален по порядку на классе Мг .

Доказательство. Так как и0 є Мг, то

и% и

<А1 (а) + |\Ла\ | • 8 .

Из теоремы 4 и формул (8), (39) следует, что если и0 є Мг, то

1

„  -9I 1 1 1 -

и 8 и.

4 8 и0

При а (8 ) = 8 из (40) следует, что

< І4 • г 1п ч \ — | + — а 2 • 8 . а 0 2

а

иг

. 8 0 2

Из (41) следует существование 81 > 0 такого, что 81 < 80 и для любого 8 < 81

а

ир

где І7 - некоторая константа.

Если а и г ‘ 1, то из (42) следует, что З 52 > 0 такое, что 52 < 5! и для любого 5 < 52

а

иг

< 2« 17 • г 1п“«| аг I.

Из оценки снизу, приведенной в теореме 1, формул (28) и (43) следует оптимальность по порядку метода {л*5) :0 <5 < 52} , определяемого (7) при а (5) = 5 на классе решений Мг.

Из (12) и (43) следует оптимальность по порядку на логарифмическом классе Мг метода 5 ) :0 <5 < 50}, определяемого формулами (7) и (10).

Заключение. В работе установлена оптимальность по порядку метода регуляризации Тихонова в зависимости от структуры множеств решения уравнения (1).

Литература

1. Тихонов, А.И. О регуляризации некорректно поставленных задач / А.И. Тихонов // Докл. АН СССР. - 1963. - Т. 153, № 1. - С. 49-52.

1

1

2. Лифшиц, И.М. Об определении энергетического спектра Бозе-системы по ее теплоемкости / И.М. Лифшиц // ЖЭТФ. - 1954. - Т. 26, № 5. - С. 551-556.

3. Коршунов, В. А. Определение плотности состояний по термодинамическим функциям / В.А. Коршунов, В.П. Танана // ФММ. - 1979. - Т. 48, № 5. - С. 908-915.

4. Иванов, В.К. О линейных некорректных задачах / В.К. Иванов // Докл. АН СССР. - 1962. -Т. 145, № 2. - С. 270-272.

5. Иванов, В.К. Теория линейных некорректных задач и ее приложения / В.К. Иванов, В.В. Васин, В.П. Танана. - М.: Наука, 1978. - 208 с.

6. Иванов, В.К. Об оценке погрешности при решении некорректных задач / В.К. Иванов, Т. И. Королюк // Журн. вычисл. матем. и матем. физ. - 1969. - Т. 9, № 1. - С. 30-41.

7. Морозов, В.А. Оптимальная регуляризация некорректных нормально разрешимых операторных уравнений / В.А. Морозов, С.Ф. Гилязов. - М.: Изд-во МГУ, 1982. - 270 с.

8. Танана, В.П. Оптимизация методов решения операторных уравнений / В.П. Танана, М.А. Рекант, С.И. Янченко. - Свердловск.: Уральск. ун-т, 1987. - 200 с.

9. Танана, В.П. Методы решения операторных уравнений / В.П. Танана. - М.: Наука, 1981. -180 с.

10. Рудин, У. Функциональный анализ / У. Рудин. - М.: Мир, 1975. - 448 с.

Поступила в редакцию 18 мая 2007 г

ON ACCURACY OF TIKHONOV REGULARIZATION METHOD OF ZEROTH ORDER AT SOME CLASSES OF CORRECTNESS

The article presents the research of accuracy of Tikhonov regularization method of zeroth order. The optimality according to the method at some classes of correctness is proved.

Keywords: solution method of operator equations, regularization method, module of continuity, error of method.

Yaparova Natalia Mikhailovna - Candidate of Science (Physics and Mathematics), associate professor of the Calculus Mathematics department, the Mechanical-Technological Faculty, South Ural State University.

Япарова Наталья Михайловна - кандидат физико-математических наук, доцент, кафедра вычислительной математики, механико-математический факультет, Южно-Уральский государственный университет.

http: //www .math.susu.ac.ru