用MATLAB进行系统分析.doc

实验三 用MATLAB进行系统分析1、 求系统的零极点 b=2 3; a=1 0.4 1; b,a=eqtflength(b,a); z,p,k=tf2zp(b,a)z = 0 -1.5000p = -0.2000 + 0.9798i -0.2000 - 0.9798ik = 2写成零极点形式： 2、 系统零极点图的绘制 b=2 3; a=1 0.4 1; b,a=eqtflength(b,a); z,p,k=tf2zp(b,a)z = 0 -1.5000p = -0.2000 + 0.9798i -0.2000 - 0.9798ik = 2 zplane(z,p);3、 有理传输函数的部分分式展开residuez例1 b=0 -3; a=2 -5 2; r,p,k=residuez(b,a)r = -1 1p = 2.0000 0.5000k = 部分分式展开为：例2 b=18; a=18 3 -4 -1; r,p,k=residuez(b,a)r = 0.3600 0.2400 0.4000p = 0.5000 -0.3333 -0.3333k = 部分分式展开为：4、 数字滤波器的单位脉冲响应impz b=0.2 0.1 0.3 0.1 0.2; a=1 -1.1 1.5 -0.7 0.3; impz(b,a,50)5、 求解差分方程 N=41; b=0.8 -0.44 0.36 0.22; a=1 0.7 -0.45 -0.6; x=1 zeros(1,N-1); k=0:1:N-1; y=filter(b,a,x); stem(k,y) xlabel(n); ylabel(幅度);练习题：1、 用MATLAB将函数展开成零点极点乘积的形式。写出系统函数的零极点表示式。 b=1 0 2; a=1 -1 -2; b,a=eqtflength(b,a); z,p,k=tf2zp(b,a)z = 0 + 1.4142i 0 - 1.4142ip = 2 -1k = 1零极点表示式：2、 画出系统的零极点图。 b=1 0 2; a=1 -1 -2; b,a=eqtflength(b,a); z,p,k=tf2zp(b,a); zplane(z,p);3、 将系统函数展开成部分分式的形式。写出展开的结果。 b=1 0 2; a=1 -1 -2; r,p,k=residuez(b,a)r = 1 1p = 2 -1k =-1部分分式的结果：4、 求出系统的单位脉冲响应。（0=n b=1 0 2; a=1 -1 -2; impz(b,a,5)5、 若系统的输入激励是,求此时系统的输出。其中。 N=5; b=1 0 2; a=1 -1 -2; x=ones(1,N); k=0:1:N-1; y=filter(b,a,x); stem(k,y) xlabel(n); ylabel(幅度);6、 若系统的初始条件为y(-1)=2 y(-2)=-0.5,求系统全响应。 N=5; b=1 0 2; a=1 -1 -2; x=ones(1,N); y=2,-0.5; z=filtic(