由传递函数转换成状态空间模型(1)

由传递函数转换成状况空间模子——办法多!!!SISO线性定常体系高阶微分方程化为状况空间表达式外部描写←—实现问题：有了内部构造—→模仿体系内部描写实现问题解决有多种办法,办法不合时成果不合.一、直接分化法因为对上式取拉氏反变换,则按下列纪律选择状况变量,于是有写成矩阵情势式中,把这种尺度型中的A系数阵称之为友阵.只要体系状况方程的系数阵A和输入阵b具有上式的情势,c 阵的情势可以随意率性,则称之为能控尺度型.则输出方程写成矩阵情势.在须要对现实体系进行数学模子转换时,不必进行盘算就可以便利地写出状况空间模子的A.b.c矩阵的所有元素.例：已知SISO体系的传递函数如下,试求体系的能控尺度型状况空间模子.解：直接得到体系进行能控尺度型的转换,即若选（若何斟酌？）斟酌式依次对第一式求导,并带入第二式;对第二式求导,并带入第三式;依次类推,便得到写成矩阵情势式中.只要体系状况空间表达式的A阵和c 阵具有上式的情势,b阵的情势可以随意率性,则称之为能不雅尺度型从情势上看,能控尺度型和能不雅尺度型的系数阵A是互为转置,能控尺度型输入阵b和能不雅尺度型输出阵c互为转置,这种互为转置的关系被称为对偶关系.将在第六章进一步评论辩论.经由过程以上对传递函数阵的能控尺度型或能不雅尺度型转换的评论辩论,对单输入体系而言,应留意如下问题：（1）传递函数转化成能控尺度型的状况空间表达式,状况方程的构造只由传递函数阵的顶点（特点）多项式肯定,而与其零点多项式无关,零点多项式只影响输出方程的构造.（2）从能不雅尺度型的转换可以看出,系数阵A的元素仅决议于传递函数顶点多项式系数,而其零点多项式则肯定输入阵B的元素.（3）只有当传递函数零点和顶点多项式同阶时,状况,.例：求前例的能不雅尺度型的状况空间模子解：直接得到能不雅尺度型的状况空间模子,即二、串联分化法若SISO体系的传递函数顶点互异,体系传递函数分子分母写成因式相乘情势图示！！三、并联分化法（对角尺度型/约旦尺度型——特点值尺度型）（一）若SISO体系的传递函数顶点互异,则可求得对角尺度型的模子.当体系的顶点互异时,体系传递函数分子分母写成因式相乘情势写成部分分式个中,其值为选择状况变量为（绘图示意状况变量的取法）即对上式拉氏反变换,得即写成矩阵情势式中,系数矩阵A为对角阵.对角线上的元素是传递函数G（s）的顶点,即体系的特点值.b阵是元素全为1的n×1矩阵.求对角尺度型模子的输出方程中c的构造对上式拉氏反变换,得假如体系的状况方程的A阵是对角阵,暗示体系的各个变量之间是解耦的.多变量的体系解耦是庞杂体系实现准确掌握的症结问题,关于若何实现解耦掌握将在第五章评论辩论.体系的状况构造图如图所示.例：设体系的闭环传递函数如下,试求体系对角尺度型的转换,得对角尺度型的转换为（二）对SISO体系式,当其有重特点值时,可以得到约当尺度型的状况空间模子.此时模子的系数矩阵A中与重特点值对应的那些子块都是与这些特点值相对应的约当块,即其重数为j,而其余为互异的特点值,则传递函数可以用部分分式睁开成式中,,其值为.绘图示意状况变量的取法：例：设体系的闭环传递函数如下,试求体系对约当准型的状况空间模子,该体系为四阶,有一个重顶点,重数为j=2,有两个互异的顶点,按部分分式睁开可得约当尺度型的模子为。

MIMO传递函数转化为状态空间模型Matlab代码1. 介绍MIMO（多输入多输出）系统是指系统具有多个输入和多个输出的特性。

在控制系统领域中，MIMO系统的建模和分析是非常重要的。

传递函数和状态空间模型是两种常用的系统建模方法。

本文将介绍如何将MIMO系统的传递函数转化为状态空间模型，并给出相应的Matlab代码实现。

2. MIMO系统的传递函数表示MIMO系统的传递函数通常表示为一个矩阵，每个元素对应一个输入到一个输出的传递函数。

假设有n个输入、m个输出，则MIMO系统的传递函数可以表示为一个m×n的传递函数矩阵G(s)。

传递函数矩阵的元素可以用s表示，如G11(s)、G12(s)等。

3. MIMO系统传递函数到状态空间模型的转化方法MIMO系统的传递函数可以通过状态空间模型来表示。

状态空间模型的基本形式如下：\[ \dot{x}(t) = Ax(t) + Bu(t) \]\[ y(t) = Cx(t) + Du(t) \]其中，A是状态矩阵，B是输入矩阵，C是输出矩阵，D是传递函数零极点对应的矩阵。

