发明名称 基于向量方法的分数阶PI<sup>λ</sup>控制器的参数整定方法
摘要 本发明涉及一种基于向量方法的分数阶PI<sup>λ</sup>控制器的参数整定方法,属于分数阶自动控制技术领域。给定被控对象系统的数学模型传递函数,将被控对象写成被控对象向量形式,写出性能指标向量形式,利用一个校正向量L去校正被控对象向量P达到性能指标向量G,校正向量和控制器向量在复平面内构成三角形,求出控制器传递函数。优点是:可以减少计算量,而且整定过程简单易懂。当被控对象发生变化时,我们只需要计算出被控对象的模值、相角和角速度,后面的公式可以套用,从而减少了反复推导方程的过程,解决了matlab求解中存在的多解问题。
申请公布号 CN103631135B 申请公布日期 2016.01.20
申请号 CN201310615818.X 申请日期 2013.11.28
申请人 长春理工大学 发明人 王春阳;蔡年春;尹高峰;樊斌;辛瑞昊
分类号 G05B11/42(2006.01)I 主分类号 G05B11/42(2006.01)I
代理机构 吉林长春新纪元专利代理有限责任公司 22100 代理人 魏征骥
主权项 一种基于向量方法的分数阶PI<sup>λ</sup>控制器的参数整定方法,其特征在于包括下列步骤:(一)、给定被控对象系统的数学模型传递函数<img file="FDA0000838310880000011.GIF" wi="302" he="123" />并给定设计指标带宽ω<sub>c</sub>和需保持稳定的相位裕度φ<sub>m</sub>,待整定控制器传递函数形式<img file="FDA0000838310880000012.GIF" wi="352" he="122" />其中,所述T为正实数,s为拉普拉斯算子,K<sub>p</sub>表示待整定的比例系数,K<sub>i</sub>表示待整定的积分系数,λ表示待整定的积分阶次;(二)、将被控对象<img file="FDA0000838310880000013.GIF" wi="286" he="137" />写成被控对象向量形式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>j&omega;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msqrt><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></msqrt></mfrac><mo>&angle;</mo><mo>-</mo><msub><mi>arctanT&omega;</mi><mi>c</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000838310880000014.GIF" wi="789" he="170" /></maths>被控对象向量角速度<maths num="0002" id="cmaths0002"><math><![CDATA[<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><img file="FDA0000838310880000015.GIF" wi="393" he="151" /></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="FDA0000838310880000016.GIF" wi="955" he="109" />同时求得校正向量在ω<sub>c</sub>处的角速度<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&psi;</mi><mi>L</mi></msub><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><img file="FDA0000838310880000017.GIF" wi="357" he="156" /></maths>(五)将控制器<img file="FDA0000838310880000018.GIF" wi="330" he="125" />映射到复平面内得到由相角为0且模值为K<sub>p</sub>的向量加上一个相角为<img file="FDA0000838310880000019.GIF" wi="134" he="117" />且模值为<img file="FDA00008383108800000110.GIF" wi="82" he="126" />的向量的和,其中K<sub>p</sub>∠0°、<img file="FDA00008383108800000111.GIF" wi="277" he="131" />即控制器向量,利用控制器向量C去逼近校正向量L,即lim|C‑L|=0,从而校正向量和控制器向量在复平面内构成三角形;(六)令<img file="FDA00008383108800000112.GIF" wi="1104" he="108" />即<img file="FDA00008383108800000113.GIF" wi="330" he="108" />和θ=φ<sub>m</sub>+arctanTω<sub>c</sub>‑180°,利用三角形的余弦定理可以得到:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfrac><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&lambda;</mi></msubsup></mfrac><mo>=</mo><msqrt><mrow><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>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00008383108800000114.GIF" wi="1238" he="142" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mi>&pi;</mi><mo>-</mo><mfrac><mi>&pi;</mi><mn>2</mn></mfrac><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>p</mi></msub><mo>-</mo><mi>A</mi><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi></mrow><msqrt><mrow><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>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi></mrow></msqrt></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000838310880000021.GIF" wi="1231" he="173" /></maths>同时利用lim|C‑L|=0角速度相等,可得:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfrac><mrow><msub><mi>&lambda;K</mi><mi>p</mi></msub><mfrac><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mrow><mi>&lambda;</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mfrac><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>&lambda;</mi><mi>&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><mfrac><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&lambda;</mi></msubsup></mfrac><mi>c</mi><mi>o</mi><mi>s</mi><mo>(</mo><mi>&lambda;</mi><mi>&pi;</mi><mo>/</mo><mn>2</mn><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&lambda;</mi></msubsup></mfrac><mi>s</mi><mi>i</mi><mi>n</mi><mo>(</mo><mi>&lambda;</mi><mi>&pi;</mi><mo>/</mo><mn>2</mn><mo>)</mo><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><img file="FDA0000838310880000022.GIF" wi="1256" he="300" /></maths>注意到存在几何关系为<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><msub><mi>K</mi><mi>p</mi></msub><mo>+</mo><mfrac><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&lambda;</mi></msubsup></mfrac><mi>c</mi><mi>o</mi><mi>s</mi><mo>(</mo><mrow><mi>&lambda;</mi><mi>&pi;</mi><mo>/</mo><mn>2</mn></mrow><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&lambda;</mi></msubsup></mfrac><mi>s</mi><mi>i</mi><mi>n</mi><mo>(</mo><mrow><mi>&lambda;</mi><mi>&pi;</mi><mo>/</mo><mn>2</mn></mrow><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><mi>A</mi><mn>2</mn></msup></mrow>]]></math><img file="FDA0000838310880000023.GIF" wi="1011" he="175" /></maths>和<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfrac><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>&omega;</mi><mi>c</mi><mi>&lambda;</mi></msubsup></mfrac><mi>sin</mi><mrow><mo>(</mo><mi>&lambda;</mi><mi>&pi;</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>A</mi><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi><mo>,</mo></mrow>]]></math><img file="FDA0000838310880000024.GIF" wi="556" he="143" /></maths>所以可以得到,<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mo>-</mo><mfrac><mrow><msub><mi>&lambda;K</mi><mi>p</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi></mrow><mrow><msub><mi>A&omega;</mi><mi>c</mi></msub></mrow></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><img file="FDA0000838310880000025.GIF" wi="1261" he="142" /></maths>利用式(2)和(4)可求得K<sub>p</sub>和λ,代入式(1)可求得K<sub>i</sub>,求出控制器传递函数<img file="FDA0000838310880000026.GIF" wi="340" he="124" />即完成了一种鲁棒性的分数阶PI<sup>λ</sup>控制器参数整定。
地址 130022 吉林省长春市卫星路7989号