Zaloguj się
Blog na Matlablog
Forum polskich użytkowników
 
UŻYTKOWNICY GRUPY PROFIL Zaloguj się, by sprawdzić wiadomości FAQ
 



Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Rysowanie wykresu filtru cyfrowego typu IIR
Forum MATLAB Strona Główna-> Przetwarzanie sygnałów
Post Wysłany: 1 Sierpnia 2007, Sro 10:54 am Temat postu: Rysowanie wykresu filtru cyfrowego typu IIR Odpowiedz z cytatem
 
AUTOR:
BIG
Początki pisania


Dołączył: 18 Cze 2007
Posty: 10


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam mam problem z narysowaniem odpowiedzi impulsowej filtru cyfrowego typu IIR.
Odpowiedz impulsowa filtru cyfrowego zaprojektowanego metodą niezmienności odpowiedzi impulsowej równa jest :



gdzie s1,s2 to pierwiastki zespolone, zaś T jest to okres próbkowania

np. s1=-0.5000 + 0.8660i; s2=-0.5000 - 0.8660i ; T= 0.8;

Chciałbym się dowiedzieć w jaki sposób mogłbym narysowac taki filtr cyfrowy ? Ja próbowałem w następujący sposób ale coś mi nie wychodzi.W kodzie programu umieściłem jeszcze filtr analogowy o nastęującej transmitancji operatorowej HU(s) , który wychodzi dobrze.

HU(s)=(1/LC)/(s^2+s*(R/L)+1/(L*C))

Kod:
% odp. impulsowa układu II rzedu HU(s)=(1/LC)/(s^2+s*(R/L)+1/(L*C))
L=1;
C=1;
R=1;

T=0.2;
N=50;
tobs=N*T;

    
s1=-R/(2*L)+((R/(2*L))^2-1/(L*C))^(1/2);
    
s2=-R/(2*L)-((R/(2*L))^2-1/(L*C))^(1/2);
            
t=(0:T:tobs-T);
    
    
tc=(0:T:tobs-T);

a = 1/(L*C);
a1 = R/L;
a2 = 1/(L*C);

La = [a];
Ma = [1 a1 a2];

ya=impulse(La,Ma,t);

    
b = 1/(L*C*(-s1+s2));
    
c = (exp(s2*T)-exp(s1*T))*b;
    
c1= -exp(s1*T)-exp(s2*T);
    
c2 = exp(s1*T)*exp(s2*T);

    
Lc = [0 T*c 0];   
    
Mc = [1 c1 c2];

    
yc = Dimpulse(Lc,Mc,N); 
 
subplot(211);
plot(t,ya,'-',tc,yc,'.r');


Dziękuję za wszelką pomoc.
Pozdrawiam...


 

Forum MATLAB Strona Główna-> Przetwarzanie sygnałów
Wyświetl posty z ostatnich:   

Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Wszystkie czasy w strefie CET (Europa)

Skocz do:  

Statystyki forum:



Od dnia 08.06.2006 forum odwiedzano 42050647
Najwięcej użytkowników 266 było obecnych 19 Lutego 2015, Czw 7:03 pm

Aktualnie online:




Najnowsze posty na forum:
Jaskiniowiec 2018 online PL dubbing/lektor CAŁY FILM CDA  (19 Lutego 2018, Pon 8:11 pm)
Nić widmo 2018 online PL lektor/napisy CAŁY FILM CDA  (19 Lutego 2018, Pon 8:08 pm)
W ulamku sekundy 2018 online PL lektor/napisy CAŁY FILM CDA  (19 Lutego 2018, Pon 8:04 pm)
Wicher – dzikie konie online PL lektor/napisy CAŁY FILM C  (19 Lutego 2018, Pon 8:01 pm)
Maz, czy nie maz online PL napisy/lektor CAŁY FILM CDA  (19 Lutego 2018, Pon 7:55 pm)
Sludzy Diabla 2018 online PL lektor/napisy CAŁY FILM CDA  (19 Lutego 2018, Pon 7:52 pm)
Jeszcze nie koniec 2018 online PL napisy/lektor CAŁY FILM  (19 Lutego 2018, Pon 7:48 pm)
Zona czy maz 2018 online PL napisy/lektor CDA CAŁY FILM  (19 Lutego 2018, Pon 7:44 pm)
Szron 2018 online PL lektor/napisy CAŁY FILM CDA  (19 Lutego 2018, Pon 7:41 pm)
Dopasowanie danych z dwĂłch 'cell aray'  (19 Lutego 2018, Pon 5:33 pm)
Twoje prawa:
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach
Nie możesz załączać plików na tym forum
Nie możesz ściągać plików na tym forum