现代控制一级倒立摆.docx
- 文档编号:7603604
- 上传时间:2023-01-25
- 格式:DOCX
- 页数:10
- 大小:177.11KB
现代控制一级倒立摆.docx
《现代控制一级倒立摆.docx》由会员分享,可在线阅读,更多相关《现代控制一级倒立摆.docx(10页珍藏版)》请在冰豆网上搜索。
现代控制一级倒立摆
现代控制一级倒立摆
倒立摆实验
电子工程学院
自动化
学号:
1实验设备简介4
1.1倒立摆介绍4
1.2直线一级倒立摆5
2倒立摆建模6
2.1直线一阶倒立摆数学模型的推导6
2.1.1受力分析6
2.1.2微分方程建模8
2.1.3状态空间数学模型9
2.2实际系统模型建立10
3系统定性、定量分析11
3.1系统稳定性与可控性分析11
3.1.1稳定性分析11
3.1.2能控性分析.13
4极点配置的设计步骤13
4.1极点配置的计算13
4.2用MATLAB讲行极点配置的计算
15
4.3极点配置的综合分析16
5小结.17
1实验设备简介
1.1倒立摆介绍
倒立摆是处于倒置不稳定状态,人为控制使其处于动态平衡的一种摆。
如杂技演员顶杆的物理机制可简化为一级倒立摆系统,是一个
复杂,多变量,存在严重非线性,非自制不稳定系统。
常见的倒立摆一般由小车和摆杆两部分组成,其中摆杆可能是一级,二级或多级,在复杂的倒立摆系统中,摆杆的长度和质量均可变化。
1.2直线一级倒立摆
根据自控原理实验书上相关资料,直线一级倒立摆在建模时,一
般忽略系统中的一些次要因素.例如空气阻力、伺服电机的静摩擦力、系统连接处的松弛程度等,之后可将直线一级倒立摆系统抽象成小车和匀质的杆组成的系统。
倒立摆系统是典型的机电一体化系统其机械部分遵循牛顿的力学定律其电气部分遵守电磁学的基本定理.无论哪
种类型的倒立摆系统,都具有3个特性,即:
不确定性、耦合性、开环不稳定性.直线型倒立摆系统,是由沿直线导轨运动的小车以及一端固定于小车上的匀质长杆组成的系统.小车可以通过传动装置由交流伺服电机驱动•小车导轨一般有固定的行程,因而小车的运动范围是
受到限制的
2倒立摆建模
2.1直线一阶倒立摆数学模型的推导
对于忽略各种摩擦参数和空气阻力之后,直线一即倒立摆抽象为小车和均质杆组成的系统。
本系统的参数定义如下:
M
小车质量
m
摆杆质量
b
小车摩擦系数
l
摆杆转动轴心到杆质心的长度
I
摆杆惯量
F
加在小车上的力
x
小车位置(变量)
©
摆杆与垂直向上方向的夹角(输出)
0
摆杆与垂直向下方向的夹角(考虑到摆杆初始位置为竖直向下)
2.1.1受力分析
F面我们对这个系统作一下受力分析。
N和p为小车与摆杆相互作用
力的水平和垂直方向的分量
图1:
倒立摆系统小车和摆杆的受力分析
应用牛顿第二定律方法来建立系统的动力学方程过程如下:
分析小车水平方向所受的合力,可以得到以下方程:
(2-1)
Mx=F-bx-N
由摆杆水平方向的受力进行分析可以得到下面等式:
把这个等式代入上式中,就得到系统的第一个运动方程:
(2-4)
(Mm)xbxmhcos^-mh2sin-F
为了推出系统的第二个运动方程,我们对摆杆垂直方向上的合力进行
分析,可以得到下面方程:
力矩平衡方程如下:
此方程中力矩的方向,由于
合并这两个方程,约去P和N,得到第二个运动方程:
2・・
(2-8)
(Iml户mglsin^二-mlxcos-
2.1.2微分方程建模
设"=7:
'",当摆杆与垂直向上方向之间的夹角••与1(单位是弧度)
相比很小,即V1时,则可以进行近似处理:
cos=…1,Sin八-,
量,用u来代表被控对象的输入力F,线性化后得到该系统数学模型
的微分方程表达式:
<2・・■■
(I+ml3-mg忡=mix
..
(2-9)
(Mm)xbx-mi二u
2.1.3状态空间数学模型
由现代控制理论原理可知,控制系统的状态空间方程可写成如下形
式:
X二AXBu
Y二CXDu(2-10)
方程组(2-9)对X;•解代数方程,得到如下解:
一(|m|2)b2X—(Im|2)2U
I(Mm)MmlI(Mm)MmlI(Mm)Mml
整理后得到系统状态空间方程:
2.2实际系统模型建立
实际系统参数如下,求系统的传递函数、状态空间方程,并进行脉冲响应和阶跃响应的Matlab仿真。
M
小车质量
1.096kg
m
摆杆质量
0.109kg
b
小车摩擦系数
0.1N/m/sec
l
摆杆转动轴心到杆质心
的长度(本实验为0.3m)
0.25m
I
摆杆惯量
0.006kg*m*m
F
加在小车上的力
x
小车位置
T
米样频率
0.005秒
0
摆杆与垂直向下方向的
夹角
1)以外界作用力作为输入的系统状态方程:
xlo
1
0
0[刃
■0]
x
0
-0.0891
0.7167
0
x
0.8906
=
0
0
0
1
+
u
0
k
1
0
-1.9235
31.6926
0一
1
2.6838一
■xl
。
]x+『]
0一©[o一2」
3系统定性、定量分析
3.1系统稳定性与可控性分析3.1.1稳定性分析
先分析系统的稳定性,将数据代入状态方程中,利用matlab程序可
以求出系统的零极点。
程序段如下:
[num,den]=ss2tf(A,B,C,D)%状态空间表达式>传递函数
为-2.6838/(sT-31.6926)
结果如下:
零^点.
二^八、、•
z=
5.43440.5492
-5.4344-0.0000
极点:
P=
0
-5.6516
5.6081
-0.0456
由得到的p(极点)可知,有的极点在单位圆外,所以可知原系统是不稳定。
3.1.2能控性分析
我们可以利用matlab来得到系统的能控性,源代码如下:
uc二ctrb(A,B)%判断能控性
r=rank(uc)滅秩
r=4
由得到的rank(ud)的值可知,原系统的能控性矩阵为4,所以我们
可知原系统是能控的。
4极点配置的设计步骤
4.1极点配置的计算
对于如上所述的系统,设计控制器,要求系统具有较短的调整时间和
合适的阻尼。
倒立摆极点配置原理图如图所示:
®倒|立拦is貞翫蛊原陛[Z]
极点配置步骤如下:
(1)检验系统能控性(以证)
(2)计算特征值
选取期望的闭环极点,=「(i=1,2,3,4),其中:
S=-10=-10%=-23.3jS=-2-3.3j
期望的特征方程为:
-亠」2-叮=10102-3.3j23.3j
=4243194.892697.81489
因此可以得到:
S=1489
口2=697.8
«3=194,89
a4=24
由系统的特征方程:
扎
-1
00
0
人+0.0891
-0.71670
432
人1-A
=
0
0
丸-1
=k—0.089仏-31.6926k-2.8238扎
0
1.9235
-31.6926九
因此有:
a〔=0a2=2.8238a3=-31.6926a^■-0.0891
系统的反馈增益矩阵为:
K=bi-"a2_〉2a3-3a4-「4T'
(3)状态反馈增益矩阵K为:
K=bi-a2_〉2a3-3a4-4T'
珂-56.6141-32.9817114.117719.8540]
程序段如下:
[num,den]=ss2tf(A,B,C,D)
p二[-10-10-2+3.3j-2-3.3j];
K=acker(A,B,p)
4.2用MATLA进行极点配置的计算
%犬态空间表达式——>传递函
数为-2.6838/(sA2-31.6926)
%期望极点配置
%求状态增益K
结果相同:
K=[-56.6141-32.9817114.117719.8540]
系统输出响应
图3系统输出响应图
4.3极点配置的综合分析
极点配置法成功实现了同时对倒立摆摆角和小车的位置的控制,但是
在极点配置时,期望极点的选取,需要考虑、研究它们对系统品质的影响以及它们与零点分布状况的关系,还需要顾及抗干扰性能方面的要求。
状态反馈系统的主要优点是极点的任意配置,无论开环极点和零点在什么位置,都可以任意配置期望的闭环极点。
这为我们提供了控制系统的手段,假如系统的所有状态都可以被测量和反馈的话,状态反馈可以提供简单而适用的设计。
5小结
通过本次读书工程,我对状态空间有了更深的了解,从建立系统传递函数到状态空间方程再到对系统的稳定性能控性分析,都充分的将理论与实践结合了起来。
极点配置法做得很成功,输入参数后摆杆能成功倒立,当有干扰的时候,可以迅速的恢复稳定,也学会了将不稳定的发散的系统通过状态反馈和观测,进行极点配置,从而得到稳定的系统。
最后通过直线一级倒立摆系统的设计使我掌握了用matlab来进行系
统设计仿真的方法,并熟悉了现代控制理论的相关知识,学以致用,学会一般系统的建模方法和控制器设计方法。
很好的将课本上的知识运用到了生活中
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 现代 控制 一级 倒立