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

книги / Теория функций комплексной переменной

..pdf
Скачиваний:
3
Добавлен:
19.11.2023
Размер:
28.91 Mб
Скачать

МЕТОД ВИНЕРА-ХОПФА

291

где функции L+(k) и L-(k) являются

аналитическими при

Im к > р и Im к < т+ соответственно. Кроме того, предполо­ жим, что в областях своей аналитичности эти функции на бес­ конечности растут не быстрее, чем кп, где п — некоторое поло­ жительное целое число. Разбиение (51) аналитической функции L(k) обычно называется факторизацией. Возможность факто­ ризации заданной аналитической функции комплексной пере­ менной будет обоснована ниже (см. леммы 1 и 2 на с. 294, 296).

Итак, в результате факторизации исходное уравнение приве­

дено к виду

 

L+(k)U+(k) = -L-(k)U-(k).

(49)

Из предыдущих рассмотрений следует, что оно определяет некоторую целую функцию комплексной переменной к.

Так как U±(k) -> 0 при |&| —> оо, a L±(k) растут на беско­ нечности, как конечная степень к71, то данная целая функция может быть лишь полиномом Pn-i(k) степени не выше п 1 .

Если функции L±(k) растут на бесконечности, лишь как пер­ вая степень переменной к, то из соотношений (50) в силу теоре­ мы Лиувилля следует, чтр соответствующая целая функция есть постоянная С. Тогда для неизвестных U+(k) и U-(k) получим выражения

и+ м = 1 ж г = (52)

определяющие преобразования Фурье искомого решения с точ­ ностью до постоянного множителя, который может быть найден хотя бы из условий нормировки. В общем случае выражения

и+ м = Ш у и- {к)=

(53)

определяют преобразования Фурье искомого решения инте­ грального уравнения (38) с точностью до неопределенных по­ стоянных, которые можно найти из дополнительных условий задачи. Само решение определяется с помощью обратного пре­ образования Фурье (31) и (35).

Рассмотрим пример применения изложенного метода. П р и м е р 2. Рассмотрим уравнение

ОО

 

и(х) = A J e“ lx -slu(s) ds,

(54)

о

ядро которого имеет вид i?(f) = е 1^1.

Найдем преобразование Фурье функции v(£):

ОО

(55)

т - Т В £ « * * * « - Т щ к г у

292

ПРИЛОЖЕНИЕ 2

Функция V (к) (55) является аналитической функцией комплекс­ ной переменной к в полосе —1 < Im к < 1. Представим выраже­ ние