MIMO系统的传递函数可以通过以下步骤转化为状态空间模型：1）将传递函数矩阵分解为多个SISO（单输入单输出）系统的传递函数；2）针对每个SISO系统，可以将其转化为状态空间模型；3）将各个SISO系统的状态空间模型组合成一个整体的MIMO系统的状态空间模型。

4. Matlab代码实现下面我们通过一个实例来演示如何用Matlab将MIMO系统的传递函数转化为状态空间模型。

假设传递函数矩阵为：\[ G(s) = \begin{bmatrix} \frac{2s+1}{s^2+3s+2}\frac{3s+2}{s^2+4s+3} \\ \frac{4s+1}{s^2+2s+1}\frac{5s+2}{s^2+3s+2} \end{bmatrix} \]我们需要将传递函数矩阵分解为四个SISO系统的传递函数：\[ G11(s) = \frac{2s+1}{s^2+3s+2} \]\[ G12(s) = \frac{3s+2}{s^2+4s+3} \]\[ G21(s) = \frac{4s+1}{s^2+2s+1} \]\[ G22(s) = \frac{5s+2}{s^2+3s+2} \]针对每个SISO系统，我们可以将其转化为状态空间模型，以G11(s)为例：```Matlab将传递函数G11(s)转化为状态空间模型num = [2, 1]; 分子系数den = [1, 3, 2]; 分母系数[A11, B11, C11, D11] = tf2ss(num, den); 转化为状态空间模型```将各个SISO系统的状态空间模型组合成整体的MIMO系统的状态空间模型：```Matlab对四个SISO系统的状态空间模型进行组合A = [A11, A12; A21, A22];B = [B11, B12; B21, B22];C = [C11, C12; C21, C22];D = [D11, D12; D21, D22];```至此，我们成功地将MIMO系统的传递函数转化为状态空间模型，并通过Matlab代码实现了这一过程。

实验题目：传递函数到状态空间的实现 课程名称：计算机仿真 一、实验目的1、 理解并掌握传递函数转换为状态空间方程的方法2、 理解状态初值的计算方法二、 实验内容1、 应用MATLAB 编写一个可以实现传递函数到状态空间方程的可控可观规范型的ni 文件。

并用相应例题验证程序的止确性。

2、 完善该程序使具可以用來计算状态初值。

并用相应的例题验证程序 的正确性。

3、 程序中需要考虑分子分母同阶以及分母首系数不为1的两种情况。

三、 报告内容1、 给出m 文件的程序框图，及验证结果，并记录出现的错误，并给出 解决的方案。

若没有得到解决，请说清楚你的问题2、 如呆做了程序的状态初值得求解，请给岀相应的验证结杲，及程序 编写过程中出现的问题，若已经解决，给出貝体方法。

能观标准型为:2、计算状态变量初值：(1)不含u 的导数项时，则冇:A= • 0 0 •■1 0• •0 1■ ■… 0 ■…• • ••B=O' 0 ■ ■~a n~a n-l~a n-l…一如・丄Z?o s n +b 1s n "1+•••+d n ^1s+c n+…+01八一］s+a 八那么其状态空间模型能控标准型为:C=［（b n — bo （z n ） （&n _i — …@1 —加血）］ D=b n!1!实验理论传递函数为G(s)=1、 力能观=B 能观D 能观和0)X?(O) 1(0)」 yj(o)（2）系统微分方程不仅包含u 的输入项，而口包含u 的导数项，则:五、程序检验（1）输入一个分母首系数为1月.分子分母不同阶传递函数:2S 3+ 4S 2+ 3S + 5 G = -------------------------------S 4 + 2S 3 + 5S 2 + 4S + 2程序运行结果： 能控标准型:A 二0 1 0 00 1 00 0 0 1-2-4-5-2B =兀 1(0)a n-l an-2 …兀2（0）a n-2%一3…七(0) ■ • = an-3•■ • • • • • • ■^-1(0)■ 1 … _ 兀“(0)..10 (x)n xa x 1 y （o ）~Cn-l1 0 y （o ）一 Cn-2 • •… •■ y(0) ■ •+ _ Cn-3■ ■ ■…0 严)(0)_C]…0 严(()) ■ ■_ 0nxl /ix(n -1)一 C] w(O)〃(()):M(O)•• • ••• :宀(0)0 ]“"-2)(0)(/?-l)xly(0) y(0)5 342D 二能观标准型：A =0 00-21 00-40 10-50 01-2B =5342C =0 001D 二初值部分：请输入系统输出的初值二[1 ;1;1;1]请输入系统输入的初值二[0; 0; 0] x0 二12831运行结果正确（2）输入一个分母首系数为2 口分子分母同阶传递函数:S 2 + 2S + 3G =2S 2 + 5S + 3程序运行结果: 能控标准型:0. 5000初值部分：请输入系统输出的初值二[1；1] 请输入系统输入的初值二[0]xO 二A =0 -1. 5000 B =0 1 C 二1. 5000 D =0. 5000能观标准型:A 二0 1.0000 B =1. 5000 1.5000 C 二1.5000 D =1.0000 -2. 50001. 5000-1. 5000 -2. 50001. 50003. 50001.0000运行结果正确六.流程图七、实验小结通过木次实验我了解了如何通过matlab的编程来实现传递函数转化为状态空间方程的能控和能观性，并掌握了程序的状态初值的求解。

