细化fft算法

%细化FFT 
clear 
close all hidden 
format long 
 
fni=input('细化FFT处理-输入数据文件名 :','s'); 
fid=fopen(fni,'r');         %以只读方式打开数据文件 
sf=fscanf(fid,'%f',1);      %读入采样频率值 
fi=fscanf(fid,'%f',1);      %最小细化截至频率 
np=fscanf(fid,'%d',1);      %放大倍数 
nfft=fscanf(fid,'%d',1);    %FFT长度 
x=fscanf(fid,'%f',[1,inf]); %读入时程数据存成行向量 
status=fclose(fid);         %关闭数据文件 
 
nt=length(x); 
%最大细化截至频率( 
fa=fi+0.5*sf/np; 
%取大于n且最接近n的2的整数次方为FFT的长度 
nf=2^nextpow2(nt); 
%确定细化带宽的数据长度 
na=round(0.5*nf/np+1); 
 
%频移 
%建立一个按1递增的向量 
n=0:nt-1; 
%确定单位旋转因子向量(b=2*pi*ft*t,ft为细化频带中心频率，即ft=(fi+fa)/2) 
b=n*pi*(fi+fa)/sf; 
%乘以单位旋转因子进行频移 
y=x.*exp(-i*b);