Mihai Sprinceana
Un forum de programare cu de toate. Va astept sa va inscrieti si sa deveniti moderatori. Oricine este binevenit aici sa se inscrie si sa aiba acces la informatie free! Fiecare este liber sa adauge proiecte programe free etc. Ajutati acest forum sa devina o comunitate puternica unde fiecare invata de la fiecare! Tot ce trebuie sa faceti este sa va inregistrati si fiecare contributie se poate dovedi utila in timp! Forumul este free informatia free dk aveti timp liber ajutati si pe ceilalti si invatati si voi in acelasi timp! Haideti sa facem ceva pt.a ne ajuta intre noi!
Cititi regulament postare forum inainte de a posta!
|
Lista Forumurilor Pe Tematici
|
Mihai Sprinceana | Inregistrare | Login
POZE MIHAI SPRINCEANA
Nu sunteti logat.
|
Nou pe simpatie:
BrianaS pe Simpatie
 | Femeie
24 ani
Bucuresti
cauta Barbat
26 - 57 ani |
|
mihaispr
Administrator
 Inregistrat: acum 17 ani
Postari: 2142
|
|
Script matlab pentru filtru trece banda (ceva mai mult decat un filtru trece jos-lowpass filter sau highpass filter-filtru trece sus).
%este un filtru trece banda in domeniul [0.12*pi, 0.24*pi] radians/esantion(sample)
%Script Fir.m -lab5 PNS -Caracteristicile filtrelor cu raspuns impulsional finit(FIR)
%Functia bode reprezinta raspunsul in frecventa Bode pentru modele LTI
clear
n=6;
a=-pi:pi/n:pi;
%h=[1 2 3 4 5 4 3 2 1];
h=cos(a);
h=h/sum(abs(h));
sis=tf(h,1,0.1);
bode(sis)
figure
[mag,phase,w]=bode(sis);
[nf,i]=size(w);
for i=1:nf
m(i)=mag(1,1,i);
ph(i)=phase(1,1,i);
end
subplot(2,1,1)
plot(w,m)
subplot(2,1,2)
plot(w,ph)
Pentru rularea scriptului Fir.m apasati tasta F5 sau scrieti in consola:
>>Fir %numele scriptului fara extensie!!!!
Efectul rularii scriptului in command window:
Raspuns pt.maginitudine script:
%Script magnitudine.m
n = 6;
a = -pi:pi/n:pi;
h = cos(a);
h = h/sum(abs(h));
fvtool(h,'normalizemagnitudeto1','on');
Efect:
|
|
| pus acum 15 ani |
|
mihaispr
Administrator
Inregistrat: acum 17 ani
Postari: 2142
|
|
Cum definim filtrele in matlab:
METODA 1:
FILTRU TRECE JOS:
>>d = fdesign.lowpass('N,Fc',5,0.1) ;
Hd = design(d, 'butter') ;
z = filter(Hd,y)
FILTRU TRECE SUS:
>>d = fdesign.highpass('N,Fc',5,0.1) ;
Hd = design(d, 'butter') ;
z = filter(Hd,y)
FILTRU TRECE BANDA:
>>d = fdesign.bandpass('N,Fc',5,0.1) ;
Hd = design(d, 'butter') ;
z = filter(Hd,y)
METODA 2:
[b,a]=butter(5,0.1);
Hdnew = dfilt.df2t(b,a);
Hdnew.arithmetic = 'fixed';
Hd = sos(Hdnew)
Pentru mai multe detalii despre folosirea filtrelor in matlab (lista completa functii dedicate) tastati in command window:
>>help filterdesign
|
|
| pus acum 15 ani |
|