同步电机系统传递函数模型和状态空间模型同步电机系统是工业中非常常见的一种电机系统，广泛应用于风力发电、输电输能、轨道交通、工业生产等领域。

在电机运行过程中，了解同步电机系统传递函数模型和状态空间模型对于对其理解和进行系统控制非常重要。

同步电机系统的传递函数模型可以被表示为： G(s)= K / (T1s + 1)(T2s + 1)，其中K是增益，T1和T2是时间常数。

传递函数模型的意义是对输入和输出之间的关系进行建模。

在同步电机系统中，输入通常是电压，输出是转速。

如传递函数模型所示，输入的电压可以被分解成两个滞后时间常数的二阶系统，其中每个常数都代表了电机系统的动态转速响应。

这样的建模方法可以使我们更好地理解同步电机系统的转速响应，并从中得出控制策略。

除了传递函数模型，同步电机系统也可以由状态空间模型来描述，即将系统表示为一系列宏观的状态和它们之间的转换。

状态空间模型可以被表示为：x’=Ax+Bu, y=Cx+Du。

状态空间模型通过宏观的描述电机系统内部的变化，阐述了系统的动态响应方式和控制之间的关系。

在一个转速控制的例子中，输入电压被视为操作量u(t)，并且与电机系统的输出y(t)通过状态方程x’(t)=Ax(t)+Bu(t)和输出方程y(t)=Cx(t)+Du(t)进行联系。

这样的状态空间模型可以帮助我们设计和分析电机控制器。

总的来说，同步电机系统传递函数模型和状态空间模型都有各自的用途。

传递函数模型可以为我们提供关于电机系统动态特性的优越的定量分析方法，而状态空间模型则更适合于控制器的设计以及在运行实际系统时的仿真与测试。

或者，这两种模型也可以结合起来，并被用于对同步电机系统进行整体观察和控制。

在电机系统中，应用这样的模型可以使我们更好地理解它们的动态特性，并对其进行控制。

这样的导向是成为我们更有资格和信心地应对电机工业领域中的挑战和机遇的关键因素。

已知传递函数求状态空间表达式在控制系统理论中，常常需要将已知的传递函数转换为状态空间表达式。

这是因为状态空间形式更加直观，便于进行控制器设计和系统分析。

首先，我们需要将传递函数化简为标准形式：$$G(s) = frac{b_0 s^n + b_1 s^{n-1} + cdots + b_{n-1} s + b_n}{s^n + a_1 s^{n-1} + cdots + a_{n-1} s + a_n}$$其中 $n$ 为传递函数的阶数，$b_i$ 和 $a_i$ 是系数。

接下来，我们可以通过状态空间的基本方程来表示传递函数： $$begin{aligned}dot{x} &= Ax + Buy &= Cx + Duend{aligned}$$其中，$x$ 是 $n$ 维状态向量，$u$ 是 $m$ 维输入向量，$y$ 是$p$ 维输出向量。

$A$、$B$、$C$、$D$ 是系数矩阵，它们的维度分别为 $n times n$、$n times m$、$p times n$ 和 $p times m$。

我们可以通过下列步骤获得$A$、$B$、$C$ 和 $D$：1. 首先，将传递函数分解为零极点形式：$$G(s) =kfrac{(s-z_1)(s-z_2)cdots(s-z_n)}{(s-p_1)(s-p_2)cdots(s-p_n )}$$其中，$k$ 是比例系数，$z_i$ 和 $p_i$ 是零点和极点。

2. 利用零极点分解结果，构造传递函数的控制分式表达式：$$G(s) = kfrac{(s-z_1)}{(s-p_1)} cdot frac{(s-z_2)}{(s-p_2)} cdots frac{(s-z_n)}{(s-p_n)}$$3. 对每个控制分式，构造对应的状态空间模型：$$begin{aligned}dot{x_i} &= p_i x_i + uy_i &= z_i x_iend{aligned}$$其中，$i$ 取值为 $1$ 到 $n$。

实验题目：传递函数到状态空间的实现课程名称：计算机仿真一、实验目的1、理解并掌握传递函数转换为状态空间方程的方法2、理解状态初值的计算方法二、实验内容1、应用MATLAB®写一个可以实现传递函数到状态空间方程的可控可观规范型的m文件。

并用相应例题验证程序的正确性。

2、完善该程序使其可以用来计算状态初值。

并用相应的例题验证程序的正确性。

3、程序中需要考虑分子分母同阶以及分母首系数不为1的两种情况＜三、报告内容1、给出m文件的程序框图，及验证结果，并记录出现的错误，并给出解决的方案。

