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 45737527
Najwięcej użytkowników 266 było obecnych 19 Lutego 2015, Czw 7:03 pm

Aktualnie online:




Najnowsze posty na forum:
Telewizja online za darmo  (24 Maj 2018, Czw 11:00 pm)
Ja, Godard 2017 Obejrzyj Cały Film Online Lektor PL  (23 Maj 2018, Sro 7:44 pm)
Wykres i widmo zerowej składowej  (23 Maj 2018, Sro 3:09 pm)
(s01e10) Nawiedzona panna młoda 1 Sezon 10 Odcinek Online  (22 Maj 2018, Wto 10:17 pm)
Nawiedzona panna młoda s01e09 transmisja na Ĺźywo  (22 Maj 2018, Wto 10:17 pm)
Nawiedzona panna młoda s01e09 gdzie obejrzeć online?  (22 Maj 2018, Wto 10:16 pm)
(s01e09) Nawiedzona panna młoda 1 Sezon 9 Odcinek Online  (22 Maj 2018, Wto 10:16 pm)
Nawiedzona panna młoda s01e08 gdzie obejrzeć online?  (22 Maj 2018, Wto 10:16 pm)
(s01e08) Nawiedzona panna młoda 1 Sezon 8 Odcinek Online  (22 Maj 2018, Wto 10:15 pm)
(s01e07) Nawiedzona panna młoda 1 Sezon 7 Odcinek Online  (22 Maj 2018, Wto 10:15 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