fftseq.m

matlab中对于语音编码的叙述和代码演示

function [M,m,df]=fftseq(m,ts,df)
% [M,m,df]=fftseq(m,ts,df)
% [M,m,df]=fftseq(m,ts)
%FFTSEQ 生成M, 它是时间序列m的FFT
% 对序列填充零，以满足所要求的频率分辨率df
% ts是采样间隔. 输出df是最终的频率分辨率
% 输出m是输入m的补过零的版本，M是FFT.
fs=1/ts;
if nargin == 2 %函数变量个数为2，表示 [M,m,df]=fftseq(m,ts)
n1=0;
else %否则函数变量个数为3，表示[M,m,df]=fftseq(m,ts,df)
n1=fs/df;
end
n2=length(m); %求序列长度
n=2^(max(nextpow2(n1),nextpow2(n2))); %将n1或n2值按2的幂次度量，取大者
M=fft(m,n); %得到FFT
m=[m,zeros(1,n-n2)]; %时间序列补零
df=fs/n; %得到最终的频率分辨率