若没有得到解决，请说清楚你的问题2、如果做了程序的状态初值得求解，请给出相应的验证结果，及程序编写过程中出现的问题，若已经解决，给出具体方法。

四、实验理论1、传递函数为 --------------------那么其状态空间模型能控标准型为：A= B=C=能观标准型为：能观能控能观能控能观2、计算状态变量初值：D= 能控能观能控(1)不含u的导数项时，则有:y (o ) y (o )I = I]y (n 」(0)-五、程序检验程序运行结果: 能控标准型：A = 0 1 0 0 0 0 1 0 00 1-2 -4 -5 -2 B = 0、1(0)--a na n■…a1们 ■ y(0) ■ [ _ C n_ C n_2 —C 2 — Ci -u(0)-X 2(0)an/ an & (1)y(0)_皿 ......... —G 0u(0) X 3(0)a =a n(-aay(0)a+ _C n ( ....................................aaaU(0)■X n ±(0)a 1......... 0 y (n「0)a一 C jU2)(0) 川(0) _1 1 0 .................... 0 一yj (0)一I0 … 0屮2)(0)一 u 的输入项，而且包含u 的导数项，则: n 1n n n 1n (n-1)(n-1) 1x 1( 0)X 2(0)I : I I IX n ( 0)(2)系统微分方程不仅包含 (1)输入一个分母首系数为 1且分子分母不同阶传递函数:5 3 4 2D =能观标准型：A =0 0 0 -21 0 0 -40 1 0 -50 0 1 -2B =5342C =0 0 0 1D =初值部分：请输入系统输出的初值=[1；1；1；1]请输入系统输入的初值=[0；0；0]x0 =12831运行结果正确（2）输入一个分母首系数为2且分子分母同阶传递函数: 程序运行结果：能控标准型：A =0 1.0000-1.5000 -2.5000B =1C =1.5000 1.5000D =0.5000能观标准型：A =0 -1.50001.0000 -2.5000B =1.50001.5000C =1.5000 1.5000D =0.5000初值部分：请输入系统输出的初值=[1；1]请输入系统输入的初值=[0]x0 =3.50001.0000运行结果正确六、流程图七、实验小结通过本次实验我了解了如何通过matlab的编程来实现传递函数转化为状态空间方程的能控和能观性，并掌握了程序的状态初值的求解。

自动控制理论自动控制第二章周立芳徐正国连续时间控制系统的数学模型浙江大学控制科学与工程学系1第二章要点✓引言✓电路及组成✓线性代数与状态的基本概念✓传递函数及方块图✓机械传递系统✓其他的数学建模实例✓系统传递函数的计算✓非线性系统的线性化✓系统整体传递函数的确定✓仿真图✓信号流图从传函数到状间模的转换✓从传递函数到状态空间模型的转换2从传递函数到状态空间模型的转换◆从传递函数到并联状态图◆并联状态图◆A 矩阵的对角化◆利用状态变换求解状态方程◆状态方程的标准形式可控标准型◆◆可观标准型◆从方块图到状态空间模型控制科学与工程学系并联状态图由下面微分方程描述的SISO 系统可以由相应的传递函数表示并联状态图)()()( ;)()())(()(1210111i ii i i ni i n n n n n n s f s U s Z f s G s G c s s s c s c s c s c s G λλλλ-==+=---++++=∑=--并联状态图系统的状态转移信号流图如下图所示，图中省略了状态变量的初始值z i (t 0)。

Z 1(s)λ1f 1前馈通道Z 2(s)f 2U(s)Y(s)λ2:f n())()(1∑=+=ni i n s G c s G f λnc nZ n (s))()()( iii i i s f s U s Z s G λ-==图5.31 式(*) 的并联解耦仿真图（w=n ）并联状态图于是系统的状态空间模型为：所有元素均为1⎥⎤⎢⎡⎥⎤⎢⎡ 2111000λλnn +Λ=⎥⎥⎥⎢⎢⎢+⎥⎥⎥⎢⎢⎢=ub z u z z1000λw=n, d n ≠0, 否则d n =0[]ud u c f f f y n n n n +=+=⎦⎣⎦⎣z c z21A 是对角阵此时系统动态方程称为状态空间模型系统矩阵A 是对角阵，此时系统动态方程称为正则标准型状态空间模型，系统矩阵可表示为Λ(or A*)，相应的状态变量称为规范变量（canonical variables ）。

传递函数形式转换传递函数模型可转化成零极点模型，也可转换成状态空间模型。

例一、模型转换22235(56)52530()(1)(2)32s s s s G s s s s s -+-+==-+-+% name:model_transform_01clear,clc % 清除此前运行程序的变量，清屏num1=5*[1,-5,6] % 传函分子多项式 den1=conv(conv([1,-1],[1,-1]),[1,2]) % 传函分母多项式 tfGs=tf(num1,den1) % 传函的分式形式 zpGs=zpk(tfGs) % 转换成传函的零极点形式[z,p,k]=zpkdata(zpGs,'v') % 列出零极点和比例系数 ssGs=ss(tfGs) % 转换成状态空间形式运行结果：num1 = 5 -25 30den1 = 1 0 -3 2Transfer function:5 s^2 - 25 s + 30-----------------s^3 - 3 s + 2Zero/pole/gain:5 (s-3) (s-2)-------------(s+2) (s-1)^2z = 3.00002.0000p = -2.00001.00001.0000k = 5a = x1 x2 x3x1 0 0.75 -0.25x2 4 0 0x3 0 2 0b = u1x1 4x2 0x3 0c = x1 x2 x3y1 1.25 -1.563 0.9375d = u1y1 0Continuous-time model.状态方程：XaX bu =+ 11223300.750.25440000200xx x x u x x -⎛⎫⎡⎤⎛⎫⎛⎫ ⎪ ⎪ ⎪⎢⎥=+ ⎪ ⎪ ⎪⎢⎥ ⎪ ⎪ ⎪⎢⎥⎝⎭⎣⎦⎝⎭⎝⎭输出方程：T Y c X du =+[]1231.25 1.5630.93750x y x u x ⎛⎫ ⎪=-+⋅ ⎪ ⎪⎝⎭例二、传递函数模型的串联% series exampleclear,clcnum1=[0,1];den1=[1,1];num2=[0,2];den2=[1,2];num3=[0,3];den3=[1,3];Gs1=tf(num1,den1)Gs2=tf(num2,den2)Gs3=tf(num3,den3)Gs4=series(Gs1,Gs2) % 不允许三连串Gs4=series(Gs1,Gs2,Gs3) Gs5=series(Gs4,Gs3)Gs6=Gs1*Gs2*Gs3 % 允许三连乘运行结果：Transfer function ：1-----Transfer function:2-----s + 2Transfer function:3-----s + 3Transfer function:2-------------s^2 + 3 s + 2Transfer function:6----------------------s^3 + 6 s^2 + 11 s + 6Transfer function:6----------------------s^3 + 6 s^2 + 11 s + 6例三、传递函数模型的并联% parallel exampleclear,clcnum1=[0,1];den1=[1,1];num2=[0,2];den2=[1,2];num3=[0,3];den3=[1,3];Gs1=tf(num1,den1)Gs2=tf(num2,den2)Gs3=tf(num3,den3)Gs4= parallel (Gs1,Gs2) % 不允许三连并Gs4= parallel (Gs1,Gs2,Gs3) Gs5= parallel (Gs4,Gs3)Gs6=Gs1+Gs2+Gs3 % 允许三连加运行结果：Transfer function:1-----Transfer function:2-----s + 2Transfer function:3-----s + 3Transfer function:3 s + 4-------------s^2 + 3 s + 2Transfer function:6 s^2 + 22 s + 18----------------------s^3 + 6 s^2 + 11 s + 6Transfer function:6 s^2 + 22 s + 18----------------------s^3 + 6 s^2 + 11 s + 6例四、传递函数模型的反馈% feedback_example clear,clcnum1=[0,1];den1=[1,1];num2=[0,2];den2=[1,2]; num3=[0,3];den3=[1,3];Gs1=tf(num1,den1) Gs2=tf(num2,den2) Gs3=tf(num3,den3)Gs4= series (Gs1,Gs2)Gs5= feedback (Gs4,Gs3) % 建立反馈闭路的闭环传函 Gs6= feedback (Gs1*Gs2,Gs3) % 建立反馈闭路的闭环传函运行结果：Transfer function:1-----12123()1G G G s G G G =+Transfer function:2-----s + 2Transfer function:3-----s + 3Transfer function:2-------------s^2 + 3 s + 2Transfer function:2 s + 6-----------------------s^3 + 6 s^2 + 11 s + 12Transfer function:2 s + 6-----------------------s^3 + 6 s^2 + 11 s + 12例五、传递函数模型的单位反馈% unit_feedback_exampleclear,clcnum1=[0,1];den1=[1,1];num2=[0,2];den2=[1,2];Gs1=tf(num1,den1)Gs2=tf(num2,den2)[num3,den3]=series(num1,den1,num2,den2);[num4,den4]=cloop(num3,den3) % 求出单位反馈闭环传函的分子分母Gs4=tf(num4,den4) % 求出单位反馈闭环传函运行结果：Transfer function:1-----s + 1Transfer function:2-----s + 2num4 = 0 0 2 den4 = 1 3 4Transfer function:2-------------s^2 + 3 s + 4。

X n =_a n X ia n 4X 2-a 1X n u由传递函数转换成状态空间模型一一方法多！!!SISO 线性定常系统高阶微分方程化为状态空间表达式SISOy (n j+a y D+azy W )+…+a n y =b 0u 俨)+b 1u (m _L )+…+b m u(n ^m )b °s mb,sm 4s nys n」 a 2s n^ ■ a n外部描述＜实现问题：有了内部结构一-模拟系统 内部描述‘X = Ax +bu y =cx+ du实现冋题解决有多种方法，方法不同时结果不同直接分解法 因为Y(s) Z(s) _ Z(s) Y(s) U(s) Z(s) U(s) Z(s)n―1b )s m bs m bmQ ss ys 亠 亠a n 」s a n:丫(s) =(b °s m +bs m '+…+b m 」s + b m )Z(s)iU (s) = (s n+a 1s n，十■八 +a n/S + a n)Z(s)对上式取拉氏反变换，则jy =b o Z (m)+32^)+…+b m'Z + b m Z<(n )丄(n 4) I ■ ■ ■.u=z +az ++a n 』z+a n zX 2 = X 3G(s)二SISO按下列规律选择状态变量,即设x 1 二 z, X 2 二乙 ,X n(nd),于是有_x ；l - 0ir x j 「0] X 2■01—4y 二[b 2 b 1 b °] X 2 =[30] f uX 1X ；式中，|心为n -1 A 系数阵称之为友阵。

