用MATLAB进行系统分析.docVIP

实验三 用MATLAB进行系统分析 求系统的零极点 b=[2 3]; a=[1 0.4 1]; [b,a]=eqtflength(b,a); [z,p,k]=tf2zp(b,a) z = 0 -1.5000 p = -0.2000 + 0.9798i -0.2000 - 0.9798i k = 2 写成零极点形式： 系统零极点图的绘制 b=[2 3]; a=[1 0.4 1]; [b,a]=eqtflength(b,a); [z,p,k]=tf2zp(b,a) z = 0 -1.5000 p = -0.2000 + 0.9798i -0.2000 - 0.9798i k = 2 zplane(z,p); 有理传输函数的部分分式展开residuez 例1 b=[0 -3]; a=[2 -5 2]; [r,p,k]=residuez(b,a) r = -1 1 p = 2.0000 0.5000 k = [] 部分分式展开为： 例2 b=[18]; a=[18 3 -4 -1]; [r,p,k]=residuez(b,a) r = 0.3600 0.2400 0.4000 p = 0.5000 -0.3333 -0.3333 k = [] 部分分式展开为： 数字滤波器的单位脉冲响应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) 求解差分方程 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(幅度); 练习题： 用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.4142i p = 2 -1 k = 1 零极点表示式： 画出系统的零极点图。 b=[1 0 2]; a=[1 -1 -2]; [b,a]=eqtflength(b,a); [z,p,k]=tf2zp(b,a); zplane(z,p); 将系统函数展开成部分分式的形式。写出展开的结果。 b=[1 0 2]; a=[1 -1 -2]; [r,p,k]=residuez(b,a) r = 1 1 p = 2 -1 k = -1 部分分式的结果： 求出系统的单位脉冲响应。（0=n=4） b=[1 0 2]; a=[1 -1 -2]; impz(b,a,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(幅度); 若系统的初始条件为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(b,a,y); k=0:1:N-1; y=filter(b,a,x,z); stem(k,y) xlabel(n);ylabel(幅度)