Цк) = 1 - V S XV(к) = ** ~ г(2+А ~

(56)

в виде (51), где

£ + (к) =

£ _ (*) = * - < •

(57)

Функция £+(&) в (57) является аналитической и отличной от нуля функцией /с в области Im к > Im у/2 Х — 1. При 0 < Л < < - эта область определяется условием Im к > л/1 2А, причем

у/ 1 2А ^ jj, < 1 . При Л > - функция L+(k) является аналити-

а

ческой и отличной от нуля в области Im к > 0. Функция £_(& ), очевидно, представляет собой отличную от нуля аналитическую

функцию в области Im к < 1 . Поэтому при 0 < Л < -> обе функции удовлетворяют требуемым условиям в полосе ц < Im к < 1 .

При iZ < Л общей областью аналитичности функций £+(&) и

L-(k) является полоса 0 < Im А: < 1 . Таким образом, необходи­ мая факторизация функции (56) произведена.

Рассмотрим выражения U±(k)L±(k). Так как U±(k) —> 0 при |&| —> оо, a L±(A;), согласно (57), растут на бесконечности, как первая степень к, то целая функция Рп(к), совпадающая с

U+(k)L+(k) при Im к > у. и с U-(k)L-(k) при Im к < 1 , может быть лишь полиномом нулевой степени. Поэтому

 

Щ(к)Ь+(к) =

С.

 

 

(58)

Отсюда

 

 

 

 

 

 

 

 

U+W =

C » - Z

+ l

 

 

(59)

 

 

 

 

и, согласно (31),

 

 

 

 

 

 

 

 

 

оо+гт

 

 

 

 

 

С

А

k + i

-ikx

 

 

и+(х) =

I

JU

(60)

,/2 ?

1

к2 (2A — 1)6

 

aK’

 

 

 

—оо+гт

где fJL < Т < 1.

Для вычисления интеграла (60) можно применить методы гл. 5. Замкнув контур интегрирования при х > 0 дугой полу­ окружности в нижней полуплоскости и оценив интеграл по этой дуге с помощью леммы Жордана, после элементарных вычисле-

МЕТОД ВИНЕРА-ХОПФА

293

ний получим

и+(х) = D cos

(61)

 

л/2А — 1

где D — новая постоянная. При 0 < А <

^ это решение экспо-

ненциально возрастает с ростом ж, при |

< А < оо — ограничено

на бесконечности.

Итак, уже пример решения однородного интегрального урав­ нения (38) выявляет основную идею метода Винера-Хопфа, за­ ключающуюся в представлении с помощью факторизации ис­ ходного функционального уравнения (47) в виде целой функ­ ции (49). Дадим теперь обоснование возможности факторизации аналитической функции комплексной переменной, причем бу­ дем исходить из несколько более общего функционального урав­ нения, чем уравнение (47).

4. О бщ ая схем а метода В и н ер а-Х оп ф а. В общем случае задача, решаемая методом Винера-Хопфа, сводится к следую­ щей.

Требуется определить функции Ф+(&) и Ф_(&) комплексной переменной к, аналитические соответственно в полуплоскостях Im к > т_ и Im к < т+ (т_ < т+), стремящиеся к нулю при |&| —> оо в своих областях аналитичности и удовлетворяющие в

полосе (т_ < Im к < т+) функциональному уравнению

 

A(k)V+(k) + В (*)Ф _ (*) + С(к) = 0.

(62)

Здесь А{к), В (к), С (к) — заданные функции комплексной переменной к, аналитические в полосе т_ < Im к < т+, причем А(к) и В (к) отличны от нуля в этой полосе.

Основная идея решения этой задачи основана на возмож­ ности факторизации выражения А(к)/В(к), т. е. возможности представить его в виде

А(к) _ Ь+(к) В (к) L - (к)

где функции L+ (к) и L - (к) являются аналитическими и отлич­ ными от нуля соответственно в полуплоскостях Im к > т_ и Im к < т+, причем полосы т_ < Im к < т+ и т'_ < Im к < т’+

имеют общую часть. Тогда с помощью (63) уравнение (62) мож­ но переписать в виде

L + № + ( k ) + £_(*)*_(*) + £ _(*)§jg = 0.

(64)

Если последнее слагаемое в (64) можно представить в виде

L- Ю Щ = D+W + D-(k),

(65)

294 ПРИЛОЖЕНИЕ 2

где функции D+(k) и D-(k) являются аналитическими в полу­

плоскостях Im к > T!L и Im к < т"

соответственно, и все три

полосы т_ < Im к < т+, т!_ < Im к < т+ и т" < Im к <

имеют

общую часть — полосу

< Im к <

, то в этой полосе имеет

место функциональное уравнение

 

 

L+ (k)V+ {k) +

D+{k) = -Ь-{к)Ъ -(к) - D-(k).

(66)

Левая часть (66) представляет собой функцию, аналитическую в полуплоскости т£ < Im к, правая — функцию, аналитиче­ скую в области Im к < т\. Из равенства этих функций в полосе т£ < Im к < т+ следует, что существует единственная целая функция Р{к), совпадающая соответственно с левой и правой частями (66) в областях их аналитичности. Если все функции, входящие в правые части (63) и (65), растут на бесконечности

в своих областях аналитичности не быстрее, чем кп+1, то из условия Ф±(&) —)• 0 при \к\ —> оо следует, что Р{к) является

полиномом Рп(к) степени не выше п. Тем самым равенства

« + (* ) =

ь ф )

(67)

И

