A Comparative Study on Algorithms of Robust and Stochastic MPC for Uncertain Systems
-
摘要: 近几十年来,不确定系统模型预测控制的理论和应用得到了飞速发展.本文简要地回顾了不确定系统中鲁棒模型预测控制和随机模型预测控制的发展历史,总结了它们的相关应用,并较为细致地分析了线性不确定系统模型预测控制的各种主要算法.通过总结各种算法的通用模型、运作方式、问题规模,以及它们保证递归可行性、稳定性的方法,分析了部分算法可行域间的关系,揭示了各种算法的主要特点、适用场合和未来可发展方向,并通过仿真实例直观地分析了各种算法的性能和可靠性.Abstract: In recent years, the development of model predictive control (MPC) of uncertain systems has been remarkable. This paper briefly reviews the development of robust MPC and stochastic MPC, summarizes their applications, and expounds and discusses the main algorithms of linear uncertain systems in these two fields. By summarizing their general models, the ways they work, computational complexities, and the ideas they use to ensure recursive feasibility and stability, we reveal the relationship of feasible sets among some of them and unravel the main features of these algorithms and their application situations. Finally, we demonstrate the performance of all the algorithms through certain simulation cases and give some indication on the future development of these two fields.
-
Key words:
- Uncertain system /
- robust MPC /
- stochastic MPC /
- stability /
- feasibility
-
在控制与决策等领域中, 往往会遇到许多优化问题, 有时我们不仅需要找到全局最优解, 还需要找到局部最优解进行辅助决策, 这类优化问题称为多峰搜索 (Multi-peaksearching) 或多峰优化 (Multimodal optimization) 问题.随着问题的复杂性不断提升, 传统的数学规划理论很难解决.
大规模可并行的群智能算法的出现为解决多峰优化问题带来了新的思路.其中粒子群算法 (Particleswarm optimization, PSO) 是一种基于群体协作的随机搜索算法, 但由于其对局部最优解搜索能力不好并且需要较多的先验知识, 许多学者对其进行了改进[1-4].同样, 差分进化算法 (Differential evolution, DE) 作为进化算法的一个重要分支, 收敛速度比PSO快, 但由于贪婪选择策略, 在求解具有大量局部优点的问题时, 种群极易陷入局部最优.有学者使用临近衰减策略[5]和双目标[6]等策略来解决这个问题.这些算法能够在局部找到接近最优的极值点, 但无法给出所有的极值点.也有学者尝试使用蚁群算法[7]和在极值点附近循环扰动的Nelder-Mead单纯形法[8]来找出所有极值, 但采用的测试函数极值点均较少, 而且搜索时间较长, 不大适合用于搜索有极多极值的函数.
量子算法以其完备的理论和强大的并行能力[9], 被很多学者用于群智能算法中.文献[10]在2010年提出了模拟量子谐振子的概念, 文献[11]中给出了算法实现, 即多尺度量子谐振子算法 (Multi-scale quantum harmonic oscillatoralgorithm, MQHOA).该算法数学原理简单物理含义明确且无需先验知识.完备的采样能力和精确的搜索聚焦能力使MQHOA能较为准确地在多维空间实现对最优中把MQHOA与模拟退火算法 (Simula-ted annealing, SA)、量子解的搜索.文献[12]将退火算法 (Quantumannealing, QA)、量子粒子群算法 (Quantum-behaved particle swarmoptimization, QPSO) 进行对比, 对15种不同的标准测试函数求最小值, 实验结果进一步说明了其具有更好的计算速度和计算精度.
MQHOA集中处理单峰优化问题, 主要根据量子隧道效应[12]保证全局采样点数目足够大以避免陷入局部较优解位置.算法收敛的依据是这些采样点都收敛到一个最优解上, 但对于多峰函数, 采样点将分散在各个峰上而得不到聚集, 导致算法无法收敛.为扩展MQHOA对多峰优化问题的应用, 本文提出了基于划分的多尺度量子谐振子算法 (Partition-based MQHOA, P-MQHOA), 它把求解域均匀划分成N个局部域, 把MQHOA的收敛特性应用在局部域内, 每个局部域中心内放置一个采样点进行谐振子迭代操作, 那么这些采样中心点的平衡位置即为所求的解.实验使用可调节参数的双势阱量子系统函数来验证算法对于等值极值点和不等极值点寻优准确性.对于全部极值点的搜寻效果, 本算法对比了蚁群算法对同一个函数的全极值搜寻, 实验结果超出了蚁群算法.此外, 通过对14个基准函数的全局多峰优化, 对比遗传算法、粒子群算法和差分进化算法效果都要好.当极值点极多的情况下也可以采用多级分辨率策略求解, 求解精度高, 查全率高, 速度快.
1. 多峰优化问题描述
工程上的多峰优化是把这个工程问题映射到一个多维函数上, 要得到多个最优方案就是在这个多维函数上求多个极大值 (或者是求这个多维函数取相反数的多个极小值).
为了全文描述方便, 本文给出以下四个定义:
定义 1.对于某 $D$ 维N极值函数 ( $N\ge2$ ), 第 $i$ 个极值点和第 $j$ 个极值点相邻, 他们坐标之间第 $k$ 维差的绝对值记作 $\lambda_{i, j, k}$ .取 $\lambda_{i, j}=\min\nolimits_{1\le k\le D}\lambda_{i, j, k}, \lambda_{\rm min}=\min\nolimits_{1\le i < j\le N}\lambda_{i, j}$ , 把 $\lambda_{\rm min}$ 称为该函数实际波长分辨率.
定义 2.工程上需要识别某 $D$ 维多极值函数上的极值点满足如下要求: $\lambda_{i, j}\ge \lambda_0\ge \lambda_{\rm min}$ .把 $\lambda_0$ 称为需求波长分辨率, 这些极值点称之为满足需求波长分辨率 $\lambda_0$ 的极值点.
定义 3.把算法求得的极值点坐标与对应的实际极值点坐标之间允许的各维最大差的绝对值称作该算法求解的位置精度, 记为 $A_l$ .
定义 4.某 $D$ 维函数定义域上界为 ${\pmb u}$ , 下界为 ${\pmb s}$ , 那么称 ${\pmb l}={\pmb u}-{\pmb s}=(l_1, l_2, \cdots, l_D)$ 为该函数的长度域, 其中 $l_k$ 为该函数在第 $k$ 维上定义域的长度.
多峰优化问题的数学描述即寻找 $D$ 维函数 $f (x_1, x_2, \cdots, x_D)$ 在定义域上某点 $X_i (x_{i, 1}, x_{i, 2}, \cdots, x_{i, D})$ 处取得极值. $X_i^-=(x_{i, 1}-dx_1, x_{i, 2}-dx_2, \cdots, x_{i, D}-dx_D), X_i^+=(x_{i, 1}+dx_1, x_{i, 2}+dx_2, \cdots, x_{i, D}+dx_D)$ , $f_i^-$ 为函数在点 $X_i^-$ 求得的值, $f_i^+$ 为函数在点 $X_i^+$ 求得的值, $f_i$ 为函数在点 $X_i$ 取得的函数值, $X_i^*$ 为实际极值点, 那么判断该函数在 $X_i$ 处取得极小值的依据:
$\begin{align}\left\{\begin{array}{ll} f_i^->f_i, f_i^+>f_i \\ X_i\in\{X\mid\left|X-X_i^*\right|\le A_l \} \\ \end{array} \right.\end{align}$
(1) 判断该函数在 $X_0$ 处取得极大值的依据:
$\begin{align}\left\{ \begin{array}{ll} f_i^-<f_i, f_i^+<f_i \\ X_i\in\{X\mid\left|X-X_i^*\right|\le A_l \} \\ \end{array} \right.\end{align}$
(2) 2. 基于划分的MQHOA原理
量子谐振子是物理上描述微粒动态振动的物理量, 可以近似描述分子的运动规律.MQHOA算法把量子谐振子群作为智能粒子群来解决优化问题, 振子的平衡位置即势能最低位置, 就是所求函数的最优解.量子谐振子从高能态到低能态是一个逐渐收敛的过程[13], 直到收敛到基能态: $\psi_0(x)=\frac{\sqrt{\alpha}}{\sqrt[4]{\pi}}{\rm e}^{-\alpha^2\frac{(x-x_i)^2}{2}}$ , 定义MQHOA算法在尺度 $\sigma_s$ 下的波函数为 $k$ 个以 $x_i$ 为中心的高斯概率密度函数的迭加:
$\begin{align}\psi_{QHO}(x)=\sum\limits_{i =1}^k{\frac{1}{\sqrt{2\pi}\sigma_s}}{\rm e}^{-\frac{(x-x_i)^2}{2\sigma_s^2}}\end{align}$
(3) $\psi_{QHO}(x)$ 近似代表了以 $k$ 个最优解为中心的目标函数可行域上的概率分布, 通过引入多尺度优化函数二进信息高斯采样模型, 二进小波尺度函数采用高斯函数按照不同尺度和不同精度要求进行聚焦搜索[9], 可以保证信息的不遗漏采样, 从而获得全局最优解, 详细算法描述参见文献[11].
P-MQHOA扩展了MQHOA对多峰优化问题的处理能力, 并且具有对所有最优解的记忆能力和更快的收敛速度.其主要思想是把MQHOA的全局优化特性局限在了一个较小划分区域内, 在这个区域内只设置一个高斯采样中心点, 并且让这个采样中心点在反复迭代过程中不断收敛到精度范围内.P-MQHOA对整个函数定义域划分的N个区域与N个采样中心点一一对应, 最终找到N个解 (包括重复解).由于这N个采样中心点之间的采样计算是独立进行的, 并不会由于某几个采样中心点收敛不了而影响其余采样中心点, 因此其收敛速度要比MQHOA要快得多, 并且算法通过对记忆的N个解进行整理可以得到全局最优解.
算法对函数定义域进行均匀划分, 划分区域内如何选择要生成的高斯采样点的标准差直接影响算法的性能和执行结果.各维度划分成的各均匀区段中心位置处随机生成 $M$ 个标准差为 $\sigma_s$ 的高斯分布点, 在同一个区段内, $\sigma_s$ 与这些分布点有如下关系: $\sigma_s$ 越大分布点越接近均匀分布, $\sigma_s$ 越小分布点越集中在中心点, 通过 $\sigma_s$ 的衰减可以使得这些分布点最终聚集在搜索精度范围内.算法初始的 $\sigma_s$ 大小很重要, 分布点太均匀的话很大概率会搜索到临近的其他区段, 分布点太集中的话很大概率会忽略区域边界附近的极值点.假如某个区段内初次生成的高斯分布点超出边界的概率为 $p$ , $M$ 次伯努利实验后可得到超出边界的分布点数目平均值为 $Mp$ .算法一般选择 $\sigma_s$ 为该区段长度的四分之一, 通过查找正态分布表可知, $p=0.0456$ , 即每个分布点超出区段的概率为0.0456.这样设置不仅可以保证对本区段内的极值搜索, 又不失高斯分布的中心点集中特性.
对定义域下界为 ${\pmb s}$ 长度域为 ${\pmb l}$ 的 $D$ 维函数在第 $j$ 维上进行均匀划分, 划分成 $k$ 个区段后高斯概率密度叠加可得第 $j$ 维上采样点的概率分布:
$\begin{align}\psi_{jQHO}(x)=\sum\limits_{i =1}^k{\frac{1}{\sqrt{2\pi}\sigma_s}}{\rm e}^{-\frac{\left(x-\frac{l_i(i-\frac{1}{2})}{k}-s_i\right)^2}{2\sigma_s^2}}\end{align}$
(4) 其对应的密度图如图 1所示, 虚线为原来的高斯概率分布, 实线为叠加后的概率分布.
叠加后某一区段的高斯分布由于临近区段分布点的 $p$ 概率跳入而弥补了本区段 $p$ 概率跳出流失的分布点, 这样就保证了 $M$ 个采样点在各区段数目相等, 并且有 $Mp$ 个点作为区段之间的交流媒介保证了各区段的联系.这种联系使得在某区段内找不到极值点的情况下可以跳出该区域协助临近区段寻找极值点.把各个维度上的均匀划分区段组成一个个均匀划分区域, 在各个划分区域内寻找极值是相互独立的操作, 故算法可以多线程并行执行.
为了保证在需求波长分辨率 $\lambda_0$ 下识别该分辨率下的所有极值点, 本算法需要在可行域上划分 $K$ 个区域并设置 $K$ 个采样中心点, $K$ 的大小决定算法能准确寻找极值的数目.以寻找一维函数极小值为例, 定义搜索区域长度为函数定义域上均匀划分的区段长度.如图 2所示, $\psi (x)$ 为初始采样点在区域内分布的概率密度函数.当搜索区域长度超出该一维函数的波长分辨率 $\lambda_0$ (两个波谷之间的水平距离) 时, 算法在当前搜索区域内至少能找到一个极值点, 如图中函数 $f_1 (x)$ 所示; 当搜索区域长度等于 $\lambda_0$ 时, 极值点或者在该区域被找到或者在邻近区域被找到, 如图中函数 $f_2 (x)$ 所示; 当搜索区域长度小于 $\lambda_0$ 时, 当前搜索区域无法找到极值点, 只能在邻近区域找到, 如图中函数 $f_3(x)$ 所示.因此, 本算法需要通过设置合理的采样中心点数使得搜索区域长度大于等于 $\lambda_0$ .
$K$ 的大小也影响算法的执行速度.以某 $D$ 维函数为例, 其可行域 ${\pmb l}=(l_1, l_2, \cdots, l_D)$ , 如果需要找到所有极值, 那么搜索区域长度正好为该函数的实际波长分辨率 $\lambda_{\rm min}$ , 以此算出合理的采样中心点数.在第 $i$ 维上置采样中心点个数 $k_i$ 为 $\lceil{ \frac{l_i}{\lambda_{\rm min}}}\rceil $ , 即在该维度上划分出 $k_i$ 个区段, 要找到 $D$ 维函数所有极值点则需初始划分区域数 $K=\prod\nolimits_{i=1}^Dk_i$ .当 $\lambda_{\rm min}$ 很小而 $D$ 很大时, 算法需要花费极多CPU时间在所有区域内进行计算.类似多分辨率的小波分析[14], 本文提出多级分辨率策略, 其基本思想是先在大波长分辨率下进行区域划分, 在找到极值点的区域内再进行小波长分辨率的划分, 而在不存在极值点的区域就不再进行划分计算, 如此重复可以使搜索区域长度逼近 $\lambda_{\rm min}$ , 从而找到满足需求分辨率极值点.这种策略避免了原来固定分辨率下对无极值区域的计算浪费, 而把更多的CPU计算集中在极值较多的区域.实际应用中需求波长分辨率 $\lambda_0$ 一般比 $\lambda_{\rm min}$ 大得多, $K$ 取 $\prod\nolimits_{i=1}^D\left\lceil\frac{\lambda_0}{l_i}\right\rceil$ .当 $K$ 不是很大时, 可以直接采用固定分辨率策略, 反之采用多级分辨率策略.
3. P-MQHOA算法过程描述
3.1 P-MQHOA算法的三个阶段
P-MQHOA主要有三个阶段.初始划分阶段:该阶段主要是对 $D$ 维度函数定义域进行区域均匀划分, 然后在区域中心放置一个采样中心点; 量子谐振子收敛阶段:该阶段采用高斯函数来进行采样最终定位中心点; 归并阶段:该阶段把位置精度范围内的中心点归为一类, 并从该类中选出最好的中心点.
3.1.1 初始划分阶段
输入需求波长分辨率 $\lambda_0$ 下所求 $D$ 维度函数定义域上界为 ${\pmb u}$ , 下界为 ${\pmb s}$ , 由定义4可求得长度域 ${\pmb l}$ .对于第 $j$ 维划分区段数记作 $k_j$ , 若算法采取固定分辨率策略, 初始第 $j$ 维上划分区段数 $k_j$ 为 $\lceil\frac{l_i}{\lambda_0}\rceil$ ; 若算法采取N级分辨率策略, 初始第 $j$ 维上划分区段数为 $\lceil\sqrt[N]{\frac{l_j}{\lambda_0}}\rceil$ .总采样点数 $K=\prod\limits_{j=1}^Dk_i$ , 初始化采样中心点, 某个中心点在第 $i$ 个区段第 $j$ 维的坐标值可以记作 $c_{i, j}=\frac{l_j (i-1/2)}{k_i}i+s_j, i\in[1, k_j], j\in[1, D]$ , $i$ 与 $j$ 均为整数.
3.1.2 量子谐振子收敛阶段
基于划分的MQHOA的采样中心点之间互不影响, 在每个划分区域内的收敛过程都是独立的, 因此各中心点可以并行执行.如图 3所示为寻找某一维函数 $f (x)$ 极大值的收敛过程.初始过程的采样点概率分布为6个高斯分布函数的叠加, 见式 (4).其区域长度正好为这个叠加函数的周期长度.按照前面说明, 该区域的高斯函数提供 $M (1-p)$ 个采样点, 由于在可行域外产生的采样点很少可以忽略不计, 那么很容易证明该区域内约有 $Mp$ 个采样点由其他区域的高斯采样提供, 而越靠近的区域提供越多.在图 3第一步中左起第二个区域内约有 $Mp/2$ 个采样点落在第三个区域, 而第三个区域内的采样点效果显然比原区域好, 因此下一轮采样将在第三个区域进行.在找到当前区域最大值采样点后, 以这个点为中心, 搜索区域长度衰减, 重新在新的搜索区域寻找最大值, 如图 3中第二第三步, 以此类推, 直到所有中心点要么收敛在极值点附近, 要么找不到极值点.图 3最后一步中所有中心点都收敛在极值点附近.
3.1.3 归并阶段
图 3最后收敛阶段有6个中心点分别收敛于4个波峰, 其中第2、3、4个中心点聚集在一个极值点附近, 需要对这些点进行归并操作.若任两个点各维坐标最大差的绝对值小于位置精度 $A_l$ , 按照定义3, 把这两个点归为一类, 最终从此类中取出对应函数值最大 (最小) 的点作为最终的极值点.
3.2 P-MQHOA算法流程
综上所述, 整个算法流程可以简要描述如下:
步骤 1.初始化:
步骤 1.1.输入 $\lambda_0, M, A_l, {\pmb u}, {\pmb s}$ .
步骤 1.2.根据公式初始化长度域 ${\pmb l}$ , 标准差 $\sigma_s=\frac{\lambda_0}{4}$ , 在各维度上均匀划分产生 $K$ 个中心点, 其中第 $i$ 个点 $center_i=(c_{i, 1}, c_{i, 2}, \cdots, c_{i, D})$ .
步骤 2.收敛过程:
步骤 2.1.在各维上以 $c_i$ 为中心点, $\sigma_s$ 为标准差进行 $M$ 次高斯采样, 当采样点超出可行域边界时以边界点为准.
步骤 2.2.从 $M$ 个采样点中选择对应函数值最大 (最小) 的一个点作为新的中心点, $\sigma_s$ 衰减为原来一半.如果满足 $\sigma_s\geq A_l$ , 重复执行步骤2.1, 否则继续执行步骤3.1.
步骤 3.归并:
步骤 3.1.对最后生成的 $K$ 个解根据式 (1) 或式 (2) 作极值判定, 去除不满足式子的点.
步骤 3.2.计算各相邻点之间各维度差值的最小值, 若小于 $A_l$ 则合并这些点并取其中最大 (最小) 的作为最终解.
3.3 P-MQHOA算法两种策略
针对维度较低、需求波长分辨率不高的函数进行多峰寻优, 采用固定分辨率策略即可解决.但维度较高或者需求波长分辨率很小情况下, 应当采用多级分辨率策略.
对于N级分辨率, 算法经历初始划分阶段和量子谐振子收敛阶段后, 找到输出中心点不为空所对应的区域, 对这些区域执行下一级的初始划分阶段的操作并对划分的子区域执行量子谐振子收敛阶段的操作, 递归搜寻到第N级后对所有输出的中心点执行归并阶段的操作, 可输出最终的极值点.
以图 4为例, 对 $f (x, y)$ 在需求分辨率为0.1下求极值, 长度域 $l_x=l_y=5$ , 采用三级分辨率策略初始划分区域数为 $\lceil\sqrt[3]{\frac{5}{0.1}}\times\sqrt[3]{\frac{5}{0.1}}\rceil=16$ , 如子图 $4 (a)$ 所示; 经过收敛和归并过程后确认极值点在6号与11号区域, 重新在这两个区域进行第二级划分, 如子图 4(b)所示; 重复以上过程第三级划分后实际分辨精度可达到 $\frac{5}{4^{3}}=0.078$ , 满足需求分辨率终止迭代, 如子图 4(c)所示.若每个中心点的采样次数为 $M$ , 可以计算总的采样次数为 $16(1+2+3) M=96M$ , 而采取固定分辨率需要次数为 $50\times50M=2 500M$ , 可见其省去很多计算时间, 但多级数情况下也容易丢失部分极值点.因为 $M$ 在粗分辨率下相对于细分辨率就要小, 就好比把黄豆随机撒在球场上和撒在球门里的区别.如子图 4(b)所示若11号区域未找到极值点那么在子图 4(c)中就会丢失很多极值点.因此 $M$ 的选择和级数的选择应该根据实际情况而定.
此外, 多级分辨率策略下不仅在同级区域内横向可并行执行, 不同级区域也可纵向并行执行, 计算过程不会相互干涉, 在保持原QHOA算法低耗时、快收敛、高精度的基础上该算法有很好的多峰寻优能力.
4. 算法实例分析
本节通过寻优精确性、全极值点寻优和全局多峰优化三个角度来测试该算法的综合性能, 实验使用的计算机采用Mac系统, 主频29 GHz, 在JDK7环境下的浮点计算精度为小数点后16位, 实验寻找的极值点默认为极小值点, 设置位置精度 $A_l$ 为0.000001.
4.1 寻优准确性
准确性是衡量一个算法好坏的重要标准, 为了测试本算法在需求波长分辨率下寻找极值点的准确度, 实验使用双势阱量子系统函数来验证.双势阱量子系统是研究微观物理学的一个重要物理模型, 其函数表达式为
$\begin{align}V_{sys}=V_0\frac{(x^2-a^2)^2}{a^4}+\delta x\end{align}$
(5) 该函数有两个极小值, 其中 $V_0$ 表示函数的值域长度, $2|a|$ 表示两个极小值点之间的水平距离, $2|\delta||a|$ 表示两个极小值点之间的垂直距离.置 $V_0=10$ , $a$ 分别取 $1, 0.5, 0.1, 0.05, \cdots, 0.0000001$ , 对应函数实际分辨率 $\lambda_{\rm min}$ 为 $2, 1, 0.5, 0.1, \cdots, 0.0000005$ .设置函数定义域为 $[-1.5, 1.5]$ , 初始采样中心点数 $K=30$ , 固定分辨率下可识别波长分辨率为 $\lceil\frac{1.5-(-1.5)}{30}\rceil=0.1$ , 实验设置需求波长分辨率 $\lambda_0=0.1$ , 测试算法在实际波长分辨率 $\lambda_{\rm min}$ 取不同值时识别极值的准确度.此外, 为了测试极值点间的垂直距离对寻找极值点准确度的影响, 实验分别取 $\delta=1$ 与0, 前者对应不等值极值双势阱函数, 后者对应等值极值双势阱函数.经过1 000次试验, 统计寻找极值点的次数, 不等值极值点函数统计结果如图 5所示, 等值极值点函数统计结果如图 6所示.其中横坐标为 $a$ , 三条曲线分别代表找到第一个极值点的次数、找到第二个极值点的次数和全部找到极值点次数.
从两张图可以发现本算法可以找出需求分辨率下几乎所有极值.令 $\lambda_0=2|a|$ , 解出 $a_0=0.05$ ; 令 $A_l=2|a|$ , 解出 $a_1=0.0000005. a_0$ 与 $a_1$ 正好是两张图曲线趋势变化的分界点.当 $a < a_0$ 时均能找到两个极值点; 当 $a=a_0$ 时极大概率找到两个极值点; 当 $a_0 < a\leq a_1$ 时, 找到两个极值点的概率趋于稳定; 当 $a>a_1$ 时找到两个极值点的概率逐渐减少.实验结果与前面图 2说明相符合, 证明了需求波长分辨率下寻优准确性.
此外可以发现:
1) 不管对于等极值函数和不等极值函数, 该算法在定义域内均能至少找到一个极值点, 并且随着 $\lambda_{\rm min}$ 减小, 三条曲线变化趋势几乎一致.
2) 图 4中 $a$ 取0.05与0.0005之间时算法更容易找到极值更小的第一个极值点, 也就是说对于不等值极值点, 该算法容易找到极值更明显的点, 这个特性使得该算法适合用于全局优化.
4.2 全极值点寻优
4.2.1 全极值点寻优性能
函数的全部极值点往往可以反映函数整体变化趋势, 大部分算法对全极值点的搜索消耗了很多时间, 而且精度不高, 例如文献[7]中的Ant算法与文献[8]中RePAMO算法全极值点的搜寻成为多峰优化问题上的一个难题, 这里采用两种策略的P-MQHOA作全极值寻优.
实验使用简单一维震荡函数 $y=x{\rm sin}\frac{1}{x}, x\in[-0.5, 0.5]$ , 如图 7.该函数在原点附近有无穷极值, 越靠近原点极值点分布越紧密, 直至波长分辨率趋近于0.故不可能找出理论上的所有极值点, 实验指定需求波长分辨率, 然后找出该波长分辨率下的极值点数.
指定需求分辨率下的理论极值点数目可以这样求.对 $y$ 求导, $y'=-\frac{1}{x}\cos{\frac{1}{x}}+\sin{\frac{1}{x}}$ , 令 $y'=0$ , 有 $\tan{\frac{1}{x}}=\frac{1}{x}$ , 当 $\frac{1}{x}\in[-0.5, 0.5]$ 时可以得到近似解, 其中正近似解 $x_+=\frac{2}{(2k+1)\pi}, k\in N^*$ .相邻两个波谷水平距离 $\lambda=\frac{2}{(2k+1)\pi}-\frac{2}{(2k+5)\pi}$ , 在 $\lambda_0 < \lambda$ 条件下可解得
$\begin{align}k<\sqrt{1+\frac{2}{\lambda_0\pi}}-1.5, k\in N^*\end{align}$
(6) 这是正近似解数, 由于函数对称性, 负近似解数等于正近似解数, 所以满足式 (6) 的最大正整数的两倍记作 $K_0$ , 即为满足需求分辨率的极值点个数.
为验证该算法对多维函数的多极值寻优性能分析, 实验使用Griewank函数, 函数表达形式: $f (x_i)=\sum\limits_{i=1}^D\frac{x_i^2}{4000}-\prod\limits_{i=1}^D\cos\left (\frac{x_i}{\sqrt{i}}\right)+1$ 这里选择维度为2, 如图 8所示, $x_i$ 在[ $-$ 100, 100]内, 该函数有固定等分辨率的极值点数记为 $K_0$ .
实验设置不同需求波长分辨率 $\lambda_0$ , 进行50次独立实验.对比固定分辨率和二级分辨率下算法的性能, 前者找到的极值点个数用 $K_1$ 表示, 后者用 $K_2$ 表示, 平均耗时 (Timeconsuming) 分别为 $TC_1$ 、 $TC_2$ , 单位为秒, 采样中心点数 (Center-point numbers) 分别用 $CN_1$ 、 $CN_2$ 表示.震荡函数结果见表 1.
表 1 一维震荡函数两种策略寻优性能比较Table 1 Comparison of two strategies for one dimensional shock functionK0 λ0 K1 CN1 TC1 K2 CN2 TC2 2 0.1 3.9 10 0.0030 4.0 4 0.0038 12 0.01 13.9 100 0.0180 11.9 10 0.0084 46 0.001 47.7 1 000 0.1499 41.6 31 0.0348 156 0.0001 157.0 10 000 1.312 144.2 100 0.167 500 0.00001 493.2 100 000 10.374 461.0 316 0.767 1592 0.000001 1515.0 1 000 000 95.614 1 419.6 1 000 3.067 $K_0$ 用来表示满足需求波长分辨率 $\lambda_0$ 的极值点数, 实际上算法可以在原点附近找到两个比需求波长分辨率还小的极值点, 如 $\lambda_0=0.1$ 时算法还可以找到分别在 ( $-$ 0.1, 0) 与 (0, 0.1) 域内的两个极值点, 但这两个极值点距离最近极值点的坐标距离远小于 $\lambda_0$ , 因此找到满足 $\lambda_0$ 的极值点数实际为 $K_1-2$ 和 $K_2-2$ .比较 $K_0$ 与 $K_1-2$ 和 $K_2-2 $ 的关系可知:随着 $\lambda_0$ 减小, 算法可以找到更多的极值点, 而且相比多级分辨率策略固定分辨率策略下找到满足 $\lambda_0$ 的极值点数更接近 $K_0$ , 但多级分辨率策略下算法花费时间较短.如取 $\lambda_0$ 为0.000001时固定分辨率与二级分辨率时间消耗相差约三十倍, 其原因是单个采样中心点迭代次数恒定的前提下两种策略算法所用的采样中心点数相差很大, 采用固定分辨率策略共需要 $CN_1=\lceil (\frac{L}{\lambda_0})^D\rceil=\lceil (\frac{(0.5-(-0.5))}{0.000001})^1\rceil=1 000 000$ 个采样中心点, 而采用二级分辨率策略初始只需 $CN_2=\lceil (\frac{L}{\lambda_0})^{D/N}\rceil=\lceil (\frac{(0.5-(-0.5))}{0.000001})^{1/2}\rceil=1 000$ 个采样中心点, 初始平均划分1 000个区域, 这1 000个区域中存在很多无极值的区域, 在第二级划分时不会对这些无极值区域划分, 这样总的采样中心点相对来说就变很少了.同样把函数换成二维Griewank函数进行50次独立实验, 对结果取平均值, 见表 2所示.
表 2 二维Griewank函数两种策略寻优性能比较Table 2 Comparison of two strategies for two dimensional Griewank functionK0 λ0 K1 CN1 TC1 K2 CN2 TC2 1 369 20 75.7 10 0.051 40.3 4 0.036 1 369 10 307.8 400 0.183 231.2 10 0.291 1 369 5 1 049.2 1 600 0.726 617.5 31 1.472 1 369 2 1 333.4 10 000 4.078 670.2 100 1.657 1 369 1 1 340.5 40 000 15.312 783.0 316 9.64 1 369 0.5 1 353.10 160 000 57.272 1 924.1 1 000 21.04 随着 $\lambda_0$ 减小算法找到的极值点变多, 到 $\lambda_0 < 5$ 后两种策略找到的极值点数均开始变得稳定.这是由于Griewank函数极值点坐标均匀分布在各维度上, 各个相邻极值点在各维度上的坐标差均相等, 等于实际波长分辨率, 当需求分辨率小于实际波长分辨率后算法找到的极值点数就趋于稳定了.此外, 采用固定分辨率策略划分的MQHOA几乎能够找到需求波长分辨率下所有极值点, 而二级分辨率策略下查全率较低.如上表所示, 二维Griewank函数在 $\lambda_0$ =2时采用固定分辨率策略几乎能够找到所有极值点, 而采用二级分辨率策略下只能找到实际数目的一半, 这是因为在粗分辨率下忽略的极值点所在区域不会参与到第二级的划分.但整体来说多级分辨率策略可以用较短的时间完成算法搜索, 尤其当需求波长分辨率很小的情况下这种时间差别更明显.
与其他智能算法比较, 本算法在时间和精度上都有很大优势.如文献[5]中利用蚁群算法对函数 $y=x_1^2+x_2^2-\cos (18x_1)-\cos (18x_2)$ , $x_1, x_2\in[-1, 1]$ 求所有极大值点.与利用P-MQHOA算法得出的结果比较如表 3所示.对36个极值点按照各维坐标进行基数排序, No.为找到极值的序号, $((x_1, x_2, f (x_1, x_2))$ 分别为坐标点和对应的值.
表 3 P-MQHOA与蚁群算法全极值搜索比较Table 3 Comparison of total-extremum-searching of P-MQHOA and ant colony algorithmNO. P-MQHOA
(x1, x2, f/(x1, x2))ANT 1 (-0.878093608, -0.878093599, 3.53255484) (-0.8781, -0.8781, 3.5326) 2 (-0.878093550, -0.526852857, 3.042136418) -0.8737, -0.5310, 3.0362 3 (-0.878093617, -0.17561706, 2.796928371) (-0.8725, -0.1810, 2.7873) 4 (-0.878093568, 0.175617049, 2.796928371) (-0.8725, -0.5310, 3.0362) 5 (-0.878093628, 0.526852825, 3.042136418) (-0.8700, 0.5175, 3.0176) 6 (-0.878093611, 0.878093569, 3.53255484) (-0.8675, 0.8675, 3.4966) 7 (-0.52685281, -0.878093573, 3.042136418) (-0.5310, -0.8737, 3.0362) 8 (-0.526852791, -0.526852812, 2.551717996) (-0.5268, -0.5268, 2.5517) 9 (-0.526852794, -0.175617069, 2.306509949) (-0.5352, -0.1773, 2.3056) 10 (-0.526852803, 0.17561705, 2.306509949) (-0.5215, 0.1700, 2.2969) 11 (-0.526852866, 0.526852858, 2.551717996) (-0.5200, 0.5200, 2.5367) 12 (-0.526852827, 0.878093628, 3.042136418) (-0.5175, 0.8700, 3.0176) 13 (-0.175617027, -0.878093592, 2.796928371) (-0.1810, -0.8725, 2.7873) 14 (-0.175617075, -0.526852805, 2.306509949) (-0.1773, -0.5252, 2.3056) 15 (-0.175617028, -0.175617046, 2.061301903) (-0.1756, -0.1756, 2.0613) 16 (-0.175617006, 0.175617028, 2.061301903) (-0.1756, 0.1756, 2.0559) 17 (-0.175617067, 0.526852847, 2.306509949) (-0.1700, 0.5215, 2.2969) 18 (-0.175617056, 0.878093620, 2.796928371) (-0.1675, 0.8700, 2.7759) 19 (0.175617041, -0.878093592, 2.796928371) (0.1675, -0.8700, 2.7759) 20 (0.175617068, -0.526852815, 2.306509949) (0.1700, -0.5215, 2.2969) 21 (0.175617061, -0.175617041, 2.061301903) (0.1715, -0.1715, 2.0559) 22 (0.175617051, 0.175617057, 2.061301903) (0.1756, 0.1756, 2.0613) 23 (0.175617054, 0.526852783, 2.306509949) (0.1773, 0.5252, 2.3056) 24 (0.175617039, 0.878093590, 2.796928371) (0.1810, 0.8725, 2.7873) 25 (0.526852762, -0.878093542, 3.042136418) (0.5175, -0.8700, 3.0176) 26 (0.526852876, -0.526852885, 2.551717996) (0.5200, -0.5200, 2.5367) 27 (0.526852785, -0.175617064, 2.306509949) (0.5215, -0.1700, 2.2969) 28 (0.526852823, 0.175617064, 2.306509949) (0.5252, 0.1773, 2.3056) 29 (0.52685278, 0.52685282, 2.551717996) (0.5268, 0.5268, 2.5517) 30 (0.526852781, 0.878093598, 3.042136418) (0.5310, 0.8737, 3.0362) 31 (0.878093552, -0.878093598, 3.53255484) (0.8675, -0.8675, 3.4966) 32 (0.87809359, -0.526852814, 3.042136418) (0.8700, -0.5175, 3.0176) 33 (0.878093528, -0.175617038, 2.796928371) (0.8700, -0.1675, 2.7759) 34 (0.878093604, 0.175617041, 2.796928371) (0.8725, 0.1810, 2.7873) 35 (0.878093617, 0.526852807, 3.042136418) (0.8737, 0.5310, 3.0362) 36 (0.878093597, 0.878093581, 3.53255484) (0.8781, 0.8781, 3.5326) 以上实验数据是在固定分辨率策略下设置中心点数为200, 可以算出可识别波长分辨率为0.142, 小于函数实际的波长分辨率, 因此可找到所有极值点.另外, 相对于蚁群算法求出所有极值点花费的3 203.297 s, P-MQHOA只需0.123 s, 并且其精度相对于蚁群算法的0.01精确到0.0000005, 无论是时间还是精度上P-MQHOA都占很大优势.实际上, 这里每个中心点只生成20个高斯采样点也能找出所有极值点, 时间花费可以更少, 因此P-MQHOA性能上还有进一步提高的空间.
综上讨论, P-MQHOA适合函数的全极值点寻优, 并且可以根据工程具体需求采用不同的策略:当需要查找所有极值点时可以采用固定分辨率策略, 当对时间要求较高时可以采用多级分辨率策略.实际上实验还发现这两种策略所找到的极值都包含全局最优解, 因此在不考虑求所有极值的情况下, 对于全局优化问题可以考虑用多分辨率策略.
4.3 全局多峰优化
在控制与决策领域中, 例如无人机的紧密编队飞行控制应用、航空发动机的燃烧主动控制应用、导弹控制参数的自动调节应用等[15], 系统的代价函数普遍地存在多个全局极值, 如何能够使系统的代价函数稳定地工作在全局极值点, 最大地发挥系统的性能, 是工业控制界一个亟待解决的问题.大部分多峰优化智能算法致力于解决这个问题, 而基于划分的MQHOA是适合全极值点寻优的算法, 可以从众多极值点中筛选出最佳极值点.
本节实验选用14个函数[6]作为测试集合, 其中 $f_1\sim f_7$ 和 $f_{12}$ 为一维函数, $f_{14}$ 为三维函数, 其余都是二维函数.比较P-MQHOA与文献[6]中其他智能算法的的区别.这些算法包括:基于双目标差分进化的多峰优化算法 (Multimodal optimization using abiobjective differential evolution algorithm, MOBi-DE)、拥挤度差分进化算法 (Crowding DE, CDE)、基于形态的差分进化算法 (Speciation-based DE, SDE)、欧氏适应度距离半径的粒子群算法 (Fitness-Euclidean distance ratio PSO, FER-PSO)、基于物种形成原理的粒子群算法 (Speciation-based PSO, SPSO)、协方差矩阵适应进化算法 (Covariance matrix adaptation evolutionstrategy, CMA-ES)、固定小生境半径的协方差矩阵适应进化算法 (CMA-ESwith fixed niche radius, CMA)、双目标多群体遗传算法 (Biobjective multipopulation genetic algorithm, BMPGA)、基于小生境的带有精英策略的非支配排序的遗传算法 (Niching-basedNSGA-Ⅱ) 和环型拓扑粒子群优化算法 (Particle swarm optimization using aring topology, rpso).其中r2pso中每个粒子只与环形拓扑中右边的粒子交互, r3pso与环形拓扑中的左右两边的粒子均交互.如表 4所示, 所有算法都经过50次独立重复实验, 统计找到的全局峰数, 对结果取平均值, 然后对每个算法进行排名记录到括号中去, 最终统计总排名.其中ε用来控制极值点的精度, No.表示对应函数应有全局峰数. P-MQHOA算法对一维函数采用固定分辨率策略, 二维函数采用二级分辨率策略, 三维函数采用三级分辨率策略.
表 4 二维Griewank函数极值分布拟合结果Table 4 Fitting results of extremum values0 distribution of Griewank functionFunc. No. ε P-MQHOA MOBi-DE CDE SDE S-CMA CMA SPSO PER-PSO r2pso r3pso Niching-based BNPGA
NSGA-Ⅱf1 1 0.05 1(4) 1(4) 1(4) 1(4) 1(4) 1(4) 0.44(12) 0.82(10) 0.76(11) 0.84(9) 1(4) 1(4) f2 1 0.05 1(5) 1(5) 1(5) 1(5) 1(5) 1(5) 0.40(12) 1(5) 0.88(11) 0.96(10) 1(5) 1(5) f3 2 1E-6 2(2.5) 2(2.5) 2(2.5) 1.96(5) 1.92(7) 1.95(6) 1.38(9) 0.8(10) 0.48(12) 0.6(11) 2(2.5) 1.5(8) f4 5 1E-6 5(1.5) 5(1.5) 3.84(10) 4.70(6) 0.04(12) 0.6(11) 4.88(3) 4.84(4) 4.68(7) 4.74(5) 4.37(9) 4.64(8) f5 1 1E-6 1(6) 1(6) 0.72(12) 1(6) 1(6) 1(6) 1(6) 1(6) 1(6) 1(6) 1(6) 1(6) f6 5 1E-6 5(1.5) 5(1.5) 3.96(9) 4.6(7) 0(12) 0.64(11) 4.92(4) 4.96(3) 4.88(5) 4.72(6) 4.52(8) 3.56(10) f7 1 1E-6 1(5.5) 1(5.5) 0.8(11) 1(5.5) 1(5.5) 1(5.5) 1(5.5) 1(5.5) 1(5.5) 1(5.5) 0.6(12) 1(5.5) f8 4 5E-4 4(1.5) 4(1.5) 0.32(12) 3.72(4.5) 3.72(4.5) 3.43(7) 0.84(11) 3.68(6) 2.92(9) 2.76(10) 3.74(3) 3.34(8) f9 2 1E-6 2(3) 2(3) 0.04(12) 2(3) 1.6(7) 2(3) 0.08(11) 1.96(6) 1.44(9) 1.56(8) 2(3) 1.24(10) f10 1 1E-6 1(2) 1(2) 0.52(9) 0.32(10) 0.06(12) 0.18(11) 0.56(8) 1(2) 0.88(6) 0.76(7) 0.94(5) 0.96(4) f11 18 5 E-2 17.5(1) 17.4(2) 11.39(12) 12.78(9) 12.04(10) 12(11) 14.36(7) 15.61(6) 15.95(5) 16.45(4) 16.94(3) 13.72(8) f12 6 0.05 6(1.5) 6(1.5) 5.56(6.5) 4.88(12) 5.56(6.5) 5.81(3) 5.6(5) 5.28(10) 5.52(8) 5.16(11) 5.36(9) 5.71(4) f13 36 1E-3 35.6(1) 35.4(2) 33.8(3) 23.8(7) 24.6(6) 23.6(8) 25.72(5) 21.8(12) 22.4(9) 22.2(10) 22.05(11) 31.76(4) f14 218 1E-3 163.56(2) 175.88(1) 152(3) 50.6(8) 0.6(12) 32.5(11) 70.12(6) 68.6(7) 40.6(10) 45.4(9) 92.76(5) 121.54(4) Total ranks 38 39 111 92 109.5 102.5 104.5 92.5 113.5 111.5 85.5 88.5 参与比较的算法主要是对差分进化算法 (DE)、遗传算法 (GA) 和粒子群优化算法 (PSO) 的各种改进算法, 属于进化算法的分枝, 也都是基于群体智能的随机并行优化算法, 通过模仿生物群体内个体间的合作与竞争产生的启发式群体智能来指导优化搜索.这些算法使用小生境来提高种群多样性避免早熟, 可以同时使用多个子种群来定位和跟踪多个最优解.小生境的半径和小生境的数量影响整体的运算速度和精度, 算法性能对该参数很敏感, 需要调试多次, 其中大部分改进算法需要较多控制参数才能达到比较满意的结果, 而P-MQHOA只需要设置与需求波长分辨率相对应的采样中心点数即可.根据表 4统计的排名数据绘制热力图, 如图 9所示, 其中同一行颜色越淡表示找到的全局最优点越多, 反之越少.比如对于 $f_1$ 函数, 由颜色深浅可见SPSO算法的寻优效果最差, 其他PSO算法表现也不佳, 而P-MQHOA与其他遗传算法均能找到对应的全局最优点.
从上图可看出, P-MQHOA与MOBi-DE对14个函数测试下所找到的最优解数目差不多, 而其他算法相比之下较为逊色.单独比较两种算法对这14个函数测试所花费的平均时间, 如表 5所示, 时间单位为秒.
表 5 P-MQHOA算法与MOBi-DE算法时间消耗 (s)Table 5 Time consumption of P-MQHOA algorithm and MOBi-DE algorithm (s)Func. f1 f2 f3 f4 f5 f6 f7 f8 f9 f10 f11 f12 f13 f14 P-MQHOA 0.01 0.01 0.01 0.01 0.01 0.01 0.02 0.1 0.1 0.1 1.76 0.28 4.32 12.14 MOBi-DE 1.56 1.46 2.82 2.64 1.54 2.04 1.68 18.34 15.28 41.24 46.28 36.49 51.74 72.36 可见采取固定分辨率策略的P-MQHOA的一维函数全局寻优时间消耗比MOBi-DE算法低两个数量级, 采取多级分辨率策略的P-MQHOA的多维函数全局寻优时间消耗比MOBi-DE算法也低很多.这是因为P-MQHOA局部收敛速度要比MOBi-DE等大部分智能算法要快的多.例如P-MQHOA对 $f_{14}$ 全局寻优过程中, 能够在13秒内找到约30 000个极值点, 然后筛选出全局极值点.
通过以上实验对比进行分析, 归纳如下:
GA采用二进制编码, 需要进行交叉变异选择遗传操作, 比较繁杂, 并且以前的知识随着种群的改变被改变, 没有记忆性, 有较多敏感的控制参数, 收敛速度慢; PSO一般采用实数编码, 粒子通过内部速度变化进行更新, 没有GA复杂的遗传操作, 此外通过当前搜索得到最优点进行信息共享, 较优解被所有粒子保存, 具有记忆能力, 收敛速度相对GA较快但很容易陷入局部最优解, 算法不稳定; DE采用实数编码、基于差分的简单变异操作和一对一的竞争生存策略, 降低了遗传操作的复杂性, 其独特的记忆能力使其可以动态跟踪当前的搜索情况以调整其搜索策略, 控制参数较少, 收敛速度相对于GA与PSO都要快. P-MQHOA无需编码, 拥有固定分辨率与多级分辨率两种策略, 每个区域的搜索相互独立, 每个采样中心点具有对所在区域极值的记忆能力, 操作简单, 控制参数较少, 收敛速度最快.
5. 结论
为解决多峰优化问题, 本文改进了MQHOA, 提出P-MQHOA算法, 把MQHOA的全局优化特性分散在多个均匀划分区域, 可以快速在各个区域找到极值点, 从而解决多峰优化问题.当实际问题除了最优解还需要备择解时, 该算法为全极值点优化中提供次优解提供了条件.此外针对具体问题本文提出了固定分辨率和多级分辨率两种策略, 前者有很好的全极值搜索能力, 后者有较强的实时性, 两种策略都能用于解决全局多峰优化问题.后面的研究工作将集中在该算法在多级分辨率策略下的全极值寻优性能研究, 提高其处理高维问题的能力, 扩展该算法在工业上的应用.
-
表 1 算法主要参数
Table 1 Main parameters of algorithms
说明 问题规模 变量数目:指最终OCP决策变量的数量, 包括所有的松弛变量.
约束数目:是指最终OCP的约束数目.
平均CPU时间:指每一次求解最终OCP所花费的平均时间.本质属性 可行域范围:假设OCP $^M$ 的所有决策变量为 $\Theta^M(x_k)$ , 其中 $x_k$ 是算法 $M$ 下的初始状态, OCP $^M$ 的可行域可定义为: ${F}^{M}=\{x| \exists \Theta^M(x)$ 使得OCP $^M$ 有可行解 表 2 Min-max算法的问题规模
Table 2 Problem scale of algorithms
算法 变量数目 约束数目 OL2M $Nm+Nr+1 $ $1LMI+Nn_{hx}+Nn_{hu}+n_S $ OL2M* $Nm+Nr+1 $ $ (q^{ N}-1)(1L+Nn_{hx}+n_S)+Nn_{hu} $ FF2M $Nm+Nr+1 $ $1LMI+N n_{hx}+Nn_{hu}+(2Nr+1)n_S $ FF2M* $Nm+Nr+1 $ $ (q^{N}-1)(N n_{hx}+Nn_{hu}+(2Nr+1)n_S+1) $ DF2M $N(N-1)/2+N m+N^2(n_{hx}+n_{hu})r+N r+1$ $1LMI+N n_{hx}+2N^2(n_{hx}+n_{hu})r+N n_{hu}+(2Nr+1)n_S $ DF2M* $N(N-1)/2+N m+N^2(n_{hx}+n_{hu})r+N r+1$ $ (q^{ N}-1)(N n_{hx}+2N^2(n_{hx}+n_{hu})r+N n_{hu}+(2Nr+1)n_S+1) $ DME2M* $ (q^{ N}-1)m/(q-1)+1$ $q^{ N}(1+ Nn_{hx} +N n_{hu} +n_S)$ 表 3 PT RMPC的Tube和参数化
Table 3 Tubes and parameterization of PT RMPC
输入Tube 状态Tube 参数化 $\left.\begin{array}{*{20}ll} \mathfrak{U}_k=\{{U}_k, {U}_{k+1},\cdots,{U}_{k+N-1}\}\\ \text{其中}\\ {U}_{k+j}=\bigoplus_{i=0}^{j} {U}_{k+j}^{i}, \forall j\in \textbf{N}_{[0,N-1]}\\ {U}_{k+j}^{i=0}=\{u_{k+j}^{i=0} \in \textbf{R}^m\}, \forall j\in \textbf{N}_{[0,N-1]}\\ {U}_{k+j}^{i}=co\{u_{k+j}^{i,l} \in \textbf{R}^m,\forall l \in \textbf{N}_{[1,q]},\\ \forall i\in \textbf{N}_{[1,j]},\forall j\in \textbf{N}_{[1,N-1]}\} \end{array}\right.$ $\left.\begin{array}{*{20}ll} \mathfrak{X}_k=\{{X}_k, {X}_{k+1},\cdots,{X}_{k+N-1}\}\\ \text{其中}\\ {X}_{k+j}=\bigoplus_{i=0}^{j} {X}_{k+j}^{i}, \forall j\in \textbf{N}_{[0,N]}\\ {X}_{k+j}^{i=0}=\{x_{k+j}^{i=0} \in \textbf{R}^n, \forall j\in \textbf{N}_{[0,N]}\}\\ {X}_{k+j}^{i}=co\{x_{k+j}^{i,l} \in \textbf{R}^n,l \in \textbf{N}_{[1,q]},\\ \forall i\in \textbf{N}_{[1,j-1]}, \forall j\in \textbf{N}_{[2,N]}\}\\ {X}_{k+j}^{i=j}=G{W}=co\{x_{k+j}^{i=j,l}=G\widetilde{w}_l,\\ \forall l\in \textbf{N}_{[1,q]}, \forall j\in \textbf{N}_{[1,N]}\}\end{array}\right.$ $\left.\begin{array}{*{20}ll} x_{k+j}=\sum_{i=0}^jx_{k+j}^i\\ \text{其中,}\forall i\in \textbf{N}_{[1,j]}, \forall j\in \textbf{N}_{[1,N]}\\ x_{k+j}^i=\sum_{l=1}^{q}\lambda^l_{k+j}x_{k+j}^{i,l}\\ \text{ 且 } \sum_{l=1}^{q}\lambda^l_{k+j}=1\\ u_{k+j}=\sum_{i=0}^ju_{k+j}^i\\ \text{其中 }~\forall i\in \textbf{N}_{[1,j]}, \forall j\in \textbf{N}_{[1,N-1]}\\ u_{k+j}^i=\sum_{l=1}^{q}\lambda^l_{k+j}u_{k+j}^{i,l},\\ \lambda^l_{k+j} \text{ 同 }~ x_{k+j}^i \end{array}\right.$ 表 4 基于Tube的RMPC算法的问题规模
Table 4 Problem scale tube-based RMPC
算法 变量数目 约束数目 RT $Nm+n$ $Nn_{hx}+Nn_{hu}+n_f+n_{sr}$ HT $Nm+ n +1$ $N(n_{hx}+n_{hu}+1)+n_f+n_{sr} $ PT $ (N-1)^2/2(mq+n_{hx}+n_{hu} )+ $ $N(N+1)/2+3N+Nn_s+N(N-1)(3(1+q)q/2+2)/2+ $ $ N(n+n_{h}+1)+N^2$ $n_{hx}+n_{hu}+N(1+q)qn_s/2+n_s$ 表 5 RMPC算法对比
Table 5 Comparison of RMPC algorithms
算法 变量数目 约束数目 平均CPU时间(s) MPC $5 $ $30 $ 0.2884 OL2M $16 $ $1LMI+35$ 0.3884 OL2M* $16 $ $25585 $ 49.4144 FF2M $16 $ $1LMI+135$ 1.1977 FF2M* $16 $ $139128 $ 267.5556 DF2M $326$ $1LMI+735 $ 1.9504 DF2M* $326$ $752928 $ 1444.8 DME2M* $ 342$ 36864 88.9460 RT $ 7$ $81$ 0.3247 HT 8 $86$ 0.3329 PT 140 636 0.7974 表 6 主要SMPC算法
Table 6 Main SMPC algorithms
表 7 RMPC和SMPC应用
Table 7 Applications of RMPC and SMPC algorithms
随机模型预测控制 鲁棒模型预测控制 无人汽车 驾驶转向控制[75] (DE)车辆导航[76] (ST)变道辅助[77] (SG)
巡航控制[78] (SG)车道保持与避障[79] (ST) [80] (DE)
自动驾驶控制[81] (SG)驾驶员行为建模[82] (DE)轨迹跟踪[83] (TB)传动动力系统控制[84] (MM)
车道保持与避障[85] (TB)智能家居 房屋气候控制[86] (SG)房屋遮阳镜片控制[87] (DE)
房屋能量控制[88] (SG)暖通空调系统建模与能量控制[89] (SG)温度控制[90] (MM)变风量空调系统控制[91] (MM) 电子电路 网络直流电机控制系统[92] (DE) 永磁同步电机驱动器[93] (TB)电子电路[94] (MM)热力电路[95] (MM) 飞行器 能源管理系统[96] (DE) 起落架系统振动抑制[97] (MM)
无人直升机轨迹跟踪控制[98] (MM)无人机飞行控制[99] (MM)
高超声速飞行器飞行控制[100] (MM)机器人 机器人导航与避障[101] (DE) 轨迹跟踪[102] (TB) 医疗 常压等离子射流控制[74] (ST)急救车辆调度[103] (SG) 药房库存管理[104](DE) 静脉麻醉控制[105] (MM) 过程工业 过程控制(四罐过程)[106] (DE) 浮选过程控制[107] (MM)热交换网控制[108] (MM)
连续搅拌反应釜[109] (MM)燃煤电站锅炉燃烧系统[110] (MM)
蒸馏塔控制[111] (MM)热轧带钢自动厚度控制[112] (MM)电网 微电网能量管理系统[113] (DE)微电网操控[114] (SG)
能源储存与生产的优化调度[115] (SG)电力调度[116] (SG)能源局域网优化调度[117] (SG)可再生能源微电网控制[118] (MM) 风力发电 风力涡轮机机械疲劳抑制[119] (DE)电池储能系统控制[120] (DE)风力发电系统控制[121] (DE) 风力涡轮机阻尼控制[122] (MM)风力涡轮机控制[123] (MM) 水资源 饮用水网络[124] (DE)水资源管理与利用[125] (DE) 城市交通 铁路货运车规模与分配控制[126](DE)能源管理[127] (SG) 城市道路交通网控制[128] (MM)地铁列车调度[129] (MM) 金融 动态套期保值和期权定价[130] (SG)欧式期权动态对冲[131] (SG) 其他 多层供应链管理[132] (SG)智能热网中的能量平衡[133](SG)云副本放置技术[134] (DE) 内燃机的热量管理[135] (MM)基于图像的视觉伺服控制[136] (MM) -
[1] Kassmann D E, Badgwell T A, Hawkins R B. Robust steady-state target calculation for model predictive control. AIChE Journal, 2000, 46(5): 1007-1024 doi: 10.1002/(ISSN)1547-5905 [2] Qin S J, Badgwell T A. A survey of industrial model predictive control technology. Control Engineering Practice, 2003, 11(7): 733-764 doi: 10.1016/S0967-0661(02)00186-7 [3] Hrovat D, Di Cairano S, Tseng H E, Kolmanovsky I V. The development of model predictive control in automotive industry: a survey. In: Proceedings of the 2012 IEEE International Conference on Control Applications (CCA). Dubrovnik, Croatia: IEEE, 2012. 295-302 [4] 席裕庚, 李德伟, 林姝.模型预测控制--现状与挑战.自动化学报, 2013, 39(3): 222-236 http://www.aas.net.cn/CN/abstract/abstract17874.shtmlXi Yu-Geng, Li De-Wei, Lin Shu. Model predictive control-status and challenges. Acta Automatica Sinica, 2013, 39(3): 222-236 http://www.aas.net.cn/CN/abstract/abstract17874.shtml [5] Mesbah A. Stochastic model predictive control: an overview and perspectives for future research. IEEE Control Systems, 2016, 36{6: 30-44 doi: 10.1109/MCS.2016.2602087 [6] Wang Y J, Rawlings J B. A new robust model predictive control method I: theory and computation. Journal of Process Control, 2004, 14(3): 231-247 doi: 10.1016/S0959-1524(02)00133-6 [7] Mayne D Q. Model predictive control: recent developments and future promise. Automatica, 2014, 50(12): 2967-2986 doi: 10.1016/j.automatica.2014.10.128 [8] Campo P J, Morari M. Robust model predictive control. In: Proceedings of the 1987 American Control Conference. Minneapolis, MN, USA: IEEE, 1987. 1021-1026 [9] Allwright J C, Papavasiliou G C. On linear programming and robust model-predictive control using impulse-responses. Systems & Control Letters, 1992, 18(2): 159-164 https://www.researchgate.net/publication/222725809_On_Linear_Programming_and_Robust_Model-Predictive_Control_Using_Impulse-Responses [10] Zheng Z Q, Morari M. Robust stability of constrained model predictive control. In: Proceedings of the 1993 American Control Conference. San Francisco, California, USA: IEEE, 1993. 379-383 [11] Lee J H, Yu Z H. Worst-case formulations of model predictive control for systems with bounded parameters. Automatica, 1997, 33(5): 763-781 doi: 10.1016/S0005-1098(96)00255-5 [12] Scokaert P O M, Mayne D Q. Min-max feedback model predictive control for constrained linear systems. IEEE Transactions on Automatic Control, 1998, 43(8): 1136-1142 doi: 10.1109/9.704989 [13] Blanchini F. Control synthesis for discrete time systems with control and state bounds in the presence of disturbances. Journal of Optimization Theory and Applications, 1990, 65(1): 29-40 doi: 10.1007/BF00941157 [14] Gossner J R, Kouvaritakis B, Rossiter J A. Stable generalized predictive control with constraints and bounded disturbances. Automatica, 1997, 33(4): 551-568 doi: 10.1016/S0005-1098(96)00214-2 [15] Lee Y I, Kouvaritakis B. Constrained receding horizon predictive control for systems with disturbances. International Journal of Control, 1999, 72(11): 1027-1032 doi: 10.1080/002071799220579 [16] Raković S V, Kerrigan E C, Kouramas K I, Mayne D Q. Invariant Approximations of Robustly Positively Invariant Sets for Constrained Linear Discrete-Time Systems Subject to Bounded Disturbances, Technical Report CUED/F-INFENG/TR.473, Department of Engineering, University of Cambridge, Cambridge, UK, 2004. [17] Limon D, Alvarado I, Alamo T, Camacho E F. Robust tube-based MPC for tracking of constrained linear systems with additive disturbances. Journal of Process Control, 2010, 20(3): 248-260 doi: 10.1016/j.jprocont.2009.11.007 [18] Cannon M, Kouvaritakis B, Rakovic S V, Cheng Q F. Stochastic tubes in model predictive control with probabilistic constraints. IEEE Transactions on Automatic Control, 2011, 56(1): 194-200 doi: 10.1109/TAC.2010.2086553 [19] Raković S V, Kouvaritakis B, Findeisen R, Cannon M. Homothetic tube model predictive control. Automatica, 2012, 48(8): 1631-1638 doi: 10.1016/j.automatica.2012.05.003 [20] Raković S V, Kouvaritakis B, Cannon M, Panos C, Findeisen R. Parameterized tube model predictive control. IEEE Transactions on Automatic Control, 2012, 57(11): 2746-2761 doi: 10.1109/TAC.2012.2191174 [21] Muñoz-Carpintero D, Kouvaritakis B, Cannon M. Striped parameterized tube model predictive control. In: Proceedings of the 19th IFAC World Congress. Cape Town, South Africa: IFAC, 2014. 11998-12003 [22] Farina M, Giulioni L, Magni L, Scattolini R. A probabilistic approach to model predictive control. In: Proceedings of the 52nd Annual Conference on Decision and Control (CDC). Florence, Italy: IEEE, 2013. 7734-7739 [23] Farina M, Giulioni L, Magni L, Scattolini R. An approach to output-feedback MPC of stochastic linear discrete-time systems. Automatica, 2015, 55: 140-149 doi: 10.1016/j.automatica.2015.02.039 [24] Hokayem P, Chatterjee D, Lygeros J. On stochastic receding horizon control with bounded control inputs. In: Proceedings of the 48th IEEE Conference on Decision and Control, Held Jointly with the 28th Chinese Control Conference. Shanghai, China: IEEE, 2009. 6359-6364 [25] Hokayem P, Chatterjee D, Ramponi F, Chaloulos G, Lygeros J. Stable stochastic receding horizon control of linear systems with bounded control inputs. In: Proceedings of the 19th International Symposium on Mathematical Theory of Networks and Systems. Budapest, Hungary: MTNS, 2010. 31-36 [26] Primbs J A, Sung C H. Stochastic receding horizon control of constrained linear systems with state and control multiplicative noise. IEEE Transactions on Automatic Control, 2009, 54(2): 221-230 doi: 10.1109/TAC.2008.2010886 [27] Hokayem P, Cinquemani E, Chatterjee D, Ramponi F, Lygeros J. Stochastic receding horizon control with output feedback and bounded controls. Automatica, 2012, 48(1): 77-88 doi: 10.1016/j.automatica.2011.09.048 [28] Lucia S, Zometa P, Kögel M, Findeisen R. Efficient stochastic model predictive control based on polynomial chaos expansions for embedded applications. In: Proceedings of the 54th Annual Conference on Decision and Control (CDC). Osaka, Japan: IEEE, 2015. 3006-3012 [29] Kouvaritakis B, Cannon M, Raković S V, Cheng Q F. Explicit use of probabilistic distributions in linear predictive control. Automatica, 2010, 46(10): 1719-1724 doi: 10.1016/j.automatica.2010.06.034 [30] Kothare M V, Balakrishnan V, Morari M. Robust constrained model predictive control using linear matrix inequalities. Automatica, 1996, 32(10): 1361-1379 doi: 10.1016/0005-1098(96)00063-5 [31] Kouvaritakis B, Cannon M. Model Predictive Control: Classical, Robust and Stochastic. Switzerland: Springer, 2015. [32] Mayne D. Robust and stochastic model predictive control: are we going in the right direction? Annual Reviews in Control, 2016, 41: 184-192 doi: 10.1016/j.arcontrol.2016.04.006 [33] Bemporad A, Morari M. Robust model predictive control: a survey, Robustness in Identification and Control. London: Springer, 1999. 207-226 [34] Löfberg J. Approximations of closed-loop minimax MPC. In: Proceedings of the 42nd IEEE Conference on Decision and Control. Maui, HI, USA: IEEE, 2003, 2: 1438-1442 [35] Alamo T, Ramirez D R, de La Peña D M, Camacho E F. Min-max MPC using a tractable QP problem. Automatica, 2007, 43(4): 693-700 doi: 10.1016/j.automatica.2006.10.006 [36] Gao Y, Chong K T. The explicit constrained min-max model predictive control of a discrete-time linear system with uncertain disturbances. IEEE Transactions on Automatic Control, 2012, 57(9): 2373-2378 doi: 10.1109/TAC.2012.2186090 [37] de la Peña D M, Ramírez D R, Camacho E F, Alamo T. Explicit solution of min-max MPC with additive uncertainties and quadratic criterion. Systems & control letters, 2006, 55(4): 266-274 https://www.researchgate.net/publication/257012284_Explicit_solution_of_min-max_MPC_with_additive_uncertainties_and_quadratic_criterion [38] Bertsekas D P. Nonlinear Programming (Second edition). Belmont, Massachusetts: Athena Scientific, 1999 [39] Löfberg J. Minimax Approaches to Robust Model Predictive Control [Ph.D. dissertation], Linköping University, Sweden, 2003. [40] 王振英. 不确定扰动情形下离散系统约束鲁棒模型预测控制[硕士学位论文], 哈尔滨理工大学, 中国, 2005.Wang Zhen-Ying. Constrained Robust Model Predictive Control of Discrete with Uncertaint Disturbance Case [Master dissertation], Harbin University of Science and Technology, China, 2005. [41] Bemporad A. Reducing conservativeness in predictive control of constrained systems with disturbances. In: Proceedings of the 37th IEEE Conference on Decision and Control. Tampa, FL, USA: IEEE, 1998, 2: 1384-1389 [42] Goulart P J, Kerrigan E C, Maciejowski J M. Optimization over state feedback policies for robust control with constraints. Automatica, 2006, 42(4): 523-533 doi: 10.1016/j.automatica.2005.08.023 [43] Langson W, Chryssochoos I, Raković S V, Mayne D Q. Robust model predictive control using tubes. Automatica, 2004, 40(1): 125-133 doi: 10.1016/j.automatica.2003.08.009 [44] Gautam A, Chu Y C, Soh Y C. Optimized dynamic policy for receding horizon control of linear time-varying systems with bounded disturbances. IEEE Transactions on Automatic Control, 2012, 57(4): 973-988 doi: 10.1109/TAC.2011.2170109 [45] Muñoz-Carpintero D, Cannon M, Kouvaritakis B. Robust MPC strategy with optimized polytopic dynamics for linear systems with additive and multiplicative uncertainty. Systems & Control Letters, 2015, 81: 34-41 https://www.researchgate.net/publication/271549753_Recursively_feasible_Robust_MPC_for_linear_systems_with_additive_and_multiplicative_uncertainty_using_optimized_polytopic_dynamics [46] Blanchini F. Survey paper: set invariance in control. Automatica, 1999, 35(11): 1747-1767 doi: 10.1016/S0005-1098(99)00113-2 [47] Herceg M, Kvasnica M, Jones C N, Morari M. Multi-parametric toolbox 3.0. In: Proceedings of the 2013 European Control Conference (ECC). Zurich, Switzerland: IEEE, 2013. 502-510 [48] Mayne D Q, Rawlings J B, Rao C V, Scokaert P O M. Constrained model predictive control: stability and optimality. Automatica, 2000, 36(6): 789-814 doi: 10.1016/S0005-1098(99)00214-9 [49] Goodwin G C, Kong H, Mirzaeva G, Seron M M. Robust model predictive control: reflections and opportunities. Journal of Control and Decision, 2014, 1(2): 115-148 doi: 10.1080/23307706.2014.913837 [50] Zheng P Y, Li D W, Xi Y G, Zhang J. Improved model prediction and RMPC design for LPV systems with bounded parameter changes. Automatica, 2013, 49(12): 3695-3699 doi: 10.1016/j.automatica.2013.09.024 [51] 平续斌, 丁宝苍, 韩崇昭.动态输出反馈鲁棒模型预测控制.自动化学报, 2012, 38(1): 31-37 http://www.aas.net.cn/CN/abstract/abstract17542.shtmlPing Xu-Bin, Ding Bao-Cang, Han Chong-Zhao. Dynamic output feedback robust model predictive control. Acta Automatica Sinica, 2012, 38(1): 31-37 http://www.aas.net.cn/CN/abstract/abstract17542.shtml [52] Cuzzola F A, Geromel J C, Morari M. An improved approach for constrained robust model predictive control. Automatica, 2002, 38(7): 1183-1189 doi: 10.1016/S0005-1098(02)00012-2 [53] Kouvaritakis B, Rossiter J A, Schuurmans J. Efficient robust predictive control. IEEE Transactions on Automatic Control, 2000, 45(8): 1545-1549 doi: 10.1109/9.871769 [54] Mayne D Q, Langson W. Robustifying model predictive control of constrained linear systems. Electronics Letters, 2001, 37(23): 1422-1423 doi: 10.1049/el:20010951 [55] Mayne D Q, Seron M M, Raković S V. Robust model predictive control of constrained linear systems with bounded disturbances. Automatica, 2005, 41(2): 219-224 doi: 10.1016/j.automatica.2004.08.019 [56] Raković S V, Kouvaritakis B, Findeisen R, Cannon M. Simple homothetic tube model predictive control. In: Proceedings of the 19th International Symposium on Mathematical Theory of Networks and Systems-MTNS. Budapest, Hungary: MTNS, 2010. 1411-1418 [57] Kolmanovsky I, Gilbert E G. Theory and computation of disturbance invariant sets for discrete-time linear systems. Mathematical Problems in Engineering, 1998, 4(4): 317-367 doi: 10.1155/S1024123X98000866 [58] Calafiore G, Campi M C. Uncertain convex programs: randomized solutions and confidence levels. Mathematical Programming, 2005, 102(1): 25-46 doi: 10.1007/s10107-003-0499-y [59] Campi M C, Garatti S, Prandini M. The scenario approach for systems and control design. Annual Reviews in Control, 2009, 33(2): 149-157 doi: 10.1016/j.arcontrol.2009.07.001 [60] Schildbach G, Fagiano L, Frei C, Morari M. The scenario approach for stochastic model predictive control with bounds on closed-loop constraint violations. Automatica, 2014, 50(12): 3009-3018 doi: 10.1016/j.automatica.2014.10.035 [61] Cannon M, Cheng Q F, Kouvaritakis B, Raković S V. Stochastic tube MPC with state estimation. Automatica, 2012, 48(3): 536-541 doi: 10.1016/j.automatica.2011.08.058 [62] Schwarm A T, Nikolaou M. Chance-constrained model predictive control. AIChE Journal, 1999, 45(8): 1743-1752 doi: 10.1002/(ISSN)1547-5905 [63] Magni L, Pala D, Scattolini R. Stochastic model predictive control of constrained linear systems with additive uncertainty. In: Proceedings of the 2009 European Control Conference (ECC). Budapest, Hungary: IEEE, 2009. 2235-2240 [64] Farina M, Scattolini R. Model predictive control of linear systems with multiplicative unbounded uncertainty and chance constraints. Automatica, 2016, 70: 258-265 doi: 10.1016/j.automatica.2016.04.008 [65] Farina M, Giulioni L, Scattolini R. Stochastic linear model predictive control with chance constraints——a review. Journal of Process Control, 2016, 44: 53-67 doi: 10.1016/j.jprocont.2016.03.005 [66] Batina I. Model Predictive Control for Stochastic Systems By Randomized Algorithms [Ph.D. dissertation], Technische Universiteit Eindhoven, Netherlands, 2004. [67] Blackmore L, Ono M, Bektassov A, Williams B C. A probabilistic particle-control approximation of chance-constrained stochastic predictive control. IEEE Transactions on Robotics, 2010, 26(3): 502-517 doi: 10.1109/TRO.2010.2044948 [68] Campi M C, Garatti S. The exact feasibility of randomized solutions of uncertain convex programs. SIAM Journal on Optimization, 2008, 19(3): 1211-1230 doi: 10.1137/07069821X [69] Campi M C, Garatti S. A sampling-and-discarding approach to chance-constrained optimization: feasibility and optimality. Journal of Optimization Theory and Applications, 2011, 148(2): 257-280 doi: 10.1007/s10957-010-9754-6 [70] Calafiore G C, Campi M C. The scenario approach to robust control design. IEEE Transactions on Automatic Control, 2006, 51(5): 742-753 doi: 10.1109/TAC.2006.875041 [71] Calafiore G C. Random convex programs. SIAM Journal on Optimization, 2010, 20(6): 3427-3464 doi: 10.1137/090773490 [72] Alamo T, Tempo R, Luque A, Ramirez D R. Randomized methods for design of uncertain systems: sample complexity and sequential algorithms. Automatica, 2015, 52: 160-172 doi: 10.1016/j.automatica.2014.11.004 [73] Cannon M, Kouvaritakis B, Wu X J. Probabilistic constrained MPC for multiplicative and additive stochastic uncertainty. IEEE Transactions on Automatic Control, 2009, 54(7): 1626-1632 doi: 10.1109/TAC.2009.2017970 [74] Gidon D, Graves D B, Mesbah A. Model predictive control of thermal effects of an atmospheric pressure plasma jet for biomedical applications. In: Proceedings of the 2016 American Control Conference (ACC). Boston, MA, USA: IEEE, 2016. 4889-4894 [75] 王秋, 曲婷, 陈虹.基于随机模型预测控制的自主车辆转向控制.信息与控制, 2015, 44(4): 499-506 http://www.cnki.com.cn/Article/CJFDTOTAL-XXYK201504021.htmWang Qiu, Qu Ting, Chen Hong. Steering control of autonomous vehicles based on stochastic model predictive control. Information and Control, 2015, 44(4): 499-506 http://www.cnki.com.cn/Article/CJFDTOTAL-XXYK201504021.htm [76] Carvalho A, Gao Y Q, Lefevre S, Borrelli F. Stochastic predictive control of autonomous vehicles in uncertain environments. In: Proceedings of the 12th International Symposium on Advanced Vehicle Control. Tuat, Japan: AVEC, 2014. [77] Schildbach G, Borrelli F. Scenario model predictive control for lane change assistance on highways. In: Proceedings of the 2015 IEEE Intelligent Vehicles Symposium (Ⅳ). Seoul, Korea: IEEE, 2015. 611-616 [78] Moser D, Waschl H, Kirchsteiger H, Schmied R, del Re L. Cooperative adaptive cruise control applying stochastic linear model predictive control strategies. In: Proceedings of the 2015 European Control Conference (ECC). Linz, Austria: IEEE, 2015. 3383-3388 [79] Liu C C, Carvalho A, Schildbach G, Hedrick J K. Stochastic predictive control for lane keeping assistance systems using a linear time-varying model. In: Proceedings of the 2015 American Control Conference (ACC). Chicago, IL, USA: ACC, 2015. 3355-3360 [80] Gray A, Gao Y Q, Lin T, Hedrick J K, Borrelli F. Stochastic predictive control for semi-autonomous vehicles with an uncertain driver model. In: Proceedings of the 16th International IEEE Annual Conference on Intelligent Transportation Systems-(ITSC). The Hague, the Netherlands: IEEE, 2013. 2329-2334 [81] Carvalho A, Lefévre S, Schildbach G, Kong J, Borrelli F. Automated driving: the role of forecasts and uncertainty-a control perspective. European Journal of Control, 2015, 24: 14-32 doi: 10.1016/j.ejcon.2015.04.007 [82] 曲婷. 基于随机模型预测控制的驾驶员行为建模[博士学位论文], 吉林大学, 中国, 2015.Qu Ting. Driver Behavior Modeling Based on Stochastic Model Predictive Control [Ph.D. dissertation], Jilin University, China, 2015. [83] Kayacan E, Ramon H, Saeys W. Robust trajectory tracking error model-based predictive control for unmanned ground vehicles. IEEE/ASME Transactions on Mechatronics, 2016, 21(2): 806-814 doi: 10.1109/TMECH.2015.2492984 [84] Zhu X Y, Zhang H, Cao D P, Fang Z D. Robust control of integrated motor-transmission powertrain system over controller area network for automotive applications. Mechanical Systems and Signal Processing, 2015, 58-59: 15-28 http://d.scholar.cnki.net/detail/SJESTEMP_U/SJES14123100018385 [85] Gao Y Q, Gray A, Tseng H E, Borrelli F. A tube-based robust nonlinear predictive control approach to semiautonomous ground vehicles. Vehicle System Dynamics, 2014, 52(6): 802-823 doi: 10.1080/00423114.2014.902537 [86] Zhang X J, Grammatico S, Schildbach G, Goulart P, Lygeros J. On the sample size of randomized MPC for chance-constrained systems with application to building climate control. In: Proceedings of the 2014 European Control Conference (ECC). Strasbourg, France: IEEE, 2014. 478-483 [87] Yao J. Determining the energy performance of manually controlled solar shades: a stochastic model based co-simulation analysis. Applied Energy, 2014, 127: 64-80 doi: 10.1016/j.apenergy.2014.04.046 [88] Tanner R A, Henze G P. Stochastic control optimization for a mixed mode building considering occupant window opening behaviour. Journal of Building Performance Simulation, 2014, 7(6): 427-444 doi: 10.1080/19401493.2013.863384 [89] Ma Y D, Matuško J, Borrelli F. Stochastic model predictive control for building HVAC systems: complexity and conservatism. IEEE Transactions on Control Systems Technology, 2015, 23(1): 101-116 doi: 10.1109/TCST.2014.2313736 [90] Xu X H, Wang S W, Huang G S. Robust MPC for temperature control of air-conditioning systems concerning on constraints and multitype uncertainties. Building Services Engineering Research and Technology, 2010, 31(1): 39-55 doi: 10.1177/0143624409352420 [91] 杨世忠, 任庆昌.鲁棒模型预测控制在变风量空调系统中的应用.计算机工程与应用, 2013, 49(6): 8-11, 104 http://www.cnki.com.cn/Article/CJFDTOTAL-JSGG201306004.htmYang Shi-Zhong, Ren Qing-Chang. Robust model predictive control for variable air volume air conditioning system. Computer Engineering and Applications, 2013, 49(6): 8-11, 104 http://www.cnki.com.cn/Article/CJFDTOTAL-JSGG201306004.htm [92] Ren M F, Zhang J H, Jiang M, Yu M, Xu J L. Minimum (h, Φ) entropy control for non-Gaussian stochastic networked control systems and its application to a networked dc motor control system. IEEE Transactions on Control Systems Technology, 2015, 23(1): 406-411 doi: 10.1109/TCST.2014.2324978 [93] Preindl M. Robust control invariant sets and Lyapunov-based mpc for ipm synchronous motor drives. IEEE Transactions on Industrial Electronics, 2016, 63(6): 3925-3933 doi: 10.1109/TIE.2016.2527722 [94] Zhao X D, Zhang L X, Shi P, Karimi H R. Robust control of continuous-time systems with state-dependent uncertainties and its application to electronic circuits. IEEE Transactions on Industrial Electronics, 2014, 61(8): 4161-4170 doi: 10.1109/TIE.2013.2286568 [95] Rhouma A, Bouani F. Robust model predictive control of uncertain fractional systems: a thermal application. IET Control Theory & Applications, 2014, 8(17): 1986-1994 https://www.researchgate.net/publication/270579793_Robust_model_predictive_control_of_uncertain_fractional_systems_A_thermal_application [96] Shahsavari B, Maasoumy M, Sangiovanni-Vincentelli A, Horowitz R. Stochastic model predictive control design for load management system of aircraft electrical power distribution. In: Proceedings of the 2015 American Control Conference (ACC). Chicago, IL, USA: IEEE, 2015. 3649-3655 [97] Hajiloo A, Xie W F. The stochastic robust model predictive control of shimmy vibration in aircraft landing gears. Asian Journal of Control, 2015, 17(2): 476-485 doi: 10.1002/asjc.v17.2 [98] Marantos P, Bechlioulis C P, Kyriakopoulos K J. Robust trajectory tracking control for small-scale unmanned helicopters with model uncertainties. IEEE Transactions on Control Systems Technology, 2017, doi: 10.1109/TCST.2016.2642160 [99] Alexis K, Papachristos C, Siegwart R, Tzes A. Robust model predictive flight control of unmanned rotorcrafts. Journal of Intelligent & Robotic Systems, 2016, 81(3-4): 443-469 https://www.researchgate.net/publication/279232323_Robust_Model_Predictive_Flight_Control_of_Unmanned_Rotorcrafts [100] 张少康. 基于LPV模型的高超声速飞行器鲁棒预测控制方法研究[硕士学位论文], 天津大学, 中国, 2014.Zhang Shao-Kang. Research on Robust Model Predictive Control Based on LPV Model for Hypersonic Vehicles [Master dissertation], Tianjin University, China, 2014. [101] Blackmore L, Ono M, Williams B C. Chance-constrained optimal path planning with obstacles. IEEE Transactions on Robotics, 2011, 27(6): 1080-1094 doi: 10.1109/TRO.2011.2161160 [102] González R, Fiacchini M, Guzmán J L, Á lamo T, Rodríguez F. Robust tube-based predictive control for mobile robots in off-road conditions. Robotics and Autonomous Systems, 2011, 59(10): 711-726 doi: 10.1016/j.robot.2011.05.006 [103] Goodwin G C, Medioli A M. Scenario-based, closed-loop model predictive control with application to emergency vehicle scheduling. International Journal of Control, 2013, 86(8): 1338-1348 doi: 10.1080/00207179.2013.788215 [104] Jurado I, Maestre J M, Velarde P, Ocampo-Martinez C, Fern{andez I, Tejera B I, del Prado J R. Stock management in hospital pharmacy using chance-constrained model predictive control. Computers in Biology and Medicine, 2016, 72: 248-255 doi: 10.1016/j.compbiomed.2015.11.011 [105] Chang H, Krieger A, Astolfi A, Pistikopoulos E N. Robust multi-parametric model predictive control for LPV systems with application to anaesthesia. Journal of Process Control, 2014, 24(10): 1538-1547 doi: 10.1016/j.jprocont.2014.07.005 [106] Jurado I, Millán P, Quevedo D, Rubio F R. Stochastic MPC with applications to process control. International Journal of Control, 2015, 88(4): 792-800 doi: 10.1080/00207179.2014.975845 [107] Tang Y, Peng C, Yin S, Qiu J B, Gao H J, Kaynak O. Robust model predictive control under saturations and packet dropouts with application to networked flotation processes. IEEE Transactions on Automation Science and Engineering, 2014, 11(4): 1056-1064 doi: 10.1109/TASE.2013.2283304 [108] Bakošová M, Oravec J. Robust model predictive control for heat exchanger network. Applied Thermal Engineering, 2014, 73(1): 924-930 doi: 10.1016/j.applthermaleng.2014.08.023 [109] Ghaffari V, Naghavi S V, Safavi A A. Robust model predictive control of a class of uncertain nonlinear systems with application to typical CSTR problems. Journal of Process Control, 2013, 23(4): 493-499 doi: 10.1016/j.jprocont.2013.01.009 [110] 刘正峰, 王景成, 史元浩, 张浪文. 燃煤电站锅炉燃烧系统的鲁棒分布式模型预测控制. 见: 第25届中国过程控制会议. 大连, 中国: 中国自动化学会过程控制专业委员会, 2014. 6Liu Zheng-Feng, Wang Jing-Cheng, Shi Yuan-Hao, Zhang Lang-Wen. Robust distributed model prediction control for coal-fired power plant boiler combustion system. In: Proceedings of the 25th Chinese Process Control Conference. Dalian, China, 2014. 6 [111] Yorgashov Yu A. Robust model predictive control of industrial reactive distillation column. Journal of Shenyang University of Chemical Technology, 2009, 23(1): 92-96 http://en.cnki.com.cn/Article_en/CJFDTOTAL-SYHY200901022.htm [112] 杨斌虎, 杨卫东, 曲蕾.热轧带钢自动厚度控制鲁棒预测模型与控制策略.控制理论与应用, 2011, 28(6): 799-804 http://www.cnki.com.cn/Article/CJFDTOTAL-KZLY201106007.htmYang Bin-Hu, Yang Wei-Dong, Qu Lei. Robust predictive-model-control strategy for automatic gauge control in hot strip finishing milling. Control Theory & Applications, 2011, 28(6): 799-804 http://www.cnki.com.cn/Article/CJFDTOTAL-KZLY201106007.htm [113] Cominesi S R, Farina M, Giulioni L, Picasso B, Scattolini R. Two-layer predictive control of a micro-grid including stochastic energy sources. In: Proceedings of the 2015 American Control Conference (ACC). Chicago, IL, USA: IEEE, 2015. 918-923 [114] Hans C A, Sopasakis P, Bemporad A, Raisch J, Reincke-Collon C. Scenario-based model predictive operation control of islanded microgrids. In: Proceedings of the 54th Annual Conference on Decision and Control (CDC). Osaka, Japan: IEEE, 2015. 3272-3277 [115] Zhu D H, Hug G. Decomposed stochastic model predictive control for optimal dispatch of storage and generation. IEEE Transactions on Smart Grid, 2014, 5(4): 2044-2053 doi: 10.1109/TSG.2014.2321762 [116] Patrinos P, Trimboli S, Bemporad A. Stochastic MPC for real-time market-based optimal power dispatch. In: Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC). Orlando, FL, USA: IEEE, 2011. 7111-7116 [117] 张彦, 张涛, 刘亚杰, 郭波.基于随机模型预测控制的能源局域网优化调度研究.中国电机工程学报, 2016, 36(13): 3451-3462 http://www.cnki.com.cn/Article/CJFDTOTAL-ZGDC201613007.htmZhang Yan, Zhang Tao, Liu Ya-Jie, Guo Bo. Stochastic model predictive control for energy management optimization of an energy local network. Proceedings of the CSEE, 2016, 36(13): 3451-3462 http://www.cnki.com.cn/Article/CJFDTOTAL-ZGDC201613007.htm [118] Velarde P, Maestre J M, Ocampo-Martinez C, Bordons C. Application of robust model predictive control to a renewable hydrogen-based microgrid. In: Proceedings of the 15th European Control Conference (ECC). Aalborg, Denmark: ECC, 2016. 1209-1214 [119] Riverso S, Mancini S, Sarzo F, Ferrari-Trecate G. Model predictive controllers for reduction of mechanical fatigue in wind farms. IEEE Transactions on Control Systems Technology, 2017, 25(2): 535-549 doi: 10.1109/TCST.2016.2572170 [120] Kou P, Gao F, Guan X H. Stochastic predictive control of battery energy storage for wind farm dispatching: using probabilistic wind power forecasts. Renewable Energy, 2015, 80: 286-300 doi: 10.1016/j.renene.2015.02.001 [121] 胡婧. 随机模型预测控制在风力发电系统的应用[硕士学位论文], 华北电力大学, 中国, 2015.Hu Jing. Research on application of stochastic model predictive control in wind power generation systems [Master dissertation], North China Electric Power University, China, 2015. [122] Evans M A, Cannon M, Kouvaritakis B. Robust MPC tower damping for variable speed wind turbines. IEEE Transactions on Control Systems Technology, 2015, 23(1): 290-296 doi: 10.1109/TCST.2014.2310513 [123] Mirzaei M, Poulsen N K, Niemann H H. Robust model predictive control of a wind turbine. In: Proceedings of the 2012 American Control Conference (ACC). Montréal, Canada: IEEE, 2012. 4393-4398 [124] Grosso J M, Ocampo-Martínez C, Puig V, Joseph B. Chance-constrained model predictive control for drinking water networks. Journal of Process Control, 2014, 24(5): 504-516 doi: 10.1016/j.jprocont.2014.01.010 [125] 黄翀胤, 罗灿, 赵勇, 韩守东.随机预测控制模型及其在水资源可持续中的应用.水电能源科学, 2011, 29(12): 16-20 doi: 10.3969/j.issn.1000-7709.2011.12.005Huang Chong-Yin, Luo Can, Zhao Yong, Han Shou-Dong. Application of stochastic predictive control model to sustainable utilization of water resources. Water Resources and Power, 2011, 29(12): 16-20 doi: 10.3969/j.issn.1000-7709.2011.12.005 [126] Milenković M S, Bojović N J, Švadlenka L, Melichar V. A stochastic model predictive control to heterogeneous rail freight car fleet sizing problem. Transportation Research Part E: Logistics and Transportation Review, 2015, 82: 162-198 doi: 10.1016/j.tre.2015.07.009 [127] Di Cairano S, Bernardini D, Bemporad A, Kolmanovsky I V. Stochastic MPC with learning for driver-predictive vehicle control and its application to HEV energy management. IEEE Transactions on Control Systems Technology, 2014, 22(3): 1018-1031 doi: 10.1109/TCST.2013.2272179 [128] Tettamanti T, Luspay T, Kulcsár B, Péni T, Varga I. Robust control for urban road traffic networks. IEEE Transactions on Intelligent Transportation Systems, 2014, 15(1): 385-398 doi: 10.1109/TITS.2013.2281666 [129] Li S K, De Schutter B, Yang L X, Gao Z Y. Robust model predictive control for train regulation in underground railway transportation. IEEE Transactions on Control Systems Technology, 2016, 24(3): 1075-1083 doi: 10.1109/TCST.2015.2480839 [130] Bemporad A, Bellucci L, Gabbriellini T. Dynamic option hedging via stochastic model predictive control based on scenario simulation. Quantitative Finance, 2014, 14(10): 1739-1751 doi: 10.1080/14697688.2011.649780 [131] 张卫国, 杜谦.基于随机模型预测控制的欧式期权动态对冲研究.华南理工大学学报(社会科学版), 2016, 18(4): 1-9 http://www.cnki.com.cn/Article/CJFDTOTAL-HNLS201604002.htmZhang Wei-Guo, Du Qian. Dynamic hedging european option based on stochastic model predictive controldynamic hedging European option based on stochastic model predictive control. Journal of South China University of Technology (Social Science Edition), 2016, 18(4): 1-9 http://www.cnki.com.cn/Article/CJFDTOTAL-HNLS201604002.htm [132] Schildbach G, Morari M. Scenario-based model predictive control for multi-echelon supply chain management. European Journal of Operational Research, 2016, 252(2): 540-549 doi: 10.1016/j.ejor.2016.01.051 [133] Rostampour V, Keviczky T. Robust randomized model predictive control for energy balance in smart thermal grids. In: Proceedings of the 2016 European Control Conference (ECC). Aalborg, Denmark: IEEE, 2016. 1201-1208 [134] Ghanbari H, Litoiu M, Pawluk P, Barna C. Replica placement in cloud through simple stochastic model predictive control. In: Proceedings of the 2014 IEEE 7th International Conference on Cloud Computing (CLOUD). Anchorage, AK, USA: IEEE, 2014. 80-87 [135] Pizzonia F, Castiglione T, Bova S. A Robust Model Predictive Control for efficient thermal management of internal combustion engines. Applied Energy, 2016, 169: 555-566 doi: 10.1016/j.apenergy.2016.02.063 [136] Hajiloo A, Keshmiri M, Xie W F, Wang T T. Robust online model predictive control for a constrained image-based visual servoing. IEEE Transactions on Industrial Electronics, 2016, 63(4): 2242-2250 https://www.researchgate.net/profile/Amir_Hajiloo/publication/287965502_Robust_On-Line_Model_Predictive_Control_for_a_Constrained_Image_Based_Visual_Servoing/links/56e1926608aec09a8bc08154.pdf [137] Carrasco D S, Goodwin G C. Feedforward model predictive control. Annual Reviews in Control, 2011, 35(2): 199-206 doi: 10.1016/j.arcontrol.2011.10.007 [138] Paulson J A, Mesbah A, Streif S, Findeisen R, Braatz R D. Fast stochastic model predictive control of high-dimensional systems. In: Proceedings of the 53rd Conference on Decision and Control. Los Angeles, California, USA: IEEE, 2014. 2802-2809 [139] Terejanu G, Singla P, Singh T, Scott P D. Approximate interval method for epistemic uncertainty propagation using polynomial chaos and evidence theory. In: Proceedings of the 2010 American Control Conference. Baltimore, MD, USA: IEEE, 2010. 349-354 [140] Bayer F A, Müller M A, Allgöwer F. Tube-based robust economic model predictive control. Journal of Process Control, 2014, 24(8): 1237-1246 doi: 10.1016/j.jprocont.2014.06.006 [141] Hovgaard T G, Larsen L F S, Jorgensen J B. Robust economic MPC for a power management scenario with uncertainties. In: Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC). Orlando, FL, USA: IEEE, 2011. 1515-1520 [142] Bayer F A, Lorenzen M, Müller M A, Allgöwer F. Robust economic model predictive control using stochastic information. Automatica, 2016, 74: 151-161 doi: 10.1016/j.automatica.2016.08.008 [143] Feller C, Ebenbauer C. A stabilizing iteration scheme for model predictive control based on relaxed barrier functions. Automatica, 2017, 80: 328-339 doi: 10.1016/j.automatica.2017.02.001 期刊类型引用(9)
1. 张浩,解艺旋,杨晓楠,岳晨阳,马明,夏立坤. 融合迭代自组织数据分析的肺结节自动检测算法. 计算机与数字工程. 2024(05): 1336-1340+1358 . 百度学术
2. 崔冬梅,魏国辉,张喜科,马志庆. 肺部肿瘤影像辅助诊断分类算法研究进展. 中国胸心血管外科临床杂志. 2022(07): 892-898 . 百度学术
3. 张战成,张大龙,罗晓清. 带有语义特征得分的肺结节检测方法. 计算机应用. 2020(03): 925-930 . 百度学术
4. 于胡飞,温景熙,辛江,唐艳. 基于生成对抗网络的医学数据域适应研究. 大数据. 2020(05): 45-54 . 百度学术
5. 杨举平. CT与普通放射检查对肺部孤立性球形病变的诊断对比. 影像研究与医学应用. 2019(02): 85-86 . 百度学术
6. 霍其润,李建武,陆耀,秦明. 基于变分的CT图像环形伪影校正. 自动化学报. 2019(09): 1713-1726 . 本站查看
7. 顾广海. 螺旋CT扫描技术在肺结节患者临床诊断中的应用效果分析. 中国医药指南. 2019(35): 28-29 . 百度学术
8. 李祥霞,李彬,田联房,张莉,朱文博. 基于稀疏表示和随机游走的磨玻璃型肺结节分割. 自动化学报. 2018(09): 1637-1647 . 本站查看
9. 于超,王璐,程道文. 基于图像降阶的纹理特征挖掘方法在医学上的应用. 东北师大学报(自然科学版). 2018(03): 53-57 . 百度学术
其他类型引用(31)
-