zfft.m

ZFFT轨道信号上应用的一些文档

%%%%%%%%%%%%%%%%%%ZFFT过程%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%{
ZFFT过程
（1）将关心频段搬至零点处；
（2）进行>频段长度(M)点的IFFT；
（3）将IFFT后序列做N点的FFT。放大倍数：N/M 
IFFT过程
（1）X(m)取共轭;(2)FFT;(3)取共轭;(4）除以长度length(X)
%}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ff,dd,mag]=zfft(yy,zfft_len,fmin,fs);
%yy为待ZFFT信号；
%zfft_len为ZFFT长度；
%fmin为ZFFT起始点；
%fs为fft点数
%ff为恢复的横坐标(单位HZ)；
%dd为ZFFT频谱
%mag为放大倍数
d=1:1:zfft_len;
for i=1:zfft_len        %%%yy中从fmin开始的点，去共轭后搬到d中
    d(i)=conj(yy(i+fmin));
end
fft_num=2^nextpow2(zfft_len);  %%%ifft点数，取接近zfft_len的，大于zfft_len的2的指数。
mag=fs/fft_num;
d=fft(d,fft_num);              
d=conj(d)/fft_num;     
dd=fft(d,fs); %%%%fs为fft点数
figure
plot(abs(dd));
f=(0:length(dd)-1)*(length(dd)-1)/length(dd);  %%%恢复ZFFT横坐标 
ff=f/mag+fmin;