ф -(*) = P- (L _ ( ^ ~ W

(68)

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

Итак, применение метода Винера-Хопфа основано на пред­ ставлениях (63) и (65). Возможность этих представлений обес­ печивается следующими леммами.

Лемма 1 . Пусть функция F(k) является аналитической

в полосе т_ < Im к < т+, причем в этой полосе F(k)

равно­

мерно стремится к нулю при \к\

оо. Тогда в данной полосе

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

 

 

F(k) = F+ (k) +

F-(k),

(69)

где функция F+(k) аналитическая в полуплоскости Im к >

>т_, а функция F-(k) — в полуплоскости Im k < r + .

До к а з а т е л ь с т в о . Рассмотрим произвольную

точку ко, лежащую в данной полосе, и. построим прямоуголь­ ник abed, содержащий точку ко, внутри и ограниченный отрез­ ками прямых Im к — т'_, Im к = т'+, Re к = —A, Re к = А,

МЕТОД ВИНЕРА-ХОПФА

295

где т<т'_ <т'+ <т+ (рис. 2 ). По формуле Коши

A+ir'_ А+гт^.

**>-£ /

С -к о

+ A / Г

^

+

 

—A+ir'

 

A+ir'

 

ко

 

 

 

 

—A+ir!_

 

 

 

—A+ir+

 

 

 

+

J _

[

F(0

dC +

J L

/

n o

<K- (70)

 

2m

J

£ —ko

^

2m

J

С —ко

 

 

 

A+ir+'

 

 

 

—A + I T + 1

 

Перейдем в (70) к пределу при А

оо. Так как по условию лем­

мы F(k) равномерно стремится к нулю при \к\ —> оо, то предел второго и четвертого слагаемых в правой части (70) равен нулю,

и мы получим

 

 

F(k0) = F+ (k0) + F-(k0),

(71)

где

 

 

oo+ir'

 

 

)

<F- l d<>

(72)

 

oo+iT '_

оо-H T +'

(73)

oo+ir^.

Интегралы (72) и (73), как интегралы, зависящие от парамет­

р а 1),

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

функции комплексной пере­

менной ко при условии,

 

что точка ко не лежит

 

на контуре

интегриро­

 

вания.

 

F+ (ko)

 

В

частности

 

является

аналитиче­

 

ской функцией в полу­

 

плоскости Im ко > r l , а

 

.Р_(А;о) — в полуплоско­

 

сти Im ко <

7+. В силу

 

произвольности выбора

 

точки

ко

и

прямых

 

т!_ и

соотношения

Рис. 2

(71)-(73) и доказывают

лемму.)*

 

 

 

*) См. с. 69.

296

ПРИЛОЖЕНИЕ 2

За м е ч а н и е 1 . Заметим, что из сходимости интегралов

(72)и (73) следует ограниченность построенных таким образом функций F+ (к) и F-(k) при \к\ -> оо в данной полосе.

Лемма 2 . Пусть функция Ф (к) является аналитической и отличной от нуля в полосе т_ < Im к < т+, причем Ф (к) равномерно в этой полосе стремится к единице при |/?| —> оо. Тогда в данной полосе имеет место представление

Ф(/с) = Ф+(&) • Ф_(&),

(74)

где функции Ф+(к) и Ф-(к) являются аналитическими и от­ личными от нуля соответственно в полуплоскостях Im к > г- и Im к < т+ .

Д о к а з а т е л ь с т в о . Рассмотрим функцию F(k) = = 1п Ф (к), которая, очевидно, удовлетворяет всем условиям лем­ мы 1 . Поэтому для функции F(k) возможно представление (71)—(73). Полагая

Ф+(к) = exp [F+(k)],

Ф-(к) = exp [ Л (*)],

(75)

где функции F+ (к) и F-(k) определены формулами (72),

(73),

получаем

 

 

1пФ+(1Ь) = F+(k),

1пФ -(к) = F_(fc).

(76)

Тогда формула (71) дает

 

 

1пФ(/г) = \ п Ф + ( к ) + 1пФ_(/г),

(77)

откуда и следует соотношение (74). Так как функции * + (* ) и F-(k), согласно лемме 1, являются аналитическими в полуплос­ костях Im к > т- и Im к < т+ соответственно, то и функции Ф+(&) и Ф-(&), определенные по формулам (75), будут обладать требуемыми свойствами. Лемма доказана.

З а м е ч а н и е 2 . Возможность факторизации (74) сохра­ няется в том случае, когда функция Ф(&) имеет конечное число нулей к{ в полосе т_ < Im к<т+.

Для доказательства леммы 2 в этом случае достаточно ввести

вспомогательную функцию

 

 

F(k) = In

2 + Ь2)N/2 Ф(к)

(78)

 

П(к ~~ ki)a*

 

