发明名称 单输入单输出时滞系统的分数阶PID控制方法
摘要 本发明公开了一种单输入单输出时滞系统的分数阶PID控制方法,步骤如下:1)利用工控系统对控制对象进行自动继电反馈辨识,得到被控对象的传递函数模型;2)适当地选取控制器参数积分阶次λ和微分阶次μ,使其满足条件0<λ<2和0<μ<2;3)对于给定的λ,μ,通过遍历一个控制参数,计算二维平面上关于另外两个控制参数的稳定域边界线;4)通过判断边界线的哪一侧具有更少的不稳定极点,确定被这些边界线所分割的哪个区域是控制参数的稳定域;5)建立分数阶PID控制器监控模块,通过鼠标点击图形界面上的控制参数稳定域选取控制参数,使其满足预定的性能指标值。本发明适用于具有任意传递函数模型的单输入单输出系统。
申请公布号 CN102073270B 申请公布日期 2015.08.05
申请号 CN201110034727.8 申请日期 2011.01.27
申请人 浙江工业大学 发明人 禹鑫燚;欧林林;周佩冬;陈宣光;冯远静;俞立
分类号 G05B11/42(2006.01)I 主分类号 G05B11/42(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵
主权项 一种单输入单输出时滞系统的分数阶PID控制方法,其特征在于,具体步骤如下:(1)在系统进入分数阶PID控制器的设计之前,先利用工控系统对控制对象进行自动继电反馈辨识,得到具有如下传递函数的被控对象模型<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>N</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>D</mi><mi>s</mi></mrow></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mi>&theta;s</mi></mrow></msup></mrow>]]></math><img file="FDA0000676929920000011.GIF" wi="348" he="144" /></maths>其中,<img file="FDA0000676929920000012.GIF" wi="671" he="136" />a<sub>i</sub>,b<sub>i</sub>,α<sub>i</sub>,β<sub>i</sub>,i=0,1,2,…,n为任意实数且满足β<sub>n</sub>&gt;…&gt;β<sub>1</sub>&gt;β<sub>0</sub>≥0,α<sub>n</sub>&gt;…&gt;α<sub>1</sub>&gt;α<sub>0</sub>≥0和α<sub>n</sub>&gt;β<sub>n</sub>;(2)适当地选取分数阶PID控制器C(s)=k<sub>p</sub>+k<sub>i</sub>/s<sup>λ</sup>+k<sub>d</sub>s<sup>μ</sup>中的积分阶次λ和微分阶次μ,使其满足条件0&lt;λ&lt;2和0&lt;μ&lt;2;(3)对于给定的λ,μ,根据被控对象的模型参数,通过遍历一个控制参数,计算二维平面上关于另外两个控制参数的稳定域边界线,参数空间的边界线由下述三部分组成:(a)实根边界(RRB):k<sub>i</sub>=0(b)无限边界(IRB):存在三种无限边界线当β<sub>n</sub>+μ&gt;α<sub>n</sub>时,IRB曲线为k<sub>d</sub>=0;当β<sub>n</sub>+μ=α<sub>n</sub>时,IRB曲线可以表示为‑a<sub>n</sub>/b<sub>n</sub>≤k<sub>d</sub>≤a<sub>n</sub>/b<sub>n</sub>;当β<sub>n</sub>+μ&lt;α<sub>n</sub>时,IRB曲线不存在;(c)复根边界(CRB):当λ+μ≠2时,采用3‑D法计算复根边界线;当λ+μ=2时,采用奇异频率法来计算复根边界线;(4)通过判断实根边界(RRB)、无限边界(IRB)和复根边界(CRB)的哪一侧具有更少的不稳定极点,确定被这些边界线所分割的哪个区域是控制参数的稳定域,令<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>-</mo><msub><mi>R</mi><mi>A</mi></msub><msub><mi>R</mi><mi>B</mi></msub><mo>-</mo><msub><mi>I</mi><mi>A</mi></msub><msub><mi>I</mi><mi>B</mi></msub></mrow><mrow><msup><mi>&omega;</mi><mi>&lambda;</mi></msup><mi>sin</mi><mrow><mo>(</mo><mi>&lambda;</mi><mfrac><mi>&pi;</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>R</mi><mi>A</mi></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>I</mi><mi>A</mi></msub><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>I</mi><mi>A</mi></msub><msub><mi>R</mi><mi>B</mi></msub><mo>-</mo><msub><mi>R</mi><mi>A</mi></msub><msub><mi>I</mi><mi>B</mi></msub></mrow><mrow><msup><mi>&omega;</mi><mi>&lambda;</mi></msup><mi>sin</mi><mrow><mo>(</mo><mi>&lambda;</mi><mfrac><mi>&pi;</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><msup><msub><mi>R</mi><mi>A</mi></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>I</mi><mi>A</mi></msub><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000676929920000021.GIF" wi="1326" he="209" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>g</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>&omega;</mi><mi>&lambda;</mi></msup><mi>sin</mi><mrow><mo>(</mo><mi>&lambda;</mi><mfrac><mi>&pi;</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>&omega;&theta;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>&omega;&theta;</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000676929920000022.GIF" wi="1078" he="122" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>h</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msup><mi>&omega;</mi><mi>&lambda;</mi></msup><msub><mi>k</mi><mi>p</mi></msub><mi>cos</mi><mrow><mo>(</mo><mi>&lambda;</mi><mfrac><mi>&pi;</mi><mn>2</mn></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000676929920000023.GIF" wi="509" he="128" /></maths>其中,R<sub>A</sub>和I<sub>A</sub>分别表示在s=jω时N(s)的实部和虚部,R<sub>B</sub>和I<sub>B</sub>分别表示在s=jω时s<sup>λ</sup>D(s)的实部和虚部,本发明根据以下规则确定稳定域位于边界线的哪一侧:(a)当λ+μ&lt;2时,沿着ω增大的方向,稳定域位于边界线的左侧;(b)当λ+μ&gt;2时,沿着ω增大的方向,稳定域位于边界线的右侧;(c)当λ+μ=2时,按以下规则判别:(I)若k<sub>i</sub>&gt;0,则k<sub>p</sub>&gt;[f<sub>1</sub>(ω)sin(ωθ)+f<sub>2</sub>cos(ωθ)]<sub>ω=0</sub>;(II)若k<sub>i</sub>&lt;0,则k<sub>p</sub>&lt;[f<sub>1</sub>(ω)sin(ωθ)+f<sub>2</sub>cos(ωθ)]<sub>ω=0</sub>;(III)当所有ω<sub>η</sub>∈Ω<sup>+</sup>时,k<sub>i</sub>&lt;ω<sub>η</sub><sup>2</sup>k<sub>d</sub>+g(ω<sub>η</sub>)+h(ω<sub>η</sub>);(Ⅳ)当所有ω<sub>η</sub>∈Ω<sup>‑</sup>时,k<sub>i</sub>&gt;ω<sub>η</sub><sup>2</sup>k<sub>d</sub>+g(ω<sub>η</sub>)+h(ω<sub>η</sub>);其中,Ω<sup>+</sup>={ω<sub>η</sub>∈R<sup>+</sup>|f(ω<sub>η</sub>)‑k<sub>p</sub>=0∧f'(ω<sub>η</sub>)&gt;0}Ω<sup>‑</sup>={ω<sub>η</sub>∈R<sup>+</sup>|f(ω<sub>η</sub>)‑k<sub>p</sub>=0∧f'(ω<sub>η</sub>)&lt;0}f(ω)=f<sub>1</sub>(ω)sin(ωθ)+f<sub>2</sub>cos(ωθ)且f′(ω)为f(ω)关于ω的导数;(5)基于步骤(3)和(4)中求解控制参数稳定域的算法和GUI开发软件在工控机中实现分数阶PID控制器的调节、仿真,建立分数阶PID控制器监控模块,用户界面能够进行被控对象模型的参数输入,超调、幅值裕度和相位裕度等性能指标的设置,控制参数稳定域和系统输出响应曲线的显示,并通过鼠标点击所获得的控制参数区域中各组不同的控制参数值,给出所对应的输出响应曲线及系统各性能指标的值;若不满足预定的性能指标值,则重新选取控制参数,若满足,则可切换到在线控制状态,直接对被控对象进行控制。
地址 310014 浙江省杭州市潮王路18号