![](/user_photo/_userpic.png)
книги из ГПНТБ / Романов В.Г. Обратные задачи для дифференциальных уравнений спецкурс для студентов НГУ
.pdf
|
S K « S « S K ^ , |
K = 0, 1 , 2, |
Jf-i |
|
||||||
является в пространстве |
|
|
f [ 5 M , s K ^ J |
оператором сжатия.Дей- |
||||||
сгвительно, для любых |
cCj, эс2 е |
CcsK, S K H ] |
|
|
|
|
||||
|
|
|
|
s |
|
|
|
|
|
|
J^x, |
- 4^г« |
= |
I J X(s.l)[xt{l) |
- |
xz(l)]di |
I < |
|
|||
< maoc |
J i^f(5,f)Ul |
| x - |
x |
|| < |
|
|||||
« / 5 ^ , - s J maoc |
I T f ^ l - f ^ - x j |
= |
d-ls^-sj-ixt-xj |
. |
||||||
Но по условию разбиения |
• |
^ ( s ^ - . s j |
<i |
|
для всех к . Поэто |
|||||
му каждый из операторов |
Лк |
действительно является |
оператором |
|
||||||
сжатия. Решая последовательно |
уравнения ( I I ) , |
найдем единственное |
||||||||
непрерывное решение уравнения |
( 8 ) . |
|
|
|
|
|
||||
Заметим, что |
если к уравнению |
(8) |
применить непосредственно |
|
метод последовательных приближений, то можно установить, что он
сходится сразу |
на всем интервале [o,iJ |
(см.[104,110,137]). |
|||
|
3 . Рассмотрим нелинейное |
интегральное уравнение |
|
||
|
|
осМ = А- J3T(s. i x(l))dl, |
(13) |
||
|
|
о |
|
|
|
где |
функция |
I, х{1)) |
непрерывна |
по совокупности |
аргумен |
тов |
в области |
OiS,{ $ i, |
i l l , |
и пусть в этой |
области |
|
|
11^ эьлх) | ^ Q . |
(14) |
Покажем, что уравнение (13) имеет единственное непрерывное на от
резке |
[o,i] |
решение, |
удовлетворяющее неравенству |
I x u i / , |
если |
|
|
|
< |
ПЫП | X ; |
— |
— | |
(15) |
|
|
|
max j max |
№(s,t,x)!tW |
|
|
Пусть |
^ - |
оператор: |
|
|
|
|
30
Возьмем в пространстве |
Clo,i] |
любой |
элемент, принадлежащий |
||||||||||||
шару |
5fo, М) |
|
радиуса |
М |
с центром в точке |
о |
. Покажем, |
||||||||
что оператор |
с# |
при условии |
(15) переводит |
этот |
элемент |
в |
эле |
||||||||
мент .также принадлежащий шару |
S(o, М). |
Действительно, |
|
|
|||||||||||
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
MacL |
« Ш- |
та-х j |ЯМ,а:^Ы<й |
4 |
|
|
|
|||||||
|
|
|
|
|
|
О |
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
|
|
|
|
|
|
|
|
|
|
|
< |
Ш- |
max |
} |
max |
lX(s,l,x)l |
ctt |
< J/. |
|
|
||||
Покажем теперь, |
что при выполнении условия (15) |
оператор |
Л |
на |
|||||||||||
шаре |
S(o,JUj |
|
является |
оператором |
сжатия. Пусть |
х±,хг-любые |
|||||||||
два элемента |
шара |
|
S(o,JU). |
Тогда |
|
|
|
|
|
|
|
||||
Щ - i o r J |
=Ш- max | |
feftU*.tftf- J t l s , l . x , a i ) ] di |
|
I 4 |
|||||||||||
|
=ё lAl-tnax |
|
j lf=. 3¥(^,x)l |
|
I • /эс,(1) - |
a: |
N)| eH |
^ |
|||||||
Здесь |
cc*(l) |
- |
фиксированное число, при каждом фиксированном i , |
||||||||||||
заключенное между |
числами |
|
x^l) |
и |
хг(1) |
(средняя точка |
в тео |
реме о конечных приращениях). Полученное неравенство показывает, что оператор Л на шаре 3(о,М) есть оператор сжатия. Из ус тановленных фактов и следует высказанное ранее утверждение о су
ществовании единственного решения уравнения |
( 1 3 ) . Отметим, |
что |
вместо интегральных уравнений ( 5 ) , ( 8 ) , (13) |
мы могли бы рассмат |
ривать соответствующие им системы интегральных уравнений. Измене ния, которые следует при этом внести в рассуждения, чтобы устано вить аналогичные теоремы единственности и существования, доволь но очевидны.
§ 4. О корректности прямых и обратных задач для дифференциальных уравнений
В этой лекции мы рассмотрим очень важный вопрос, связанный с корректностью задач для дифференциальных уравнений. Обычно он рассматривается в курсах математической физики (.см., например, [73,98,112,138]) . Мы приведем здесь краткое изложение этого воп роса. Связано это с тем, что при рассмотрении обратных задач для дифференциальных уравнении оказывается целесообразным несколько
31
изменить понятие корректности задачи.
Из теории обыкновенных дифференциальных уравнений и уравнений в частных производных известно, что дифференциальное уравнение оп ределяет целое семейство решений, зависящее либо от некоторого числа произвольных постоянных, либо от некоторого числа произволь ных функций. Чтобы задача имела определенный физический смысл, не обходимо выделить из всего класса решений дифференциального урав нения то единственное решение, которое описывает данный физичес кий процесс. Такое выделение осуществляется обычно заданием до полнительной информации о решении, чаще всего это задание началь ных и граничных условий для дифференциального уравнения. Поясним это на примере. Уравнение
|
|
|
• « и - с Ч * х |
|
|
|
( I ) |
описывает |
процесс малых колеоаний струны. При этом функция |
||||||
и(х, 1} |
|
равна смещению струны в точке |
х |
в момент |
времени i от |
||
положения равновесия, совпадающего с осью |
х . Числовой коэффици |
||||||
ент с |
характеризует скорость передачи |
сигналов в струне (ско |
|||||
рость |
звука). В случае конечной струны |
ахаХ, |
закрепленной |
||||
в точках |
х=о |
и х=Х, естественные |
граничные условия для урав |
||||
нения |
( I ) |
имеют вид: |
|
|
|
|
|
|
|
|
n(o,i)^o, |
u[X,L}=o. |
(2) |
Однако задание условии (2) не достаточно для выделения единствен ного решения уравнения ( I ) , Чтобы оно было полностью определено, необходимо задать еще начальные условия, которые сводятся в дан ном случае к заданию смещений струны и начальных скоростей ее то чек:
и (х, б) = Cf>(X), ыг (х, ol = (j)lx). ' (3)
В курсах математической физики показывается, что задание началь
ных условии (3) и граничных условий (2) однозначно |
определяет ре |
|||
шение уравнения ( I ) . |
|
|
|
|
В случае уравнения теплопроводности |
|
|
||
|
щ=ли, |
|
(4) |
|
рассматриваемого в конечной области пространства |
х |
=[xa,xz,x3)t |
||
ограниченной поверхностью |
S, |
для выделения однозначного решения |
||
достаточно задать тепловой |
режим на поверхности 5 , например: |
|||
nlxDJ |
= / М ) , |
|
(5) |
32
и начальное распределение температуры внутри .5: |
|
nivz, о) = ср(ос). |
(6) |
Итак, задание начальных и граничных условий преследует цель выде ления из всего класса решений дифференциального уравнения единст венного решения. Вместе с тем таких условий должно быть минималь ное количество, иначе они могут вступить в противоречие друг с другом. Другими словами, задаваемые условия должны обеспечивать существование решения дифференциального уравнения при произволь ном выборе функций, входящих в эти условия, из некоторого функцио нального множества.
Подчеркнем еще одно существенное обстоятельство. Рассматривая данное дифференциальное уравнение с заданными граничными и началь ными условиями, мы можем ставить задачу отыскания его решений.при надлежащих, вообще говоря, различным функциональным пространствам. Выбор их диктуется, как правило, физическим смыслом задачи. Напри-
мер, можно рассматривать задачу |
отыскания решения уравнения |
( I ) |
||
при условиях ( 2 ) , |
(3) в классе |
дважды |
непрерывно-дифференцируемых, |
|
в любой конечной |
области пространства |
эс, i функций, я можно |
рас |
сматривать задачу об отыскании решении того же уравнения, при тех же самых условиях, принадлежащего множеству дважды дифференцируе мых функций, суммируемых со степенью р . Другими словами, мы можем выбирать достаточно произвольным образом функциональное простран ство решений дифференциального уравнения. При этом функции, входя щие в граничные и начальные условия, не могут быть произвольными, они должны обеспечивать принадлежность решения дифференциального уравнения выбранному функциональному пространству. Для этого они сами должны принадлежать некоторому функциональному пространству (вообще говоря, не совпадающему с пространством отыскиваемых ре шений). Это становится особенно понятным, если мы посмотрим на за дачи для дифференциальных уравнений с точки зрения функционально го анализа. Действительно, выберем функциональное пространство U решений дифференциального уравнения. Дифференциальное уравнение в совокупности с начальными и граничными условиями определяет опе
ратор Л, который ставит в соответствие любому решению ие U на
бор функций, входящих в дополнительные (начальные |
и граничные)ус- |
||||
ловия. Для уравнения ( I ) |
это функции . <pfccj, <р(рс), |
а для уравне |
|||
ния |
(4) — функции f(x,i), |
(floe). |
Рассматривая |
этот набор функ |
|
ций |
как элемент ]• функционального пространства F , |
мы получаем, |
|||
что |
решение задачи для дифференциального уравнения |
|
эквивалентно. |
33
решению операторного |
уравнения |
|
||||
|
|
|
|
|
du = / |
(7) |
при условии, что |
« е |
U |
. Чтобы решение этого уравнения |
существова |
||
л о , |
необходимо, |
чтобы |
/ |
являлось образом некоторого элемента |
||
ueU, |
т . е . принадлежало |
множествуi являющемуся областью значений |
оператора d . Таким образом, множество данных задачи вполне опре
деляется заданием пространства решений U. |
|
|
Условия существования |
и единственности |
решения уравнения (7) |
гарантируют существование |
обратного операторам? : |
|
|
u = d"f, |
(8) |
который и решает задачу, ставя в соответствие данным задачи,т.е. элементу / , решение ueU, Необходимо, однако .вспомнить, что ос новная цель, для которой решаются математические задачи, заключа ется в математическом описании физических процессов. При этом данные задачи берутся из опыта и они не могут быть измерены аб солютно точно, то есть в данных задачи всегда присутствуют ошибки измерения. Поэтому, чтобы математическая постановка физической за дачи описывала реальный физический процесс, к задаче необходимо предъявить еще некоторые требования, которые отражают физический факт малого измерения решения при небольшом изменении данных, или, как принято говорить, устойчивости решения к малым возмущениям входных данных. Сформулируем сказанное выше на принятом в матема
тике языке. Поскольку |
мы видели, |
что прямые задачи для дифферен |
|||||
циального |
уравнения приводят к решению уравнения ( 7 ) , |
то мы сфор |
|||||
мулируем понятие устойчивости решения для этого уравнения. |
|||||||
Итак, |
пусть рассматривается |
уравнение (7), где |
Л |
- |
оператор, |
||
действующий из некоторого метрического пространства |
U в метричес |
||||||
кое пространство F. |
Говорят, что |
решение уравнения |
(7) |
устойчиво |
|||
к малым изменениям правой часта |
fe |
F уравнения ( 7 ) , |
если для лю |
||||
бого е>о |
найдется такое. 8>о , |
что для любого элемента |
F для |
||||
которого |
|
•Mb |
|
< *. |
|
|
|
выполнено |
неравенство |
|
|
|
|
||
jXjl-u. й)< е. |
|
|
|
||||
|
|
|
|
|
|||
При этом |
|
Для нелинейного оператора |
А |
устойчи |
вость зависит, вообще говоря, от элемента / . На одних элементах оператор может быть устойчив, на других нет. В случае линейного оператора А устойчивость или неустойчивость имеет место сразу для
34
всех элементов {е. F.
Математическая задача решения уравнения ( 7 ) , удовлетворяющая требованиям существования, единственности и устойчивости к малым изменениям входных данных носит название корректной задачи. Сфор
мулируем это понятие более точно. |
|
|
|
|||||
|
Пусть Л - |
оператор, действующий из метричеокого пространства |
||||||
U |
в метрическое пространство F. |
Задача решения уравнения |
(7) на |
|||||
зывается поставленной |
корректно |
на паре пространств ( J , F, |
если |
|||||
она |
удовлетворяет следующим требованиям: |
|
|
|
||||
|
1) |
Решение |
задачи |
существует |
при любом |
/е F |
и принадлежит |
|
пространству |
U, |
|
|
|
|
|
||
|
2) |
Решение |
задачи |
единственно в пространстве |
V. |
|
||
|
3) |
Решение |
задачи |
устойчиво |
на любом элементе |
-f^F. |
|
|
|
Отметим, что задача может быть корректна |
на одной паре |
прост |
ранств и некорректна на другой. Ясно, например, что при расшире
нии пространства F можно прийти в противоречие |
с требованием су |
||||||
ществования решения при любом /е F |
. В том случае, когда |
оператор |
|||||
d линеен, для корректности задачи |
(7) |
на паре |
банаховых |
пространс |
|||
ранств U, |
F необходимо |
и достаточно, |
чтобы для |
оператора |
Л |
су |
|
ществовал |
ограниченный |
обратный Л |
, действующий из F |
в |
U,при |
чем область определения обратного оператора совпадала с простран ством F (см. [98J. стр.507).
Приведем пример корректной задачи. Рассмотрим для простоты не
ограниченную струну - о о < ос < + = > о . |
В этом случае достаточно |
для |
||||
уравнения |
( I ) |
задать начальные условия в виде ( 3 ) . Будем считать, |
||||
что функция |
if(-x) |
является дважды непрерывно дифференцируемой, |
||||
а функция |
ф(ос) |
один раз непрерывно дифференцируемой в любой ко |
||||
нечной области, т . е . при любых конечных а и 6 |
срсх) е |
Cz[a,S]t |
||||
c^fcje С |
|
Пусть решение u(ocj) |
уравнения ( I ) |
разыскивается |
в классе функций, имеющих непрерывные частные производные до вто
рого порядка включительно в области eD={-X<ix<X, |
o^laT), |
где Z |
и Т - произвольные положительные числа, т . е . |
и&€2аЬ). |
Покажем, |
что задача поставлена в этом случае корректно. Действительно, не трудно убедиться непосредственной проверкой, что решение уравне
ния ( I ) , при условиях (3)., в |
случае безграничной |
струны имеет |
вид: |
||
ip(x-cl) |
+ cp(-x+ci) |
i |
Г , |
,., |
(9) |
«ftr,«= г |
|
+ _ |
ffmidl. |
x-el
35
При этом принадлежность |
<ре С2 1-{'Л+сТ), (Х+сТ)] |
, |
|
фе Ci[-lj£->-cT), |
(Х+сТ}] |
гарантирует существование непрерывных |
|
частных производных до второго порядка включительно в области |
|||
Единственность решения задачи ( I ) , (3) в классе С2(Ъ) |
обычно |
доказывается в курсах математической физики. Покажем теперь, что
решение задачи устойчиво к малым изменениям функций |
(fix) |
и ф(х) |
||
в соответствующих функциональных |
пространствах С г и |
С. |
Введем |
|
для этого нормы соответствующих пространств: |
|
|
||
Icpll "max |
^[\(f{x)\ |
+ \(p'm\ + \ц>"1Х)\] , |
|
|
Проводя оценки для уравнения ( 9 ) , находим
Ы(-х, III « max Шх)\ + Т max / ф(х) \ « ШЦ+ Т 1ф1.
Дифференцируя равенство (9) и делая аналогичные оценки, получаем
1и{(х.1)\ < с-ЩЧ + «ф№ •
\u^x.i)\ « 1Кр!1 + -£-1фП;
lvH(x,l)l |
<сг1срЦ + сйф$. |
Складывая все полученные неравенства, получаем, что
HuH^(3+2c+cz) |
llcpil +(Т+2+с + -*-Иф11. |
(Ю) |
Возьмем теперь произвольное е>о . Тогда, если
то |
|
|
|
|
и, следовательно, .чалым по норме функциям |
(f и ф |
соответствуют |
||
малые решения |
и[х,1) |
задачи ( I ) , ( 3 ) . |
Так как |
задача линейна, |
36
го отсюда следует, что малым по норме изменениям функций у я ф соответствуют малые изменения решения.
Заметим, что та же самая задача, если считать <f и ср просто непрерывными функциями, а для решения сохранить прежнее простран
ство |
С I'D), будет некорректной, так как не будут выполнены ус |
ловия |
корректности I ) и 3 ) . |
Приведем теперь пример задачи, некорректной по существу, ког да никакое требование на гладкость исходной информации не приво дит к корректности задачи. Пример такой задачи был построен Адамаром, чтобы подчеркнуть важность третьего из условий коррект ности.
Пусть |
а = и (ж, ц] - |
решение уравнения Далласа |
|
||
в долуполосе |
у>о , |
-7Г<х<тг, удовлетворяющее условиям |
|||
ut-зг. у) = и 1я, у) = |
tt(x,o) = о, |
и^х,о)=ётлт |
пх. (12) |
||
Нетрудно проверить, |
что решением этой задачи является функция |
||||
|
|
•utx,p |
=~e'm-suinx-sk |
nij. |
(13) |
Можно показать, что решение поставленной задачи единственно. При
п-^~= |
функция |
emsinn.x, |
представляющая из себя данные |
за |
||
дачи ( I I ) , (12) , стремится равномерно к нулю вместе со всеми про |
||||||
изводными. Тем не менее решение задачи |
( П ) , (12) , как видно |
из |
||||
формулы |
(13), при любом фиксированном |
у>о |
и достаточно боль |
|||
шом п |
имеет вид синусоиды со сколь угодно |
большой амплитудой. |
||||
И, следовательно, |
сколь угодно |
малому |
изменению данной задачи со |
ответствуют отнюдь не малые вариации решения. Итак, рассмотрен ная задача Коши для уравнения Лапласа поставлена некорректно.
Из приведенных примеров видно, что существуют, вообще говоря, два сорта некорректных задач. Задачи, которые являются некоррект ными в одних пространствах, но могут быть сделаны корректными при другом выборе пространств. Такие некорректные задачи естественно назвать слабо некорректными. Наряду с этим существуют некоррект ные задачи (как,например, рассмотренная задача Копш для уравне ния Лапласа), когда ни в каких функциональных дространсгБах.норма которых использует конечное число производных, задача не явля ется корректной. Задачи подобного-сорта будем называть сильно не-
37
корректными задачами. До сравнительно недавнего времени считалось, что некорректные задачи не имеют физического смысла. Это мнение, однако, опровергается развитием прикладных наук. В настоящее вре мя известно большое число практически важных задач, которые не являются корректно поставленными в смысле данного определения. К числу таких задач относится и большинство многомерных обратных задач для дифференциальных уравнений, на важность которых было указано во введении. В связи с этим А.Н.Тихоновым было предложено [153] новое понятие корректности, которое является физически оп равданным для многих прикладных задач,некорректных в классичеоком смысле.
для подхода А.Н.Тихонова к вопросу корректности характерно,
что рассматривается некоторое множество JllcU, |
существенно более |
|||||||||
узкое, чем все |
пространство U. |
Чаще всего |
это |
некоторое |
компакт |
|||||
ное множество. Пусть образ множества |
Ж при отображении с помо |
|||||||||
щью оператора |
d |
в пространстве |
F |
есть |
множество |
й. , |
то есть |
|||
1& = ДШ]. |
Задача |
(7) |
называется условно |
корректной |
(или |
коррект |
||||
ной по А.Н.Тихонову), |
если для него |
выполнены следующие условия: |
1)априори известно, что решение задачи существует и принад лежит некоторому заданному множеству Ж функционального прост ранства U;
2)решение задачи единственно на множестве Ж;
3) бесконечно малым вариациям / « Я , не выводящим решение за пределы множества М , соответствуют бесконечно малые вариации решения.
Множество Ж называется множеством корректности задачи.
Отметим некоторые особенности данного определения корректнос ти, отличающие его от классического. При классическом определении корректности предполагается, что данные задачи / принадлежат некоторому метрическому пространству F. Первое из требований классичеокого понятия корректности требует доказательства, что
при любом {е F |
уравнение |
имеет решение « е U. |
Это обычно делается установлением теоремы |
существования. Соответствующее требование корректности задачи по
Тихонову постулирует |
существование решения задачи, принадлежаще |
||
го некоторому заданному множеству JILc |
U , на некотором множестве |
||
данных. Относительно |
этого последнего |
множества не делается |
ни |
каких предположений, |
структура его не |
описывается. Тем самьм |
на- |
38
добность в доказательстве теоремы существования снимается. |
|
||||||
Требование единственности решения задачи на множестве Л |
со |
||||||
впадает о классическим требованием единственности |
задачи в |
про |
|||||
странстве решений |
U. |
|
|
|
|
|
|
Перейдем н третьему требованию. При классическом понятии кор |
|||||||
ректности решение |
задачи существует |
при любом |
/е.Р |
и принадлежит |
|||
пространству U. Естественно |
поэтому |
потребовать, |
чтобы малому |
||||
изменению элемента |
/ по метрике пространства |
F |
|
соответствова |
|||
ло малое изменение |
решения по |
норме пространства |
U. |
Корректность |
|||
задачи по Тихонову |
совсем не |
требует |
знания множества данных, но |
постулирует принадлежность его решения к некоторому множеству^/. Естественно при этом, что малое изменение данных задачи по метри ке любого наперед выбранного пространства может привести к тому,
что решение задачи |
вообще не существует или существует, но не при |
|
надлежит заданному |
множеству JU. В связи с |
этим разумно модифици |
ровать требование |
непрерывной зависимости решения от данных зада |
|
чи, потребовав эту |
непрерывную зависимость |
только для таких ва |
риаций данных, которые не выводят решение |
за пределы множества Ж. |
|
Это и отражено в третьем требовании корректности. |
Посмотрим, почему же понятие корректности по А.Н.Тихонову оказывается довольно естественным при рассмотрении многих задач и особенно при изучении обратных задач для дифференциальных урав нений. Расомотрим, например, обратную задачу о колебаниях струны. В том случае, когда плотность струны меняется от точки к точке, процесс колебаний струны описывается уравнением ( I ) , в котором
с = сбс). Пусть мы рассматриваем задачу определения функции |
cfoc) |
по некоторым известным функционалам от решений прямой задачи |
(I)- |
( 3 ) . В качестве таких функционалов здесь естественно рассматри вать смещения отдельных точек струны в различные моменты времени. При этом физика задачи подсказывает, что решение следует искать в классе функций положительных, ограниченных сверху и снизу неко торыми константами, которые можно указать, зная возможный набор материалов, использованных при конструировании струны. Возможно, что известны и другие свойства струны, например, что плотность ее меняется непрерывно, либо куоочно-непрерывно. В результате мы
заранее можем указать множество Л |
, которому должна принадлежать |
|||
функция |
c(xj . Если теперь данные |
задачи взяты из конкретного фи |
||
зического |
эксперимента «о |
струной, |
у.которой с(х)&М, |
то у нао |
не может быть сомнений в |
том, что |
существует реальная струна, ха- |
||
|
|
3.9. |
|