mydft.m

digital signal processing-MATLAB

function y=mydft(x)
%《数字信号处理教程——MATLAB释义与实现》
% 用简单的速度很低的方法计算DFT
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% y=mydft(x)
% ------------------------------------------------------------
% x为给定时间序列
% y为x的离散傅立叶变换
%
N=length(x);            % 输入序列的长度
n=0:N-1; k=n;           % 确定时域位置序列n和频域位置序列k
WN=exp(-j*2*pi/N);      % 计算DFT所需的旋转因子
nk=n'*k; WNnk=WN.^nk;   % 构成旋转因子矩阵
Xk=x*WNnk;y=Xk;         % 按DFT定义计算x的傅立叶变换