只要系统状态方程的系数阵 A 和输入阵b 具有上式的形式，c 阵的形式可以任意, 则称之为能控标准型。

则输出方程y 二 b °X n b i X n 」b m 」X 2 b m X i写成矩阵形式_X L IX 2y = [b m b m」b 1 b 0 ]'X n 」」n 一分析A,b,c 阵的构成与传递函数系数的关系。

现代控制理论试题B 卷及答案一、1 系统[]210,01021x x u y x ⎡⎤⎡⎤=+=⎢⎥⎢⎥-⎣⎦⎣⎦能控的状态变量个数是，能观测的状态变量个数是cvcvx 。

2试从高阶微分方程385y y y u ++=求得系统的状态方程和输出方程（4分/个） 解 12。

…..233118x x x x y x ==--=010080x ⎡⎢=⎢⎢-⎣分） 00⎣(5分)解 1．答：若存在控制向量序列(),(1),,(1)u k u k u k N ++-，时系统从第k 步的状态()x k 开始，在第N 步达到零状态，即()0x N =，其中N 是大于0的有限数，那么就称此系统在第k 步上是能控的。

若对每一个k ，系统的所有状态都是能控的，就称系统是状态完全能控的，简称能控。

…..….…….（3分）2.[][]320300020012 110-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=CA ………..……….（1分） [][]940300020012 3202=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=CA ……..……….（1分） ⎤⎡⎤⎡110C 1分）0140x ⎡⎤=⎢⎥⎣⎦ ()⎥⎦⎢⎢⎢⎣-=-8181881C U ……..…………..…….…….（1分） 11188P ⎡⎤=-⎢⎥⎣⎦……..………….…..…….…….（1分） ⎦⎤⎢⎣⎡=43412P ……..………….…...…….…….（1分）1314881148P -⎡⎤-⎢⎥=⎢⎥--⎢⎥⎣⎦..………….…...…….…….（1分） 101105C A PAP -⎡⎤==⎢⎥-⎣⎦………….…...…….…….（1分） ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-==1011 43418181Pb b C ……….…...…….…….（1分）1分） 解（3分） 3分）2分）（81分）11121112221222420261p p p p p ⎪-+=⎨⎪-=-⎩………...……....…….…….（1分） 112212743858p p p ⎧=⎪⎪=⎨⎪=⎪⎩………...…………....…….…….（1分）1112122275485388p p P p p ⎡⎤⎡⎤⎢⎥==⎢⎥⎢⎥⎣⎦⎣⎦...…………....…….…….（1分） 111211122275717480 det det 05346488p p P p p ⎡⎤⎡⎤⎢⎥=>==>⎢⎥⎢⎥⎣⎦⎣⎦………...（1分） P 正定，因此系统在原点处是大范围渐近稳定的.………（1分）八、给定系统的状态空间表达式为1010x --⎡⎢=-⎢⎢⎣2322213332223321(21)3313332(3)(26)64E E E E E E E E E E E λλλλλλλλλλ=+++++++++++++=+++++++++ -- 2分 又因为 *32()331f λλλλ=+++ ------- 1分列方程32123264126333E E E E E E +++=++=+= ----- 2分1232,0,3E k E =-==- ----------- 1分观测器为10312ˆˆ0110010113x x u y ---⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=-++⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦------- 1分 方法 2λ⋅分 分分分10ˆ0110x -⎡⎢=-⎢⎢⎣九 分） 1200A tAt A t e e e ⎛⎫= ⎪⎝⎭1A t t e e =…………………………..……….（1分） 11210()12s sI A s ---⎛⎫-= ⎪--⎝⎭101111212s s s s ⎛⎫ ⎪-= ⎪ ⎪- ⎪---⎝⎭………..……….（1分）(){}2112220t A t t t t e e L sI A e ee --⎛⎫=-= ⎪-⎝⎭……….…（1分）()112200000t At tt tt e e L sI A e e e e --⎛⎫ ⎪⎡⎤=-= ⎪⎣⎦ ⎪-⎝⎭……….……….（2分） 222001000001t t tt t t t e e e e e e e ⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪== ⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭……………..……….（2分）一、（（ × （ × （ √ （ √二、（的能控标准型、能观标准型和对角线标准型，并画出能控标准型的状态变量图。

