一类位置随动系统的测速反馈控制

1、一类位置随动系统的测速反馈控制1位置随动系统原理1.1位置随动系统工作原理图1-1位置随动系统原理图该系统为一自整角机位置随动系统，用一对自整角机作为位置检测元件，并形成比较 电路。发送自整角机的转子与给定轴相连：接收自整角机的转子与负载轴（从动轴）相连。 TX与TR组成角差测量线路。若发送自整角机的转子离开平衡位置转过一个角度0，则在 接收自整角机的单相绕组转子的单相绕组上将感应出一个偏差电压、，它是一个振幅为 uem、频率与发送自整角机激励相同的交流调制电压。即ue = um sinOt在一定范围内，uem 正比于0 -0，即u = k 0 -0 ，所以可得u = k 0 -0 sin o

2、t这就是随动系统中接收自r cem er ce er c整角机所产生的偏差电压的表达式，它是一个振幅随偏差（0r-0c）的改变而改变的交流电 压。因此，ue经过交流放大器放大，放大后的交流信号作用在两相伺服电动机两端。电动 机带动负载和接收自整角机的转子旋转，实现0,=0c，以达到跟随的目的。为了使电动机 转速恒定、平稳，引入了测速负反馈。系统的被控对象是负载轴，被控量使负载轴转角0，电动机是执行机构，功率放大器 器信号放大作用，调制器负责将交流调制为直流电供给直流测速发电机工作电压，测速电 动机是检测反馈元件。1.2单元电路模块分析1.2.1自整角机自整角机是常用的位置检测装置，将角位移或者

3、直线位移转换成模拟电压信号的幅值 或相位。自整角机作为角位移传感器，在位置随动系统中是成对使用的。与指令轴相连的 是发送机，与系统输出轴相连的是接收机。则自整角机的表达式为u(t) = K 0 (t)-0 (t) = KAO(t)在零初始条件下，拉氏变换为u(s) = KgA0 (s)，则自整角机的传递函数为G1( s)=u (s)AO (s)自整角机的结构图如图1-2所示u (s)图1-2自整角机1.2.2功率放大器由于运算放大器具有输入阻抗很大，输出阻抗小的特点，在工程上被广泛用来作信号 放大器。其输出电压与输入电压成正比，即有u (t) = K u(t) - u, (t)G2(s)=在零

4、初始条件下，拉氏变换为u (s) = K u(s)-(s) = KAu(s)，则传递函数为 七(s) = K u (s)a式中ua(s)为输出电压，u(s)为输入电压，Ka为放大倍数。图1-3功率放大器1.2.3两相伺服电机由伺服电机的性质可知，则表达式为T tM + = Ku (t)m dt 2 dt m a在零初始条件下，拉氏变换为(T s2 + s)e (s) = K u (s)，则伺服电机的传递函数为G (s) =mM = K3 u (s)s (T s +1)伺服电机的结构图如图1-4所示图1-4两相伺服电动机1.2.4直流测速发电机与调制器测速发电机的输出电压ut与其转速O成正比，即

5、有 于是可得测速发电机的微分方程在零初始条件下，拉普拉斯变换u (s) = K舶(s)，则直流测速发电机的传递函数为 tt cG (S)=尝=Ks4 e (s) t直流测速电机的结构图如图1-5所示。Kse Jis)图1-5直流测速发电机与调制器1.2.5减速器由系统的减速装置可知，减速器的表达式为e (t) = i-o (t)在零初始条件下，拉普拉斯变换为e(s) = i-e(s)。则减速器的传递函数为 cmG (s) = c(s) = i5 em(s)其中i为减速比。减速器的结构图如图1-6所示9图1-6减速器91.3各部分元件传递函数2)交流放大器G (s) = & 3)= K2 u (

6、s)a3)两相伺服电机G (s)=七里)3 u (s)s(T s +1)其中T = R J (Rf + C C )是电动机机电时间常数;m a m a m m eK广CJ (RJm + Cfe)是电动机传递系数;4)直流测速发动机G (S)=竺立=K s49 (s)t5)减速器 G (s) = (，= i59m(s)1.4位置随动系统的结构图从与系统输入量9有关的比较点开始，依据电路单元模块结构图中的信号流向，把各 r环节的结构图连接起来，将系统的输入量9于结构图的最左端，系统输出量9于结构图的最右端，便得到系统结构图，如图1-7所示图1-7位置随动系统的结构图则系统的信号流图如图1.8所示1

7、.5相关函数的计算1)开环传递函数：G(s)KKJJT s 2 + (KK K + 1)s2)闭环传递函数：中(s)=KKKmT s2 + (KK K + 1)s + K K K imt a m& a m电动机机电常数：T = RJ /(R f + C C )m a L a L m e电动机传递系数：K =C /(R f + C C )m m a m m e初始条件：放大增益为K =50，电桥增益气=3，测速电机增益K = 0.18V-s ,R = 7.5。，L = 14.25mH，J = 0.007kg - m2，Ce = Cm = 0.3N - m / A，f = 0.2N - m - s

