sy3.txt

几个Matlab编程示例

%3、离散信号的DTFT和DFT
n=0:15;
x=cos(5*pi/16*n);
y1=fft(x,16);
absy1=abs(y1);
subplot(3,1,1);stem(absy1);
title('16点DFT幅度');
y2=fft(x,32);
absy2=abs(y2);
subplot(3,1,2);stem(absy2);
title('32点DFT幅度');
k=0:500;w=(pi/500)*k;
y3=x*(exp(-j*pi/500)).^(n'*k);
absy3=abs(y3);
subplot(3,1,3);plot(w/pi,absy3);
title('DTFT幅度');