tf2ss函数TF2SS函数是一种将传递函数转换为状态空间模型的MATLAB函数。

在控制系统中，状态空间模型是一种常见的表示方法，它使用一组状态变量和一个控制输入来描述系统的动态行为。

TF2SS函数将传递函数转换为这种状态空间模型，并可以在控制系统分析和设计中使用。

TF2SS函数的语法为：[A,B,C,D]=tf2ss(num,den)，其中num是传递函数的分子多项式，den是传递函数的分母多项式。

返回的A、B、C、D是状态空间模型的矩阵，其中A是状态转移矩阵，B是输入矩阵，C是输出矩阵，D是传递矩阵。

下面将详细介绍TF2SS函数及其应用。

一、TF2SS函数的原理在信号处理和控制系统中，传递函数通常被用来描述系统的输入和输出之间的关系。

传递函数是一个复数函数，表示输出与输入的比值。

状态空间模型是一种更一般的描述系统的方法，它包含状态变量和控制输入，描述了过程的完整的物理和数学过程，包括状态、输入和输出之间的所有关系。

状态空间模型是一组微分方程，它可以从初始状态开始通过控制输入来预测未来的系统行为，通常还包括噪声。

TF2SS函数可以将传递函数转换为状态空间形式，由此得到了系统的状态空间表达式，这有利于控制系统的建模和分析。

