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

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

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

Глава 1. Решения типа Вейля

Доказательство. Пусть k 2 f2;:::;n 1g (случаи k = 1, k = n требуют незначительных изменений). Имеем

Dk = jFk 1 ^Tkj = r k jFk 1^ Fk ^Tk+1j =

= r k dkjFk 1 ^Yk ^Tk+1j = r k dkDk+1:

 

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

 

dk(q;r) = r k

Dk(q;r)

 

(1.122)

 

Dk+1(q;r)

 

(при k = 1 и k = n следует положить D1 = 1, Dn+1 = 1 соответственно). Используя утверждение 1.2, выводим из (1.122)

dk(q;r) = r k

Dk(q;r)

 

= r k

D0k + o(1)

 

= dk(0;r)(1 + o(1))

Dk+1(q;r)

 

 

D0;k+1 + o(1)

 

при r ! ¥. Аналогично получается и второе утверждение следствия.

2

˜ : 1

Всюду далее будем обозначать Y(q;x;r) = Y(q;x;r)(W (rx)) .

Теорема 1.10. Пусть p > 2. Тогда для элементов матрицы

: ˜ 1 ˜

b(q;x;r) = (Y0(x;r)) Y(q;x;r)

справедливы представления

bjk(q;x;r) =

dj;k + djk(q;x;r)

;

1 + dk(q;x;r)

где dj;k – символ Кронекера, djk;dk 2C(Xp;BC([0;¥);C0(S)) и для любого луча l = fr = zt;t 2[0;¥)g, z2S nf0g ограничения djkjl ; dkjl принадлежат

C(Xp;BC([0;¥);H (l))).

Доказательство. Определение матрицы b(q;x;r) эквивалентно представлениям

˜

 

n

 

˜

 

 

 

 

Yk(q;x;r) = åbjk(q;x;r)Y0 j(x;r)

 

 

j=1

 

 

для каждого фиксированного k =

 

:

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

1;n

в соотношения

 

 

 

 

 

˜

˜

˜

˜

˜

Fk 1

^Yk = Fk; Yk ^Tk = 0;

91

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

получаем

 

 

 

 

 

8 j=1

 

 

 

^

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

n

 

 

˜

 

 

 

 

˜

 

 

 

˜

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

å bjkFk

 

 

Y0 j = Fk;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

˜

 

 

 

 

˜

= 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

å bjkY0 j

^

Tk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:

вытекаетj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда,

 

 

как

 

следствие,>

 

следующая

СЛАУ

относительно

fbjkgkj=1

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(здесь и далее произвольное k фиксировано):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

åmi jbjk = ui;

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.123)

mi j =

 

 

 

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

 

 

 

 

 

 

 

T˜k

 

 

 

 

 

 

 

 

˜

 

 

 

˜

 

˜

 

Y˜

 

˜

 

 

 

^

˜

Y˜

 

 

 

 

˜

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0;k

 

 

 

 

; i = 1;k 1;

 

Y˜ 01^

Y0 j

^Y˜ 0;i ^1

 

0;i 1

 

 

 

^1 Y˜^0; j

 

 

mi j =

 

Fk 1

 

 

Y0k

^Y0;i 1

 

 

Y0;i+1

 

 

Y0n

;i

= k;n;

(1.124)

 

 

 

 

 

 

^ ^

^

 

 

^ ^

 

 

^

 

 

^

 

 

 

 

 

 

 

 

 

 

 

 

 

˜

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

˜

 

 

˜

 

 

 

 

˜

 

 

 

 

 

 

 

˜

 

 

 

 

 

 

 

 

 

 

(1.125)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ui

= 0; i = 1;k

 

1:

 

 

 

 

 

 

 

 

 

 

(1.127)

ui =

 

Fk ^Y0k ^ ^Y0;i 1 ^Y0;i+1

^ ^Y0n

 

; i = k;n

(1.126)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

˜

 

 

 

 

 

 

 

 

С учётом теоремы 1.6, леммы 1.2 и свойств матрицы

 

( ; ) из (1.123)

Y0

следует,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mi j(q;x;r) =

