Файл: Романов В.Г. Обратные задачи для дифференциальных уравнений спецкурс для студентов НГУ.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.06.2024
Просмотров: 215
Скачиваний: 1
да по ранее доказанному (см. § 5 главы П) вектор
представляет собой единичный вектор касательной в точке х к лучу Ц(х'х). Поэтому равенство (5) можно вдоль луча Г„(х°х) запи сать в виде
|
|
T J a ^ x ) |
= |
rt±ix). |
|
|
|
|
Интегрируя |
его вдоль луча |
Ц(х°х) |
от точки х° |
до точки |
х\ |
|||
находим |
|
|
|
|
|
|
|
|
|
zL(x° х 1 ) = j |
nt№) |
ds. |
|
(7) |
|||
|
|
Cfx;x') |
|
|
|
|
|
|
В этой формуле через ds |
обозначен элемент |
евклидовой длины ду |
||||||
ги кривой |
£(х°х% Лучи |
£(х°х') |
зависят |
только |
от |
п„(х) |
и мо |
|
гут считаться поэтому известными. Линеаризированная |
постановка |
обратной кинематической задачи приводится таким образом к следую
щей задаче интегрально;! геометрии: от функции f i j x ; |
известны |
|||
интегралы по лучам |
[^(х'х1), |
зависящим от функции |
njx) |
и свя |
зывающим пары точек |
сс° х ' е 5 ; |
требуется найти через эти интег |
||
ралы функцию n.t(x). |
Исследование этой задачи мы можем провести, |
опираясь на изложенные в предыдущем параграфе результаты. |
|
||
Рассмотрим вначале частный случай, |
когда исследование |
прово |
|
дится наиболее просто. А именно, пусть |
функция njx) |
зависит |
только от одной координаты. Выделим ее, обозначив через у , а за
остальными (n-i) |
координатой сохраним прежнее |
обозначение |
х . |
||
Таким образом, |
теперь л = п(х,у), |
a |
na~njyi. |
Лучи Га в этом |
|
случае являются плоскими^кривыми. Каждый луч £ |
лежит в плоскости |
||||
размерности 2, параллельной оси ц. |
Поэтому в каждом сечении |
обла |
|||
сти £ ) такой плоскостью мы получаем |
задачу интегральной геомет |
рии на плоскости, В связи с этим мы можем в дальнейшем считать х одномерной координатой й рассмотреть задачу в двумерной плоскос
ти х,у. |
Пусть |
"гО - |
полуплоскость у&о , 5 - ее граница, |
то есть |
||
прямая у=о . Лучи Ц |
соединяют пары точек границы 5. |
Чтобы |
||||
найти функцию |
njx,у) |
в области 50 , нужно, конечно, чтобы лу |
||||
чи Ц |
лежали в области 50 . Ранее нами было выяснено условие,при |
|||||
котором лучи |
J7 представляют собой вид дуг, лежащих в полуплос |
|||||
кости цго. |
Условие |
это сводится к тому, что функция |
па<ц) долж |
|||
на быть монотонно убывающей функцией у . Если функция |
njy) |
име- |
170
ет по у |
отрицательную производную на некотором отрезке [о,ИЗ,то |
||||||
в полосе |
o<,ys:H |
пространства |
х,у |
из каждой |
точки |
[$,-ф мож |
|
но провести луч |
I2($,vK |
имеющий в этой точке |
свою вершину и |
||||
опирающийся концами на ось у=о. |
Именно такое свойство |
мы требо |
вали у семейства кривых в задачах интегральной геометрии. Пока
жем, что условие |
n'jyxo |
при у^[о,НЗ в совокупности с усло |
||||||||
вием ограниченности |
на этом же отрезке |
п"(у) |
являются |
достаточ |
||||||
ными, чтобы |
задача |
|
(7) |
имела |
в полосе |
<Он=[(х,у)-. |
оауаН} |
|||
единственное |
решение. |
Заметим прежде всего, что лучи |
О*,??) ин |
|||||||
вариантны относительно |
параллельного переноса вдоль оси у=о. Это |
|||||||||
является следствием |
того |
факта, что функция п . |
зависит |
только |
||||||
от ц и не зависит |
от х . Проверим, что дифференциальные |
свойст |
||||||||
ва семейства |
кривых |
Е($,т?) |
удовлетворяют требованиям |
теоремы1 |
предыдущего параграфа. Воспользуемся для этого уравнением сейсми
ческого луча в одномерной среде. Из формул |
(28) , (29) § 5 главы П |
в принятых нами здесь обозначениях следует, |
что уравнение луча |
имеет вид: |
|
|
|
|
j |
-fillip |
|
|
|
|
Из |
сравнения этой формулы с формулой (15) предыдущего параграфа |
|||||||
следует, что в обозначениях § 3 , |
|
|
|
|
|
|||
|
|
|
7•с |
|
|
|
|
|
Используя для |
frpy |
прежнее обозначение |
р , находим |
|
||||
|
|
|
|
|
|
|
|
О ) |
|
|
|
П-Р' |
|
|
|
(p.(p,i2) |
|
Нам нужно изучить дифференциальные |
свойства |
функции |
в |
|||||
области G={(p7Ti)-- |
О б р = а ^ «//Т } . |
Условие, что п'о(г)<0 |
||||||
при |
ze[o,H] |
приводит к тому, что подинтегральное выражение |
в |
|||||
формуле (9) имеет интегрируемую особенность |
f^-z)"^ |
и интеграл |
||||||
(9) |
всегда конечен. Из формулы (9) |
сразу |
следует, что в области |
|||||
G: |
|
|
|
|
|
|
|
|
|
|
%(Р,Ч)>°, |
<pj(0,4)=o. |
|
|
|||
Дифференцируя равенство (9) по переменной |
р , находим |
|
|
171
Из формулы (10)следует, |
|
что |
t ^ > o |
внутри области Or и непрерыв- |
||||||||||||||||
но дифференцируема по |
р |
и 7? . Особо надо |
исследовать |
позедение |
||||||||||||||||
этих производных при |
р —- о. |
Легко,однако, |
показать, |
используя |
||||||||||||||||
условие |
п'/укр, |
|
уе[о,Н], |
|
|
и условие |
ограниченности |
п'су), |
||||||||||||
что производные -|^, |
|
|
|
|
эрЗт? |
|
существуют |
при р^-о |
и огра |
|||||||||||
ничены. Действительно, используя формулу Тейлора, |
находим |
|
|
|||||||||||||||||
Здесь 0(р") |
означает величину |
порядка р" |
при |
р — о . |
Подстав |
|||||||||||||||
ляя |
это |
разложение |
|
в формулу |
( 1 0 ) , |
получаем |
|
|
|
|
|
|
||||||||
|
|
|
|
|
э р |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отсюда легко следует, что предельные значения при |
р—о |
указанных |
||||||||||||||||||
производных конечны, и в частности |
|
|
>о. |
Таким |
образом, |
|||||||||||||||
функции |
(р^1р,т1) |
удовлетворяют |
всем |
требованиям |
|
теоремы I |
§ |
3 , |
||||||||||||
Нам остается проверить, |
что |
весовая функция, порождаемая |
за |
|||||||||||||||||
дачей ( 7 ) , также удовлетворяет |
условия//! теоремы, |
|
h данном |
случае |
||||||||||||||||
р= |
|
Используя формулу |
( 8 ) , |
можно налти янное |
|
вгрокение |
д л я ^ : |
|||||||||||||
|
|
|
^ = J 7 7 ] X F = |
|
* ш . |
|
|
|
|
|
|
i n ) |
||||||||
Так |
как |
у=у/хr -$, |
|
тр, |
|
тоv |
весовая функция |
j> =j)(x-$,•>?). |
|
Из фор- |
||||||||||
|
|
|
V |
|
|
\x l |
|
|
|
п.iv) |
|
|
|
|
|
|
|
|||
мулы |
( i l ) следует, |
|
что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
filO, |
V) |
= i. |
|
|
|
|
|
|
|
|
Дифференцируя |
соотношение |
( I I ) |
по |
зс |
и по |
, находим |
|
|
|
|
||||||||||
|
|
э х |
|
п.гч» |
|
э х |
|
п„ г (^ |
|
" ° , 7 Z ' > |
|
|
|
|
||||||
|
|
Э£_ |
|
|
njyjnjfy) |
|
|
n'jy) |
Э£ |
|
|
|
|
|
(12) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
Производную § ^ мы можем вычислять, используя равенство
x = i + r-i>* (ft Щ-у, |
v)- |
172 |
|
Дифференцируя его по переменной , находим
|
|
|
т. |
2Е. |
JLSL + |
Ж |
|
= 0 |
|
|
|
|
|
||
Отсюда |
|
Эр |
Эу |
Эт? |
Эт? |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
эу |
_ |
р |
эт? |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
э р |
|
|
|
|
|
|
|
Из непрерывности |
^ |
Р |
и условия |
- 1 ^ > о следует |
непрерыв- |
||||||||||
|
|
|
|
|
|
|
эр |
|
|
|
|
|
|
||
ность |
а я |
Тогда из формул (12) заключаем о непрерывности |
част |
||||||||||||
ных производных весовой |
функции. |
|
|
|
|
|
|
|
|
||||||
Итак, |
семейство |
кривых £(£,?z) |
и весовых функций |
удовлет |
|||||||||||
воряет |
условиям теоремы I |
§ 3 . Отсюда следует |
теорема. |
|
|
|
|||||||||
Т е о р е м а I . Пусть-функция |
ri(x, у) |
представима |
в виде |
||||||||||||
|
|
|
П(х,у) |
= псЩ) |
+ |
п1(х,у), |
|
|
|
|
|
||||
причем функция njy) |
задана |
и удовлетворяет |
условиям |
|
|
|
|||||||||
|
п*(у»о, |
<(yl<o, |
|
|
|
K ' t y j / < < = ~ , |
|
0*</бН, |
|||||||
а функция |
nt/x,y) |
мала по сравнению с |
i\iy) |
и финитна. Тогда в |
|||||||||||
области |
^) = {(х,у) |
.•-•=>= <ос<°°, ойу&Н} |
функция п(х,у),в |
ли |
|||||||||||
неаризированной постановке, однозначно определяется временами |
|||||||||||||||
t(x°xi) |
|
для любых пар точек |
эс°, |
x'eS. |
|
|
|
|
|
||||||
В формулировке этой теоремы есть некоторый дефект. Дело |
в |
||||||||||||||
том, что при исследовании |
задачи интегральной |
геометрии мы пред |
|||||||||||||
полагали, |
что есть от функции njx,y) |
интегралы для любых |
|
|
|||||||||||
(£,??j€50. В то же время между |
точками |
x°xi |
и параметрами |
£,7?, от |
|||||||||||
вечающим концам кривой |
ГЦ,У]), |
может не быть взаимно |
однознач |
||||||||||||
ного соответствия. Это связано с возможной неоднозначностью |
со |
||||||||||||||
ответствия паре точек |
х°х*е |
|
{у = о} |
экстремали |
Га[х°х*), |
их со |
|||||||||
единяющей. Ранее мы об этом уже говорили при исследовании |
одно |
||||||||||||||
мерной обратной задачи. Взаимно однозначное соответствие будет |
|||||||||||||||
обеспечено, если считать, что функция |
п0ш} |
удовлетворяет усло- |
|||||||||||||
вию |
j i |
|
на |
|
|
|
|
" |
|
|
|
|
|
|
|
2 _ In n„ty)^o |
[о, И]. В противном случае для справед- |
||||||||||||||
|
Ч |
|
|
|
|
|
|
|
|
|
|
|
|
%[х° х1 ), |
|
ливости теоремы мы должны считать, |
что известны времена |
||||||||||||||
отвечающие всевозможным экстремалям, соединяющим пару точек |
|
||||||||||||||
Отметим еще интересный |
случай, |
когда |
скорость |
ЧЩр^хгт воз- |
173
растает |
по линейному |
закону |
|
|
|
|
|
|
|
Vjy) |
= а+ |
, |
а>о, |
&>о. |
|
||
В этом случае лучи, соединяющие пару точек |
х°, |
оси ^=о,пред- |
||||||
ставляют собой дуги окружностей, центр которых |
находится на пря- |
|||||||
мой у=—^. Тогда, если функцию |
n.tix,y) |
доопределить в области |
||||||
t |
|
|
|
|
i |
|
|
|
-— ^у<о |
нулем и продолжить в область |
цё—— |
четным, |
относи |
||||
тельно прямой у=~-£ |
образом, го задача (7) |
приводится к |
извест |
|||||
ной задаче определения |
функции через |
ее интегралы по окружностям, |
||||||
центр которых пробегает |
прямую |
y=--g. |
Впервые внимание на этот |
факт был обращен в работе [ 8 3 ] .
Аналогично изложенному выше, легко исследуется случай, когда
функция п,(х) |
зависит только от расстояния до фиксированной |
точ |
||||||
ки х |
пространства. В этом случае |
разумно рассмотреть оферт с |
цен |
|||||
тром в этой точке. Пусть |
т = /х-х|. |
Тогда, |
если функция |
п„(г) |
име |
|||
ет конечную вторую производную и |
^ [ ^ Я ь ' " 1 ' ] > ° , т о |
сочно- |
||||||
непрерывная функция ntlx) |
-однозначно находится внутри фиксирован |
|||||||
ной сферы с центром в х |
по r . J x ' x * ) , где |
ос* х*—любая |
пара |
то |
||||
чек |
этой сферы. Доказательство этого |
факта |
содержится в § 5 гла |
|||||
вы П работы |
[ 1 2 8 ] . |
|
|
|
|
|
|
|
Отметим, что эти два довольно |
простых для исследования случая ' |
имеют существенное прикладное значение для геофизики. Действитель но, к настоящему времени существует достаточно много одномерных скоростных разрезов Земли, либо вдоль радиуса Земли (глобальных),
либо как функция глубины |
в г ".инок районе (.региональных). |
Поэтому |
распределение i\,(z) или |
\:ожет считаться известным, |
Естест |
венно тогда использовать |
линеаризированную 'постановку обратно:; |
задачи для уточнения распределения поля скоростей сейсмических волн в Земле, в частности для наделения горизонтальных неоднородностей. Это тем более разумно, что есть все основания предпола гать, что отклонения в скоростях волн от радиального распределе ния достаточно малы. На основании полученной задачи интегральной
геометрии могло сконструировать |
вычислительные |
алгоритмы опреде |
ления малых добавок к скоростям |
волн (см. [ 4 ] , |
[ 5 ] , [ 1 2 8 ] ) . |
Перейдем |
теперь |
к исследованию случая, когда функция' njx) |
|
существенно |
зависит |
от всех переменных. Нам предстоит |
" отоа слу |
чае выяснить поведение лучей п среде, характеризуемой |
скоростным |
||
|
|
174 |
|