2019基于MATLAB仿真的残障电动轮椅车速控制系统设计报告

1、泡南呼黑学院电子课程设计报告题目：基于MATLAB仿真的残隙电动轮椅车速控制系统设计课程：自动控制原理学生姓名：学生学号：年级：专业：班级：指导教师：机械与电气工程学院制20152015 年 3 3 月基于MATLAMATLA仿真的残障动轮椅车速控制系统设计1 1.课程设计的任务与要求1.11.1设计课题基于MATLAB仿真的残障电动轮椅车速控制系统设计系统结构图为:其中，控制器为Gc=K,传感器模型为G1(s)、被控对象模型为s2编程绘制系统的根轨迹，并求分离点坐标的增益K值、与虚轴交点时(4)编程当系统的速度误差系数K=8时，绘制系统单位阶跃响应和单位斜坡的曲线，利用MATLAB进行系统动

2、态特性分析(求出其性能指标o%,ts的值)和稳态误差。G2(s)=0(s1)(0.25s1)1.21.2课程设计的任务(D简述具有头盔传感器的残障电动轮椅车速控制的基本原理及应用;(2)编程求解系统的传递函数C0;R(s)(3)图1 1系统方框图(5)编程绘制当K=8时系统的Bode图、Nyquist图，求出相角裕量和幅值裕量，判断系统稳定性。(6)当控制器为P控制律Gc=4(2s+1)时，编程此时系统的Bode图、Nyquist图，求出相角裕量和幅值裕量，判断系统稳定性。1.31.3课程设计的目的(1)正确理解传递函数及根轨迹的概念；(2)掌握根轨迹的绘制法则，能熟练绘制跟轨迹；(3)根据根

3、轨迹定性分析系统指标随参数变化的趋势；(4)熟练使用MATLAB工具绘制系统的根轨迹和传递函数；(5)使用MATLAB工具对系统进行稳定性分析、稳态误差分析以及动态特性析；(6)使用MATLAB工具画出Bode图、求出相角裕量和幅量并判断系统的稳定性。2 2 . .具有头盔传感器的残障电动轮椅车速控制的基本原理及应用一种新型的电动轮椅装有一种非常实用的速度控制系统，使颈部以下有残障的人士也能自动驾驶这种电动轮椅。该系统在头盔上以间隔90安装了四个速度传感器，用来指示前、后、左、右四个方向。头盔传感系统的综合输出与头部运动的幅度成正比。国内对电动轮椅的研究较晚，尤其是智能电动轮椅，研究还不完善，

4、但近几年发展很快。国内厂商生产的电动轮椅大部分为四轮式和六轮式，一般都具有调速、翻越简单路障和防倾倒等功能。虽然国内电动轮椅研究还不太完善，但在一定的基础上还是有所提高的。例如有些生产商在原有轮椅的研究上，发明出利用驱动左、右动力后轮的左、右电机串联连接设计，从而具有差动速度功能，使电动轮椅行驶时稳定舒适、转向可靠。近几年，还出现了手扶电动、可爬梯以及站立式电动轮椅。随着机器人技术、人工智能技术和传感器技术的进步，电动轮椅的研究朝着高性能、多功能、 智能化和人性化的方向发展。 智能轮椅不但可以为老年人和残疾人提供一种良好的代步工具，而且可以具有自主导航、自主避障、人机对话等服务机器人所具有的各

5、种功能，因而可以帮助残疾人和老年人提高自己的生活自理能力和工作能力，使他们更好地融入社。初期的研究，赋予轮椅的功能一般都是低级控制，如简单的运动、速度控制及避障等。随着机器人控制技术的发展，移动机器人大量技术用于轮椅，电动轮椅在更现实的基础上，有更好的交互性、适应性、自主性。应用领域随着科学技术的发展，电动轮椅的强大功能不仅适用于年老体弱的老年人和重度残疾的伤患，同样的，它也适合于大型车间工人的代步工具。3 3 . .控制系统设计原理3.13.1反馈控制系统的基本知识反馈控制系统又称闭环控制系统，是在闭环控制系统中，把输出量检测出来，经过物理量的转换，再反馈到输入端去与给定值进行比较，并利用比

6、较后的偏差信号，以一定的控制规律产生控制作用，抑制内部或外部扰动对输出量的影响，逐步减少以致消除这一偏差，从而实现要求的控制性能。闭环控制的特点是：在控制器和被控对象之间，不仅存在着正向作用，而且存在反馈作用，既系统的输出量对控制量有直接影响，将检测出来的输出量送回到系统的输入端，并与信号比较的过程称为反馈，若反馈信号与输入信号想减，则称负反馈。反之，若相加，则称正反馈，输入信号与反馈信号之差称为偏差信号， 偏差信号作用于控制器上， 控制器对偏差信号进行某种运算，产生一个控制作用，是系统的输出量趋向于给定数值，闭环的实质就是利用负反馈的作用来减小系统的误差，因此闭环控制又称为反馈控制。图2 2

7、反馈控制系统方框图3.23.2控制系统时域分析法对于线性系统，常用的分析方法有三种：时域分析法，根轨迹法和频域分析法。时域分析法，具有直观准确的优点，并且可以提供系统时间响应的全部信息，尤其适用于低阶阶段。时域分析法是根据微分方程，利用拉氏变换直接求出系统的时间响应，然后按照响应的曲线来分析系统的性能。3.33.3控制系统根轨迹原理(1)所谓根轨迹是指，开环系统的每一个参数从零变化到无穷大时，闭环系统特征方程根在S平面上的轨迹称为根轨迹；(2)根轨迹的分离点与分离角是两条或两条以上根轨迹分支在S平面上相遇又立即分开的点，该点与实轴正方向的夹角即为分离角；(3)在理论分析中，往往只能画出根轨迹草

8、图，而利用MATLAB,则可以迅速绘制出精确的根轨迹图形。MATLAB绘制根轨迹的函数为rlocus,常用格式为rlocus(sys),sys为系统开环传递函数模型名称：rlocus(num,den,k),num为开环传递函数分子多项式，den为分母多项式，k为根轨迹增益，k的范围可以指定，若k未给出，则默认为k从0T+oo,绘制出完整的根轨迹。利用函数rlocfind可以显示根轨迹上任意一点的相关数值， 以此判断对应根轨迹增益下闭环系统的稳定性。3.43.4频率特性及图解法传统的频率分析是绘制频率特性曲线的渐近线，或通过人工计算数据，绘制较为详细的伯德图，奈氏图、对数幅相频率特性图，方法复杂

9、还不一定能保证绘制的精度。而应用MATLAB提供的相关函数，可以快速、精确地绘制出这三种图形的准确曲线，并计算出频域性能指标，对系统进行分析与设计。3.53.5稳态误差分析(1)频域法以控制系统特性作为数学模型，不必求解系统的微分方程或动态方程，而是做出系统频率特性的图形，然后通过频域和时域之间的关系来分析系统的性能，因而比较方便。频域的相对稳定性即稳定裕度常用相角裕度,和幅值裕度h来度量。相角裕度：设8c为系统的截止频率，则AM)=Gjc)Hjc)=1定义相角裕度为=18. G(j-c)H(2相角裕度的含义是，对于闭环稳定系统，如果系统开环相频特性再滞后,度，则系统处于临界稳定状态。幅值裕度

10、：设”为系统的穿越频率，则系统在、处的相角(x)G(j，x)H(jx)=(2k1)二；k=0,_T幅值裕度h的含义是，对于闭环稳定系统，如果系统开环幅频特性在增大倍,则系统处于临界稳定状态。(2)所谓时域分析法根据描述系统的微分方程的性能或传递函数，直接解出控制系统的时间响应，然后依据响应的表达式或描述曲线来分析系统的性能时域分析法包括稳定性分析、稳定性能分析(稳态误差)、动态性能分析三方面。控制系统的稳定性是由系统的闭环极点唯一确定的，而控制系统的动态性能则由该系统的闭环零、极定义相角裕度为1G(jx)H(jx)点所决定。因此，可以根据闭环的零、极点间接的研究控制系统的性能斜坡输入、阶跃输入

11、、加速度输入作用下的稳态误差与稳态误差系数：a、斜坡输入作用下的稳态误差与稳态误差系数2若r(t)=Rt,其中R表示速度输入函数的斜率，则R(s)=R/s当用静态速2度误差系数表示系统在斜坡输入作用下的稳态误差可将R(s)=R/S代入sR(s)ess(二)=limsE(s)=limssQ1G(s)H(s)KKV=limsG(s)H(s)=lims-71式中sT7称为静态速度误差系数，其单位为s0在斜坡输入作用下，速度误差的含义并不是指系统稳态输出与输入之间存在速度上的误差，而是指系统系统稳态输出与输入之间存在位置上的误差。b、阶跃输入作用下的稳态误差与稳态误差系数若r(t)=R*1，R为输入阶

12、跃函数的幅值，则R(s)=R/s。当R(s)=R/s时/、RRe5s(:)二上1limsG(s)H(s)1KP-Kp=limsG(s)H(s)华上0*八3口由 T，式中 T 称为静态位置块差系数。C、加速度输入作用下的稳态误差与稳态加速度系数23若r(t)=Rt/2,其中为加速度输入函数的速度变化率，则R(s)=R/s0如果用静态加速度误差系数表示系统在加速度输入作用下的稳态误差可得R2_lims2G(s)H(s)2_Ka=limsG(s)H(s)=lims式中s：0sQ3.6MATLAB3.6MATLAB软件的介绍MATLAB系统五个主要部分组成，下面分别加以介绍ess(二)=limsG(s

13、)H(s)s0RKVess(二)s0MATLAB是高层次的矩阵/数组语言， 具有条件控制、函数调用、数据结构、输入输出、面向对象等程序语言特点。利用它既可以进行小规模编程，完成算法设计和算法实验的基本任务，也可以进行大规模编程，开发复杂的应用程序。(2)MATLAB工作环境这是对MATLAB提供给用户使用的管理功能的总称，包括管理工作空间中的变量输入输出的方式和方法，以及开发、调试、管理M文件的各种工具。(3)图形图像系统这是MATLAB图形系统的基础，包括完成2D和3D数据图示、图像处理、动画生成、图形显示等功能的高层MATLAB命令，也包括用户对图形图像等对象进行特性控制的底层MATLAB

14、命令，以及开发GUI应用程序的各种工具。(4)MATLAB数学函数库这是对MATLAB使用的各种数学算法的总称，包括各种初等函数的算法，也包括矩阵运算、矩阵分析等高层次数学算法。4 4 . .课题设计的分析与计算4.14.1闭环传递函数的求解中(s)二3闭环传递函数公式为：1+H(s)G(s),由此公式计算可得3210.25s1.75s3.5s22.5k4.24.2闭环传递函数的仿真(1)MATLAB语言体系称为静态加速度误差系数。文件编嘱(i)宜看 M 播入工具受面爸口 MD 犁励dH曰国口864202468864202468SO.O.SSO.O.S(7(73PUO3PUO口由fl)m-xv

15、ajeu-gew-m-xvajeu-gew-3-2.5-2-1.5-0.5D3-2.5-2-1.5-0.5DRealAxis(seconds1)- -%TTAJJIkM1TTAJJIkM12.52.5J J二dexCO,2S,I.7S,dexCO,2S,I.7S,3.3.S,S,盒】:G=fGtun*dwrOG=f p=p=f f一亡dbdbulitG#11ulitG#11A=2.52.5_4-=_4-=_，一，=_4_4Q.N5sQ.N5s 3 3+1.+1.VSmFVSmF+ +3.5s3.5s+ +4.54.5C Conon工rruouwrruouw_ _*im*im - -七zrans

16、f-zransf mug25mug25: : d.en=d.en=O-25j,O-25j,1 1- -7 755S S3 322: : sys=sys=f(niinijf(niinijden)den); ; f fiaureiaure(1)(1) pzn.ap(sys)pzn.ap(sys) -fxgure-fxgure(2)(2) xJLocus(sys) K,PK,P=r=rlocfindde-n)Selectapointintheeraphicswindowlocfindde-n)Selectapointintheeraphicswindowse1eed_polot:=0,03790,0

