exzfft_m.m

来自「在谱分析中」· M 代码 · 共 20 行

function y=exzfft_m(x,fi,fs,nfft,np)
% x 被测信号，被测信号长度要>=nfft*np
% fi  细化的最低频率
% fs  采样频率
% nfft 作细化FFT长
% np 放大倍数
% y  细化FFT输出
nt=length(x);     %计算读入数据长度
fa=fi+0.5 * fs/np;     %最大细化截止频率
%nf= 2^nextpow2(nt);     %取大于nt且最接近nt的整数次方为FFT长度
%na=round(0.5 * nf/np+1);     %确定细化带宽的数据长度
% 频移
n=0: nt-1;     %建一个递增向量
b=n*pi* (fi+fa)/fs;     %乘单位旋转因子进行频移
y=x.*exp(-i*b);
%  滤波和下采样
c=decimate(y,np); 
% FFT
y=fft(c, nfft) * 2/nfft;