Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Романов В.Г. Обратные задачи для дифференциальных уравнений спецкурс для студентов НГУ

.pdf
Скачиваний:
43
Добавлен:
23.10.2023
Размер:
8.33 Mб
Скачать

 

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)\ < с-ЩЧ + «ф№ •

\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.

 

Соседние файлы в папке книги из ГПНТБ