function lydfilavspilling

m = 'Sound (WAV) file';
c = 'lydfil3.wav';
[m d] = wavfinfo(c);
% d % For utskrift om ønskelig
Fs = 44100;  % Samplingsfrekvens
type = 'int16';  % Ikke nødvendig?
[y, Fs, type] = wavread(c, [2001 12000]);
wavplay(y,Fs);  % NB: Velg sound(y,Fs); dersom du er på en Linux-maskin!

N = 10000;          % Velg et mindre tall under uttesting av programmet
                    % og dersom programmet stopper HELT opp (minst 2-3 min
                    % for N = 10 000.
T = 10000/Fs        % Sløyfer her semikolon for å få ut en testutskrift
t = linspace(0,T*(N-1)/N,N);
g = zeros(N,1);
delta_f = 1/T;
f_max = N/T
f = linspace(0,f_max*(N-1)/N,N);

g = y(:,1);  % MERK: y som vi får fra wavread kan være av typen y(N,2)
             % dersom lydfilen er i stereo. Denne programlinjen plukker 
             % ut bare en av kanalene slik at vi får et mono-signal i 
             % fortsettelsen.
% Plotter opprinnelig funksjon
plot(t,g,'-k');
%xlim([-0.2,6.5]);    % Kan begrense plottet dersom man ønsker det, men
%ylim([-1.2,1.2]);    % pass på at du ikke mister oversikten!
title('Opprinnelig tids-funksjon');
xlabel('Tid (sek)');
ylabel('Relativt signal');
figure;

% Fourier transform starter her

an = zeros(N,1);
bn = zeros(N,1);

for j=1:N
    sum_a = 0;
    sum_b = 0;
    delta_omega = (j-1)*2*pi/N;
    for i=1:N
        sum_a = sum_a + g(i)*cos(delta_omega*(i-1));
        sum_b = sum_b + g(i)*sin(delta_omega*(i-1));
    end;
    an(j) = sum_a.*(1/N);
    bn(j) = sum_b.*(1/N);
end;

plot(f,an,'.-r');
%xlim([-0.2,6.5]);
hold on;
plot(f,bn,'.-b');
title('Fourier-komponenter (cos-ledd: røde, sin-ledd: blå)');
xlabel('frekvens (Hz)');
ylabel('Relativt signal');