17、379+ +0.049510.04951K K= =0-&554|0-&554|F F二4.00264.0026十0000100001- -1 1. .0&370&37+ +1 1. .4966n.4966n.T.T.09370937- -1.496611.49661当K=8时，由公式c(QC),结合下列仿真图所得数据，可求出超* *FigureFigure2 2,文件?看M M蕤入工具CDCD身面窗口他)帮助出)n na a$一二|入一发片|曰困|一图6MATLAB6MATLAB根轨迹仿真图分析：通过图像可以看出理论求出的分离点与图形根轨迹上的分离点的位置很相

18、近。从图中分离点坐标的增益K为0.85,而理论分析值K为0.67误差不大。6.6.分析系统动态特性分析和稳态误差6.16.1当K=8K=8时求和仃和上对于输入信号为单位阶跃信号时r(t)=1,则R(s)=1/s,对于I型具稳态误差:essG)=0对于输入信号为单位斜坡信号时r(t)=t,则R(s)=12，其稳态误差：S21会(二户厂0.125如)-c(二)二t2c/c/t=3.07sRealAxis(secondsRealAxis(seconds1 1) )864202468864202468(spu口口邕要 num=8;SK=3 den=O,25j1.75,3.5,2; gh=tf(nuji

19、ijden):石开环传函 sys=f=dback(gh,1):%1 表示单位负反馈，闭环传函t=0:0,1:10:u=t 淳输入为斜坡函数 y=lsim(sysJu,t)制系统占 ys 对由 u 和 t 描述的输入信号白 subplot2j1);%分割区域 plot3%?jtjUjb-): holdon； title(theResponsetoarampinput; xlabel(Time(sec); ylabel(Amplitude);subplot&2); step(sys)图7 7输入单位斜坡、阶跃信号程序图9 9K=8K=8时单位阶跃曲线分析:（1）速度误差随着时间的增大而减小

20、，所以它的速度误差系数为零，和计算结由仿真图我们可以得出结论：（2）位置误差随着时间的增大而逐渐减小，最后只存在微小误差，由计算可知，这个误差仅有0.18。分析：由单位阶跃响应曲线，可得出以下数据：ts=3.07sc（tp）=1.25c=0.785再由公式可得超调量为：c(tp)一。(二)1.25-0.785工=100%=100%=37.2%调节时间由响应曲线直接可得:ts=3.07s此项值与理论值一致。7 7.判断系统的稳定性7.17.1当K=8时，求出相角裕量和幅值裕量1）相角裕量的分析计算：在系统幅值穿越频率6c处，使系统达到临界稳定状态所需的附加相位滞后角称为相角裕量，用丫表示。即干1

21、80（+（）ac）若系统相角裕量为正值，闭环系统稳定；反之,如果相角裕量为负值，闭环系统不稳定。由Bode图可知，可以看到所有环节都作用后，系统依然没有穿越1800。年3.69oh=1.12由于60o,所以该系统稳定。2）幅值裕量的分析计算由于系统无右半平面的开环极点，从图11可以看出，奈奎斯特曲线不包围（-1,j0）点，系统稳定。7.27.2当K=8K=8的BodeBode和NyquistNyquist图c（二）1.251 k=8;ima=(2.5*k|：:den=0(25L：5152);figurebod*(0;gndInyquist(G)图10K=810K=8时的BodeBode和Nyq

22、uistNyquist程序G*=G*=J2S3P*P*= =3_3_999999产Vf=3.74213.7421Wc=3.53.5；1616A.菌Ewmtuc笠*俯r3s目-Ed-S)rtnGS)rtnGFreuenciiralsFreuenciirals；C7BC7BUagngUagng画173173S S陛nGnGFrwiMncy(ridFrwiMncy(rid1 108950895PMMidwJ:-6PMMidwJ:-61S1ST TFrequencyIraFrequencyIraiwiw图12K=812K=8时BodeBode图帆营网即4 4a a22i1010劭G GM3Jlnwg6

23、B2M3Jlnwg6B2：FHtqMhCFHtqMhC (ri9s(ri9s力倍ReaRea：AxeAxe分析：(1)相角裕量的分析计算：在系统幅值穿越频率0c处， 使系统达到临界稳定状态所需的附加相位滞后角称为相角裕量，用丫表示。即.180(,c)若系统相角裕量为正值，闭环系统稳定；反之,如果相角裕量为负值，闭环系统不稳定。对于高阶系统，一般难以准确计算截止频率，由图可知：Afec)=G(jo)+H(jo)=1.c=0.6244=3.66计算结果与MATLAB仿真结果相同，同时，该系统稳定。(2)幅值裕量的分析计算：由于系统无右半平面的开环极点， 从图11可以看出， 奈奎斯特曲线不包围(-1

24、,j0)点，系统稳定。8 8.当控制器为P P控制律G。二”如+1)时，系统的BodeNyquistBodeNyquist图。)皿也平。1/潟温力；仁怕b戒心(ndfiurt(2):ny(uisvG;。电/刊立墀诩G)InfInfPM=42.446942.4469g=InfInfu=8.35478.3547图1212当Gc=4(2s+1)时的BodeBode和NyquistNyquist的程序图1313当Gc=4(2s+1)+1)时的BodeBode图BodDiagramFfequencyraMFfequencyraM) )K K) )图1414当Gc= =4(4(2s2s+ +1)1)时的N