0 Wi(rx)

+ mˆ i j(q;x;r)

 

 

1

 

 

mi j

 

 

 

 

 

 

 

;

Wj(rx)

!

n

r

 

 

 

 

 

 

 

 

x)

 

 

 

 

 

 

W

 

(

ui(q;x;r) =

0 Wi(rx)

+ uˆi(q;x;r)

 

 

1

 

 

ui

 

 

 

 

 

 

 

;

Wk(rx)

!

n

r

 

 

 

 

 

 

 

 

x)

 

 

 

 

 

 

W

 

(

где скалярные функции

i j = mˆ i j(q;x;r), i = uˆi(q;x;r) таковы, что

i j;uˆi 2C(Xp;BC([0;¥);C0(S)) и для любого луча l = fr = zt;t 2[0;¥)g,

z 2S nf0g ограничения i jjl ; uˆijl 2C(Xp;BC([0;¥);H (l))).

Поскольку скалярные функции с такими свойствами, как у i j, i образуют алгебру относительно поточечного умножения, применяя к системе (1.123) правило Крамера и учитывая, что

m0i j = di; j( 1)i k+1D0k; i = 1;k 1;

m0i j = di; j( 1)i kjfj; i = k;n;

u0i = di;kjfj

(di; j – символ Кронекера), получаем требуемые представления.

2

92

Глава 1. Решения типа Вейля

Полученный результат допускает естественное уточнение в том случае, если характеристические функции не обращаются в нуль на некотором (заданном) множестве.

Определение 1.2. Пусть L – некоторое (вообще говоря, неограниченное) подмножество сектора S. Будем говорить, что q 2Xp принадле-

n

жит классу G0p(L), если Õ Dk(r) 6= 0 для всех r 2L.

k=1

Лемма 1.12. В представлении

ˆ

Yk(q;x;r) = Y0k(x;r) +Wk(rx)Yk(q;x;r)

для любого отрезка [0;T ];0 < T < ¥ и любого замкнутого множества

 

 

ˆ

p

(L);C([0;T ];C0(L))).

 

 

L S справедливо Yk 2C(G0

Для доказательства леммы нам понадобится следующее обобщение результата о непрерывной зависимости решений невырожденных СЛАУ от параметров.

Лемма 1.13. Рассмотрим СЛАУ Ax = b, где A N n, n N – матрица, b 2 CN . Пусть A = A(p;t), b = b(p;t), p 2 P, t 2 T , где P;T – метрические пространства, причём T – компакт. Пусть СЛАУ зависит от параметров непрерывно в том смысле, что A(p; ) 2C(T ), b(p; ) 2C(T ) для любого фиксированного p2P и A2C(P;C(T )), b2C(P;C(T )). Кроме того предположим, что СЛАУ однозначно разрешима при любых p 2P, t 2T . Обозначим её решение x(p;t). Тогда x 2C(P;C(T )).

Доказательство. 1. Случай n = N очевиден. Далее полагаем N > n. 2. Зафиксируем произвольное p 2 P. Пусть an(p;t)(2CN ) обозначает n-й столбец матрицы A(p;t). Из однозначной разрешимости СЛАУ следует, что a1(p ;t)^ ^an(p ;t) 6= 0 для всех t. Отсюда для каждого t0 найдутся векторы a(n0+)1;:::;a(N0) d0 > 0, такие что a1(p ;t)^ ^an(p ;t)^

n+1 ^ ^

 

N

6 0

при всех

t :

 

t

m0

)

d0

. В силу компактности T

a(0)

 

a(0)

=

 

dist( ;t

 

 

Tj

– компакт;

 

 

 

 

 

 

 

 

jS1

 

 

 

 

 

 

 

 

 

существует конечное разбиение T =

=

Tj, такое что:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

– для каждого

j =

 

найдутся

векторы an(+j)

1;:::;aN( j),

такие что

1;m

a

(p ;t)

^ ^

a

(p ;t)

^

a( j)

 

^ ^

a( j)

= 0 при всех t

2

T

.

1

 

 

n

 

 

 

n+1

 

N

6

 

 

 

j

 

Определим матрицу A˜( j)(t) := (a1(p ;t);:::;an(p ;t);an(

+j)

1;:::;aN( j)) и век-

тор j(t) := (x1(p ;t);:::;xn(p ;t);0;:::;0)T 2CN . Тогда при каждом t 2Tj

˜( j)( ) есть единственное решение СЛАУ ˜( j)( ) ( j)( ) = ( ; ). В силу x t A t x t b p t

пункта 1 имеем ( j) 2C(Tj), j = 1;m, следовательно, x(p ; ) 2C(T ).

93

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

3.Пусть pk ! p в P. Покажем, что для каждого j = 1;m x(pk;t) !

!x(p ;t) равномерно по t 2 Tj.

Определим матрицу

˜

( j)

(p;t)

 

 

 

 

 

 

 

 

 

 

 

( j)

( j)

)

A

 

:= (a1(p;t);:::;an(p;t);an+1;:::;aN

и вектор

j

(

;

 

) := (

x1

(

;

t

);:::;

xn

(

 

;

t

);

0

;:::;

0

)T

N . Тогда при каждом

t 2 Tj,

 

 

 

p

t( j)

 

 

p

 

 

 

 

p

 

 

˜

( j2) C

 

( j)

(p;t) = b(p;t).

p 2 P x˜

 

(p;t)

 

есть решение СЛАУ

 

A

(p;t)x˜

 

Поскольку

 

функция

 

 

 

˜

( j)

(p;t)

 

 

принадлежит

классу C(P;C(Tj))

 

 

detA

 

 

 

t Tj

j

 

 

 

 

 

 

 

j

>

0, то найдется окрестность

P

j

точки

p , такая

и inf

 

detA˜( j)(p ;t)

 

2

 

˜( j)

(p;t) 6= 0 для всех p 2 Pj, t 2 Tj. Но тогда в силу пункта 1

 

 

что detA

 

 

имеем ( j) 2C(Pj;C(Tj)), что даёт требуемую равномерную сходимость x(pk;t) ! x(p ;t).

2

Доказательство леммы 1.12. 1. Прежде всего заметим, что для

любого q

2

Gp

(L) функция Y˜

k

(q;x;r)

:= (W (rx)) 1Y

(q;x;r) допускает

 

0

 

 

 

k

k

 

 

непрерывное продолжение на [0;¥) L и всюду на указанном множе-

 

 

 

 

 

˜

˜

˜

˜

˜

˜

стве имеют место разложения Fk = Y1 ^ ^Yk, Tk = Fk ^ ^Fn, k =

=1;n (обоснование этого утверждения дословно повторяет доказательство теоремы 1.9). В силу указанной разложимости тензоров Fk, Tk, k =

=1;n условие Dk(r) 6= 0, r2L гарантирует однозначную разрешимость

для всех (x;r) 2 [0;¥) L системы

 

 

 

 

˜

˜

˜

˜

 

˜

 

Fk 1

^Yk = Fk; Yk ^Tk = 0:

 

Применяя лемму 1.13, заключаем, что

˜

p

(L);C([0;T ] K)) для

Yk 2C(G0

любого отрезка [0;T ], 0 < T < ¥ и любого pкомпакта K L. Заметим,

что такое включение эквивалентно

˜

(L);C([0;T ];C(K))). В си-

Yk 2C(G0

 

 

 

 

 

ˆ

 

 

лу произвольности компакта K отсюда следует также, что Yk(q; ; ) 2

2C([0;T ] L).

 

 

 

 

 

ˆ

2. Из теоремы 1.10 следует, что для любого фиксированного q 2Xp

Yk(q;x;r) ! 0 при r ! ¥ равномерно по x 2[0;T ] для любого отрезка

[0;T ];0 < T < ¥. Таким

p

Yk(q; ; ) 2C([0;T ];C0

(L))

для

 

 

 

образом, имеем ˆ

 

 

любого фиксированного q 2 G0 (L).

p

, и пусть последователь-

 

3.

Зафиксируем произвольный

 

 

p

 

q0 2 G0 (L)

 

 

ность

qm 2 G0 (L) такова, что qm ! q0 по норме Xp.

 

 

В силу теоремы 1.10 имеем (в тех же обозначениях)

dk(q0; ; ) 2BC([0;¥);C0(S)):

Следовательно, найдется R > 0, такое что jdk(q0;x;r)j 1=2 при всех x 2 [0;T ], r 2 L : jrj > R. Тогда в силу

sup sup jdk(qm;x;r) dk(q0;x;r)j ! 0;

x2[0;T ] r2S;jrj>R

94

Глава 1. Решения типа Вейля

sup sup jdjk(qm;x;r) djk(q0;x;r)j ! 0

x2[0;T ] r2S;jrj>R

(что следует из теоремы 1.10), получаем

sup sup jbjk(qm;x;r) bjk(q0;x;r)j ! 0

x2[0;T ] r2S;jrj>R

при m ! ¥. С учетом непрерывности и равномерной ограниченности

˜

Y0(x;r) отсюда заключаем:

 

sup

 

 

sup

 

ˆ

;x;r)

ˆ

;x;r)k ! 0; m ! ¥:

 

 

 

 