TF2SS函数可以通过一系列步骤来使用。

1. 定义传递函数传递函数可以定义为一个分子和一个分母多项式的比值。

在MATLAB中，这可以通过一组向量来表示。

num=[1 2] %分子多项式den=[1 3 2] %分母多项式G=tf(num,den) %定义一个传递函数变量2. 调用TF2SS函数G变量包含了一个传递函数，可以通过调用TF2SS函数将其转换为状态空间模型。

[A,B,C,D]=tf2ss(num,den) %将传递函数转换为状态空间模型3. 检查结果转换后的结果由四个矩阵组成：状态转移矩阵A，输入矩阵B，输出矩阵C和传递矩阵D。

可以使用MATLAB命令disp来查看这四个输出。

现代控制理论试题B 卷及答案一、1 系统[]210,01021x x u y x ⎡⎤⎡⎤=+=⎢⎥⎢⎥-⎣⎦⎣⎦能控的状态变量个数是，能观测的状态变量个数是cvcvx 。

2试从高阶微分方程385y y y u ++=求得系统的状态方程和输出方程（4分/个）解 1． 能控的状态变量个数是2，能观测的状态变量个数是1。

状态变量个数是2。

…..（4分）2．选取状态变量1x y =，2x y =，3x y =，可得 …..….…….（1分）12233131835x x x x x x x u y x ===--+= …..….…….（1分）写成010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦…..….…….（1分）[]100y x = …..….…….（1分）二、1给出线性定常系统(1)()(),()()x k Ax k Bu k y k Cx k +=+=能控的定义。

（3分）2已知系统[]210 020,011003x x y x ⎡⎤⎢⎥==⎢⎥⎢⎥-⎣⎦，判定该系统是否完全能观？(5分)解 1．答：若存在控制向量序列(),(1),,(1)u k u k u k N ++-，时系统从第k 步的状态()x k 开始，在第N 步达到零状态，即()0x N =，其中N 是大于0的有限数，那么就称此系统在第k 步上是能控的。

若对每一个k ，系统的所有状态都是能控的，就称系统是状态完全能控的，简称能控。

…..….…….（3分） 2.[][]320300020012 110-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=CA ………..……….（1分） [][]940300020012 3202=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=CA ……..……….（1分） ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=940320110 2CA CA C U O ………………..……….（1分） rank 2O U n =<，所以该系统不完全能观……..….…….（2分）三、已知系统1、2的传递函数分别为2122211(),()3232s s g s g s s s s s -+==++-+求两系统串联后系统的最小实现。

传递函数和状态空间方程引言传递函数和状态空间方程是控制系统工程中常用的数学模型和分析工具。

它们用于描述和分析动态系统的行为和性能，对于控制系统的设计和优化起着关键作用。

传递函数定义在控制系统中，传递函数是一个描述输入和输出之间关系的数学函数。

传递函数通常用G(s)表示，其中s是复数变量，表示系统的复频域特性。

传递函数描述了一个线性、时不变系统对输入信号的响应。

传递函数的一般形式如下：b0*s^n + b1*s^(n-1) + ... + bnG(s) = ---------------------------------------s^m + a1*s^(m-1) + ... + am其中n和m分别是传递函数的分子和分母的最高次幂。