25、yquistNyquist图分析:这是高阶开环传递函数，可以直接从图中得出 一=8.35h=Inf=42.44因为俨0o,所以系统处于稳定的状态9 9.总结及体会9.19.1课程设计小结MATLAB不失为一款非常好的软件，但是我们所涉及的东西还是太少，就这次设计来说，可以使用MATLAB自带的控制系统工具箱来简单进行课设的各个环节。就我个人认为，在本次课设中遇到的困难不少，主要还是对课程即超前校正均里审DogrwnDogrwnRealRealSrtniGSrtniG% %1414I I卿4 4而SysmGSysmG刎二牺鹏3防Few7的不熟悉，在熟悉这个过程中花了不少的时间，实际在对MATLA

26、B进行制作的过程中，没有遇到很大的困难，无非是一些代码以及格式的学习。我认为，在这学期的实验中，不仅培养了独立思考、动手操作的能力，在各种其它能力上也都有了提高。更重要的是，在实验课上，我们学会了很多学习的方法。而这是以后最实用的，真的是受益匪浅。要面对社会的挑战，只有不断的学习、实践，冉学习、冉实践。这对于我们的将来也有很大的帮助。以后，不管有多苦，我想我们都能变苦为乐，找寻有趣的事情，发现其中珍贵的事情。就像中国提倡的艰苦奋斗一样，我们都可以在实验结束后变的更加成熟，会面对需要面对的事情。回顾起此课程设计，至今我仍感慨颇多，从理论到实践，在这段日子里，可以说得是苦多于甜，但是可以学到很多很

