发明名称 基于期望闭环传递函数的伺服系统控制器优化方法
摘要 基于期望闭环传递函数的伺服系统控制器优化方法,本发明涉及基于期望闭环传递函数的伺服系统控制器优化方法。本发明的目的是为了解决目前提供的控制器没有建立频响指标和控制器参数之间的定量关系,设计结果的好坏很大程度上取决于设计者的经验,不能保证闭环系统性能是最优的,对于双十控制器优化方法较少报道。基于期望闭环传递函数的伺服系统控制器优化方法具体为:步骤一、确定被控对象的传递函数;步骤二、构造期望闭环传递函数的结构形式;步骤三、建立系统频响指标、剪切频率和稳定裕度约束条件与期望闭环传递函数参数之间的数学关系;步骤四、得到优化的期望闭环传递函数的参数;步骤五、求解出控制器。本发明属于控制器优化技术领域。
申请公布号 CN104181905B 申请公布日期 2016.09.14
申请号 CN201410471644.9 申请日期 2014.09.16
申请人 哈尔滨工业大学 发明人 陈松林;陈婷;李明
分类号 G05B19/418(2006.01)I 主分类号 G05B19/418(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 基于期望闭环传递函数的伺服系统控制器优化方法,其特征在于:基于期望闭环传递函数的伺服系统控制器优化方法具体是按照以下步骤进行的:步骤一、系统辨识,确定被控对象的传递函数;被控对象的传递函数选取为<img file="FDA0001025633160000011.GIF" wi="547" he="135" />式中,k<sub>G</sub>为系统增益,τ<sub>e</sub>为电气时间常数,τ<sub>m</sub>为机电时间常数,s为Laplace算子;通过代入标称参数或者参数辨识的方法测得实际被控对象传递函数中的参数k<sub>G</sub>、τ<sub>e</sub>和τ<sub>m</sub>,从而确定被控对象的传递函数;步骤二、构造期望闭环传递函数的结构形式;期望闭环传递函数的结构形式为<img file="FDA0001025633160000012.GIF" wi="731" he="150" />式中Φ<sub>E</sub>(s)为构造的期望闭环传递函数,ξ为阻尼比,w<sub>n</sub>为自然角频率,s为Laplace算子,T为辅助惯性环节时间常数;其中,ξ和w<sub>n</sub>为待优化的参数,而T用于保证控制器的物理可实现性,将其取为常值,通常其取值要小于<img file="FDA0001025633160000013.GIF" wi="136" he="127" />w<sub>c</sub>为剪切频率;步骤三、建立系统频响指标、剪切频率和稳定裕度约束条件与期望闭环传递函数参数之间的数学关系;分别由下列解析关系式给出:1)双十频响指标i)幅值误差对于期望闭环传递函数Φ<sub>E</sub>(s),定义w<sub>A1.1</sub>为其幅频特性首次大于1.1的频率点,定义w<sub>A0.9</sub>为其幅频特性首次小于0.9的频率点;对于期望闭环传递函数Φ<sub>E</sub>(s),随参数阻尼比ξ,自然角频率w<sub>n</sub>,辅助惯性时间常数T的变化,从单调性方面考虑幅频特性会出现先增后减、先减后增再减以及单调递减三种情况,相应地,幅值误差首次超过10%的频率点可分以下2种情况进行求取,一种为由w<sub>A1.1</sub>决定的情况,另一种为由w<sub>A0.9</sub>决定的情况;首先给出几个关键判别式的定义a.<img file="FDA0001025633160000021.GIF" wi="910" he="234" />其中,<img file="FDA0001025633160000022.GIF" wi="899" he="507" />b.<img file="FDA0001025633160000023.GIF" wi="1038" he="239" />其中,<img file="FDA0001025633160000024.GIF" wi="362" he="150" />c.Δ<sub>A0.9</sub>=B<sub>s</sub><sup>2</sup>‑4A<sub>s</sub>C<sub>s</sub>,其中,<img file="FDA0001025633160000025.GIF" wi="859" he="502" />和<img file="FDA0001025633160000026.GIF" wi="433" he="485" />A.当A<sub>min</sub>≤0.9;A<sub>min</sub>&gt;0.9且M<sub>r</sub>&lt;1.1;<img file="FDA0001025633160000027.GIF" wi="893" he="143" /><img file="FDA0001025633160000028.GIF" wi="614" he="486" />或<img file="FDA0001025633160000029.GIF" wi="632" he="614" />时,幅值误差小于10%的最高频率点由幅值首次变为0.9时的频率w<sub>A0.9</sub>决定,具体解析表达式分三种情况,分别为:a)若<img file="FDA00010256331600000210.GIF" wi="206" he="150" />或<img file="FDA00010256331600000211.GIF" wi="246" he="151" />则<img file="FDA00010256331600000212.GIF" wi="827" he="206" />b)若Δ<sub>A0.9</sub>≥0,A<sub>s</sub>&lt;0,则<img file="FDA0001025633160000031.GIF" wi="819" he="205" />c)若Δ<sub>A0.9</sub>&lt;0,A<sub>min</sub>≤0.9,则<img file="FDA0001025633160000032.GIF" wi="625" he="220" />B.A<sub>min</sub>&gt;0.9且M<sub>r</sub>≥1.1时,其幅值误差小于10%的最高频率点由幅值首次变为1.1时的频率点w<sub>A1.1</sub>决定,计算公式为:<img file="FDA0001025633160000033.GIF" wi="846" he="203" />ii)相角误差对于期望闭环传递函数Φ<sub>E</sub>(s),其相频特性的数学表达式为<img file="FDA0001025633160000034.GIF" wi="750" he="135" />由于其具有单调递减特性,只定义相位首次变为‑10°的频率为w<sub>P10</sub>;由于w<sub>P10</sub>&lt;&lt;w<sub>n</sub>,忽略w<sup>2</sup>,令<img file="FDA0001025633160000035.GIF" wi="331" he="70" />得w<sub>P10</sub>的解析表达式为<maths num="0001"><math><![CDATA[<mrow><msub><mi>w</mi><mrow><mi>P</mi><mn>10</mn></mrow></msub><mo>=</mo><mfrac><mrow><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>Tw</mi><mi>n</mi></msub><mo>+</mo><mn>2</mn><mi>&xi;</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>8</mn><msub><mi>T&xi;w</mi><mi>n</mi></msub><msup><mi>tan</mi><mn>2</mn></msup><mrow><mo>(</mo><mfrac><mi>&pi;</mi><mn>18</mn></mfrac><mo>)</mo></mrow></mrow></msqrt><mo>-</mo><mrow><mo>(</mo><msub><mi>Tw</mi><mi>n</mi></msub><mo>+</mo><mn>2</mn><mi>&xi;</mi><mo>)</mo></mrow></mrow><mrow><mn>4</mn><mi>T</mi><mi>&xi;</mi><mi>tan</mi><mrow><mo>(</mo><mfrac><mi>&pi;</mi><mn>18</mn></mfrac><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0001025633160000036.GIF" wi="1006" he="271" /></maths>综合幅值误差与相角误差的分析结果可知,若要满足双十频响指标,则必须满足min{w<sub>A1.1</sub>,w<sub>P10</sub>}≥w<sub>fr</sub>或min{w<sub>A0.9</sub>,w<sub>P10</sub>}≥w<sub>fr</sub>,w<sub>fr</sub>为满足双十频响指标要求的最高频率;2)剪切频率剪切频率也是控制系统的约束条件之一,下面给出了剪切频率与期望闭环传递函数参数之间的数学关系:定义<img file="FDA0001025633160000037.GIF" wi="590" he="335" />令<img file="FDA0001025633160000038.GIF" wi="451" he="479" />并定义判别式<img file="FDA0001025633160000039.GIF" wi="355" he="67" />则剪切频率w<sub>c</sub>的解析表达式分为三种情况,分别为:i)若Δ<sub>c</sub>≥0,A<sub>c</sub>≥0,则<img file="FDA0001025633160000041.GIF" wi="579" he="221" />ii)若Δ<sub>c</sub>≥0,A<sub>c</sub>&lt;0,则<img file="FDA0001025633160000042.GIF" wi="819" he="209" />iii)若Δ<sub>c</sub>&lt;0,则<img file="FDA0001025633160000043.GIF" wi="843" he="221" />3)稳定裕度i)相位裕度根据相位裕度的定义,代入剪切频率w<sub>c</sub>的表达式,即可得到相位裕度γ与期望闭环传递函数参数的关系为:<img file="FDA0001025633160000044.GIF" wi="670" he="143" />ii)幅值裕度根据幅值裕度的物理意义,并利用劳斯稳定判据,可得幅值裕度k<sub>g</sub>与期望闭环传递函数参数的解析表达式为:<maths num="0002"><math><![CDATA[<mrow><msub><mi>k</mi><mi>g</mi></msub><mo>=</mo><mfrac><mrow><mo>(</mo><mn>2</mn><msub><mi>T&xi;w</mi><mi>n</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo><mo>(</mo><msub><mi>Tw</mi><mi>n</mi></msub><mo>+</mo><mn>2</mn><mi>&xi;</mi><mo>)</mo></mrow><mrow><msub><mi>Tw</mi><mi>n</mi></msub></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0001025633160000045.GIF" wi="590" he="142" /></maths>步骤四、运用MATLAB工具箱求解优化问题,得到优化的期望闭环传递函数的参数;步骤五、运用被控对象传递函数及期望闭环传递函数求解出控制器。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号