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

Обратная задача рассеяния для систем дифференциальных уравнений с особенностью

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
767 Кб
Скачать

Глава 2. Решение обратной задачи рассеяния

Для заданной матрицы-функции v( ) 2 V, удовлетворяющей условию 2 теоремы 2.7, определим:

v](r) := (d (r)) 1v(r)d+(r);r2S0;

 

 

 

 

ˆ

]( ) :=

v

]

( )(

] ( )) 1

I

;

 

 

 

 

 

 

v

 

r

 

r

v0

r

 

 

 

V

]

(v;x;r)

 

 

 

 

]

 

 

]

1

 

1

;

 

:=F0 (x;r)v

(r)(v0(r))

(F0 (x;r))

 

 

ˆ

]( ;

;

r

) :=

 

( ;

]( )(

( ; )) 1

:

 

 

 

V

v

x

F0

 

x

 

r v

r F0

x

r

 

 

Заметим, что по построению матрица v](r) – нижнетреугольная

P-диагональная, нетривиальные диагональные блоки имеют вид

 

 

 

 

 

 

 

vk]+1;k

1

!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

0

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где vk]

+1;k(r) = dk+(r)=dk+1(r), k 2 I , а нетривиальные блоки матрицы

](r)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

:

 

 

 

 

 

 

 

 

 

 

 

vk]+1;k (v0] )k+1;k

 

 

 

 

 

 

 

 

 

 

 

 

0 !

 

 

 

 

 

Определим следующий зависящий от параметра x 0 оператор:

A]( )

( ) := ( +

f

)( )

(

f

)( ) ]( ;

r

) =

( )

(

f

)( ) ˆ ]( ;

r

);

x

f r

C

r

C

r V x

 

f r

C

r V x

 

где V ](x;r) := V ](v;x;r).

Лемма 2.13. Предположим, что v( ) 2 V удовлетворяет условию 2 теоремы 2.7. Тогда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ ]

 

nL1

 

 

 

N

для каждого x 2 [0;¥); более того,

 

1.

V (x; ) 2

 

 

 

 

= C0

(Sn) \ L2(S)

 

 

 

Vˆ ]( ; ) 2C [0;¥); n=1C0(

 

n)

\L2(S) , причём Vˆ ](0;r) 0.

 

 

 

S

 

2.

]

 

 

 

 

 

]L

 

 

 

 

]

(0) = Id.

 

A

(x) 2L(L2(S)), A ( ) 2C([0;¥);L(L2(S))), причём A

 

 

 

Доказательство. Зафиксируем произвольное n2f1;:::;Ng. При r2

2

S

n

d (r) = d

r k

, где d – некоторые ненулевые константы.

 

0k

k

 

 

 

k

 

 

 

 

 

 

 

 

 

 

Из свойств функций dk( ) следует

 

 

d+

 

 

 

 

 

]

 

 

(r)=dk+1(r) = r

 

k+1

k

 

 

 

 

 

 

vk+1;k(r) = dk

 

 

dk+1 + uˆk+1;k(r) ;

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

k

 

141

Игнатьев М. Ю. Обратная задача рассеяния

где k+1;k( ) 2C0(Sn) \ L2(Sn). Тогда

v]k+1;k(r) (v]0)k+1;k(r) = r k+1 k k+1;k(r);

или в матричной форме

](r) = r uˆ(r)r ;

где матрица-функция uˆ( ) 2C0(Sn) \ L2(Sn) P-диагональная, строго нижнетреугольная. Теперь из представления

 

 

ˆ ]( ;

r

) =

( ;

)

( )

(

( ;

r

)) 1

=

 

 

 

 

V x

F0 x

r r

uˆ r r

F0 x

 

 

 

 

 

 

= cˆ(rx)x lx (rx) uˆ(r)(rx) cˆ(rx)x lx 1 ;

где l – постоянная нижнетреугольная

обратимая матрица, следует, что

 

 

 

 

 

 

ˆ

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x;r) продолжается до непрерывной функции от (x;r) 2 [0;¥) Sn.

V

 

Поскольку uˆ(r) строго нижнетреугольная, из Re k < Re k+1, k = 1;n 1 вытекает, что

 

(rx) uˆ(r)(rx) ! 0; x lx ! l0 := diag(l11;:::;lnn)

 

 

 

 

 

ˆ

]

(0;r) 0.

 

при x ! 0. Таким образом, получаем V

 

 

 

Далее, имеем

 

 

 

 

 

 

Vˆ

](x;r) =

