sinifft.m

% Demonstroidaan kuinka Fourier-muunnos purkaa signaalin sen
% sinimuotoisiin komponentteihin

function sinifft
    % Luodaan periodinen signaali, jossa kaksi eritaajuista
    % sinimuotoista komponenttia
    t = 0:0.05:9.95;
    f1 = 1.2;
    f2 = 5.5;
    x = 2.3*sin(2*pi*f1*t) + 0.5*sin(2*pi*f2*t);

    % Näytteistysaika ja maksimitaajuus
    dt = 0.05;
    fmax = 1/(2*dt);

    % Esitetään signaali
    figure(1);
    L = length(x);
    plot( t, x );

    % Lasketaan Fourier-muunnos ja esitetään komponenttien suuruudet
    y = fft( x );
    freq = 2*fmax*(0:L/2)/L;
    figure(2);
    plot( freq, 2*abs(y(1:L/2+1))/L )
end