%Czyszczenie pozostalosci
clc
clear
%Wprowadzenie krótkiego ciągu liczb
a = abs((floor(rand(1,5)*100))-60);

%wczytanie dzwieku
%a = waveread('sygnal1.wav');

N = length(a);
Wn=exp(j*2*pi/N);
A=zeros(1,N);

%Wprowadzona macierz
disp('Wprowadzone dane')
disp('a = ');
disp(a);

disp('1) Wyznaczenie DFT z definicji')
for n=1:N
	w=Wn^(n-1);
    for k=1:N
		A(n)=A(n)+a(k)*w^(k-1);
	end;
end;
disp('A = ');
disp(A);

disp('2) Wyznaczenie DFT z postaci macierzowej');

for i=0:N-1;
    for k=0:N-1;
        A1(i+1,k+1)=Wn^(i*k);
    end;
end;
A2 = [A1*(a')]';
disp('A1 = ');
disp(A2);

disp('3) Sprawdzenie rezultatow za pomoca funkcji fft()');
A3=[fft(a)];
disp('A2 = ');
disp(A3);