„проектиране на филтри
Аналогово-цифрово билинейно
Изтегляне 0.94 Mb.
|
| 3.Аналогово-цифрово билинейно z-преобразуване.
Формулата, по която се изчисляват цифровите честоти от дадени аналогови е: , където: - е кръговата честота на аналоговия филтър; - – период на дискретизация, съответно е честотата на дискретизация; - – кръгова честота на цифровия филтър; - Fa – аналогова честота ; - f – цифрова честота. След извършената апроксимация намираме цифровите честоти които имат стойности: Wp3=[60 2885]/Fn; %гранични честоти на цифров режекторен филтър Ws3=[299 2624]/Fn; %честоти на задържане на цифров режекторен филтър След билинейно Z-преобразуване чрез функцията bilinear(), при честота на дискретизация Fn=10.195kHz, трансформираме аналоговата лентова ПФ в цифрова. Тук получаваме стойности на N(z) и D(z), подреждаме ги по степени на z и HBS(z) следва да бъде: 4.ПНД на аналогов и цифров РФ. За изчертаването на ПНД на аналогов РФ се използва командата zplane (Ns,Ds), a за изчертаването на ПНД на цифров РФ – zplane(Nz,Dz). За извеждането на числените стойности на полюсите и нулите се използва командата tf2zp; 5.Криви на затихване на филтрите. 6.Амплитудно-честотни характеристики на аналогов и цифров РФ. За начертаване графиките на затихванията на двата филтъра ще са ни нужни следните команди[2]: [T1,w]=freqs(Ns2,Ds2); //за изчисляване съответно предавателната функция и кръговата честота на аналоговия РФ; [T2,w2] = freqz(Nz,Dz,Fs);// за изчисляване съответно предавателната функция и кръговата честота на цифровия РФ; m=abs(H); //за изчисляване на амплитудно-честотната характеристика в линейни единици; md=20*log10(m); //за изчисляване на АЧХ в dB; plot(w2*Fs/(2*pi) ,md2);// изчертаване на цифровата АЧХ; plot(w1/(2*pi),md2);//изчертава аналоговата АЧХ; Изтегляне 0.94 Mb. Сподели с приятели:
|