(数字信号处理实验.docxVIP

实验一信号、系统及系统响应1、实验目的认真复习采样理论、离散信号与系统、线性卷积、序列的z 变换及性质等有关内容；掌握离散时间序列的产生与基本运算，理解离散时间系统的时域特性与差分方程的求解方法，掌握离散信号的绘图方法；熟悉序列的z 变换及性质，理解理想采样前后信号频谱的变化。2、实验内容a. 产生长度为500 的在[0,1]之间均匀分布的随机序列，产生长度为500 的均值为0 单位方差的高斯分布序列。b. 线性时不变系统单位脉冲响应为h(n)=(0.9)nu(n)，当系统输入为x(n)=R10(n)时，求系统的零状态响应，并绘制波形图。c. 描述系统的差分方程为：y(n)-y(n-1)+0.9y(n-2)=x(n)，其中x(n)为激励，y(n)为响应。计算并绘制n=20,30,40,50,60,70,80,90,100 时的系统单位脉冲响应h(n)；计算并绘制n=20,30,40,50,60,70,80,90,100 时的系统单位阶跃响应s(n)；由h(n)表征的这个系统是稳定系统吗？d. 序列x(n)=(0.8)nu(n)，求DTFT[x(n)]，并画出它幅度、相位，实部、虚部的波形图。观察它是否具有周期性？e. 线性时不变系统的差分方程为y(n)=0.7y(n-1)+x(n)，求系统的频率响应H(ejω)，如果系统输入为x(n)=cos(0.05πn)u(n)，求系统的稳态响应并绘图。f. 设连续时间信号x(t)=e-1000|t|，计算并绘制它的傅立叶变换；如果用采样频率为每秒5000 样本对x(t)进行采样得到x1(n)，计算并绘制X1(ejω)，用x1(n)重建连续信号x(t)，并对结果进行讨论；如果用采样频率为每秒1000 样本对x(t)进行采样得到x2(n)，计算并绘制X2(ejω)，用x2(n)重建连续信号x(t)，并对结果进行讨论。加深对采样定理的理解。g. 设X1(z)=z+2+3z-1，X2(z)=2z2+4z+3+5z-1，用卷积方法计算X1(z)X2(z)。h. 已知系统方程为y(n)=0.9y(n-1)+x(n)，求系统函数H(z)并绘制其零极点图，求系统的频率响应H(ejω)并绘制其幅度和相位波形，求系统的单位脉冲响应h(n)并绘图。i. 系统方程为：y(n)-0.4y(n-1)+0.75y(n-2)=2.2403x(n)2.4908x(n-1)+2.2403x(n-2)，验证系统是否为线性系统、是否为时不变系统。3、实验报告要求a. 根据实验目的与内容，写出相应的计算过程，绘制各种波形图。b. 简述上机调试程序的方法，写出实验中遇到的问题及解决方法。c. 总结实验体会和感想。?实验二 DFT 和FFT1、实验目的认真复习周期序列DFS、有限长序列DFT 的概念、旋转因子的定义、以及DFS 和DFT的性质等有关内容；复习基2-FFT 的基本算法，混合基-FFT 的基本算法、Chirp-Z 变换的算法等快速傅立叶变换的方法。掌握有限长序列的循环移位、循环卷积的方法，对序列共轭对称性的含义和相关内容加深理解和掌握，掌握利用DFT 分析序列的频谱特性的基本方法。掌握FFT 算法的基本原理和方法、Chirp-Z 变换的基本原理和方法，掌握利用FFT 分析序列的频谱特性的方法。熟悉利用MATLAB 进行序列的DFT、FFT 的分析方法。2、实验内容a. 设周期序列( ) { ,0,1, 2,3,0,1, 2,3,0,1, 2,3, x?? n =?? ??}，求该序列的离散傅立叶级数X??(k) =DFS[x??(n)]，并画出DFS 的幅度特性。b. 设周期方波序列为1 ( 1)( ) ( 0, 1, 2, )0 ( ( 1) 1mN n mN Lx n mmN L n m N? ≤ ≤ + ?=?? + ≤ ≤ + ? = ± ±?? ??其中N 为基波周期，L/N 是占空比。(1) 用L 和N 求|X??(k)|的表达式；(2) 当L 和N 分别为：L=5，N=20；L=5，N=40；L=5，N=60 以及L=7，N=60 时画出DFS 的幅度谱；(3) 对以上结果进行讨论，总结其特点和规律。c. 设有限长序列x(n) ={0,1,2,3}，计算DTFT[x(n)]=X(ejω)，并画出它的幅度谱；然后利用12 , 0,1,2,34k kk πω = = 对X(ejω)进行采样，并证明它等于实验a 中的X??(k)。d. 序列x(n)=R4(n)，计算DTFT[x(n)]=X(ejω)，并绘制其幅度和相位谱。(1) 计算x(n)的4 点DFT，并绘制DFT 的幅度与相位谱；(2) 将x(n)补零形成8 点序列，计算8 点DFT，并绘制幅度与相位谱，求频率分辨率；(3) 将x(n)补零形成16 点序列