版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、Good is good, but better carries it.精益求精,善益求善。matlab大作业-平衡杆小球位置控制系统-华中科技大学-MATLAB语言、控制系统分析与设计大作业平衡杆小球位置控制系统设计与仿真专业:电气工程及其自动化班级:设计者:学号:华中科技大学电气与电子工程学院2008年1月平衡杆小球位置控制系统设计与仿真一、问题描述与实验要求Aballisplacedonabeam,seefigurebelow,whereitisallowedtorollwith1degreeoffreedomalongthelengthofthebeam.Aleverarmisatta
2、chedtothebeamatoneendandaservogearattheother.Astheservogearturnsbyanangletheta,theleverchangestheangleofthebeambyalpha.Whentheangleischangedfromtheverticalposition,gravitycausestheballtorollalongthebeam.Acontrollerwillbedesignedforthissystemsothattheballspositioncanbemanipulated.Forthisproblem,wewil
3、lassumethattheballrollswithoutslippingandfrictionbetweenthebeamandballisnegligible.Theconstantsandvariablesforthisexamplearedefinedasfollows:Mmassoftheball0.11kgRradiusoftheball0.015mDleverarmoffset0.03mGgravitationalacceleration9.8m/s2Llengthofthebeam1.0mJballsmomentofinertia9.99e-6kgm2Rballpositio
4、ncoordinateAlphabeamanglecoordinateThetaservogearangleSystemEquationsTheLagrangianequationofmotionfortheballisgivenbythefollowing:Linearizationofthisequationaboutthebeamangle,alpha=0,givesusthefollowinglinearapproximationofthesystem:Theequationwhichrelatesthebeamangletotheangleofthegearcanbeapproxim
5、atedaslinearbytheequationbelow:Substitutingthisintothepreviousequation,weget:DesignrequirementsThedesigncriteriaforthisproblemare:Settlingtimelessthan3secondsOvershootlessthan5%1数学模型的建立将上面推导的简化式做拉普拉斯变换,得到:变化后得到:化简为:2设计目标希望能精确小球的位置,即要求小球的稳态误差为零,同时希望因扰动引起的稳态误差也能为零。对于动态性能的要求,希望小球能较快且平稳,期望调节时间Ts为3s,超调量小
6、于5%。根据时域和频域指标的关系,可将时域性能指标转换为频率响应的约束条件。如系统的带宽与闭环系统自然振荡频率n和阻尼比有关,而n与调节时间有关。相角裕度PM和阻尼比有关,进而与超调量相关。(1-4)zeta=-log(.05)/sqrt(pi2+(log(.05)2);PM=100*zeta;wbw=(4/(3*zeta)*sqrt(1-2*zeta2)+sqrt(4*zeta4-4*zeta2+2);得0.6901,PM69.0107deg,wbw1.9785rad/sec3开环响应首先用MATLAB描述上述模型,并观察开环系统阶跃响应。J=9.99e-6;m=0.11;R=0.015;g
7、=9.8;d=0.03;l=1.0;num=m*g*d*R*R;den=(J+m*R*R)*l00;ball=tf(num,den)bode(ball)margin(ball)Bode图如下:图1-1小球位置开环频率特性曲线由图1-1知原系统的相角裕度为0deg,截止频率为0.458rad/sec,由此判断闭环系统的调节时间不能满足设计要求。期望的截止频率为1.9785rad/sec,相角裕度为70deg,而校正前系统在1.9785rad/sec处的幅度为-20dB,相角为-180deg。需增加系统截止频率和相角裕度。二、闭环系统控制器设计概述将小球的位置做完控制量,构成闭环,系统闭环控制框图
8、如图2-1所示,控制器设计采用串联校正,设计思想采用开环频率特性设计方法,即根据系统的开环频率性能指标确定校正器参数。图2-1闭环系统控制框图设计步骤:(1)根据稳态误差要求确定系统控制器结构;根据调节时间的要求确定截止频率。(2)确定未校正系统需增加的最大相位超前角及新的幅值穿越频率;(3)确定控制器的传递函数;(4)校验,验算校正后系统性能。控制器设计(1)确定控制器的结构由上图可知,为了达到象角要求,需要增加一个积分环节并添加一个零点,构成PI控制器,零点的位置以低于截止频率10倍频程即可。由期望的截止频率wc所对应的幅值Mag(wc),确定PI的增益Kp,Ki。()()由此确定PI控制
9、器的传递函数为:()J=9.99e-6;m=0.11;R=0.015;g=9.8;d=0.03;l=1.0;num=m*g*d*R*R;den=(J+m*R*R)*l00;ball=tf(num,den);w=logspace(0,2,101);mag,phase,w=bode(num,den,w)L=find(w5.5);wc=w(L(1)Mag=20*log10(mag(L(1)Kp=10(-Mag/20)Ki=0.1*wc*Kpsysc=tf(KpKi,10)bode(sysc*ball);margin(sysc*ball);图2-2PI校正后系统BODE图根据图2-2所示,加PI校正后
10、系统的截止频率为2rad/sec满足要求,但相角裕度PM为-5.7deg,需增加滞后校正环节,欲减小的相角由超调量的要求确定。(2)确定超前校正环节参数根据超前校正步骤,已知待补偿相角和计算校正器参数和。()()由此确定超前校正控制器的传递函数为:从上面的图可以看出,相角裕度应该再提高96。但是可知,超前校正控制器的加入会影响截止频率,从而减小相角裕度,所以在这里,我们选用大于96的超前角100。PM=100;a=(1+sin(PM*pi/180)/(1-sin(PM*pi/180);mag0,phase0,w0=bode(sysc*ball);adb=20*log10(mag0);am=-1
11、0*log10(a);wgc=spline(adb,w0,am);T=1/(wgc*sqrt(a);con2=tf(a*T1,T1);w=logspace(2,3,101);bode(con2*sysc*ball,w)margin(con2*sysc*ball)图形如下:图2-3加超前校正的系统BODE图可以看到相角裕度为78.3,满足相角裕度要求,截止频率比期望的大了一些。3)确定控制器的零极点和增益通过调整增益K1使系统的截止频率在指定范围。mag1,phase1,w1=bode(con2*sysc*ball,w)L1=find(w15);wc1=w1(L1(1)Mag1=20*log10
12、(mag1(L1(1)K1=10(-Mag1/20)bode(K1*con2*sysc*ball,w)margin(con2*sysc*ball);图2-4校正后的系统BODE图(4)校验,验算校正后系统性能。这时候的时域响应情况为:figure(2)sys_cl=feedback(K1*sysc*con2*ball,1);%CLOSED-LOOPSYSTEMt=0:0.01:1;step(sys_cl,t)图2-5校正后闭环系统阶跃响应由上图可知,响应时间非常的短,远远低于期望的3sec,但是系统的超调量为33%,也远大于需要的5%。在这里,反复的调整参数,PM于截止频率Wc。最后当PM=8
13、5,Wc=2rad/sec时,得到如下响应曲线:图2-6校正后闭环系统阶跃响应响应的超调量为4%,响应时间为0.66sec,符合设计要求。校验扰动下的稳态误差:d_cl=feedback(ball,K1*sysc*con2);%DISTRUBANCESTEPsys_d=zpk(d_cl)figure(3)step(sys_d,t)图形如下:图2-7校正后扰动输入系统阶跃响应三、结论由以上的仿真和设计可以看到,设计出的系统完全符合要求:控制器为:controller=K1*sysc*con2系统传递函数如下:0.20946s(s+229)-(s+0.44)(s+0.2086)(s2+228.4s
14、+2.484e004)性能指标为超调量4%,调节时间0.66s,单位阶跃输入和扰动输入的稳态误差均为0。四、课程学习体会控制理论作为一门以应用为主的学科,其灵魂在于工业应用。因此,在长期逗留在书本上的学习中,我们虽然表面上掌握了许多控制理论知识,但是感性认识的缺乏以及控制中许多约定俗成的规范,使得我们脑海中始终没有形成一幅控制理念的图画,感觉课程深奥,难懂。非常幸运的,我们能够开设一门matlab课程设计的课程,它让我们在一台简单的计算机上,就能够看到控制的过程;简单的换掉几个参数,就能很快仿真出新的控制结果,给我们理解控制理论以及其他种种缺乏感性认识的课程提供了强有力的工具。让人欣喜!感谢老
15、师在教学过程中的指导!五、附录本文所设计程序清单:J=9.99e-6;m=0.11;R=0.015;g=9.8;d=0.03;l=1.0;num=m*g*d*R*R;den=(J+m*R*R)*l00;ball=tf(num,den);bode(num,den)margin(num,den)figure(2)w=logspace(0,2,101);mag,phase,w=bode(num,den,w)L=find(w2);wc=w(L(1)Mag=20*log10(mag(L(1)Kp=10(-Mag/20)Ki=0.1*wc*Kpsysc=tf(KpKi,10)bode(sysc*ball,
16、w)margin(sysc*ball)figure(3)PM=85;a=(1+sin(PM*pi/180)/(1-sin(PM*pi/180);mag0,phase0,w0=bode(sysc*ball);adb=20*log10(mag0);am=-10*log10(a);wgc=spline(adb,w0,am);T=1/(wgc*sqrt(a);con2=tf(a*T1,T1);w=logspace(2,3,101);bode(con2*sysc*ball,w)margin(con2*sysc*ball)figure(4)mag1,phase1,w1=bode(con2*sysc*ball,w)L1=find(w12);wc1=w1(L1(1)Mag1=20*log10(mag1(L1
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 电网劳务分包投标方案(技术方案)
- 2024年重庆市消防操作员《中级技能-维保方向》科目真题冲刺卷9月份B卷
- 陶瓷机械产业深度调研分析及投资机会研究报告模板
- 球体类产品(新)
- 小学教学小论文13篇
- 2025版高考英语总复习高考语法与写作专题8定语从句教学案外研版
- 心系健康网创业项目策划书
- 信息必刷卷01(湖南长沙专用)(参考答案)
- 考点10 记叙文插叙的作用-2024年中考语文现代文阅读高频考点解题指津(全国)(解析版)
- 全球传播导论 课件 高金萍 第1-5章 全球传播的概念与要素-全球传播与国家形象
- 不搞恶作剧(小学安全教育课件)
- 囊性肾癌(课件)
- 厂房钢结构屋面及彩钢板施工方案
- (整理)五十音图字帖
- 山东电力物资部各岗位职责(共5篇)
- 藤野先生教学设计(朱则光)
- 基站电源规划配置原则-V3
- 第一课-笔墨千秋
- 实验小学新苏教版四年级科学上册10.《摩擦力》第1课时教案
- 高压电缆直埋技术要求
- 最新人教版三年级上册数学《4.2.2 三位数减两位数(连续退位)》PPT课件
评论
0/150
提交评论