当前位置：文档库 › 现代控制理论矩阵的复习

现代控制理论矩阵的复习

矩阵的复习(现代控制理论用)

(线性代数的复习)

一．矩阵

矩阵定义为矩形阵列表，它的元素可能是实数、复数、函数或算子，表中每一个元素表示一定的数学概念或工程信息，则此矩形表整体称为矩阵。

???

?? ??=nn n n n n a a a a a a a a a A 2

22221

11211 ???

??=jk j k k b b b b b b B 1221111 其中A 具有n 行与m 列，ij a 表示矩阵A 的第i 行，第j 列元素。矩阵B 为第j 行，k 列，即=A （ij a ），=B （jk b ）当n=m 时称为方阵，可表示为

??=?nn n n n n a a a a a a a a a 2

2222111211

＝det(A)

称为方阵A 的n 阶行列式。 1．子式或子行列式

给定n 阶行列式及其任一元素jk a ，划去△的j 行第k 列的全部元素而得到一个（n-1）阶行列式，称为原行列式的一个子行列式，简称子式，以jk ?表示。

即:33

2322

131211

a a a a a a a a a =? 划掉23a 的子行列式 32

3112

1123a a a a =

? 23?=?jk j=2 k=3 23a a jk =

n 阶行列式则有n 行n 列共2n 个子式。

2．余因子或代数余子式

如果用()

j +-1去乘jk a 的子式jk ?得到的结果称为jk a 的余因子或代数余子式。

()

2231+-=A -

=?2332

3112

11a a a a 一般为()

jk k

j jk A ?-=+1 ，其余类推。

3．行列式的值行列式的拉普拉斯展开式为任一行或任一列的各元素与对应的余因子的乘积之和。

232221

1211

a a a a a a a a a =?＝()111+-()()23

131231

333

131221

233

2322

1111a a a a a a a a a a a a a a

a ++-+-+

上式为第一列各元素展开的拉普拉斯展开式，任一行、任一列的元素展开结果之△值都相同。一般为 ()∑==

=?n

k jk jk

A a

a 1

det

4．行列式的基本性质

○

1将行列式的行与列交换其值不变。 ○

2将行列式的任两行（行与行）或两列交换，其值不变但符号相反。 ○

3将行列式任一行（任一列）的元素，乘以一个给定的数以后再加到另外的任一行（或任一列）的对应元素上去，行列式的值不变。

○

4若行列式中任一行或任一列的元素为零，则该行列式的值为零。 ○

5若有两行或两列的元素相同或成比例，即任一行或列的元素，刚好是与一行或列的元素的k 倍，则行

列式的值亦为零。

○

6用常数λ乘行列式任一行（或任一列）的诸元素，等于用λ乘这行列式。即：i 在行列式中任一行（或列）元素的公因子可以提到行列式之外。

ii 在n ×n 阵A 乘λ与每个元素，等于乘n

λ ，与上边是有区别的。A A n

λλ= 二．对角线矩阵

方阵A 中，由j=k 的全部元素jk a 所组成的对角线称为主对角线，除主对角线上各元素外，其余元素均为零者，称为对角线矩阵。

()nn nn a a a diag a a a A

1211

00=????

三．单位矩阵或么阵

当对角线矩阵中的主对角线各元素均为1时，其余的元素都为零。12211====nn a a a

这种特殊的对角线矩阵称为单位矩阵或么阵。()11110101

diag I =????

? ??=

四．向量

只有一列的矩阵????

? ??n x x x 21称为列向量，n 维称为n 维列向量。

只有一行的矩阵()n x x x 21

称为n 维行向量。

五．零阵所有元素均为零的矩阵。六．奇异矩阵

如果方阵相关行列式之值为零时，此方阵称为奇异矩阵。

七．非奇异矩阵

方阵相关行列式之值不为零者。八．转置矩阵

如果n ×m 矩阵A 的行与列交换，所得m ×n 矩阵称为该矩阵的转置矩阵，矩阵A 的转置矩阵用T

A 表示。

??????? ??=nm n n m m a a a a a a a a a A 2

22221

11211

??=nm m

n n T

a a a a a a a a a A 21222

1212111

()

A A T

T =

九．对称矩阵

如果方阵A 等于它的转置矩阵，即A =T

A 则方阵称A 为对称矩阵。十．正定矩阵

如果方阵A 的所有主子行列式为正，则矩阵A 称为正定矩阵。

??=nm n n m m a a a a a a a a a A 21

22221

11211

当满足 11a >0 ,

1211a a a a >0 ， … ，

n n m

m a a a a a a a a a

222

11211

>0 则A 为正定阵，这是判别正定矩阵的赛尔维斯特（sylverster ）准则。

还有共轭矩阵，赫末特矩阵（复数矩阵），逆矩阵，等等。

十一．矩阵的运算与变换

