发明名称 一种基于RBF-ARX模型稳定参数估计的非线性系统建模方法
摘要 本发明公开了一种基于RBF‑ARX模型稳定参数估计的非线性系统建模方法,针对实际工业过程中普遍存在的非线性动态和物理模型难以获取的情况,采用RBF‑ARX模型描述对象的动态特性,并针对该类模型参数较多及潜在的病态可能导致的模型参数过大或不稳定的问题,将待估计的模型参数分类成线性参数和非线性参数分开优化,设计了一种融合了正则化技术和非线性最小二乘法算法的可保证该类模型参数稳定的离线参数估计方法。与现有技术相比,本发明可大幅提高RBF‑ARX模型的长期预测精度和鲁棒性,并适用于所有能将参数进行分类成线性参数部分和非线性参数部分的非线性模型参数优化问题。对于基于计算机数值仿真分析的工程设计和优化问题具有很高的实用价值。
申请公布号 CN106021829A 申请公布日期 2016.10.12
申请号 CN201610569906.4 申请日期 2016.07.19
申请人 中南大学 发明人 彭辉;曾小勇;周峰;田晓盈;孙亚鹏;覃业梅
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 长沙正奇专利事务所有限责任公司 43113 代理人 马强;王娟
主权项 一种基于RBF‑ARX模型稳定参数估计的非线性系统建模方法,其特征在于,包括以下步骤:1)构建RBF‑ARX模型:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&phi;</mi><mn>0</mn></msub><mrow><mo>(</mo><mrow><mi>w</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>&phi;</mi><mrow><mi>y</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>w</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mi>y</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mi>i</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>q</mi></munderover><msub><mi>&phi;</mi><mrow><mi>u</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>w</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mi>u</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mi>j</mi></mrow><mo>)</mo></mrow><mo>+</mo><mi>e</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&phi;</mi><mn>0</mn></msub><mrow><mo>(</mo><mrow><mi>w</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>=</mo><msubsup><mi>c</mi><mn>0</mn><mn>0</mn></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>c</mi><mi>k</mi><mn>0</mn></msubsup><mi>exp</mi><mrow><mo>{</mo><mrow><mo>-</mo><msubsup><mi>&lambda;</mi><mi>k</mi><mi>y</mi></msubsup><mo>|</mo><mo>|</mo><mi>w</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mo>-</mo><msubsup><mi>z</mi><mi>k</mi><mi>y</mi></msubsup><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup></mrow><mo>}</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&phi;</mi><mrow><mi>y</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>w</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>=</mo><msubsup><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mn>0</mn></mrow><mi>y</mi></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi></mrow><mi>y</mi></msubsup><mi>exp</mi><mrow><mo>{</mo><mrow><mo>-</mo><msubsup><mi>&lambda;</mi><mi>k</mi><mi>y</mi></msubsup><mo>|</mo><mo>|</mo><mi>w</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mo>-</mo><msubsup><mi>z</mi><mi>k</mi><mi>y</mi></msubsup><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup></mrow><mo>}</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&phi;</mi><mrow><mi>u</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>w</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>=</mo><msubsup><mi>c</mi><mrow><mi>j</mi><mo>,</mo><mn>0</mn></mrow><mi>u</mi></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>c</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow><mi>u</mi></msubsup><mi>exp</mi><mrow><mo>{</mo><mrow><mo>-</mo><msubsup><mi>&lambda;</mi><mi>k</mi><mi>u</mi></msubsup><mo>|</mo><mo>|</mo><mi>w</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mo>-</mo><msubsup><mi>z</mi><mi>k</mi><mi>u</mi></msubsup><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup></mrow><mo>}</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>z</mi><mi>k</mi><mi>v</mi></msubsup><mo>=</mo><msup><mrow><mo>(</mo><mrow><msubsup><mi>z</mi><mrow><mi>k</mi><mo>,</mo><mn>1</mn></mrow><mi>v</mi></msubsup><mo>,</mo><msubsup><mi>z</mi><mrow><mi>k</mi><mo>,</mo><mn>2</mn></mrow><mi>v</mi></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>z</mi><mrow><mi>k</mi><mo>,</mo><mi>d</mi></mrow><mi>v</mi></msubsup></mrow><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><mi>v</mi><mo>=</mo><mi>y</mi><mo>,</mo><mi>u</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001053363030000011.GIF" wi="1013" he="662" /></maths>其中:<img file="FDA0001053363030000012.GIF" wi="157" he="54" />为t时刻的非线性系统的输出;<img file="FDA0001053363030000013.GIF" wi="158" he="53" />为t时刻的非线性系统的输入;<img file="FDA0001053363030000014.GIF" wi="153" he="55" />为高斯白噪声;p和q分别为非线性系统模型的输出和输入阶次;φ<sub>0</sub>,φ<sub>y</sub>,<sub>i</sub>和φ<sub>u,j</sub>为不同的RBF神经网络;w(t‑1)为RBF神经网络的输入向量,用以引导RBF‑ARX模型描述非线性系统的动态特性;<img file="FDA0001053363030000015.GIF" wi="69" he="63" /><img file="FDA0001053363030000016.GIF" wi="507" he="62" />为RBF神经网络的线性权重;<img file="FDA0001053363030000017.GIF" wi="50" he="54" />和<img file="FDA0001053363030000018.GIF" wi="51" he="55" />为RBF神经网络的缩放因子;<img file="FDA0001053363030000019.GIF" wi="45" he="55" />和<img file="FDA00010533630300000110.GIF" wi="46" he="60" />为RBF神经网络的中心向量;<img file="FDA00010533630300000111.GIF" wi="69" he="70" />为向量的2范数;m为RBF神经网络中间层的结点个数;<img file="FDA00010533630300000112.GIF" wi="270" he="62" />为RBF神经网络的中心向量的具体数值,其中v表示非线性系统的输入或输出,d为RBF神经网络输入向量w(t‑1)的维数;j=1,2,…,q;i=1,2,…,p;k=1,2,…,m。2)将所述RBF‑ARX模型表述为如下非线性最小二乘形式:<img file="FDA00010533630300000113.GIF" wi="540" he="60" />其中,<img file="FDA00010533630300000114.GIF" wi="1054" he="261" /><img file="FDA0001053363030000021.GIF" wi="1101" he="255" /><img file="FDA0001053363030000022.GIF" wi="964" he="534" />3)构造如下带线性参数惩罚的优化目标函数:<maths num="0002"><math><![CDATA[<mrow><mo>(</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>N</mi></msub><mo>,</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>L</mi></msub><mo>)</mo><mo>=</mo><mi>arg</mi><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><msub><mi>&theta;</mi><mi>N</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>L</mi></msub></mrow></munder><mo>|</mo><mo>|</mo><mi>Y</mi><mo>-</mo><mi>&Phi;</mi><mo>(</mo><msub><mi>&theta;</mi><mi>N</mi></msub><mo>)</mo><msub><mi>&theta;</mi><mi>L</mi></msub><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&gamma;&theta;</mi><mi>L</mi><mi>T</mi></msubsup><msup><mi>P</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&theta;</mi><mi>L</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0001053363030000023.GIF" wi="797" he="87" /></maths>其中Y为非线性系统输出y(t)的向量;Φ(θ<sub>N</sub>)为<img file="FDA0001053363030000024.GIF" wi="266" he="56" />的向量;γ和P为需要调整的正则化参数。
地址 410083 湖南省长沙市岳麓区麓山南路932号