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

книги / Операционное исчисление и обобщенные ряды Лагерра

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
9.58 Mб
Скачать

f - f - } ^ ( 0 - ^ 5 5 r { [ l n 2 + W(n)p + 3[В Д + I n z i m ’Ü ) -

Ф'(»)] + W{n)} = A n,i{z).

Замечание. Применяя теорему сдвига (при а > 0 ) к выра­ жению (1), получаем

{(f - a)-» lnfi(f - а)} r,(t - «) - «Г«* А а, н (г).

(4.3.3)

Этим соотношением воспользуемся ниже.

п. 2. Функции с неинтегрируемыми степенными особен­ ностями в произвольной точке. Очевидно, (ЕеЯ>—1):

а

txvi(О -/) (а — t) ■+■ J етЧч*- dt = Ф (я, X, z), а > 0.

(4.3.4)

0

Так как изображение Ф(а, Я, z) аналитично для всех к ф —1,

—2 , ... (в точках к = —1, —2 , ... эта функция имеет про­ стые полюсы), то аналитическим продолжением соответст­ вия (4) по к имеем уже для всех к ф —1, —2 , ...

{*}

(«. К г), о > 0.

<4.3.5)

Аналогичным приемом, изложенным в § 1, разложим функцию Ф(а, Я, г) в ряд Лорана по степеням к + п . В резуль­ тате получим

Ф(а, X, г) =

 

(—1)п

^

 

(Х+л)*

(4.3.6)

(п—i)!(À+»)z"

+ 2d -Р*0 («>г» 1п«) • —

---- ,

 

 

 

 

 

й-»0

 

 

 

где

 

 

 

 

Р(ь ](а, г, 1па) =

 

 

 

 

 

 

 

 

 

 

( — 1 ) “

1 f

n l n k + 1 <*

- щ { — А

я ’к '

(

s l l n ft

___

n!z»

L

A+l

n , \ * /

(U

{i-l)\ Z à\ sj

(u-0*+1

a

k

/ft\

 

 

 

.....1; 2,...,2,

n + 1;

-

~ " T

2

( S

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.3.7)

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

 

- 0

=

+

2

{!?} £т г -]ч(*М«—»).

 

 

 

 

 

 

(4.3.8).,

Сравнивая (6) и (8), учитывая при этом (5), находим:

 

{f—eln*t},>j(f)Tj(a—t)

Pj^(a, z, lna),

(4.3.9)

где

(a, z, lna) определяется по формуле (7).

 

Аналогично тому, как было получено соответствие (9),

устанавливаем

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

следующих соответствий

(ниже к = 0 , 1 , 2 , . . . ;

п=1, 2 , ... )

 

 

 

{(t -

a)- "In* (t -

a)} 7j(f)Y|(«-

f) -

e -“/2 Й 0 (a, A

(4.3.10)

здесь и ниже ln(f—a) понимается в смысле главного значе­ ния

{(а - 0 _в1п*(а -

0} < t) т)(а -

г) -*-e~*l* Q<»>(a, z),

(4.3.11)

где

 

 

 

 

Ф &

z) =

- Р [ п) (a,

z, 0), (а, г) =

 

 

= ( - 1 ) " - 1 Р Р [ « , г , 1 п ( а е « ) ] .

(4.3.12)

 

 

Суммируя попарно соответствия (3), (10) и (3), (11), по­

лучаем следующие соответствия:

 

 

{(t - a)~» In* (t -

a)}vj(f) с-*!-[An, hW + Й и)(а, *)],

(4.3.13)

{\t — «|~*ln* \t —a|}nj(f)-s- e~alz[An, * (z)+ Q(ftn)(a, z)].

(4.3.14)

n. 3. Регуляризация некоторых расходящихся интегра­ лов. В силу обобщенного интеграла Лапласа — Карсона со­ ответствие (9) принимает вид

a

 

Р р (Я, г, lna) = [ €Г*1‘ [t~n lnkt) dt.

(4.3.15)

b

 

Полагая здесь z = l, получим регуляризованное

значение

следующего интеграла :

 

а

 

Jе~* {t_nlnfe£} dt = P (An)(«, 1, lna)

(4.3.16)

О

 

(при л= 1 второе слагаемое в (7) должно быть положено рав­ ным нулю).

Переходя к пределу в обеих частях (15) при z-*-+ оо, по­ лучаем

*+1 lnft+1 a

при

n — 1

hi_^

In8a

(4.3.17)

с”- 1 2 d

s!(n-l)*-8+1 при

71 = 2 ,3 ,...

s=0

v

 

Совершенно аналогично с помощью соответствия (10) вычисляются следующие регуляризованные значения инте­ гралов:

Jе~* {(t — a)-» In* (t — a)} dt =

e~aQlbn) (a, 1)

0

 

 

 

(4.3.18)

k = 0,1 ,

2 , . . .

 

 

 

_ è ï ln№(ae‘K)’

71 =

1

 

 

( - 1)n _ 1A !

V

,

» n

•■

n-l

£d~

 

e-o

sl(n—1) '

 

(4.3.19)

 

 

 

 

Замечание. Равенства (18) и (19) тажже могут быть по­ лучены из (16) и .(17) соответствующей заменой переменных интегрирования, причем при выводе (18) необходимо учесть второе из равенств (12).

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

Как известно, понятие конечной части расходящегося интеграла было впервые введено Коши [132], назвавшим его «необычным интегра­ лом». Он также указал на возможность дифференцирования и интегри­ рования по параметру таких интегралов. Конечная часть расходящегося интеграла была использована в теории дифференциальных уравнений к частных производных Р. Адемаром [114].

Однако большая заслуга принадлежит Ж. Адамару [157, 158], который не только изучил свойства введенного им независимо от Коши понятия конечной части, но и, распространив эту концепцию на кратные интегралы, применил эти понятия для решения задачи Коши уравнения гиперболического типа второго порядка. Обобщение его идеи было осуществлено Ван-дер-Корпутом [218].

Впоследствии появился ряд работ, где упомянутая концепция была исследована с различных точек зрения. Так, в теории распределений Шварца [202] конечная часть трактуется как разновидность псевдофунк­ ции и играет связующую роль между обычными функциями и распреде­ лениями.

Понятия конечной части и логарифмической части расходящихся интегралов были применены к решению задачи Коши для уравнений вто­ рого порядка всех трех типов и системы уравнений в частных производ­ ных в работе [128].

Лавуан [172, 173] и Рюхе [192] с помощью понятия конечной части расширили операционное исчисление.

Джилли [153—155] использовал понятие конечной части для регу­ ляризации степенных особенностей любого отрицательного, но нецелогопорядка при обобщении преобразования Лапласа.

Следует также указать на исследование Бохме [126], где дана регу­ ляризация особенностей типа {£*1п£}, здесь г может быть любым, в том числе целым отрицательным числом. При этом он опирается на работу Бутсера [129], в которой операторное исчисление Микусинского [76]

применено к исследованию конечной части расходящихся

интегралов-

типа свертки.

|т)(г), t - n \ n 4

О.

о. вида (ï—a)“n lnfe (f—a)ij(*), |t—a| -nlnft |t—a

(a—i)

впервые были изучены Ж. Н&урзбаевым [80].

 

Г л а в а 5

ЧИСЛЕННОЕ ОБРАЩЕНИЕ ПРЕОБРАЗОВАНИЯ ЛАПЛАСА МЕТОДАМИ ГАРМОНИЧЕСКОГО АНАЛИЗА

Ниже речь будет идти о приближенном (численном) ре­ шении интегрального уравнения Фредгольма первого рода

(1 )

0

где F(p) — заданная функция; /(f) — искомая функция. Известно [54, 64, 120], что указанная задача относится

к классу некорректных задач. Одной из существенных ха­ рактеристик этой некорректности является следующая осо­ бенность: сколь угодно малые погрешности, вносимые в вы­ числение значений функции F(p), могут существенно влиять на значения вычисляемой функции /(f). Этот факт можно проиллюстрировать следующим простым примером.

Пусть решается задача обращения (1). Известно,

Нетрудно видеть, что для любого е > 0 существует такое А, что при а > А х

Таким образом, изображения F(p) и i|)(p)=Г(Р) + рз_|_дг

численно отличаются друг от друга на величину, модуль которой не превышает е. Тогда как соответствующие им ори-

гиналы разнятся величиной |flsinaï [, которая при а > А характеризуется высокочастотными колебаниями.

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

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

Пусть изображение F{p) разлагается в сумму вида

*ЧР) = 2ь «*Ф*(Р).

(2)

Пусть далее известны оригиналы <pft(i), отвечающие изобра­ жениям Ф*Ср). Тогда ряд (2) порождает ряд

