mydft.m

用与语音信号处理中的预处理中MATLAB原代码

function y=mydft(x)
% 按DFT定义计算x的傅立叶变换
N=length(x);

n=0:N-1;%确定时域位置序列n

k=n;%确定频域位置序列k
WN=exp(-j*2*pi/N); %设定WN因子  
nk=n'*k; %产生一个含nk值的整数矩阵
WNnk=WN.^nk;   %求出W矩阵
Xk=x*WNnk;y=Xk; %    求出离散傅立叶级数系数