где щ — кратность нуля k{\N полное число нулей с учетом их кратности; положительная постоянная Ь > |т_|, |т+ | выбирается из условия, чтобы функция, стоящая под знаком логарифма, не имела дополнительных нулей в полосе т_ < Im к < т+ . Послед­ няя функция, очевидно, стремится на бесконечности к единице. Построенная таким образом функция F(k) по-прежнему удовле­ творяет всем условиям леммы 1 .

Доказанные леммы и определяют возможность представле­ ний (63), (65), составляющих основу метода Винера-Хопфа.

МЕТОД ВИНЕРА-ХОПФА

297

Мы рассмотрели применение метода Винера-Хопфа для ре­ шения функционального уравнения (62). Легко видеть, что к

этому уравнению сводится и

н е о д н о р о д н о е

и н т е ­

г р а л ь н о е у р а в н е н и е

на полубесконечном промежутке

с ядром, зависящим от разности аргументов:

 

оо

 

 

гх(я) = Л J v(x —s)u(s) ds + f(x).

(79)

о

 

 

Будем предполагать, что ядро уравнения (79) и функция f(x) удовлетворяют условиям (43), и будем искать решение

уравнения (79), удовлетворяющее условию1)

 

|гх+(ж)| < Мхё1* при х ->• оо (ii<T+).

(80)

Тогда, проводя рассуждения, аналогичные рассуждениям при выводе функционального уравнения (47) для однородно­ го интегрального уравнения, получаем, что в случае уравнения (79) в полосе /х < 1ш к < т+ должно удовлетворяться функцио­ нальное уравнение

U+{k) + U~{k) = A^ V (k )U + (k ) + F+ (k) + F-{k)

(81)

или

 

Цк)17+(к) + U-(k) - F(k) = 0,

(82)

где

 

L(k) = 1 - VbrXV(k).

(83)

Уравнение (82) является частным случаем уравнения (62). Функция L(k) в полосе т_ < Im к < т+ является аналитиче­ ской и равномерно стремится к единице при \к\ -» оо, так как \V(k)\ —> 0 при |&| —¥ оо. Если, кроме того, функция L(k) име­ ет конечное число нулей в этой полосе, то все условия леммы 2 выполнены и функцию L(k) можно представить в виде

m =

(84)

где L+(k) является аналитической функцией в верхней полу­ плоскости Im к > т_, a L-(k) в нижней полуплоскости Im к < < т+ . Тогда уравнение (82) принимает вид

L+(k)U +(k)+L-(k)U -(k)-L-(k)F+{k)-F+ (k)L.(k) = 0. (85)

Для приведения последнего уравнения к виду (66) достаточно разложить последнее слагаемое:

F + (k)L. (к) = D+ (к) + D_ ),

(86)

х) Мы опять не останавливаемся на обосновании существования решения уравнения (79), удовлетворяющего условию (80).

298 ПРИЛОЖЕНИЕ 2

на сумму функций D+(k) и D -(k), являющихся аналитическими в полуплоскостях Im к > р и 1ш к <т+ соответственно.

Для обоснования возможности такого представления заме­ тим, что в силу условия (43) функция (к) является аналити­ ческой в верхней полуплоскости Im к > т_ и равномерно стре­ мится к нулю при \к\ —> оо. Функция L-(k) является анали­ тической в нижней полуплоскости Im к < т+, и по способу ее построения в силу леммы 2 и замечания к лемме 1 можно так провести факторизацию (84), чтобы L -(k ) оставалась ограни­ ченной в полосе г < Im к < т+ при |/г| -> оо. Отсюда следует, что для функции (к)Ь-(к) в полосе т_ < Im к < т+ выпол­ нены все условия леммы 1 , что и достаточно для обоснования представления (86).

Проведенные рассмотрения позволяют при дополнительных условиях, что функции L±(k) растут на бесконечности не бы­ стрее, чем кп, представить преобразования Фурье решения неод­ нородного интегрального уравнения (79) в виде

U+ (k) = **«(*)+ д +(*)>

= -Pn(k) + L-(k)F-(k) + D-(k)

 

 

L-(k)

 

, ч

(87)

Само решение может быть получено из (87) с помощью формул

(31)и (35) обратного преобразования Фурье.

5.Задачи, приводящие к интегральным уравнениям

сядром, зависящим от разности аргументов.

5.1.Вывод уравнения Милна. К интегральным уравнениям

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

Пусть в полупространстве х > 0, заполненном однородным веществом, плотность которого определяется числом щ частиц в единице объема, распространяется поток нейтронов. Будем счи­ тать, что частицы вещества являются тяжелыми атомами, рас­ сеивающими нейтроны так, что абсолютная величина скорости нейтронов остается постоянной, а меняется лишь ее направле­ ние. Рассмотрим стационарный процесс и предположим, что все нейтроны имеют одну и ту же абсолютную величину скорости «о = 1 и плотность их распределения зависит лишь от одной ко­ ординаты х. Введем функцию f(x,p), характеризующую плот­ ность нейтронов в сечении х, скорость которых составляет с по­