(x;r)r uˆ(r)r (

 

(x;r)) 1 = Y˜

(x;r)uˆ](x;r)Y˜ 1

(x;r);

 

F0

 

F0

 

 

0

0

 

(2.57)

где

](x;r) = W (rx)duˆ(r)d 1(W (rx)) 1:

Заметим, что в силу диагональности матриц W (rx) и d все ненулевые элементы ](x;r) находятся в позициях (k + 1;k), k 2I и равны

]

(x;r) =

dk+1

 

Wk+1(rx)

k+1;k(r):

(2.58)

 

dk

Wk(rx)

 

Напомним, что для k 2I имеем, в частности, Re(rRk+1) = Re(rRk) (r2 2Sn), что с учетом Re k+1 Re k означает, что дробь

Wk+1(rx)

Wk(rx)

 

 

 

 

 

 

допускает непрерывное ограниченное]

продолжение на

[0;¥) Sn. Итак,

из представления (2.58)]

вытекает uˆ ( ; )2C([0;¥);C0(Sn)\L2(Sn)). Тог-

ˆ

 

 

 

(Sn) \L2(Sn)).

да из (2.57) следует V ( ; ) 2C([0;¥);C0

2

 

 

 

142

Глава 2. Решение обратной задачи рассеяния

Зафиксируем теперь произвольную v( ) 2 V0, удовлетворяющую условию 2 теоремы 2.7. Определим

 

]

 

 

 

j(x;r) := (I + (Cp

(x; ))(r)) F0 (x;r);

где p](x;

) := (A](x)) 1Vˆ ](x; ). Отметим, что в силу пункта 2 леммы

 

 

]

(x)) 1 существует при всех

2.13 найдётся x0 > 0 такое, что оператор (A

x 2 [0;x0].

 

 

 

Лемма 2.14. Для всех x 2(0;x0] и r 2CnS справедливо равенство j(x;r) = y(x;r)d(r) (здесь, как и ранее в леммах 2.11, 2.12, y(x;r) =

= y(v;x;r)).

Доказательство. 1. Определим матрицу-функцию:

 

1

 

D(x;r) := Y0(x;r)d(r) F0 (x;r)

; r2CnS:

(2.59)

Заметим, что в силу свойств функции d( ) при r 2S0 существуют предельные значения D (x;r) и справедливо равенство

V ](x;r) = (D (x;r)) 1V (x;r)D+(x;r):

(2.60)

Переписав D(x;r) в виде

D(x;r) = Y0(x;r)d(r)(d0(r)) 1 (Y0(x;r)) 1 ;

воспользовавшись диагональностью и асимптотикой матрицы-функции d( ), получим

D(x;r) = I + o(1); r ! ¥;r2CnS:

(2.61)

Рассуждая как при доказательстве пункта 2 теоремы 2.5, из асимптотики (2.61) выводим

D(x;r) I = C(D+(x; ) D (x; )) (r);r2CnS;

(2.62)

а также равенство

 

 

 

 

 

 

 

 

 

 

 

1

w

 

dx

 

 

(D+(x)

 

D (x)) =

 

2pi

 

)(

x

)

 

(

 

 

 

S

r x

 

 

 

 

 

 

11

=r D( ) r D(r); r; 2CnS;

что эквивалентно представлению

CFr ( ) =

1

1

D(r); r; 2CnS;

 

D( )

 

r

r

143

Игнатьев М. Ю. Обратная задача рассеяния

где

Fr(x) := r 1 x(D+(x) D (x)); x2S:

Здесь и далее в данном доказательстве для краткости в списке аргументов опущен произвольный x 2 (0;x0].

Повторяя дальнейшие рассуждения из доказательства пункта 2 теоремы 2.5, из полученных равенств получаем следующее (где f ( )2L2(S) произвольна):

+ +

C C f D D (r) = (C f )(r)D(r) C f D (r): (2.63)

2. Введем функцию

h(r) := D+(r) D (r) + (C p)(r)(D+(r) D (r)) + p(r)D+(r):

Из (2.62), (2.63) следует

 

 

 

 

 

 

(Ch)(r) = D(r) I + (Cp)(r)D(r); r2CnS:

 

(2.64)

Из (2.64), переходя к соответствующим пределам, вычисляем

A]h = (D+ I) D V ] +V ] + (C+p)D+ (C p)D V ]:

Учитывая, что D V ] = V D+ и C+p + (C p)V = Ap = V

 

I, получаем

A]h = V ] I, откуда вытекает, что при x 2[0;x0] p] = h.

 

