8-位置随动系统解读

1、指导教师评定成绩： 审定成绩： 重 庆 邮 电 大 学自 动 化 学 院自动控制原理课程设计报告 设计题目：位置随动系统单位（二级学院）： 学 生 姓 名： 专 业： 班 级： 学 号： 指 导 教 师： 设计时间： 重庆邮电大学自动化学院制目 录一、设计题目2二报告正文3摘要321 问题一的分析与求解422 问题二的分析与求解523 问题三的分析与求解1024 问题四的分析与求解14三、设计总结18四、参考文献19五、附录20附录一20附录二20一、 设计题目自动控制原理课程设计任务书1某位置随动系统原理如下图所示。输入量为转角，输出量为转角，为圆盘式滑动电位器，SM为伺服电动机，TG为测速

2、发电机。功放k3+15v-15v+15v-15v-k1-k2SMTG+30k20k10k10k10k10kuiuauoutu2u1RpRp要求：（） 查阅相关资料，分析系统的工作原理，指出被控对象、被控量和给定量，画出系统方框图。（） 分析系统每个环节的输入输出关系，代入相关参数求取系统传递函数。（） 分析系统时域性能和频域性能。（） 运用根轨迹法或频率法校正系统，使之满足给定性能指标要求。（已知条件和性能要求待定）二、 设计报告正文摘要随动系统是指系统的输出以一定的精度和速度跟踪输入的自动控制系统，并且输入量是随机的，不可预知的，主要解决有一定精度的位置跟随问题，如数控机床的刀具给进和工作台

3、的定位控制，工业机器人的工作动作，导弹制导、火炮瞄准等。控制技术的发展，使随动系统得到了广泛的应用。位置随动系统是反馈控制系统，是闭环控制，其位置指令是经常变化的，要求输出量准确跟随给定量的变化，输出响应的快速性、灵活性和准确性成了位置随动系统的主要特征。本次课程设计研究的是一类位置随动系统的滞后校正，首先通过分析原理求出传递函数，并利用主导极点进行降阶，得出一个二阶系统传递函数，并通过MATLAB分析时域和频域的各个性能，得出相角裕度太小和超调量太大，然后设计PD控制装置，改善系统的阻尼比，来使系统的各个性能达到要求。关键词：随动系统 主导极点 滞后校正 PD控制 阻尼比2.1 问题一的分析

4、与求解系统工作原理本系统利用圆盘式电位器作为转角输入量，在圆盘式电位器两端加有15V电源，每当圆盘式转动一定角度，圆盘式点位器上输出电压也跟随变化，该变化的电压与被控制圆盘式电位器反馈回来的电压值进行比较，得到一个电压差。圆盘式电位器与圆盘式电位器工作原理相同，电压差输入运算放大器，运算放大器上加有反馈电阻。根据运算放大器“虚短”与“虚亏”方法，运算放大器的放大系数为反馈电阻与放大器前端电阻10k的比值，即放大系数为3，电压差经过放大器放大后输出电压，后面测速发电机反馈回来电压值，电机转速不同，则输出电压也会不同，运算放大器输出电压与测速反馈电压进行比较，得到差值电压，运算放大器与工作原理及分

5、析方法相同，比例放大器放大系数为，运算放大器构成比例调节器，可快速响应电压变化。比例调节器输出电压，经过“功率放大器”放大功率，既可以放大输入电压又可以放大输出电压，以驱动直流伺服电机SM转动，直流伺服电机SM同时带动测速发电机TG。测速发电机是用来测量角速度并可将它转换为电压量，即反馈电压。直流伺服电机SM上连接减速器，减速器和齿轮转化系统与相连接，即可就能够调节位置。功放输出电压大，对应电动机SM转速高，相应电位器也会转动较大角度。相反。电压小，电动机转速低，则电位器转动角度较小。同时电位器也时刻反馈回来其所在位置，即反应为反馈电压，电位器每转动一个角度，对应的控制电位器转动一个角度，这样

6、便实现了位置随动的控制。其中，控制系统的被控对象是电动机；被控量是电压；系统的输入量为转角，输出量为转角；给定量是圆盘式滑动变位器，运算放大器，运算放大器，功率放大器，伺服电动机，测速发电机。系统方框图如下：2.2 问题二的分析与求解2.21 电位器：电位器是一种把线位移变换为电压量的装置。在控制系统中，单个电位器用作为信号变换装置。+15v-15v圆盘式滑动电位器方框图如下： 可得传递函数为：，得。其中为电位器最大工作角。2.22 放大环节 实物图如下：放大环节的微分方程为c (t)=Kr（t） ，式中，K 为常数，称放大系数或增益。放大环节的传递函数为 G（s ）=K 。 放大环节的方框图

7、如图所示。在一定的频率范围内，放大器、减速器、解调器和调制器都可以看成比例环节。 2.23 功率放大器：实物图如下：功率放大器：功放(功率放大器)的原理就是利用三极管的电流控制作用或场效应管的电压控制作用将电源的功率转换为按照输入信号变化的电流。 故传递函数为 2.24 电机 实物图如下：电机的传递函数求解如下：电枢回路电压平衡方程 (2-4-1)式中是电枢旋转时铲射的反电势,其大小与激磁磁通成正比,方向一样电枢电压相反,即=,是反电势系数。电磁转矩方程 （2-4-2）式中，是电机转矩系数；是电枢电流产生的电磁转矩。电动机轴上的转矩平衡方程 （2-4-3）式中，是电动机和负载折合到电动机轴上的