1．矩阵的加减法

)(jk jk b a B A +=+ )(jk jk b a B A -=-

矩阵的加减法满足交换律，结合律。C B A C B A ++=++)( 2．矩阵的乘法

数与矩阵相乘，等于数与矩阵各个元素相乘。)(jk ka kA =

矩阵与矩阵相乘，A 的各行元素分布与B 的各列元素对应乘积之和。

)(l j a A =是n ×m 矩阵 )(l k b B =是n ×r 矩阵 ∑==n

l lk jl jk b a C 1

所得矩阵C 的阶是m ×r

结合律、分配律适用于矩阵的乘法。

C AB BC A )()(= CA BA A C B +=+)(

3．两矩阵和与乘积的转置法则

T T T B A B A +=+)(

T T T A B AB =)( 注意 T T T B A AB ≠)(

4．逆矩阵

○

1定义：若对于给定的方阵A ，存在一个矩阵B ，使得AB=I ，则称B 是A 的逆矩阵，并以符号1

-A 表示之。

○

2矩阵的秩在一个n ×m 的矩阵A 中，任取k 行k 列，位于这些行列相交处的元素构成k 阶行列式，此行列式叫做A 的k 阶子式，而矩阵A 中不等于零的那个子矩阵的最大阶数，叫做矩阵A 的秩。若A 为方阵，而且其所有相当的行列式之值不为零，即满足条件:

)det(A ≠0

则称A 是满秩的n 阶方阵。

○

3余因子转置矩阵（或称伴随矩阵）将A 的余因子矩阵转置以后，（即行与列交换）所得之新矩阵，并以符号T

jk A )(或)(kj A 或adjA . ○

4逆矩阵的求法若A 是满秩的n 阶方阵，（即)det(A ≠0）则存在唯一的逆1

-A 。

A 1-A ＝1-A A ≡I

1-A 由下面公式确定

-A ＝

a a a a a a a a a A adjA

A A A A nn

n n

n n kj T jk

212

1212111)

det()

det()()