kYk(qm

Yk(q0

 

x2[0;T ] r2

S

;jrj>R

 

 

 

 

 

 

 

В то же время из пункта 1 следует, что

 

 

 

sup

 

 

sup

 

ˆ

 

ˆ

;x;r)k ! 0; m ! ¥:

 

 

 

 

kYk(qm;x;r)

Yk(q0

 

x2[0;T ] r2L;jrj R

 

 

 

 

 

 

 

Таким образом, имеем

 

 

 

 

 

 

 

 

lim

ˆ

(q

; ; )

ˆ

(q

; ; )

 

= 0:

 

 

m

!

¥ kYk

 

m Yk

0

kC([0;T ];C0(L))

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 1.11. Пусть p > 2. Тогда справедливо представление

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

Yk(q;x;r) = Y0k(x;r) +Wk(rx)Yk(q;x;r);

ˆ

p

(l);C([0;T ];H (l))) для любого луча l = fr = zt;t 2[0;¥)g,

где Yk 2C(G0

z 2 S n f0g и любого отрезка [0;T ];0 < T < ¥.

Доказательство. Подставляя в СЛАУ

Fk 1 ^Yk = Fk; Yk ^Tk = 0

фундаментальные тензоры в виде (здесь и далее в рамках данного доказательства W обозначает W (rx))

F

= F0

!

 

; F = F0

!

 

= T 0

 

 

 

+ W k 1Fˆ

+ W kFˆ ; T

+ W kTˆ

 

k 1

k 1

 

k 1

k

k

k

k

k

 

k

 

 

 

 

ˆ

 

 

 

 

 

 

 

и, выполнив замену Yk = Y0k +WkYk, преобразуем исходную систему

к следующему виду:

 

 

 

 

 

 

 

 

 

 

 

˜0

ˆ

= fk;

ˆ

˜ 0

 

 

 

 

(1.128)

где

 

Fk 1

^Yk

Yk ^Tk = gk;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ ˆ

˜

ˆ

ˆ

 

˜

 

ˆ

ˆ

ˆ

 

fk = Fk Fk 1 ^Y0k Fk 1 ^Yk;

gk = Y0k

^Tk Yk ^Tk:

Из теоремы 1.6 известно, что функции

ˆ

 

ˆ

 

ˆ

; ) при-

Fk 1( ; ; );Fk( ; ; );Tk( ;

надлежат C(Xp;BC([0;¥);H (l))).

95

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

 

ˆ

p

(l))). С уче-

