fft_arb.m

来自「beamlet变化的工具箱」· M 代码 · 共 35 行

function fhat = fft_arb(x)
% fft_arb -- FFT for arbitrary length sequence
%	fhat = fft_arb(x)
% Input
%	x		array(1,p)
% Output
%	fhat 	FFT of length p
%
p = length(x);
l = log2(p);
l1 = ceiling(l);
n = 2^l1;
y = zero(1,n);
y(1:p) = x(:);
fy = fft(y);
del = p/n;
fhat = ufrint(fy,1,del);
fhat = fhat(1:p);

% Copyright 2000 David L. Donoho

%
% Part of BeamLab Version:200
% Built:Friday,23-Aug-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%
%
% Part of BeamLab Version:200
% Built:Saturday,14-Sep-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%