О проблеме потери точности при преобразовании информации Текст научной статьи по специальности «Математика»

УДК 517.948

О ПРОБЛЕМЕ ПОТЕРИ ТОЧНОСТИ ПРИ ПРЕОБРАЗОВАНИИ ИНФОРМАЦИИ

Н.Ю. Колесникова, Т.Н. Рудакова, A.B. Танана

THE PROBLEMS OF THE ACCURACY LOST DURING INFORMATION TRANSFER

N.Y. Kolesnikova, T.N. Rudakova, A.V. Tanana

Рассматривается метод M.M. Лаврентьева решения операторных уравнений в гильбертовых пространствах. Получены точные оценки погрешности, доказывающие его оптимальность. Этот метод применён к решению обратной задачи Коши для уравнения теплопроводности. Показано, что при преобразовании информации этим методом происходит минимальная потеря точности.

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

In this paper, we proved the optimality of M.M. Lavrentiev method for solving of the ill-posed problems. The precise estimations of the accuracy of this method were obtain. This method has been applied to the inverse problem of Cauchy for heat transfer equation. When we used this method the transfer information will be with minimum loss of accuracy.

Keywords: operator equations, regularization, the optimal method, error estimate, ill-posed problem.

Пусть ^ и и - некоторые гильбертовы пространства, а Т - неограниченный оператор, действующий из Р в и . Предположим, что нам требуется определить значение и0 -Т/0, но вместо /0 в качестве исходной информации нам известна пара (/5,5), в которой /5 е ^, а ||/5 - /01| < 8. Требу-ется эту информацию преобразовать в пару (иг, г), в которой иг е и, а \иг -щ|| < 8 .

