如果您无法下载资料，请参考说明：
1、部分资料下载需要金币，请确保您的账户上有足够的金币
2、已购买过的文档，再次下载不重复扣费
3、资料包下载后请先用软件解压，在使用对应软件打开
1.有限长序列的DTFT>>w=[pi/8pi/4pi*3/8pi/2pi*5/8pi*3/4pi*7/8pi];>>num=[13579111315171921];>>h=freqz(num,1,w);>>subplot(2,2,1);>>stem(w/pi,abs(h));>>title('幅度谱')>>xlabel('数字频率');>>ylabel('振幅');>>subplot(2,2,2);>>stem(w/pi,real(h));>>title('虚部');>>xlabel('数字频率');>>ylabel('振幅');>>subplot(2,2,3);>>stem(w/pi,imag(h));>>title('实部');>>xlabel('数字频率');>>ylabel('振幅');>>subplot(2,2,4);>>stem(w/pi,angle(h));>>title('相位谱');>>xlabel('数字频率');>>ylabel('相位');2.有限长序列的FFT和DTFT>>w=[pi/8pi/4pi*3/8pi/2pi*5/8pi*3/4pi*7/8pi];>>n=0:15;%16点的算法>>num=cos(5*pi*n/6);>>h=freqz(num,1,w);>>subplot(2,2,1);>>stem(w/pi,abs(h));>>title('DTFT幅度谱')>>xlabel('数字频率');>>ylabel('振幅');grid;>>n=0:15;>>x=cos(5*pi*n/6);>>k=fft(x);>>subplot(2,2,2);>>stem(n,abs(k));>>title('FFT幅度谱')>>xlabel('数字频率');>>ylabel('振幅');grid;%DTFT和FFT算出的幅值轮廓一致，不过FFT的数字频率0~2π，DTFT只是0~π%32点的算法同上
立即下载
