弹簧质量阻尼模型Word文档下载推荐.docx
《弹簧质量阻尼模型Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《弹簧质量阻尼模型Word文档下载推荐.docx（12页珍藏版）》请在冰豆网上搜索。
设缓冲器的摩擦力与活塞的速度成正比，其中
=1kg，
=2kg，
=
=100N/cm，
=300N/cm，
=3N
s/cm，
=6N
s/cm。
由图2.1，根据牛顿第二定律,，建立系统的动力学模型如下：
对
有：
（2-1）
（2-2）
3建立状态空间表达式
令
则原式可化为：
化简得：
（2-3）
（2-4）
整理得：
（2-5）
代入数据得：
则系统的状态空间表达式为
4化为对角标准型
当系统矩阵A有n个不相等的特征根
时，相应的有n个不相等的特征向量
，所以有矩阵A的特征矩阵
根据矩阵论
线性变换得：
可以使用matlab进行对角标准型的运算，matlab作为一种数学运算工具，很大程度的方便了了我们的计算，对于这个弹簧-质量-阻尼系统是一个四阶的状态空间表达式，所以可以用matlab简化计算。
（1）求特征值与特征向量
A=[0010;
0001;
-40030096;
150-2003-4.5]
B=[00;
00;
10;
00.5]
C=[1000;
0100]
[P,J]=eig（A）
求得结果：
P=
0.0007-0.0402i0.0007+0.0402i0.0401-0.0698i0.0401+0.0698i
-0.0171+0.0157i-0.0171-0.0157i0.0176-0.0792i0.0176+0.0792i
0.86500.86500.6682+0.2084i0.6682-0.2084i
-0.3442-0.3621i-0.3442+0.3621i0.70500.7050
J=
0.3667+21.5183i000
00.3667-21.5183i00
001.8833+8.4864i0
0001.8833-8.4864i
（2）P矩阵求逆
PN=inv（P）
PN=
3.4167+9.7803i-2.1017-9.2399i0.3466-0.2323i-0.4703-0.1054i
3.4167-9.7803i-2.1017+9.2399i0.3466+0.2323i-0.4703+0.1054i
-3.3554+3.4224i3.7199+3.2032i0.2886-0.0353i0.5337-0.2409i
-3.3554-3.4224i3.7199-3.2032i0.2886+0.0353i0.5337+0.2409i
（3）带入公式
解得对角标准型为：
5求状态空间表达式的解
（1）求状态转移矩阵
其中，T为特征向量
状态转移矩阵为：
5可控性与可观性
不同于经典控制理论，能控性和能观性，是一个具有实际意义的概念，经典控制理论中用传递函数描述系统的输入-输出特性，输出量即被控量，只要系统是因果系统并且稳定，输出量便可以受控，且输出量总是可以被测量的，因而不需要能控能观性的提出。
但是现代控制理论是建立在状态空间表达式描述系统的基础上的，状态方程描述输入u（t）引起状态x（t）的变化过程，输出方程描述有状态变化引起的输出y（t）的变化。
能控能观便是定性的描述输入u（t）对状态x（t）的控制能力，输出y（t）对状态x（t）的反应能力，他们分别回答了“输入能否控制状态的变化”------可控性
“状态的变化能否有输出反映出来”----------可观性
另外在工程上常用状态变量作为反馈信息，可是状态x（t）的值通常是难测的，往往需要从测量到的y（t）中估计出状态，如果输出y（t）不能完全反映出系统的状态x（t），那么就无法实现对状态的估计。
能控性定义：
当系统用状态方程描述时，给定系统的任意初始状态，可以找到允许的输入量，在有限的时间内使系统的所有状态达到任一终止状态，则称系统是完全可控的。
有状态方程x’（t）=Ax（t）+Bu（t）其解为：
如果有限的时间内0<
t<
t1内通过输入量u（t）的作用把系统的所有状态引向状态x（t1）[设x（t1）=0]，则应有：
即在给定x（0-）和A、B的条件下求可以使x（t）=x（t1）的u（t）。
换言之：
上述方程有解则系统能控。
根据凯莱－哈米尔顿定理，e-At、eAt可写成有限级数：
如果方程有解，等式右边左侧矩阵应满秩=n，此时系统是可控的。
求可控性：
n=4满秩所以系统是可控的
可观性定义：
当系统用状态方程描述时，给定控制后，如果系统的每一个初始状态x（0-）都可以在有限的时间内通过系统的输出y（t）唯一确定，则称系统完全可观。
若只能确定部分初始状态，则称系统部分可观。
有状态方程x’（t）=Ax（t）+Bu（t）y（t）=Cx（t）其解为
由于在讨论能观性问题时，输入是给定的，上式右侧第二项是确知的，设u（t）=0。
y（t）=CeAtx（0-）。
根据凯莱－哈米尔顿定理，e-At、eAt可写成有限级数：
如果方程有解，等式右侧中间侧矩阵应满秩。
其中，秩=n（系统的阶数）
求可观性：
n=4满秩所以系统是可观的
6求系统的输入输出传递函数
对于两输入两输出的系统求得的传递函数是一个二阶的传递函数阵,其中包含四个传递函数
Transferfunctionfrominput1tooutput...
s^2+4.5s+200
#1:
--------------------------------------------
s^4-4.5s^3+541.5s^2-1800s+3.5e004
3s+150
#2:
Transferfunctionfrominput2tooutput...
0.5s^2-4.5s+200
矩阵函数阵：
7分析开环稳定性
稳定性定义是系统在受到小的外界扰动后，系统状态方程解的收敛性。
系统正常工作要求是系统在受到外界扰动后，虽然其原平衡状态被打破，但在扰动消失后，仍然能恢复到原来的平衡状态，或者趋于另一平衡状态继续工作，且线性系统稳定性与输入作用无关。
研究系统的稳定性对于研究系统能否正常工作具有很重要的意义，稳定性是自动控制系统正常工作的必要条件，是系统的重要特征。
我们不仅要分析一个系统是否稳定还要解决的问题便是怎样使一个系统稳定。
经典控制理论稳定性判别方法有很多，例如代数判据，niquist判据，根轨迹判据等。
而现代控制理论经常用李雅普诺夫第二法求稳定性。
（1）利用特征根的方法
根据上述结果求得的特征根为0.3667+21.5183i0.3667+21.5183i1.8833+8.4864i1.8833-8.4864i，四个特征值全部都在坐标轴的右半平面，所以系统是不稳定的。
（2）利用利亚普诺夫第二法求解
其中：
将其换成矩阵形式
可以看出A不是正定的，所以系统不稳定。
8利用状态反馈将系统的闭环极点配置到合适的值
状态反馈是将系统每一个状态变量乘以相应的反馈系数馈送到输入端与参考输入相加，其和作为受控系统的控制输入。
原受控对象为经过状态反馈后得到的闭环系统为
闭环系统期望极点选取原则为以下几点：
1）n维控制系统有n个期望极点；
2）期望极点是物理上可实现的，为实数或共轭复数对；
3）期望极点的位置的选取，需考虑它们对系统品质的影响（离虚轴的位置），及与零点分布状况的关系。
4）离虚轴距离较近的主导极点收敛慢，对系统性能影响最大，远极点收敛快，对系统只有极小的影响。
闭环极点0.3667+21.5183i0.3667-21.5183i1.8833+8.4864i1.8833-8.4864i
配置状态反馈后，系统应稳定，所以期望极点应在虚轴左侧，所以期望闭环极点-1.1001+64.5549i-1.1001-64.5549i-5.6499+25.4952-5.6499-25.4952
得到极点配置矩阵
K=
1.0e+003*
1.2556-0.03750.0157-0.0332
0.97182.89690.08390.004
验证极点配置结果是正确的：
ans=
-1.1001-64.5549i-1.1001+64.5549i-5.6499-25.4952i-5.6499+25.4952i
求得开环传递函数阶跃响应曲线（没有经过状态反馈的）：
没有上升时间
经过状态反馈的传递函数：
状态空间表达式为
Matlab解得闭环传递函数：
s^2+6.758s+1648
--------------------------------------------------
s^4+13.5s^3+4875s^2+4.86e004s+2.843e006
-38.95s-335.9
19.59s+168.7
0.5s^2+3.371s+827.8
反馈后的阶跃响应：
阶跃响应上升时间是0.034s，配置后系统最终稳定。
9设计全维状态观测器
不是所有的系统状态物理上都能够直接测量得到。
需要从系统的可量测参量，如输入u和输出y来估计系统状态。
状态观测器基于可直接量测的输出变量y和控制变量u来估计状态变量，是一个物理可实现的模拟动力学系统。
如果系统状态是完全能观测的，那么根据输出y的测量，可以唯一的确定系统的初始状态，所以只要满足一定的条件，便可以从可测量y和u中把x间接重构出来。
全维渐进状态结构图：
系统的原极点为0.3667+21.5183i0.3667+21.5183i1.8833+8.4864i1.8833-8.4864i
期望极点应是原极点的2-5倍，并不是越快越好。
期望极点为-1.4668+86.0732i-1.4668-86.0732i-7.5332+33.9456i-7.5332-33.9456i
相应的全维观测器是：
10带观测器的输出状态空间表达式：
D=0
分别得到输出和观测状态的传递函数：
Transferfunctionfrominput"
y1"
tooutput...
17.51s^3+5899s^2+1.23e005s+7.911e006
y1_e:
------------------------------------------------
s^4+18s^3+8664s^2+1.152e005s+8.96e006
-52.62s^3+359.1s^2-6844s+4.369e005
y2_e:
x1_e:
x2_e:
3200s^3+9.295e004s^2+7.698e006s+3.967e006
x3_e:
122.3s^3+1.405e005s^2-8571s+5.737e007
x4_e:
y2"
51.13s^3+1379s^2+2.945e004s+8.583e005
4.985s^3+5001s^2+1.152e005s+8.303e006
1609s^3-1.337e005s^2+1.373e006s-1.853e007
2246s^3+2.437e004s^2+6.754e006s-4.844e007
Inputgroups:
NameChannels
Measurement1,2
Outputgroups:
OutputEstimate1,2
StateEstimate3,4,5,6
分析其稳定性
根据阶跃响应图来看系统是稳定的，并且上升时间为0.0015s
附录：
D=0;
G=ss（A,B,C,D）
Gc=tf（G）
Po=[-1.1001+64.5549i-1.1001-64.5549i-5.6499+25.4952i-5.6499-25.4952i]
K=place（A,B,Po）
a=A-B*K
g=ss（a,B,C,D）
gc=tf（g）
eig（a）'
P1=[-1.4668+86.0732i-1.4668-86.0732i-7.5332+33.9456i-7.5332-33.9456i]
L=place（A'
C'
P1）'
est=estim（G,L）
step（est）
tf（est）
结束语：
这次大作业是把正本现代控制理论书的内容串联起来，现代控制理论是从状态空间表达式出发研究一个系统的稳定性和控制品质，而自动控制原理则是从输入输出出发来研究一个系统的，这份大作业的最终目的是要设计一个状态观测器去改变系统的控制品质，通过这次大作业，我锻炼了自己查找信息的能力，无论是从网上或者是书籍上，培养了自己自主学习的能力和做一项任务的能力，我通过去图书馆查阅相关的书籍，上搜索引擎去搜索相关的内容，解答做大作业是的困惑。
还有就是学习了matlab软件的使用。
判断系统的能控能观性以及稳定性，状态空间表达式和传递函数之间的化简，都是为世纪一个状态观测器所服务的。
总之，通过这次大作业，收获不少的，不同于传统考核，让我书上所学和知识和实践相结合的机会，真正锻炼了自己各方面的能力，提高自己的实践能力。
华北电力大学
升级会员 