kuolemat.m

% Käytetään Fourier-muunnosta datan analysointiin. Tiedosto
% kuolemat.txt sisältää tapaturmaisesti vuosina 1973-1978
% Yhdysvalloissa kuolleiden määrät kuukausittain.
function kuolemat
    % Ladataan data ja esitetään se graafisesti
    load kuolemat.txt
    figure(1)
    plot( kuolemat(:,1), kuolemat(:,2) )

    % Määritetään aikaväli ja maksimitaajuus
    dt = kuolemat(2,1)-kuolemat(1,1)
    fmax = 1/(2*dt)

    % Datam pituus ja keskiarvo
    L = size(kuolemat,1)
    ave = sum( kuolemat(:,2) ) / L;

    % Fourier-muunnos datasta, josta keskiarvo on vähennetty pois.
    y = fft( kuolemat(:,2)-ave );
    freq = 2*fmax*(0:L/2)/L;
    figure(2);
    plot( freq, 2*abs(y(1:L/2+1))/L )
end