Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék

10. fejezet
Fourier analízis
Fourier analízist és frekvenciaspektrumot az fft(y,N) paranccsal végezhetünk, ahol N 2 egész számú hatványa
(mivel a Matlab gyors Fourier transzformációt alkalmaz). Inverz transzormáció az ifft paranccsal
számítható.
Például vegyünk egy 1000Hz-el mintavételezett jelet, mely egy 50 és egy 120 Hz-es komponenset tartalmaz,
de el van rontva egy véletlenszerű, zérusátlagú zajjal. Végezzük el a Fourier-transzformációt és számítsuk ki
a komplex amplitúdó abszolúrtékét (fft demo.m ).
t = 0:0.001:0.6;
x = sin(2*pi*50*t)+sin(2*pi*120*t);
y = x + 2*randn(size(t));
n=8;
Y = fft(y,2^n);
Pyy = Y.* conj(Y) / 2^n;
f = 1000*(0:2^(n-1))/2^n;
hossz=min(length(Pyy),length(f));
subplot(2,1,1)
plot(t,x), xlabel(’t’), ylabel(’x(t)’)
grid on, axis([0 0.6 -2.5 2.5])
subplot(2,1,2)
plot(f(1:hossz),Pyy(1:hossz))
xlabel(’f [Hz]’), ylabel(’amplitudo’)
grid on
fft demo.m
38