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

книги из ГПНТБ / Земанян А.Г. Интегральные преобразования обобщенных функций

.pdf
Скачиваний:
37
Добавлен:
25.10.2023
Размер:
13.27 Mб
Скачать

І квадратично интегрируема на 0 < І°о , то поставлен­

ная задача решается с

помощью подходящего ортоиор-

мального разложения;

см.,

например,

Армстронг [1],

Хаггинс [1J, Каутц [1],

Ли

[1], Лернер

[1], Янг и Х аг­

гинс [1].

эту

классическую задачу, допу­

Мы обобщим теперь

стив, что сигнал / может быть обобщенной функцией, и бу­ дем искать в качестве решения такую последовательность

{/N

 

= . каждый элемент которой является обычной дей­

(0}

n

 

o

 

функцией,

обладающей

указанными

ствительнозначной

выше свойствами 1) и 2), причем /,ѵ ->- / на 0 <

t

< оо в не­

котором обобщенном смысле.

 

 

 

 

Если

/ ЕЕ

А ' ,

причем

Л '

соответствует

 

системе 2

(п. 9.8),

то решение можно найти с помощью преобразова­

ния Лагерра. Выберем порядок

а

для функций фп, соот­

 

ветствующих преобразованию Лагерра обобщенных функ­ ций, равным некоторому неотрицательному четному

числу. Обрежем затем ортонормальное разложение обоб­

щенной функции / на номере

п =

N

и умножим резуль­

t

 

тат на функцию Хевисайда 1+ ( ). Тогда

 

(3)

 

 

 

N

 

 

 

 

 

 

 

/.V (о = и (о S (/> Фп) Фп («)•

 

 

 

 

71— 0

 

 

 

 

(s) функ­

Легко показать, что преобразование Лапласа

 

ции

(t)

является рациональной

функцией,

имеющей

/іѵ s

нуль

при

s =

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

N

+ 1 + а/2

 

 

в точке = — 1/2. Таким образом, свойства 1) и 2) выпол­

нены.

 

 

t

 

А '

Кроме того, мы знаем, что /іѵ ( ) сходится в

 

к /

(t)

при

N

-> оо. Поэтому последовательность элементов

вида (3) является решением.

преобра­

 

Другое решение можно получить, применив

зование Якоби и сделав замену переменных. Предположим

теперь, что

/ €Е

(/),

 

где / — интервалt (0, оо). Пусть

X

(£) — действительная непрерывная неубывающая функ­

ция, тождественно равная —1 на —оо <

<

оо и такая,

что

dx/dt

0 на 0

< Ц

<

о о , а 1ітл;(і) =

1.

Далее, пусть

 

 

/ — ► С О

преобразование Лапласа функции х (t) + 1 является ра­ циональной функцией, имеющей нуль на бесконечности, все полюсы которой лежат в левой половине 5-плоскости (Re s < 0 ) , исключая простой полюс в начале координат. Существует множество способов, в помощью которых мож­ но построить такую функцию х (t) -f 1, подходящим обра­

340

зом выбирая точки расположения полюсов и нулей ее пре­ образования Лапласа (Земанян 113], [14]). Пусть t (ж) обозначает функцию, обратную функции ж (t). Делая за­ мену переменных в формуле для распределений на интер­

вале

О

< it

с ,о о

(Земанян [1],

стр.

30),

мы можем пост­

роить распределение / U (ж)]. На оси ж носитель функции

/ U (ж)]

является компактным

подмножеством

интервала

J

=

(— 1,

1); поэтому /

[t

(ж)] е=

Ш'

(/) d

Л '

, где

Л '

обоз­

начает

теперь пространство

обобщенных

функций,

соот­

ветствующее преобразованию tЯкоби в случае,

когда

а

и

ß — неотрицательные четные целые числа. При таком вы­

боре

а

и

ß

функции

 

фп ( ) — полиномы.

Разложим

/ U (ж)]

в

ряд

/ [*(*)!=

2

ЯпФпОг),

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

п=0

 

 

 

 

 

 

 

 

 

 

а п = (/

[* (*)]»

 

(*))

=

