[发明专利]一种基于函数型权RBF-ARX模型的双回路水箱液位预测控制方法在审

专利信息
申请号: 201610139368.5 申请日: 2016-03-11
公开(公告)号: CN105676645A 公开(公告)日: 2016-06-15
发明(设计)人: 彭辉;周锋;覃业梅 申请(专利权)人: 中南大学
主分类号: G05B13/04 分类号: G05B13/04;G05D9/12
代理公司: 长沙正奇专利事务所有限责任公司 43113 代理人: 马强;王娟
地址: 410083 湖南*** 国省代码: 湖南;43
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明公开了一种基于函数型权RBF-ARX模型的双回路水箱液位预测控制方法。该方法采用一种数据驱动的系统辨识技术,设计了双回路水箱液位系统的一种基于函数型权RBF-ARX模型的建模方法。本发明采用的一种基于函数型权的RBF-ARX模型,可有效减少RBF网络隐层结点的数目,与其它一般化非线性ARX模型相比,该模型有更好的预测精度,是一种适用性广的、可有效描述双回路水箱液位系统全局非线性动态特性的建模方法。本发明基于双回路水箱系统函数型权RBF-ARX模型的全局非线性特性,采用带约束的二次型性能指标,设计了蕴含模型全部非线性信息的预测控制算法,可进一步提高双回路水箱液位控制系统的动静态性能指标,具有较高的实用价值和较好的应用前景。
搜索关键词: 一种 基于 函数 rbf arx 模型 回路 水箱 预测 控制 方法
【主权项】:
一种基于函数型权RBF‑ARX模型的双回路水箱液位预测控制方法,其特征在于,包括以下步骤:1)对双回路水箱系统离线建立如下结构函数型权RBF‑ARX模型:<mrow><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&phi;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>w</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>=</mo><mn>1</mn></mrow><mn>22</mn></munderover><msub><mi>&phi;</mi><mrow><mi>f</mi><mo>,</mo><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>w</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>k</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><mn>1</mn></mrow><mn>20</mn></munderover><msub><mi>&phi;</mi><mrow><mi>e</mi><mo>,</mo><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>w</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mi>e</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>k</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>&xi;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>;</mo></mrow><mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&phi;</mi><mi>0</mi></msub><mrow><mo>(</mo><mi>w</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mo>=</mo><msubsup><mi>c</mi><mn>0</mn><mn>0</mn></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msubsup><mi>v</mi><mi>m</mi><mn>0</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msubsup><mi>&lambda;</mi><mi>m</mi><mi>f</mi></msubsup><mo>|</mo><mo>|</mo><mi>w</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>-</mo><msubsup><mi>z</mi><mi>m</mi><mi>f</mi></msubsup><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&phi;</mi><mrow><mi>f</mi><mo>,</mo><msub><mi>k</mi><mn>1</mn></msub></mrow></msub><mrow><mo>(</mo><mi>w</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mo>=</mo><msubsup><mi>c</mi><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><mn>0</mn></mrow><mn>0</mn></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>1</mn></msub></mrow><mn>0</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msubsup><mi>&lambda;</mi><mi>m</mi><mi>f</mi></msubsup><mo>|</mo><mo>|</mo><mi>w</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>-</mo><msubsup><mi>z</mi><mi>m</mi><mi>f</mi></msubsup><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&phi;</mi><mrow><mi>e</mi><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub></mrow></msub><mrow><mo>(</mo><mi>w</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mo>=</mo><msubsup><mi>c</mi><mrow><msub><mi>k</mi><mn>2</mn></msub><mo>,</mo><mn>0</mn></mrow><mi>e</mi></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub></mrow><mi>e</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msubsup><mi>&lambda;</mi><mi>m</mi><mi>e</mi></msubsup><mo>|</mo><mo>|</mo><mi>w</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>-</mo><msubsup><mi>z</mi><mi>m</mi><mi>e</mi></msubsup><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>v</mi><mi>m</mi><mn>0</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><mn>0</mn></mrow><mn>0</mn></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><mn>1</mn></mrow><mn>0</mn></msubsup><msubsup><mi>w</mi><mn>1</mn><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><mn>2</mn></mrow><mn>0</mn></msubsup><msubsup><mi>w</mi><mn>2</mn><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>1</mn></msub></mrow><mi>f</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><mn>0</mn></mrow><mi>f</mi></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><mn>1</mn></mrow><mi>f</mi></msubsup><msubsup><mi>w</mi><mn>1</mn><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><mn>2</mn></mrow><mi>f</mi></msubsup><msubsup><mi>w</mi><mn>2</mn><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub></mrow><mi>e</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub><mo>,</mo><mn>0</mn></mrow><mi>e</mi></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub><mo>,</mo><mn>1</mn></mrow><mi>e</mi></msubsup><msubsup><mi>w</mi><mn>1</mn><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>m</mi><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub><mo>,</mo><mn>2</mn></mrow><mi>e</mi></msubsup><msubsup><mi>w</mi><mn>2</mn><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>w</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>w</mi><mn>1</mn><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mtd><mtd><msubsup><mi>w</mi><mn>2</mn><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>.</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>其中:f(t)为t时刻的水箱液位输出向量;e(t)为t时刻的水箱电磁调节阀门开度输入向量;ξ(t)为高斯白噪声;状态向量为2范数;为RBF神经网络的中心向量和缩放因子;为RBF神经网络的函数型权系数,是水箱状态量w(t‑1)的线性函数;是RBF神经网络相应的权重系数;为常数系数;非线性参数和线性参数均通过SNPOM优化方法离线优化辨识获得;2)基于水箱系统函数型权RBF‑ARX模型的全局非线性特性设计预测控制器结构如下:<mfenced open = "{" close = ""><mtable><mtr><mtd><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mover><mi>V</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></munder><mi>J</mi><mo>=</mo><mo>|</mo><mo>|</mo><mover><mi>F</mi><mo>^</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><msub><mi>F</mi><mi>r</mi></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>|</mo><msubsup><mo>|</mo><mi>Q</mi><mn>2</mn></msubsup><mo>+</mo><mo>|</mo><mo>|</mo><mover><mi>E</mi><mo>^</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>|</mo><msubsup><mo>|</mo><msub><mi>R</mi><mn>1</mn></msub><mn>2</mn></msubsup><mo>+</mo><mo>|</mo><mo>|</mo><mi>&Delta;</mi><mover><mi>E</mi><mo>^</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>|</mo><msubsup><mo>|</mo><msub><mi>R</mi><mn>2</mn></msub><mn>2</mn></msubsup></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><msub><mi>E</mi><mi>min</mi></msub><mo>&le;</mo><mover><mi>E</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>E</mi><mi>max</mi></msub><mo>,</mo><msub><mi>&Delta;E</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>&le;</mo><mi>&Delta;</mi><mover><mi>E</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>&Delta;E</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></mtd></mtr></mtable></mtd></mtr></mtable></mfenced>其中,Fr(t)分别表示t时刻的预测输出序列和期望输出序列,为t时刻模型预测的t+g时刻的输出液位,f(t+g|t)为t时刻输出液位的t+g时刻的期望值;分别表示t时刻的预测控制输入序列和输入增量序列,Δe(t)=e(t)‑e(t‑1)为t时刻的输入增量;e(t+q)=e(t+5)为t+q时刻的输入信号,30≥q≥6;R1=diag{R1,…,R1}6、R2=diag{R2,…,R2}6为控制加权矩阵,R1=[0.0001 0.0、00R2=[0.28 0.32];Q=diag{Q,…,Q}30为误差加权矩阵,Q=[1 1];Emin、Emax为控制输入量约束序列,ΔEmin、ΔEmax为控制输入增量约束序列,其中,Emin=[0,0;…;0,0;]6×2、Emax=[100,100;…;100,100;]6×2、ΔEmin=[‑10,‑10;…;‑10,‑10]6×2、ΔEmax=[10,10;…;10,10]6×2;基于水箱系统函数型权RBF‑ARX模型的全局非线性特性设计的预测控制算法,通过实时优化得到电磁调节阀的开发控制信号量e(t),最终达到精确控制水箱中液位的目的。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中南大学,未经中南大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201610139368.5/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top