用途传递函数可用于描述系统的频率响应和稳定性特性。

传递函数可以反映系统对不同频率的输入信号的放大或衰减情况，帮助工程师了解系统的动态特性。

传递函数还可以用于控制系统的设计和分析。

通过对传递函数进行数学运算和变换，可以获得系统的稳定性、动态响应以及频域特性等关键性能指标。

工作方式传递函数的输入是一个复数变量s，代表系统的频域特性。

通过将s带入传递函数的表达式中，可以得到系统的输出。

传递函数的输出代表了系统对输入信号的响应。

通过对传递函数表达式进行分析和计算，可以获得系统的稳定性、频率响应和动态响应等关键性能指标。

状态空间方程定义在控制系统中，状态空间方程是一种用状态变量表示系统状态的数学模型。

状态空间方程描述了系统的状态和状态变化随时间的规律。

状态空间方程的一般形式如下：dx/dt = Ax + Buy = Cx + Du其中，x是系统的状态向量，表示系统的状态变量；u是系统的输入向量，表示系统的输入信号；y是系统的输出向量，表示系统的输出信号；A、B、C和D是系统的系数矩阵。

用途状态空间方程可以用于描述和分析系统的动态行为和稳定性特性。

状态空间方程是一种直观、物理意义明确的模型，可以帮助工程师理解系统的内部状态和相互关系。

由传递函数转换成状态空间模型——方法多!!! SISO 线性定常系统高阶微分方程化为状态空间表达式SISO ()()()()()()m n u b u b u b y a y a y a y m m m n n n n ≥+++=++++--- 1102211)(2211110nn n n mm m a s a s a s b s b s b s G +++++++=--- 假设1+=m n外部描述←—实现问题：有了部结构—→模拟系统部描述SISO ⎩⎨⎧+=+=ducx y bu Ax x实现问题解决有多种方法，方法不同时结果不同。

一、直接分解法因为1011111()()()()()()()()1m m m mn n n nY s Z s Z s Y s U s Z s U s Z s b s b s b s b s a s a s a ----⨯=⨯=⨯++++++++⎩⎨⎧++++=++++=----)()()()()()(1111110s Z a s a s a s s U s Z b s b s b s b s Y n n n n m m m m 对上式取拉氏反变换，则⎩⎨⎧++++=++++=----z a z a za z u zb z b z b z b y n n n n m m m m 1)1(1)(1)1(1)(0 按下列规律选择状态变量，即设)1(21,,,-===n n z x zx z x ，于是有 ⎪⎪⎩⎪⎪⎨⎧+----===-u x a x a x a xx xx xn n n n 12113221写成矩阵形式式中，1-n I 为1-n 阶单位矩阵，把这种标准型中的A 系数阵称之为友阵。

只要系统状态方程的系数阵A 和输入阵b 具有上式的形式，c 阵的形式可以任意，则称之为能控标准型。

则输出方程121110x b x b x b x b y m m n n ++++=--写成矩阵形式⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=--n n m m x x x x b b b b y 121011][ 分析c b A ,,阵的构成与传递函数系数的关系。

系统建模与仿真习题二及答案1. 考虑如图所示的典型反馈控制系统框图(1)假设各个子传递函数模型为66.031.05.02)(232++-+=s s s s s G ，s s s G c 610)(+=，21)(+=s s H 分别用feedback （）函数以及G*Gc/(1+G*Gc*H)（要最小实现）方法求该系统的传递函数模型。

(2) 假设系统的受控对象模型为s e s s s G 23)1(12)(-+=，控制器模型为 ss s G c 32)(+=，并假设系统是单位负反馈，分别用feedback （）函数以及G*Gc/(1+G*Gc*H)（要最小实现）方法能求出该系统的传递函数模型？如果不能，请近似该模型。

解：（1）clc;clear;G=tf([2 0 0.5],[1 -0.1 3 0.66]);Gc=tf([10 6],[1 0]);H=tf(1,[1 2]);G1=feedback(G*Gc,H)G2=G*Gc/(1+G*Gc*H)Gmin=minreal(G2)结果：Transfer function:20 s^4 + 52 s^3 + 29 s^2 + 13 s + 6s^5 + 1.9 s^4 + 22.8 s^3 + 18.66 s^2 + 6.32 s + 3Transfer function:20 s^8 + 50 s^7 + 83.8 s^6 + 179.3 s^5 + 126 s^4 + 57.54 s^3 + 26.58 s^2 + 3.96 ss^9 + 1.8 s^8 + 25.61 s^7 + 22.74 s^6 + 74.11 s^5 + 73.4 s^4 + 30.98 s^3+ 13.17 s^2 + 1.98 s Transfer function:20 s^4 + 52 s^3 + 29 s^2 + 13 s + 6s^5 + 1.9 s^4 + 22.8 s^3 + 18.66 s^2 + 6.32 s + 3(2)由于s c e s s s s G s G 232)1(3624)(*)(-++= 方法1：将s e 2-转换为近似多项式。