Файл: Романов В.Г. Обратные задачи для дифференциальных уравнений спецкурс для студентов НГУ.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.06.2024
Просмотров: 171
Скачиваний: 0
n.(x,y) = n„(oc,</) + nt Lx,y)t
где функция nc(xty) имеет в области £> ограниченные частные производные до второго порядка и удовлетворяет условиям
о<а ^ пв(х,у) |
< |
к о о , |
0<а^~Ъ^ |
n . ( 3 C >V ) 4 |
1 |
^ < ^ |
||||
а функция |
п.(х,у) |
мала |
по |
норме |
и пред ставима в виде |
|
||||
где (ркеС2[0,Ю, |
<РкеС*1- |
—,<=-). |
|
|
|
|
||||
Тогда |
в области |
£> |
|
функция п(х,у), |
в линеаризированной по |
|||||
становке, |
однозначно |
определяется |
временами |
г fx" ос*) |
|
для любых |
||||
пар точек |
х°х1е{у |
= |
о}. |
|
|
|
|
|
|
|
Для доказательства |
достаточно |
область •£) |
разбить |
|
на полоски |
|||||
толщиной не более h„ |
и воспользоваться |
следствием I |
§ |
3. Заме |
тим только, что на самом деле для единственности решения задачи
(7) нам достаточно знать функцию г , ( х ° х ' ) |
только |
для тех |
пар |
||
точек x°xie{y |
= o}, |
которым соответствуют |
лучи |
/7(£>??), |
их |
соединяющие и принадлежащие целиком области £>. |
|
|
|||
Можно получить аналогичные теоремы единственности для линеа |
|||||
ризированной |
обратной |
задачи, связанные с классами функций |
|
п 1( х .у). отвечающих следствиям 2 и 3 § 3 . Формулировка их и доказа тельство очевидны.
В процессе |
линеаризации задачи мы предполагали, |
что функция |
|
njx,y), |
мала и при этом условии получили приближенное равенство |
||
( 7 ) . |
Возникает, |
однако, вопрос, в каком смысле п± |
должна быть |
мала, чтобы решение линеаризированной задачи приближало нас к ис тинному решению? Вопрос этот очень важный в прикладном отношении. Нужно знать, в каком случае, решая приближенное уравнение (7),мы получим решение, заслуживающее доверия. Для случая, когда функ
ция п0 зависит только от расстояния до |
фиксированной точки, |
это |
|||
проделано в работе [1031 Получено достаточное условие, чтобы |
ре |
||||
шение уравнения (7) |
было близко к истинной функции |
гц. Мы приве |
|||
дем эти достаточные |
условия, |
перефразируя на случай |
п0=п0(у). |
Ус |
|
ловия эти следующие: |
функции |
п, и. п |
должны быть дважды непре |
рывно дифференцируемы, иметь ограниченные вторые производные, а
функции |
п. и [n-nj |
- удовлетворять |
условиям (17) . |
По-видимому, |
они являются достаточными и в случае, |
когда функция |
зависит |
||
от ас и |
у. |
|
|
|
184
В практическом отношении при решении обратной кинематической задачи линеаризацию задачи можно использовать для построения ме тода последовательных приближений. Так как уравнение (7) являет ся приближенным, то решив его, мы найдем только некоторое прибли
жение к функции |
п, . Объявим его |
первым приближением и обозначим |
||||||||||||||
через |
|
|
По формуле ( I ) найдем первое |
приближение функции |
п, |
|||||||||||
которое |
обозначим через |
п"': |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
п"'(х) |
= п„{х) |
+ |
n f ( x ) . |
|
|
|
|
|
|||
По функции |
|
n!''lx) |
построим лучи |
^ ( х ° х Ч |
и времена |
г*(х°х?). |
||||||||||
Тогда для |
разности |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
W X ) - п„м ) (х) |
= |
П ^ ( Х ) |
|
|
|
|
||||
получим приближенное равенство . |
|
|
|
|
|
|
|
|
|
|||||||
в котором |
|
|
|
|
|
|
1 |
' |
|
|
|
|
|
|
|
|
|
|
|
|
т*(х;х') |
= |
т(х;х*> |
- |
г*и° х*>. |
|
|
|
|
||||
Решая уравнение ( 7 ' ) , |
мы получим |
|
|
, по которому можно постро |
||||||||||||
ить далее |
п?=п0+п™, |
|
затем лучи |
[^(х'х*), |
|
времена пробега |
||||||||||
по ним |
т0'х>1х.°х') |
|
и,используя |
линеаризацию, получить |
уравнение |
|||||||||||
для третьего |
приближения |
п1* |
и т . д . Вообще, |
на |
«.-том шаге мы |
|||||||||||
получаем |
уравнение |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
< " * Л |
Й : Х * ) = |
5 |
с ^ ) |
|
^х, |
|
|
|
|
|||
в котором лучи |
Гк(х°х') |
строятся по функции |
|
|
|
|
||||||||||
|
|
|
|
|
|
п'0я1Сх)=Па1х)ч-гг^(х), |
|
|
|
|
|
|||||
а времена |
|
т^**" |
находятся по |
формуле |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
•rf'Vrx1 ) = Ttx'x1) |
г- |
T^lxXx1), |
|
|
|
||||||
причем |
r . ' ^ f x ' x 1 ) — времена пробега |
сигнала |
в среде, |
характери |
||||||||||||
зуемой функцией |
nl*(x). |
|
|
|
|
|
|
|
|
|
|
|
||||
Вопрос |
сходимости |
n^lxi |
к |
n&d |
требует |
своего |
выяснения:. |
|||||||||
Для сходимости, конечно, |
существенно, чтобы функции |
п!*Ш |
удон- |
|||||||||||||
летноряли условиям типа условий на- |
nix) |
леммы I . |
Можно наделать |
|||||||||||||
ся, что |
при |
«„(х), |
достаточно |
близком к |
л ( х ) , |
достаточно |
не |
|||||||||
скольких итераций, чтобы получить удовлетворительную |
точность. |
|||||||||||||||
4 . О единственности обратной кинематической задачи. Пусть |
об |
|||||||||||||||
ратная кинематическая, задача рассматривается в области |
Ц^о про- |
185
странства ос, у. |
Естественно было бы ожидать, |
что функция п(х,у), |
|
удовлетворяющая |
условиям леммы I |
однозначно определяется в облас |
|
ти £> по временам tlx'x'), |
х ° х * е {у=о}. |
К сожалению,это |
является пока гипотезой. Теоремы единственности задачи получены пока только в частных классах функций п(х,у), практически экви валентных классу аналитических функций. Различные исследования в
этом направлении |
были |
проведены |
Ю.Е.Аниконовым [9] ,[10] ,[12] , [ 1 3 ] , |
||||||||
[15] . Мы остановимся здесь только на одном результате. |
|
||||||||||
Т е о р е м а |
3 . |
Пусть п(х,у) — |
аналитическая |
функция в об |
|||||||
ласти цъо, и |
|
щх,о)>о, |
Пу(х,о)<о. Тогда она однозначно |
определя |
|||||||
ется в области |
у > о |
функцией |
т(х°х') |
для различных |
|
||||||
х°,х'е{ц |
= о}. |
|
|
|
|
|
|
|
|
|
|
Для доказательства |
теоремы достаточно показать, что различным |
||||||||||
функциям |
п(х,у) |
|
из указанного |
класса |
отвечают различные функции |
||||||
Т7х°х*). |
Рассмотрим две различных функции |
п±(х,у) |
и |
f\lx,y). Из |
|||||||
их аналитичности следует, что найдется |
такая |
точка |
(х,о) |
и такая |
|||||||
ее открытая окрестность |
5(х,а) |
={(х,у): f x - x ) I + ^ 2 < £ i , |
у>о}, |
||||||||
что для точек |
(х,у)е |
Six,г) |
одна из функций |
« ^ « ^ |
строго |
больше другой. Чтобы убедиться в этом, достаточно рассмотреть их
разность. |
Пусть, |
например, |
п^х,^) |
> nzlx,y). |
|
За счет |
умень |
||||||||
шения |
е |
всегда |
можно добиться, чтобы в замкнутой |
области |
|
5(х,£) |
|||||||||
были выполнены неравенства |
|
nilx,y)>o, |
^-nL(x,y)<o |
|
U=1,2). Тогда |
||||||||||
лучи |
/^(1°х'), |
отвечающие |
средам |
п€1х,у), |
в области |
|
5 ( х , е ) |
||||||||
имеют вид дуг, опирающихся на ось- х . |
|
|
|
_ |
|
|
|
|
|||||||
Покажем, что для лучей |
|
1?(х'хл), |
содержащихся в |
( |
|
|
), отве |
||||||||
чающие им времена |
удовлетворяют |
неравенству |
|
|
|
S |
X,E |
|
|
||||||
|
|
|
Tt(x°x') > |
тг1х°х'). |
|
|
|
|
|
|
|
||||
Действительно, из неравенства |
njx,y) |
> nz(x,y) |
следует, что |
||||||||||||
|
|
|
r , ( x ° x ' J = |
^ n1[x,y)ds |
|
>§nz(x,y)ds, |
|
|
|
|
|||||
|
|
|
|
£(х°х*) |
ГУх'х') |
|
|
|
|
|
|
||||
но так как из всех лучей, |
соединяющих точки |
ос° х* |
|
в среде .от |
|||||||||||
вечающей функции |
пг(х,у), |
минимум времени реализуется на кривой |
|||||||||||||
Г1(х°х*), |
то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
^ nzlx,y)ds |
> |
j |
i\(x,y)ds |
=тл(х°х'). |
|
|
|
|
|
||||
Отсюда и получаем требуемое неравенство. Но это неравенство |
озна |
||||||||||||||
чает, что функции |
Tiix'x'), |
|
отвечающие различным |
п^у), |
|
|
не |
186
могут совпадать между собой. Именно это нам и требовалось дока зать.
Слегка изменив рассуждения, проведенные здесь, можно устано вить теорему единственности в классе функций п(о^у), представимых в виде
|
п(х,у)=щх,у)-1>(х,у), |
П(Х,0)>0, |
tlytX,0)<0, |
|
где |
aix,y) — фиксированная положительная |
непрерывно-дифференци |
||
руемая функция, а |
&(х,у) |
- аналитическая |
в области у&о функ |
|
ция [ |
1 5 ] . |
|
|
|
5. О постановке |
задачи при наличии внутренних источников воз |
мущенна. Как мк уже видели при исследовании одномерной задачи и при исследовании линеаризированной постановки многомерной задачи условие монотонности скорости распространения волн ( п у < о ) очень существенно для однозначности решения задачи. При его нарушении однозначного решения обратная кинематическая задача не имеет. В го хе время в приложениях случай, когда скорость не монотонна, представляет особый интерес. Геофизиками высказываются довольно обоснованные к настоящему времени предположения о немонотонном поведении с глубиной скоростей сейсмических волн. Вопрос, есть ли такие понижения скорости, оказывается принципиальным при характе ристике вещества Земли. Есть основания считать из. модельных экс периментов, что в зонах понижения скорости вещество обладает осо быми свойствами. Поэтому естественен интерес, проявляемый к вопро су построения поля скоростей сейсмических волн при возможной не монотонности скорости. Однако, для этого, как уже говорилось при рассмотрении одномерной задачи, необходимо иметь источники, распо ложенные в зоне понижения скорости. Практически с такой ситуацией мы встречаемся в отдельных районах земного шара. Существуют такие районы, в которых естественные источники возмущений - землетрясе ния происходят на различных глубинах и даже заполняют целые прост ранственные области. Например, в районе Японские оотрова-Сахалин- Курильские острова-Камчатка такие землетрясения заполняют большой объем .уходящий в глубину до 600 км. Подобная же ситуация (правда, со значительно меньшей областью) характерна для района Средней Азии. Естественно поэтому рассмотреть постановку обратной кинема тической задачи с внутренними источниками, проведя некоторую мате матическую идеализацию. Заключается она в следующем: предполагает ся, что некоторая область 50 пространства я » ^ , ^ , . . . , ^ ) заполнена
187
источниками возмущении х°, от каждого из которых на поверхности
шара 5 , содержащего |
область |
50, измеряются времена пробега |
т(х'х% |
||
x°e5D, |
х'е S. |
Требуется найти скорость передачи сигналов |
|||
v(x)=-^~- |
для точек are© |
( с м . [ П , 1 4 ) ) . Своеобразие этой |
за |
||
дачи заключается в том, что не зная скорости передачи сигналов |
вну |
||||
три сферы 5 , мы не вправе считать |
известными координаты источни |
||||
ков х \ |
Более того, |
времена пробега |
т(х° х*) для каждого фикси |
рованного х° известны только с точностью до момента отсчета (на
чала |
возмущения), т . е . с точностью до аддитивной |
|
постоянной. |
В |
|||
предположении, что семейство функций Пх'х1), |
отвечающее различ |
||||||
ным х* допускает |
параметризацию (так что между |
х° и параметром |
|||||
u = (u1,...,uj |
существует взаимно однозначное и дифференцируемое |
||||||
соответствие),Ю.Е.Аяиконовым показано (см . [II] ) , |
что скорость |
v(x) |
|||||
находится по функции |
r(x°x J ) в области £> с точностью доконформ |
||||||
ных преобразований этой области. Отсюда следует, |
|
что в случае .ко |
|||||
гда |
область |
50 примыкает к поверхности сферы S, |
скорость внутри |
||||
£> |
определяется |
однозначно. |
|
|
|
||
|
6. О частных |
решениях обратной кинематической |
задачи. В целом |
||||
решение обратной кинематической задачи представляет довольно |
зна |
||||||
чительные трудности, |
как в формулировке удачного |
|
вычислительного |
||||
алгоритма, |
так и в его реализации. Поэтому немаловажную роль |
здесь |
играют различные решения обратной кинематической задачи. Существенные успехи в этом направления получены Ю.Е.Аниконовым. В рабо
те [12] им предложен довольно |
общий метод разыскания частных реше |
||||||||||||
ний. Идея метода |
|
заключается в том, что из некоторых дифференци |
|||||||||||
альных связей на заданную функцию |
т(х°х') |
следуют |
(при условии |
||||||||||
единственности решения задачи) соответствующие дифференциальные |
|||||||||||||
связи для искомой функции |
п(х). |
Например, для области ц>о |
про |
||||||||||
странства х,ц, |
если функция |
п(х,у) |
удовлетворяет условиям |
|
|||||||||
|
п |
= с + |
а |
п ч |
+Рп = |
0 ' |
|
|
П ( Х , 0 ) = СЦХ), |
( 3 9 ) |
|
||
то функция |
Tlx'x'), |
x°xie{y=o} |
|
удовлетворяет |
уравнению |
||||||||
|
|
|
|
+ |
a-[№(x°)-i£. |
+Va*(x<;-E^] +jbx=o. |
(40) |
||||||
|
|
|
|
|
|||||||||
При условии |
единственности |
решения задачи определения |
п(х,у) |
|
в |
||||||||
области |
у го |
по r(x°xJ ) |
верно и обратное: из уравнения (40) |
|
|||||||||
следует |
уравнение |
( |
3 9 ) . |
Этот факт можно использовать для отыска |
|||||||||
ния функции |
Шх,ц). |
Найдя из (40) функцию |
асх) и константы |
а |
|||||||||
и J J , можно, решив задачу |
( 3 |
9 ) |
(если <хФо |
), найти и функцию |
п[хф. |
188