实验4离散系统的分析报告.docVIP

. . 实验四 离散系统分析 实验目的 深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义.熟练掌握利用MATLAB分析离散系统的时域响应、频响特性和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法。 二、实验原理 可以在时域、复频域（Z域）及频域分析系统.在以上三种域表征系统固有特性的量分别为： ? 单位冲激响应 h(n)（时域表征）； ? 系统函数 H(z) （ Z域表征）； ? 频率响应 H(ejw)（频域表征）。 MATLAB主要从以上三方面提供了许多可用于分析线性时不变系统的函数.包含系统时域响应、系统函数、系统频域响应等分析函数。 本实验通过调用各种系统预置函数来求系统的以上几个表征量以及零极点图。 三、实验内容 1.已知某LTI系统的差分方程为： （1）初始状态 .输入 .计算系统的全响应。 程序段： N=40; b=[0.0675,0.1349,0.0675]; a=[1,-1.143,0.412]; x=ones(1,N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi); stem(y) xlabel(k); title(y[k]); 结果： 2）当以下三个信号分别通过系统时.分别计算离散系统 的零状态响应： 程序 N=30; k=0: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）计算系统的单位冲激响应。 程序： 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(y1) 结果为： （2）当信号 通过系统时.计算系统的零状态响应。 程序为： 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=ones(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) 结果为： 3.已知LTI系统的输入输出序列分别为 （1）利用解析方法分别求解系统的单位取样响应。 程序： N=20; a1=[1,-0.75]; b1=[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) 结果为： （2）利用系统辨识原理确定 并求出系统的单位脉冲响应h[k]。比较解析方法与系统辨识方法得到的系统单位冲激响应.分析误差原因。 程序： N=50; k=0:N; x1=0.5.^k; y1=0.25*0.5.^k+0.25.^k; X1=fft(x1); Y1=fft(y1); H1=Y1/X1; h1=ifft(H1) n=1:50; x2=0.25.^n; y2=0.25.^n-4*0.25.^n; X2=fft(x2); Y2=fft(y2); H2=Y2/X2; h2=ifft(H2) 结果： 误差原因：取样点数不够多引起的误差.而且不易与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