紫金学院计算机系实验报告
现代控制理论基础
实验报告
专业：______自动化_____________
年级：_______2011级______________
姓名：__________孙青山_________________
学号：____________110603152_______________ 提交日期：_____5.29__________________
实验一 系统能控性与能观性分析
1． 实验目的：
1．通过本实验加深对系统状态的能控性和能观性的理解；
2．验证实验结果所得系统能控能观的条件与由它们的判据求得的结果完全一致。 2． 实验内容：
1．线性系统能控性实验； 2．线性系统能观性实验 3． 实验原理：
系统的能控性是指输入信号u 对各状态变量x 的控制能力。如果对于系统任意的初始状态，可以找到一个容许的输入量，在有限的时间内把系统所有的状态变量转移到状态空间的坐标原点。则称系统是能控的。
系统的能观性是指由系统的输出量确定系统所有初始状态的能力。如果在有限的时间内，根据系统的输出能唯一地确定系统的初始状态，则称系统能观。
对于图10-1所示的电路系统，设i L 和u c 分别为系统的两个状态变量，如果电桥中
4
3
21R R R R ≠， 则输入电压u 能控制i L 和u c 状态变量的变化，此时，状态是能控的；状态变量i L 与u c 有耦合关系，输出u c 中含有i L 的信息，因此对u c 的检测能确定i L 。即系统能观的。
反之，当4321R R
＝R R 时，电桥中的c 点和d 点的电位始终相等， u c 不受输入u 的控制，u 只
能改变i L 的大小，故系统不能控；由于输出u c 和状态变量i L 没有耦合关系，故u c 的检测不能确定i L ，即系统不能观。
1.1 当4
321R R
R R ≠时
r u ????????+????????????
???????????? ??+-+-???? ??+++-???? ??+-+-???? ??+++-=??????0L 1u i R R 1R R 1C 1R R R R R R R R C 1R R R R R R R R L 1R R R R R R R R L 1u i c L 4321434321214343212
143432121c L (10-1) []??
?
???==c L c u i u y 10 (10-2)
由上式可简写为
bu Ax x
+= cx y = 式中??????=C L u i x ??????
???????????? ??+-+-???? ??+++-???? ??+-+-???? ??+++-=432143432121434321
2
143432121R R 1R R 1C 1R R R R R R R R C 1R R R R R R R R L 1R R R R R R R R L 1A
???
?
????=01L b 1] [0=c
由系统能控能观性判据得
][Ab b rank =2 2=???
???cA c rank
故系统既能控又能观。
1.2 当4
321R R
＝R R 时，式(10-1)变为
r u ????????+????????????
???????????? ??+-+????
??+++-=??????000L 1u i R R 1R R 1C 1-R R R R R R R R L 1u i c L 432143432121c L (10-3) y=u c =[0 1] ???
?
? ??c L u i (10-4)
由系统能控能观性判据得 ][Ab b
rank =1<2 1=??
?
???cA c rank <2
故系统既不能控又不能观，若把式(10-3)展开则有
u L
i R R R R R R R R L i
L L 1)(143432121++++-= (10-5) c c u R R R R C u
)11(14
321+-+-= (10-6) 这是两个独立的方程。第二个方程中的c u 既不受输入u 的控制，也与状态变量L i 没有任何耦合关系，故电路的状态为不能控。同时输出u c 中不含有i L 的信息，因此对u c 的检测不能确定i L ，即系统不能观。
图10-1系统能控性与能观性实验电路图
4；实验结论：
以表格形式记录测得的各项实验数据，并据此分析系统的能控性与能观性。 R3 1V 2V
Uab ucd uab ucd
1K 0.92 0.19 1.68 0.37 2K 0.92 0.24 1.68 0.46 3K 0.92 0.26 1.68 0.52
实验二 控制系统极点的任意配置
1 实验目的：
1．掌握用全状态反馈的设计方法实现控制系统极点的任意配置； 2．用电路模拟的方法，研究参数的变化对系统性能的影响。 2 实验内容：
1．用全状态反馈实现二阶系统极点的任意配置，并用电路模拟的方法予予以实现； 2．用全状态反馈实现三阶系统极点的任意配置，并通过电路模拟的方法予以实现。 3 实验原理：
由于控制系统的动态性能主要取决于它的闭环极点在S 平面上的位置，因而人们常把对系统动态性能的要求转化为一组希望的闭环极点。一个单输入单输出的N 阶系统，如果仅靠系统的输出量进行反馈，显然不能使系统的n 个极点位于所希望的位置。基于一个N 阶系统有N 个状态变量，如果把它们作为系统的反馈信号，则在满足一定的条件下就能实现对系统极点任意配置，这个条件就是系统能控。理论证明，通过状态反馈的系统，其动态性能一定会优于只有输出反馈的系统。
设系统受控系统的动态方程为
bu Ax x
+= cx y =
图11-1为其状态变量图。
图11-1 状态变量图
令Kx r u -=，其中]...[21
n k k k K =，r 为系统的给定量，x 为1?n 系统状态变量，u 为
1?1控制量。则引入状态反馈后系统的状态方程变为
bu x bK A x
+-=)( 相应的特征多项式为
)](det[bK A SI --，调节状态反馈阵K 的元素]...[21
n k k k ，就能实现闭环系统极点的任
意配置。图11-2为引入状态反馈后系统的方框图。
图11-2 引入状态变量后系统的方框图 1．典型二阶系统全状态反馈的极点配置 二阶系统方框图如11-3所示。
图11-3 二阶系统的方框图
1.1 由图得
)
15.0(10
)(+=
S S S G ，然后求得：223.0=ζ，%48≈p δ
同时由框图可得：
211
5.01)
(X S X R =+- ，2110X X = 所以：R X X X 2222
12+--= R X X ??
?
???+????
??--=2022100
[]X X y 011==
1.2 系统能控性 []242200=??
?
???-=rank Ab b
rank
所以系统完全能控，即能实现极点任意配置。 1.3 由性能指标确定希望的闭环极点 令性能指标： 20.0≤p δ，s 5.0Tp ≤
由 20.02
1≤--ζζπ
δe
p ＝，选择707.02
1＝＝
ζ (%3.4=p δ)
s n 5.01Tp 2
≤-=
ζ
ωπ，选择10＝n ω 1/S
于是求得希望的闭环极点为 07.707.721j S ±＝－，
)07.707.72
1
11007.71(2
2,1j j j S n n ±=-
±-=-±-=－ζ
ωξω 希望的闭环特征多项式为
10014.14)07.707.7)(07.707.7()(*2++=++-+=S S j S j S S ? (11-1)
1.4 确定状态反馈系数K1和K2
引入状态反馈后系统的方框图如图11-4所示。
图11-4 引入状态反馈后的二阶系统方框图
其特征方程式为
2
1
222210)(K S K S bK A SI +++-=
--
2020)22(122++++=K S K S (11-2)
由式(11-1)、 (11-2)解得
1.6,421==K K
根据以上计算可知，二阶系统在引入状态反馈前后的理论曲线如图11-5的a)、b)所示。
a) 引入状态反馈前 b) 引入状态反馈后 图11-5 引入状态反馈前后二阶系统的单位阶跃响应曲线 2．典型三阶系统全状态反馈的极点配置 2.1 系统的方框图
三阶系统方框图如11-6所示。
图11-6 三阶系统的方框图
2.2 状态方程
由图得： 21X X = []X X y C 0011===
3
2222X X X +-= R X X X 5553
13+--= 其动态方程为： R X X
????
?
?????+??????????---=500505220010 2.3能控性
由动态方程可得： [
]
3125255701001000
2=??
??
??????--=Rank b A Ab A
Rank
所以系统能控，其极点能任意配置。
设一组理想的极点为：10-=P ，223,2j P ±-= 则由它们组成希望的特征多项式为
804814)22)(22)(10(23+++=++-++=*S S S j s j S S ? (11-3)
2.4 确定状态反馈矩阵K
引入状态反馈后的三阶系统方框图如11-7所示。
图11-7 引入状态反馈后的三阶系统方框图
由图11-7可得
21310)55(2)55)(2()](det[SK K K S S S BK A SI ++++++=--
1322331010)101010()57(K S K K S K S +++++++= (11-4)
由式(11-3)、(11-4)得
14573=+K 4.13=K 4810101032=++K K 4.22=K 8010101=+K 71=K
图11-7对应的模拟电路图如图11-12所示。图中电阻R X1、R X2、R X3按下列关系式确定。
72001=Rx k ，4.22002
=Rx k
，4.12003=Rx k 根据以上计算可知，三阶系统在引入状态反馈前后的理论曲线如图11-8的a)、b)所示。
a) 引入状态反馈前 b) 引入状态反馈后
4 实验结论：
1 根据状态空间表达式，分别画出二阶和三阶系统引入状态反馈后的模拟结构图。
2 画出实验所得曲线图。包括二阶和三阶系统反馈前及反馈后的响应曲线
实验三 具有内部模型的状态反馈控制系统
1 实验目的：
1．通过实验了解内模控制的原理；
2．掌握具有内部模型的状态反馈设计的方法 2 实验内容：
1．不引入内部模型，按要求设计系统的模拟电路，并由实验求取其阶跃响应和稳态输出； 2．设计该系统引入内部模型后系统的模拟电路，并由实验观测其阶跃响应和稳态输出。 3 实验原理：
系统极点任意配置（状态反馈），仅从系统获得满意的动态性能考虑，即系统具有一组希望的闭环极点，但不能实现系统无误差。为此，本实验在上一实验的基础上，增加了系统内部模型控制。
经典控制理论告诉我们，系统的开环传递函数中，若含有某控制信号的极点，则该系统对此输入信号就无稳态误差产生。据此，在具有状态反馈系统的前向通道中引入R(s)的模型，这样，系统既具有理想的动态性能，又有对该系统无稳态误差产生。
1．内模控制实验原理 设受控系统的动态方程为
bu Ax x
+= cx y = 令参考输入为阶跃信号r ，则有： 0=r
令系统的输出与输入间的跟踪误差为
r y e -=，则有：x c r y e
=-= (12-1) 若令 x
Z =，u =ω为两个中间变量，则得 ωb Az u b x A Z
+=+= (12-2) 把(12-1)、(12-2)写成矩阵形式：
ω??
?
???+??????????
?
?=????????b z e A c z e 000.. (12-3) 若(12-3)能控，则可求得如下形式的状态反馈 z K e k 21--=ω ([]32
2k k K =) (12-4)
这不仅使系统稳定，而且实现稳态误差为零。对式(12-4)积分得 ?
--=)()(21t x k dt t e k u
引入参考输入的内部模型后系统的方框图如下图所示：
图12-1 具有内部模型系统的方框图
2．系统的极点配置 已知给定电路的动态方程为
u X X ?
?
?
???+????
??--=101410
，[]X y 01= (12-5) 或写成
21x x
= u x x x
+--=2124 由于Rank[b Ab]=2，故系统能控。 系统的方框图如图12-2所示。
图12-2 引入状态反馈前的二阶系统方框图
由动态方程可得
[]44
411011401)()(21
1++?=??
???????+-?
??=-=--S S S S b A SI C S T
由于12=n ζω，故25.0=ζ 此时超调量较大(约为50%左右)，当单位阶跃输入时，75.0=ss e 。 令状态反馈阵][21
k k K =， Kx r u -=，式中r 为系统的给定量。
则引入状态反馈后的状态方程为
br x bK A x
+-=)( 相应的特征多项式为
12221
4)1(11
4)()(k S k S k S k S
bK A SI S ++++=??
?
++-?
??+=--=ψ (12-6) 设闭环系统的希望极点为j S 112,1±-= 则由它们组成希望的特征多项式为
22)11)(11(2++=++-+=*S S j s j S ? (12-7)
对比式(12-6)和式(12-7)得
21-=k 12=k
此时 []2
22
2110212
01)]([)(2
1
1++?=?????????+-?
??=--=--S S S S
b bK A SI C S T 由于222=ζ
，故ζ 此时超调量约为4.3%，5.0=ss e 。
引入状态反馈后系统的方框图如12-3所示。
图12-3 引入状态反馈后的二阶系统方框图
3．内模控制器的设计
为使校正后的系统不仅具有良好的动态性能，而且要以零稳态误差跟踪输入，因此需在状态反馈的基础上引入内模控制。根据式（12-3）和(12-5)得
ω????
?
?????+????????????????--=??????100140100010z e z e 设闭环系统的希望极点为10,132,1-=±-=S j S ，则得希望的闭环特征方程式为：
202212)10)(1)(1()(*23+++=+++-+=S S S S j S j S S ? (12-8)
引入状态反馈后系统的特征多项式为
??
??
?
?????+++--=--32
1
1410
1det )](1det[k s k k s s
bK A S 12233)4()1(k S k S k S +++++= (12-9) 对比(5)、(6)式得
201=k ，182=k ，113=k 故校正后系统的方框图如图12-4所示。
图12-4 校正后系统的方框图
根据以上计算可知，二阶系统在引入状态反馈前后的理论曲线如图12-5的a)、b)、c)所示。
a) 引入极点配置前 b) 引入极点配置后 c)引入内模控制后
图12-5 内模控制引入前后的阶跃响应曲线
4实验结论：
根据实验结果，画出引入极点配置前、引入极点配置后、引入内模控制后的系统响应曲线。
《大学计算机基础》 上机实验报告 班级： 姓名： 学号： 授课教师： 日期：年月日
目录 一、Windows操作系统基本操作......................................................... - 1 - 二、Word文字处理基本操作 .............................................................. - 4 - 三、Excel电子表格基本操作 ............................................................ - 6 - 四、PowerPoint幻灯片基本操作....................................................... - 8 - 五、网页设计基本操作 ...................................................................... - 9 - 六、Access数据库基本操作 ............................................................ - 10 - 上机实验作业要求： ○1在实验报告纸上手写并粘贴实验结果； ○2每人将所有作业装订在一起（要包封面）； ○3全部上机实验结束后全班统一上交； ○4作业内容不得重复、输入的数据需要有差别。
实验名称一、Windows操作系统基本操作 实验目的1、掌握Windows的基本操作方法。 2、学会使用“画图”和PrntScr快捷键。 3、学会使用“计算器”和Word基本操作。 实验内容1、日历标注 利用“画图”和Word软件，截取计算机上日历的图片并用文字、颜色、图框等标注出近期的节假日及其名称，并将结果显示保存在下面（参考下面样图）。 运行结果是： 主要操作步骤是： 2、科学计算 利用“计算器”和Word软件，计算下列题目，并将结果截图保存在下面（参考样图）。 ○1使用科学型计算器，求8!、sin(8)、90、74、20、67、39、400、50.23、ln(785)的平均值、和值，并用科学计数法显示。 运行结果是： ②将以下十、八、十六进制数转换为二进制数：(894.8125)10、(37.5)8、(2C.4B)16 运行结果是：（需要下载使用“唯美计算器”） ○3计算下列二进制数的加法与乘法：101.1+11.11；1101*1011 运行结果是：（参考样图） 写出主要操作步骤： 3、实验心得体会
实验报告 ( 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
即 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 --------?? -+-+=??--??
第 1 页 共 1 页 西 安 科 技 大 学2004—2005 学 年 第2 学 期 期 末 考 试 试 题（卷） 电控 院系： 班级： 姓名： 学号： 装 订 线 装 订 线 以 内 不 准 作 任 何 标 记 装 订 线 10x =?-?
第 2 页 共 1 页 现代控制理论C 卷答案 2 解：Bode 图略 解得：开环截止频率：)/(1.2s rad c =ω； 相角裕量：)(40rad r ≈ 3 解： 1）系统的传递函数阵为： 2231231))((1 ))()((1 ][)(du a s a s a s a s a s Du B A sI C s G +?? ?? ??-----=+-=- 2）系统的状态结构图，现以图中标记的321,,x x x 为 u 2u 1
第 3 页 共 1 页 4解： 1）列写电枢电压u 为输入，以电流i 和旋转速度n 为输出的状态空间表达式。由于ω.πωn 559260==，可得 dt dn J dt d J 55.9=ω， 22)2(D g G mR J == 式中， m 为一个旋转体上的一个质点的质量，质量m 为该质量的重量G 和重力加速度g 之比，R 和D 分别为旋转体的半径和直径，综合上两 式可推得 dt dn GD dt dn D G dt d J 37548.955.922=???=ω 2）从而可得到电机电枢回路电压平衡和电机运动平衡的一组微分方程式 ???????=+=++i C n K dt dn GD u n C Ri dt di L m b e 3752 式中，摩擦系数55.9/B K b =。选择状态变量n x i x ==21,，则系统得状态空间表达式为 u L x x GD K GD C L C L R x x b m e ????????+??????????? ????? -- -=? ?????01375375212 2 21 ????????????=211001x x y
华北电力大学 实验报告| | 实验名称状态空间模型分析 课程名称现代控制理论 | | 专业班级：自动化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 为状态初值。
现代控制理论实验报告
实验一系统能控性与能观性分析 一、实验目的 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的关系 Uab Ucd
现代控制理论基础考试题 西北工业大学考试题（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. 对于任意的初始状态)(0t x 和输入向量)(t u ，系统状态方程的解存在并且 惟 一 。 （√） 5. 传递函数矩阵也能描述系统方程中能控不能观测部分的特性。 （×）
《大学计算机基础Ⅰ》课程 实验报告手册 \ 实验教师（签字） 西南大学计算机与信息科学学院 计算机基础教育系 年月日
一、实验说明 本课程实验分为一般性实验（验证和简单设计）和综合性实验（课程设计）两部分。从第3周开始参考实验任务书（本报告中的五部分）完成每周规定的实验，并根据进度按要求认真填写本实验报告中的六、七部分，此实验报告将作为实验成绩评定的依据之一。 本课程实验从开课学期第3周开始实习，每周2学时，16周结束，共28学时。除统一安排的时间外，学生还可根据自己的实际适当安排课余时间上机。上机内容参见本报告中的“五、实验任务书”部分。 二、实验目的 通过本实验，让学生掌握计算机的基本操作和基本技能，能够学会知识的运用与积累，能够举一反三，具备一定的独立解决问题的能力和信心，培养学生熟练地使用常用软件的能力及严肃认真的科学作风，为今后的学习和工作打下良好的基础。 三、实验要求 1、每次实验课将考勤，并作为实验成绩的重要依据。 2、每次实验前学生必须充分准备每次的实验内容，以保证每次上机实验的效果。实验过程中必须独立完成。 3、学期结束时，每位同学应将自己的《实验报告》交各专业班长或学习委员，由班长或学习委员以专业为单位、按学号从小到大排列好统一交给实验指导老师，否则无实验成绩。 四、实验报告要求 一共要求填写3个阶段性实验报告、1个综合性实验报告和1份学期总结，与每份实验报告对应产生的电子文档交由实验老师指定的位置，该电子文档也将作为实验成绩评定的依据之一。 五、实验任务书 教材：《大学计算机基础》第五版高等教育出版社 实验参考书：《大学计算机基础实践教程》高等教育出版社 实验一：指法练习、汉字录入 实验目的： 1．掌握鼠标和键盘的使用及正确的操作指法。 2．掌握微型计算机的打开和关闭操作 3．熟悉键盘指法和文字录入 4．了解中英文切换，全半角的切换 实验任务： 1．参见实验参考书中的实验1-1-1中的[任务1]（7页） 2．参见实验参考书中的实验1-1-1中的[任务3]（7页） 实验二：Windows的基本操作和文件管理操作 实验目的： 1．掌握Windows的基本知识和基本操作 2．掌握“Windows资源管理器”和“我的电脑”的使用 实验任务： 1．参见实验参考书中的实验1-2-1中的全部任务（14页） 2．参见实验参考书中的实验1-2-2中的全部任务（18页）
河南工业大学 现代控制理论实验报告姓名：朱建勇 班级：自动1306 学号：201323020601
现代控制理论 实验报告 专业: 自动化 班级: 自动1306 姓名: 朱建勇 学号: 201323020601 成绩评定: 一、实验题目： 线性系统状态空间表达式的建立以及线性变换 二、实验目的 1. 掌握线性定常系统的状态空间表达式。学会在MATLAB 中建立状态空间模型的方法。 2. 掌握传递函数与状态空间表达式之间相互转换的方法。学会用MATLAB 实现不同模型之 间的相互转换。 3. 熟悉系统的连接。学会用MATLAB 确定整个系统的状态空间表达式和传递函数。 4. 掌握状态空间表达式的相似变换。掌握将状态空间表达式转换为对角标准型、约当标准 型、能控标准型和能观测标准型的方法。学会用MATLAB 进行线性变换。 三、实验仪器 个人笔记本电脑 Matlab R2014a 软件 四、实验内容 1. 已知系统的传递函数 (a) ) 3()1(4)(2++=s s s s G
(b) 3486)(22++++=s s s s s G
(c) 6 1161)(232+++++=z z z z z z G （1）建立系统的TF 或ZPK 模型。 （2）将给定传递函数用函数ss( )转换为状态空间表达式。再将得到的状态空间表达式用函 数tf( )转换为传递函数，并与原传递函数进行比较。 （3）将给定传递函数用函数jordants( )转换为对角标准型或约当标准型。再将得到的对角 标准型或约当标准型用函数tf( )转换为传递函数，并与原传递函数进行比较。 （4）将给定传递函数用函数ctrlts( )转换为能控标准型和能观测标准型。再将得到的能控标 准型和能观测标准型用函数tf( )转换为传递函数，并与原传递函数进行比较。
现代控制理论课程总结 学习心得 从经典控制论发展到现代控制论,是人类对控制技术认识上的一次飞跃。现代控制论是用状态空间方法表示,概念抽象,不易掌握。对于《现代控制理论》这门课程，在刚拿到课本的时候，没上张老师的课之前，咋一看，会认为开课的内容会是上学期学的控制理论基础的累赘或者简单的重复，更甚至我还以为是线性代数的复现呢！根本没有和现代控制论联系到一起。但后面随着老师讲课的风格的深入浅出，循循善诱，发现和自己想象的恰恰相反，张老师以她特有的讲课风格，精心准备的ppt 课件，向我们展示了现代控制理论发展过程，以及该掌握内容的方方面面，个人觉得，我们不仅掌握了现代控制理论的理论知识，更重要的是学会了掌握这门知识的严谨的逻辑思维和科学的学习方法，对以后学习其他知识及在工作上的需要大有裨益，总之学习了这门课让我受益匪浅。 由于我们学习这门课的课时不是很多，并结合我们学生学习的需求及所要掌握的课程深入程度，张老师根据我们教学安排需要，我们这学期学习的内容主要有：1.绪论；2.控制系统的状态表达式；3.控制系统状态表达式的解；4.线性系统的能空性和能观性；5.线性定常系统的综合。而状态变量和状态空间表达式、状态转移矩阵、系统的能控性与能观性以及线性定常系统的综合是本门课程的主要学习内容。当然学习的内容还包括老师根据多年教学经验及对该学科的研究的一些深入见解。 在现代科学技术飞速发展中,伴随着学科的高度分化和高度综合,各学科之间相互交叉、相互渗透,出现了横向科学。作为跨接于自然科学和社会科学的具有横向科学特点的现代控制理论已成为我国理工科大学高年级的必修课。 经典控制理论的特点 经典控制理论以拉氏变换为数学工具，以单输入－单输出的线性定常系统为主要的研究对象。将描述系统的微分方程或差分方程变换到复数域中，得到系统的传递函数，并以此作为基础在频率域中对系统进行分析和设计，确定控制器的结构和参数。通常是采用反馈控制，构成所谓闭环控制系统。经典控制理论具有明显的局限性，突出的是难以有效地应用于时变系统、多变量系统，也难以揭示系统更为深刻的特性。当把这种理论推广到更为复杂的系统时，经典控制理论就显得无能为力了，这是因为它的以下几个特点所决定。 1．经典控制理论只限于研究线性定常系统，即使对最简单的非线性系统也是无法处理的；这就从本质上忽略了系统结构的内在特性，也不能处理输入和输出皆大于1的系统。实际上，大多数工程对象都是多输入－多输出系统，尽管人们做了很多尝试，但是，用经典控制理论设计这类系统都没有得到满意的结果；2．经典控制理论采用试探法设计系统。即根据经验选用合适的、简单的、工程上易于实现的控制器，然后对系统进行分析，直至找到满意的结果为止。虽然这种设计方法具有实用等很多完整，从而促使现代控制理论的发展：对经典理论的精确化、数学化及理论化。优点，但是，在推理上却是不能令人满意的，效果也
一．（本题满分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 +=+ 所以状态空间表达式为
实验三 利用MATLAB 导出连续状态空间模型的离散化模型 实验目的： 1、基于对象的一个连续时间状态空间模型，导出其相应的离散化状态空间模型； 2、通过编程、上机调试，掌握离散系统运动分析方法。 实验原理： 给定一个连续时间系统的状态空间模型： ()()()()()() x t Ax t Bu t y t Cx t Du t =+=+ （3.1） 状态空间模型（3.1）的输入信号()u t 具有以下特性： ()(),u t u kT kT t kT T =≤≤+ （3.2） 已知第k 个采样时刻的状态()x kT 和第k 个采样时刻到第1k +个采样时刻间的输入()()u t u kT =，可得第1k +个采样时刻(1)k T +处的状态 (1)((1))((1))()((1))()k T kT x k T k T kT x kT k T Bu d τττ++=Φ+-+Φ+-? （3.3） 其中： ((1))((1))A k T kT AT k T kT e e +-Φ+-== ((1))((1))A k T k T e ττ+-Φ+-= 由于输入信号在两个采样时刻之间都取常值，故对式（3.3）中的积分式进行一个时间变量替换(1)k T στ=+-后，可得 0((1))()()()AT A x k T e x kT e d Bu kT τ σσ+=+? （3.4） 另一方面，以周期T 对输出方程进行采样，得到 ()()()y kT Cx kT Du kT =+ 在周期采样的情况下，用k 来表示第k 个采样时刻kT 。因此，连续时间状态空间模型
（3.1）的离散化方程可以写成 (1)()()()()()()() x k G T x k H T u k y k Cx k Du k +=+=+ （3.5） 其中： 0()()()AT A G T e H T e d B τσσ==? （3.6） 已知系统的连续时间状态空间模型，MATLAB 提供了计算离散化状态空间模型中状态矩阵和输入矩阵的函数： [G ,H]=c2d(A,B,T) 其中的T 是离散化模型的采样周期。 实验步骤 1、导出连续状态空间模型的离散化模型，采用MA TLAB 的m-文件编程； 2、在MA TLAB 界面下调试程序，并检查是否运行正确。 例3.1 已知一个连续系统的状态方程是 010()()()2541x t x t u t ????=+????--???? 若取采样周期0.05T =秒，试求相应的离散化状态空间模型。 编写和执行以下的m-文件： A=[0 1;-25 –4]; B=[0;1]; [G ,H]=c2d(A,B,0.05) 得到 G= 0.9709 0.0448 -1.1212 0.7915 H= 0.0012 0.0448 因此，所求的离散化状态空间模型是 0.97090.04480.0012(1)()()1.12120.79150.0448x k x k u k ????+=+????-????
现代控制理论实验报告 组员： 院系：信息工程学院 专业： 指导老师： 年月日
实验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
实验六利用MATLAB设计状态观测器 ******* 学号 1121*****
实验目的： 1、学习观测器设计算法； 2、通过编程、上机调试，掌握基于观测器的输出反馈控制系统设计方法。 实验原理： 1、全阶观测器模型： () ()x Ax Bu L y Cx A LC x Bu Ly =++-=-++ 由极点配置和观测器设计问题的对偶关系，也可以应用MATLAB 中极点配置的函数来确定所需要的观测器增益矩阵。例如，对于单输入单输出系统，观测器的增益矩阵可以由函数 L=(acker(A ’,C ’,V))’ 得到。其中的V 是由期望的观测器极点所构成的向量。类似的，也可以用 L=(place(A ’,C ’,V))’ 来确定一般系统的观测器矩阵，但这里要求V 不包含相同的极点。 2、降阶观测器模型： ???w Aw By Fu =++ b x w Ly =+ 基于降阶观测器的输出反馈控制器是： ????()[()]()b a b b a b w A FK w B F K K L y u K w K K L y =-+-+=--+ 对于降阶观测器的设计，使用MATLAB 软件中的函数 L=(acker(Abb’,Aab’,V))’ 或 L=(place(Abb’,Aab’,V))’ 可以得到观测器的增益矩阵L 。其中的V 是由降阶观测器的期望极点所组成的向量。 实验要求 1．在运行以上例程序的基础上，考虑图6.3所示的调节器系统，试针对被控对象设计基于全阶观测器和降 阶观测器的输出反馈控制器。设极点配置部分希望的闭环极点是1,22j λ=-± （a ） 对于全阶观测器，1 8μ=-和 28μ=-； （b ） 对于降阶观测器，8μ=-。 比较系统对下列指定初始条件的响应： （a ） 对于全阶观测器： 1212(0)1,(0)0,(0)1,(0)0x x e e ==== （b ） 对于降阶观测器： 121(0)1,(0)0,(0)1x x e === 进一步比较两个系统的带宽。
《现代控制理论》实验指导书 俞立徐建明编 浙江工业大学信息工程学院 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 &&&
. 现代控制理论实验报告 组员： 院系：信息工程学院 专业： 指导老师： 年月日
实验1 系统的传递函数阵和状态空间表达式的转换 [实验要求] 应用MATLAB 对系统仿照[例1.2]编程，求系统的A 、B 、C 、阵；然后再仿照[例1.3]进行验证。并写出实验报告。 [实验目的] 1、学习多变量系统状态空间表达式的建立方法、了解系统状态空间表达式与传递函数相互转换的方法； 2、通过编程、上机调试，掌握多变量系统状态空间表达式与传递函数相互转换方法。 [实验内容] 1 设系统的模型如式(1.1)示。 p m n R y R u R x D Cx y Bu Ax x ∈∈∈?? ?+=+= (1.1) 其中A 为n ×n 维系数矩阵、B 为n ×m 维输入矩阵 C 为p ×n 维输出矩阵，D 为传递阵，一般情况下为0，只有n 和m 维数相同时，D=1。系统的传递函数阵和状态空间表达式之间的关系如式(1.2)示。 D B A SI C s den s num s G +-== -1)() () (()( (1.2) 式(1.2)中，)(s num 表示传递函数阵的分子阵，其维数是p ×m ；)(s den 表示传递函数阵的按s 降幂排列的分母。 2 实验步骤 ① 根据所给系统的传递函数或（A 、B 、C 阵），依据系统的传递函数阵和状态空间表达式之间的关系如式(1.2)，采用MATLA 的file.m 编程。注意：ss2tf 和tf2ss 是互为逆转换的指令； ② 在MATLA 界面下调试程序，并检查是否运行正确。 ③ [1.1] 已知SISO 系统的状态空间表达式为(1.3)，求系统的传递函数。
现代控制理论基础课程总结 学院:__机械与车辆学院_ 学号:____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])