实验四 离散系统分析 一、 实验目的 深刻理解离散时间系统的系统函数在分析离散系统的时域特性、 频域特性以 及稳定性中的重要作用及意义，熟练掌握利用MATLA分析离散系统的时域响应、 频响特性和零极点的方法。掌握利用 DTFT和DFT确定系统特性的原理和方法。 二、 实验原理 可以在时域、复频域（Z域）及频域分析系统，在以上三种域表征系统固有特性 的量分别为： 单位冲激响应h（n）（时域表征）； 系统函数H（z） （ Z域表征）； 频率响应H（ep）（频域表征）。 MATLA吐要从以上三方面提供了许多可用于分析线性时不变系统的函数，包含 系统时域响应、系统函数、系统频域响应等分析函数。 本实验通过调用各种系统预置函数来求系统的以上几个表征量以及零极点图。 三、 实验内容 1. 已知某LTI系统的差分方程为： yk -1.143yk -1 0.412yk -2 =0.0675xk 0.1349xk -1 0.0675x(2 (1)初始状态 y-1=1, y-2=2 ，输入 xk uk，计算系统的全响应。 程序段： N=40; b=0.0675,0.1349,0.0675; a=1,-1.143,0.412; x=on es(1,N); zi=filtic(b,a,1,2); y=filter(b,a,x,zi); stem(y) xlabel(k); title(yk); 结果: y.k 二 二 7 二 xdk =cos( k)uk; X2【k =cosUk)uk;X3k =cos( k)uk 10 5 10 程序 N=30; k=O:N; b=0.0675,0.1394,0.0675; a=1,-1.143,0.412; x1=cos(pi*0.1.*k); x2=cos(pi*0.2*k); x3=cos(pi*0.7*k); y1=filter(b,a,x1); y2=filter(b,a,x2); y3=filter(b,a,x3); subplot(3,1,1); stem(y1) subplot(3,1,2); stem(y2) subplot(3,1,3); stem(y3): 结果： (3)该系统具有什么特性? 答：因果稳定。 2. 已知某因果LTI系统的系统函数为: 1 _2 _3 _4 0.03571 0.1428Z 0.2143z 0.1428z 0.03571Z 1 -1.035Z0.8264Z，-0.2605Z； 0.04033Z* (1)计算系统的单位冲激响应 程序: N=50; k=1:N; b=0.0357,0.1428,0.2143,0.1428,0.0357; a=1,-1.035,0.8264,-0.2605,0.04033; y1=impz(b,a,N); stem(yl) 结果为： H(z)二 (2)当信号xkHuk cos(k)uk cos(2k)uk通过系统时，计算系统的零状 态响应。 程序为： N=50; k=1:N; b=0.0357,0.1428,0.2143,0.1428,0.0357; a=1,-1.035,0.8264,-0.2605,0.04033; x=o nes(1,N)+cos(pi*0.25*k)+cos(0.5*pi*k); y1=impz(b,a,N); y2=filter(b,a,x); subplot(2,1,1); stem(y1) subplot(2,1,2); stem(y2) 结果为：O.B 3. 已知LTI系统的输入输出序列分别为 xk二uk cos-k)uk cos(? k)uk (a) xk =(*)kuk, yk =*(1)kuk讨)kuk (b) xk =(kuk, yk =(4)kuk -(-l)kluk -1 (1)利用解析方法分别求解系统的单位取样响应。 程序： N=20; a仁1,-0.75; b仁1.25,-19/16; h1=impz(b1,a1,N); subplot(2,1,1); stem(h1) b2=1,-1; a2=1; h2=impz(b2,a2,N); subplot(2,1,2); stem(h2) 结果为：0 4 0 5 1Q 15 20 25 3D 40 45 5D a 6 (2)利用系统辨识原理确定 并求出系统的单位脉冲响应hk。比较解析方法 与系统辨识方法得到的系统单位冲激响应，分析误差原因。 程序： N=50; k=0:N; x1=0.5.Ak; y1=0.25*0.5.Ak+0.25.Ak; X1=fft(x1); 丫仁fft(y1); H仁 Y1/X1; h1=ifft(H1) n=1:50; X2=0.25.A n; y2=0.25.A n-4*0.254 n; X2=fft(x2); Y2=fft(y2); H2=Y2/X2; h2=ifft(H2) 结果： hi 二 H(ej ) 1.OF1 十 0,00001 h2 二 -3.0000 - 0.000015o iE I-； 10 I 12 lb I.； 10 误差原因：取样点数不够多引起的误差，而且不易与 k对应 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 , -0.7564 , 5.4872 , -11.7557 , 19.7533 , -29.6298 , 41.4666 , -55.2433 , 70.7979, -87.7810 (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,zeros(1,16); 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; subplot(2,1,1); stem(x) subplot(2,1,2); stem(y) 结果： ia0|- SO- f 50 - n 存匸一rE. 世预骨性 plot(k,x); title( 处理前); X=fft(x,150); h=ifft(H); H=fft(H,150); Y=X.*H; y=ifft(Y); subplot(1,2,2); plot(1：N),y); title(处理后); 结果： 实验思考题 1. 系统函数的零极点对系统频率特性有何影响？ ANS系统零点越接近1,贝肿畐频响应的波谷越低，即滤波器的滤波衰减最低 点越低；系统极点越接近1,贝肿畐频响应的波峰越高，即滤波器的滤波通带最高 点越高。极点主要影响频率响应的峰值，极点愈靠近单位圆，峰值愈尖锐；零点 主要影响频率特性的谷值，零点愈靠近单位圆，谷值愈深，当零点在单位圆上时， 频率特性为零，一个传递函数有几个极点幅度响应就有几个峰值， 对应出现一些 谷值。 2. 对于因果稳定实系数的低通、高通、带通、带阻数字滤波器，零极点分布有何 特点？ ANS:因为是因果稳定系统，所以极点都在单位圆内。若为最小相位系统，其 零点也在单位圆内。 3. 离散系统的系统函数的零极点对系统脉冲响应有何影响？ ANS系统函数的极点位置决定序列包络的变化趋势和变化频率， 极点的半径 决定了序列包络的变化趋势，而极点的幅角决定序列包络的变化频率， 而零点位 置只影响冲激响应的幅度大小和相位。 4. 若某因果系统不稳定，有哪些主要措施可使之稳定？ ANS改变参数，让极点在单位圆内。 5. 从频域利用DFT确定离散LTI系统的特性，一般会产生哪些误差，如何改善？ ANS:频谱混叠，对于带限连续信号，只要提高抽样频率使之满足时域抽样 定理；对于非带限信号，更具实际情况对其进行低通滤波，使之成为带限信号； 频谱泄露，时域加窗使之成为有限长序列；栅栏现象，在序列后补零，构成新序 列后再求频谱 6. 若使用DFT对连续LTI系统进行辨识，需要解决哪些问题？ ANS循环卷积长度的确定。 单纯的课本内容，并不能满足学生的需要，通过补充，达到内容的完善 教育之通病是教用脑的人不用手，不教用手的人用脑，所以一无所能。教育革命的对策是手脑联盟，结果是手与脑的力量都可以大到不可思议。