dtft.m

给出了数字信号处理课程设计的基本程序

function [w,xn]=DTFT(a)
p=length(a);
n=[0:p-1];i=1;
for w=-2*pi:2*pi/p:2*pi;
w0(i,:)=exp(-j*w.*n);
i=i+1;
xn(i)=dot(a,wo(i,:));
end
stem([-2*pi:2*pi/n:2*pi],xn,'.')