基于输入趋势补偿的快速模型预测控制方法

基于输入趋势补偿的快速模型预测控制方法
【专利摘要】本发明提供一种提高模型预测控制动态响应速度的控制方法，属于工业自动控制领域；该方法用于进一步提高模型预测控制系统的响应快速性，缩短系统达到稳态时间。该控制方法主要利用每一步的最优控制序列，计算出基于控制输入的变化趋势补偿当下时刻的控制量，从而使系统达到提高快速性的目的。仿真实验表明该方法可行，至少可以使系统的响应快速性提高21％。
【专利说明】
基于输入趋势补偿的快速模型预测控制方法
技术领域
[0001] 本发明涉及工业自动控制领域，具体涉及一种基于输入趋势补偿的模型预测控制 新方法。
【背景技术】
[0002] 模型预测控制(Model predictive control，MPC)是近年来广泛研究的一种反馈 控制策略，是一种除PID控制算法以外，工业界应用最多的一种新型智能控制方法。模型预 测控制的机理可以描述为:在每一采样时刻，根据获取系统当前的状态和输入信息，在线求 解一个有限时域开环优化问题，并且将得到的控制序列的第一个元素作用于系统的控制对 象，在下一采样时刻，重复上述过程，不断用新的预测值刷新优化问题并求得最新解。在线 求解开环优化问题，获得开环优化序列是模型预测控制与其他计算控制率算法主要的区 另IJ。随着计算机计算性能的不断提升，模型预测控制实现更加简单，应用更加广泛。
[0003] 模型预测控制是一种基于模型的闭环优化控制策略，其算法具有控制效果好、鲁 棒性强等优点，可以有效地克服过程的不确定性、非线性和耦合性，并能方便的处理过程被 控变量和操纵变量中的各种约束。模型控制方法可具体分为三类：1)模型算法控制(Model Algorithm Control,MAC)或者是模型预测启发控制（Model Predictive Heuristic Control)。该方法是Richalet等人提出的一种控制算法，算法的实现主要利用系统脉冲响 应，并且只适用于渐进稳定的线性系统。2)动态矩阵控制（Dynamic Matrix Control，DMC) 是一种基于对象阶跃响应的预测控制算法，通过在线求解线性规划(Linear Programming， LP)或者二次规划(Quadratic Programming，QP)问题得到满足约束条件的最优性能，该算 法适用于处理渐进稳定的线性系统，对于弱非线性系统，由于系统在稳定工作点附近具有 近似的线性特征，因而在稳态工作点处作单位阶跃响应实验，得到系统的单位阶跃响应模 型。3)广义预测控制方法，该方法由Clarke等人在保持最小方差自校正控制原理的基础上， 汲取DMC和MAC中的多步预测优化策略，提出广义预测控制（Generalized Predictive Control，GPC)是一种基于辨识被控过程参数模型，并且带有自适应机制的预测控制算法。 4)基于状态空间的模型预测控制方法，该方法基于现代控制理论的状态空间模型设计，通 过动态优化可以解决线性系统、非线性系统的控制问题，其使用范围更为广泛，是当前研究 的一个热点。各类模型预测控制算法虽然在模型、控制和性能上存在许多差异，但其核心都 是基于滚动时域原理，算法中包含了预测模型、滚动优化和反馈校正三个基本原理。
[0004] 模型预测控制(MPC)是最早应用于工业控制领域，理论上的深入研究直到近年才 开展起来。以Morar i的著作为标志，MPC理论研究逐渐成为控制理论研究的热点。Morar i和 Garcia、Rawling、Muske和Mayne等代表学者基于状态空间对于模型预测控制提供了基本概 念和理论框架。MPC存在的稳定性和计算负担的问题也得到了很好的解决。模型预测控制一 直被应用于工业控制等慢变的系统，在目前研究中对于控制系统的响应特性还没有引起足 够的重视，特别随着计算机计算性能的不断提高，模型预测控制算法能够被扩展处理一些 快变系统的控制。比如对于飞行器姿态、水下航行器舵等对于快速性要求很高的场合，还缺 乏相应的理论研究。本发明为了解决某一些快变系统控制时响应慢的问题，将提供一种利 用输入趋势补偿控制输入的模型预测控制方法，用于提高控制系统响应快速性，缩短系统 达到稳态时间。

