FIR filtro matematinio modelio gavimas

Susipažinti su FIR filtrais, pademonstruoti, kaip skaitmeniniai filtrai gali būti panaudoti įvairaus dažnio signalų filtravimui. Darbo metu bus nustatomi FIR filtro koeficientai ir perdavimo charakteristikos, filtruojami signalai.
K
filtruoti
2
1
650
1/4π
1
260
1/4π
40
Žemesnio dažnio dedamąją
1. Sugeneruoti amplitudžių A1 ir A2, dažnių f1 ir f2 fazės postūmių Θ1 ir Θ2 sinusinių sumos signalą. Diskretizuoti dažniu fs, naudoti N atskaitų. Naudojant funkciją fft() ir abs() apskaičiuoti sugeneruoto signalo spektrą. Viename lange atvaizduoti pirmas N1 šio signalo atskaitų naudodami funkciją plot() taip, kad horizontalioje ašyje būtų vaizduojamas laikas, ir signalo spektrą naudodami funkciją plot() taip, kad horizontalioje ašyje būtų atvaizduojami dažniai, kurių rėžiai yra nuo 0 iki fs/2.
x=linspace(0,1,8000);
a=1*sin(2*pi*(650)*x+(1/4)*pi);
b=1*sin(2*pi*(260)*x+(1/4)*pi);
c=a+b;
d=abs(fft(c,8000));
subplot(2,1,1)
plot(x(1:100),c(1:100));
title('Sinusiniu sumos sigalas')
e=(1:8000/2);
subplot(2,1,2)
plot(e,d(1:8000/2));
title('Sinusiniu sumos spektras')
Sinusinis sumos signalas
2. Naudojant funkciją fir1() surasti K-tos eiles FIR filtro koeficientus žemesnio dažnio dedamosios filtravimui. Naudojant funkcijas fregz() ir fregzplot() atvaizduoti sudarytos filtro amplitudinę ir fazinę perdavimo charakteristikas.
b=fir1(40,650/4000,'low');
[h,fn,s]=freqz(b,1,256,8000);
s.plot='both';
s.xunits='hz';
s.yunits='linear';
freqzplot(h,fn,s);
title('Filtro amplitudes perdavimo charakteristika')
xlabel('f')
ylabel('')
subplot(2,1,2)
title('Filtro fazine perdavimo charakteristika')
xlabel('f')
ylabel('faze')
3. Naudojant funkciją grpdelay() atvaizduoti filtro vėlinimo charakteristiką.
b=fir1(40,650/4000,'low');
grpdelay(b,1,256,8000);
title('Filtro velinimas')
xlabel('f')
ylabel('velinimas')
4. Atvaizduoti filtro koeficientus naudojant funkcija stem().
b=fir1(40,650/4000,'low');
stem(b)
title('Filtro koeficientai')
5. Apskaičiuoti filtro koeficientų N atskaitų Furjė transformaciją naudojant funkciją fft(). Viename lange atvaizduoti koeficientų amplitudinę dažninę charakteristika (abs(fft(x)) ir fazinę dažninę charakteristiką (unwrap(angle(fft(x)))*180/П) naudojant funkciją plot(), kad horizontalioje ašyje būtų atvaizduojami dažniai, kurių rėžiai yra nuo 0 iki fs/2.
b=fir1(40,650/4000,'low');
c=abs(fft(b,4000));
subplot(2,1,1);
plot((1:100),c(1:100));
title('Filtro koeficientu amplitudine daznine charakteristika')
xlabel('f')
ylabel('A')
d=unwrap(angle(fft(b,4000))*(180/pi));
subplot(2,1,2);
plot((1:100),d(1:100))
title('Filtro koeficientu fazine daznine charakteristika')
xlabel('f')
ylabel('faze')
6. Sugeneruoti N1 ilgio vienetinio impulso signalą. Viename lange atvaizduoti šį signalą ir jo spektrą. Signalą vaizduoti naudojant funkciją stairs(). Skaičiuojant spektrą naudoti N atskaitų, jį atvaizduoti naudojant funkciją plot() taip, kad horizontalioje ašyje būtų atvaizduojami...

Daugiau informacijos...

★ Klientai rekomenduoja

Šį rašto darbą rekomenduoja mūsų klientai. Ką tai reiškia?

Mūsų svetainėje pateikiama dešimtys tūkstančių skirtingų rašto darbų, kuriuos įkėlė daugybė moksleivių ir studentų su skirtingais gabumais. Būtent šis rašto darbas yra patikrintas specialistų ir rekomenduojamas kitų klientų, kurie po atsisiuntimo įvertino šį mokslo darbą teigiamai. Todėl galite būti tikri, kad šis pasirinkimas geriausias!