Файл: Амербаев, В. М. Операционное исчисление и обобщенные ряды Лагерра.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.10.2024
Просмотров: 83
Скачиваний: 0
а также В. М. Амербаева [1], которым рассмотрены алгоритмы разложе ния оригиналов в ряды по многочленам, ортогональным на отрезке [0, 1 ]
весом ы(х) (в частности, случаи Якоби, Лежандра, Чебышева). К этому циклу следует отнести работы Армстронга [116] и Мендела [177].
В*. В. Солодовниковым, А. Н. Дмитриевым, Н. Д. Егуповым [99] эти методы развиты применительно к задачам анализа и синтеза систем ав томатического управления. Ими показана высокая практическая цен ность методов разложения оригиналов в классические ортогональные ряды.
Практическое значение численного восстановления оригинала пос редством ортогональных рядов во многом определяется скоростью сходи
мости рядов. В этой связи представляет |
интерес работа [62], в ко |
торой формулируются некоторые критерии |
сходимости ортогональных |
рядов в терминах функции-изображения. |
|
Отметим, что метод разложения оригинала посредством ортогональ ных рядов в случае двумерного преобразования Лапласа был обсужден Бергером [122] и Р. Т. Джаембаевым [37].
Из ранних работ, посвященных проблеме обращения, необходимо также упомянуть исследование Бейтмена [118]. В отличие от предыду щих в ней аппроксимации ортогональными рядами подвергается не ори гинал, а само изображение. Это приводит к новым типам разложения оригиналов по семействам функций, отличных от классических ортого нальных функций и вместе с тем сохраняющих некоторые свойства ортогональных функций. Этот подход использован в данной книге при раз работке БГЛ-алгоритма, а также в других [14, 17, 19]. Кроме того, сле дует указать на работы П. И. Кузнецова [66, 67], широко использующе
го функции Ломмеля двух переменных в качестве аппарата восстановле ния оригиналов, соответствующих специальному классу изображений.
В. М. Амербаев [4, 7] рассматривал общие принципы разложения оригиналов в ряды Неймана. Им, в частности, получены новые функци ональные соотношения, представляющие независимый интерес для тео рии рядов Неймана.
Несколько ранее ряды Неймана по функциям Бесселя с целыми ин
дексами использовал Кэмби [130] в качестве аппарата численного обра щения преобразования Лапласа.
Г л а в а 6
ОБ ОДНОМ КЛАССЕ ЯДЕР УРАВНЕНИЯ ВОССТАНОВЛЕНИЯ
Рассмотрим уравнение восстановления t
U(t) + X § K ( t — х)U(x)dx = f(t).
О
Представляет интерес изучение следующего вопроса: ука зать приемы построения таких ядер K(t—т), для которых резольвента R(t—г) выражается посредством той же функ
ции, что и само ядро.
Примером такого ядра может служить функция Работнова7
00
э.(M )= t-2 Г[.('+1Д ■ <о<«<1). |
а) |
|
71=0 |
|
|
Для такого ядра справедливо свойство: если |
|
|
г7(0 + ц Jt |
Э« (A, t—x)U(x)dx = f(t), |
|
0 |
|
|
t |
|
|
то U(t) = f(t)—ц JiJa (К, |
t —x)f(x)dx, где Ra( l , t —т) = Эа(^+ [х» |
|
О |
|
|
t —т).
Следует отметить, что на поведение ядра K(t—т) могут налагаться некоторые дополнительные ограничения. Так, в наследственной теории упругости ядра K(t—т) достаточно хорошо описываются функциями вида
7 Если функцию (2, 3. 2) обозначим через Эа (X, f), то функция (1)
связана сЭ а (X, t) соотношением Эа (X, 4)= Э1_ а(Х, t).
142
K(t—т)= A(t—x) aexp[—(it—x)i_K] — ядро Вронского,
K(t—x)=A(t—x)~“exp[— —x)] — ядро Ржаницына.
Поведение функции Работнова качественно близко к по ведению ядер указанного типа. Вместе с тем, как отмеча лось выше, достоинством этой функции является то, что резольвента ядра, представленного функцией Работнова, также выражается через функцию Работнова. Это свойство явилось причиной большой популярности функции Работ нова в прикладных задачах наследственной теории упру гости.
В настоящей главе методами операционного исчисления (главным образом методами разложения) изучаются неко торые новые свойства функции Работнова, а также приво дится общий принцип конструирования функций, обладаю щих важнейшими функциональными свойствами функции Работнова. Рассмотрены частные реализации этого принци па. Уравнение восстановления ниже интерпретируется в основном в терминах наследственной теории упругости. Так, ядро этого уравнения называется ядром последействия, а резольвента этого ядра — ядром релаксации.
§1. О представлении функции Работнова
взамкнутой форме
Изучению свойств функции Работнова и построению различных аппроксимаций для нее посвящены работы [10, 18, 24, 85, 86, 88]; исследование ее асимптотических
t
свойств, а также функции Эя (Я, x)dx проведено в работе
0
[18].
Ниже изучается вопрос о представлении рассматривае мой функции в замкнутой форме, удобной для решения за дач табулирования. Под термином «замкнутая форма» по нимается такое представление функции Работнова, которое обеспечивает возможность применения более эффективных способов вычисления значений функции, чем простое сум мирование медленно сходящегося ряда (1.1). Очевидно, что для упрощения задачи табулирования функции Работнова требуется выделить имеющуюся особенность в точке t = 0, иначе говоря, представить функцию Эа (Я., t ) в виде суммы Э« (Я, t) = фа (Я, г) + Фа(Я, it), где Фа (Я, t ) — замкнутая форма суммы слабо сходящихся компонент ряда (2.3.2), причем
143
фа (X, 0) = 0, а другое слагаемое qv(A, t) — конечная часть выделенной особенности функции Работнова в точке (= 0 и имеет форму, удобную для вычислений.
Отметим некоторые ограничения, налагаемые на пара метры а и А функции Эа (А,, £), представленной рядом
|
\п +ап |
(6.1.1) |
|
0 < а < 1. |
|
|
71=0 |
|
1. |
При пользовании различного рода таблицами функ |
|
ций, содержащих много параметров, известное неудобство представляет наличие большого числа входов в таблицу. Ясно, что выписывание параметра А независимо от основ ного аргумента t в ряде (1) диктуется как соображениями компактности записи, так и некоторыми функциональными
свойствами функцииЭа (Я, £). |
|
Я = | Я | signA |
|
Однако, так как при любом вещественном |
|||
справедливо равенство |
|
|
|
Эа (A, t) = I А| |
l—i |
_i_ |
(6.1.2) |
* Эа (signA, | А| а £), |
|||
то анализу достаточно |
подвергнуть |
лишь |
функции |
Эа ( + 1, t) |
и Эа (—1, £), чем мы и воспользуемся в этой главе. |
|
2. Будем считать, что индекс а принимает рациональные |
||
значения |
( а = — , где 0 < — < 1 и т и п — вообще говоря, |
|
|
71 |
71 |
взаимно простые числа).
Остановимся вначале на интегральном функции Работнова.
Случай Я = + 1 .
Имеем Э* (+ 1, t) •*- m/ |
. |
ПР
Так как
— |
м 71““Х/ ^ »^ ^ |
^ »_ |
|
И |
* = 0 4-7' “ г2= 0 т |
||
|
|
п— 1 |
, ч , mk |
p m (N+ Un) |
г г ~, 20 |
\№Р I . |
|
то |
|
|
|
|
п — 1 |
N |
|
Э т (1, £)= 22
*= 0 г - 0
представлении
+
J
144
л—1 |
|
|
|
|
|
- t ) |
|
|
|
|
|
э т(1, х) dx. |
||
+2 Г [mN + m(k + 1)/л] j V |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
m l f + m ( k ± l ) _ 1 |
|
|||||
А=О |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В частности, при N = О имеем |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
п—1 |
m (k + l) |
|
|
|
|
|||
|
|
Эт (1, t) — ~ |
'S |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(6.1.3) |
, |
Л — 1 |
|
1 |
|
(* |
|
У*П!Ч |
|
|
|
|
|||
V |
|
|
|
|
|
|
|
|||||||
+ |
2d г I |
т (А + 1)1 |
] V — *) |
в |
|
Bm{l,x)dx. |
||||||||
|
*“ ° |
[ |
л |
Jo |
|
|
|
|
|
|
|
|
|
|
Случай |
л = —1. |
|
|
|
|
|
|
|
|
|
|
|
||
Имеем Эт (—1, f)-*- |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
Р т >п + 1 |
|
|
|
|
|
|
|
|
|
|
Так как при «-четном |
|
|
|
|
|
|
|
|
|
|
||||
|
1 |
|
»—* |
, „ |
, т(*+1) |
ЛГ , _ |
|
|
||||||
|
4 - 2 |
H |
) * |
7 |
r 2 |
|
7 |
' |
+ |
|||||
|
|
|
|
*=0 |
|
|
|
|
|
г—0 |
|
|
||
|
I 1 |
\тя |
|
1 |
^ |
|
|
|
/1 |
-,™<*+1> |
||||
|
+ |
т |
|
|
|
|
|
|
|
Т |
/ |
|
" |
, |
то |
|
|
|
|
|
|
|
|
|
\ х |
|
|
|
|
|
|
|
|
|
хг |
|
ш(ft+1) . |
|
|
|
||||
|
|
|
Л — 1 |
|
|
|
|
|
||||||
|
|
|
|
|
|
К |
|
— —— + т г —1 |
|
|||||
|
|
|
|
|
|
|
|
t |
» |
|
|
|
|
|
Э т ( |
1» 0 = 2 ( ~ 1 ) 4 2 И-(¥+-)] + |
|||||||||||||
У |
|
( - D * |
|
г |
|
m (ft+ l) |
|
|
|
|
|
|||
|
|
|
v--------- - +7rciV-*l |
9m(l,^)dx |
||||||||||
+ £ г [ „ |
( i t ! + „ ) ] ] < * |
|
") |
* |
|
|
|
|||||||
и в частности, при 2V= О |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
п—1 |
|
|
t |
m(k+1 ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
п |
|
|
|
|
|
Э » ( - 1, й - |
4 |
2 ( - » * |
-г^< *+ я'1 |
+ |
(6.1.4) |
|||||||||
|
|
|
|
|
& |
|
г (— s - ^ J |
|
|
|
||||
1 0 - 5
145
”v |
- |
|
( - 1 ) " |
C |
mUt+1) —^ |
|
|
|
( - Ц * |
j (f |
*0 * |
3m(l,t)dT |
|
+ 2 j V [m(t+l) j |
||||||
2d v |
[ |
n |
|
|
|
|
ft=0 |
|
|
|
|||
Аналогично при и-нечетном получаем
э „ ( - 1 , t ) = 2 ( - D * 2 |
t |
n |
+ В Г-1 |
+ |
||||||
В | Ш +Г |
||||||||||
|
|
|
* = 0 |
|
r = 0 |
|
||||
|
___(-1У\i+iV—1 |
|
|
|
|
|
|
|||
+ 2 |
^ * ± 1+ ^ ] I |
|
71 +* |
|
|
|||||
*=o |
Г [ |
\ n |
1 1 |
о |
|
|
|
|
|
|
и в частности при N = О |
|
|
|
|
|
|
||||
|
|
|
|
|
n—3. |
„ |
"L * * ” |
|
||
|
|
Э r - 1 |
rt== J _ y |
<-»** |
" ■+ |
|
||||
|
|
3f ( h t ) |
* j " r p ^ ± i > ] + |
|
||||||
|
|
|
|
|
|
|
|
|
(6.1.5) |
|
|
^ |
(-1)*-1 |
* |
”»(*+!) |
|
|
|
|||
+ |
2 d |
V r » ( * + l ) l |
] (* — |
x) |
" |
Э щ ( — 1» T)d x - |
||||
|
* = “ |
[ n |
J o |
|
|
|
|
|
|
|
Таким образом, как видно из формул (3)—(5), функция |
||||||||||
Работнова рационального индекса а = |
|
выражается через |
||||||||
ту же функцию целого индекса т.
В связи с этим представляет интерес рассмотреть свой ства функции Работнова целого индекса.
Покажем, что функция Эт (А, t) целого индекса выража
ется через элементарные функции. |
урав |
|||
Действительно, пусть га — первообразный корень |
||||
нения 2 т-—1 = 0. |
|
|
|
|
Тогда, суммируя по р ( О ^ т ^ р —1) следующие ряды |
|
|||
|
00 |
(xt)k = topxt exp (txmP), |
|
|
|
У. ~ |
|
||
получим |
к - 1 l w |
|
|
|
|
m —1 |
т —1 |
|
|
ce |
|
(6.1.6) |
||
2 |
"гШ |
2 |
фТк— xt 2 “'’expi^xt). |
|
к- 1 |
1 ' |
р - 0 |
г - 0 |
|
146