【发明内容】

[0005] 要解决的技术问题
[0006] 为了提高工业模型预测控制方法的响应特性，本发明提供一种利用输入趋势补偿 控制输入的模型预测控制方法，用于提高控制系统响应快速性，缩短系统达到稳态时间。
[0007] 技术方案
[0008] 本发明的技术构思为:根据k时刻的最优控制序列中预测输入f(k+l |k)对当前输 入lAklk)的影响趋势，对当前输入lAklk)进行补偿，提高控制输入的强度和精准性，进而 提高系统的响应速度。
[0009] -种基于输入趋势补偿的快速模型预测控制方法，其特征在于步骤如下：
[0010] 步骤1:建立系统的离散状态空间模型：
[0011] x(k+l)=Ax(k)+Bu(k) (1)
[0012] 其中x(k)表示k时刻系统的状态矢量，长度为n;x(k+l)表示k+1时刻系统的状态矢 量;u(k)表不k时刻系统的控制矢量，长度为m; A为系统矩阵，B为系统的控制输入矩阵；
[0013] 步骤2:构建系统的二次性能指标：
[0014]
[0015] 其中x(k+i|k)表示在k时刻对系统状态在k+i时刻的预测;u(k+i |k)表示在k时刻 对系统输入在k+ i时刻的预测；Q表不系统状态量的加权矩阵，R表不系统控制量的加权矩 阵，P表示系统终端的加权矩阵，Q、R、P矩阵均为正定矩阵，N表示预测长度;x(k+N)表示在k 时刻对系统状态在k+N时刻到无穷时刻的预测量；
[0016]步骤3:求解最优的控制序列：
[0018] x(k+i+l | k) =Ax(k+i | k)+Bu(k+i | k), (4)
[0019] umin<u(k+i |kXumax，i = 0, · · ·，N-1， (5)
[0020] 其中lf(k) = lVT(k|k)，. . .，ιι,Ι?+Ν-Ι |1〇]1'表示控制输入最优序列;111^11表示控制 输入允许的最小值输入量，u max表示控制输入允许的最大值输入量；
[0021 ]步骤4:根据最优序列趋势，计算输入补偿量：
[0022] Au(k) =y(u*(k| k)-u*(k+l | k)) (6)
[0023] 其中参数μ(.)是趋势补偿函数；
[0024] 1^仏|1〇为步骤3中所求得最优的控制序列第一个控制量，1^仏+1|1〇为步骤3中所 求得最优的控制序列第二个控制量；
[0025] 步骤5:根据补偿量，计算k时刻的最终控制输入：
[0026] u(k)=u*(k| k)+Au(k) (7)
[0027] 步骤6:将最终控制输入u(k)代入步骤1中得到x(k+l)。
[0028] 有益效果
[0029] 本发明提出的一种基于输入趋势补偿的快速模型预测控制方法，实时性强，不增 加模型预测控制算法的复杂度和计算量;补偿操作灵活简单，方便实用;可以显著提高模型 预测控制方法的响应速度，缩短调控时间。在实现算例中，快速性可至少提高21%。
【附图说明】
[0030] 图1为原始优化算法流程图与本发明改进的流程图
[0031] 图2为不加额外影响因子量的系统状态图
[0032] 图3为加入额外影响因子量的系统状态图
[0033] 图4为不加额外影响因子量的系统输入图
[0034] 图5为加入额外影响因子量的系统输入图
[0035] 图6为系统引入的额外影响因子量示意图
【具体实施方式】
[0036]现结合实施例、附图对本发明作进一步描述：
[0037] 1)通过对系统(水下航行器舵系统)进行数学建模分析，建立系统离散状态空间模 型。离散时间周期T取0.1s。
[0038] x(k+l) =Ax(k)+Bu(k) (1)
[0040] 2)对上述系统构造二次性能指标：
[0041]
[0042] 其中状态初始时刻值为
，状态的加权矩阵Q = diag([l，l，l，l])，输入 的加权矩阵为R = diag( [ 1，1，1，1 ])，终端状态的加权矩阵为P = diag( [0 · 5，0 · 5，0 · 5， 0.5])，Q、R、P矩阵均为正定矩阵，在此实例中取预测长度N=9。
[0043] 3)求解最优的控制序列：
[0045] x(k+i+l | k) =Ax(k+i | k)+Bu(k+i | k), (4)
[0046] Umin<U(k+i |kXUmax，i=0, · · ·，N-1， (5)
[0047] 其中U*(k) = |VT(k|k)，. . ·，1^(1^+8110]1??^控制输入最优序列;Umin表示控制输 入允许的最小值输入量，^^表示控制输入允许的最大值输入量；（5)式按元素满足不等式， 在此实例中对(5)式元素不作额外约束要求。
[0048] 4)根据最优序列趋势，计算输入补偿量：
[0049] Au(k)=y(u*(k|k)-u*(k+l|k)) (6)
[0050] 在此实例中取趋势补偿值参数μ = 2。
[0051] 5)根据补偿量，设计k时刻的最终控制输入：
[0052] u(k)=u*(k| k)+Au(k) (7)
[0053] 6)将公式(7)的优化结果带入步骤1)中，计算得到下一时刻的状态量。重复滚动迭 代计算100步后，结束仿真程序。
[0054]表1不加额外影响因子量和加入额外影响因子量的系统性能分析图表
[0055]
[0056] 通过对上述表格分析，可以看出加入额外影响因子量的系统达到稳态5%的响应 时间要比不加额外影响因子量的系统达到稳态5%的响应时间短，继而显示了加入额外影 响因子量的系统可以缩短系统达到稳态的时间，则证明此方法可以解决某一些快变系统控 制时响应慢的问题。
[0057] 通过对图2和图3进行比较，可以得出在相同时间内，加入额外影响因子量的系统 状态达到稳态的收敛速度比不加额外影响因子量的系统达到稳态的收敛速度要快，从而证 明了利用输入趋势补偿控制输入的模型预测控制方法，对于解决某一些快变系统控制时响 应慢的问题是可行的。
【主权项】
1. 一种基于输入趋势补偿的快速模型预测控制方法，其特征在于步骤如下： 步骤1:建立系统的离散状态空间模型： x(k+l)=Ax 化)+Bu 化） （1) 其中X化)表示k时刻系统的状态矢量，长度为n;x化+1)表示k+1时刻系统的状态矢量;U 化)表示k时刻系统的控制矢量，长度为m;A为系统矩阵，B为系统的控制输入矩阵； 步骤2:构建系统的^次性能指柄：其中X化+i I k)表示在k时刻对系统状态在k+i时刻的预测;u(k+i I k)表示在k时刻对系 统输入在k+i时刻的预测;Q表示系统状态量的加权矩阵，R表示系统控制量的加权矩阵，P表 示系统终端的加权矩阵，Q、R、P矩阵均为正定矩阵，N表示预测长度;X化+N)表示在k时刻对 系统状态在k+N时刻到无穷时刻的预测量； 步骤3:求解最优的巧制序列：其中U*化）= [U*T化|k), . . .，U*T化+N-1 |k)]T表示控制输入最优序列;Umin表示控制输入 允许的最小值输入量，Umax表示控制输入允许的最大值输入量； 步骤4 :根据最优序列趋势，计算输入补偿量： Au(k) =μ(ιι*化 I k)-u*(k+l I k)) (6) 其中参数μ(.)是趋势补偿函数； lAklk)为步骤3中所求得最优的控制序列第一个控制量，iAk+l|k)为步骤3中所求得 最优的控制序列第二个控制量； 步骤5:根据补偿量，计算k时刻的最终控制输入： u(k)=u* 化 |k)+Au(k) (7) 步骤6:将最终控制输入U化)代入步骤1中得到X化+1)。
【文档编号】G05B13/04GK106094524SQ201610532905
【公开日】2016年11月9日
【申请日】2016年7月7日
【发明人】李慧平, 范真元, 王婉青
【申请人】西北工业大学