数字信号处理实验一离散系统的时域分析.docVIP

实验一离散系统的时域分析 实验目的: （1） 熟悉并掌握离散系统的差分方程表示法； （2） 加深对冲激响应和卷积分析方法的理解。 二、实验原理与方法 离散系统 x[n\ Discrete-time systme 其输入、输出关系可用以下差分方程描述： /V M y[n-k] = ^pkx[n-k] 々=() k=G 输入信号分解力冲激信号， OO x[n] = x[m]8[n - m] o 州=一OO 记系统单位冲激响应 S[n] h[n], 则系统响应为如下的卷积计算式： OO y[n] = x[n]^h[n] = ^x[m]h[n- m] 当 =0，々 = l，2，..JV时，h[nl是有限长度的（n: [0, Mp,称系统为FIR系统; 反之，称系统为IIR系统。 三、实验内容及要求 1、 要求分别用filter、conv、impz三种函数完成。 （1） 用conV（x，y）函数计算序列｛-2 0 1-1 3｝和序列｛1 2 0 -1｝的离散卷积 程序如下： xn=[-2 01-13]; hn=[l 20-1]; yn=conv(xn,hn); 结果如下: 1 -3 1 -3 CD File Edit Debug Desktop Window Help yn -2 -4 A? (2)用yn=filter(B, A, xn)函数计算差分方程 沖)+0.7Xi-Q -⑽沖-2) - 0.6X*-X =aSz(R) - 0.44冲-1) *0.36 场-2)+102jc(n-3) 当输入序列为《(|| =逐(#时的输出结果，C3x340。 程序如下： xn=[l，zeros(l，40)]; B=[l 0.7 -0.45 -0.6];A=[0.8 -0.44 0.36 0.22]; yn=filter(A,B,xn) n=O:length(yn)-l; subplot( 1,1,1 );stem(n,yn,7) title((a),);xlable(,n,);ylable(y(n)) 结果如下： File Edit Debug Desktop Window hdelp yn = Columns 1 through 5 0.8000 -1.0000 1. 4200 -0.9440 0. 6998 Columns 6 through 10 -0.0627 -0.2076 0.5370 -0.5069 0.4719 Colujnns 11 through 15 -0. 2363 0.0736 0. 1253 -0.1964 0. 2380 Columns 16 through 20 -0.1798 0. 1151 -0.0187 -0. 0430 0. 0908 Columns 21 through 25 -0. 0941 0.0809 -0. 0445 0. 0111 0. 0207 Columns 26 through 30 -0.0362 0.0414 -0.0328 0.0198 -0. 0038 Columns 31 through 35 -0.0081 0.0158 - 0.0170 0. 0142 -0. 0081 Coliunns 36 through 40 0.0018 0.0036 -0.0065 0.0073 -0. 0059 Column 41 0.0035 A ? A Command Window 口 回 图如下: (3)用impz(B,A,N)绘制出差分方程 ^?)4-0.7X-Q-QL45Xk-2)-0.6X*-X = Q.Sz(ji)-0.44j(M-D4-O.Kz(?-2)4-102je(?-3) 所对应的单位脉冲响应/?(/!), 0 zt 19 o 程序如下： N=0:19; A=[l 0.7 -0.45 -0. 6];B=[0.8 -0.44 0.36 0.22]; hn=impz (B, A, N) n=0:length(hn) -1; subplot (1, 1, 1) ;stem(n, hn, ’. ’) title(’⑹’);xlabel (’ n’)；ylabel (’ y(n)’) 结果如下： 人 Command Window File Edit Debug Desktop Window Help 0.8000 -1.0000 1. 4200 -0.9440 0.6998 -0.0627 -0.2076 0.5370 -0.5069 0.4719 -0.2363 0.0736 0. 1253 -0.1964 0.2380 -0.1798 0. 1151 -0.0187 -0.0430 0.0908 图如下: 四、实验报告要求 简述实验原理，画出程序框图，列出实验程序淸