CC BY
26
УДК 517.948
ОБ ОЦЕНКЕ ПОГРЕШНОСТИ В ТОЧКЕ ПРИ РЕШЕНИИ ОБРАТНЫХ ЗАДАЧ
В.П. Танана, Т.С. Камалтдинова
Предложен новый подход к оценке погрешности нелинейного метода проекционной регуляризации в точке при решении некорректных задач. Приведено сравнение этой погрешности на множестве.
Ключевые слова: некорректная задача, метод проекционной регуляризации, оценка погрешности, модуль непрерывности, гильбертово пространство.
Введение
При оценке погрешности методов решения некорректно поставленных задач приходится сталкиваться с трудностью, связанной с неопределенностью точного решения. Выход из этого положения заключается в том, что дополнительно вводится класс корректности и предполагается, что точное решение задачи принадлежит этому классу. Затем (см. [1]) оценка погрешности определяется на этом классе. Естественно, что при таком подходе, для получения оценки, из класса корректности выбирается «худший» элемент. Этот подход приводит к тому, что реальная погрешность оказывается меньше теоретической.
В настоящей работе, при условии кусочной гладкости точного решения, сделана попытка сравнить оценку погрешности в точке с оценкой на классе. Для этой цели использован нелинейный метод проекционной регуляризации [2].
1. Постановка задачи
Пусть Ь2 (-да,да) - комплексное пространство, а Z с Ь2 (-да,да) . Элемент и(х) е Z тогда и только тогда, когда
П
и ( X ) = ( X ) + ^ ( X ),
где
I а (х2 - х2), |х < xi щ (х) = ^ V >
0, |х| > xi
аг и х1 > 0, а щ (-х) = щ (х) и
у (х) е Ж/2 (-да, да).
Рассмотрим линейные, неограниченные, замкнутые операторы Т и G с областями определения В (Т) и В (С) с Ь2 (-да, да) и множествами значений
Я (Т) и Я (в) с Ь2 (-да, да) такие, что
Я (Т) I В (G) = 12 (-да, да),
где Я (Т) IВ (G) - замыкание множества Я (Т) IВ (О) .
Поставим задачу вычисления значения 7/0 оператора Т в точке /0 & В (Т)
Т/ = и. (1)
Предположим, что при / = /0 элемент и0 = ТУ0 принадлежит множеству 2 , но точное значение /0 нам не известно, а вместо него даны & L2 (-да,да) и 6 > 0 такие, что
||Л -/о|| ^ 3. (2)
Требуется, используя априорную информацию / и 8, определить приближенное значение иг е Ь2 (-да,да) задачи (1), (2) и оценить уклонение - ы0||.
2. Основные понятия
Определение 1. Семейство {Т :0 < 8 < 80} операторов Т8 , непрерывно отображающих Ь2 (-да,да) в себя, будем называть методом решения задачи (5), (6), если для любого и0 £ 2
sup
{||TS fs - Uo|| : fs e L2 (-«>)>1 fs - T"Ч|| ^ 5} ^ 0 nPU 5 ^ 0 ■
fs
Введем оценку погрешности As (u0, S) метода {Ts :0 < S < S0} в точке u0 e Z
AS (u0,S) = sup{||Ts/s - uo\\: fs e L2 (-да>)> fs - T~4 < S} ■
fs 1 ’
Обозначим через Mr множество, определяемое формулой Mr = {u : u e R(T) П D(G) ,||Gu\\ < r} .
На нем определим модуль непрерывности ю (т, г ) (см. [1]), формулой
ю{т,г) = sup{\Tf\| : f e T-1 (Mr),||/|| < t}; t,r > 0, (3)
f 1 J
а также оценку погрешности As (Mr, S) метода {7 :0 < S < S0} на множестве Mr
As (Mr ,8) = sup [||Ts fs -u0\\ : u0 e Mr, /y e L2 (-да,да), /у -Г-1м0 < 8}.
uo,fs
Известно, что для любого метода {7 :0<S<S0} справедлива оценка
As (Mr, S) > ® (S, r) ; 0 < S < 4) ■ (4)
3. Нелинейный метод проекционной регуляризации
Для решения задачи (1), (2) перейдем от функций и ( x ) и f ( x ) к Фурье-образам
и (4 ) = F [и (x )] и f {4 ) = F \_f (x )] , где F - преобразование Фурье в L2 (-да, да) .
При этом переходе операторы T и G перейдут в операторы T и G, действующие в пространстве L2 ( -да, да ) и определяемые действительными функциями t (4 ) и g (4 ).
Tf ц)=t а) f а) ; f ц) G d (t)
D (T) = {f (4) : f (4) 6 L2 (-да, да), t (4) f (4) e L2 (-да, да )}
и
Gû(4) = g(4)u(4); u(4) e D(iî),
D(<5) = {û(4) : il(4) e L2 (-<^да),g(4)(4) e L2 (-да,да)}
При этом множество Mr перейдет в Mr = F [Mr ], определяемое Mr = {w (4 ) : û (4 ) e R (T ) П D (GG ) ,|| Gù (4 )| < r}.
На Mr, аналогично (3), определим модуль непрерывности со(т,r) формулой
а(т, r) = sup{ Tf : f є T 1 (Ыг), f < т}.
Из теоремы Планшереля следует, что задача (1), (2) перейдет в метрически эквивалентную задачу
/ = и. (5)
Предположим, что при / = /0 элемент й0 = ТТО принадлежит множеству2 = ^ [ * ] , но вместо /0 нам известны е Ь2 (-да, да) и 3 > 0 такие, что
fs /о
(6)
В дальнейшем предположим, что функции ? ) и g ) непрерывные, четные, поло-
жительные и строго возрастающие на [0,да). Кроме того, ?(4) ^ и g^ да при £ ^ да
Нелинейный метод проекционной регуляризации [2] зададим семейством {Т :0 < д < д0} операторов Тв , непрерывно отображающих Ь2 (-да, да) в себя и определяемых формулами
fa (f ,*)f (4); fs (4) > 33
0; fs (4) < 3S
определяется уравнением
т8 и (4) =
■Г Ш ^ 11* ", « > 0 , а * (Л, *)
[0; ||^| > а ^ ’
* 2
!\/3 ц) 14=9£2.
а
Таким образом, приближенное решение йд ) задачи (5), (6) определим формулой
и (4) = 7 Л (4) ■
4. Оценка погрешности \й8 - й0||
Лемма
g2 (4)
1. Если Z с U nMMr, а функция —убывает на [0,да) , то n=i 1 + 4
г«2^^
Теорема 1. Предположим, что выполнены условия леммы 1, й0 е 21 Мг и для любого k > 0 найдутся числа dk, &к и 4к > 0 такие, что для любого 4 , удовлетворяющего условию Щ > 4к справедливы соотношения
£ (4)
dk dk•
g (к4)
со (8, r )
Теорема 2. Пусть выполнены все условия теоремы 1. Тогда As («о)
А, (Mr ,S)
5. Решение обратной задачи Коши для уравнения теплопроводности
Рассмотрим уравнение
ди (х, г) д2и (х, г) . ,
—|—'- = —;-да< х <да, г е ( 0, ?о ], ?о > 0, (7)
, ^ ^ (°, , 10
дг дх
где и (х,I) е С{(- да,да)х[0,?0]} IС21 {(-да,да)х(0,?0]} для любого ^ е(0,?0], и(х,?), и'х (х,?), иXX (х,?) е Х2 (-да,да) I Ц (-да,да) и существует х(х) е L1 (-да,да) такая, что почти для любого ? е (0, ?0 ]
, , ди (х, ?)
|и (х, ?)| ч-д— ^ X (х), х е (-да, да).
Пусть нам дано распределение р (х) е Ь2 (-да, да) I 7 (-да, да) в момент времени /0
и(х,?0) = ф(х); -да< х<да, (8)
а начальное распределение и0 (х)
и(х,0) = и0 (х); -да< х <да (9)
требуется определить.
Предположим, что при ф(х) = р0 существует и0 (х) е Z , такое, что решение задачи (7),
(9) удовлетворяет условию и(х,?0 ) = ^>0 (х) , но точное значение ^>0 (х) нам не известно, а вместо него даны некоторое приближение ф3 (х) е Ь2 (-да,да) 17 (-да,да) и 6 > 0 такие, что
\<Ра (х) - Ро (х)||^ 3• (10)
Требуется, используя исходные данные ,8) задачи (7), (8), (10), определить приближенное решение иу (х) е Ь2 (-да,да) и оценить величину |и8 (х) - и0 (х)||£ .
Применив к задаче (7), (8), (10) преобразование Фурье ^, сведем ее к задаче вычисления значений неограниченного оператора Т
ТФз (^^и^)• (11)
|р (х ) - ро (х )||^ 3• (12)
Приближенное решение й3 (Л) задачи (11), (12) определим формулой
й (х\Л еАоф*(я); И ^а(ф*,д)
* ( ) = [0; И > а (Ф*,дУ
где а (ф8,8) удовлетворяет условию
да
| \фв (X)2аX = 982
а (Ра )
Пусть Х = F [ ^ ], Мг определим формулой
М = {и (Л): и (Л) 6 Я (т) ,|| g (Л) и (Л )|| < г},
где g(Л) = g(—Л), g & С (-да,да), g(Л) > 0, строго возрастая, стремится к бесконечности.
¿Ъ(S, г) = sup{ Тф : ф еТ 1 (Mr) ,|\ф\ < S} ,
ч>
Ag (йо,8) = sup{|iis (X) -uo (X)||: (ps є L2 (-да>)> (s -T~4 ^ 8}
<Ps
a Ag (Mr ,8) = sup {||ig - i0||: i0 є Mr ,ф8 є L2 (-да,да)
«0 $5
- T i0
< 8
}■
Теорема 3. Предположим, что Z с UnMr , —убывает на [0,да) , й0 (X) є Z I Mr
n=l 1 + ^
и I'M >33 ■
Тогда
As (u0,8)
со (8, г )
Теорема 4. При выполнении условий теоремы 3 справедливо соотношение
А3 (и0 )
— 0 при 8 — 0.
Заключение
Введена оценка погрешности в точке при решении некорректных задач. Приведено сравнение этой погрешности с погрешностью на множестве. Другие методы решения задачи (7)-(9) и ее нелинейного варианта рассмотрены, например, в работе [3].
да
Литература
1. Ivanov, V.K. Theory of Linear Ill-Posed Problems and its Applications / V.V. Vasin, V.P. Tanana. - VSP, 2002
2. Танана, В.П. Об оценке погрешности приближенного решения одной обратной задачи в классе кусочно-гладких функций / А.Б. Бредихина, Т.С. Камалтдинова // Труды Института математики и механики УрО РАН. - 2012. - Т. 18, №1. - С. 281288.
3. Табаринцева Е.В. Об оценке погрешности метода квазиобращения при решении задачи Коши для полулинейного дифференциального уравнения / Е.В. Табаринцева / / Сиб. журнал вычисл. матем. - 2005. - Т. 8, № 3. - С. 259-271.
Танана Виталий Павлович, заведующий кафедрой вычислительной математики, Южно-Уральский государственный университет, [email protected]
Камалтдинова Татьяна Сергеевна, старший преподаватель, Южно-Уральский государственный университет, KamaltdinovaTS@ mail.ru.
ON POINT-WISE ERROR ESTIMATE IN SOLVING INVERSE PROBLEMS
V.P. Tanana, South Ural State University (Chelyabinsk, Russian Federation), T.S. Kamaltdinova, South Ural State University (Chelyabinsk, Russian Federation)
A new approach to the point-wise error estimation for the projection ill posed problems is suggested in the article. We compare point-wise error estimate with the error estimate on a set.
Keywords: ill-posed problem, method of projective regularization, error estimation, module of coninuity.
References
1. Ivanov V.K., Vasin V.V., Tanana V.P. Theory of Linear Ill-Posed Problems and its Applications. VSP, 2002.
2. Tanana V.P., Bredikhina A.B., Kamaltdinova T.S. On an error estimate for an approximate solution for an inverse problem in the class of piecewise smooth functions. Supplement to Proceedings of the Steklov Institute of Mathematics and Mechanics Ural Branch of RAS. 2012. Vol. 18, No. 1. P. 281-288.
3. Tabarintseva E.V. On error estimate of the quasi-reversibility method to solve the Cauchy problem for a semielinear differential equation. Siberian journal of numerical mathematics. 2005. Vol. № 3. P. 259-271.
Поступила в редакцию 14 января 2013 г.
