книги из ГПНТБ / Амербаев, В. М. Операционное исчисление и обобщенные ряды Лагерра
.pdfСравнивая приближенную формулу (9) с точной формулой (5), заключаем, что погрешность Й(7) формулы трапеции
7+i<* |
|
|
2 Д7+1т ^ |
|
2nk |
|
Ш I |
ePtF(p)dp — -у- |
|
||||
Т + |
* Т “ + 2(0 |
|||||
Y—1с |
|
|
к——со |
|
|
|
|
|
оо |
|
|
|
|
имеет вид |
2(1) = |
— 2 |
/(t + |
kl). |
|
|
Поскольку |/(1)| |
к=1 |
|
|
|
|
|
(у>уо), то |
|
|
||||
|
|
| 2(0 | < |
Ме**- |
е—1(т—То) |
|
(5.6.10) |
Неравенство (10) оценивает величину погрешности фор мулы трапеции в зависимости от выбора параметра I, т. е. в зависимости от выбора шага h = 2л/7 формулы трапеции.
§ 7. «Смешанные» производящие функции
Из вышеизложенного со всей ясностью вытекает, что многочлены Лагерра играют особую роль в операционном исчислении и они могут с успехом применяться в самых различных целях. Ниже будет показано, что многочлены Лагерра наиболее удобны при выводе так называемых «смешанных» производящих функций.
Вначале вкратце покажем вывод формулы для одной «смешанной» производящей функции, остальные приведем без выкладок.
Используя производящую функцию [21, 93]
* |
-L |
2 |
P n(x)hn = ( l - 2 x h + h2) з , | * | < 1 |
для многочленов Лежандра Р„(л;), находим равенство
оо |
3 |
2 (2га + 1)hnPn (ж) = (1 - Л2) (1 - 2xh + k2)T .
п—0
Отсюда с помощью замены й = о>(1—р)/р, где р — операци онная переменная, получаем:
V I |
/1—р \ п |
р[(1—ма)рг+2о)р—0)8] |
(5.7.1) |
^ ( 2 г а + |
1)а> ^ р J — |
lip2_ 2ta{x+ti3)p+u>^ i z . |
в—0
130
Выше и в дальнейшем в этом параграфе использовано обозначение
£ = £(л;, со)= 1 -j- 2ясо + со2. |
(5.7.2) |
Применяя к равенству (1) обратное преобразование Лап ласа по р, находим «смешанную» производящую функцию для многочленов Лежандра и Лагерра:
00 |
|
|
|
|
|
1 |
|
|
|
|
|
2 |
(2п + l)u>nLn(t) Рп (х) = -г? - ei {[1 — со2 -f 2со (tx + |
||||||||||
л-0 |
|
|
|
|
|
У |
|
|
|
|
(5.7.3) |
|
|
+ 2 b V 1—x2)]Jо(о)— 28(1 - со2) } |
(8), |
||||||||
|
|
|
|||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
7 = |
7(t, л, со)= М*.+-т).; |
Ь= |
S(t, х, со) = |
* . (5.7.4) |
|||||||
Аналогично |
для |
многочленов |
Гегенбауэра P„(v+i/z) и |
||||||||
обобщенных многочленов Лагерра имеем: |
|
|
|
||||||||
" |
|
|
|
|
|
|
2vr(o+l) el J v(5) |
|
|||
2 r(A+2v+l) P(b+ll?)(x )H p (t)= |
jv + l/2 5vr(2v+1) =(М*. X , со). |
||||||||||
ft=0 |
|
|
|
|
|
|
|
|
|
|
(5.7.5) |
|
|
|
|
|
|
|
|
|
|
|
|
Отсюда, |
в частности при v = 0, получаем формулу для мно |
||||||||||
гочленов Лежандра и Лагерра: |
|
|
|
|
|
||||||
|
00 |
соkPk(x)Lk(t) = |
1 |
elJ о(&) — М*, х, со). |
(5.7.6) |
||||||
|
2 |
y f |
|||||||||
|
ft—о |
|
|
|
' |
|
|
|
|
|
|
Ниже в целях краткости принята условная запись типа |
|||||||||||
|
|
|
^ - |
А ^ |
а |
+ |
Аг ап Р „ |
cos |
|
||
|
|
|
тс |
|
\cos |
|
|
|
|||
которую |
надо |
трактовать |
как |
объединение двух |
формул |
||||||
cps = |
A (since -{- AtsinP), |
<ре = |
A (cosa + |
AiCosP). |
|
||||||
Для многочленов Лагерра и синусов (косинусов) имеем
2 «>*£*(t)sin |
kQ = (t, 0, со), |
(5.7.7) |
|
*—о |
cos |
|
|
131'
(t, е, |
со) = |
eTe[rtn §8 + |
0)sin (8e - 0)1 , |
(5.7.8) |
||
|
|
*“9 |
[COS |
cos |
J |
|
где $e = £(cos0, w), |
т е = |
y(£, |
cos0 , w), |
89 = |
8(t, cos0, |
«). |
Для многочленов Чебышева первого и второго типов:
2 шкТк(х) Lk (f) = |
et[«j(l — *2)sin8 + (1 + лгсо) cos8], (5.7.9) |
fc=0 |
|
Разложения (3)—(10) могут рассматриваться как «сме шанные» производящие функции для многочленов Лагерра. Они позволяют установить интегральные представления многочленов Лагерра, Гегенбауэра, Лежандра, а также ко синусов и синусов.
Действительно, умножая обе части разложения (5) на (1—a:2'v.p^+i/2) интегрируя полученное равенство по х
в пределах от —1 до 1 и учитывая ортогональность многочленрв Гегенбауэра, имеем:
1
-1
(5.7.11)
/ |
(2 у + 1 ) „ Г ( у + 1 ) |
<»nL l* 4 t) |
( л - И + 1 /2 ) Г ( ^ + 1 /2 ) |
r ( n + 2 v + l ) * |
|
Таким же образом из разложения (5) находим:
00
о
В частности, при v = 0 соответственно получим:
1
- 1
132
00
J e~*Ut, x,u>)Ln(t)dt = co"P„ (л:), |
(5.7.14) |
||||||
0 |
|
|
|
|
|
|
|
где фо(?» x, со) определяется из равенства (6). |
|
|
|||||
Аналогично из разложения (7) имеем: |
|
|
|||||
|
7Z |
|
|
|
|
|
|
v - P - \ r e(t,9, со) ” ” пШ9 = |
\ <оnL n(f), |
(5.7.15) |
|||||
|
о |
|
|
|
|
|
|
00 |
|
|
|
|
|
|
|
j |
ср* (t, 0, ш) Ln (t) dt = |
со" ®o” п0. |
|
(5.7.16) |
|||
О |
|
|
|
|
|
|
|
Наибольший интерес имеют интегральные представле |
|||||||
ния: |
|
|
|
|
|
|
|
Ь [ W r S ) |
ш |
т р ■ ы |
d x - |
S T ! е~ ’,! |
(f)- |
(5 -7Л 7> |
|
JJ 0( | Y |
|
(f)df = |
«0"/2(1+*) Pn(*), |
(5.7.18) |
|||
|
|
- -i-) |
• sin nQdQ = |
Tie~ti2 L n(t), |
|||
|
|
|
* J |
COS |
|
|
|
|
|
|
|
|
|
|
(5.7.19) |
|
|
•e-«2L n(t)dt = 2 sinra0, (5.7.20) |
|||||
|
|
|
|
|
|
cos |
|
которые являются частными случаями (13)—(16) при со = 1. Соотношения (17) и (18) примут более изящный вид, если положить в них a;= cos0:
1Z |
|
J Jo |
tg -f-) P n (COS0) sin -f -d0 = z ir i e~m Ln{t)' (5-7-21) |
133
00 |
|
|
|
J |
|
e - ^ L „ (O ^ = 2 co s-|-^(co s0 ). |
(5.7.22) |
0 |
' |
' |
|
Полагая в (17) £ = 0 и делая замену х —— и в полученном равенстве, имеем
1
^ P n(u)du |
2/2 |
(5.7.23) |
|
2 л +1 » |
что совпадает с формулой (10.10.49) работы [21]. Вышеприведенные интегральные соотношения являются
дуальными, что позволяет их использовать для решения интегральных уравнений 1 рода со специальными ядрами. Например, пара соотношений (21) и (22) дает возможность утверждать, что формальным решением интегрального урав нения
f «Л> ( y tg -§-) / (6)sin y |
dQ = Ш |
(5-7.24) |
||
b |
' |
|
|
|
является функция |
|
|
|
|
|
СО |
|
|
|
А в ) = 2 2- ¥ - 1сл ( |
со80)’ |
<5-7-25> |
||
где сп — лагерровский |
спектр |
известной функции е(/2ф(1), |
||
а также, что формальным решением уравнения |
|
|||
00 |
|
|
|
|
|
|
= |
|
(5.7.26) |
служит функция |
00 |
|
|
|
|
|
|
|
|
ф (f) = |
е-*/» 2 |
ТГ £»(*), |
(5-7-27) |
|
|
< 71=0 |
|
|
|
где ап — коэффициенты |
разложения |
известной |
функции |
|
£ |
|
|
|
|
sec-к- - /(0) по многочленам Лежандра P n(cos0).
А
134
§8. Улучшение сходимости рядов Лагерра
Вэтом разделе изучается вопрос об улучшении сходи мости рядов вида
= |
P(k) |
, |
|
|
(s.8.d |
|
|
|
|
|
|
||
k=0 |
|
|
|
|
|
|
со |
|
|
|
|
|
|
в» = 2 |
а ( 4 |
)*% *(*), |
(5.8.2) |
|||
ъ= о |
v |
' |
|
|
|
|
где L%[t) есть многочлены Лагерра |
с положительными сво- |
|||||
|
2 |
1 k \ |
(—ty |
_ |
и Q(k) явля |
|
|
[ . |
J |
~ц~ |
и Р(к) |
||
ются произвольными взаимно простыми многочленами от к, причем Q(k) не имеет целых положительных корней, а степень многочлена Р(к) ниже степени многочлена Q(k); далее функция а(и) аналитична в окрестности и = 0 и а(0) = 0.
Следуя общей идее А. Н. Крылова [59] об улучшении сходимости рядов, необходимо из рассматриваемых рядов выделить слабо сходящиеся составляющие и просуммиро вать их независимо от оставшихся сумм. Как показал Г. С. Салехов [91, 92], эти слабо сходящиеся ряды удобно суммировать, выделяя из дроби P(k)/Q(k) составляющие ви-
да(Н-1)в^ = 1’
Им же разработан алгоритм этого процесса. Следовательно, функции, улучшающие сходимость рядов (1) и (2), могут быть представлены разложением
00 |
^ |
2 |
( F ХКF I j s |
и наша задача состоит в том, чтобы найти аналитический
вид функции ts (х, t).
Интегрируя повторно s раз по х известное разложение
00 |
x ^ ( t ) = i — |
_ f X |
|
|
2 |
|
|
||
k=0 |
|
|
|
|
получаем |
оо |
|
|
|
|
|
tz |
||
|
1 |
|
||
(Я> t) |
( £1z-^' |
Z ie 1~xdT. |
||
(e-l)f* О |
||||
|
|
|
135
Так как
s —1 |
— t x |
|
(x—‘)s |
ат = 2 ( - i ) v - 1-* |
dx |
||
l — z |
|
|||
|
ft=0 |
|
1 — X |
|
|
|
|
|
|
и |
X |
|
|
|
|
—tx |
|
|
|
|
.f |
T^e1-T |
r ( ft)^r+ i(x, t), |
(5.8.3) |
|
1—x dx — 2 ( - D |
|||
|
|
r —O |
\ r 1 |
|
где A T+i(x, |
t) |
вычисляются из |
рекуррентных соотношений |
|
А т+1(х, t) = |
— |
1 — (1 — x f е |
— t A r (г = |
1 ,2 , .. .), |
причем А г(х, t) = е* Ei Г |
- E i ( - t ) |
E i ( - t ) |
||
ТО |
|
|
|
|
|
|
|
|
|
|
|
|
s —1 |
|
|
|
|
Лх,Я |
* |
+ (,_ !)! 2 |
0(xJ i |
( а 1) ^ |
1)Г(^ ) Ar+l(*’ *)* |
|||
|
|
|
|
|
|
|
|
(5.8.4) |
В частности, при х = 1 отсюда находим |
|
|||||||
|
|
Cs(l, t) — — —- + |
* ^ A s (1, t). |
(5.8.5) |
||||
Здесь величины A s (1, |
t) |
вычисляются по рекуррентным со |
||||||
отношениям |
A s+1(1, t) = ~ [1 — t A s (1, f)], |
А х(1, t) = |
||||||
= — e*Ei (—t) |
(s = |
1,2,. .. ) . |
|
вычислениями инте |
||||
Приведем |
расчеты, |
связанные с |
||||||
грала (3): |
|
|
|
|
|
|
|
|
00 |
— t x |
|
|
|
|
1/1—гс |
|
|
X* |
|
l |
|
|
|
|
|
|
e ^ d x |
= |
и |
|
|
|
|||
I1 ^ |
|
|
1—X |
|
|
|
|
|
136
1/1—X
- 2 |
J |
Далее, интегрируя но частям, имеем:
1 / 1 - х
|
|
1/ 1— X |
|
|
1 |
1 / 1 - х |
= |
]. |
-f I |
Таким образом построены функции (4) и (5), улучшаю щие сходимость рядов (1) и (2).
Для иллюстрации воспользуемся примером, аналогич ным приведенному в работе [91] :
nX* |
(t) |
<г(1, *) + |
С2(1, *)+Са(!,*)“ |
|
л2+ 1 |
||||
71=1 |
|
|
|
|
- 1 0 |
2 |
^ |
(О |
|
(л + 1) (л +2) (л+ З) (л3+ 1) |
||||
|
71=1 |
|||
Улучшающие функции £г (1, t), (г= 1, 2, 3) вычисляются по формуле (5), используя их, окончательно имеем:
2 |
nL* (t ) |
-з- + t + (f — 2fefEi (— f) |
|
л*3-}-! |
|||
71= 1 |
|
_ V _______10-Mt)_______
^ (л + 1 ) ( л + 2) (л +3) (л2+ 1 ) ’
71 = 1
где полученный ряд имеет лучшую сходимость, чем первоначальный ряд.
137
* * *
Задача о построении решения интегрального уравнения
00 |
|
F(p) = С e-Pt f(t) dt |
(1) |
о
ставилась задолго до того, как возникла теория операционного исчисле ния. К ранним работам этого направления следует отнести работу Мор ф и 5, в которой обосновывается возможность построения решения уравнения (1) в виде ряда по многочленам Лагерра.
Последующие исследования уравнения (1) (до сороковых годов) сти мулировались главным образом развитием теории ортогональных рядов, постановкой и решением задач проблемы моментов, изучением свойств аналитических функций класса Нг и их обобщением. Теория операцион ного исчисления, возникнув как часть прикладных наук, в пятидесятых годах приобрела формы самостоятельной математической дисциплины, имеющей широкий диапазон приложений.
Расцвет операционного исчисления совпал с бурным развитием вы числительной техники и ее внедрением в научно-технические расчеты.
В связи с этим вновь возродился интерес к решению уравнения (1), но уже в рамках разработки вычислительных (приближенных) методов.
Исследования, посвященные проблеме численного обращения, можно условно разбить на четыре взаимосвязанных направления.
Первое направление характеризуется конструированием различных б-образных последовательностей, приспособленных к решению проблемы восстановления оригинала по его изображению. Сюда следует отнести формулу обращения Поста—Уиддера [190, 223], различные обобщения которой были даны Поллардом [187, 188]; формулу обращения Вине ра—Шли [28], а также формулу Боас—Уиддера [125]. Кроме того, ука жем на работу [163], в которой формула обращения содержит функции Бесселя. Все эти формулы имеют скорее теоретический интерес и мало пригодны для практических расчетов.
Идея использования б-образных последовательностей для прибли женного вычисления значений оригинала использована И. И. Рябцевым [90], который сконструировал формулу обращения на базе значений изображения, исчисленных в равноотстоящих точках действительной оси. Это обстоятельство делает формулу И. И. Рябцева практически бо лее удобной г. сравнении с формулой Поста—Уиддера, в которой приме нено значение производных высшего порядка изображения, исчисленных в различных точках действительной оси. Практические приложения формул подобного типа рассмотрены в работах [156, 199]. В частности, Костом [138] они использованы для решения задач вязкоупругих сред.
Второе направление основано на сведении интегрального уравнения (1 ) к системе линейных алгебраических уравнений посредством различ
ных квадратурных формул, а также на применении метода регуля ризации решения некорректных задач, предложенного А. Н. Тихоновым
[102, 103].
Последование этого плана принадлежит Нордену [179], который за дачу обращения связывает с решением некоторой системы алгебраиче ских уравнений, при этом значения изображения исчисляются им в последовательности действительных точек, образующих геометрическую прогрессию.
5 Цитируется по [118].
138
Метод регуляризации использован в работах [104, 54]. Этот подход к проблеме обращения представляется наиболее перспективным.
Большое и интересное исследование проведено Р. Веллманом с сот рудниками [119, 120]. Ими использованы результаты Филлипса [184], Тумея [217], а также Шанкса [206], вскрыты трудности, обусловленные некорректностью задачи численного обращения, указаны ограничения на оригиналы, обеспечивающие сходимость к точному решению приближен ного решения, полученного из системы алгебраических уравнений, приб лиженно описывающих уравнение (1 ).
Однако замена точного равенства (1) приближенной системой неточ ных равенств представляется нам наиболее уязвимым местом этого приема в сравнении со всеми методами остальных направлений решения проблемы численного обращения.
Третье направление связано с применением различных типов квад ратурных формул к интегралу Меллина—Римана, являющегося, как из вестно, основой аналитического аппарата обращения преобразования Лапласа.
Первыми в этом направлении были работы Солзера [194—198]. Так, автор [197] предлагает применить к интегралу Меллина — Римана квадратурную формулу наивысшей алгебраической точности в одном частном варианте весовой функции. В последующем эта идея изучалась в широком классе весовых функций [207, 117, 5, 8].
Достаточно полное и систематическое применение квадратурных формул к задаче приближенного вычисления интеграла Меллина—Рима на осуществлено В. И. Крыловым и Н. С. Скоблей [61—64, 95—97]. Ими разработаны таблицы [64, 96], пригодные для практических расчетов.
К. Ланцош [70] ставил вопрос о применении к интегралу Меллина— Гимана формулы трапеций.
В. М. Амербаевым [8] показано, что оценка точности формулы тра
пеции тесно связана с формулой суммирования Пуассона для односто ронних интегралов Лапласа. Автором также к приближенному вычислению интеграла Меллина—Римана применены методы приближенного вычисления интегралов с быстроосциллируюгцими функциями [109], в частности, построены квадратурные формулы с чебышевскими узлами интерполирования и указаны приемы счета по этим формулам. Следует отметить, что формула суммирования Пуассона для одностороннего ин теграла Лапласа приводит к функциональному соотношению, которое было получено ранее [83, 84] в частном варианте и рекомендовано за тем для задачи численного обращения. Формула суммирования Пуассо на служит основой связи между г-преобразованием и преобразованием Лапласа. Отметим, что Вихом [219—221] к проблеме обращения приме нено г-преобразование. Им рассмотрен случай, когда изображение являет ся дробно-рациональной функцией.
Четвертое направление исторически возникло раньше всех осталь ных и связано с разработкой приемов разложения изображения на сум му (конечную или бесконечную) изображений, оригиналы которых вос станавливаются точно.
Сюда прежде всего следует отнести так называемые первую и вто рую теоремы разложения Хевисайда [161]. Обобщение второй теоремы разложения на класс мероморфных функций осуществлено Черчилем
[133, 134].
Практически ценным и достаточно универсальным аппаратом проб лемы численного обращения является теория ортогональных рядов Лагерра. Это обусловлено тем, что изображения многочленов Лагерра кон формным отображением расширенной комплексной плоскости на себя сводятся к степенной функции. Этот факт лежит в основе многих иссле дований рядов Лагерра. В частности, он использован в данной книге для
139