2

(*)•

(3)

к

 

 

Естественно, при этом возникают вопросы:

а) в каких случаях и в какой топологии ряд (3) представ­ ляет оригинал /(t), отвечающий изображению F(p);

б) каким образом вычислить коэффициенты а* ряда (3) (или что то же ряда (2)).

Для задачи численного обращения важное значение также приобретает ответ на такие вопросы, как: скорость сходимости ряда (3); в каких случаях появляется возмож­ ность улучшить сходимость ряда (3); с какой точностью могут быть определены коэффициенты а к ; какова степень влияния точности вычисления коэффициентов а * на общее представление функции f(t) рядом (3).

Очевидно, что вид разложения (3) зависит от выбора функций Ф*(р). В данной главе используется простейшая ситуация, когда в качестве <p*(f) избираются функции Лагерра 9k (0 = e-t;2 L k (t). Таким образом, здесь ряды Лагерра выступают как аппарат численного обращения преобразо­ вания Лапласа. Для вычисления коэффициентов a*(fc=0, 1, 2,...) ряда Лагерра (называемых иногда лагерровским спектром) используются формулы, являющиеся аналогами формул тригонометрического интерполирования. Этот факт позволяет широко применять методы гармонического ана-

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

§ 1. Обзор важнейших задач тригонометрического интерполирования

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

Здесь в отличие от того, как это делается в известных работах по гармоническому анализу, алгоритмы сопостав­ ляются между собой на основе связи между «дискретными»

и «интегральными» коэффициентами Фурье.

и I

Пусть функция f(x) определена на отрезке

— некоторое множество индексов.

многочлен

 

Задача

tr: Построить тригонометрический

у{х) = 2

о*, п е' 1кх, такой, что y(xk)= f(x k) (kQl),

где

AG I

 

 

xk (kQ l) — узлы, распределенные на [— 1,1].

 

Ниже приведены наиболее распространенные варианты задачи tr, обусловленные специальным выбором узлов.

Задача t r

1 п

t r 2

t r 3

 

 

 

N — число узлов

N = 2 n + 1

JV-2H+1

N — 2 n

Вид узлов

x k ~ ± -

2 k

X k -----—

 

п

Хк = 1 Г

n

Множество измене­

 

—л+1 <А <л иля

ния к

 

 

—л<А< л—1

В задаче tr{ концевые точки промежутка [—1, 1 ] вклю­ чаются в состав узлов, в задаче tr2— исключаются, в зада­ че fгз — включается одна из них.

Основой приведенной классификации служат два отно­ шения ортогональности :

Задача try :

j 2 в 71 £

п == I®*

^п kQi

Ц, если s=q»2n-\-m

J

Задачи tr2, tr3

 

 

 

 

1 2

QTi~iï g~2ri^~ _если s Ф qN -f- m

 

N bel

 

(5.1.2)

 

 

Ц, если s = g/V -f- m

 

Здесь g, s — любые целые числа; me/; символ 2 " озна­

чает, что

два крайних

члена суммы взяты с весом, рав­

ным '/г*

 

 

 

аап

Используя (1) и (2), получаем искомые коэффициенты

задачи tr :

Задачи tr2, tr3

Задача tri

 

 

2

лк

л к

 

 

" .

~ТЛП

± \ .

-

с“" ш

fk&

 

a* * = * m fhe

Периодичность по переменной s соотношений (1), (2) с периодом 2п и N соответственно позволяет установить связь коэффициентов a sn задач t r i,2, г с коэффициентами аа раз­ ложения /(г) в ряд Фурье

 

 

 

 

 

 

00

 

 

 

 

 

 

 

f(x) =

2 ake«ikx.

 

(5.1.3)

 

 

 

 

 

 

*— »

 

 

Задача tri

 

 

 

 

 

 

 

X

чт-ч/у

 

л т

x i

1

.km .km

 

J _ V "

e

- * i —

 

 

 

 

a sn =

2n2â

я

 

2 ? к 9 2

e ^ T e ^ T .

 

п

 

** n

2

 

 

me/l

 

 

«

A=—oo

jnel

 

Отсюда в силу (1),

 

 

 

 

 

 

 

 

 

(QT2B+S ~b Я—r2n+s )

(sG-0»

(5.1.4)

 

 

 

 

 

r-1

 

 

 

Задачи tr2tr3

 

 

 

 

 

Аналогично проверяется справедливость

 

 

 

 

atn — ae+ 2

+« “b Q—TN +«) (sG-0.

(5.1.5)

 

 

 

 

 

r-l

 

 

 

Соотношения (4) и (5) устанавливают связь между «дискрет­ ными» и «интегральными» коэффициентами Фурье и пред­ ставляют удобный аппарат для анализа вопросов точности тригонометрической аппроксимации.

В таблице 1 приведена более детальная классификация задач tr 1,2в связи с переходом от «комплексной формы» записи тригонометрических многочленов к «вещественной форме». Эта таблица составлялась в предположении, что ряд (3) приведен к виду

СО СО

f(x) = 2

c*cosfejt* + 2 dksink~x,

(5.1.6)

k = 0

ft=l

 

где Cft = a * + <x_fc; dh = i{a k —

аотвечающий ему тригонометрический многочлен приведен

квиду

П

71—1

 

у(х) = 2 сЬпcosArae + 2

(5.1.7)

ft=0

ft=l

 

Здесь символ S 7означает, что

первое слагаемое берется с

весом Уз-

 

 

Таблица 1 в каждом конкретном случае содержит так* же связи между «дискретными» и «интегральными» коэф­ фициентами Фурье, при этом для сокращения записи вве­

дены обозначения

во

оо

Ск Н = 2 ( с г Я + к + C r N - k ) ,

= 2 № r N + k d r N - к ) .

7*=1

Г*=1

Замечания. 1 . К. Ланцош [70] рекомендует видоизме* нить задачу tr\ и строить интерполяционный многочлен в форме (задача tr'\) :

Пгг л—1

у(х) = 2 ckncosvkx + 2 ft=0 *=1

Практическая ценность указанной модификации находит простое объяснение, если воспользоваться связью между «дискретными» и «интегральными» коэффициентами Фурье. Для задачи tri имеем с„„ = с„ + о+ 2в=2(с„ + с3|»+е5п+ ... ) .

Отсюда следует, что величина ~спп приближает величину

сп с точностью до порядка малости коэффициента с Ьп. Этот факт собственно объясняет целесообразность введения в за­ даче tri весового множителя ‘/г при коэффициенте с„п (т. е. модификации по.К. Лагацошу).

2. В нЬследующем под задачей tri будем понимать зада­ чу tr'u Это соглашение, кроме упомянутого в замечании 1 достоинства, полезно в следующем отношении. Бели в

Г Г ' ^ 9 9

Индекс

Ограниче­

Формулы

задачи

ние

 

tTiC

tris

tr2

Узлы:

k_

* k = - n

/■(*)—ве- щесTвенная функция

Узлы: h

*А = —

n —n<k<n

/(*)—ве- щественная четная функция

Узлы:

*А = ' n

— n*Ck<n

/(ас)—ве­ щественная нечетная функция

Узлы:

2k

Xk = —-

2п+1

—n<k<n

/(ас)—ве­ щественная функция

п—1

У ( * ) в

2

Ckn C037:kx + 2

<*Ал sinnkx

 

А=0

 

 

А=0

 

С А п = “

2

 

( / r o + / - m ) C O S

 

 

771=0

 

 

 

 

 

l ”—1

 

 

 

*

d ftn =

 

ifm —f —m) Sin^m —

 

177=1

 

 

 

 

c*„=cfc + c+2ll; d*B=d*+o^2n

ÿ(* )= 2

Ckn еозг^-х

 

 

 

k=o

 

 

 

 

 

2 v i"

 

*

 

Cbn =

~ n 2 à

f m

C0S1!m 'й

 

 

m =0

 

 

 

 

сАл =» Ck +

2n

 

 

 

 

л — 1

 

 

 

 

 

2/(*)^2

^Алs’in7C**

 

 

 

A=1

 

 

 

 

d*n== n ^ f m

s i n n m

 

 

ГО=1

 

 

 

 

iftn =

rfft + a^2n

 

 

 

71

 

 

 

 

 

ÿ(* )= 2

cftn cosrcfcc+^]

 

 

*=o

 

 

fc=l

 

 

2

 

 

s

77Î&

Ckn =

2 Й + Г ^

^ m + f—m) cos2n 2TtJr i

 

 

 

m=Q

 

 

d k n = 2n + l

m=1

 

m) sin2n «j^+T

 

 

 

 

 

 

 

 

_+ '

 

^Ал — ^А+°^, 2л+1

САл — ^A +

2л+1»

1D0

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