zfft的matlab实现

function [f y] = zfft（x f_start f_end fs）  
% x为采集的数据（待分析数据）  
% f_start为分析的起始频率  
% f_end为分析的截止频率  
% fs为采集数据的采样频率  
% f为输出的频率序列  
% y为输出的幅值序列（实数）  
  
fe  =  （f_start+f_end） / 2;                             % fe表示中心频率  
N   =  length（x）;                                        % 取数据长度  
r   =  0:N-1;                                            % 共N个点 
x1  =  x .* exp（-1j .* 2*pi*fe.*r ./ fs）;               % 复调制移频
bw  =  f_end-f_start;                                   % 所需细化频带带宽