f_zfft.m

利用频率细化方位估计,水声信号处理中的一项技术

function [y,zf]=f_zfft(xn,N,fs,f_step,f_down,f_up)
% ZFFT 对xn在某一频率段内进行细化傅立叶变换
% 语法:
% y=f_zfft(xn,N,f_step,f_down,f_up)
%
% 其中，XN是一个离散信号，F_DOWN是频率下限，
% F_STEP是步长，F_UP是频率上限，fs是采样频率，N是采样点数
% Y是得到XN在ZF内的细化傅立叶变换

%代码行
f_length=(f_up-f_down)/f_step+1;                           %设定频带长度
zf=f_down:f_step:f_up;

%对a(f),b(f)表达式进行循环
for f=1:f_length                                            
    a(f)=0; b(f)=0;
    for k=0:N-1
        a(f)=a(f)+xn(k+1)*cos(2*pi*k*zf(f)/fs);
        b(f)=b(f)+xn(k+1)*sin(2*pi*k*zf(f)/fs);
    end
    a(f)=a(f)*2/N;
    b(f)=b(f)*2/N;
end
%输出幅值谱矢量表达式
y=a-i*b;