( / (0 .

’IV. [ *

(01 ^ j f )

 

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f ( t ) =

7

21 = 0

Mv[z(0b

t

 

t

 

и

для любого

 

0 функция

Так как фл. — полиномы

 

ж ( ) является конечной линейной комбинацией действи­

тельных постоянных чиселt,

и действительнозначных затуха­

ющих

экспоненциальных

функций,

умноженных на

не­

отрицательные степени

 

 

N

 

 

 

то функция

 

(4)

 

 

 

 

J n (0 = 1+ (0

7

21 = 0

(<)]

 

удовлетворяет

Л '

 

1)(t) и 2). (t)В силуЮ'

свойства(I) NI -п*■.

9.4

свойствам

 

I

=

 

 

влечет сходимость в 25' (/). Отсюда не­

сходимость в

 

медленно следует, что /л

 

 

-ѵ /

в

при

oo,

где

 

 

(0, oo). Таким

 

образом,

последовательность

{/jv}, где /іѵ задано выражением (4), также дает решение нашей задачи.

Отметим, что если мы возьмем

ж (t) — 1 — 2е~р1, 0 < £ < ° ° і

где р — действительное положительное число, то преобра­

зование Лапласа

s

(()

будет иметь только

( ) функции

 

341

простые действительные полюсы, расположенные в точ­ ках вида — 7??.р, где т — неотрицательное целое число. Справедливо п более общее утверждение, а именно; если в X (t) входят только действительные экспоненты, то F n (s) может иметь только действительные полюсы. Если же мы позволим X (t) содержать затухающие синусоидальные экспоненты, то окажется, что Fn (s) может иметь как дей­ ствительные, так и комплексные полюсы.

З а д а ч а 9.10.1. Показать, что двустороннее пребрааованпе Лапласа выражения (3) есть рациональная функция, имеющая нуль при s = оо и только один полюс кратности /V + 1 + а/2, расположенный в точке s — 1/2. Как и раньше, порядок а обоб­

щенных функции Лагерра фп считается равным неотрицательному четному целому числу.

9.11. Применение преобразования Лежандра: внутренняя задача Дирихле для единичной сферы

J Іреобразоваппе Лежандра может быть использовано для решения следующей внутренней задачи Дирихле для единичной сферы. Выберем сферическую систему коор­ динат

{(г, Ѳ, а) : О ^ г < 1 , О ^ 0 ^ л , О ^ о с < 2л}

и предположим, что искомая функция ѵ не зависит от коор­ динаты а (т. е. V — V (г, 0)). Полагая р — cos Ö, мы можем переписать уравнение Лапласа в виде

rD\ rv - f Dp

(1 — р2) Dy.v =

0,

(1)

V = V ( г , arccos р),

0 < ; г < 1, — 1

<

р < 1 .

Далее, потребуем, чтобы функция ѵ (?■ , 0) оставалась огра­ ниченной в некоторой окрестности начала координат г = 0. Пусть, наконец, Л ' — пространство обобщенных функций, соответствующее системе 4а (см. п. 9.8) для пре­ образования Лежандра, и / (р) £Е Л '. Наложим гранич­ ное условие

V (г, arccos р)

/ (р) в Л ' при 7-

1 — 0 .

(2)

Чтобы формально получить решение, мы сначала приме­ ним преобразование Лежандра 11 по р к дифференциаль­

ному уравнению (1). Переставляя операции U и г D\ и полагая

V — V (г, //) = U V (г, arccos р), п = 0, 1, 2, . . .,

342

приходим

к уравнению

п (п

 

 

 

 

 

V =

 

 

 

 

 

(3)

 

 

 

 

 

rD lrV -

 

 

 

 

 

+

 

 

1 )

 

 

 

0 .

 

 

 

Решение уравнения (3) имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

V (г, п)

 

=

А (п) гп

 

+

В

(и) г-""1. г —

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как функция конечна в окрестности точки

 

0, то мы

должны положить

В (п)

=

 

0,

 

чтобы исключить слагаемое,

расходящееся при

г —>- 0.

Фупкции

А

(

п

) определяются

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

из преобразования граничного условия (2); при этом мы

предполагаем,

что

оператор

 

U

и

предельный переход

при

 

->■

1

А п)

 

 

 

 

 

V

 

 

 

п)

 

 

 

 

(х).

 

 

 

 

— 0 перестановочны. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

(

=

 

 

1—0

 

(г,

 

 

=

(U/)

 

 

 

 

 

(5)

 

 

 

 

 

 

lim

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

г—>

 

 

 

 

 

 

 

 

 

г11,

 

 

 

 

 

 

V(г, п) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вытекает

равенство

 

 

фп (р))

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

(/ (р),

 

 

 

 

 

 

 

 

 

 

п

 

Ѵ п

 

+ Ѵа

Р

,I (р),

п =

0,

 

1, 2, . . .

 

 

V

 

Р Фп (р) =

 

 

 

 

 

 

 

 

 

 

 

 

п)

(р) обозначены полиномы Лежандра.

Применяя

Через

 

 

к

(г,

 

обратное преобразование Лежандра,

мы полу­

чаем формальное

решение

 

 

 

 

оо

(/> Фп) гпФп (cos Ѳ).

(6)

 

 

V (г, Ѳ) =

V (г, arccos р) =

 

2

Теперь мы докажем,

 

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

что функция (6) — обычная и удов­

 

 

 

 

 

 

 

D fD l

 

р

 

q

 

 

 

 

 

 

 

 

 

 

 

летворяет уравнению (1) в обычном смысле. Формально

применяя

оператор

 

 

 

 

 

( ,

 

 

=

0, 1, 2, . . .)

под

зна­

ком интеграла в (6), приходим к выражению

 

 

 

D?Dfa> =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

2

 

(/» фи) Л (и—1)... (и—р - И ) / я

+

 

(І-О- (7)

T L ~ p

Из теоремы 9.6.1 вытекает, что (/, фп) имеет медленный рост; это означает, что существует постоянная С и целое число к, для которых

I (/. Фп) I < Спк, п = 1, 2, 3, . . .

Кроме того, известно, что

I D l Р п (р) I < и»«, - 1 < р < 1, » = 1, 2, 3, . . .

343

(см. Черчилл [2], стр.ап211—212)

. Поэтому ряд (7) сходится

равномерно в области

0

а < 1 ,

— 1

< 1 ,

так как функция

экспоненциально

убывает

при

n -V оа, а другие множители ограничены полиномами от п. Отсюда мы видим, что (6) — обычная, причем гладкая, функция и что оператор Лапласа можно применить к (6) почленно. Поскольку функция гпф„(р) удовлетворяет урав­ нению (1) в обычном смысле, то, как легко показать, то же самое справедливо и относительно функции (6).

Всилу сказанного мы можем перейти к пределу при

г— 0 под знаком суммирования и прийти к заключению,

что выражение (6)

стремится к (/, ф0)/|/"2 равномерно на

О

0 ^ я .

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

выражение (6) остается огра­

ниченным в

окрестности началаЛ координат.

.

ср — произ­

 

Проверим, наконецЛ ,

, условие

(2)Ь 2. Пусть

вольный элемент пространства

Согласно

лемме

9.3.3

ряд (6) сходится в

 

а потому и в

 

(/), для любого фик­

сированного

г

1.

Следовательно, мы можем вычислить

его

скалярное произведение

с ср почленно и написать

 

(у (г, arccos р), ср (р)) = 2

(/> Фп)

гп

(Фи, ф),

 

г < 1.

(8)

 

 

 

 

 

 

 

 

71— 0

 

 

 

 

 

 

Функция (/, ф„) имеет медленный рост (теоремы 9.5.1 и 9.6.1), а (фп, ср) быстро убывает (леммы 9.3.2 и 9.3.3). Поэтому правая часть равенства (8) равномерно сходится на 0 < г < 1, и мы можем перейти к пределу при г 1—0 под знаком суммирования. Тогда

С О

lim (V (г, arccos р), ср (р)) =

2 (/, Фп) (Фп,

ф) =

(/, ф),

 

,'-*1—0

 

п=0

 

 

 

 

 

 

чтокотораяи доказываетудовлетворяетформулудифференциальному(2).

уравнению (1) во внеш­

В а д а ч а 9.11.1. Найти

обычную

функцию

ѵ (г,

arccos

р),

ности единичной сферы (т.

е.

при 1 < г < о о ,

1

< р <

1),

поточечно сходится к нулю при г —> оо и сходится

в

 

к / 6Е с&'

при г —>1 + 0. Здесь через

 

снова

обозначено

пространство

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

З а д а ч а 9.11.2. Рассмотрим задачу о распространении тепла в одномерном неоднородном стержне. Этот процесс описывает­ ся дифференциальным уравнением

[aDx (1 — X2) D x + (р.т +

ѵ) D x] V (х , С) = cD tv (х, t),

(9)

— 1 < ж < 1,

0 < г < о о ,

 

344

1'де я, с,

р и V действительнее постоянные, прячем а > О, с > Ü.

Через V

(X, і) обозначена температура в стержне. Все поверхности

стержня считаются обсолютно изолированными. Отметим, что функция а (1 — х3) представляет собой теплопроводность стерж­ ня, с — постоянная, зависящая от плотности и удельной теплоем­ кости стержня, а слагаемое (рг + ѵ) D xv описывает распределен­ ный источник тепла внутри стержня. Показать, что дифференциаль­

ное уравнение (9)

можно привести к виду

 

 

 

 

 

 

а 91л и (X , I) = с Diu (х,

t),

 

 

где 91а; — дифференциальный оператор

(по

х),

соответствующий

системе 4 (п.

9.8)

для преобразования Якоби, и =

V w ѵ, w (х) =

= (1 — х)а (1 +

a;)ß. Какие ограничения на р и" ѵ

налагают усло­

вия а > - 1

и Р > - 1 ?

 

 

 

 

Предположим,

далее, что при / -> -(-

0 функция ѵ (х , t) сходится

в 3)' (/), I = (— 1,

1) к g {x)/V w (X), где g Er •?#',

Л-' — соответст­

вующее пространство обобщенных функций для

 

преобразования

Якобн.

Найти

формально решение этой задачи Коши. (Этот пример

взят из

работы Дебната [3].)

 

 

 

 

9.12, Применение первой формы конечного преобразования Ганкеля: задача Дирихле для полубесконечного цилиндра

Конечное преобразование Ганкеля оказывается полезным при решении различных граничных задач для областей с цилиндрическими границами. В этом и следующем пунктах мы приведем два соответствующих примера; первый посвящен задаче Дирихле.

Мы хотим найти решение уравнения ЛапласаRво вну­

тренней части

R

полубесконечного цилиндра единичного

радиуса. В цилиндрических координатах область

можно

описать

следующим

образом:

 

 

 

 

< 2

 

 

 

 

R

=

{(г,

Ѳ, z):

0

< г

< 1 ,

0 < 0 < 2 л ,

0

 

 

 

< оо}.

Предполагая,

 

что

искомая

функция

ѵ

не зависит от Ѳ

(т. Уе.г v(r,

V

(г,

z))

и

 

вводя

новую функцию

 

и

(г,

z ) —

=

 

V

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z),

 

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

где,

как

и всегда,

R о,ги

л-

DUi

=

0,

 

 

 

 

 

 

 

(1)

0іГ

=

 

 

 

 

 

 

 

 

 

Мы

наложим

 

 

S

r-'-''DrrDTr-'\

 

 

условия.

 

 

также

следующие граничные

Л '.

1) Если z

—> + 0 . то

и

(г, z) сходится в

Л ’

к

g

(г) SE

Здесь

Л '

обозначает пространство

обобщенных

 

функций,

соответствующее

системе 5а

(см.

п. 9.8).

 

 

 

 

 

 

345

па

2)

Если z —ѵ оо, то

и {г,

z) равномерно сходится it нулю

 

0 <

г

СZ

1.

 

 

 

 

 

 

 

Z

 

 

 

 

3)

Если

т

—>1 — 0,

то

и

(г, z)

 

равномерно сходится

к нулю на

г

 

 

z <

со

при любом

 

>

 

0.

Z

4)

Если

 

-*

+ 0 ,

то

и

(г, z) =

О ( У

г)

равномерно на

 

z < с»

при любом

Z ]> 0.

 

 

 

 

 

Отметим, что эта задача полностью аналогична рассмот­ ренной в п. 5.8. Однако там областью изменения был

полубесконечный интервал

(0, оо),

и мы использовали

преобразование

Ганкеля

(не конечное) нулевого по­

рядка. Теперь

же г изменяется в

конечном интервале

(0, 1) и мы используем первую форму конечного преоб­

разования Ганкеля

U

нулевого

 

порядка.

х

 

г

 

Уо,пМы

Упвведем-

 

обозначения, несколько отличающиеся от

указанных для системы 5я (п.

 

9.8),

а именно

 

=

 

и

 

=

 

Формально

 

примененное

преобразование

U

переводит (1)

 

в уравнение

 

 

= 0,

 

 

 

 

(2)

 

U

— U

 

п

 

y lU

 

D\U

 

 

 

 

где

(

, z) —

(и (г,+

 

 

 

 

 

 

 

 

 

 

 

 

 

z), ф„

 

(г)).

 

 

 

 

 

 

В силу граничного условия 2) мы выбираем в качестве

