《信号与系统分析》课件第8章.pptVIP

图8.15例8.22系统经对角化变换后的模拟方框图(3)求系统函数通过以上分析可见，当系统具有不可控、不可观的状态变量时，则系统函数的零、极点会有相互抵消的现象，这表明仅用系统函数描述一个系统有时并不全面。系统函数只能描述系统的外部特性，故只能反映系统中可控和可观那部分的运动规律，而不能反映不可控和不可观那部分的运动规律。由图8.15可知，在本例中系统函数只描述了系统中既连接输入又连接了输出的那部分，所以整个系统的系统函数为H(s)=。由于状态方程描述了系统内部状态，所以用状态方程和输出方程来描述一个系统的运动更全面、更详尽。可以证明，如果一个系统既是可控制的,又是可观测的，则系统的内部描述和外部描述才是等价的。8.5系统的状态变量分析的MATLAB实现例8.23已知系统的状态方程和输出方程为试求状态转移矩阵eAt、系统函数矩阵H(s)和单位冲激响应矩阵h(t)。解求eAt要用到MATLAB符号工具箱中的expm函数,它的用法很简单,见下面的MATLAB程序。单位冲激响应矩阵可用公式h(t)=CeAtB+Dδ(t)实现,再用拉普拉斯反变换就可以求得系统函数矩阵H(s)［MATLAB程序］symst;A=［01;-2-3］;B=［10;11］;C=［10;11;02］;D=［00;10;01］;expm(t*A)ht=C*expm(t*A)*B+D*sym(′Dirac(t)′)Hs=laplace(ht)［程序运行结果］ans=［-exp(-2*t)+2*exp(-t),exp(-t)-exp(-2*t)］［-2*exp(-t)+2*exp(-2*t),2*exp(-2*t)-exp(-t)］ht=［-2*exp(-2*t)+3*exp(-t),exp(-t)-exp(-2*t)］［2*exp(-2*t)+Dirac(t),exp(-2*t)］［-6*exp(-t)+8*exp(-2*t),4*exp(-2*t)-2*exp(-t)+Dirac(t)］Hs=［-2/(s+2)+3/(1+s),1/(1+s)-1/(s+2)］［2/(s+2)+1,1/(s+2)］［-6/(1+s)+8/(s+2),4/(s+2)-2/(1+s)+1］例8.24已知连续时间系统的状态方程、输出方程、激励信号和系统的起始状态为求系统的零输入响应和零状态响应。解利用关系y(t)=CeAtλ(0-)+［CeAtB+Dδ(t)］x(t)来求解系统的完全响应。其中,第一部分为系统的零输入响应,用时域法求解即可,而表达式中的第二部分为系统的零状态响应,需要借助于拉普拉斯变换。［MATLAB程序］A=［-21;-30］;B=［0;1］;C=［11］;D=0;ramda0=［1;0］;yzir=C*expm(t*A)*ramda0yzsr=ilaplace(laplace(C*expm(t*A)*B+D*sym(′Dirac(t)′))*laplace(sym(′Heaviside(t)′)))［程序运行结果］yzir=exp(-t)*cos(t*2^(1/2))-2*2^(1/2)*exp(-t)*sin(t*2^(1/2))yzsr=-exp(-t)*cos(t*2^(1/2))+1例8.25给定离散时间系统状态方程、输出方程、激励信号和系统的初始条件分别为求系统输出y(n)。解计算Ak,系统函数H(z)和输出信号的Z变换Y(z)的公式分别如下式所示:其中,输出信号的Z变换的第一项对应零输入响应的Z变换,即零输入响应为yzir(n)=CAkλ(0)。在MATLAB环境中,利用符号法进行响应的分析时,零输入响应可以用时域法进行求解,也可以用Z变换法进行求解,但零状态响应就只能用Z变换法进行求解,否则要用到时域的卷积和,MATLAB的符号工具箱没有提供符号卷积函数。注意,在符号工具箱中,单位取样信号δ(n)使用charfcn［0］(n)表示,而ε(n)是用Heaviside(n)表示。而且要注意,不是字符串,而是符号变换来表示。