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

Aktualnie online:




Najnowsze posty na forum:
Zlecę wykonanie przenośnika taśmowego w MATLAB simulink  (21 Lutego 2018, Sro 2:00 pm)
BOTOKS - 4 ODCINKI SERIALU [ZA DARMO]  (21 Lutego 2018, Sro 8:40 am)
Top Gear {S25E01} online napisy PL 25x01 CDA/Zalukaj  (20 Lutego 2018, Wto 9:54 pm)
Falling Water {S02E08} online napisy PL 2x08 CDA/Chomikuj  (20 Lutego 2018, Wto 9:50 pm)
The Seven Deadly Sins {S02E07} online napisy PL 2x07 CDA  (20 Lutego 2018, Wto 9:46 pm)
Waco {S01E05} online napisy PL 1x05 CDA/Zalukaj/Chomikuj  (20 Lutego 2018, Wto 9:41 pm)
The Magicians (Magicy) {S03E07} online napisy PL 3x07 CDA  (20 Lutego 2018, Wto 9:37 pm)
Stargate: Origins {S01E04} online napisy PL 1x4 CDA/Zalukaj  (20 Lutego 2018, Wto 9:32 pm)
Wojny magazynowe {S11E21} online napisy PL 11x21 CDA  (20 Lutego 2018, Wto 9:29 pm)
Katastrofy w przestworzach {S18E02} online napisy PL 18x02  (20 Lutego 2018, Wto 9:23 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