det()(==

○

5逆矩阵变换 111)(---=A B AB 注意111)(---≠B A AB

A A =--11)(

只有方阵才有逆矩阵，并且此方阵应是非奇异矩阵。（不等于零） 5．例1：2)(B A +＝？

解： 2)(B A +＝22))((B BA AB A B A B A +++=++

一般矩阵乘法不满足交换律，即A B ≠B A 当 A B ＝B A 时 2)(B A +＝222B AB A ++

例2：求A ＝???

??212104321的逆矩阵1-A 。

解：)det(A ＝－2

2214＋02231－11

1=－2(8－2)－1(1－12)=－12+11＝－1

伴随矩阵adjA ??????????------8341146211＝???

???

--04211

2211

204

1431223

122

141032213

1-A ＝????

??????----=????? ??-----?-=8341146211

834114621111det A A T

例3：求下列三个联立方程式的解：

1x ＋32x ＋3x ＝2

1x －22x ＋3x ＝3 21x +2x +33x ＝1

解：A =????? ??-312121131 ?

???

? ??=321x x x x ????? ??=132b

即b Ax = b A x 1-=

)det(A ＝

112-－3113＋21213-＝（－6－1）－（9－1）＋2×（3＋2）＝－5

伴随矩阵adjA ＝5

011

7----

1-A ＝

A A

det ???????

?----=1110515

115857 ????? ??321x x x ???????? ??----=1110515115857?????

??--=???????? ??+----+=?

???? ??45153313253521524514132 即：1x ＝5

33，2x ＝－51

，3x ＝－4

例4：确定下面矩阵的秩

A ＝??

??????-151431020211 解：??????????-151431020211→??????

???170310021→??????????-200031002

1→??????????--2000310621 该矩阵维数为3

现代控制理论基础考试题A卷及答案

即 112442k g k f M L M ML θθθ??=-+++ ??? && 212 44k k g M M L θθθ??=-+ ??? && （2）定义状态变量 11x θ=，21x θ=&，32 x θ=，42x θ=& 则一．（本题满分10分）如图所示为一个摆杆系统，两摆杆长度均为L ，摆杆的质量忽略不计，摆杆末端两个质量块（质量均为M ）视为质点，两摆杆中点处连接一条弹簧，1θ与2θ分别为两摆杆与竖直方向的夹角。当12θθ=时，弹簧没有伸长和压缩。水平向右的外力()f t 作用在左杆中点处，假设摆杆与支点之间没有摩擦与阻尼，而且位移足够小，满足近似式sin θθ=，cos 1θ=。（1）写出系统的运动微分方程；（2）写出系统的状态方程。【解】（1）对左边的质量块，有 ()2111211 cos sin sin cos sin 222 L L L ML f k MgL θθθθθθ=?-?-?-&& 对右边的质量块，有 ()221222 sin sin cos sin 22 L L ML k MgL θθθθθ=?-?-&& 在位移足够小的条件下，近似写成： ()1121 24f kL ML Mg θθθθ=---&& ()2122 4kL ML Mg θθθθ=--&&

2 / 7 1221 334413 44244x x k g k f x x x M L M ML x x k k g x x x M M L =?? ???=-+++ ???? ? =????=-+? ????? &&&& 或写成 11 223 34401 000014420001000044x x k g k x x M L M f ML x x x x k k g M M L ? ? ?? ?????????? ??-+???? ???????????=+???? ????? ??????????????????? ????-+?? ? ? ?????? ? &&&& 二．（本题满分10分）设一个线性定常系统的状态方程为=x Ax &，其中22R ?∈A 。若1(0)1?? =??-??x 时，状态响应为22()t t e t e --??=??-?? x ；2(0)1??=??-??x 时，状态响应为 2()t t e t e --?? =??-?? x 。试求当1(0)3??=????x 时的状态响应()t x 。【解答】系统的状态转移矩阵为()t t e =A Φ，根据题意有 221()1t t t e t e e --????==????--???? A x 22()1t t t e t e e --????==????--???? A x 合并得 2212211t t t t t e e e e e ----????=????----?? ??A 求得状态转移矩阵为 1 22221212221111t t t t t t t t t e e e e e e e e e -----------?????? ?? ==????????------???? ????A 22222222t t t t t t t t e e e e e e e e --------?? -+-+=??--??

(完整word版)现代控制理论基础试卷及答案,推荐文档

现代控制理论基础考试题西北工业大学考试题（A卷）（考试时间120分钟）学院：专业：姓名：学号：一．填空题（共27分，每空1.5分） 1.现代控制理论基础的系统分析包括___________和___________。 2._______是系统松弛时，输出量、输入量的拉普拉斯变换之比。 3.线性定常系统齐次状态方程是指系统___________时的状态方程。 4.推导离散化系统方程时在被控对象上串接一个开关，该开关以T为周期进行开和关。这个开关称为_______。 5.离散系统的能______和能______是有条件的等价。 6.在所有可能的实现中，维数最小的实现称为最小实现，也称为__________。 7.构造一个与系统状态x有关的标量函数V(x, t)来表征系统的广义能量， V(x, t)称为___________。 8.单输入-单输出线性定常系统，其BIBO稳定的充要条件是传递函数的所有极点具有______。 9.控制系统的综合目的在于通过系统的综合保证系统稳定，有满意的 _________、_________和较强的_________。 10.所谓系统镇定问题就是一个李亚普诺夫意义下非渐近稳定的系统通过引入_______，以实现系统在李亚普诺夫意义下渐近稳定的问题。 11.实际的物理系统中，控制向量总是受到限制的，只能在r维控制空间中某一个控制域内取值，这个控制域称为_______。 12._________和_________是两个相并行的求解最优控制问题的重要方法。二．判断题（共20分，每空2分） 1.一个系统，状态变量的数目和选取都是惟一的。（×） 2.传递函数矩阵的描述与状态变量选择无关。（√） 3.状态方程是矩阵代数方程，输出方程是矩阵微分方程。（×） 4.对于任意的初始状态) ( t x和输入向量)(t u，系统状态方程的解存在并且惟一。（√） 5.传递函数矩阵也能描述系统方程中能控不能观测部分的特性。（×） 6.BIBO 稳定的系统是平衡状态渐近稳定。（×） 7.一个系统能正常工作，稳定性是最基本的要求。（√） 8.如果系统的状态不能测得，只要系统能观测，可以采用状态观测器实现状

现代控制理论实验报告

实验报告 ( 2016-2017年度第二学期) 名称：《现代控制理论基础》题目：状态空间模型分析院系：控制科学与工程学院班级： ___ 学号： __ 学生姓名： ______ 指导教师： _______ 成绩：日期： 2017年 4月 15日

线控实验报告一、实验目的： l.加强对现代控制理论相关知识的理解； 2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析；二、实验内容 1 第一题：已知某系统的传递函数为G (s) S23S2 求解下列问题：（1）用 matlab 表示系统传递函数 num=[1]; den=[1 3 2]; sys=tf(num,den); sys1=zpk([],[-1 -2],1); 结果： sys = 1 ------------- s^2 + 3 s + 2 sys1 = 1 ----------- (s+1) (s+2) （2）求该系统状态空间表达式： [A1,B1,C1,D1]=tf2ss(num,den); A = -3-2 10 B = 1 C = 0 1

第二题：已知某系统的状态空间表达式为： 321 A ,B,C 01：10 求解下列问题：（1）求该系统的传递函数矩阵：（2）该系统的能观性和能空性：（3）求该系统的对角标准型：（4）求该系统能控标准型：（5）求该系统能观标准型：（6）求该系统的单位阶跃状态响应以及零输入响应：解题过程：程序： A=[-3 -2;1 0];B=[1 0]';C=[0 1];D=0; [num,den]=ss2tf(A,B,C,D); co=ctrb(A,B); t1=rank(co); ob=obsv(A,C); t2=rank(ob); [At,Bt,Ct,Dt,T]=canon(A,B,C,D, 'modal' ); [Ac,Bc,Cc,Dc,Tc]=canon(A,B,C,D, 'companion' ); Ao=Ac'; Bo=Cc'; Co=Bc'; 结果：（1） num = 0 01 den = 1 32 （2）能控判别矩阵为： co = 1-3 0 1 能控判别矩阵的秩为： t1 = 2 故系统能控。（3）能观判别矩阵为： ob = 0 1

现代控制理论课后习题答案

绪论为了帮助大家在期末复习中能更全面地掌握书中知识点，并且在以后参加考研考博考试直到工作中，为大家提供一个理论参考依据，我们11级自动化二班的同学们在王整风教授的带领下合力编写了这本《现代控制理论习题集》（刘豹第三版），希望大家好好利用这本辅助工具。根据老师要求，本次任务分组化，责任到个人。我们班整体分为五大组，每组负责整理一章习题，每个人的任务由组长具体分配，一个人大概分1~2道题，每个人任务虽然不算多，但也给同学们提出了要求：1.写清题号，抄题，画图（用CAD或word画）。2.题解详略得当，老师要求的步骤必须写上。3.遇到一题多解，要尽量写出多种方法。本习题集贯穿全书，为大家展示了控制理论的基础、性质和控制一个动态系统的四个基本步骤，即建模、系统辨识、信号处理、综合控制输入。我们紧贴原课本，强调运用统一、联系的方法分析处理每一道题，将各章节的知识点都有机地整合在一起，力争做到了对控制理论概念阐述明确，给每道题的解析赋予了较强的物理概念及工程背景。在课后题中出现的本章节重难点部分，我们加上了必要的文字和图例说明，让读者感觉每一题都思路清晰，简单明了，由于我们给习题配以多种解法，更有助于发散大家的思维，做到举一反三！

这本书是由11级自动化二班《现代控制理论》授课老师王整风教授全程监管，魏琳琳同学负责分组和发布任务书，由五个小组组组长李卓钰、程俊辉、林玉松、王亚楠、张宝峰负责自己章节的初步审核，然后汇总到胡玉皓同学那里，并由他做最后的总审核工作，绪论是段培龙同学和付博同学共同编写的。本书耗时两周，在同学的共同努力下完成，是二班大家庭里又一份智慧和努力的结晶，望大家能够合理使用，如发现错误请及时通知，欢迎大家的批评指正！ 2014年6月2日

现代控制理论实验

华北电力大学实验报告| | 实验名称状态空间模型分析课程名称现代控制理论 | | 专业班级：自动化1201 学生姓名：马铭远学号：2 成绩：指导教师：刘鑫屏实验日期：4月25日

状态空间模型分析一、实验目的 1.加强对现代控制理论相关知识的理解； 2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析；二、实验仪器与软件 1. MATLAB7.6 环境三、实验内容 1 、模型转换图 1、模型转换示意图及所用命令传递函数一般形式： MATLAB 表示为： G=tf(num,den)，，其中 num,den 分别是上式中分子，分母系数矩阵。零极点形式： MATLAB 表示为：G=zpk(Z,P,K) ，其中 Z，P ，K 分别表示上式中的零点矩阵，极点矩阵和增益。传递函数向状态空间转换：[A,B,C,D] = TF2SS(NUM,DEN)；状态空间转换向传递函数：[NUM,DEN] = SS2TF(A,B,C,D,iu)---iu 表示对系统的第 iu 个输入量求传递函数；对单输入 iu 为 1。

例1：已知系统的传递函数为G（S）= 2 2 3 24 11611 s s s s s ++ +++ ,利用matlab将传递函数和状态空间相互转换。解：1.传递函数转换为状态空间模型： NUM=[1 2 4];DEN=[1 11 6 11]; [A,B,C,D] = tf2ss(NUM,DEN) 2.状态空间模型转换为传递函数： A=[-11 -6 -11;1 0 0;0 1 0];B=[1;0;0];C=[1 2 4];D=[0];iu=1; [NUM,DEN] = ss2tf(A,B,C,D,iu); G=tf(NUM,DEN) 2 、状态方程状态解和输出解单位阶跃输入作用下的状态响应： G=ss(A,B,C,D);[y,t,x]=step(G);plot(t,x). 零输入响应 [y,t,x]=initial(G,x0)其中，x0 为状态初值。

现代控制理论考试试卷A

北京航空航天大学 2019－2020 学年第二学期期末《现代控制理论》 A卷班级______________学号 _________ 姓名______________成绩 _________ 2020年6月22日

班号学号姓名成绩《现代控制理论》期末考试卷一、（本题10分）某RLC 电路如题一图所示，其中u 为输入信号、y 为输出信号、i 为流过网络的电流。若令状态x 1=i ，x 2=y ，建立系统的动态方程，并判断系统的可控性和可观测性（所有参数非零）。题一图二、（本题10分）系统的动态方程为 010*********???? ????=+????-???????? x x u ， []001=y x 若[](0)001=-T x ，()()δ=u t t （单位脉冲信号），求()x t 和()y t 。三、（本题15分）已知系统具有如下形式： []111122********* a b x Ax bu a x b u b y cx c c c x l l l éù éùêúêúêúêú=+=+êúêú êúêú???? == (1). 若12=l l ，给出系统可控并且可观测的充分必要条件；若12≠l l ，20=b ，

给出系统可控的充分必要条件（即参数12123123,,,,,,,a a b b b c c c 需满足的条件）； (2). 若11=-l ，11=a ，[][]12123301,1000b b c c c b éùéù êúêú êúêú==êúêúêúêú??? ?，计算系统的传递函数()G s ，并给出该传递函数的可观标准型最小阶实现。四、（本题20分）已知系统具有如下形式： []1112212200 n n A A x Ax bu x u A A b y cx c x éùéù êúêú=+=+êúêú????== 其中， 11A 为(1)(1)-?-n n 的方阵，22A 为11?的方阵，12A 为(1)-n 维列向量，21A 为(1)-n 维行向量，n b 和n c 分别为非零实数。 (1). 证明系统既可控又可观测的充分必要条件是：1112(,)A A 可控且1121(,)A A 可观测； (2). 若A 的特征多项式为()p s ，而 110100001000011000 A éù êúêúêúêú=êúêúêú êú?? 求系统的传递函数，并证明若系统既可控又可观测，则有(1)0≠p 。五、（本题15分）已知系统动态方程如下： 210431x x u éùéù êúêú=+êúêúêúêú???? , 11y x éù=êú?? (1). 判断系统的可控性。若系统可控，将系统化为可控标准型； (2). 是否可以用状态反馈将A bk -的特征值配置到{}2,3--？若可以，求出状态反馈增益阵k 。

现代控制理论实验报告

实验一系统能控性与能观性分析一、实验目的 1.理解系统的能控和可观性。二、实验设备 1.THBCC-1型信号与系统·控制理论及计算机控制技术实验平台；三、实验内容二阶系统能控性和能观性的分析四、实验原理系统的能控性是指输入信号u对各状态变量x的控制能力，如果对于系统任意的初始状态，可以找到一个容许的输入量，在有限的时间内把系统所有的状态引向状态空间的坐标原点，则称系统是能控的。对于图21-1所示的电路系统，设iL和uc分别为系统的两个状态变量，如果电桥中则输入电压ur能控制iL和uc状态变量的变化，此时，状态是能控的。反之，当时，电桥中的A点和B点的电位始终相等，因而uc不受输入ur的控制，ur只能改变iL的大小，故系统不能控。系统的能观性是指由系统的输出量确定所有初始状态的能力，如果在有限的时间内根据系统的输出能唯一地确定系统的初始状态，则称系统能观。为了说明图21-1所示电路的能观性，分别列出电桥不平衡和平衡时的状态空间表达式：平衡时：

由式（2）可知，状态变量iL和uc没有耦合关系，外施信号u只能控制iL的变化，不会改变uc的大小，所以uc不能控。基于输出是uc，而uc与iL无关连，即输出uc中不含有iL的信息，因此对uc的检测不能确定iL。反之式（1）中iL与uc有耦合关系，即ur的改变将同时控制iL和uc的大小。由于iL与uc的耦合关系，因而输出uc的检测，能得到iL的信息，即根据uc的观测能确定iL（ω）五、实验步骤 1．用2号导线将该单元中的一端接到阶跃信号发生器中输出2上，另一端接到地上。将阶跃信号发生器选择负输出。 2．将短路帽接到2K处，调节RP2，将Uab和Ucd的数据填在下面的表格中。然后将阶跃信号发生器选择正输出使调节RP1，记录Uab和Ucd。此时为非能控系统，Uab和Ucd没有关系（Ucd始终为0）。 3．将短路帽分别接到1K、3K处，重复上面的实验。六、实验结果表20－1Uab与Ucd的关系

(完整版)现代控制理论考试卷及答案

西北工业大学考试试题（卷）2008 －2009 学年第2 学期

2009年《现代控制理论》试卷A 评分标准及答案第一题（10分，每个小题答对1分，答错0分）（1）对（2）错（3）对（4）错（5）对（6）对（7）对（8）对（9）对（10）错第二题（15分）（1）)(t Φ（7分）：公式正确3分，计算过程及结果正确4分 ? ? ? ???+-+---=-=Φ?? ?? ??????+- +-+- +-+- ++-+=??????-+++=-??? ???+-=------------t t t t t t t t e e e e e e e e A sI L t s s s s s s s s s s s s A sI s s A sI 22221 11 2222}){()(22112 21221112112 213)2)(1(1 )(321 （2）状态方程有两种解法（8分）：公式正确4分，计算过程及结果正确4分 ??????-+-+-=????? ???????+-+++-+++-++??????+--=??????????? ???????++-++++-=-+-=??????---+-=????? ?+--+??? ???+--=??????-Φ+Φ=------------------------------??t t t t t t t t t t t t t t t t t t t t t e e te e e te s s s s s s L e e e e t x t x s s s s s L x A sI L t x s BU A sI x A sI s X e e t e e t d e e e e e e e e e t x t x d t Bu x t t x 222 21 22212 21111122)(02222210 2344}2414)1(42212)1(4 {2)()(} )2()1(4) 2()1()3(2{)}0(){()() ()()0()()(2)34()14(22222)()()()()0()()(或者 ττ τττττττ 第三题（15分，答案不唯一，这里仅给出可控标准型的结果）（1）系统动态方程（3分） []x y u x x 0010 1003201 00010=???? ??????+??????????--=&

《现代控制理论基础》考试题B卷及答案

一．（本题满分10分）请写出如图所示电路当开关闭合后系统的状态方程和输出方程。其中状态变量的设置如图所示，系统的输出变量为流经电感2L 的电流强度。【解答】根据基尔霍夫定律得： 1113222332 1L x Rx x u L x Rx x Cx x x ++=?? +=??+=? 改写为1 13111 22 322 312 11111R x x x u L L L R x x x L L x x x C C ? =--+?? ?=-+???=-?? ，输出方程为2y x = 写成矩阵形式为

[]11 111222 2 331231011000110010R L L x x L R x x u L L x x C C x y x x ??? --???????????????? ???????=-+???? ??????? ??????????????? ? ???-?????? ? ? ??? ?? ?=??? ?????? 二．（本题满分10分）单输入单输出离散时间系统的差分方程为 (2)5(1)3()(1)2()y k y k y k r k r k ++++=++ 回答下列问题：（1）求系统的脉冲传递函数；（2）分析系统的稳定性；（3）取状态变量为1()()x k y k =，21()(1)()x k x k r k =+-，求系统的状态空间表达式；（4）分析系统的状态能观性。【解答】（1）在零初始条件下进行z 变换有： ()()253()2()z z Y z z R z ++=+ 系统的脉冲传递函数： 2()2 ()53 Y z z R z z z +=++ （2）系统的特征方程为 2()530D z z z =++= 特征根为1 4.3z =-，20.7z =-，11z >，所以离散系统不稳定。（3）由1()()x k y k =，21()(1)()x k x k r k =+-，可以得到 21(1)(2)(1)(2)(1)x k x k r k y k r k +=+-+=+-+ 由已知得 (2)(1)2()5(1)3()y k r k r k y k y k +-+=-+-112()5(1)3()r k x k x k =-+- []212()5()()3()r k x k r k x k =-+-123()5()3()x k x k r k =--- 于是有： 212(1)3()5()3()x k x k x k r k +=--- 又因为 12(1)()()x k x k r k +=+ 所以状态空间表达式为

现代控制理论基础考试题B卷及答案

-----好资料学习分）一．（本题满分10请写出如图所示电路当开关闭合后系统的状态方程和输出方程。其中状L态变量的设置如图所示，系统的输出变量为流经电感的电流强度。2

【解答】根据基尔霍夫定律得：uLx?Rx?x??3111 ?x?Lx?Rx?3222 ?xx?Cx??213 1R1?ux?x??x?? 311LLL ?1111R?x??x?x?232x?y，输出方程为改写为LL ?222?11x?x?x? 123CC? 写成矩阵形式为更多精品文档．学习-----好资料 ?R1??1??0?????LL?????11Lxx??????111 ?1R??????u?x?0?x0?????????22LL ?????22????0xx???????33??11

???0????? ?CC???x???1?????x1y?00???2???x???3 10分）二．（本题满分单输入单输出离散时间系统的差分方程为)k2r(r?3y(k)?(k?1)??y(k2)?5y(k?1) 回答下列问题：）求系统的脉冲传递函数；（1 ）分析系统的稳定性；（2)y?(kx(k))r(kx(k)?x(k?1)?，，（3）取状态变量为求系统的状态空间表达式；112（4）分析系统的状态能观性。【解答】z变换有：1（）在零初始条件下进行????2 )z?2)?zRz(?5z?3zY(2?(Yz)z?系统的脉冲传递函数： 23R(z)z?5z?（2）系统的特征方程为20?5?z?3zD(z)? 1z?0.7?z?4.3??z，，所以离散系统不稳定。，特征根为211)(k1)?rx)?y(k)x(k)?(k?(xk 3）由，，可以得到（1211)(k??(k?1)y(k?2)?r?kx(?1)?x(k2)?r12由已知得)?1)?3x(kk(?2rk)?5x()k3?1)?y(k?r??(yk?2)r(k1)?2(k)5y(11??)x(k?5)x(k?r(k)3?)2?r(k)(x ?5(k)?3rkk3??x()2112于是有：)k3(?(?1)?3xk)5xk)?r(?(xk221又因为)?k(??(xk1)x)r(k21所以状态空间表达式为更多精品文档．学习-----好资料 ?x(k?1)x(k)101????????11??r(k)?????????x?3?3?5(x(k)k?1)?????????22 ?x(k)?????101y(k)????x(k)???2（4）系统矩阵为0101??????????，输出矩阵为0c?110?0G?cG?1，?????3?5?3?5????c10????能观性矩阵为，，系统完全能观。2Q?rank??Q????oo cG01???? 三．（本题满分10分）回答下列问题：（1）简述线性系统的对偶原理；（2）简述线性定常系统的状态稳定性与输出稳定性的相互关系； r?2rr阶线性解耦系统等效于多少个独立的单输入单输出系统？输出（3）输入【解答】（1）若线性系统1与线性系统2互为对偶，则系统1的能控性等价于系统2的能观性，系统1的能观性等价于系统2的能控性。（2）若线性定常系统的状态稳定，则输出必稳定，反之，若线性定常系统的输出稳定，则状态未必稳定。当且仅当线性定常系统的传递函数没有零极点对消现象时，其状态稳定性和输出稳定性才是等价的。 r?2rrr个独立的单输入单输出系统。输入）输出阶线性解耦系统等效于（3 四．（本题满分10分） x?x?x cos x?2211?，判

哈尔滨工业大学《现代控制理论基础》考试题B卷及答案

哈工大2010 年春季学期现代控制理论基础试题B 答案题号一二三四五六七八卷面分作业分实验分总分满分值 10 10 10 10 10 10 10 10 80 10 10 100 得分值第 1 页（共 8 页）班号姓名一．（本题满分10分）请写出如图所示电路当开关闭合后系统的状态方程和输出方程。其中状态变量的设置如图所示，系统的输出变量为流经电感2L 的电流强度。【解答】根据基尔霍夫定律得： 1113222332 1L x Rx x u L x Rx x Cx x x ++=?? +=??+=? 改写为1 13111 22 322 31 211111R x x x u L L L R x x x L L x x x C C ? =--+?? ?=-+???=-?? ，输出方程为2y x = 写成矩阵形式为

现代控制理论实验报告

现代控制理论实验报告组员：院系：信息工程学院专业：指导老师：年月日

实验1 系统的传递函数阵和状态空间表达式的转换 [实验要求] 应用MATLAB 对系统仿照[例]编程，求系统的A 、B 、C 、阵；然后再仿照[例]进行验证。并写出实验报告。 [实验目的] 1、学习多变量系统状态空间表达式的建立方法、了解系统状态空间表达式与传递函数相互转换的方法； 2、通过编程、上机调试，掌握多变量系统状态空间表达式与传递函数相互转换方法。 [实验内容] 1 设系统的模型如式示。 p m n R y R u R x D Cx y Bu Ax x ∈∈∈?? ?+=+=& 其中A 为n ×n 维系数矩阵、B 为n ×m 维输入矩阵 C 为p ×n 维输出矩阵，D 为传递阵，一般情况下为0，只有n 和m 维数相同时，D=1。系统的传递函数阵和状态空间表达式之间的关系如式示。 D B A SI C s den s num s G +-== -1)() () (()( 式中，)(s num 表示传递函数阵的分子阵，其维数是p ×m ；)(s den 表示传递函数阵的按s 降幂排列的分母。 2 实验步骤 ① 根据所给系统的传递函数或（A 、B 、C 阵），依据系统的传递函数阵和状态空间表达式之间的关系如式，采用MATLA 的编程。注意：ss2tf 和tf2ss 是互为逆转换的指令； ② 在MATLA 界面下调试程序，并检查是否运行正确。 ③ [] 已知SISO 系统的状态空间表达式为，求系统的传递函数。

, 2010050010000100001 0432143 21u x x x x x x x x ? ? ??? ? ??????-+????????????????????????-=????????????&&&&[]??? ? ? ???????=43210001x x x x y 程序: A=[0 1 0 0;0 0 -1 0;0 0 0 1;0 0 5 0]; B=[0;1;0;-2]; C=[1 0 0 0]; D=0; [num,den]=ss2tf(A,B,C,D,1) 程序运行结果: num = 0 den = 0 0 0 从程序运行结果得到:系统的传递函数为: 2 4253 )(s s s S G --= ④ [] 从系统的传递函数式求状态空间表达式。程序: num =[0 0 1 0 -3]; den =[1 0 -5 0 0]; [A,B,C,D]=tf2ss(num,den) 程序运行结果: A = 0 5 0 0 1 0 0 0 0 1 0 0

《现代控制理论》.

《现代控制理论》实验指导书俞立徐建明编浙江工业大学信息工程学院 2007年4月

实验1 利用MATLAB 进行传递函数和状态空间模型间的转换 1.1 实验设备 PC 计算机1台（要求P4-1.8G 以上），MATLAB6.X 或MATLAB7.X 软件1套。 1.2 实验目的 1、学习系统状态空间模型的建立方法、了解状态空间模型与传递函数相互转换的方法； 2、通过编程、上机调试，掌握系统状态空间模型与传递函数相互转换的方法。 1.3 实验原理说明设系统的状态空间模型是 x Ax Bu y Cx Du =+?? =+?& (1.1) p y R ∈其中：n x R ∈是系统的状态向量，是控制输入，m u R ∈是测量输出，A 是维状态矩阵、是维输入矩阵、是n n ×m n ×n p ×B D C 维输出矩阵、是直接转移矩阵。系统传递函数和状态空间模型之间的关系如式（1.2）所示。 1()()G s C sI A B D ?=?+ (1.2) 表示状态空间模型和传递函数的MATLAB 函数。函数ss （state space 的首字母）给出了状态空间模型，其一般形式是 SYS = ss(A,B,C,D) 函数tf （transfer function 的首字母）给出了传递函数，其一般形式是 G=tf(num,den) 其中的num 表示传递函数中分子多项式的系数向量（单输入单输出系统），den 表示传递函数中分母多项式的系数向量。函数tf2ss 给出了传递函数的一个状态空间实现，其一般形式是 [A,B,C,D]=tf2ss(num,den) 函数ss2tf 给出了状态空间模型所描述系统的传递函数，其一般形式是 [num,den]=ss2tf(A,B,C,D,iu) 其中对多输入系统，必须确定iu 的值。例如，若系统有三个输入和，则iu 必须是1、2或3，其中1表示，2表示，3表示。该函数的结果是第iu 个输入到所有输出的传递函数。 21,u u 3u 1u 2u 3u 1.4 实验步骤 1、根据所给系统的传递函数或（A 、B 、C 、D ），依据系统的传递函数阵和状态空间模型之间的关系（1.2），采用MATLAB 的相关函数编写m-文件。 2、在MATLAB 界面下调试程序。例1.1 求由以下状态空间模型所表示系统的传递函数， ?? ? ? ? ?????=?????? ?????+???????????????????????=??????????321321321]001[1202505255100010x x x y u x x x x x x &&&

现代控制理论试题(详细答案)

现代控制理论试题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分） 12233131 835x 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 (),()3232 s s g s g s s s s s -+==++-+ 求两系统串联后系统的最小实现。（8分）解 112(1)(1)11 ()()()(1)(2)(1)(2)4 s s s s g s g s g s s s s s s -+++== ?=++--- …..….……. （5分）最小实现为

现代控制理论课程学习心得.

现代控制理论基础课程总结学院:__机械与车辆学院_ 学号:____2120120536___ 姓名:_____王文硕______ 专业:___交通运输工程__ 《现代控制理论》学习心得摘要:从经典控制论发展到现代控制论,是人类对控制技术认识上的一次飞跃。现代控制论是用状态空间方法表示,概念抽象,不易掌握。对于《现代控制理论》这门课程,本人选择了最为感兴趣的几个知识点进行分析,并谈一下对于学习这么课程的一点心得体会。关键词:现代控制理论;学习策略;学习方法;学习心得在现代科学技术飞速发展中,伴随着学科的高度分化和高度综合,各学科之间相互交叉、相互渗透,出现了横向科学。作为跨接于自然科学和社会科学的具有横向科学特点的现代控制理论已成为我国理工科大学高年级的选修课和研究生的学位课。从经典控制论发展到现代控制论,是人类对控制技术认识上的一次飞跃。经典控制论限于处理单变量的线性定常问题,在数学上可归结为单变量的常系数微分方程问题。现代控制论面向多变量控制系统的问题,它是以矩阵论和线性空间理论作为主要数学工具,并用计算机来实现。现代控制论来源于工程实际,具有明显的工程技术特点,但它又属于系统论范畴。系统论的特点是在数学描述的基础上,充分利用现有的强有力的数学工具,对系统进行分析和综合。系统特性的度量,即表现为状态;系统状态的变化,即为动态过程。状态和过程在自然界、社会和思维中普遍存在。现代控制论是在引入状态和状态空间的概念基础上发展起来的。状态和状态空间早在古典动力学中得到了广泛的应用。在5O年代Mesarovic教授曾提出“结构不确定

性原理”,指出经典理论对于多变量系统不能确切描述系统的内在结构。后来采用状态变量的描述方法,才完全表达出系统的动力学性质。6O年代初,卡尔曼(Kalman从外界输入对状态的控制能力以及输出对状态的反映能力这两方面提出能控制性和能观性的概念。这些概念深入揭示了系统的内在特性。实际上,现代控制论中所研究的许多基本问题,诸如最优控制和最佳估计等,都是以能能控性和能观性作为“解”的存在条件的。现代控制理论是一门工程理论性强的课程,在自学这门课程时,深感概念抽象,不易掌握;学完之后,从工程实际抽象出一个控制论方面的课题很难,如何用现代控制论的基本原理去解决生产实际问题则更困难,这是一个比较突出的矛盾。对现代控制理论来说,首先遇到的问题是将实际系统抽象为数学模型,有了数学模型,才能有效地去研究系统的各个方面。许多机电系统、经济系统、管理系统常可近似概括为线性系统。线性系统和力学中质点系统一样,是一个理想模型,理想模型是研究复杂事物的主要方法,是对客观事物及其变化过程的一种近似反映。现代控制论从自然和社会现象中抽象出的理想模型,用状态空间方法表示,再作理论上的探讨。线性系统理论是一门严谨的科学。抽象严谨是其本质的属性,一旦体会到数学抽象的丰富含义,再不会感到枯燥乏味。线性系统理论是建立在线性空间的基础上的,它大量使用矩阵论中深奥的内容,比如线性变换、子空间等,是分析中最常用的核心的内容,要深入理解,才能体会其物理意义。比如,状态空间分解就是一种数学分析方法。在控制论中把实际系统按能控性和能观性化分成四个子空间,它们有着确切的物理概念。线性变换的核心思想在于:线性系统的基本性质(如能控性、能观性、极点、传递函数等在线性变换下都不改变,从而可将系统化为特定形式,使问题的研究变得简单而透彻。在学习现代控制理论教材时,发现不少“引而未发”的问题。由于作者有丰富的教学经验与学术造诣,能深入浅出阐述问题,发人深省。因此,通过自己反复阅读教材,就能理解这些内容。比如,在探讨线性系统的传递函数的零极点相消时,如果潜伏着

现代控制理论实验报告

现代控制理论实验指导书实验一：线性系统状态空间分析 1、模型转换图1、模型转换示意图及所用命令传递函数一般形式： )()(11101110n m a s a s a s a b s b s b s b s G n n n n m m m m ≤++++++++= ---- MATLAB 表示为：G=tf(num,den)，其中num,den 分别是上式中分子，分母系数矩阵。零极点形式： ∏∏==--= n i j m i i p s z s K s G 1 1 ) () ()( MATLAB 表示为：G=zpk(Z,P,K)，其中 Z ，P ，K 分别表示上式中的零点矩阵，极点矩阵和增益。传递函数向状态空间转换：[A,B,C,D] = TF2SS(NUM,DEN)；状态空间转换向传递函数：[NUM,DEN] = SS2TF(A,B,C,D,iu)---iu 表示对系统的第iu 个输入量求传递函数；对单输入iu 为1; 验证教材P438页的例9-6。求P512的9-6题的状态空间描述。 >> A=[0 1;0 -2];

>> B=[1 0;0 1]; >> C=[1 0;0 1]; >> D=[0 0;0 0]; >> [NUM,DEN] = ss2tf(A,B,C,D,1) NUM = 0 1 2 0 0 0 DEN = 1 2 0 >> [NUM,DEN] = ss2tf(A,B,C,D,2) NUM = 0 0 1 0 1 0 DEN = 1 2 0 给出的结果是正确的，是没有约分过的形式 P512 9-6 >> [A,B,C,D]=tf2ss([1 6 8],[1 4 3])

现代控制理论期末试卷

一、（10分，每小题1分） 1、任一线性连续定常系统的系统矩阵均可对角形化。(×) 2、对SISO 线性连续定常系统，传递函数存在零极点对消，则系统一定不能观且不能控制。(×) 3、对线性连续定常系统，非奇异变换后的系统特征值不变。(√) 4、对于线性连续定常系统的最小实现是唯一的。(×) 5、稳定性问题是相对于某个平衡状态而言的。(√) 6、Lyapunov 第二法只给出了判定稳定性的充分条件。(√) 7、对于SISO 线性连续定常系统，状态反馈后形成的闭环系统零点与原系统一样。(√) 8、对于一个系统，只能选取一组状态变量。(×) 9、对于一个n 维的线性定常连续系统，若其完全能观，则利用状态观测器实现的状态反馈闭环系统是2n 维的。(√) 10、对线性定常系统，其Lyapunov 意义下的渐近稳定性和矩阵特征值都具有负实部是一致的。(√) 二（10分，每小题5分）（1）简述平衡状态及平衡点的定义。（2）简述状态方程解的意义。解：（1）状态空间中状态变量的导数向量为零向量的点。由平衡状态在状态空间中所确定的点称之为平衡点。（2）线性连续定常系统状态方程的解由两部分组成，一部分是由初始状态所引起的自由运动即零输入响应，第二部分是由输入所引起的系统强迫运动，与输入有关称为零状态响应。三、（10分）考虑如图的质量弹簧系统。其中，m 为运动物体的质量，k 为弹簧的弹性系数，h 为阻尼器的阻尼系数，f 为系统所受外力。取物体位移为状态变量x 1，速度为状态变量x 2，并取位移为系统输出y ，外力为系统输入u ，试建立系统的状态空间表达式。解： f ma =……………………………….……1分令位移变量为x 1，速度变量为x 2，外力为输入u ，有 122u kx kx mx --=………………………………2分于是有 12x x =………………………………..……………1分 2121k h x x x u m m m =--+……….….……………….2分再令位移为系统的输出y ，有