function?[xxyyy2]=zoom_fft（fifpnpnfftfsx）

if（（fp-fi）>（fs/np/2））
warning（‘頻帶范圍?is?larger?than?fs/np/2‘）;
end

%fifp?????????細化起始頻率，截止頻率
%np????????????放大倍數
%nfft??????????fft長度
%?fs????????????采樣頻率
%?x?????????????數據
nt=length（x）;??
fa=（fp+fi）/2-0.25*fs/np;???%最大細化起始頻率
fb=（fp+fi）/2+0.25*fs/np;???%最大細化截止頻率
nf=2^nextpow2（nt）;?????????%取大于nt的且最接近于nt的2的整數次方為暫時fft長度
na=round（0.5*nf/np+1）;

%頻移
n=0:nt-1;
b=n*pi*（fa+fb）/fs;?%確定單位旋轉因子向量（b=2*pi*ft*tft為細化頻帶中心頻率，即ft=（fa+fb）/2）
y=x.*exp（-i*b）;????%乘以單位旋轉因子進行頻移

%頻移后的低通濾波????????????????????
b=fft（ynf）;???????%fft變換
a（1:na）=b（1:na）;???%正頻率帶通內的元素賦值
a（nf-na+1:nf）=b（nf-na+1:nf）;%負頻率帶通內的元素賦值
b=ifft（anf）;??????%IFFT變換??

%重采樣
c=b（1:np:nt）;

%進行ZOOM-FFT
a=fft（cnfft）*2/nfft;?%進行點數為nfft的fft
%變換結果重新排序
??????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2133??2012-09-13?18:14??細化譜.m

-----------?---------??----------?-----??----

?????????????????2133????????????????????1