www.pudn.com > PCMsheji.rar > FFTSEQ.M

function [M,m,df]=fftseq(m,ts,df)  
%		[M,m,df]=fftseq(m,ts,df) 
%		[M,m,df]=fftseq(m,ts) 
%FFTSEQ	产生m的FFT变换M。 
%		序列补零以达到需要的频率分辨率df 
%		ts是抽样间隔. 输出df是最终的频率分辨率 
%		输出m是输入m补零后的结果。 
fs=1/ts; 
if nargin == 2 
  n1=0; 
else 
  n1=fs/df; 
end 
n2=length(m); 
n=2^(max(nextpow2(n1),nextpow2(n2))); 
M=fft(m,n); 
m=[m,zeros(1,n-n2)]; 
df=fs/n;