8、，减速比 i = 0.1 。在MATLAB中调用tf()函数和feedback()函数，求系统的开、闭环传递函数。代码如下：Ka=50; Ke=3;Kt=0.18;Ra=7.5; La=0.01425;J=0.007*9.8; Cm=0.3;Ce=0.3; f=0.2; i=0.1;Tm=Ra*J/(Ra*f+Cm*Ce);Km=Cm/(Ra*f+Cm*Ce);%开环传递函数分子系数，按s降幕排列%开环传递函数分母系数，按s降幕排列%调用tf()函数，求出开环传递函数%调用feedback()函数，求出单位反馈闭环传递函数num=Ka*Km*Ke*i;den=Tm,Ka*Km*Kt+1,0;s

9、1=tf(num,den)sys=feedback(s1,1)程序运行结果:开环传递函数2.830.3236 s 2 + 2.698 s闭环传递函数2.832.830.3236s 2 + 2.698s0.3236 s 2 + 2.698 s + 2.83即有开环传递函数为G(s) =气 T s 2 + (KK K + 1)s即有闭环传递函数为0(s)=KKKJT s2 + (KK K + 1)s + K K K imt a m& a m2.830.3236s 2 + 2.698s + 2.831.6开环系统频域特性求解求系统的幅值裕度和相角裕度，可直接调用margin()函数。margin()

10、函数可以从频率 响应数据中计算出幅值裕度、相位裕度及其对应的角频率。调用格式为margin(sys)其中 sys为系统的开环传递函数。代码如下：s1=tf(num,den);margin(sl);%调用margin()函数，求校正前系统的相角裕度和幅值裕度grid on;Matlab运行结果如下图所示。校正前，截止频率七=1.04rad/ s ；相角裕度y = 829 ； 幅值裕度为+sdB。Bode DiagramGm = Inf dB (at Inf rad/sec) , Pm = 82.9 deg (at 1.04 rad/sec)50I I I I_IIrnI I I I IImi i

