发明名称 基于向量方法的鲁棒分数阶PD<sup>μ</sup>控制器的参数整定方法
摘要 本发明涉及一种基于向量方法的鲁棒分数阶PD<sup>μ</sup>控制器的参数整定方法,属于分数阶自动控制技术领域。给定被控对象系统的数学模型传递函数,将被控对象写成被控对象向量形式,写出性能指标向量形式,利用一个校正向量L去校正被控对象向量P达到性能指标向量G,校正向量和控制器向量在复平面内构成三角形,求得待整定的微分系数。本方法可以减少计算量,而且整定过程简单易懂,当被控对象发生变化时,我们只需要计算出被控对象的模值、相角和角速度,后面的公式可以套用,从而减少了反复推导方程的过程,解决了matlab求解中存在的多解问题。
申请公布号 CN103605290A 申请公布日期 2014.02.26
申请号 CN201310624499.9 申请日期 2013.11.28
申请人 长春理工大学 发明人 王春阳;蔡年春;辛瑞昊;樊斌;尹高峰
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 吉林长春新纪元专利代理有限责任公司 22100 代理人 魏征骥
主权项 1.一种基于向量方法的鲁棒分数阶PD<sup>μ</sup>控制器的参数整定方法,其特征在于包括下列步骤:(一)、给定被控对象系统的数学模型传递函数<img file="FDA0000423866640000011.GIF" wi="375" he="146" />并给定设计指标带宽ω<sub>c</sub>和需保持稳定的相位裕度φ<sub>m</sub>,待整定控制器传递函数形式C(s)=K<sub>p</sub>+K<sub>d</sub>s<sup>μ</sup>;其中,所述T为正实数,s为拉普拉斯算子,K<sub>p</sub>表示待整定的比例系数,K<sub>d</sub>表示待整定的微分系数,μ表示待整定的微分阶次;(二)、将被控对象<img file="FDA0000423866640000012.GIF" wi="349" he="144" />写成被控对象向量形式:<img file="FDA0000423866640000013.GIF" wi="975" he="164" />被控对象向量角速度<maths num="0001"><![CDATA[<math><mrow><msub><mi>&psi;</mi><mi>P</mi></msub><mo>=</mo><mo>-</mo><mfrac><mi>T</mi><mrow><msup><mrow><mo>(</mo><msub><mi>T&omega;</mi><mi>c</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>1</mn></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>(三)、将设计指标带宽ω<sub>c</sub>和需保持稳定的相位裕度φ<sub>m</sub>,写出性能指标向量形式:G(jω<sub>c</sub>)=1∠φ<sub>m</sub>-180°,角速度为0;(四)、利用一个校正向量L去校正被控对象向量P达到性能指标向量G,即L·P=G,根据步骤(二)和步骤(三),可以求得在ω<sub>c</sub>处校正向量L,<img file="FDA0000423866640000015.GIF" wi="983" he="109" />同时求得校正向量在ω<sub>c</sub>处的角速度<maths num="0002"><![CDATA[<math><mrow><msub><mi>&psi;</mi><mi>L</mi></msub><mo>=</mo><mo>-</mo><mfrac><mi>T</mi><mrow><msup><mrow><mo>(</mo><msub><mi>T&omega;</mi><mi>c</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>1</mn></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>(五)、将控制器C(s)=K<sub>p</sub>+K<sub>d</sub>s<sup>μ</sup>映射到复平面内得到由相角为0且模值为K<sub>p</sub>的向量(K<sub>p</sub>∠0°)加上一个相角为<img file="FDA0000423866640000017.GIF" wi="57" he="116" />μ且模值为K<sub>d</sub>ω<sup>μ</sup>的向量<img file="FDA0000423866640000018.GIF" wi="320" he="114" />的和,即控制器向量。利用控制器向量C去逼近校正向量L,即lim|C-L|=0,从而校正向量和控制器向量在复平面内构成三角形;(六)、令<img file="FDA0000423866640000019.GIF" wi="1143" he="108" />即<img file="FDA0000423866640000021.GIF" wi="389" he="107" />和θ=φ<sub>m</sub>+arctanTω<sub>c</sub>-90°,利用三角形的余弦定理可以得到:<maths num="0003"><![CDATA[<math><mrow><msub><mi>K</mi><mi>d</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&mu;</mi></msubsup><mo>=</mo><msqrt><msubsup><mi>K</mi><mi>p</mi><mn>2</mn></msubsup><mo>+</mo><msup><mi>A</mi><mn>2</mn></msup><mo>-</mo><mn>2</mn><msub><mi>K</mi><mi>p</mi></msub><mi>A</mi><mi>cos</mi><mi>&theta;</mi></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>cos</mi><mrow><mo>(</mo><mi>&pi;</mi><mo>-</mo><mfrac><mi>&pi;</mi><mn>2</mn></mfrac><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>p</mi></msub><mo>-</mo><mi>A</mi><mi>cos</mi></mrow><msqrt><msubsup><mi>K</mi><mi>p</mi><mn>2</mn></msubsup><mo>+</mo><msup><mi>A</mi><mn>2</mn></msup><mo>-</mo><mn>2</mn><msub><mi>K</mi><mi>p</mi></msub><mi>A</mi><mi>cos</mi><mi>&theta;</mi></msqrt></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>同时利用lim|C-L|=0角速度相等,可得:<maths num="0005"><![CDATA[<math><mrow><mfrac><mrow><mi>&mu;</mi><msub><mi>K</mi><mi>p</mi></msub><msub><mi>K</mi><mi>d</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mrow><mi>&mu;</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mi>sin</mi><mrow><mo>(</mo><mi>&mu;&pi;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow></mrow><mrow><msup><mrow><mo>(</mo><msub><mi>K</mi><mi>p</mi></msub><mo>+</mo><msub><mi>K</mi><mi>d</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&mu;</mi></msubsup><mi>cos</mi><mrow><mo>(</mo><mi>&mu;&pi;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>K</mi><mi>d</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&mu;</mi></msubsup><mi>sin</mi><mrow><mo>(</mo><mi>&mu;&pi;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>=</mo><msub><mi>&psi;</mi><mi>L</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>注意到存在几何关系为<maths num="0006"><![CDATA[<math><mrow><msup><mrow><mo>(</mo><msub><mi>K</mi><mi>p</mi></msub><mo>+</mo><msub><mi>K</mi><mi>d</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&mu;</mi></msubsup><mi>cos</mi><mrow><mo>(</mo><mi>&mu;&pi;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>K</mi><mi>d</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&mu;</mi></msubsup><mi>sin</mi><mrow><mo>(</mo><mi>&mu;&pi;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><mi>A</mi><mn>2</mn></msup></mrow></math>]]></maths>和<maths num="0007"><![CDATA[<math><mrow><msub><mi>K</mi><mi>d</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&mu;</mi></msubsup><mi>sin</mi><mrow><mo>(</mo><mi>&mu;&pi;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mi>sin</mi><mi>&theta;</mi><mo>,</mo></mrow></math>]]></maths>所以可以得到,<maths num="0008"><![CDATA[<math><mrow><mfrac><mrow><mi>&mu;</mi><msub><mi>K</mi><mi>p</mi></msub><mi>sin</mi><mi>&theta;</mi></mrow><msub><mi>A&omega;</mi><mi>c</mi></msub></mfrac><mo>=</mo><msub><mi>&psi;</mi><mi>L</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>利用式(2)和(4)可求得K<sub>p</sub>和μ,代入式(1)可求得K<sub>d</sub>,即完成了一种鲁棒性的分数阶控制器参数整定。
地址 130022 吉林省长春市卫星路7989号