решения уравнения

(2)

функцию

 

 

 

 

 

 

 

 

 

 

 

U

(

п

, z) =

A

(n)e~VnZ,

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где А (п) не зависит от z. Используя преобразование гра­ ничного условия 1), получаем

A (n) = G (и) ^ (g, ф„).

И наконец, вычисляя обратное преобразование для функ­ ции (3), мы получаем желаемое решение

и {г, z ) = 2 G (,l) e ""ЧЛ»').

(4)

П—1

 

где

У 2rJo(V .

J M

Удостоверимся теперь в том, что и (г, z) — действительно решение уравнения (1). Прежде всего, справедливы фор­ мулы

У п ~ Я ( п

п оо,

(5)

346

и

І А Ы Г Ч ^ .

» = 1 , 2 , 3 _______

где К — некоторая постоянная (см. Янке, Эмдѳ и Леш [1]). Кроме того, из соотношения (5) и теорем 9.5.1 и 9.6.1 следует, что G (п) — О (пр) при п —> оо для некоторого целого числа р , тогда как е~упг = О (e~nnZ) равномерно на Z z < оо при любом фиксированном Z > 0. Заметим, наконец, что

D rJ О(УпГ) = —IJnJl (УпГ),

D p -Ji (УпГ) = УпГJo (УпГ).

а функции

J

о (г/„г) и

J x (упг)

равномерно ограничены

на

0 < г <

1,

 

» =

1, 2,

3, . . . В силу

всего сказанного

ряд (4)

равномерно сходится при 0 < г

< 1 ,

Z ^

z с

оо

Z

0).

Это утверждение останется верным и после при­

менения

любого

из

операторов

D , r ~

D ,rD rr~'^,

D .,

D i

 

 

 

 

 

в формуле (4) под знаком суммирования. Таким образом,

можно применить оператор

S

0, +

D i

к (4) почленно,

 

 

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

Установленная равномерная сходимость ряда (4) по­ зволяет нам перейти к пределу под знаком суммирования при z со, г -> 1 — 0 или г —► -{-0; в действительности это можно сделать даже для ряда, состоящего из абсолют­ ных величин членов ряда (4). Отсюда непосредственно

следует выполнимость

граничных

условий

1),

3)

