发明名称 复杂滞后过程的先进控制方法及其系统
摘要 一种复杂滞后过程的先进控制方法及实现该方法的系统,至少包括:预测该复杂滞后过程输出的未来值;对该过程当前的输出进行预测;对该过程进行在线反馈校正;对该过程进行最优控制;依据约束限对计算得到的最优控制率判断,得到满足约束条件的控制增量,并将该控制增量施加到被控过程中去。本发明依据被控过程具体情况和测点情况,选取相关的过程变量为状态变量,使用实测的状态变量计算最优控制率;并通过选取可测干扰变量设计模型预测前馈来计算最优控制率,实现了实测状态变量反馈设计,增大了控制系统对被控过程工艺变化、操作控制变化的适应范围,提高了控制系统的鲁棒稳定性及抗干扰能力。
申请公布号 CN1285016C 申请公布日期 2006.11.15
申请号 CN03123876.9 申请日期 2003.06.09
申请人 石油大学(北京) 发明人 胡品慧;闫峰;袁璞
分类号 G05B13/04(2006.01);G05B11/01(2006.01) 主分类号 G05B13/04(2006.01)
代理机构 北京同立钧成知识产权代理有限公司 代理人 刘芳;刘国伟
主权项 1、一种复杂滞后过程的先进控制方法,其特征在于:该方法至少包括如下的步骤:步骤1:选择第h个输出y<sub>h</sub>的预测时域为p<sub>h</sub>,预测y<sub>h</sub>未来第p<sub>h</sub>步的值,即k+p<sub>h</sub>时刻的未来输出值,由下式确定:<maths num="001"><![CDATA[ <math><mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mi>h</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><msub><mi>p</mi><mi>h</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>c</mi><mi>h</mi></msub><msup><mi>A</mi><msub><mi>p</mi><mi>h</mi></msub></msup><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>l</mi><mi>a</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>p</mi><mi>h</mi></msub></munderover><msub><mi>c</mi><mi>h</mi></msub><msup><mi>A</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>A</mi><mi>i</mi></msub><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><msub><mi>p</mi><mi>h</mi></msub><mo>-</mo><msub><mi>&tau;</mi><mi>i</mi></msub><mo>-</mo><mi>j</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="002"><![CDATA[ <math><mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>l</mi><mi>b</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>p</mi><mi>h</mi></msub></munderover><msub><mi>c</mi><mi>h</mi></msub><msup><mi>A</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>B</mi><mi>i</mi></msub><mi>U</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><msub><mi>p</mi><mi>h</mi></msub><mo>-</mo><msub><mi>d</mi><mi>i</mi></msub><mo>-</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>p</mi><mi>h</mi></msub></munderover><msub><mi>c</mi><mi>h</mi></msub><msup><mi>A</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msup><mi>FV</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><msub><mi>p</mi><mi>h</mi></msub><mo>-</mo><mi>j</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="C031238760002C3.GIF" wi="215" he="50" />为对y<sub>h</sub>(k+p<sub>h</sub>)的预测值;y<sub>h</sub>(k+p<sub>h</sub>)为第h个过程输出y<sub>h</sub>的预测时域为p<sub>h</sub>的k+p<sub>h</sub>时刻的未来输出值;c<sub>h</sub>为输出矩阵C的第h行向量,且C=[c<sub>1</sub><sup>T</sup>,…,c<sub>h</sub><sup>T</sup>,…,c<sub>r</sub><sup>T</sup>]<sup>T</sup>;τ<sub>i</sub>为状态变量纯滞后时间,d<sub>i</sub>为控制纯滞后时间,j代表未来j时刻,V为干扰变量,V(k+p<sub>h</sub>-j)为干扰变量未来p<sub>h</sub>-j时刻的值,U为操纵变量,U(k+p<sub>h</sub>-d<sub>i</sub>-j)为操纵变量未来p<sub>h</sub>-d<sub>i</sub>-j时刻的值,X为状态变量,X(k)为状态变量当前时刻的值,X(k+p<sub>h</sub>-τ<sub>i</sub>-j)为状态变量p<sub>h</sub>-τ<sub>i</sub>-j时刻的值;A∈R<sup>n×n</sup>,A<sub>i</sub>∈R<sup>n×n</sup>,i=1,2,…l<sub>a</sub>,B<sub>i</sub>∈R<sup>n×r</sup>,i=1,2,…l<sub>b</sub>,C∈R<sup>r×n</sup>,F∈R<sup>n×q</sup>,均为已知的常数矩阵,用于计算最优控制率;然后处理未来状态及未来干扰的预测值;步骤2:对该复杂滞后过程当前时刻的输出值进行预测;步骤3:对该复杂滞后过程进行在线反馈校正;步骤4:对该复杂滞后过程进行最优控制;步骤5:依据规则库的判断逻辑规则关系,对计算得到的各种变量和最优控制率进行约束限判断,得到满足约束条件的控制增量,并将该控制增量施加到被控过程中去。
地址 102249北京市昌平区府学路20号