okna, sita in viri

19.6 Uporaba programa MATLAB pri načrtovanju sit s FIR 159
izračunamo frekvenčni odziv pri frekvencah med 0 Hz in f s Hz. Pri tem je
vektor f dolžine l.
Zapisali smo le najpomembnejše možnosti funkcije freqz. Podrobnejši
opis vseh možnosti najdemo v “Pomoči”, ki je sestavni del paketa Matlab.
Opis freqz zaključimo s preprostim zgledom njene uporabe (zgled 19.2 na
strani 162).
MATLAB 19.1: Izračun in izris amplitudne karakteristike sita s FIR. Frekvenčna karakteristika naj bo izračunana v
512 točka za izbrano frekvenco vzorčenja.
...
% IZRAČUN FREKVENČNE KARAKTERISTIKE
[H,f]=freqz(h,1,512,fs);
% IZRIS AMPLITUDNE KARAKTERISTIKE
mag=20*log10(abs(H));
plot(f,mag),grid on;
xlabel(’frekvenca [Hz]’);
ylabel(’amplitudni odziv [dB]’)
% b = h: impulzni odziv sita
% a = a_0 = 1: sito s FIR
% 512: število točk računanja FFT
% fs: frekvenca vzorčenja
% izračun amplitude v dB
% izris grafa na podlago z mrežo
% oznaka horizontalne osi
% oznaka vertikalne osi
Okenska metoda
Postopek načrtovanja standardnih, frekvenčno selektivnih sit s FIR z linearno
fazno karakteristiko z uporabo okenske metode smo širše opisali v razdelku
na strani . Iz opisa sledi, da ga izvedemo v naslednjih korakih:
1. izberemo frekvenčno karakteristiko prototipnega sita,
2. izračunamo idealni odziv prototipnega sita (za idealna frekvenčno selektivna
sita so d[n] zbrani v tabeli 18.1 na strani 109),
3. izberemo okensko funkcijo w[n],
4. ocenimo število koeficientov sita N,
5. izračunamo koeficiente okenske funkcije w[n],
6. izračunamo koeficiente sita,
Za te korake ima Matlab v orodnem kovčku “Signal Processing Toolbox”
pripravljeno funkcijo fir1. Njena osnovna sintaksa je naslednja:
h = fir1(r,fm)
Vektor h vsebuje N koeficientov impulznega odziva sita z mejno frekvenco
f m. Parameter r določa red sita – ta je pri sitih s FIR za ena manjši od števila
datoteka: signal_C