и

4)

соответственно.

 

 

 

 

 

 

 

 

Л ,

Наконец, в силу леммы 9.3.3 и наших предыдущих

рассуждений

правая часть

равенства (4) сходится

в

 

а поэтому

и

в

Ь 2 {1),

при

всех

фиксированных

z

>

0.

Мы можем,

следовательно,

вычислить скалярное

произ­

ведение

и

(г,

z)

с любой функцией

ср (г) 6Е

Л

почленно

и

получить

формулу

со

 

 

 

 

 

 

 

-

 

 

 

е~Ѵп\(Фп, ф).

 

 

(6)

 

(и {г, z), ф (г))= 2

 

 

 

 

 

 

 

П=і . . .

 

 

 

 

 

Но функция G (и) имеет медленный рост, в то время как (Фш ф) быстро убывает при п -> оо. Поэтому ряд в (6) сходится равномерно на 0. z <Со°, и мы можем перейти

Ö47

к пределу под знаком суммирования при z -> -|-0; тогда

со

(г, Z), Ф (г)) - » 2 G («) (Фи. Ф) = (g. ф)-

7Ѵ=1

Таким образом, граничное условие 1) выполнено. Этим завершается проверка того, что (4) действительно является

решениемв конечной. области вида

 

дифференциальное уравнение