В силу леммы 1.12 имеем Yk( ; ; )2C(G0 (l)˜;C([0;T ];C0

том непрерывности и ограниченности функции Y0k(x;r) заключаем, что

функции fk( ; ; ) и gk( ; ; ) принадлежат C(G0p(l);C([0;T ];H (l)))

 

ˆ

 

 

Далее будем искать Yk(q;x;r) в виде

 

 

ˆ

n

˜

 

ˆ

(1.129)

Yk(q;x;r) = åbjk(q;x;r)Y0 j(x;r):

j=1

Подстановка (1.129) в систему (1.128) приводит к следующим представлениям для коэффициентов разложения:

 

 

 

 

 

 

 

 

ˆ

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bjk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (

)j k

W n

 

1

 

˜

˜

˜

 

˜

^gk

 

;

j

<

k

;

1

 

D0k!

 

Y01

^ ^Y0; j 1 ^Y0; j+1 ^ ^Y0;k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bjk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)j k f

W n

1

 

 

˜

˜

˜

˜

 

 

 

 

 

 

 

 

 

= (

 

 

 

 

;

j

= ;

n

;

 

1

j j

!

 

 

fk ^Y0k ^ ^Y0; j 1

^Y0; j+1

^ ^Y0n

 

 

n k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

!

 

 

r

x)

из которых с учетом доказанного выше и свойств функции

W

(

 

 

 

 

 

 

 

 

ˆ

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

(см. парагр. 1.1) заключаем, что bjk 2C(G0 (l);C([0;T ];H (l))). Теперь,

используя снова свойства функций

˜

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y0 j(x;r), из (1.129) получаем тре-

буемое утверждение.

2

Дальнейшее уточнение полученных результатов также возможно при условии дополнительной гладкости и убывания потенциала при x ! 0.

Обозначим через P(S) множество функций, представимых в виде

F(r) = å f (l)exp(lr);

l2L

где множество L (своё для каждой функции F 2 P(S)) таково, что Re(lr) < 0 для всех l 2 L, r 2 S. Заметим, что множество скалярных функций из P(S) образует алгебру относительно поточечных арифметических операций.

Теорема 1.12. Пусть q( ) удовлетворяет условиям теоремы 1.7. Тогда при каждом фиксированном x > 0 и r ! ¥, r2S справедлива

асимптотика

r(Y(q;x;r) Y0(x;r))exp( rxR) = fГ(x) + qˆ(x)f + E(x;r) + o(1);

96

Глава 1. Решения типа Вейля

где Г(x) – некоторая диагональная матрица-функция, E(x; ) 2P(S).

Доказательство. В силу теоремы 1.7 справедливы асимптотики

 

˜

˜0

 

 

(x)fa + E(x;r) + o(1);

 

rFk(q;x;r) = rFk (x;r) + å fk;a

 

 

 

 

 

a2Ak

 

 

 

 

 

˜

 

˜ 0

˜ 0

 

 

 

 

rTk(q;x;r) = rTk

(x;r) + d0kTk (x;r)+

 

 

 

+

 

å Tk0;a (k)gk;a;a (k)(x)fa + E(x;r) + o(1);

(1.130)

 

 

a2An k+1

 

 

 

 

 

Y˜ 0k(x;r) = fk + E(x;r) + O r 1

;

(1.131)

где

a (k) := (k;:::;n)

,

a (k) := (1;:::;k). Здесь и

 

 

 

 

 

далее одним и тем же

символом E( ; ) обозначаются различные функции, такие что E(x; ) 2

2 P(S).

Дальнейшие рассмотрения проводятся при (произвольном) фиксированном x, который для краткости обозначений не будет указываться

в списках аргументов.

 

 

 

 

 

 

 

 

Перепишем

соотношения,

определяющие решения типа Вейля,

в виде

 

 

 

 

 

 

 

 

 

 

 

 

˜

˜

˜

˜

˜

 

 

 

 

Fk 1

^Yk

= Fk;

Yk ^Tk = 0:

 

 

Подставив в эти соотношения

˜

 

 

 

 

 

 

Yk в виде

 

 

 

 

 

 

 

 

˜

˜

ˆ

;

 

 

 

(1.132)

 

 

 

Yk

= Y0k

+ Yk

 

 

 

получим

 

 

 

 

 

 

 

 

 

 

˜

ˆ

˜

˜

˜

ˆ

 

˜

˜

˜

Fk 1

^Yk = Fk Fk 1 ^Y0k; Yk

^Tk =

Y0k ^Tk:

Полученные соотношения преобразуем в СЛАУ

 

 

 

 

n

 

 

 

 

 

 

 

 

 

åmi jgjk = ui; i =

 

 

 

(1.133)

 

 

1;n

 

 

 

j=1

 

 

 

 

 

 

 

относительно коэффициентов fgjkg разложения

 

 

 

 

 

ˆ

n

 

 

 

 

 

(1.134)

 

 

 

 

 

 

 

 

 

 

 

 

Yk = ågjkf j:

 

 

 

 

j=1

Коэффициенты в системе (1.133) определяются следующими соотношениями:

 

 

 

mi j

=

F˜k 1 ^f j ^fa

 

;

 

 

ui =

 

˜ ˜

˜

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Fk Fk 1

^Y0k) ^fa

; a = a (k) ni; i = k;n;

 

 

 

 

 

 

 

 

 

 

 

 

97

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

 

 

 

˜

 

 

 

 

 

˜

˜

 

 

 

 

 

 

mi j =

 

 

; ui =

 

 

; a

= a (k 1) ni; i = 1;k 1:

fa ^f j ^Tk

fa ^Y0k ^Tk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пользуясь

(1.130),

(1.131)

и учитывая, что

 

 

 

 

 

 

 

 

 

 

˜0

 

˜

˜0

 

˜

 

˜ 0

= 0;

 

 

 

 

 

 

 

 

Fk 1 ^Y0k = Fk

; Y0k ^Tk

 

 

получим асимптотики для коэффициентов СЛАУ (1.133) при r ! ¥

( ) =

O r

1

 

; j = i; m ( ) = m0 + O

r

1

;

 

0

= (

 

 

)k i

 

f ;

 

=

k;n

;

mi j r

 

 

 

6

ii r

 

 

ii

 

 

mii

 

 

 

 

 

1

 

 

 

 

j j

 

i

(1.135)

и

 

mi j(r) = mi0j + O r 1 ; i =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

; j =

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

1;k 1

1;k 1

 

 

 

 

 

 

 

 

 

 

 

mi j(r) = mi0j + E(r) + O r 1 ; i =

 

;

j =

 

;

 

 

 

 

 

 

 

где

1;k 1

k;n

 

 

 

 

 

 

 

 

 

 

 

mi0j

= T 0

fa

^

f j

^

fa (k)

; a = a

(k

 

1)

n

i;

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

k;a (k)j

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

; j = i; j < k;

 

 

 

 

 

 

1

 

 

 

=

 

;

( ) =

 

( ) =

( ) +

 

 

 

 

;

 

 

k;n

mi j r

 

O

 

 

 

6

 

 

 

mi j r

 

E r

 

O

 

 

 

 

 

 

j

 

 

(1.136)

для i =

 

mii(r) = mii0 + O r 1 ; mii0 = ( 1)k 1 ijfjTk0;a (k)

 

 

 

 

 

 

 

 

 

 

(1.137)

1;k 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично вычисляем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rui(r) = ui1 + E(r) + o(1);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.138)

