книги / Операционное исчисление и обобщенные ряды Лагерра
..pdfа также В. М. Амербаева [1], которым рассмотрены алгоритмы разложе ния оригиналов в ряды по многочленам, ортогональным на отрезке fO, 1J весом со(я) (в частности, случаи Якоби, Лежандра, Чебышева). К этому циклу следует отнести работы Армстронга [116] и Мендела [177].
В. В. Солодовниковым, А. Н. Дмитриевым, Н. Д. Егуповым [99] эти методы развиты применительно к задачам анализа и синтеза систем ав томатического управления. Ими показана высокая практическая цен ность методов -разложения оригиналов в классические ортогональные
ряды.
Практическое значение численного восстановления оригинала пос редством ортогональных рядов во многом определяется скоростью сходи
мости рядов. В этой овязи представляет |
интерес (работа [62], в ко |
торой формулируются некоторые критерии |
сходимости ортогональных |
рядов -в терминах функции-изображения.
Отметим, что метод разложения оригинала посредством ортогональ ных рядов в случае двумерного преобразования Лапласа был обсужден Бергером [122] и P. Т. Джаембаевым [37].
Из ранних работ, посвященных проблеме обращения, необходимо также упомянуть исследование Бейтмена [118]. В отличие от предыду щих в ней аппроксимации ортогональными рядами подвергается не ори гинал, а само изображение. Это приводит к новым типам разложения оригиналов по семействам функций, отличных от классических ортого нальных функций и вместе с тем сохраняющих некоторые свойства ортогональных функций. Этот подход использован в данной книге при раз работке ЫГЛ-алгоритма, а также в других [14, 17, 19]. Кроме того, сле дует указать на работы П. И. Кузнецова [66, 67], широко использующе го функции Ломмеля двух переменных в качестве аппарата восстановле ния оригиналов, соответствующих специальному классу изображений.
В. М. Амербаев [4, 7] рассматривал общие принципы разложения оригиналов в ряды Неймана. Им, в частности, получены новые функци ональные соотношения, представляющие независимый интерес для тео рии рядов Неймана.
Несколько ранее ряды Неймана по функциям Бесселя с целыми ин дексами использовал Кэмби [130] в качестве аппарата численного обра щения преобразования Лапласа.
Г л а в а 6
ОБ ОДНОМ КЛАССЕ ЯДЕР УРАВНЕНИЯ ВОССТАНОВЛЕНИЯ
Рассмотрим уравнение восстановления
■t
U (t)+X §K (t —x)U(x)dx=№ .
О
Представляет интерес изучение следующего вопроса: ука зать приемы построения таких ядер K (t—т), для которых резольвента R (t—т) выражается посредством той же функ ции, что и само ядро.
Примером такого ядра может служить функция Работнова7
00
Э . ( М ) - г - 2 гД + ц -р (•><“ < « • |
< 1 > |
71—0 |
|
Для такого ядра справедливо свойство: если t
17(f)+ p j Э . (Я, f—x)U(x)dx=f(t),
О
t
то U (t)=f(t)— р, J-Rot (Я, f—x)f(x)dx, где R* (Я, f—т) = Э«(Я+рг
О
t—т).
Следует отметить, что на поведение ядра K (t—т) могут налагаться некоторые дополнительные ограничения. Так, в наследственной теории упругости ядра К{$—х) достаточно хорошо описываются функциями вида
7 Коли функцию (2. 3. 2) обозначим через Эа |
(Я, t), то функция (1) |
связана сЭа (Я, t) соотношением Эа (К |
t). |
142 |
|
K (t—T) = A ( t — T ) 1*exp[— Я(?—T)] — ядро Ржаницьша.
Поведение функции Работнова качественно близко к по ведению ядер указанного тина. Вместе с тем, как отмеча лось выше, достоинством этой функции является то, что резольвента ядра, представленного функцией Работнова, также выражается через функцию Работнова. Это свойство явилось причиной большой популярности функции Работ нова в прикладных задачах наследственной теории упру гости.
В настоящей главе методами операционного исчисления (главным образом методами разложения) изучаются неко торые новые свойства функции Работнова, а также приво дится общий принцип конструирования функций, обладаю щих важнейшими функциональными свойствами функции Работнова. Рассмотрены частные реализации этого принци па. Уравнение восстановления ниже интерпретируется в основном в терминах наследственной теории упругости. Так, ядро этого уравнения называется ядром последействия, а резольвента этого ядра — ядром релаксации.
§1. О представлении функции Работнова
взамкнутой форме
Изучению свойств функции Работнова и построению различных аппроксимаций для нее посвящены работы [10, 18, 24, 85, 86, 88] ; исследование ее асимптотических
t
свойств, а также функции |* Эа (Я,, т)dx проведено в работе
Ь
[18].
Ниже изучается вопрос о представлении рассматривае мой функции в замкнутой форме, удобной для решения за дач табулирования. Под термином «замкнутая форма» по нимается такое представление функции Работнова, которое обеспечивает возможность применения более эффективных способов вычисления значений функции, чем простое сум мирование медленно сходящегося ряда (1.1). Очевидно, что для упрощения задачи табулирования функции Работнова требуется выделить имеющуюся особенность в точке f = 0, иначе говоря, представить функцию Эв(Я, t) в виде суммы Э. (Я, f ) — фа (Я, *)+Фа(Я, t), где Ф« (Я, t ) — замкнутая форма суммы слабо сходящихся компонент ряда (2.8.2), причем
•Фа (Я, 0)’= 0, a другое слагаемое фа(Я, *)— конечная часть выделенной особенности функции Работнова в точке f= 0 и имеет форму, удобную для вычислений.
Отметим некоторые ограничения, налагаемые на пара метры а и Я функции Э« (Я, t), представленной рядом
|
со |
|
|
Эа (К t) = Г-1 21 l ^ V + ij T ’ 0 < « < 1 . |
( 6 . 1 . 1 ) |
1. |
При пользовании различного рода таблицами функ |
|
ций, |
содержащих много параметров, известное неудобство |
представляет наличие большого числа входов в таблицу. Ясно, что выписывание параметра Я независимо от основ ного аргумента t в ряде (1 ) диктуется как соображениями •компактности записи, так и некоторыми функциональными свойствами функцииЭа (Я, t).
Однако, так как при любом вещественном Я= |Я|э^пЯ
справедливо равенство |
|
|
|
1—« |
1 |
|
|
Эа (М ) = I * I “ |
Эа (signX, I Я |“ t), |
(6.1.2) |
|
то анализу достаточно |
подвергнуть |
лишь |
функции |
Эа ( + 1 , t) и Эа (—1, f), чем мы и воспользуемся в этой главе. 2. Будем считать, что индекс а принимает рациональные
значения (а = — , где 0 < -^ -< 1 и т и п — вообще говоря,
взаимно простые числа).
Остановимся вначале на интегральном представлении
•функции Работнова. Случай Я = + 1 .
Имеем Э_т (+ 1, t) -*• _т /п_1 •
Так как
|
n—1 |
. m k N , + . |
|
1_______ |
|
|
+ |
l) |
TO |
|
|
|
TJ—1 ДГ |
*т ( т 1+г) - 1 |
эт(1>t) =2 2 |
г [ „ ( Ш + .)] + |
|
Т |
»-о г-0 |
. V _________-_________ |
L |
|
|
|
|
|
|
|
л |
Эт(1 ,z)dx |
||||
+ 2 iT [m N + m(k + l)/n] |
Q |
|
|
|
|
|
|
|
||||||
k=0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В частности, при N = 0 имеем |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
n—i |
|
”*<*+1) |
|
|
|
|||
|
э= |
(î, |
|
~ |
|
2 |
j ï |
|
|
|
+ |
|
||
|
|
|
|
|
|
â r [ |
2 <|±.«| |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(6.1.3) |
^ |
1 |
l / r e |
C |
+ l ) ! |
m(k+l) |
— |
t |
■s) |
я dz. |
3 . ( 1 , |
||||
" b ^ |
p |
C A |
] |
( f |
|
|||||||||
A“ ° r [ — F — Jo |
|
|
|
|
|
|
|
|
|
|
|
|||
Случай À = — 1 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Имеем |
( - 1 , f) - |
|
|
|
|
. |
|
|
|
|
|
|
|
|
Так как при я-четном |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
п —'1 |
|
|
, |
|
* |
/гс(А + 1) |
N |
. ^ |
чг |
|
||
|
|
2 |
н |
)‘ |
4 |
т |
Г |
|
2 |
|
т |
+ |
|
|
|
|
А- 0 |
|
|
* |
|
|
|
|
|
|
|||
|
, _ . |
|
^ |
|
л—1 |
|
|
, |
w , тп(А+1 |
|
||||
н - ( т Н ^ 2 < - » * ( т Г ^ |
|
|||||||||||||
то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
я- |
|
SÎSÜ Î |
+т г - 1 |
|
||||
|
|
|
|
|
|
— |
|
f |
» |
|
|
|
|
э ? ( - 1л , = | ( - 1)‘ | 7 [£^ ± г Гт | +
S ? |
(-D * |
I |
t частности, при N = 0 |
|
|
|
|
ЖЧ |
|
э т ( - 1 , t) = 4 “i |
|
|
в* |
г *=о |
4! » + ! ! +в д ч
’э “ ( , т )
|
|
ЩЯт{к+1)4 |
и |
t |
I* |
п |
||
7 |
|
Ть+тг + |
ПЕШ |
||
|
I |
п |
|
(—1)* |
n |
m(4+l) |
|
|
4 " 2 |
T»Ги*№+1) 1 I (t — т) |
n |
(1, т) dx, |
||
|
I » |
J o |
|
|
|
Аналогично при л-нечетном получаем |
|
||||
|
Л — 1 |
N |
Ш^±Ктт-Х |
||
|
|
|
t |
л |
|
^ ( - i . * ) - 2 < « * 2 г [ „ ( i ± î + r ) ] + |
|||||
п-х |
|
* |
( Ж + * ) _ а |
||
+?0 r[“(Ç}+Jf)] I (t-,) |
|
|
* 4 - U 4 * . |
||
и в частности при N = О |
|
|
|
|
|
|
|
|
|
mlk+1) |
|
Л—1 |
|
|
|
|
(6.1.5) |
(—1)6—1 р |
го(*+1) |
|
|
||
+2„ГлЦй-ЦЛ j (* — т) л |
|
1Эт (—1, |
|||
*=С |
|
|
|
|
|
Таким образом, как видно из формул (3)—(5), функция |
|||||
Работнова рационального индекса |
|
Ш |
выражается через |
||
а = — |
ту же функцию целого индекса тп.
В связи с этим представляет интерес рассмотреть свой ства функции Работнова целого индекса.
Покажем, «что функция Эт(К t) целого индекса выража ется через элементарные функции.
Действительно, пусть о — первообразный корень урав нения г т—1 = 0.
Тогда, суммируя по р {О ^ т ^ р — 1) следующие ряды
00 |
|
|
|
2 щ |
№)к = ®Pxt ехР |
|
|
получим |
^ |
^ |
|
2 ~v!l\ 2 |
— xt 2 |
“Гехр ("fxt). |
(6.1.6) |
4 - 1 W л - 0 |
г - 0 |
Так как |
|
|
|
|
и—1 |
(т, если Ы О (modт), |
|
|
2 «>г*= |
(6.1.7) |
|
|
г=о |
(О, если к ф О (modт), |
|
то равенство (6) преобразуется к виду |
|
||
|
f(r+ l)m —1 |
vl —771 71— 1 |
|
г-о |
Г [т(г+ 1)] |
хтТ= - — 2 <о* exp (®kxt). |
|
|
к=о |
|
Выражение, стоящее в левой части этого равенства, есть функция Эт(хт , t) , поэтому окончательно имеем
|
ж1- 7" |
71—1 |
|
|
|
|
||
9ra(*m,t) = |
2 |
u)fcexp (<Bbjct), |
(6.1.8) |
|||||
т |
|
|||||||
|
|
|
ft-0 |
|
|
|
|
|
|
|
|
|
2та |
|
|
|
|
где х произвольный параметр и а = « в . |
|
|
|
|||||
Следствие 1. При ж= + 1 |
получим для |
произвольного |
||||||
четного или нечетного т |
|
|
|
|
|
|
|
|
|
|
771—1 |
|
|
|
|
||
э т ( + м ) = - ^ - |
2 ш*ехр (*“*)♦ |
|
(6.1.9) |
|||||
|
|
Й = 0 |
|
|
|
|
||
При JC= — 1, если т — нечетно, из |
равенства |
(8) имеем |
||||||
|
|
771—1 |
|
|
|
|
||
э т (—1, t) = |
-jj- |
2 ш*ехР |
ÿco*)< |
(6.1 .10) |
||||
|
|
*-о |
|
|
|
|
||
|
|
|
|
Til |
|
для |
любого ж |
|
Следствие 2. Положим х= Х = ет , тогда |
||||||||
(четного или нечетного) справедлива формула |
|
|||||||
|
|
|
77Z— 1 |
Xü,fcexPUXü)A- |
(6.i.ii) |
|||
эм(-1,t)= --£- 2 |
mftïo
Вцелях дальнейшего изучения свойств функции Работ* нова целого индекса удобно ввести в рассмотрение функцию
1 |
771—1 |
2x1 |
^(») (г) = — |
2 o)W*exp (гш*), |
где со= еЯГ. |
|
Л—О |
|
Функция 1$Чг) обладает следующими свойствами.
Свойство 1. этфт, t) = |
|
|
|
и Р = — 1 (т — нечет |
||||
При р= + 1 имеем Эт (1, t)= li(nt)(t) |
||||||||
ное) Э«(— 1, t)=h<m>(— t). |
|
|
|
|
|
|||
ni |
при любом |
т |
имеем Эда(— 1 , t)~ |
|||||
Бели р=Я=ел», то |
||||||||
——Я/j(m) (Я, f). |
|
|
если N =0(m odm ), |
|||||
Свойство 2. |
|
(0) = I |
||||||
|
если ZV#0(mod/»). |
|||||||
Отсюда сразу следует, что |
||||||||
|
|
|
|
|
||||
|
Эт |
|
1 , если т = 1 , |
|
||||
|
|
0, если т ф 1 . |
|
|||||
|
|
|
|
|||||
Свойство 3. |
|
|
|
|
|
|
|
|
Справедливы следующие равенства |
|
|
||||||
|
|
|
chz, если ZV=0(mod2), |
|||||
|
|
|
sh2, если ZV^Q(mod2), |
|||||
1%>Z = |
|
2izi |
2*i |
|
4K Î |
4JÇÎ |
||
\ - е г + ё * |
+ e 8‘ |
+ |
e 8 + |
e 8 % |
||||
или для функции Работнова целого индекса будем иметь |
||||||||
|
|
|Э1(+ 1 , t) — в1 |
|
|
||||
|
|
(Э1(—1, О = |
в-», |
|
|
|||
|
|
|Э2(+ 1 , |
f) = |
sh f |
|
|
||
|
|
|э2(—1, t) — sinf, |
|
|
||||
Эа(+1» t) — -g-e* |
j f e |
2cos^^Ç-f |
з") |
|||||
Э8( - 1 , |
t) = |
-§-e“f + -§ - e T s |
i |
n f |
-----£-). |
|||
Свойство 4. Z^t)(2) = Zl“ )(2:), |
|
|
|
|
||||
где r=AT(modm), причем i— |
наименьший неотрицательный |
вычет числа W<modm). Это свойство следует непосредствен но из равенства ш*г*=шг*, если r=N (m odm ). Из этого свойст
ва, в частности, |
следует, |
что Эт (1, t) = Z<r>") (t), если г —Is s |
|
= 0(mod/ra). |
|
|
(г). |
Свойство 5. |
(о>4г) = |
<*>~* |
Действительно,, с одной стороны, справедливо равенство
т—1
= - ~ ~ 2 Ч к+в)ехР((в*+*г),
Щ |
к“° |
с другой стороны, при изменении ft от 0 до т— 1 последова тельность чисел S h, определяемых соотношением
S* s ft -J- s (modm), 0 < ft < т — 1 ,
пробегает полную систему наименьших неотрицательных вычетов, поэтому имеем
2 ^ +в)ехР (ü)ft+e 2) = |
2 wm ехР (ш**2) = ^ |
(2)* |
Л«=0 |
*=о |
|
откуда и следует искомое равенство. |
|
|
Свойство 6. Производная |
s-го порядка функции |
(г) |
вычисляется по формуле |
|
|
& {£’<*)}=■iTw.
где r=sn+$(modm).
Это равенство можно доказать следующим образом:
|
|
771“ 1 |
£-s |
(2М |
2 ш(п+в)* ехр <•**)= |
|
Ш—1 |
|
= |
-i- 2 |
шгй ехр (<0kz) = Z(m)(в). |
m*н>
Вчастности, из этого свойства и равенств (7), (9)— (11) следует, что
d»~i |
fl, |
если |
s = 0 (mod/n) |
d t* ~ l |
— |о^ если |
s ф о (modm) |
|
и, наконец, для любого целого т >О |
|
||
d‘~l |
(—l)i+1, если |
s = 0(modm), ft = |
|
|
Z)]f=o — |
если s^O(modnz). |
|
|
0, |
Свойство 7. Повторный интеграл s-ro порядка функции V^ixz) вычисляется по формуле
t
V
m r
ï è ï ) S ( t - x y - 4 ^ ( x 4 . ) d ï = |
W M - 2 Г[го—k + m r ] |
о |
г»0 |
где ftss л— $(mod7ft)
и
если m ~ k > p
v =
[V]* если w — & < р,
где р — наименьший неотрицательный вычет числа s— 1
(modm), а квадратные скобки означают целую часть числа, заключенного в скобки.
Последовательным интегрированием по частям можно
получить |
|
t |
d ' = - à ^ r fexp (*« *) - |
(^ïj![<*“' x)J_1exP |
|
0 |
|
откуда следует, что
|
t |
|
|
0 |
|
l r=0 |
©(n-* )rexp(a>rxf) — |
2 ш(" *+г)/ |
r«*o |
Z*0 |
Пусть к — наименьший неотрицательный вычет числа (л—s)(modm). Тогда сравнение п—a+r^O(modw) будет иметь место лишь для чисел г вида г=тп—« + тР- Границы изменения р определяются из неравенства
■ ü ri+ 1 ) < [!= !l + |
J L - „ - i . |
|
||
т |
1 * |
m J ~ |
т |
|
Из этого неравенства следует, что «ели |
то ^ долж’ |
|||
но изменяться |
в пределах 0^р^= |“пГ | |
** 60,111 же |
т—й^р, то р изменяется в пределахс 0 ^ Р ^ ( \ г }