Файл: Романов В.Г. Обратные задачи для дифференциальных уравнений спецкурс для студентов НГУ.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