27、多的东西，同时不仅可以巩固了以前所学过的知识，而且学到了很多在书本上所没有学到过的知识。通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会服务，从而提高自己的实际动手能力和独立思考的能力。在设计的过程中遇到问题，可以说得是困难重重，但可喜的是最终都得到了解决。实验过程中，也对团队精神的进行了考察，让我们在合作起来更加默契，在成功后一起体会喜悦的心情。果然是团结就是力量，只有互相之间默契融洽的配合才能换来最终完美的结果。9.29.2收获体会经过不断的努力之后，我的课程设计终于完成了。课程设计不仅是

28、对前面所学知识的一种检验， 而且也是对自己能力的一种提高。 以前只是仅仅停留在书本的知识， 感觉自己什么都不懂，这次是真实的操作一次，感觉变得充实了。通过这次课程设计，我才明白学习是一个长期积累的过程，在以后的工作、生活中都应该不断的学习，努力提高自己知识和综合素质，特别是对于我，基础比较差，一定不能太过于心急，要静下心来慢慢的研究。在这次课程设计中也使我们的同学关系更进一步了，同学之间互相帮助，有什么不懂的大家在一起商量，听听不同的看法对我们更好的理解知识，所以在这里非常感谢帮助我的同学，我也明白学习不是埋头苦读书，而是合理的利用资源，从同学那里，老师那里得到的有用的想法和信息，特别是网上有

