okna, sita in viri
okna, sita in viri
okna, sita in viri
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
19.6 Uporaba programa MATLAB pri načrtovanju sit s FIR 161<br />
Rezultati, ki jih dobimo pri načrtovanju s programom Matlab, se lahko<br />
razlikujejo od rezultatov, ki jih dobimo pri drugih podobnih programih. Razlike<br />
lahko nastanejo na primer zaradi drugačnega skaliranja. Tako MATLAB<br />
priredi amplitudi vrednost 1 na sred<strong>in</strong>i prepustnega pasu. To skaliranje lahko<br />
izključimo z dopisom besede noscale:<br />
b = fir1(r,fm,nonscale)<br />
Ker je tudi računanje koeficientov okenskih funkcij specifično, mnogokrat<br />
vmesni rezultati niso direktno prenosljivi v druge programe. Zato moramo<br />
biti pazljivi pri uporabi rezultatov oziroma moramo pred njihovo uporabo v<br />
navodilih preveriti konkretne algoritem računanja.<br />
ZGLED 19.6.1 (Hamm<strong>in</strong>govo okno)<br />
V zgledu 19.3.1 smo “peš” določili koeficiente impulznega odziva aproksimacije idealnega<br />
nizkopasovnega <strong>sita</strong>, ki izpolnjuje naslednje zahteve:<br />
mejna frekvenca idealnega <strong>sita</strong>:<br />
šir<strong>in</strong>a prehodnega pasu:<br />
slabljenje zapornega pasu:<br />
1,5 kHz<br />
0,5 kHz<br />
> 50 dB<br />
frekvenca vzorčenja: 8 kHz<br />
Ponovimo postopek izračuna s pomočjo programa Matlab <strong>in</strong> orodnega kovčka “Signal<br />
Process<strong>in</strong>g Toolbox”.<br />
REŠITEV: Izbiro <strong>okna</strong>, določitev števila koeficientov <strong>in</strong> mejne frekvence naredimo<br />
enako kot v zgledu 19.3.1. Iz tabele glede na šir<strong>in</strong>o prehodnega pasu <strong>in</strong> zahtevanega<br />
dušenja zapornega pasu izberemo Hamm<strong>in</strong>govo okno. Iz podatka za šir<strong>in</strong>o<br />
prehodnega pasu, ki je zapisan v tabeli 19.2, izračunamo število koeficientov <strong>sita</strong>:<br />
⌈<br />
N = 2π 3,3 ⌉ ⌈<br />
⌉<br />
3,3<br />
= 2π<br />
= ⌈52,8⌉ = 53 .<br />
∆Ω 2π ·0,0625<br />
Določimo še mejno frekvenco aproksimacijskega <strong>sita</strong>:<br />
f c = f c + ∆ f /2 = 1,5 + 0,25 = 1,75 kHz<br />
f m = 2 f c / f s = 2.1,75/5 = 0,7 .<br />
Koeficientov <strong>sita</strong> izračunamo s funkcijo fir1. Za primerjavo izračunamo še koeficiente<br />
<strong>sita</strong>, ki ga dobimo pri uporabi pravokotnega <strong>okna</strong> <strong>in</strong> koeficiente Hamm<strong>in</strong>govega <strong>okna</strong>.<br />
Za vse to zadostuje naslednji programom:<br />
datoteka: signal_C
Change language