%****************************************************************
% 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
|