ложительным направлением оси х угол в, где р — cos в *). Чис­ ло нейтронов в единице объема в данном сечении, направление

*) Очевидно, —1 ^ р ^ 1 при 0 ^ в ^ тг.

МЕТОД ВИНЕРА-ХОПФА

299

скорости которых лежит в пределах (р,р -J- dp), определяется величиной /(ж, р) dp.

Полная плотность нейтронов р(х) в данном сечении равна

1

р(х) = f f(x,p) dp.

(88)

- 1

 

Нашей ближайшей целью является вывод уравнения для функ­ ции распределения /(ж, у). Для этого составим соотношение пол­ ного баланса числа нейтронов, имеющих направление скорости в интервале (р, р + dp) и находящихся в слое между сечениями ж и ж + dx. В силу стационарности процесса поток нейтронов, выходящих из данного слоя

pf(x + dx,p)dp-pf(x,p)dp,

(89)

определяется разностью между числом нейтронов, приобретших скорость в заданном направлении (р, p+dp) в результате рассея­ ния на частицах вещества в данном слое, и числом нейтронов, имевших скорость в заданном направлении и изменивших это направление после рассеяния. Мы будем считать, что рассеяние нейтронов на частицах вещества является изотропным, т. е. оно равновероятно во всех направлениях, и вероятность рассеяния нейтрона на одной частице характеризуется эффективным се­ чением рассеяния Q. Тогда число нейтронов, имевших заданное направление скорости (р, р + dp) и рассеянных в данном слое, очевидно, равно

/ (ж, д) dp ’ Qno dx,

(90)

а число нейтронов, приобретших в результате рассеяния ско­ рость в требуемом направлении, равно

1

| dpQ • щ dx J /(ж, р') dp1.

(91)

- 1

 

На основании выражений (89) - (91) уравнение баланса запи­ шется в виде

p f (ж + dx, р) dp - pf(ж, р) dp =

1

= —Q - nof(x>p)dpdx + ^Y-dpdx J f(x,p')dp'. (92) - l

Разделим обе части равенства на dpdx и перейдем к пределу при dx —> 0. Учтя (88), получим уравнение для функции рас­ пределения нейтронов в виде

/*§ £ = ~Q nof(x,n) + ^р р (х).

(93)

300

ПРИЛОЖЕНИЕ 2

Это уравнение часто называется уравнением переноса или транспортным уравнением. Оно справедливо не только в случае рассмотренной конкретной физической задачи, но и для многих других физических процессов, связанных с переносом вещества

или излучения1).

Для дальнейшего удобно переписать уравнение (93) в несколько ином виде, введя безразмерную пространственную ко­

ординату £, связанную с х соотношением х = А£, где Л = 1 nQo

средняя длина свободного пробега. Тогда уравнение переноса примет вид 2)

р Ц = -/(£ > Р) + \ р(0 -

(94)

Функция /(£, д) должна быть подчинена граничным услови­ ям, вытекающим из физической постановки задачи. Будем счи­ тать, что поток нейтронов из внешнего полупространства £ < 0 равен нулю, а при £ -> оо имеется постоянный поток нейтро­ нов единичной мощности в отрицательном направлении оси £ (т.е. при £ -> оо отсутствуют нейтроны, направление скоро­ сти которых составляет с отрицательным направлением оси £ острый угол, отличный от нуля). Тогда граничные условия для функции /(£,/2) запишутся в виде

/ ( 0, д ) = 0,

 

0,

/ ( oo,/i) = 0,

1

(95)

< у < 0.

Установим важные следствия уравнения (94) и условий (9 5), для чего сначала проинтегрируем (94) по ц:

1

1

1

щ /

f(£,p)pd(i = - f

+ ip ( 0 f dn =

= - Ж ) + / » (? )= 0. (96)

Так как интеграл i(£) = / /(£,//)// d/i равен потоку нейтронов

- 1

через данное сечение, то уравнение (96) дает

 

щ = 0 или Д £) = const.

(97)

*) Подробный вывод уравнения переноса для более общих случаев см., например, в кн. : М о р с и Ф е ш б а х . Методы теоретической физики.

Т.1. - М. : ИЛ, 1958.

2)Мы сохранили для функции /(£, у) старое обозначение.

Соседние файлы в папке книги