fftseq.m

pcm编码设计和最佳量化器设计仿真程序。建议在matlab 6.1 或以上的版本

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;