(!)

 

 

З а д а ч а

9.12.1.

 

Решить

 

 

 

 

{(г,*):

0 <

г < 1,

0 <

z <

с>

 

 

 

прп

следующих граничных условиях:

 

 

 

 

 

 

 

1)

если z —> + 0 , то

и (г, z)

сходится в cé ' к

(г) е

 

 

 

2)

если z —» с — 0,

то и (г,

z) сходится в J k '

к g2 (г) е

 

 

3)

если г —* 1 — 0,

то и (г,

z) сходится к нулю равномерно на

 

<

z < Zn, если 0 <

Z x <

Zn < с;

 

 

 

 

 

 

 

 

4)

если г —» +

0, то и (г, z) = О ( V г ) равномерно на Z x < z <

<

Zn,

если 0 <

Zj < Z2 < с.

 

 

 

 

 

 

 

 

 

 

Здесь j b ' обозначает то же самое постранство обобщенных функ­

ций,

что п ранее в этом пункте. Проверьте полученный ответ.

 

 

 

З а д а ч а

9.12.2.

Получить формально

решение следующей

задачи. Пусть Я обозначает область, заключенную

между двумя

полубесконечнымп концентрическими цилиндрами с

радиусами а

и 6; в цилиндрических координатах эта область имеет вид

 

