500
19. Zastosowanie przekształcenia dyskretnego 2t
a następnie dzieląc licznik przez mianownik, otrzymujemy:
+ 13z‘2 + 40z~3 + ...
1 :(1 — 4z"‘+3z~2) = l+4z‘
1—4z~*+3z'2
4z~' —3z~2 4z-‘-16z-2 + 12z-ł
13z‘2-12z‘3 13z~2 —52z~3 + 39z~4
40z~3 —39z~4 40z-3-160z-4+120z'5
wobec tego
F(z) = 1+4z-' + 13z-2 + 40z-3+...
skąd
/o = 1. /i = 4, fi = 13, /j = 40,...
Metoda ta jest wygodna, gdy celem obliczeń jest kilka początkowych wyrazów ciągu {/„}•
19.4.2. Metoda rozkładu na ułamki proste
Zakładamy, że równanie charakterystyczne transformaty
amzm + am-1zm~1 + ... + a1z+a0 = 0 (19.27)
ma różne od zera pierwiastki jednokrotne zt, z2,zm. Jak wiadomo, funkcję wymierną (19.24) przy m > v można wyrazić wzorem
(19.28)
G(z) _ y Ak H(z) * = !Z-Z*’
gdzie Ak, k = 1, 2,..., m, są stałymi. Powyższe równanie mnożymy stronami przez z —Zj, i — 1,2,..., m, otrzymujemy
G(z)
H(z)
Z
k= 1 fc*i
(z Zj) = Ai + (z-Zj)
a po przejściu do granicy przy z->zt
A: = lim
G(z)(z-z,) H(z) ‘
Stąd po zmianie i na k oraz zastosowaniu reguły de THospitala, znajdujemy
G(z)(z—zk) G(zt)
Ak = lim
Z — Zk
H(z)
H'(zky
501
19.4. Przekształcenie odwrotne
Wykorzystując zależność (19.11), można wyrażenie (19.28) przedstawić w postaci
Akar{żl~le(n-l)},
H\Z) k= 1
a zgodnie z określeniem przekształcenia mamy 2t{Ą le(n — 1)} =
X
= V (n— wobec tego
n = 0
= Z ^ Z ** 1e(n— 1)^ " = X 2 "e(n-l) Z ^*1 *•
^(Z) k= 1 n = 0 n = 0 *=1
Wynika więc, że ogólny wyraz oryginału dyskretnego {/„} ma postać
* = 1
Ponieważ e(n — 1) = 1 dla n ^ 1 oraz e(n— 1) = 0 dla n = 0, otrzymujemy po wzięciu pod uwagę zależności (19.29) wzór przedstawiający ogólny wyraz oryginału {/„}:
(19.30)
r
f _ y G(zk) ,
przy czym /0 = 0.
Gdy wielomiany G(z) i H(z) w zależności (19.24) są tego samego stopnia (m = v), należy najpierw podzielić G(z) przez H(z), otrzymując w ten sposób /0 = bjam. Dowodzi się, że pozostałe wyrazy ciągu/, ,/2,... można obliczać na podstawie wzoru
Przykład 2. Obliczymy oryginał transformaty z przykładu 1.
Mamy tu G(z) = z2, //(z) = z2 —4z + 3 i H'(z) = 2z—4, a pierwiastkami równania W(z) = 0 są z, = 1 i :2 = 3. Ponieważ wielomiany G(z) i H(z) są tego samego stopnia, więc po podzieleniu G(z) przez H(z) znajdujemy f0 = 1. Przyjmując m = 2 we wzorze (19.30), otrzymujemy a stąd /t = 4, f2 = 13, /3 = 40,..., zgodnie z wynikiem otrzymanym w przykładzie 1.
Przypadku pierwiastków wielokrotnych nie będziemy rozpatrywać. Jest on omówiony w pracy [28],
19.5.1. Uwagi ogólne
Przy rozwiązywaniu równań różniczkowych na emc dokonuje się dyskretyzacji, rozpatrując wartości funkcji w regularnych chwilach tk, k = 0, 1, 2,..., o stałym