11、 i i_iimio-1Frequency (rad/sec)10图1-9位置随动系统的频域特性1.7绘制Bode图和Nyquist图在Matlab环境下绘制系统的Bode图和Nyquist图。则输入代码如下:G=tf(2.83,0.3236 2.698 0);figure(1)bode(G);figure(2)nyquist(G);运行程序得系统的Bode图和Nyquist图。分别如图1-10和图1-11所示Bode Diagram5Q iiiiiijI I iiiiims省mu烈巨(Eiilip(DMnsJZd10110101102103Frequency (rad/sec)图1-10位置

12、随动系统的Bode图图1-11位置随动系统的Nyquist图2设计测速反馈控制2.1控制要求要求在原来的基础上设计一测速反馈控制使其阻尼比变为0.75,即是通过改变反馈比 K来改变整个系统的开环传递函数的阻尼比。2.2测速反馈控制的原理输出量的导数同样可以改善系统的性能。通过将输出的速度信号反馈到系统输入端与 误差信号进行比较，其效果与比例一微分控制相似，可以增大系统阻尼，改善系统动态性 能。如果系统的输出量是机械位置，如此系统就是角位移，则可采用测速发电机将角位移 变换为正比于角速度的电压，从而获得速度反馈。2.3设计测速反馈控制改变阻尼比有上面计算可得，位置随动系统的闭环传递函数为中(s)

13、=气人上.T s2 + (KK K + 1)s + K K K imt a m8 a m通过设计测速反馈控制使系统的阻尼比为0.75，则只需要改变K的值，而不需要改 变系统的其它系数。有二阶系统的闭环传递函数的标准形式：0( s)=2s2 + 2。 s + 2则其中 2 = KKK：=兰=8.7454, 2S = KKaKm +1,75。即求得n T 0.3236n TO = 2.9573, K = 0.0462加入测速反馈控制后的系统结构图和信号流图与未校正之前的只有K的取值不同，加 入测速反馈控制后的传递函数中用K = 0.0462.即此时的开环传递函数G，(s)和闭环传递函 数0( s)

14、分别为g ,(s)二KKKi=283T s2 + (KK K + 1)s 0.3236s2 + 1.436s0( s)=KKKi=2T s 2 + (KK K + 1)s + KKK i 0.3236s 2 + 1.436s + 2.83mt a m8 a m2.4校正前后系统的单位阶跃反馈曲线在Matlab软件环境中输入Matlab程序代码，分别得到校正前后的单位阶跃反馈曲线。输入代码如下：num=2.83;den=0.3236 2.698 2.83;G=tf(num,den);step(G)axis(0 8 0 1.5)grid minor即得系统校正前的单位阶跃反馈曲线如图2-1所示St

15、ep Response1234567Time (sec)图2-1系统校正前的单位阶跃反馈曲线输入代码如下：num=2.83;den=0.3236 1.436 2.83;G=tf(num,den);step(G);grid minoraxis(0 5 0 1.5)即得系统校正后的单位阶跃反馈曲线如图2-2所示Step Response图2-2系统校正后的单位阶跃反馈曲线3校正前后系统的MATLAB仿真分析3.1系统的Simulink 的仿真用Matlab的Simulink进行仿真，在Simulink中新建一个文件将单位阶跃信号、负反 馈、开环传递函数和示波器连接起来，对位时随动系统进行仿真，即如

16、图3-1所示图3-1 Simulink仿真模型在Simulink中对单位阶跃信号的设置如图3-2所示图3-2单位阶跃信号在Simulink中对开环传递函数的设置如图3-3所示图3-3开环传递函数的设置将开环传递函数的参数分别设置为校正前后的，运行Simulink仿真可得校正前后的系统的单位阶跃反馈的曲线如图3-4所示图3-4校正前后系统的单位阶跃反馈曲线3.2系统校正前后的仿真分析由校正前后的仿真图比较可以看出，校正后的系统的上升时间明显变短，而上升时间 是系统响应速度的一种度量，上升时间越短，响应速度越快，所以校正后的系统的响应速 度较快。稍加分析不难得出上升时间变短的原因，因为校正后的阻尼

17、比变成了 0.75比校 正前的1.41要小，且从校正前的过阻尼系统变为校正后的欠阻尼系统，因此校正后的系 统响应的阻碍能力就相对减小，从而缩短了上升时间。其次，从图中还可以看出，系统的校正前无超调量，相对校正后的超调量很小，这也 是因为阻尼比减小且由过阻尼变为欠阻尼的缘故。心得体会经过此次为期近两星期的的自动控制理论课程设计，我对本课程又有了更进一步的认 识，特别是对于位置随动系统的测速反馈控制这一类系统来说。此次课程设计分为四个部分，一类位置随动系统的建模、传递函数的求解、测速反馈 校正装置的设计以及校正前后系统的仿真分析，只有把每一部分都做好才能保证本次课程 设计的成功；也就是说每一部分都需要认真分析。对位置测速反馈系统这类系统，我首先从书本上的介绍和例题入手，然后抽取其实分 析的套路，再对自己的课程设计进行分析建模，一步一步按照一般系统的分析方法做下来。 在用Matlab绘图的时候我又加深了对之前学过的信号与系统的实验课的理解和掌握，一举 两得。最终如愿完成了课程设计的所有任务。通过本次课程设计，我学会将课堂上所学到的理论知识运用到了实际问题中。并通过 对知识的综合利用，进行了一系列的分析，比较，提高