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