Так как оператор Т неограничен, то без дополнительной информации о значении и0 г - оо , что совершенно неудовлетворительно. Поэтому в качестве дополнительной информации можем использовать задание класса корректности Мс[/, которому должно принадлежать точное значение и0. Примером такого класса корректности может служить множество М из и, удовлетворяющее условию

Г 6 с [г-1 (м)]

Колесникова Наталья Юрьевна - старший преподаватель кафедры вычислительной математики ЮУрГУ;

1 @таП.ги

Рудакова Татьяна Николаевна - канд. физ.-мат. наук, доцент кафедры вычислительной математики ЮУрГУ, rtn@susu.ac.ru

Танана Алексей Витальевич - канд. физ.-мат. наук, доцент кафедры математики Уральского горного университета, г Екатеринбург; tvpa@susu.ac.ru

Одна из основных проблем математической информатики заключается в разработке методов преобразования информации (/5,8) в (г/£,г) и определении нижней грани 8 среди всех возможных.

В настоящей работе методом М.М. Лаврентьева [1] решено операторное уравнение первого рода Аи = /5;ме*У,/5 е/%

в предположении, что |^4-11| = °о , а ||/§ - /01| < 5 .

Доказана оптимальность этого метода и найдена нижняя грань 8 среди всех возможных. Результаты работы проиллюстрированы на обратной задаче Коши для уравнения теплопроводности.

1. Понятие метода преобразования информации

Пусть и, Т7 и V - гильбертовы пространства, А - инъективный линейный ограниченный оператор, отображающий и в ^ и имеющий неограниченный обратный, В - линейный ограни-

Kolesnikova Natalia Yurievna - senior lecturer of Computational mathematics department of SUSU; Natasha720221@ mail.ru

Rudakova Tatiana Nikolaevna - PhD, associate professor of Computational mathematics department of SUSU; rtn@susu.ac.ru

Tanana Aleksei Vitalievich - PhD, associate professor of Mathematics department of Ural Mountain University, Yekaterinburg; tvpa@susu.ac.ru

ченный оператор, отображающий V в U, а Mr = 2?£г , где Sr = |v: v g V, ||v|| < -rJ . Предположим, что

а-^с[а(мг)\ (1)

Рассмотрим операторное уравнение первого рода Au=f (2)

и предположим, что при f = fо существует точное решение и0 уравнения (2), которое принадлежит множеству Мг, но вместо /о нам известны fb е F и 5 > 0 такие, что

||/5-/о1И

Требуется, используя априорную информацию Мг5 /8 и 8, определить приближенное решение иъ уравнения (2) наименее уклоняющееся от точного решения и0 И оценить уклонение ||и5 - и01 на классе Мг.

Определение 1 Оператор Р непрерывно отображающий пространство F в U будем называть методом приближенного решения уравнения (2) на множестве Мг.

Введем количественную характеристику точности метода Р на множестве Мг

Д8 (Р) = sup {!« - Pfb I: и <= Mr, IЛи - /8II < 8}. (3)

«,/б

Таким образом, метод Р преобразует исходную информацию (Мг, /ъ, 5) в информацию [Pf§, А§ (Р))

Мы заинтересованы в том, чтобы величина Д5 (Р) была минимальна. Для этого через C[F,U] обозначим множество операторов непрерывно отображающих пространство F в U и определим

число Ag^ формулой

Af = inf {Д5 (Р):Р<= С [F,I7]}, (4)

где А§ (Р) - определена формулой (3).

Определение 2. Метод Popt будем называть оптимальным на классе Мг, если

Д5 (рор1) = АГ

Теперь следуя [2], определим функции (т, г) и со(т, г) формулами: со1(т,г) =

= sup {||м, -и2\\-щ,и2емг, I Ащ -Аи2 ||<т} (5)

И

ш (т, г) = sup {\\и\\ иеМг, \\Аи\\ < т| Л (6)

где т > О

Из леммы доказанной в [3, с. 17] следует,

что

ю1(т,г) = ю(т,2г) (7)

Лемма 1. Пусть к > 0. Тогда справедливо равенство ф(&т, кг) = £со(т, г)

Доказательство. При к = 0 лемма очевидна. Пусть к > 0 и т > г \АВ\ Тогда кх > кг |/42?||, а из (6) следует, что

а>(т,г) = г||у4В| (8)

и

ф(£т, кг) = кг ||^#|. (9)

Из (8) и (9) следует, что ш(&т, &г) = £а>(т, г)

Пусть к> 0 и х<г||у4^||. Тогда из того, что

иеМг и ||у4и||<т следует, что киеМ^ и |у4(Ь/)| < кх

Таким образом

ксо(х, г) < ю(£т, кг) . (10)

В другую сторону. Пусть и е М& и \Аи\ < кх.

Тогда -еМг и \\А к

< X, т.е.

ую(£т, *г)<ю(т,г),

или

со(А:т, кг) < А:со(т, г) . (11)

Из(10)и(11) следует утверждение леммы.

Лемма 2. Пусто Ре С [/%£/], а (т, г) определена формулой (7). Тогда справедлива оценка

Дв(/>)*1Ю1(25,г).

Доказательство. Пусть 8 - достаточно малое положительное число. Тогда, на основании (7) существуют элементы щ,и2 е Мг такие, что

-и2\ > ©! (25, г)-8, (12)

а

\Ащ-Аи2\<2Ъ (13)

- (Ащ -Аи2 )

Если положить /§ =---------------, то из (13)

будет следовать, что

\\Аи1 - 7ъ | ^6 и \Аи1 - 7& | ^5 • (14>

Так как

тах {|и, - Р/51|, ||м2 - Р/51| > 1 у»! - м21, (15)

то из (12), (14) и (15) следует, что шахЩм, -.Р/5||,||и2 -а из (3), что

Д8 (Р) > тах {¡и, - Р/ъ ||; ||м2 - Р/51||. (17)

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

А5(Р)>1со1(2§,г)-|. (18)

Ввиду произвольности е из (18) следует утверждение леммы.

Из формул (4), (7) и лемм 1 и 2 следует, что

А^>ш(5 ,г). (19)

2. Исследование на оптимальность метода М.М. Лаврентьева

Пусть и = Р = У = Я , а операторы А* А и

ВВ* положительно определены. Тогда на основании леммы, доказанной в [4], для операторов А и В имеют место полярные разложения А = ()А и

В = ВР, где <2 и Р - унитарные операторы,

А = у]А* А , а В = ^¡ВВ* . Кроме того, предположим, что спектр 8р{А) оператора А совпадаете

отрезком [о, 1^1], а

В = о(а), (20)

где С/(а) - строго возрастающая, непрерывная на отрезке [о, ||^4||] и дифференцируемая на интервале (о, ||,4||) функция, такая, что <3(0) = 0 . Рассмотрим уравнение

гв(о)а = т; а е [о, ¡А\\\ (21)

Из (21) следует, что если т<0 и т < гС? (||^4||

то это уравнение имеет единственное решение

а(т) = \|/|^- , где \|/(х) - функция обратная

G(a)a. Из известной теоремы об обратной функции следует, что у еС^О, О(|^||)И|] и 1|/(0) = 0

Таким образом

а (т) 0 при т -> 0 . (22)

Лемма 3. Если выполнены все условия на операторы А и В, сформулированные выше, а т < г ||у4#|| , то справедлива формула

ш(т, г) = /*?[а(т)],

где а(т) решение уравнения (21).

Доказательство. Пусть 8 - достаточно малое, положительное число. Тогда, выбрав натуральное число п0 таким, что

«0-1-

К/ [а(т)]-К7

"о

-а(т)

<8,

(23)

рассмотрим подпространство Я0, определяемое формулой

Но =Ецх)Н-

^о~1-

Кт)

я,

(24)

где : 0 < а < ||^||| - разложение единицы, порожденное оператором А [5, с. 336]. Легко проверить, что

МГ=ВБГ (25)

Пусть у0 еН0 и

||у0|| = г. (26)

Тогда из (25) и (26) следует, что

и0 = Ву0 е Мг (27)

Так как и0 еЯ0, то на основании (23)-(26)

||м0||>га[а(т)]-е. (28)

Ввиду того, что и0,Аи0 еЯ0, а функция б(ст) строго возрастает, из (24) следует, что

1| < Ю [а (т)] а (т) = т . (29)

Из (27) и (29) следует, что

||и0||<ю(т,г), (30)

а из (28) и (30), что

ш (т, г) > г& [а (т)] - 8.

Ввиду произвольности 8 будем иметь

ю (т, г) > К? [а (т)] . (31)

Теперь докажем неравенство в другую сторону. Для этого представим пространство Я в виде ортогональной суммы

Н = Н1+Н2, (32)

подпространство Нх = Ец^Н, а Н2 = |я.

Из теоремы [5, с. 336] следует, что подпространства #! и #2 инвариантны для операторов А и В. Из того, что и0 еМг, а рфт, (33)

следует существование элемента такого, что

Ы<>- (34)

И

и0 = Ву0 (35)

Используя (32), представим элемент у0 в виде ортогональной суммы

У0=Г1+У2, (36)

где V,. =рг(у0,Н,),1 = 1,2.

Пусть ^ = ¡V! ||, а г2 = ||у2 ||. Тогда из (34) и (36) следует, что

г,2 +г22 <г2 (37)

Из инвариантности подпространств Нх и Н2 для

оператора В и формулы (35) следует, что и0 = и\ + и2 И

щ =Щ еЯ,;/ = 1,2. (38)

Из инвариантности подпространств Я1 и Я2 для

оператора А будем иметь

Аи1 еЯ;,/ = 1,2. (39)

Из (33), (38) и (39) следует, что

Г;

Ыг/у < —т;/ = 1,2. н и г

(40)

Так как С(<т) строго возрастает, то из (38) следует, что

ЩЦ ^ (т)],

а из (40), что

и? < -

(41)

(42)

г а(т)

Ввиду того, что

г20[а( т)]а(т) = -^-т, (43)

из (42) и (43) следует, что

||м2||<г2с[а(т)]. (44)

Из (37) и (38), (41) и (44) следует, что

||«о||<К7[а(т)]. (45)

Ввиду произвольности и0 из (33) - (35) и (45) сле-

дует, что

ю(т,г)<К?[ст(т)], (46)

а из (31) и (46), что

ш(т, г) = Кх [а(т)].

Из формулы (46) следует утверждение леммы.

Из леммы 3 и формулы (22) следует, что множество Мг является классом корректности для уравнения (2), а из формулы (19) и леммы 3 следует, что

ДГ>К7[<т(т)]. (47)

Так как, используя лемму работы [4], уравнение (2) можно заменить эквивалентным

Au = g, (48)

где А = л/а* А , g = $ f, а множество Мг , с учетом формулы (25) совпадает с множеством ВБг,

где В = ^ВВ* , то задачу приближенного решения уравнения (48) поставим следующим образом.

Предположим, что при g = g0 еН существует точное решение и0 уравнения (48), которое принадлежит множеству Мг, но точное значение правой части g0 нам не известно, а вместо него даны некоторое приближение g^^H и уровень погрешности 5 > 0 такие, что |^§ - #01| < 5 .

Требуется по исходным данным Мг, g8 и 5 определить приближенное решение иъ уравнения (48), минимально уклоняющееся от и0 на классе Мг и оценить это уклонение.

Будем решать эту задачу методом М.М. Лаврентьева [1], в котором используется регуляризующее семейство операторов . 0 < а < а0}, действующих из Н в Н и определяемых формулой

7^=2?[с + а£] 1,ае(0,а0], (49)

а С = АВ .

Лемма 4. Для любого а > 0 оператор , определяемый формулой (11) ограничен и О(ст)

о<а<|И G (а) а + а Доказательство. Так как

IKII2 =,ff IMI2,

(50)

|M2=(/U,/U), (51)

то ввиду самосопряженности оператора Ra из (50) и (51) следует, что

IKII2 =sup(*ag,g). (52)

II#1

Из (49) и (52) следует, что

IKf = supijB2[c + a£] 2g,gl. (53)

||*Л L J )

Пусть {Ест • 0 < ст < ЦлЩ - спектральное разложение единицы, порожденное оператором А. Тогда из (49) следует, что

2 ¥ G2 (а)

= / г...■;.У..п2 dEcg, (54)

0 |^aG(a) + aj

а из (53) и (54), что

II ||2 И оЧа)

IKII 1 г . (55)

Ы-1 0 [aG(cj) + a]

Учитывая (17), получим, что

IKf * sup

G»

°<a<||4 j^a(7 + a J ||g||<l 0

а из (56), что

sup ¡d(Ecg, g), (56)

sup

G2(a)

Так как функция

°-a-MI [aG (a) + a J G2(a)

(57)

непрерывна на

[aG(a) + a]

отрезке [О, И|], то существует значение а е [о, ||^||] такое, что

G2(a) _ G2 (a)

• = sup

(58)

[aG (a) + a] [aG (a) + a]

Из соотношений (57) и (58) следует утверждение леммы.

Лемма 5. Для любых а и г > 0 справедливо соотношение

II - - II

sup \\lLCv-Bv <ra шах ----------—•

Цфг 11 о<с<|л|| aG(o) + a

Доказательство. Так как

В [с + оЕ]“' Cv-Bv = -аВ [с + а£]-1 v, (59) то из (49) и (59) следует, что

IK Cv - jBv|| = а р [с + аЯ] v|| Если v ^ 0, то из (60) следует, что

||^Су-Яу|| = а|Н|

в\_С + аЕ]

-1 V

IM

(60)

(61)

W

.(62)

(63)

Так как

sup||i?aCv-i?v||= sup ¡iLCy-Z^vll,

о<М^

то из (61) следует, что

sup||iLCv-i?v|| = ra sup в\с + аЕ~\

IHN'' 11 ми L J

Из (49) и (62) следует, что

SUP IK^v - 5v|| < m I.

Ъ¥г

Из (63) и леммы 4 следует утверждение леммы. Пусть и0 е Mr, a ||gs - Аи01| < 8 Тогда

К - К gb I ^ ||мо - Ra Äu01 + 5 ¡7^ I. (64)

Так как

||и0 - Ra Äu01 < sup 1^ Cv0 - Bv01|, (65)

то из лемм 3, 4 и соотношений (64), (65) следует, что

<66)

Обозначим через G'(cj) производную от функции G(g), а через cj(S) решение уравнения (21) при г - 8.

Теорема 1. Пусть для любого сге(о, |

/ ч °2(а)

G'(a)>0, — j возрастает, 8<rG

а

(8) =

g2(^(S))

G'(ä(S))

,а

(67)

Тогда

А5 ( ^«(5) ) ^ ^ (§)] .

Доказательство. Пусть щ произвольный элемент из множества Мг, а ||^§ - Аи01| < 8 . Тогда из формулы (66) следует, что

К• <б8>

О(о)

Теперь вычислим шах<

aG(a) + a(S)

Для этого продифференцируем функцию G(cj)

aG(a) + a(8)

•0<a< И||

oG(ct) +сс(8) ä(8)G'(a)-G2 (a)

aG(a) + a(8)]2

(69)

и исследуем поведение числителя в правой части равенства (69). Получаем, что при а < а (8)

G'(ä(S))

При a = ä(8)

G2^(a)-G2(a) = 0,

(70)

(71)

С'(а(8)) и при а>ст(8)

О2 (ст(5))

—; 1 \'о'(а)-в2(о)<0 (72)

0'(а(5)) К ’ К ’

Из соотношений (70)-(72) следует, что

б (<*) / / чч

тах -----. ; 1. = С(а(8)),

ОйафЦ стС (ст) + а (5) V V //

а из (21), (67) и (68), что

||%-^а(б)Я5||^^(а(8)). (73)

Ввиду произвольности элементов и0, gs из соотношения (73) следует утверждение теоремы.

Из теоремы 1 и формулы (47) следует:

Теорема 2. Пусть для любого сге(о, ||у4||) С1 (а\

С'\а)>0, функция — |-у возрастает, 8<гО(|^|)|<4||, G2Гa(S)l

а а(8) = —-Г . Тогда метод М.М. Лавренть-G'La(8)J

ева , определяемый (49), оптимален.

3. Обратная задача Коши для уравнения теплопроводности

Рассмотрим уравнение теплопроводности ди(х,

dt

■ = Aw(x, t),

(74)

где X =(х1,х2,..., х„)еД„, ¿е [о, г], Т> 0 и

А $д2и

А« = I—-/=1 дх;

Предположим, что температурное поле и(х, Т) в момент времени Т нам известно приближенно, а распределение температуры и0 (х) = и(х,¿0 ) в момент времени ¿0 < Т требуется определить.

Далее обозначим и(х,Т) через /(*) и предположим, ЧТО и0 И / е ¿2 (Ж„ )

Поставленную таким образом задачу называют обратной задачей Коши для уравнения теплопроводности. Так как эта задача [6] является примером некорректно поставленной задачи, то мы предположим, что при точном значении /0(х) существует начальное распределение температуры

у0 (*) = и(х,0), которое принадлежит простран-ству (К„ ) и

<г.

(75)

Пусть точное значение /0 (х) нам не известно, а вместо него дано некоторое 8 - приближение /5 (х) е ¿2 (К„ ) и уровень погрешности 8 > О такой, что

(76)

Требуется, используя исходные данные (У^,8,г) задачи, построить приближенное решение иъ (х) Е ¿2 (К„ ), наиболее близкое к точному решению Щ (х) , И оценить уклонение ||г/5 -^оЦ^ .

Для решения поставленной задачи воспользуемся преобразованием Фурье (стр. 412 [7]), которое определяется следующим образом

1 , ,_ч

О \ //

,, I и(х)б

/_ \п/2 ^ V /

(2п) к„

’с£с, (77)

где /=>/ч, А, = (А1,Х2,...,Хй)еК„, а (х,х) = £х,х/,.

/=1

Обратное преобразование Фурье Ф“1 будет иметь вид

1 с Л / —\ /(а.,л)

»V / /

ы(х) = Ф_1 [и

\п/2

| й(я,)е^</Х.(78)

(2я)"-"к.

Хорошо известно, что преобразование Фурье Ф отображает пространство 1^ (Кя) на себя и при

этом имеет место теорема Планшереля [7, с. 412], которая утверждает изометричность этого преобразования, Т.е. ДЛЯ любого и е ¿2 )

НЧ, -К, • (79)

Теперь используем преобразование Фурье Ф для решения уравнения (74). Для этого дополнительно предположим, что для любого t функции и(х^),

их (х,/) и иП2 (х,/) абсолютно интегрируемы по

всему пространству Ки, функция щ (х,^) имеет интегрируемую мажоранту \}/(х),

(80)

Бир к/(х,/) < 1|/(х) *е[0 ,Т]

Бир

/е[0,Г]

— | и(х,/)е !^*х}сБс

Таким образом, преобразование Фурье Ф переводит уравнение (74) в обыкновенное дифференциальное уравнение первого порядка

й;(м)=-1Л2й(м).

/=1

Решая уравнение (82), получим

(82)

-1#

й{м) = у0(х)е 1=1 , (83)

где у0 (X) = Ф [у0 (х)], а у0 (х) = и (х, 0) .

Из (83) следует, что

й0(х) = у0(х)е'=> , (84)

а

п

_ _

/о(х) = у0(^)е -1 . (85)

Таким образом, из (84) и (85) следует, что операторы В и С = АВ могут быть представлены формулами

Ву(^ = е '-1 ^(^-)>

п

-У х?т Су(1) = е% ур),

(86)

(87)

где у(^), и Су^к)еЬ2 (К„), а оператор-

ное уравнение, которое нам предстоит решать, будет иметь вид

Ай(к) = е& й(1) = /(х). Из (86) и (88) следует, что В = С(А), а из (86), (88) и (89) следует, что

в (а)

(88)

(89)

(90)

Используя метод М.М. Лаврентьева для решения уравнения (88), получим, что

-1

-1^

+ а

(91)

где /8(а.) = Ф[/в(х)]

Следуя (21), рассмотрим уравнение

&(ст)ст = —.

(92)

Из (90) следует, что решение а (8) уравнения (92) имеет вид

. .1^1о

(93)

с(8) =

г),(9о; о (8) =

Из (67), (90) и (93) следует, что Т~10(дЛ

'о *.г

(94)

Из (91) и (94) следует, что приближенное решение уравнения (88) будет иметь вид

ГО-

т^х}

о г=1

*0 2

/=1

/=1

-ЛМ (95)

е /=1 + а(5)е

Из (73), (90) и (93) следует, что для приближенного решения определяемого формулами

(94) и (95), справедлива оценка

II -/ ч и Тч°

-а(5) .

- и0 < г т 8 т .

(96)

Применяя к щ ^ обратное преобразование Фурье Ф-1, получим приближенное решение щ (х) задачи (74)-(76) для которого будет выполняться оценка (96).

Заключение

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

Литература

1 Лаврентьев, М.М. Об интегральных уравнениях первого рода / М.М. Лаврентьев // ДАН СССР - 1959. - Т. 127, № 1. - С 31-33.

2. Иванов, В.К. Об оценке погрешности при решении линейных некорректно поставленных задач / В.К. Иванов, Т.И Королюк // ЖВМ и МФ -1969 - Т. Я № 1. - С. 30-41

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

4. Менихес, Л.Д. Конечномерная аппроксимация в методе М.М. Лаврентьева / Л.Д. Менихес,

B.П. Танана // Сибирский журнал вычислительной математики. - 1988, - Т 1, М 1 - С. 59-66.

5. Люстерник, Л.А. Элементы функционального анализа / Л.А. Люстерник, В.И. Соболев. — М.. Наука, 1965.

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

7 Колмогоров, А.Н. Элементы теории функций и функционального анализа / А.Н. Колмогоров,

C.В. Фомин. - М.. Наука, 1972.

Поступила в редакцию 1 июня 2009 г.