Next: 9. Cvičení 9: Adaptive
Up: 8. Cvičení 8: Odhad
Previous: 8.1 Kaskáda FIR prediktorů
Charakteristika prediktorů 2. řádu není příliš selektivní.
Velké zlepšení z hlediska selektivity lze dosáhnout umístěním komplexně
sdružených pólů v úhlech nul původního filtru poblíž jednotkové kružnice.
Přenosová funkce výsledného systému má pak tvar
 |
(8.8) |
V čitateli vidíme dvě komplexně sdružené nuly
a ve jmenovateli
dva komplexně sdružené póly
na stejných
úhlech
,
s modulem
.
Na Obr. 8.8
je vidět srovnání systémů s přenosovou funkcí (8.8)
pro různé moduly pólů.
Figure 8.8:
Porovnání frekvenčních charakteristik systémů s přenosovou funkcí
(8.8) pro různé volby
(jsou vyneseny i polohy nul a pólů)
|
Modul pólů
(ovlivňuje selektivitu) budu voli fixní.
Jediná věc, která lze upravovat je opět pouze úhel nuly
případně
. Váhy zavedu tedy shodně s (8.2).
Z přenosové fce (8.8) jednoduše přejdeme ke struktuře.
Kaskáda IIR prediktorů 2. řádu je na Obr. 8.9.
Figure 8.9:
Kaskáda IIR prediktorů 2. řádu pro odhad dvou frekvencí
 |
Z Obr. 8.9 jednoduše sesvíme rovnice filtrace
(oproti (8.3) přibydou zpětné vazby)
K tomu, abychom sestavili rovnice pro úpravu
a
potřebujeme opět určit derivace
,
.
Opět začneme s
.
Zpětnou vazbou se situace značně komplikuje, ale za předpokladu,
že
,
a
lze považovat za přibližně
shodné
(dostatečně malé
)
dospějeme k rekurzivnímu vztahu pro
(
označím
)
![$\displaystyle \frac{\partial e[n]}{\partial c_2[n]} = d_2[n] = - 2 x_2[n-1] + 2 r e[n-1] + 2 r c_2[n] d_2[n-1] - r^2 d_2[n-2] .$](img642.png) |
(8.10) |
Derivaci
stanovíme opět
s použitím LTI předpokladu. Prohodíme opět oba prediktory jako na
Obr. 8.2
a napíšeme rovnice filtrace
a již obdobně jako v (8.10) (předpoklad pomalu měnících se vah)
určím problematickou derivaci
![$\displaystyle \frac{\partial e[n]}{\partial c_1[n]} = d_1[n] = - 2 y_1[n-1] + 2 r e[n-1] + 2 r c_1[n] d_1[n-1] - r^2 d_1[n-2] .$](img644.png) |
(8.12) |
Obrázek výsledné struktury lze opět převzít, viz. Obr. 8.3.
S použitím (8.12) a (8.10)
již jednoduše sestavíme rovnice pro
úpravu vah
Cvičení 8.2: Vystupní signál
modelujte shodně jako ve cvičení 8.1.
Napište skript v Matlabu implementující strukturu Obr. 8.9
(kaskádu IIR prediktorů 2. řádu pro odhad dvou frekvencí)
Porovnejte výsledky (odhady frekvencí na čase, atd.)
oproti výsledkům získaným ve cvičení 8.1
(kaskádu FIR prediktorů 2. řádu pro odhad dvou frekvencí).
Výsledky:
Figure 8.10:
Kaskáda IIR, rozptyl adidivního rušení
: (zleva doprava)
amplitudová spektra vstupu
a výstupu
,
průběh odhadu frekvencí na čase, modulové frekvenční charakteristiky
kaskády a obou prediktorů 2. řádu, polohy nul a pólů přenosové funkce
kaskády v z-rovině.
|
Next: 9. Cvičení 9: Adaptive
Up: 8. Cvičení 8: Odhad
Previous: 8.1 Kaskáda FIR prediktorů
Mirek
2006-12-12