《实验二、离散系统分析》由会员分享,可在线阅读,更多相关《实验二、离散系统分析(9页珍藏版)》请在金锄头文库上搜索。
1、本科学生验证性实验报告学号 104090459 姓名 张静 学院 物电学院 专业、班级 10 电子实验课程名称 数字信号处理实验 教师及职称 杨卫平 教授 开课学期 2013 至 2013 学年 下 学期填报时间 2013 年 5 月 23 日云南师范大学教务处编印实验序号 2 实验名称 离散系统分析实验时间 2013.5.21 实验室 同析 3-313一实验预习1实验目的深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义,熟练掌握利用 MATLAB 分析离散系统的时域响应、频响特性和零极点的方法。掌握利用 DTFT 和 DFT 确定系统特性的原理和方
2、法。2实验原理、实验流程或装置示意图MATLAB 提供了许多可用于分析线性时不变连续系统的函数,主要包含有系统函数、系统时域响应、系统频域响应等分析函数。1. 离散系统的时域响应离散时间 LTI 系统可用如下的线性常系数差分方程来描述:已知输入信号 xk以及系统初始状态 y-1,y-2 .,就可以求出系统的响应。MATLAB 提供了计算差分方程的数值解的函数,可以计算上述 差 分方程描述的离散系统的单位脉冲响应、零状态响应和完全响应等。 在调用 MATLAB 函数时,需要利用描述该离散系统的系数函数。对差分方程进行 Z变换即可得系统函数:在 MATLAB 中可使用向量 a 和向量 b 分别保存
3、分母多项式和分子多项式的系数:这些系数均从 z0 按 z 的降幂排列。 离散系统的单位冲激响应 hk的计算 h,k = impz(b, a):计算系统的单位脉冲响应 hk和相应的时间向量 k;也可简写为:h = impz(b, a)。其中: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 和初始
4、状态作用下的完全响应 yk。zi 是由系统的初始状态经过 filtic 函数转换而得到的初始条件:zi= filtic (b, a, Y0) ,Y0 为系统的初始状态,Y0= y-1,y-2,y-3,.。)1(11 NkyakyayakNL0 MxbxbbM )(1)( )1(0 zabzzbzXYHNNMML,11Na ,110L,11NaaL,110MbbL,11NaaL ,110MbbL2离散系统的系统函数零极点分析离散 LTI 系统的系统函数 H(z)可以表示为零极点形式:使用 MATLAB 提供的 roots 函数计算离散系统的零极点;使用 zplane 函数绘制离散系统的零极点分布
5、图。注意:在利用这些函数时,要求 H(z)的分子多项式和分母多项式的系数的个数相等,若不等则需要补零。3离散系统的频率响应当离散因果 LTI 系统的系统函数 H(z)的极点全部位于 z 平面单位圆内时,系统的频率响应可由 H(z)求出,即H, w=freqz(b, a, n): 计算系统的 n 点频率响应 H,w 为频率点向量。H=freqz(b, a, w) :计算系统在指定频率点向量 w 上的频响;freqz(b,a): 绘制频率响应曲线。其中:b 和 a 分别为系统函数 H(z)的分子分母系数矩阵;4利用 DTFT 和 DFT 确定离散系统的特性 在很多情况下,需要根据 LTI 系统的输
6、入和输出对系统进行辨识,即通过测量系统在已知输入 xk激励下的响应 yk来确定系统的特性。若系统的脉冲响应为 hk,由于存在 yk=xk* hk,因而可在时域通过解卷积方法求解hk。在实际应用中,进行信号解卷积比较困难。因此,通常从频域来分析系统,这样就可以将时域的卷积转变为频域的乘积,从而通过分析系统输入序列和输出序列的频率特性来确定系统的频率特性 H( ) ,再由 H( ) 得到系统的脉冲响应 hk。jeje若该 LTI 系统输入 xk的 DTFT 为 X( ) ,系统输出 yk的 DTFT 为 Y( ) ,则系统j je的频率特性可表示为:有限长序列的 DTFT 可以利用 FFT 计算出
7、其在区间内的 N 个等间隔频率点上的样点值:X = fft(x,N); Y = fft(y,N);再利用 H = Y./X 和 h = ifft(H,N) 可以得到系统的单位脉冲响应 hk3实验设备及材料计算机一台及 MATLAB 仿真软件 。4实验方法步骤及注意事项实验方法步骤:先打开电脑,然后再打开 MATLAB 仿真软件,在 Blank M-File 中输入程序,然后再编译运行程序,直到程序能编译运行为止。).(1)( 21)1(0 NMNNMM pzpzkazabbzXYH L)(jjej e)(ej z)e()(jj注意事项:(1).为了省时间以及编译的方便性,程序应该在 Blank
8、 M-File 中输入,而不应该在Command Window 中直接运行;(2).在使用 MATLAB 时应注意中英输入法的切换,在中文输入法输入程序时得到的程序是错误的;(3). MATLAB 中两个信号相乘表示为 x.*u,中间有个. ,同样两个信号相除也是如此,也就是在实验中要注意乘和点乘的区别。二实验内容1实验现象与结果1.已知某 LTI 系统的差分方程为:(1)初始状态 ,输入 计算系统的完全响应。(2)当以下三个信号分别通过系统时,分别计算离散系统的零状态响应:(3)该系统具有什么特性?(1)a=1,-1.143,0.412;b=0.0675,0.1349,0.0675;N=10
9、0;x=ones(1,N);zi=filtic(b,a,1,2);y=filter(b,a,x,zi)stem(y);(2)a=1,-1.143,0.412;b=0.0675,0.1349,0.0675;N=100;k=1:N;x1=cos(pi/10*k);y1=filter(b,a,x1)20675.1349.0675. 21kxkxkxyy2,1yu)7cos(;)cos(;)0cos( 321 ukxukx stem(y1);x2=cos(pi/5*k);y2=filter(b,a,x2)stem(y2);x3=cos(7*pi/10*k);y3=filter(b,a,x3)stem(
10、y3);4.已知某离散系统的输入输出序列。输入序列:2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0.0083,0.0041,0.002,0.001,0.0005,0.0002,0.0001,0.0001,后面的数值均趋于 0;输出序列:0.0056,-0.0259, 0.073,-0.1593 ,0.297,-0.4974,0.7711,-1.1267,1.5702,-2.1037,2.724,-3.4207,4.174,-4.9528,5.7117,-6.3889,6.9034 ,-7.1528 ,7.012,-6.3322,4.9416,-2.648
11、,-0.7564,5.4872,-11.7557,19.7533 ,-29.6298 ,41.4666, -55.2433,70.7979,-87.7810(1)绘出输入输出信号的波形。(2)计算该系统的频率响应 ,并绘出其幅频特性。(3)计算该系统的单位脉冲响应 ,并绘出其波形。(1)x=2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0.0083,0.0041,0.002,0.001,0.0005,0.0002,0.0001,0.0001;y=0.0056,-0.0259,0.073,-0.1593,0.297,-0.4974,0.7711,-1.126
12、7,1.5702,-2.1037,2.724,-3.4207,4.174,-4.9528,5.7117,-6.3889,6.9034,-7.1528,7.012,-6.3322,4.9416,-2.648,-0.7564,5.4872,-11.7557,19.7533,-29.6298,41.4666,-55.2433,70.7979,-87.7810;subplot(2,1,1);stem(x);subplot(2,1,2);)(jeHkhstem(y); (2)x=2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0.0083,0.0041,0.002,0
13、.001,0.0005,0.0002,0.0001,0.0001;y=0.0056,-0.0259,0.073,-0.1593,0.297,-0.4974,0.7711,-1.1267,1.5702,-2.1037,2.724,-3.4207,4.174,-4.9528,5.7117,-6.3889,6.9034,-7.1528,7.012,-6.3322,4.9416,-2.648,-0.7564,5.4872,-11.7557,19.7533,-29.6298,41.4666,-55.2433,70.7979,-87.7810;X=fft(x,256);Y=fft(y,256);H=Y./
14、X;H0=abs(H);plot(H0);(3)h=ifft(H);plot(h0);2对实验现象、实验结果的分析及其结论思考题:系统函数的零极点对系统频率特性有何影响?答:系统零点越接近 1,则幅频响应的波谷越低,即滤波器的滤波衰减最低点越低; 系统极点越接近 1,则幅频响应的波峰越高,即滤波器的滤波通带最高点越高;极点主要影响频率特性的峰值,极点愈靠近单位圆,峰值愈尖锐;零点主要影响频率特性的谷值,零点愈靠近单位圆,谷值愈深,当零点在单位圆上时,频率特性为零,一个传递函数有几个极点幅度响应就有几个峰值,对应出现一些谷值。结论:当离散 LTI 系统的系统函数 H(Z)的 ROC 包含单位圆时
15、,系统的频率响应是系统函数在单位圆上的 Z 变换,系统函数与系统的稳定性表现在对于离散 LTI系统,系统稳定的充分必要条件是 H(Z)的 ROC 包含 Z 平面上的单位圆。若离散 LTI 系统是因果系统,由于 hk是因果序列,所以稳定系统的 ROC 必须包含单位圆及单位圆外 Z 平面的全部区域,由于 H(Z)在其 ROC 中不能有极点,故因果的离散 LTI 系统稳定的充要条件是 H(Z)的极点全部位于 Z 平面的单位圆内。总结:通过本次试验可以得出需要注意的是,离散时间系统的系统函数可能有两种形式,一种是分子和分母多项式均按 Z 的正幂次项排列,另一种是分子和分母多项式均按 Z 的负幂次项排列。无论采用哪种排列形式,分子和分母多项式系数向量的维数一定要相同,不同的要用 0 补齐,否则 Z=0 的零点或极点就可能被漏掉。系统函数 H(Z)的零极点分布完全可以决定系统的特性,对离散系统特性的分析具有非常重要的意义。通过对系统函数零极点的分析可以获得离散系统以下几个