资源描述：
《离散时间系统时域及频域的分析》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、...页眉实验报告课程名称数字信号处理原理及实现实验名称离散时间系统的时域与频域分析班级与班级代码082510012实验室名称（或教室）3-804任课教师学号姓名：实验日期：2011年10月31日广东商学院教务处制....页脚...页眉姓名实验报告成绩评语：指导教师（签名）年月日说明：指导教师评分后，实验报告交院（系）办公室保存。....页脚...页眉实验二离散时间系统的时域与频域分析实验目的：通过实验，了解离散时间信号的傅立叶变换的数学依据；理解离散时间信号傅立叶变换在数字信号处理中的重要意义——即数字信号处理中的频域离散化的概念；掌握离散时间信号及系统的变换域分析——Z变换定义及应用；掌握
2、离散时间信号傅立叶变换的计算及应用。实验步骤和图形：1.(P66例2-27)解：MATALAB实现程序：Impseq文件：function[x,n]=impseq(n0,n1,n2)n=n1:n2;x=[(n-n0)==0];程序：>>b=[1,2,1];a=[1,-0.5,0.25];x=impseq(0,0,100);n=0:100;h=filter(b,a,x);stem(n,h,'.');title('脉冲响应波形');sum(abs(h))ans=6.5714[H,w]=freqz(b,a,200);magH=abs(H);....页脚...页眉phaH=angle(H);subpl
3、ot(2,1,1),plot(w/pi,magH);grid;xlabel('');ylabel('幅度');title('幅频响应')subplot(2,1,2);plot(w/pi,phaH/pi);grid;xlabel('频率（单位：pi）');ylabel('相位（单位：pi）');title('相频响应');n1=0:200;x1=5+3*cos(0.2*pi*n1)+4*sin(0.3*pi*n1);y1=filter(b,a,x1);subplot(2,1,1);plot(n1,x1);title('输入x1波形');subplot(2,1,2);plot(n1,y1);tit
4、le('输出y1波形');....页脚...页眉n1=0:200;x2=5+3*cos(0.2*pi*n1)+4*sin(0.8*pi*n1);y2=filter(b,a,x2);subplot(2,1,1);plot(n1,x2);title('输入x2波形');subplot(2,1,2);plot(n1,y2);title('输出y2波形');....页脚...页眉2.(P84例3-4，3-5)例3-4解：MATALAB实现程序：（1）四点序列的DTFT>>x=[1,1,1,1];w=[0:1:500]*2*pi/500;>>[H]=freqz(x,1,w);>>magH=abs(H);
5、phaH=angle(H);>>subplot(2,1,1);plot(w/pi,magH);grid>>xlabel('');ylabel('
6、X
7、');title('DTFT的幅度')>>subplot(2,1,2);plot(w/pi,phaH/pi*180);grid>>xlabel('以pi为单位的频率');ylabel('度');title('DTFT的相角')>>>>subplot(2,1,1);plot(w/pi,magH);grid实验结果：（....页脚...页眉（2）N=4;w1=2*pi/N;k=0:N-1;X=dft(x,N);magX=abs(X),phaX=ang
8、le(X)*180/pisubplot(2,1,1);plot(w*N/(2*pi),magH,'--');axis([-0.1,4.1,0,5]);holdonstem(k,magX);ylabel('
9、X(k)
10、');title('DFT的幅度：N=4');text(4.3,-1,'k')holdoffsubplot(2,1,2);plot(w*N/(2*pi),phaH*180/pi,'--');axis([-0.1,4.1,-200,200]);holdonstem(k,phaX);ylabel('度');title('DFT相角：N=4');text(4.3,-200,'k')实验结
11、果：....页脚...页眉例3-5解：MATALAB实现程序：x=[1,1,1,1];w=[0:1:500]*2*pi/500;[H]=freqz(x,1,w);magH=abs(H);phaH=angle(H);N=8;w1=2*pi/N;k=0:N-1;x=[x,zeros(1,4)];X=dft(x,N);magX=abs(X);phaX=angle(X)*180/pi;subplot(2,