29、很多很好的资料，对自己的自学能力也是很好的提高。我的心得也就这么多了，总之，不管学会的还是学不会的的确觉得困难比较多，真是万事开头难，不知道如何入手。最后终于做完了有种如释重负的感觉。此外，还得出一个结论：知识必须通过应用才能实现其价值！有些东西以为学会了，但真正到用的时候才发现是两回事，所以我认为只有到真正会用的时候才是真的学会了。也不能因为做出来就以为自己什么都懂了，在很多的地方还需要进一步的去学习和研究，就想这个课程设计虽然能实现轮椅的车速控制的功能，但是识别率不是特别高的，只是在理论上能行得通，商业也已经运用到了。以前对MATLAB也是仅局限在数学实验上学的那些知识，但是通过此次的课程

30、设计，才了解到MATLAB在图形图像处理方面的强大功能。在此要感谢我们的指导老师对我们的悉心指导，感谢老师给我们的帮助，以前我是不相信我能做出来的，是您的谆谆教诲让我重新找到自信。在设计过程中，我通过?阅大量有关资料，与同学交流经验和自学，并向老师请教等方式，使自己学到了不少知识，也经历了不少艰辛，但收获同样巨大。在整个设计中我懂得了许多东西，也培养了我独立工作的能力，树立了对自己工作能力的信心，相信会对今后的学习工作生活有非常重要的影响。而且大大提高了动手的能力，使我充分体会到了在创造过程中探索的艰难和成功时的喜悦。虽然这个设计做的也不太好，但是在设计过程中所学到的东西是这次课程设计的最大收获和财富，使我终身受益。1010.参考文献1卢京潮.自动控制原理.西北工业大学M.2004.2刘卫国.MATLAB程序设计与应用.北京：高等教育