发明名称 基于模型预测控制的集结预测控制系统及其控制方法
摘要 一种工业过程控制技术领域的基于模型预测控制的集结预测控制系统及其控制方法,包括:输入映射模块、传统预测控制模块、传感器和逆映射模块,其中:输入映射模块与传统预测控制模块相连接并将预测控制的输入控制序列进行线性映射,传统预测控制模块与控制对象相连以接接收系统对象状态值,传感器与控制对象相连且传感器的输出端与传统预测控制模块的输入端相连以传输检测信号,传统预测控制模块的输出端连接至逆映射模块,逆映射模块的输出端与控制对象相连接以输出控制策略。本发明的方法对集结预测控制器的控制性能进行优化,在降低预测控制器的在线计算量,提高其实时性的同时,保证了控制器的控制性能,具有很强的实用性和应用价值。
申请公布号 CN101930213B 申请公布日期 2013.09.11
申请号 CN200910197512.0 申请日期 2009.10.22
申请人 上海交通大学 发明人 李德伟;席裕庚;姜舒
分类号 G05B13/02(2006.01)I 主分类号 G05B13/02(2006.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 张泽纯
主权项 1.一种集结预测控制系统,其特征在于,包括:输入映射模块、传统预测控制模块、传感器和逆映射模块,其中:输入映射模块与传统预测控制模块相连接并将预测控制的输入控制序列进行线性映射,传感器与控制对象相连且传感器的输出端与传统预测控制模块的输入端相连以传输检测信号,传统预测控制模块的输出端连接至逆映射模块,逆映射模块的输出端与控制对象相连接以输出控制策略;所述的集结预测控制系统的控制方法包括以下步骤:第一步、采用状态方程<img file="FDA00002898145800011.GIF" wi="394" he="98" />描述化工过程连续搅拌釜式反应器系统,其中:状态向量x(t)为反应器内液体的浓度和温度,控制量u(t)为冷却剂的流量,<img file="FDA00002898145800012.GIF" wi="56" he="82" />和<img file="FDA00002898145800013.GIF" wi="52" he="76" />分别为系数矩阵;对状态方程进行离散化,得到离散状态方程x(k+1)=Ax(k)+Bu(k);第二步、工控机根据采样周期,在k时刻发布采样命令,对化工过程连续搅拌釜式反应器系统的状态向量进行检测,检测的信号通过模拟量输入通道经A/D转换后,输送到工控机中;第三步、工控机第二步检测得到的系统状态,经过输入映射模块映射处理U(k)=H<sub>1</sub>V(k),将优化变量变为V(k),其中:U(k)为k时刻的预测控制输入序列,H<sub>1</sub>为映射系数矩阵;第四步,逆映射模块接收到传统预测控制模块的输出后,进行逆映射计算得到优化控制序列,即当前时刻预测控制器的预测控制序列,并将控制序列的第一个控制变量输出到控制对象;第五步、工控机将预测控制输入序列的第一个控制变量输出到控制对象,优化控制完成,在下一个采样周期,重复从第二步开始执行新的优化控制周期;所述的优化变量满足以下条件:<maths num="0001"><![CDATA[<math><mrow><munder><mi>min</mi><mrow><mi>V</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></munder><mi>J</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><mi>N</mi></munderover><msup><mi>x</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>|</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>Q</mi><mi>i</mi></msub><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>|</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>u</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>|</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>R</mi><mi>i</mi></msub><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>|</mo><mi>k</mi><mo>)</mo></mrow></mrow></math>]]></maths>s.t. x(k+i+1|k)=Ax(k+i|k)+Bu(k+i|k)U(k)=H<sub>1</sub>V(k),U(k)∈Ω<sub>u</sub>,x(k+i|k)∈Ω<sub>x</sub>其中:U(k)=[u<sup>T</sup>(k),...,u<sup>T</sup>(k+N-1|k)]<sup>T</sup>为预测控制在k时刻预测控制输入序列,Q<sub>i</sub>和R<sub>i</sub>为系统状态和输入加权矩阵,N为控制时域长度;所述的优化控制序列满足逆映射计算公式:U(k)=H<sub>2</sub>V(k),其中,H<sub>2</sub>映射系数矩阵与输入映射计算公式中的映射系数矩阵H<sub>1</sub>相同,映射系数矩阵H<sub>1</sub>或H<sub>2</sub>通过以下步骤获得:1)根据控制对象差分方程和传统预测控制模块设计,通过对象状态x(k)和预测控制输入序列U(k)表述时域长度内系统状态序列X(k),具体表达式如下:X(k)=Sx(k)+GU(k),其中:G和S分别通过系数矩阵A和B进行表述:<maths num="0002"><![CDATA[<math><mrow><mi>S</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>A</mi></mtd></mtr><mtr><mtd><msup><mi>A</mi><mn>2</mn></msup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msup><mi>A</mi><mi>N</mi></msup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><img file="FDA00002898145800022.GIF" wi="587" he="314" />2)选取[A<sup>N-1</sup>B,A<sup>N-2</sup>B,…,B]中线性无关的列向量做为基向量,求解方程组[A<sup>N-1</sup>B,A<sup>N-2</sup>B,…,B]x=0的解空间的基U<sub>0</sub>,得到矩阵G<sub>1</sub>=GU<sub>0</sub>,然后计算中间矩阵<img file="FDA00002898145800023.GIF" wi="528" he="92" />将中间矩阵W按列分成W<sub>1</sub>∈R<sup>(Nm-n)×m</sup>和W<sub>2</sub>∈R<sup>(Nm-n)×(Nm-m)</sup>两个部分,并计算出列满秩的基本映射系数矩阵H<sub>10</sub>∈R<sup>(Nm-n)×m</sup>,满足<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>H</mi><mn>10</mn><mi>T</mi></msubsup><msub><mi>W</mi><mn>2</mn></msub><mo>=</mo><mn>0</mn><mo>;</mo></mrow></math>]]></maths>3)在H<sub>10</sub>矩阵中对应[A<sup>N-1</sup>B,A<sup>N-2</sup>B,…,B]中基向量所处的列位置插入单位向量的行列,其中1元素所在的行号对应于该向量所在的列号,得到基本映射系数矩阵H<sub>0</sub>,具体如下:<img file="FDA00002898145800031.GIF" wi="1338" he="779" />4)基本映射系数矩阵的自由时域长度优化,根据实际应用的操作范围,选取时域长度,将H<sub>0</sub>分成<maths num="0004"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>H</mi><mn>01</mn></msub></mtd></mtr><mtr><mtd><msub><mi>H</mi><mn>02</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>满足以下矩阵形式:<maths num="0005"><![CDATA[<math><mrow><mi>H</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>I</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>H</mi><mn>02</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>其中:I为单位矩阵,H<sub>01</sub>的行数为自由时域长度。
地址 200240 上海市闵行区东川路800号