8、粘性摩擦系数；是电动机和负载折合到电动机轴上的转动惯量。由式（2-4-1）（2-4-3）中校区中间变量，及，便可得到以为输出量，为输入量的直流电动机微分方程： （2-4-4）在工程应用中，犹豫电枢电路电感较小，通常忽略不计，因而上式可简化为 （2-4-5）式中，=是电动机的时间常数；=，=是电动机传递系数。上面我们已经求的电枢控制直流电动机简化后的微分方程为：式中可视为负载扰动转矩。根据线性系统的叠加原理，可分别求到和到的传递函数，以便研究在和分别作用下的电动机转速的性能，将他们叠加后，便是电动机转速的相应特性。为求，令=0，则有在初始条件下，即=0时，对上式各项求拉氏变换，并令，则的s的传递

9、方程，由传递函数定义，于是有 下图是它的方框图 2.25 测速发电机 由于测速发电机接有负载电阻，故测速发电机的传递函数为：式中，是测速发电机经分压后的输出电压；为测速发电机输入角度。下图为测速发电机的方框图s 故传递函数为 根据各环节的输入输出关系及系统的结构框图可求得传递函数如下： 式中为电位器的传递系数；，为运算放大器的比例系数；为功率放大器的放大系数；为电动机的转矩系数；和分别为电动机电枢绕组的电感和电阻；J和f分别为折算到电动机轴上的总转动惯量和总粘性摩擦系数；为与电动机反电动势有关的比例系数。传递函数可变为： （电枢电感很小，可知为非主导极点可舍去。） 则 （令=,K=）2.3 问

10、题三的分析与求解利用MATLAB辅助分析法分析系统的频域性能和实域性能：令 又 带入，得把 3.464由MATLAB计算（程序见附录一，各参数计算公式见附录二）可得： = 3.4640 （无阻尼振动频率） =3.3703 （阻尼振荡频率） =76.6440 （阻尼角） =0.5352 （上升时间） =0.9321 （峰值时间） =0.4743 （超调量） =3.7491 （误差带=5%时的调节时间） =4.9989 （误差带=2%时的调节时间） =2.2247 （谐峰峰值） =3.2739 （谐峰频率） =5.1777 （带宽频率） =3.2843 （截止频率） =25.9787 （相角裕度）

11、相应的图如下：图1 零极点图图2 根轨迹图图3 波特图图4 奈奎斯特图 图5 系统的单位脉冲响应图图6 系统的单位阶跃响应图图7 系统的单位斜坡响应图2.4 问题四的分析与求解由所得的数据和各响应图可知系统由于相角裕度太小（25.978730）,至使超调量太大（0.4743=47.43%）。因此可用PD控制校正相角裕度。校正后要求，求所应添加的校正框图根据校正要求可知用PD控制由 可得新的闭环特征方程为：求等效阻尼比 所以校正后的特征方程为：其中 校正后所得的参数如下： = 3.4640 （无阻尼振动频率） = 2.8260 （阻尼振荡频率） = 54.6689 （阻尼角） = 0.7740

12、（上升时间） = 1.1117 （峰值时间） = 0.1079 （超调量） = 1.4976 （误差带=5%时的调节时间） = 1.9968 （误差带=2%时的调节时间） = 1.0598 （谐峰峰值） = 1.9933 （谐峰频率） = 4.0760 （带宽频率） = 2.5318 （截止频率） = 57.7098 （相角裕度）相应的图如下： 图8 零极点图图9 根轨迹图图10 波特图。图11 奈奎斯特图图12 系统的单位脉冲响应图图13 系统的单位阶跃响应图图14 系统的单位斜坡响应图校正后相角裕度（超调量为）由以上可知满足校正要求。三、 设计总结参考文献【1】 姜启源，谢金星，叶俊.数学

13、模型（第三版）M.北京：高等教育出版社,2003.8.【2】 孙亮，杨鹏.自动控制原理M.北京：北京工业大学出版社，1999.9.【3】 毕云飞，李世平.新概念office2003三合一图解教程M.北京：清华大学出版社,2004.9.【4】 孙亮.MATLAB语言与控制系统仿真M.北京:北京工业出版社,2002.7.【5】 翁剑枫，叶志前.MATLAB LabVIEW SystemView仿真分析基础M. 北京：机械工业出版社,2005.1.【6】 胡寿松.自动控制原理M.北京：科学出版社,2007.6.附录一：MATLAB源代码如下：zeta=0.231;wn=3.464wd=wn*sqrt

14、(1-zeta*zeta)m=atan(sqrt(1-zeta*zeta)/zeta)*180/pitr=(pi-m/180*pi)/wdtp=pi/wdmp=exp(-zeta*pi/(sqrt(1-zeta*zeta)ts1=3/(zeta*wn)ts2=4/(zeta*wn)mr=1/(2*zeta*sqrt(1-zeta*zeta)wr=wn*sqrt(1-2*zeta*zeta)wb=wn*sqrt(1-2*zeta*zeta+sqrt(2-4*zeta*zeta+4*zeta*zeta*zeta*zeta)wc=wn*sqrt(sqrt(1+4*zeta*zeta*zeta*zeta)-2*zeta*zeta)r=atan(2*zeta/sqrt(sqrt(1+4*zeta*zeta*zeta*zeta)-2*zeta*zeta)*180/pinum=wn*wn;den=1 2*zeta*wn wn*wn;sys=tf(num,den);figure(1)pzmap(sys);den1=1 2*zeta*wn;sys1=tf(num,den1);figure(2)rlocus(sys1);figure(3)margin(sys);figure(4)nyquist(sys);axis quualp=roots(den)t