发明名称 一种增强型的增益鲁棒分数阶PID控制器参数整定方法
摘要 本发明涉及一种增强型的增益鲁棒分数阶PID控制器参数整定方法,属于电机位置系统的分数阶自动控制技术领域。电机系统控制器在拥有三个可调参数K<sub>p</sub>,K<sub>i</sub>,K<sub>d</sub>前提下,又添加了两个参数λ,μ,这五个参数配合调整,可大大增加系统的稳定性及动态响应,同时由于增加了约束条件Arg[G(jω<sub>b</sub>)]=Arg[C(jω<sub>b</sub>)P(jω<sub>b</sub>)]=-π+φ<sub>m</sub>和Arg[G(jω<sub>h</sub>)]=Arg[C(jω<sub>h</sub>)P(jω<sub>h</sub>)]=-π+φ<sub>m</sub>,还可以满足增益K<sub>1</sub>变化较大情况下的系统鲁棒性需求。对电机系统控制精度及鲁棒性得到了很大的提升。
申请公布号 CN104777746A 申请公布日期 2015.07.15
申请号 CN201510166587.8 申请日期 2015.04.09
申请人 长春理工大学 发明人 王春阳;刘雪莲;吴萌;蔡年春;辛瑞昊;李明秋;田成军
分类号 G05B11/42(2006.01)I;H02P29/00(2006.01)I 主分类号 G05B11/42(2006.01)I
代理机构 吉林长春新纪元专利代理有限责任公司 22100 代理人 魏征骥
主权项 一种增强型的增益鲁棒分数阶PID控制器参数整定方法,电机系统位置控制模型的传递函数形如:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>K</mi><mn>1</mn></msub><mrow><mi>Js</mi><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mi>s</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>T</mi><mn>2</mn></msub><mi>s</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000696417260000011.GIF" wi="1236" he="187" /></maths>其中J为有效负荷的惯性参数;T<sub>1</sub>,T<sub>2</sub>为系统的时间常数,s为拉普拉斯算子;K<sub>1</sub>为增益常数;其特征在于所述电机系统位置鲁棒控制器的参数整定方法,包括以下步骤:(1)对于电动位置系统被控对象的数学模型传递函数P(s),其待整定FOPID控制器传递函数形式<img file="FDA0000696417260000012.GIF" wi="501" he="128" />待整定参数为比例系数K<sub>p</sub>,积分系数K<sub>i</sub>,微分系数K<sub>d</sub>,积分阶次λ,微分阶次μ,并给定需校正穿越频率ω<sub>c</sub>和需保持稳定的相位裕度φ<sub>m</sub>;(2)利用MATLAB画出被控对象P(s)的伯德图,求得在频率ω<sub>c</sub>处的模值|P(jω<sub>c</sub>)|,及相角Arg[P(jω<sub>c</sub>)],同时给定ω<sub>c</sub>周围两个频率点ω<sub>b</sub>,ω<sub>h</sub>,其中ω<sub>b</sub>,ω<sub>h</sub>经验值取值范围30%ω<sub>c</sub>≤ω<sub>b</sub>&lt;ω<sub>c</sub>&lt;ω<sub>h</sub>≤350%ω<sub>c</sub>,并利用MATLAB所画出被控对象P(s)的伯德图,求得在频率ω<sub>b</sub>处的相角Arg[P(jω<sub>b</sub>)]及ω<sub>h</sub>处的相角Arg[P(jω<sub>h</sub>)];(3)将增益鲁棒性条件作为目标函数,令:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>V</mi><mo>=</mo><mfrac><mrow><mi>dArg</mi><mo>[</mo><mi>G</mi><mo>]</mo><mrow><mo>(</mo><mi>j&omega;</mi><mo>)</mo></mrow></mrow><mi>d&omega;</mi></mfrac><msub><mo>|</mo><mrow><mi>&omega;</mi><mo>=</mo><msub><mi>&omega;</mi><mi>c</mi></msub></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000696417260000013.GIF" wi="975" he="157" /></maths>其中G(jω)=C(jω)P(jω);(4)利用幅值裕度稳定性条件,根据穿越频率的定义,开环系统传递函数G(jω<sub>c</sub>)=C(jω<sub>c</sub>)P(jω<sub>c</sub>)在穿越频率ω<sub>c</sub>处幅值为1,于是可以得到:|G(jω<sub>c</sub>)|=|C(jω<sub>c</sub>)P(jω<sub>c</sub>)|=1     (3)将(3)等式转换为不等式,可得:||C(jω<sub>c</sub>)|‑1/|P(jω<sub>c</sub>)||≤0     (4)(5)利用相位裕度稳定性条件,系统开环传递函数在穿越频率ω<sub>c</sub>处相位裕度为φ<sub>m</sub>;Arg[G(jω<sub>c</sub>)]=Arg[C(jω<sub>c</sub>)P(jω<sub>c</sub>)]=‑π+φ<sub>m</sub>     (5)其中Arg[C(jω<sub>c</sub>)P(jω<sub>c</sub>)]=Arg[C(jω<sub>c</sub>)]+Arg[P(jω<sub>c</sub>)];将(5)等式转换为不等式,可得|Arg[C(jω<sub>c</sub>)P(jω<sub>c</sub>)]+π‑φ<sub>m</sub>|≤0     (6)(6)为扩大系统开环相频特性的平坦范围,添加两个约束条件Arg[G(jω<sub>b</sub>)]=Arg[C(jω<sub>b</sub>)P(jω<sub>b</sub>)]=‑π+φ<sub>m</sub>     (7)Arg[G(jω<sub>h</sub>)]=Arg[C(jω<sub>h</sub>)P(jω<sub>h</sub>)]=‑π+φ<sub>m</sub>     (8)将(7)(8)等式转换为不等式,可得|Arg[C(jω<sub>b</sub>)P(jω<sub>b</sub>)]+π‑φ<sub>m</sub>|≤0     (9)|Arg[C(jω<sub>h</sub>)P(jω<sub>h</sub>)]+π‑φ<sub>m</sub>|≤0     (10)(7)利用MATLAB中的Fmincon函数工具箱,将(2)作为目标函数,(4)(6)(9)(10)作为约束条件,可求解出增强型分数阶控制器五个参数最优解,即得到K<sub>p</sub>,K<sub>i</sub>,K<sub>d</sub>,λ,μ;(8)将上述求出的五个参数带入<img file="FDA0000696417260000021.GIF" wi="500" he="133" />即完成控制器参数整定。
地址 130022 吉林省长春市卫星路7989号