u1i = ( 1)k ijfjfka dikjfjf(k 1);a (k 1); a = a (k 1) [fig; i = k;n;

(1.139)

где dik – символ Кронекера,

u1i = ( 1)k i jfjTk0;a (k)gk;b;a (k); b = a0 nk;a = a (k 1) ni; i = 1;k 1:

(1.140) Из СЛАУ (1.133) следует предварительная оценка gik(r) = O r 1 при r ! ¥. Подставляя в (1.133) gik(r) виде gik(r) = r 1ik(r), где ik(r) =

= O(1), получим при i = k;n

mii(r)gˆik(r) = u1i + E(r) åmi j(r)gˆ jk(r) + o(1);

j6=i

что с учетом (1.135), (1.138) приводит к асимптотическим формулам

u1

ik(r) = g1ik + E(r) + o(1); g1ik = i0 ; (1.141)

mii

98

Глава 1. Решения типа Вейля

i = k;n.

Аналогично для i < k имеем

ˆ

1

ˆ

å

ˆ

mii(r)gik(r) = ui

+ E(r) åmi j(r)gjk(r)

mi j(r)gjk(r) + o(1):

 

 

j k

j<k; j6=i

 

В силу (1.136) данное соотношение преобразуется к виду

m0iiik(r) = u1i + E(r) åmi j(r)gˆ jk(r) + o(1):

j k

Пользуясь снова (1.136) и полученными выше асимптотиками (1.141) для i = k;n, убеждаемся, что формулы (1.141) верны и для i < k.

Далее, для i = k;n из (1.141), (1.139), (1.135) получаем

1

˜1

+ fka; a = a (k 1) [i:

(1.142)

gik

= di;kgik

 

Из теоремы 1.7 следует, что

 

(k)(x)fa (k)

^fa0

 

 

 

fka(x) = ca

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Напомним, что (произвольный) линейный оператор V , действующий в Cn, может быть продолжен на внешнюю алгебру ^Cn таким образом, что равенство

V (h1 ^ ^hm) = (V h1) ^ ^(V hm)

справедливо для любого набора векторов h1;:::;hm, m n; при этом для любого h2^nCn имеем V h = jV jh (где jV j обозначает определитель матрицы оператора V в стандартном базисе fe1;:::;eng). В последующих выкладках f обозначает соответствующее расширение оператора, порожденного матрицей перестановок f. Кроме того, будем учитывать, что

(f 1V f)(k) = f 1V (k)f

для любой n n матрицы V . С учетом сказанного

имеем

 

 

 

 

 

 

 

 

 

 

 

 

fka(x) = ca f f 1(k)(x)fea (k) ^(fea0 )

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

(k)

 

 

 

 

1

 

 

 

(k)

ea (k) ^ea0 :

= jfa ^fa0 jjfj f qˆ

 

(x)fea (k) ^ea0 = jea ^ea0 j f qˆ(x)f

 

 

 

 

 

 

 

 

1)

[

i

 

 

 

Заметим далее,

что для мультииндекса

a = a

(k

 

 

и произвольной

n n матрицы V справедливо равенство (проверяемое непосредственным вычислением)

jea

^ea0 j V (k)ea (k)

^ea0 = Vik:

 

 

 

 

 

 

 

Подставляя полученные соотношения

в (1.142),

получаем

 

gik1

= di;kik1 + f 1qˆf ik ; i =

 

:

(1.143)

k;n

99

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

Аналогично из (1.136), (1.140) получаем для i < k

 

 

 

 

 

gik1 = gk;b;a (k);

b = a0 nk;a = a (k 1) ni:

(1.144)

Из теоремы 1.7 следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gkab(x) = ca

 

(n k+1)(x)fb ^fa0

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

получим

 

 

 

 

для b = a. Повторяя рассуждения,

приведенные выше,

 

 

 

 

gkab(x) = jea ^ea0 j

f 1qˆ(x)f (n k+1) eb ^ea0 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В частности, имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gk;b;a (k) = jeb ^eb0 j f 1qˆ(x)f (n k+1) ea (k) ^eb0 :

 

 

 

 

Наконец, для b = a0

n

k;a =

 

(k

 

1)

n

i, i < k

 

 

 

n

n

 

a

 

 

 

и произвольной

 

матрицы V имеем

 

 

 

 

 

 

 

 

 

 

^eb0 = Vik:

 

 

 

 

 

jeb ^eb0 j V (n k+1)ea (k)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя полученные соотношения

в (1.144),

получаем

 

 

 

 

 

 

gik1 = f 1qˆf ik ; i =

 

:

 

 

(1.145)

 

 

1;k 1

 

 

Подставляя (1.145) и (1.143) в (1.141), приходим к равенствам

 

 

 

 

rgik = di;k1ik + f 1qˆf ik + E(r) + o(1);

что в терминах матрицы g = (gik)i;k=1;n запишется в виде

rg(x;r) = Г(x) + f 1qˆ(x)f + E(x;r) + o(1);

где Г(x) – диагональная матрица. Переписав (1.134) в матричном виде

ˆ

как Y(x;r) = fg(x;r), получаем окончательно требуемый результат.

2

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