天线方位角位置随动系统的建模与分析

2-5b天线方位角位置随动系统建摸系统的原理图如图2-7所示，其方块图如图2-8所示。系统的任务是使输出的天线方位角θ0(t)跟踪输入方位角θi(t)的变化，试建立该系统的数学模型。系统的参数值如下：电源电压V=10v；功率放大器的增益和时间常数K1=1，T1=0.01s；伺服电动机的电枢回路电阻Rd=8Ω，转动惯量Ja=0.02Kgm2，粘性摩擦系数fa=0.01Nms/rad，反电势系数Ce=0.5Vs/rad，转矩系数Cm=0.5Nm/A；减速器各齿轮的齿数为Z1=25，Z2=Z3=250；负载端的转动惯量JL=1Kgm2粘性摩擦系数fL=1Nms/rad。解：采用组合系统建摸法，根据原理图2-7可以将系统划分为六个环节：输入电位器，差分放大器，功率放大器，电动机，减速器和输出电位器。首先建立各个环节的数学模型，然后将它们组合起来则可得系统的数学摸型。1环节的数学模型(1)输入电位器与输出电位器由于输入电位器与输出电位器的线路和电位器的结构均相同，故这两个环节的传递函数是一样的。对电位器环节的输出电压与输入角位移的特性进行线性化处理则可视其为一比例环节。由图2-7可知；当动触头位于电位器中心时其输出电压为零；朝前或朝后转动5圈其输出电压变化均为10V。于是可得它们的传递函数为00()()100.318/()52()ipotiususkvradss(2)差分放大器与功率放大器放大器通常工作在放大状态，可不考虑饱和的影响。差分放大器的时间常数比起功率放大器以及系统的其他环节的时间常数要小得多，可以忽视不计。故这两个环节的输入输出传递函数分别为差分放大器ceCKsUsU)()(功率放大器1()1()10.011dociuskusTss(3)电动机在小功率伺服系统中直流电动机的结构图中，由于电动机的电枢回路电感很小，可以忽略不计。图中的J与f为折算至电动机轴上系统转动部分的等效转动惯量和等效粘性摩擦系数，其值分别为2212(/)0.021(25/250)0.03aLJJJZZ2212(/)0.011(25/250)0.02aLfffZZ将具体参数值代入，于是可求得电动机的电枢（空载）电压与转子角位移之间的传递函数为图2-7天线方位角位置随动系统原理图()/2.083()[()][1/(/)](1.708)mmmddodsememdsCCJRussRJfCCssJfCCRss(4)减速器齿轮减速器的传动比为21/izz=250/25=10，于是减速器的传递函数为()1/0.1()omsis2系统的输入输出模型将个环节的数学模型按照信号的传递关系组合起来，则可绘制系统的结构图如图2-8所示。应用梅森公式或结构图化简，由图则可求得系统的传递函数为32()6.62()101.71170.86.62ocicskssssk3-16c（天线方位角位置控制系统的时域分析）在2-5b已建立该系统的数学模型，其结构如图2-7所示，若系统采用比例控制（即其前置放大器的增益cK是可调的），试求（1）当cK＝1000时系统的时间响应特性；（2）若要求超调量20%p，cK应调整为何值？并分别应用MATLAB和近似估算的方法求系统的暂态性能，说明近似估算的适用条件。解：（1）增益cK的稳定取值范围天下女子面目可憎言语无味由系统的传递函数32()6.62()101.71170.86.62ocicsKssssK，可得系统的特征方程为32det()101.71170.86.620csIAsssK于是可构造劳斯表如下：320S1170.8S101.716.62S17372.076.620S6.62cccKKK根据劳斯判据，为确保系统稳定必须使劳斯表第一列的元素不变号，即图2-8天线方位角位置随动系统结构图θi(s)θ0(s)θm(s)减速器17372.07-6.62Kc06.62Kc0故可求得系统在比例控制时，放大器增益的稳定取值范围为0Kc2624.18。（2）当Kc＝1000时系统暂态特性的分析由题2-5b系统结构图2-7可得，系统的开、闭环传递函数为0.06620.0388()(0.011)(1.708)(0.011)(0.591)cckKKGsssssss32()6.62()()101.71170.86.62ocicsKsssssK由于Kc远低于其临界值（Ker＝2624.18），在系统的3个开环极点中功率放大器的极点Po3＝-1/0.01＝-100远离虚轴（相应的时间常数T3=0.01很小）可忽略不计，故系统的开、闭环传递函数可简化为22220.06620.06620.0662()'()(0.011)(1.708)(1.708)(1.708)'()66.2()'()1'()1.70866.22cckkknknnKKGsGssssssssGsssGsssss（3-9）式中66.28.14/nrads，1.708/(2)0.105n。于是可按二阶规范系统的表达式，估算系统的暂态性能如下：20.391pnts%77.71%100*21/ep44.68snts（3）当要求超调量20%p时系统暂态特性的分析(a)暂态性能的估算由于系统响应的超调量远低于第（2）项的相应值（71.77％），可以预计这时的增益Kc1000，故开环时间常数T3＝0.01可忽略不计。于是系统的开环传递函数可简化为'()0.0662/[(1.708)]kcGsKss，从而系统的闭环传递函数可近似为2222'()66.2'()1'()1.70866.22kcnkcnnGsKsGsssKss（3-10）根据对超调量的要求2/1100%0.2pe，则可得系统的阻尼比为22ln(1/)0.4558[ln(1/)]pp而21.708ns，20.0662ncK。故可求得1.708/(2)1.8736/nrads，2/0.066253.03cnK。于是根据二阶规范系统的性能指标表达式，则可估算系统的暂态性能如下：2'1.881pnts'20%p4'4.68snts（b）应用MATLAB进行时域分析由式（3-9）和(3-10)应用MATLAB程序A3.1，则可求得系统的实际暂态性能为pt1.9s，p＝20.8546％，st＝4.4（s）；并绘制系统的单位阶跃响应曲线，如图3-20所示。图中实线为按式（3-9）求得的准确相应曲线，用”.”标出数据点的为按式（3-10）将系统视为二阶系统的近似响应曲线。由图中可见：当Kc远比其临界值小得多时，忽略小环节时间常数的影响并进行近似估算，在工程上是可行的。MATLAB程序A3.1G1=tf(0.0662*53.03,[1,1.708,0.0662*53.03]);G=tf(6.62*53.03,[1,101.71,170.8,6.62*53.03]);T=0:0.1:8;y=step(G,t);y1=step(G1,t);plot(t,y,’-‘,t,y1,’.’),gridXlabel(‘t(sec)’)Legend(‘y’,’y1’)Y=dcgain(G);[yp,p]=max(y);peaktime=t(p)Percentovershoot=100*(yp-Y)/YK=length(t);while(y(k)0.98*Y&(y(k)1.02*Y),k=k-1;endSettingtime=t(k)(4)稳态特性分析由式（3-9）可知：系统为1型的，且开环增益K=0.0388Kc＝38.8（当Kc＝1000时）或2.06（当Kc＝53.03）。于是系统的误差系数为pKvKK0aK。因此该系统跟踪阶跃输入信号没有稳态误差；跟踪单位斜坡输入信号的稳态误差为1/38.80.026100011/2.060.48553.03cscKeKK当＝时当＝时但无法跟踪抛物线函数或其它变化更快的输入信号，若勉强使用则跟踪误差将随时间而不断地增大，当t时其终值误差将趋于无穷大。图3-20Kc=53.03时系统的单位阶跃响应曲线