zdft.m

用MATLAB实现的DSP函数

%DFT的计算
close all
t=linspace(1e-3,100e-3,10);
xn=sin(100*2*pi*t);
N=length(xn);
WNnk=dftmtx(N);Xk=xn*WNnk;
figure;
subplot(1,2,1),stem(1:N,xn),title('时域离散序列x(n)');
subplot(1,2,2),stem(1:N,abs(Xk)),title('x(n)的DFT变换结果');