发明名称 基于预测函数控制优化的非自衡对象控制方法
摘要 本发明公开了一种基于预测函数控制优化的非自衡对象控制方法。本发明方法首先基于非自衡对象的阶跃响应数据建立被控对象的模型,挖掘出基本的对象特性;然后依据预测函数控制的特性去整定相应P控制器的参数;最后对被控对象实施P控制。本发明将预测函数控制的性能赋给了P控制,有效地提高了传统控制方法的性能,同时也促进了先进控制方法的应用。
申请公布号 CN103336437B 申请公布日期 2015.11.25
申请号 CN201310310130.0 申请日期 2013.07.19
申请人 杭州电子科技大学 发明人 张日东;吴胜
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 杭州君度专利代理事务所(特殊普通合伙) 33240 代理人 杜军
主权项 基于预测函数控制优化的非自衡对象控制方法,其特征在于该方法的具体步骤是:步骤(1).通过被控对象的实时阶跃响应数据建立被控对象的模型,具体是:Ⅰ.给被控对象一个阶跃输入信号,记录被控对象的阶跃响应曲线;Ⅱ.将对应的阶跃响应曲线滤波处理成一条光滑曲线,然后将数值发生变化的曲线段拟合成一条直线,计算出模型的增益参数K<sub>m</sub>:K<sub>m</sub>=q*l其中,q为过程模型的控制量的阶跃变化幅度,K<sub>m</sub>是建立的被控对象模型的增益系数,l是拟合的直线的斜率;Ⅲ.记录Ⅱ中滤波处理后的光滑曲线上每个采样时刻对应的阶跃响应数据,第一个采样时刻为T<sub>s1</sub>,相邻两个采样时刻间隔的时间为T<sub>s1</sub>,采样时刻顺序为T<sub>s1</sub>、2T<sub>s1</sub>、3T<sub>s1</sub>……;在记录的阶跃响应数据中找出数据开始上升的起始点a<sub>I</sub>,之前的数据分别记做a<sub>1</sub>,a<sub>2</sub>,…a<sub>I‑1</sub>,模型的滞后时间参数τ为τ=(I‑1)T<sub>s1</sub>,最后得到的被控对象的传递函数模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>K</mi><mi>m</mi></msub><mi>s</mi></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mi>&tau;</mi><mi>s</mi></mrow></msup></mrow>]]></math><img file="FDA0000750054600000011.GIF" wi="317" he="130" /></maths>其中,G(s)为被控对象的传递函数,s为拉普拉斯变换算子;步骤(2).设计被控对象的P控制器,具体是:a.对得到的传递函数在采样时间T<sub>s</sub>下加一个零阶保持器离散化,得到离散模型为:y<sub>m</sub>(k)=y<sub>m</sub>(k‑1)+K<sub>m</sub>T<sub>s</sub>u(k‑1‑L)y<sub>m</sub>(k)为k时刻的模型预测输出,u(k‑1‑L)为k‑1‑L时刻的控制输入,L为离散传递函数模型的时滞,L=τ/T<sub>s</sub>;b.计算被控对象去掉纯滞后以后在预测函数控制下的第P步预测输出,形式如下:y<sub>mav</sub>(k)=y<sub>mav</sub>(k‑1)+K<sub>m</sub>T<sub>s</sub>u(k‑1)y<sub>mav</sub>(k+P)=y<sub>mav</sub>(k)+K<sub>m</sub>PT<sub>s</sub>u(k)其中,P为预测步长,y<sub>mav</sub>(k+P)为k时刻去掉纯滞后的被控对象在预测函数控制下的第P步预测输出,y<sub>mav</sub>(k)为k时刻去掉纯滞后的模型输出;c.校正k时刻的实际输出得到包含未来预测信息的新的过程输出值,形式如下:y<sub>pav</sub>(k)=y<sub>p</sub>(k)+y<sub>mav</sub>(k)‑y<sub>mav</sub>(k‑L)其中,y<sub>pav</sub>(k)为校正得到的k时刻包含未来预测信息的新的过程输出值,y<sub>p</sub>(k)为k时刻的实际输出值;d.选取预测函数控制方法的参考轨迹y<sub>r</sub>(k+P)以及目标函数J,形式如下:y<sub>r</sub>(k+P)=β<sup>P</sup>y<sub>p</sub>(k)+(1‑β<sup>P</sup>)c(k)J=min(y<sub>r</sub>(k+P)‑y<sub>mav</sub>(k+P)‑e(k))<sup>2</sup>e(k)=y<sub>pav</sub>(k)‑y<sub>mav</sub>(k)+(P+L)(y<sub>pav</sub>(k)‑y<sub>mav</sub>(k)‑y<sub>pav</sub>(k‑1)+y<sub>mav</sub>(k‑1))其中,β为参考轨迹柔化系数,c(k)为k时刻的设定值,e(k)为k时刻校正的误差值,这里采用新的误差校正方式;e.依据步骤d中的目标函数求解P控制器中的参数,这里将控制量u(k)进行变换:u(k)=K<sub>p</sub>e<sub>1</sub>(k)e<sub>1</sub>(k)=βy<sub>p</sub>(k‑1)+(1‑β)c(k‑1)‑y<sub>p</sub>(k)其中,K<sub>p</sub>为位置式P控制器的参数,e<sub>1</sub>(k)为k时刻参考轨迹值与实际输出值之间的误差;结合上述式子,可以求得:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>K</mi><mi>p</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>y</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>P</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mrow><mi>m</mi><mi>a</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>K</mi><mi>m</mi></msub><msub><mi>PT</mi><mi>s</mi></msub><msub><mi>e</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000750054600000021.GIF" wi="696" he="148" /></maths>f.得到P控制器的参数K<sub>p</sub>以后构成控制量u(k)=K<sub>p</sub>e<sub>1</sub>(k)作用于被控对象;g.在下一时刻,依照b到f中的步骤继续求解新的K<sub>p</sub>,依次循环。
地址 310018 浙江省杭州市下沙高教园区2号大街