Mechatronika II, grupa 3

Ćwiczenie 3

Analiza widmowa – FFT

Zadanie 1)

% Dane wejściowe

okno=0.1;

fs=1000;

dt=1/fs;

t=0:dt:okno;

% Sygnał wyjściowy

y=10+100*cos(2*pi*50*t)+80*cos(2*pi*150*t+pi)+50*cos(2*pi*250*t)+40*cos(2*pi*430*t+pi/2);

% Wykreślenie wykresu sygnału badanego w funkcji czasu

y1=y;

subplot(4,1,1);

plot(t,y1);

% Przetwarzanie i wykreślenie sygnału wyjściowego szybką transformatą Fouriera

yf1=fft(y1);

% Wykreślenie postaci czasowej sygnału

Freq=(0:n);

subplot(4,1,2);

bar(Freq,yf1);

% Usunięcie części ujemnej wykresu poprzez wyznaczenie wartości bezwzględnej wartości funkcji

yf2=abs(yf1);

% Wykreślenie postaci częstotliwościowej sygnału w postaci słupków (które będą miały największe wartości przy częstotliwościach składowych sygnału)

subplot(4,1,3);

% Wykreślenie postaci częstotliwościowej przy przeskalowanej osi OX (osi częstotliwości)

plot(Freq(1:n/2),yf2(1:n/2));

subplot(4,1,4);

grid;

on;

f=Freq*fs/n;

stem(f(1:n/2),yf2(1:n/2);

Zadanie 2)