Теперь, используя снова (2.64), получаем

 

 

]

 

 

 

 

 

j(x;r) = (I + (Cp

(x; ))(r)) F0 (x;r) =

 

 

 

 

 

 

 

 

 

= (D(x;r) + (Cp)(r)D(x;r)) F (x;r) =

 

 

 

 

 

 

 

 

= (I + (Cp)(r))D(x;r) F (x;r) =

 

 

 

 

 

1

 

 

 

= (I + (Cp)(r))Y0(x;r)d(r) F0 (x;r)

 

F (x;r) = y(x;r)d(r):

 

 

2

 

 

 

 

Доказательство теоремы 2.7. Необходимость условий тео-

ремы фактически доказана ранее

в параграфах 2.1,

2.2, а также

1.3 (см. следствия 1.2, 1.3). Будем доказывать достаточность. Пусть фиксирована произвольная матрица-функция v( )2V, удовлетворяющая условиям 1–3.

1. Рассмотрим функцию y(x;r) = y(v;x;r). По построению данная функция – аналитическая по r в каждом из секторов Sn.

144

Глава 2. Решение обратной задачи рассеяния

Согласно лемме 2.12, как функция переменной x y(x;r) почти всюду на (0;¥) удовлетворяет дифференциальному уравнению

y0 = (rB + x 1A + q(x))y;

(2.65)

где q(x) := q(v;x) и в силу условия 3 q( ) 2Xp, p > 2.

Далее, в силу леммы 2.11 при каждом r2CnS имеет место асимп-

тотика

 

y(x;r) = (f(r) + o(1))exp(rxR); x ! ¥:

(2.66)

Кроме того, поскольку из леммы 2.9 следует включение p(v; ; ) 2 2 BC([0;¥);L2(S)), мы имеем, в частности,

kp(v;x; )kL2(S) = O(1)

при x ! 0. Отсюда с учетом леммы 2.10 вытекает следующая оценка для k-го столбца yk(x;r):

yk(x;r) = O(x k ); x ! 0:

(2.67)

Сопоставляя (2.65), (2.66), (2.67) c определением решений типа

Вейля, убеждаемся, что yk(x;r) = Yk(q;x;r) для всех x 2 (0;¥) и r 2

n

2 C n S таких, что Õ Dk(q;r) 6= 0.

k=1

2.Рассмотрим функцию j(x;r). В силу леммы 2.13 имеем p](0;r)

I, откуда следует асимптотика для k-го столбца

jk(x;r) = x k (hk + o(1)); x ! 0:

(2.68)

Далее, из леммы 2.14 имеем jk(x;r) = dk(r)yk(x;r), что дает продолжение j(x;r) на все x > 0, после которого будет выполняться оценка

jk(x;r) = O(exp(rxRk)); x ! ¥:

(2.69)

Кроме того, из соотношения jk(x;r) = dk(r)yk(x;r)

вытекает, что

jk(x;r) (для каждого k) удовлетворяют той же системе дифференци-

альных уравнений (2.65) (и с той же матрицей q( )), что и yk(x;r).

С учетом (2.68), (2.69) это означает, что jk(x;r) =Fk (q;x;r) для всех

n

x 2 (0;¥) и r 2 C n S таких, что Õ Dk(q;r) 6= 0.

k=1

3. Пусть n 2f1;:::;Ng произвольно. Для любого r 2Sn, такого что

n

Õ Dk(q;r) 6= 0, имеем

k=1

Dk(q;r) =

145

Игнатьев М. Ю. Обратная задача рассеяния

= Y1(q;x;r) ^ ^Yk 1(q;x;r) ^(r

k

 

n

 

 

Fk (q;x;r)) ^ ^(r

 

Fn (q;x;r)) ;

 

k

2 f

1;:::;n

g

 

 

 

 

 

где

 

 

также произвольно. Пользуясь результатами пунктов

1, 2, перепишем данное представление в следующем виде:

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

Dk(q;r) = Õ(r j dj(r))jfj:

 

(2.70)

j=k

Заметим, что по принципу аналитического продолжения представление (2.70) справедливо для всех r 2 Sn. В силу свойств функций dj( ) из

(2.70) следует, что Dk(q;r) 6= 0 для всех r 2Sn, n = 1;N, k = 1;n. Таким образом, q( ) 2 G0p.

4. При п.в. r 2 S0 существуют предельные значения y (x;r), равные по построению (I +(C p(x; ))(r))Y0 (x;r). Уравнение A(x)p(x; ) =

= ˆ ( ; ), эквивалентное соотношению

