2022年实验二离散系统分析.pdf

《2022年实验二离散系统分析.pdf》由会员分享，可在线阅读，更多相关《2022年实验二离散系统分析.pdf（12页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、本科学生验证性实验报告学号 9 姓名张静学院 物电学院专业、班级 10 电子实验课程名称数字信号处理实验教师及职称杨卫平教授开课学期 2013 至 2013 学年下学期填报时间 2013 年 5 月 23 日云南师范大学教务处编印精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 1 页，共 12 页 - - - - - - - - - - 实验序号2实验名称离散系统分析实验时间实验室同析 3-313一实验预习1实验目的深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意

2、义，熟练掌握利用MATLAB 分析离散系统的时域响应、频响特性和零极点的方法。掌握利用DTFT和 DFT确定系统特性的原理和方法。2实验原理、实验流程或装置示意图MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。1. 离散系统的时域响应离散时间 LTI 系统可用如下的线性常系数差分方程来描述：已知输入信号xk 以及系统初始状态y-1， y-2 .， 就可以求出系统的响应。MATLAB提供了计算差分方程的数值解的函数，可以计算上述差分方程描述的离散系统的单位脉冲响应、零状态响应和完全响应等。在调用 MATLAB 函数时， 需要利用

3、描述该离散系统的系数函数。对差分方程进行Z 变换即可得系统函数：在 MATLAB 中可使用向量a 和向量 b 分别保存分母多项式和分子多项式的系数：这些系数均从z0 按z的降幂排列。离散系统的单位冲激响应hk 的计算 h,k = impz(b, a)：计算系统的单位脉冲响应hk 和相应的时间向量k；也可简写为： h = impz(b, a)。其中：)1( 111NkyaNkyakyakyNN)1(1110MkxbMkxbkxbkxbMM)()(1)()()()1(111)1(1110zazbzazazazbzbzbbzXzYzHNNNNMMMM, 111NNaaaa,110MMbbbbb精品资

4、料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页，共 12 页 - - - - - - - - - - h,k = impz(b, a, n)：计算 n 点单位脉冲响应hk ；也可简写为： h = impz(b, a, n)。impz(b, a)：绘制单位脉冲响应hk 的图形。离散系统响应yk 的计算 y = filter(b, a, x)：计算系统在输入x 作用下的零状态响应yk ；y = filter(b, a, x, zi)：计算系统在输入x 和初始状态作用下的完全响应yk 。zi 是由系统的初始状

5、态经过filtic函数转换而得到的初始条件：zi= filtic (b, a, Y0) ，Y0为系统的初始状态，Y0= y-1，y-2，y-3，.。2离散系统的系统函数零极点分析离散 LTI 系统的系统函数H(z) 可以表示为零极点形式：使用 MATLAB 提供的 roots函数计算离散系统的零极点；使用 zplane 函数绘制离散系统的零极点分布图。注意：在利用这些函数时，要求H(z) 的分子多项式和分母多项式的系数的个数相等，若不等则需要补零。3离散系统的频率响应当离散因果LTI 系统的系统函数H(z)的极点全部位于z 平面单位圆内时，系统的频率响应可由H(z)求出，即H, w=freqz

6、(b, a, n)：, 111NNaaaa,110MMbbbbb, 111NNaaaa,110MMbbbbb).()().()(1)()()(2121)1(111)1(1110NMNNNNMMMMpzpzpzzzzzzzkzazazazbzbzbbzXzYzH)(jjeje)e()()e(jHzHHz精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 3 页，共 12 页 - - - - - - - - - - 计算系统的n 点频率响应H，w为频率点向量。H=freqz(b, a, w)：计算系统在指定频率点

7、向量w上的频响；freqz(b,a)： 绘制频率响应曲线。其中： b 和 a 分别为系统函数H(z) 的分子分母系数矩阵；4利用 DTFT和 DFT确定离散系统的特性在很多情况下， 需要根据 LTI 系统的输入和输出对系统进行辨识，即通过测量系统在已知输入xk 激励下的响应yk 来确定系统的特性。若系统的脉冲响应为hk ，由于存在yk=xk*hk ，因而可在时域通过解卷积方法求解hk 。在实际应用中， 进行信号解卷积比较困难。因此，通常从频域来分析系统，这样就可以将时域的卷积转变为频域的乘积，从而通过分析系统输入序列和输出序列的频率特性来确定系统的频率特性H(je) ，再由 H(je) 得到系

8、统的脉冲响应hk 。若该 LTI 系统输入xk 的 DTFT为 X(je) ,系统输出yk 的 DTFT为 Y(je) ，则系统的频率特性可表示为：有限长序列的DTFT可以利用 FFT计算出其在区间内的N个等间隔频率点上的样点值： X = fft(x,N); Y = fft(y,N);再利用 H = Y./X 和 h = ifft(H,N) 可以得到系统的单位脉冲响应hk3实验设备及材料计算机一台及MATLAB 仿真软件。4实验方法步骤及注意事项实验方法步骤：先打开电脑，然后再打开MATLAB 仿真软件，在Blank M-File中输入程序，然后再编译运行程序，直到程序能编译运行为止。)e()