Я

= {(г , Ѳ, z):a <

г < 6, 0 <

а <

Ь <

оо,

0 <

0 <

2л,

0 < z <

оо).

Найти обычную функцию и =

и (г,

 

 

 

 

 

 

 

 

 

 

 

 

 

z), которая не зависит от 0 и

удовлетворяет в области Я дифференциальному уравнению (1).

Пусть

— пространство обобщенных функций, соответствующее

системе 5с для третьей формы конечного преобразования Ганкеля,

где |х = 0. Мы предполагаем, что функция и удовлетворяет следую­

щим граничным условиям:

 

 

 

 

 

 

 

 

 

 

1)

если z —* + 0, то и (г,

z)

сходится в c é ' к j

(г) Е

 

а

 

2)

если z —» со, то

 

и (г,

z)

равномерно

сходится к нулю

на

<

г

< 6;

 

 

оо

 

 

 

 

 

 

 

 

 

 

3)

если г —* Ъ — 0 или г —» а +

0, то и (г, z) равномерно схо­

дится к нулю на Z ^ z <

 

прп Z >

0.

 

 

 

 

 

9.13.

Применение второй формы

 

 

 

 

 

конечного преобразования Ганкеля:

тепловой поток

 

вбесконечном цилиндре с условием излучения

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

348

уравнению

D ;v

-f r_1Z)ri; =

D tv,

c t

V —

v{r,

z),

(l)

 

0

C r

< 1 ,

0

<

oo.

 

 

 

 

 

 

 

В качестве начального условия мы потребуем, чтобы при t —' -|-0 функция V (г, t) сходилась в некотором смысле к обобщенной функции / (г); в качестве граничного усло­ вия мы потребуем, чтобы при любом фиксированном ' > 0

lim [D Tv + au] = 0,

(2)

1-»1—о

 

где а — действительная положительная постоянная. Если V обозначает температуру внутри цилиндра, то условие (2) соответствует излучению тепла поверхностью цилиндра.

Для решения этой задачи можно использовать вторую форму конечного преобразования Ганкеля нулевого по­

рядка;^ при этом мы прежде всего положим и = )f r v ,

b = V r U

а постоянную а, входящую в соотношение (2),

будем считать равной

а,

входящей в форму­

Лпостоянной'

лу (3) п. 9.8. Мы будем также предполагать, что

g

принад­

лежит пространству

обобщенных

функций,

 

соответ­

ствующему системе 56 (см. п. 9.8). Дифференциальное

уравнение принимает теперь вид

 

 

 

 

S 0tru = D tu,

u =

u{r, t ) ,

 

0

< r

< 1 ,

0 < г <

оо, (3)

а начальноеи

 

 

 

 

 

 

 

условие переходит в

 

t

 

 

 

(г,

t)

g

(г )

в

Л '

при

->• + 0 .

(4)

Граничное условие

превращается в

 

 

 

 

lim

[D ru

+

(а — г/2) и]

=

0.

(5)

 

г->1—О

 

 

 

 

 

 

 

 

 

 

 

При этом подразумевается, что t имеет любое фиксирован­ ное положительное значение, и предел на 0 < ; t <[ оо по­ нимается в поточечном смысле.

Пусть теперь U обозначает преобразование, соответ­ ствующее системе 5Ъ (п. 9.8); положим

U (п, t) = U и (г, t), п — 1, 2, 3, . . .

Оператор U преобразует дифференциальное уравнение (3) к виду

—&.U (п, t) = D tU (re, <), 0 < t < оо,

(6)

349

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