V x

I + (C+p(x; ))(r) = (I + (C p(x; ))(r))V (x;r)

(для п.в. r 2 S0), влечет равенство y+(x;r) = y (x;r)v(r) для п.в. r 2 2S0. Но по доказанному в пункте 1 y(x;r) = Y(q;x;r). Таким образом, для п.в. r 2 S0 имеем

Y+(q;x;r) = Y (q;x;r)v(r):

(2.71)

Из результатов парагр. 3 известно, что для q( ) 2 G0p

(а это свойство

установлено выше в пункте 3) предельные значения Y (q;x;r) существуют для всех r 2S0 и представляют собой непрерывные функции от r 2S0. Следовательно, соотношение (2.71) справедливо для всех r 2S0, что означает v(r) = v(q;r).

2

Теорема 2.8 показывает, что в случае v( ) 2V00 выполнение наиболее трудно проверяемого условия 3 следует из выполнения условия 1.

Теорема 2.8. Пусть v( )2V00 такова, что выполнено условие 1 теоремы 2.7. Тогда q(v;x) непрерывна по x 2 [0;¥) и справедлива оценка q(v;x) = O(x m) при x ! ¥, где m 0 произвольно.

Доказательство. В условиях теоремы формула восстановления принимает вид (см. также доказательство леммы 2.12, пункт 1)

q(v;x) = 1 w [B;p(v;x;r)]dr:

2pi

S

Введем функции

P (x;r) := E (x;r)(I + wˆ (r))(E (x;r)) 1;

146

Глава 2. Решение обратной задачи рассеяния

где матрицы-функции wˆ ( ) – те же, что в лемме 2.7 (функции P ( ; ) теперь строятся иначе!). Построим по этим функциям зависящий от параметра x 2 [0;¥) оператор B(x):

(B(x) f )(r) =

= (C+( f ( )P+ 1(x; )))(r)P+(x;r) (C ( f ( )P+ 1(x; )))(r)P (x;r):

Тогда для p(x) = p(v;x) имеем представление

 

 

 

ˆ

(x);

(2.72)

p(x) + H(x)p(x) = B(x)V

где

 

 

 

 

H( ) := B( )A( )

Id

:

 

x

x x

 

 

Заметив, что

P 1(x;r) := E (x;r)(I wˆ (r))(E (x;r)) 1;

и повторяя вычисления, используемые при доказательстве леммы 2.7, получим для произвольной j 2 L2(S)

BAj = j+ C j C+(P 1 I) C (P+ 1 I) (P+ P )+

+(C j) C+(I P 1) +C (P+ 1 I) (P+ P )+

 

 

 

 

 

+C

(C j)V1 P C+

(C j)V1

P+:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x)

 

 

 

 

 

 

 

 

 

 

 

Таким образом, для оператора H

 

 

 

справедливо представление

 

H( ) =

C

(

ˆ)(

 

 

)

 

(

 

 

) ˆ(

 

 

 

)+

 

 

 

 

 

 

x j

 

 

jP P+ P

 

 

C

 

j

P P+ P

 

 

 

 

 

 

 

 

где

 

 

 

 

 

+C (C j)V1

(P P+) (C j)V1P+;

(2.73)

 

 

 

 

 

 

ˆ

 

 

 

 

ˆ

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(x;r) := P+(x;r) P (x;r);

 

 

 

 

 

ˆ

:=

C

+( 1

I

);

ˆ

:=

C

( 1

I

);

V1

=

1

 

1:

 

P+

 

 

P

 

P

 

 

P+

 

 

 

V P+

 

P

 

Далее, имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B ˆ

=

C

+( ˆ 1)(

 

 

 

)

C

( ˆ

1)(

P

I

) + ˆ

1

:

(2.74)

V

 

V P+

P+ I

 

 

V P+

V P+

 

 

Исходя

из представлений (2.72), (2.73), (2.74), запишем q(x) =

= q(v;x) в виде суммы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q(x) = [B;q0(x)] + [B;q1(x)] + [B;q2(x)];

 

 

 

(2.75)

147

Игнатьев М. Ю. Обратная задача рассеяния

где

 

 

 

 

 

q0(x) =

1

w drVˆ (x;r)P+ 1(x;r);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2pi

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q1

(x) =

 

1

 

 

w dr(H(x)p(x;

))(r);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2pi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

q2

(x) =

1

w d C+

Vˆ

(x; )P+ 1(x; ) ( )(P+(x; )

 

I)

 

 

 

 

