发明名称 一种基于预测补偿的网络化控制系统跟踪控制方法
摘要 本发明公开了一种基于预测补偿的网络化控制系统跟踪控制方法,网络预测控制器由预测产生器和时延补偿器两部分组成,k时刻网络预测产生器输出的控制增量序列为:<img file="DDA0001060912400000011.GIF" wi="1030" he="75" />时延补偿器t时刻的控制增量为:<img file="DDA0001060912400000012.GIF" wi="1350" he="86" />本发明以引入很小的预测误差为代价缩短预测步长,提高了预测控制器对网络时延的补偿效果,相比传统的网络预测控制更易实现。
申请公布号 CN106209474A 申请公布日期 2016.12.07
申请号 CN201610603162.3 申请日期 2016.07.27
申请人 江南大学 发明人 潘丰;刘婷
分类号 H04L12/24(2006.01)I;H04L7/00(2006.01)I 主分类号 H04L12/24(2006.01)I
代理机构 代理人
主权项 一种基于预测补偿的网络化控制系统跟踪控制方法,其特征在于,网络预测控制器由预测产生器和时延补偿器两部分组成,具体包括以下步骤:(1)一类离散时间线性时不变多输入多输出受控系统,其增广系统状态空间描述为:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mover><mi>A</mi><mo>&OverBar;</mo></mover><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mover><mi>B</mi><mo>&OverBar;</mo></mover><mi>&Delta;</mi><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>C</mi><mo>&OverBar;</mo></mover><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001060912370000011.GIF" wi="406" he="135" /></maths>其中,<img file="FDA0001060912370000012.GIF" wi="358" he="62" />Δu(t)=u(t)‑u(t‑1);x(t)∈R<sup>n</sup>、u(t)∈R<sup>m</sup>和y(t)∈R<sup>l</sup>分别为被控对象的状态向量、控制输入和控制输出,A∈R<sup>n×n</sup>,B∈R<sup>n×m</sup>,C∈R<sup>l×n</sup>为系数矩阵;系统满足(A,B)完全可控,(A,C)完全可观;系统输出y(t)跟踪外部参考输入y<sub>0</sub>(t)实现渐近跟踪,即跟踪误差e(t)满足:<img file="FDA0001060912370000013.GIF" wi="1206" he="110" />状态不可直接测,建立状态观测器:<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>|</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>A</mi><mo>&OverBar;</mo></mover><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>|</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mover><mi>B</mi><mo>&OverBar;</mo></mover><mi>&Delta;</mi><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mover><mi>L</mi><mo>&OverBar;</mo></mover><mo>&lsqb;</mo><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>|</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>|</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mover><mi>C</mi><mo>&OverBar;</mo></mover><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>|</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001060912370000014.GIF" wi="788" he="119" /></maths>其中,x(t+1|t)为当前时刻对下一时刻的状态观测值,x(t|t‑1)和y(t|t‑1)分别为前一时刻对当前时刻的状态观测值和输出观测值,L∈R<sup>n*l</sup>为观测器增益;增广观测器增益<img file="FDA0001060912370000015.GIF" wi="424" he="118" />(2)k时刻网络预测产生器输出的控制增量序列为:<maths num="0003"><math><![CDATA[<mrow><mi>&Delta;</mi><msup><mover><mi>U</mi><mo>&OverBar;</mo></mover><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><msub><mi>&tau;</mi><mrow><mi>s</mi><mi>c</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>K</mi><msub><mi>&tau;</mi><mrow><mi>s</mi><mi>c</mi><mi>k</mi></mrow></msub></msub><mo>&CenterDot;</mo><mo>&lsqb;</mo><msub><mi>&Gamma;</mi><msub><mi>&tau;</mi><mrow><mi>s</mi><mi>c</mi><mo>,</mo><mi>k</mi></mrow></msub></msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>-</mo><msub><mi>&tau;</mi><mrow><mi>s</mi><mi>c</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>+</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><msub><mi>&tau;</mi><mrow><mi>s</mi><mi>c</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>Y</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001060912370000016.GIF" wi="989" he="76" /></maths>其中,<img file="FDA0001060912370000017.GIF" wi="646" he="71" /><img file="FDA0001060912370000018.GIF" wi="85" he="53" />为1×(N<sub>u</sub>+τ<sub>sc,k</sub>‑1)的块矩阵,从第τ<sub>sc,k</sub>+1块开始为单位矩阵I∈R<sup>m×m</sup>;k时刻数据包在反馈通道中的实时时延为τ<sub>sc,k</sub>,N<sub>p</sub>为预测区间长度,N<sub>u</sub>‑1为优化区间长度;<img file="FDA0001060912370000019.GIF" wi="70" he="52" />为N<sub>p</sub>×1的块矩阵,<img file="FDA00010609123700000110.GIF" wi="542" he="71" /><img file="FDA00010609123700000111.GIF" wi="83" he="55" />为下三角矩阵,满足j‑i≤τ<sub>sc,k</sub>‑1且<img file="FDA00010609123700000112.GIF" wi="401" he="79" />Y<sub>0</sub>(k)为系统跟随目标值序列:Y<sub>0</sub>(k)=[y<sub>0</sub>(k+1) … y<sub>0</sub>(k+N<sub>p</sub>)]<sup>T</sup>;R、Q为权重系数矩阵;系统输出和控制输入分别为y(k‑τ<sub>sc,k</sub>)和u(k‑τ<sub>sc,k</sub>),<img file="FDA00010609123700000113.GIF" wi="334" he="62" />为k时刻状态观测器最新更新的状态观测值:<img file="FDA00010609123700000114.GIF" wi="1330" he="63" />(3)时延补偿器t时刻的控制增量为:<maths num="0004"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mrow><msub><msup><mi>&tau;</mi><mo>*</mo></msup><mrow><mi>c</mi><mi>a</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow></msub><mo>&CenterDot;</mo><mi>&Delta;</mi><msup><mover><mi>U</mi><mo>&OverBar;</mo></mover><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><msup><mi>&tau;</mi><mo>*</mo></msup><mrow><mi>c</mi><mi>a</mi><mo>.</mo><mi>t</mi></mrow></msub><mo>|</mo><mi>t</mi><mo>-</mo><msub><msup><mi>&tau;</mi><mo>*</mo></msup><mi>t</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>S</mi><mrow><msub><msup><mi>&tau;</mi><mo>*</mo></msup><mrow><mi>c</mi><mi>a</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow></msub><msub><mi>K</mi><mrow><msub><msup><mi>&tau;</mi><mo>*</mo></msup><mrow><mi>s</mi><mi>c</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow></msub><mo>&lsqb;</mo><msub><mi>&Gamma;</mi><mrow><msub><msup><mi>&tau;</mi><mo>*</mo></msup><mrow><mi>s</mi><mi>c</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow></msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><msup><mi>&tau;</mi><mo>*</mo></msup><mi>t</mi></msub><mo>+</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><msub><msup><mi>&tau;</mi><mo>*</mo></msup><mi>t</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>Y</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><msup><mi>&tau;</mi><mo>*</mo></msup><mrow><mi>c</mi><mi>a</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA00010609123700000115.GIF" wi="1324" he="86" /></maths>其中,<img file="FDA00010609123700000116.GIF" wi="78" he="62" />为1×(N<sub>u</sub>‑1)的块矩阵,第<img file="FDA00010609123700000117.GIF" wi="61" he="55" />块为单位矩阵I∈R<sup>m×m</sup>,<img file="FDA00010609123700000118.GIF" wi="301" he="62" />为<sub>t</sub>时刻补偿器中最新控制增量序列,τ<sup>*</sup><sub>sc,t</sub>和τ<sup>*</sup><sub>ca,t</sub>分别为时延补偿器选取的最新优化控制增量序列对应的时延;网络化跟踪控制闭环系统为:x(t+1)=Ax(t)+B[u(t‑1)+Δu(t)]。
地址 214122 江苏省无锡市滨湖区蠡湖大道1800号