Informacijos apdorojimo sistemos: FIR skaitmeniniai filtrai LD4

Darbo tikslas. 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.
Darbo užduotis. Dėstytojas turi paskirti jums variantą. Naudodami MATLAB paketą sugeneruokite nurodytų parametrų sinusinių sumos signalą. Sudarykite FIR filtrą, kuris filtruotų vieną iš sinusinio signalų dedamųjų ir paliktų kitą. Nustatykite filtro charakteristikas. Sudarytų filtrų filtruokite sinusinių sumos, vienetinio impulso ir atsitiktinį signalus. Atvaizduokite ir palyginkite pradinius ir filtruotus signalus bei jų spektrus. Baigę darbą, pasikvieskite dėstytoją ir pademonstruokite atliktus darbus. Dėstytojas pasižymės, kad atlikote darbą.
6 variantas:
Darbo eiga.
1. Sugeneruokite amplitudžių A1 ir A2, dažnių f1 ir f2 fazės postūmių Θ1ir Θ2 sinusinių sumos signalą. Diskretizuokite dažniu fs, naudokite N atskaitų. Naudodami funkcijas fft() ir abs() apskaičiuokite sugeneruoto signalo spektrą. Viename lange atvaizduokite 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.
% ----------------------------------------------------------
% Pirma uzduotis
% Atliko Ricardas Marozas Vit-4
% ------------------------------------------------------------
clear;
clc;
N=1050; % signalo atskaitu kiekis
N1=45; % vaizduojamos pirmos N1 signalo atskaitu
fs=16000; % diskretizavimo daznis
A1=3; % pirmo signalo amplitude
f1=850; % pirmo signalo daznis
Q1=(5/4)*pi; % pirmo signalo fazes postumis
A2=3; % antro signalo amplitude
f2=380; % antro signalo daznis
Q2=(5/4)*pi; % antro signalo fazes postumis
n=(0:N-1); % atskaitu indeksai
t=n/fs; % laiko momentai
% -----------------------------------------------------------
x=A1*sin(2*pi*(f1/fs)*n + Q1) + A2*sin(2*pi*(f2/fs)*n + Q2);
subplot(2,1,1); plot(t, x);
xlabel('t'), ylabel('x(t)')
title('Sinusiniu sumos signalas')
axis([t(1) t(N1) min(x) max(x)])
xsp=abs( fft(x));
f=linspace(0, fs, N); % linspace(pirmas, paskutinis, kiekis)
subplot(2,1,2); plot(f, xsp);
xlabel('f'), ylabel('|A|')
title('Sinusiniu sumos signalo spektras')
axis([f(1) fs/2 min(xsp) max(xsp)+0.1*max(xsp)])
2. Naudodami funkciją fir1() suraskite K-tos eiles FIR filtro koeficientus nurodytos dažnio dedamosios filtravimui. Naudodami funkcijas freqz() ir freqzplot() atvaizduokite sudaryto filtro amplitudinę ir fazinę perdavimo charakteristikas.
% ----------------------------------------------------------
% Antra uzduotis
% Atliko Ricardas Marozas Vit-4
% ------------------------------------------------------------
clear;
clc;
N=1050; % signalo atskaitu kiekis
N1=45; % vaizduojamos pirmos N1 signalo atskaitu
fs=16000; % diskretizavimo daznis
A1=3; % pirmo signalo amplitude
f1=850; %...

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!