2pi S

 

 

 

 

 

 

 

 

 

 

 

 

1

w d C

Vˆ

(x; )P+ 1

(x; ) ( )(P (x; )

 

I):

 

 

 

2pi S

 

 

 

 

 

 

 

 

 

Последнее представление преобразуем к виду

(2.76)

(2.77)

q2(x) =

 

1

w d (Vˆ P+ 1)(x; )Pˆ

(x; )+

 

 

 

 

 

 

 

 

 

 

2pi

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w d (Vˆ P+ 1)(x; )Pˆ+(x; ):

 

 

 

 

 

 

 

 

+

 

1

(2.78)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2pi

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя представление (2.41) и оценку (2.42), получаем

 

 

ˆ

1

 

( ; ) = +( ; ) ( )( +( ; )) 1

 

(2.79)

 

 

E

(x;r)wˆ (r)(E (x;r)) 1

+ O(exp( tjrjx)); t > 0

V P+

 

x r

 

E x r wˆ+ r E x r

 

 

 

равномерно по (x;r) : jrxj > 1, т. е., в частности, по r 2 supp(v( )v0( )) при всех достаточно больших x. Отметим, что из полученного результата вытекает, в частности, оценка

( ;

r

) = ˆ ( ;

)

1

( ;

r

) +

1

( ;

r

)

1

( ;

r

) =

V1 x

V x

r

P+

x

 

P+

x

 

P

x

 

= O(exp( tjrjx)); t > 0 (2.80)

равномерно по r2supp(v( ) v0( )) при всех достаточно больших x. Рассмотрим теперь подробнее функции E (x;r)wˆ (r)(E (x;r)) 1.

Пусть r 2 Sn, где n 2 f1;:::;ng, произвольно. Тогда справедливо представление

exp(rxR )wˆ (r)exp( rxR ) = å exp(rx(Rk Rk+1))(wˆ (r))k;k+1ek;k+1;

k2I

148

Глава 2. Решение обратной задачи рассеяния

где e jk – матрица, единственный ненулевой элемент которой равен 1 и находится в позиции ( j;k). Пользуясь асимптотиками функций E (x;r), после очевидных преобразований приходим к выражению

E (x;r)wˆ (r)(E (x;r)) 1 =

= å exp(rx(Rk

 

 

 

 

m 1

 

 

 

Rk+1))(wˆ

 

(r))k;k+1

å

(rx) sMks + O x m

;

k2I

 

 

s=0

 

 

в котором Mks некоторые постоянные матрицы и O( ) равномерно по r2supp(v( ) v0( )) при всех достаточно больших x. Полученное представление можно, в свою очередь, преобразовать к виду

E (x;r)wˆ (r)(E (x;r)) 1 =

m 1

 

 

 

= å å x s exp(rx(Rk

 

Rk+1)) fks (r) + O x m

; (2.81)

k2I s=0

 

 

где fks ( ) – некоторые бесконечно дифференцируемые матрицы-функ- ции.

Поскольку для каждого n =

 

 

 

 

 

 

1;N имеем

 

 

w drexp(rx(Rk

 

Rk+1)) fks (r) = O

x m

 

Sn

 

 

 

 

при x ! ¥ (интеграл в левой части есть фактически преобразование Фурье – Лапласа от финитной бесконечно дифференцируемой функции), из (2.81) и (2.79) следует оценка

q0(x) = O x m :

Далее, пользуясь леммой 2.6 (см. доказательство леммы 2.7), получаем оценки

где

C Fks (x; ) L¥(S) = O x m

;

 

 

 

 

откуда вытекает

Fks (x;r) = exp(rx(Rk Rk+1)) fks (r);

 

kPˆ (x; )kL¥(S) = O x m :

 

 

 

(2.82)

В силу представления (2.78) из (2.82) и (2.79) следует оценка

q2(x) = O x m :

149

Игнатьев М. Ю. Обратная задача рассеяния

Таким образом, для завершения доказательства теоремы осталось рассмотреть функцию q1( ). Прежде всего заметим, что в силу представления (2.73) и оценок (2.80), (2.82)

 

 

(Id + H(x)) 1

kH(x)k = O

x m :

(2.83)

Но тогда

= O(1) и из (2.72), (2.74) следует оценка

kp(x; )kL2(S) = O(1).

Теперь из представления (2.77) и оценки (2.83)

 

 

функции

 

 

(см. (2.73)) вытекает

с учетом

финитности по

r

(H(x)p(x; ))(r)

 

 

 

 

 

q1(x) = O x m :

2

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]