%****************************************************************
% File:  FOURIER.M
%
% Scopo: Creazione matrice di Fourier per FFT
%
% Uso:	f = fourier(m,n)
%
% Input:   m = 2^r radice m-esima unita'
%	     n	   coefficienti eventuale sviluppo
%		   se omesso n=m.
%
% Output:  f   matrice di Fourier
%
% Richiami:
% Commenti:
% Autori: mf
%*****************************************************************

function f=fourier(m,n)
if nargin == 1
   n=m;
end
i=sqrt(-1);
w=exp(2*pi*i/m);
vw=w.^[0:m-1]';
f=ones(m,n);
f(:,2)=vw;
for j=3:n;
   f(:,j)=f(:,2).*f(:,j-1);
end
return