9、e()e(jjjXYH精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 4 页，共 12 页 - - - - - - - - - - 注意事项：(1). 为了省时间以及编译的方便性，程序应该在Blank M-File中输入， 而不应该在Command Window中直接运行；(2). 在使用 MATLAB 时应注意中英输入法的切换，在中文输入法输入程序时得到的程序是错误的；(3). MATLAB 中两个信号相乘表示为x.*u,中间有个 . ，同样两个信号相除也是如此，也就是在实验中要注意乘和点乘的区别。二实

10、验内容1实验现象与结果1. 已知某 LTI 系统的差分方程为：（1）初始状态，输入计算系统的完全响应。（2）当以下三个信号分别通过系统时，分别计算离散系统的零状态响应：（3）该系统具有什么特性（1）a=1,;b=,;N=100;x=ones(1,N);zi=filtic(b,a,1,2);y=filter(b,a,x,zi)stem(y);20675.011349.00675.02412.0 1143.1kxkxkxkykyky22, 11yykukx)107cos(;)5cos(;)10cos(321kukkxkukkxkukkx精品资料 - - - 欢迎下载 - - - - - - - -

11、 - - - 欢迎下载 名师归纳 - - - - - - - - - -第 5 页，共 12 页 - - - - - - - - - - （2）a=1,;b=,;N=100;k=1:N;x1=cos(pi/10*k);y1=filter(b,a,x1)stem(y1);精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6 页，共 12 页 - - - - - - - - - - x2=cos(pi/5*k);y2=filter(b,a,x2)stem(y2);精品资料 - - - 欢迎下载 - - - -

12、- - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 7 页，共 12 页 - - - - - - - - - - x3=cos(7*pi/10*k);y3=filter(b,a,x3)stem(y3);4. 已知某离散系统的输入输出序列。输入序列： 2， ， ， ， ， ， ， ， ， ， ， ， ， ， ，后面的数值均趋于0；输出序列：， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ，（1）绘出输入输出信号的波形。（2）计算该系统的频率响应，并绘出其幅频特性。（3）计算该系统的单位脉冲响应，

13、并绘出其波形。(1)x=2,;y=,;subplot(2,1,1);)(jeHkh精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页，共 12 页 - - - - - - - - - - stem(x);subplot(2,1,2);stem(y); (2)x=2,;y=,;X=fft(x,256);Y=fft(y,256);H=Y./X;H0=abs(H);plot(H0);精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - -

14、 - -第 9 页，共 12 页 - - - - - - - - - - (3)h=ifft(H);plot(h0);2对实验现象、实验结果的分析及其结论思考题：系统函数的零极点对系统频率特性有何影响答：系统零点越接近1，则幅频响应的波谷越低，即滤波器的滤波衰减最低点越低；系统极点越接近 1，则幅频响应的波峰越高，即滤波器的滤波通带最高点越高；极点主要影响频率特性的峰值，极点愈靠近单位圆，峰值愈尖锐；零点主要影响频率特性的谷值，零点愈靠近单位圆，谷值愈深，当零点在单位圆上时，频率特性为零， 一个传递函数有几个极点幅度响应就有几个峰值，对应出现一些谷值。结论：当离散 LTI 系统的系统函数H(Z

15、)的 ROC 包含单位圆时，系统的频率响应是系统函数在单位圆上的 Z变换， 系统函数与系统的稳定性表现在对于离散LTI 系统，系统稳定的充分必要条件是H(Z)的 ROC 包含 Z 平面上的单位圆。若离散LTI 系精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 10 页，共 12 页 - - - - - - - - - - 统是因果系统，由于hk 是因果序列，所以稳定系统的ROC 必须包含单位圆及单位圆外 Z 平面的全部区域，由于H(Z)在其 ROC 中不能有极点，故因果的离散LTI 系统稳定的充要条件是H

16、(Z)的极点全部位于 Z 平面的单位圆内。总结：通过本次试验可以得出需要注意的是，离散时间系统的系统函数可能有两种形式，一种是分子和分母多项式均按Z 的正幂次项排列， 另一种是分子和分母多项式均按Z 的负幂次项排列。 无论采用哪种排列形式，分子和分母多项式系数向量的维数一定要相同，不同的要用 0 补齐，否则Z=0的零点或极点就可能被漏掉。系统函数H(Z) 的零极点分布完全可以决定系统的特性，对离散系统特性的分析具有非常重要的意义。通过对系统函数零极点的分析可以获得离散系统以下几个方面的特征：1. 系统脉冲响应h（n）的时域特性2. 离散系统的稳定性3. 离散系统的频率特性（幅频响应和相频响应）MATLAB 为我们提供了求解这些特性及绘制系统函数零极点图的相关函数，如 zplane 、impz、freqz等教师评语及评分：精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 11 页，共 12 页 - - - - - - - - - - 签名：年月日精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 12 页，共 12 页 - - - - - - - - - -