книги из ГПНТБ / Видаль П. Нелинейные импульсные системы
.pdfляемую натуральными числами, можно считать f(n) не прерывной функцией п. В таком случае имеем уравнение первой степени с конечными разностями. Для соотноше ния вида
f(n + p)=F^f(n), f ( n + 1), |
.... f(n + p— 1), |
я], |
(1-1) |
||
или |
|
. . .+ fp f(n )= g |
(1 - 2 ) |
||
f(n + p)+hf(n + p—l ) + |
|||||
знание р значений последовательности f(ri) |
позволяет |
||||
определить следующее значение. |
|
|
|
|
|
Эти соотношения являются рекуррентными последо |
|||||
вательностями р-го порядка. |
Это |
же f(n) |
можно |
рас |
|
сматривать как непрерывную |
функцию от |
п. |
Это |
соот |
ношение является разностным уравнением р-го порядка. Решить рекуррентное соотношение — это значит опре
делить функцию |f(n) |
по р |
начальным |
значениям /о, |
||
fi> •••>fp-i- |
|
|
|
|
|
|
|
б) |
|
|
РАЗНОСТИ В |
По определению первая разность, обозначаемая Д, |
|||||
определяется так: |
|
|
|
|
|
|
A f(n )= /( n + l) - f ( n ) . |
|
(1-3) |
||
Последующее применение этого определения приво |
|||||
дит к разностям высшего порядка |
|
|
|||
A2 {f (л) ] = А [А/ (п) ] = A {f (п+ 1) —/ (п) ]= |
|
||||
|
= f (н + 2 )—2f(n+1) 4-f(n); |
|
(1-4) |
||
Д’ [f (п)\= |
Д [Д2f (п)} = |
Д [f(n + |
2) - 2/ (я + |
1)+ "I |
(1-5) |
+ f(n)] = |
f(n+3)- Zf(n4 -2) + |
Sf(n+ \) - f(nj l |
|||
|
Amlf (п)} = Д”» '1 / (п + 1) - Д”»-1/ (я). |
(1-6) |
|||
Примем обозначения |
|
|
|
||
|
/( л + 1 )= :(1 + Д )/( я ); |
|
(1-7) |
||
f ^ + 2)= ( l + A ) f ( / i + l ) = (l+ A )*f(*); |
(1-8) |
/ ( « + я » ) ; = ( 1 - к д ) да/ ( « ) - |
(1-9) |
Пример 1-1. Для линейного разностного уравнения второго по рядка с постоянными коэффициентами
aof(n) + a if(n + l)+ a 2f(n+2) = g (« ), |
(ЫО) |
где f(n) — неизвестная функция, п — независимая переменная, |
коэф |
фициенты do, fli, а2 и функция g заданы. |
|
Это уравнение связывает значения f, отстоящие соответственно на один или два интервала. Без потери общности путем простой за мены переменных можно перейти к последовательностям первого порядка.
Уравнение (1-10), полностью соответствующее вышеназванным определениям, может быть представлено в другой форме, более ясно иллюстрирующей аналогию с дифференциальным уравнением.
Подстановка уравнений (1-7) и (1-8) в (1-10) |
позволяет на |
писать: |
|
W V (n ) +biAf(n) +bof(n) =g(n). |
(1-11) |
Коэффициенты b ввязаны с коэффициентами а следующими со |
|
отношениями: |
|
Ьо — й0-f- я, -f- Яг! |
|
6, = я, + 2аг\ |
( 1- 12) |
Ьг= яг.
Аналогия между уравнением (1-11) и линейным дифференциаль ным уравнением второго порядка с постоянными коэффициентами ясна. Др соответствует dP/dnP. Различие заключается в конечной величине независимой переменной п в разностных уравнениях, стре мящейся к нулю в дифференциальном уравнении.
в) ВЕКТОРНО-МАТРИЧНЫЕ РАЗНОСТНЫЕ УРАВНЕНИЯ
Рассмотрим вектор F(n), определенный через
f{n)
f ( n + p - |
1) |
и обозначим через А матрицу ( п Х р ) |
|
0 |
0 |
01
- и
14 «
й через В(п) вектор-столбец
О
В (л )= ;
о
е
Если для данного п функция fp удовлетворяет соот ношению (1-2), то F (п) удовлетворяет рекуррентному соотношению
F (/i+ l)= A (n )F (n ) + B(n) |
(1-13) |
и соответственно уравнение (1-13) называется неодно родным уравнением,а уравнение
F(n+l) =A{n) F( n) |
(1-14) |
называется однородным уравнением, соответствующим уравнению (1-13).
1-2. ПРЕОБРАЗОВАНИЕ КАРСОНА
Дискретное преобразование Карсона F*(p) |
функции |
J f(n) определяется как |
|
СО |
|
с [f(x)] = F*(p ) = p f lf{x)\e-v*dx, |
(1-15) |
9 |
|
где p= a+ jсо — комплексная величина (параметр преоб разования).
Выражение (1-15) можно сравнить с классическим преобразованием Карсона
ПР) = Р \ e-*4{t)dt. |
(1-16) |
о |
|
Сходство между этими двумя выражениями очевид но [непрерывная переменная t, соответствует дискретной переменной х, и непрерывная функция f (t) — ступенча той функции f(x)].
Интеграл, определяемый выражением (1-15), сходит ся к конечному пределу, когда р>оо, где сг0 — радиус сходимости.
15
Приведем |
примеры |
использования |
преобразования Карсона. |
||||||
1) fW = 1(*), |
|
|
|
|
|
|
|
|
|
где |
|
|
|
( |
1 при х |
0; |
|
|
|
|
1 |
(х) = |
|
|
|||||
|
< |
0 при |
х < |
(рис. 1-3) |
|||||
|
|
|
|
( |
0. |
|
|
||
Используя (1-16), |
получаем: |
|
|
|
|
||||
|
|
|
|
00 |
|
|
|
|
|
|
F* (р) = |
р J 1 e~v*dx — — [<?“ рх]§° = 1• |
|||||||
|
е**, |
|
|
о |
|
|
|
|
|
2) } (х) = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1f ( X |
х ' |
|
|
|
|
|
|
|
|
|
|
|
7,1 |
|
|
|
|
|
|
ч Г Г " |
||
|
|
|
|
|
|
|
J |
b |
<х<0 |
|
|
|
|
|
X |
|
1Д С |
||
|
|
|
|
|
|
|
0 |
1 г з ь s |
|
|
Рис. 1-3. |
|
|
|
|
|
Рис. 1-4. |
||
где а — действительное |
число (рис. |
1-4). |
|
|
|||||
|
+00 |
|
|
|
|
00 |
п+\ |
||
F * ( p ) = p | |
[ J |
вах \ e - v * d x = ^ p |
j |
еапe~Pxd x = |
|||||
оо |
О |
я+ 1 |
|
-| |
оо |
/1=0 |
u |
|
|
р |
|
|
|
|
|||||
= 2 i |
р |
j ' e~p*dx |= |
|
[ c - nP — <?-("+»)?] == |
|||||
|
я |
|
|
-I |
n=.0 |
|
|
|
|
ti—О |
|
|
|
|
|
|
|||
|
|
00 |
|
|
|
1—<?-p |
|
e? —■1 |
|
(1 — <?-J>)^' a — p ) n == |
|
||||||||
|
(P |
^ |
ev — ea |
||||||
|
|
n=0 |
|
1 — e |
|||||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
а) |
ОСНОВНЫЕ ТЕОРЕМЫ |
1. Теорема суперпозиции
J Дискретное преобразование Карсона обладает свой ством суперпозиции, т. е. преобразование произведения функции на константу можно записать:
C{af(x)] = aC[f(x)]=aF*(p), |
(1-17) |
а если |
(1-18) |
f(x)— ag(x)+bh(x), |
то получим:
C{f(x)]=F*(p)=aG*(p)+bH*(p),
16
где
б* (Р) = С [£(.V)]; И* (Р) =С{1г (х)].
2. Теорема смещения
Дискретное преобразование Карсона функции J f(x) обозначается так:
С [(/ (х)] = F* (р) = р J J / (х) e~Pxdx,
О
откуда следует, что
Clf(x + k)]= p Jj(f(x + k)e-p*dx =
О
00
= р ( J f (г)е~р{г~к) dr.
г—к
Считая r = x-\-k, получаем:
+ 00 |
n—k—1 |
C[f{x + k)) = p0* j |
J /О') e~Prdr—pevh J J f(r)e~Prdr. |
Интегрирование от 0 до бесконечности не зависит от обозначения, принятого для переменной. Учитывая это, получаем:
С ( /( * + k)] = er>* F*(p) — 2 |
J/OO^P j |
e-p’-dr'j |
||
|
k-\ |
|
|
|
■.ePhF *(p )-eP h'£l (l ~e~p)e-Pr l |
f (r); |
|
||
C [/ {x + k)\= eP*F* (p) - |
1))£ e " |
J / (r). (1-19) |
||
|
|
о |
|
|
Исходя из этого соотношения следует, что |
|
|||
C[f(x+l )]=ePF* (р) - ( е р - 1)/ (0); |
(1-20) |
|||
C{f(x+2)]=e*PF*(p) — (e2P—eP)f(Q) — (еР— 1)/(1); |
||||
|
|
|
|
( 1-21) |
С (/(х + 3)] |
еЗРр* (р) — (е3р _е2р)f (0) _ |
|
||
— (е2Р— |
е Р )/(1 > — |
( е Р 1)/( |
И&учно - |
|
2—352 |
библиотек* СС СР |
|
сКссГУоП/иц*3
3.Изображение приращений (разност
C W ( x ) ] = C [ f ( x + l ) - f ( x ) ] = C[f(x+ \)] -
—С [/(х )]= (е?— 1) {F*(p)—/(0)]; |
|
СШ( х)]= ( е р - 1) {F*(p)- f( 0)]- |
(1-23) |
В частности, |
|
F * (p )= t^ r= T -+ f(°)- |
О '24) |
Аналогичным образом докажем, что изображение второй разности будет:
С[Д2Д*)]=(еР— 1 )2F*{p)-f(0) (е*— 1]*—Д/(0) (еР— 1) (1-25)
A /( 0 )= f( l) - f (0 ) .
Таким же образом
С,[Д2/(х) ]= ( е Р - 1)3F* (р) —f (0) [еР— 1Р -
—А/(0) (еР— I)2—А2/(0) (еР— 1), |
(1-26) |
||
А2/(0) = /(2 ) —2f (1) + /(0 ). |
|
||
|
|
4. Теорема свертки |
|
Можно показать, что если F*i(p) |
и F*z(p) |
дискрет |
|
ные преобразования Карсона от fi(x) |
и /г(я)> т0 |
|
|
F*(p) = F\ (p)F\(p) = |
С S fi (m)f2 (x — m)\, |
||
|
m=.0 |
|
|
ИЛИ |
|
|
|
F*(p) = C S |
f, (■« — да) f. (да) |
(1-27) |
|
m=0 |
|
|
Эта теорема очень полезна при отыскании оригинала, соответствующего произведению изображений.
5.Теорема о конечных и начальных значениях
Применение выражения (1-24) позволяет показать,
)
Пт / (*) = |
lim e?F* (р) |
(1-28) |
х~>оо |
Р+о |
|
18
и что
|
lim f(x) = f (0) = |
lim |
CV |
|
(1-29) |
|
|
ер — 1 F(p)- |
|||||
|
|
|
р-*00 |
|
' Г |
|
|
|
|
|
|
|
|
|
|
б. Дискретное преобразование Лапласа |
||||
Можно определить дискретное преобразование Лап |
||||||
ласа ![Л. 1-6] по формуле |
|
|
|
|
||
|
|
00 |
• • |
|
|
|
|
|
F\ (Р) = V e - p * f (х) = L [f (х)}. |
(1-30) |
|||
|
|
jc=0 |
|
|
|
|
Дискретное преобразование |
Карсона |
имеет вид: |
|
|||
|
|
|
+ 00 |
|
|
|
|
F*(p) = F\(p) = |
p f |
Jf(x)e~r>*dx = |
|
||
|
оо |
лс+1 |
о |
оо |
х + \ |
|
|
|
|
||||
= |
S p |
[ I f {x)e-v*dx=Yi If(x)p |
^e--Pxdx = |
|||
|
x=zO |
x |
*= 0 |
x |
|
|
00 |
|
■X |
|
00 |
|
|
|
4 |
|
|
|
||
= — S |
I /Чх)[г~р(*+1) — e-T>x\— — S i |
f (x)e~r>x(er>~ 1). |
||||
*= 0 |
|
|
|
x = 0 |
|
|
Отсюда соотношение, связывающее эти два преоб |
||||||
разования, будет: |
|
|
|
|
||
|
F*(p) = - ( e P - l ) f L(p) = - ? £ ± F * L(p). |
(1-31) |
С помощью этого равенства производится модифика ция общей таблицы дискретных преобразований Карсо на. Дискретное преобразование Лапласа получается умножением дискретных преобразований Карсона на
е р(е Р— 1 ) - \
Дискретное преобразование Карсона имеет то пре имущество, что изображение константы есть константа. Это соотношение важно при изучении разностных урав нений системы (переходная функция системы есть: пре образование Карсона реакции на единичный входной импульс).
Доказательства, касающиеся свойств этого преобра зования, очевидно, одинаковы. Отметим, однако, что тесь рема смещения записывается как
L {f (•* + 1)] = ev [Fl(р) — 1(0)}, |
(1-32) |
2* |
№ |
и
fl(р) = |
L- ^ r + f (0) |
(1-33) |
||
L If (x + |
2)] = |
g2p [F*L (p ) - |
f (0) - |
f (1) e-p]: (1-34) |
L [f (x + |
3)] = |
e*p [P\ (p) - |
f (0) - |
f (0) - |
|
- f ( l ) e - p '- f ( 2 ) e ~ 4 |
(1-35) |
||
|
б) ИЗОБРАЖЕНИЯ НЕКОТОРЫХ ФУНКЦИИ |
|||
Преобразования |
выполняются с |
использованием |
уравнений (1-3) и (1-23):
Af ( * ) = / ( * + ! ) —f(x);
C [ A f ( x ) l = ( e P - l ) [ F * ( p ) - f m
Рассмотрим несколько примеров. Используем при этом табл. 1-1, в которой приведены преобразования
обычных функций. |
|
1. |
И |
Найдем С (а), где а — действительная величина,
f(x)=a и A f{x )= a —а = 0.
Применив уравнение (1-23), получим:
(еР— 1)С (а)—а(еР— 1) = 0 ,
С (а)=а. (1-36)
2.Изображение равномерно нарастающей функции
Пусть f(x) — х (рис. 1-5),
А(х) = х + 1—х = 1.
Согласно уравнению (1-36) (п = 1 )
(ер— l)C (x) = 1,
откуда
0-37)
Таким же образом можно показать, что если f(x) — =х?, то
Д (*2) = (х + 1) 2 _ Х2= 2х + 1;
С(х2)= 2 С ( х )+ С (1 ),
20
и с помощью вышенайденных результатов
О Д |
ер + 1 |
(1-38) |
|
(ер — I)2 ‘ |
|||
|
|
Рис. 1-5. Ступенчатая |
Рис. 1-6. Ступенчатая |
функция f ( x ) —x. |
функция f(x) = С Х. |
Аналогичным способом показывается, что
О Д |
е*р + 4ер + 1 |
|
(1-39) |
|
(ev— I)3 |
’ |
|||
|
|
|||
еьР + |
\\е*р+ \\ер + 1 |
(1-40) |
||
О Д = |
(ер — 1)* |
|
3. Изображение экспоненциальной нарастающей функции
Найдем преобразование \{x)— dx, где d — действи тельная величина. Если d= 0, то мы получаем преоб разование Лапласа (рис. 1-6) от единичной ступенчатой функции
A (dx) = dx+l—dx(d— 1) dx.
Отметим, что разность Ad* имеет тот же характер, что и функция d* (так же, как и для непрерывной экс поненциально-показательной функции):
С [Adx]= (d— l)C{d*];
(er>— 1) С [dx]—eP+1 = (d— 1) С (dx) ,
откуда
С (d*) = |
ер — 1 |
(1-41) |
|
еР— d |
|||
Если d— ea, то |
ер — 1 |
|
|
С (еах) = |
(1-42) |
||
ер — е{ |
|||
|
|
21
Таблица 1-1
Дискретные преобразования Карсона
F*(P)
- у - и ( Р )
ер _ 1 e - av----—
ер
1
1
еР — 1
(^р + 1)
(ер — I)2
(62р + 4й>р + 1)
(е е — I)3
е3Р+ 11е2Р+
(ер — I)4
+11ер+ 1
---------------
1 (ер — I)2
1
( е р - 1 ) 3
1 (еР — 1)п
ер — 1 ер — а
f(х) ДЛЯ X5: 0
и ( х ) — и(х — 1), |
Л |
или и (х) и (1 — х) |
)1 а (х) |
1 (х — а) |
|
1 |
|
X |
|
X2 |
|
Примечания
—
а—действитель ная величина
—
—
X3 —
X4 —
X (х — 1)
2! |
|
х (х — 1) (х — 2) |
|
3! |
|
х (х — 1) ...'(х — п + 1) |
п —положитель- |
п\ |
ное целое число |
ах |
а —действитель |
|
ная величина |
22