zoomfft.m

《数字信号处理：基于计算机的方法（第三版）》资源下载第三部分。 著者： （美）米特拉著 作译者： 孙洪等译 ISBN号： 7-121-02652-X 出版日期： 2006-06 出版

function [XF] = zoomfft(X,R,K,I)
% ZOOMFFT(X,R,K,I) computes K samples of the N-point of the
% vector X beginning with the I-th sample where N = R*K
%
N = R*K;
X = [X zeros(1,N - length(X))];
y = zeros(K,R);
for r = 1:R
    y(:,r) = X(r : R : N)';
end
Y = fft(y);
XF = zeros(1,K);
for n = I:I+K-1
    for r = 1:R
    XF(n-I+1) = XF(n-I+1) + exp(-i*(n-1)*(r-1)*2*pi/N)*Y(mod(n-1,K)+1,r);
    end
end