原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
继教学院 X 禹志
控制系统实践大作业
球杆控制系统设计与仿真
学号：X姓名：禹志
球杆控制系统设计与仿真
一 球杆系统的数学模型
球杆系统的物理模型如下图所示：
图1：球杆系统模型
在这个模型里，我们假设球在杆上只有纯滚动而没有滑动并且我们忽略球和杆之间的摩擦力。我们定义这个模型中的变量和常量如下：
m
球的质量
0.1 kg
R
球的半径
0.02
d
轮的半径
0.05
g
重力加速度
9.8 m/s^2
L
杆长
1.0 m
J
球的转动惯量
1.0e-5 kgm^2
r
球的位置
杆与水平的夹角
伺服机构与水平的夹角
表1：模型的变量和常量
球杆系统的数学建模：
系统的动能：
系统的势能：
该系统的拉格朗日方程为：
通过拉格朗日法得到球的运动方程：
若我们假设杆与水平的夹角是小角度且，对上面的方程进行线性化，我们得到：
杆的角度与伺服机构的角度的关系可以近似的线性化为：
综合上面几个方程我们得到：
对方程两边做拉普拉斯变换：
整理后我们得到系统线性化后的开环传递函数：
二 控制器的设计
在这里，我们采用普遍使用的PID控制器。整个系统的闭环框图可以表达为：
图2：线性化后的控制框图
令 ，得：
我们取=21，=0，=14，将球稳定在杆的中点，仿真结果如下图所示：
图3：仿真结果
通过仿真结果表明将系统线性化后，该系统的超调量小于1%，上升时间3秒。
如果我们进一步考虑该模型的非线性性，可以得到
在Matlab Simulink中我们可以用如下方框图表示和的关系：
图4：到的方框图
其中f(u)是 (m*g*sin(u[3])- m*(L-u[1])*(u[4])^2) /(J/(R^2)+m)，u[1]是r， u[2]是，u[3]是， u[4]是。
这样我们可以得到非线性情况下的系统框图：
图5：非线性系统框图
我们取=1.8，=0，=3.8，将球稳定在杆的中点，仿真结果如下图所示：
图6：仿真结果
通过仿真结果表明将系统按照非线性处理后，该系统的超调量小于1%，上升时间约8秒。
,,,的响应曲线如下图所示：
图7：,,,的响应曲线
基于GUI的动画仿真效果如下图所示：
图8：基于GUI的仿真效果
附件：
1 系统线性化后的仿真程序 ballbeam.mdl
2 系统按非线性的仿真程序 nonlinearballbeam.mdl
3 基于GUI的动画显示 ballfig.m ballfig.fig，运行时请运行ballfig.m
注：仿真程序及GUI动画均有Matlab 6.5 版本下完成并通过。
参考文献：
Derivation of the Dynamics of the Ball and Beam System ，Dr Ben Cazzolato
Control of a Ball and Beam System ，Wei Wang ，School of Mechanical Engineering The University of Adelaide ，South Australia 5005 ，AUSTRALIA
文档评论（0）