digitalny filter - uplne zaciatocnicka otazka :-)

[Daniel Valuch]

ano, ku IIR filtru som sa uz dostal, nasiel som nejake materialy na 
internete.
Tento algoritmus presne popisuje RC filter

A =RC/dt;
x=[zeros(1,10) ones(1,50) zeros(1,100)];
y_previous=0;
for n=1:length(x);
    y(n)=(A*y_previous+x(n))/(A+1);
    y_previous=y(n);
end;

snazim sa to teraz previest do standardneho zapisu pomocou zlomku 
prenosovej funkcie:

y(n) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb)
                 - a(2)*y(n-1) - ... - a(na+1)*y(n-na)

takze koeficient b(1) by mal byt 1/(A+1) a koeficient a(2) by mal byt A.

V matlabe je implementovana funkcia filter(b,a,x) ktora by to mala 
realizovat ako
           b(1) + b(2) z^-1 + b(3) z^-2 + ...
Y(z) = ----------------------------------- X(z)
           1     + a(2) z^-1 + a(3) z^-2 + ...

ale asi som trochu natvrdly lebo to nefunguje :-/ Neviem ale vytvorit 
citatel a menovatel aby filter spravne fungoval.

b.




Ondrej wrote:
>
>  Přesně RC článku může odpovídat pouze IIR (FIR má konečnou impulsovou 
> odezvu a RC nekonečnou - záludná otázka u státnic mimochodem).
> Nicméně pro návrh FIR lze použít např. fci FIR1, které se zadá jen řád 
> a zlomová frekvence. Pak je tam ještě zajímavé FIRCLS (nebo tak 
> nějak), které se zadá požadované pásmo a ona interpoluje vhodný FIR 
> filtr. Pro IIR jsou tam obdobné fce.
> No a když pak máte koeficienty, tak ty se předhodí fci FILTER spolu se 
> signálem a výstupem je vyfiltrovaný signál.
>
> OH
>