dtft1.m

matlab离散信号时域分析

%分别计算16点序列的x[k]=cos5*pi*k/16   0<=k<=15 的16点和32点的DFT，绘出幅度图形
%并绘出该序列的DTFT的图形

k=0:15;
x=cos(5*pi.*k/16);
N=16;
N2=32;
k2=0:31;
%计算DFT

h=fft(x,N);
h2=fft(x,N2);
k3=linspace(0,2*pi,1024);

h3=freqz(x,[1],k3);





subplot(3,1,1);
stem(k/16,abs(h)); 

title('DFT 16点幅度特性');

subplot(3,1,2);
stem(k2/32,abs(h2),'-b'); 

title('DFT32点T幅度特性');

subplot(3,1,3);

plot(k3/2*pi,abs(h3),'-b'); 

title('DTFT幅度特性');