基于Hammerstein-Wiener模型的广义预测控制

李 泰，侯小燕，林鹤云

（1.东南大学电气工程学院，江苏南京210096；2.江苏科技大学电子信息学院，江苏镇江212003）

基于Hammerstein-Wiener模型的广义预测控制

李 泰1，2，侯小燕2，林鹤云1

（1.东南大学电气工程学院，江苏南京210096；2.江苏科技大学电子信息学院，江苏镇江212003）

提出了一种新型的基于Hammerstein-Wiener模型的广义预测控制策略。采用基于最小二乘支持向量机的Hammerstein-Wiener模型描述非线性系统动态特性，作为被控对象预测模型。同时，针对现有遗传算法和混沌粒子群优化算法收敛速度慢和精度低等缺点，给出一种拟牛顿信赖域混沌粒子群混合优化算法，作为预测控制的滚动优化策略，函数测试和非线性对象的广义预测控制的滚动优化表明该算法的优越性。最后，对设计的预测控制器进行实例仿真，结果表明它能满足系统实时稳定运行的需求，取得了良好的控制效果。

广义预测控制；Hammerstein-Wiener模型；拟牛顿信赖域；混沌粒子群

0 引 言

广义预测控制（generalized predictive control，GPC）［1］是随着自适应技术发展起来的一种预测控制方法，它控制效果好，鲁棒性强，适用于具有时滞或开环不稳定的非线性系统，因此在工业控制中得到大量应用。然而，广义预测控制通常需要进行Diophantine方程求解、矩阵求逆、最小二乘递推等复杂计算，求解高阶复杂非线性系统控制问题时计算量较大。广义预测模型常采用CARIMA模型、状态空间模型等参数模型，但是针对复杂非线性系统，预测模型仍没有统一有效的建模方法。滚动优化采用对控制量和偏差加权的性能指标，常用最小二乘递推寻优，需要进行大量数据和复杂的矩阵计算，且寻优性能常在快速性和稳定性之间折中，其优化性能有待进一步提高。

文献［2］采用ARMAX模型作为广义预测控制的预测模型，设计了预测控制器，但是模型参数未知或存在动态不确定性时，该方法不能实现参数的一致估计和性能指标的极小化；文献［3］利用神经网络对非线性系统进行逼近和控制，但神经网络对大范围非线性系统的逼近存在逼近精度和网络复杂性的矛盾。T-S模型在非线性对象的建模和控制中也得到广泛应用，文献［4－5］采用T-S模糊模型对被控对象建模，实现了模型实时更新，提高了被控系统的逼近精度，但其结构辨识与参数辨识的混合导致计算量大。Hammerstein-Wiener模型［6－9］作为一种经典的分块结构非线性模型，包括了两个非线性模块，能更好地描述非线性系统，且结构简单，计算量少，为非线性预测模型提供了良好的模型结构。

文献［10］将遗传算法作为广义预测控制中的滚动优化来处理带约束的优化问题，但遗传算法在函数寻优时容易陷入局部极值，产生早熟且进化代数较多。粒子群优化（particle swarm optimization，PSO）算法［11－13］也被引入到预测控制的滚动优化中，解决了被控对象存在约束时寻优困难和计算复杂的问题，但PSO易产生早熟收敛；混沌粒子群优化算法［14－15］结合混沌和粒子群算法的优点，保持了全局寻优的快速性且避免陷入局部极值，但算法后期收敛精度和收敛速度仍未得到改善。文献［16］将拟牛顿法与信赖域法相结合，在信赖域求解子问题不理想时，采用拟牛法可以避免信赖域半径的盲目缩小和重新求解子问题，加快了寻优速度。

本文将Hammerstein-Wiener模型、拟牛顿信赖域（quasi-Newton trust region，QN-TR）混沌粒子群优化（chaotic particle swarm optimization，CPSO）算法和广义预测控制的思想相结合，提出一种新型的基于Hammerstein-Wiener模型的广义预测控制策略。为了提高CPSO算法效率，将拟牛顿信赖域和混沌粒子群相结合，可以避免陷入局部极值，并且加快局部极值的收敛速度。将拟牛顿信赖域混沌粒子群优化算法作为预测控制的滚动优化策略，可避免复杂的矩阵计算且提高滚动优化的性能。仿真结果验证了该方法的有效性。

1 广义预测控制的基本思想

设非线性被控对象的数学模型用如下具有随机阶跃扰动非平稳噪声的离散差分方程描述为

式中，y，u，ε为系统输出、输入和均值为零、方差为δ2的白噪声；Δ为差分因子，，其中ai，bi，ci为系统参数。

将式（1）两边同时乘以差分因子Δ得：

引入Diophantine方程：

将式（3）代入式（1）并化简得

式中，yp（（k＋j）／k）＝B（z－1）Rj（z－1）Δu（k＋j－1）＋Sj（z－1）y（k）为预测输出；εp（k＋j）＝Rj（z－1）ε（k＋j）为预测误差。

广义预测控制问题可以表达为在模型约束和控制约束下求出使性能指标达到最小的一组最优控制量序列。采用如下二次型性能指标：

式中，P为预测长度；M为控制长度（M≤P）；λj为控制加权系数；yr（k＋j）为参考轨迹。

2 Hammerstein-Wiener预测模型辨识

在非线性系统动态建模和控制中，Hammerstein-Wiener模型作为一类块结构非线性动态模型受到了高度重视。该模型同时结合了线性动态模型和静态非线性函数模型，结构简单，容易辨识，能较好地反映过程特征，已被广泛应用于系统建模［9］。

单输入单输出非线性Hammerstein-Wiener模型结构如图1所示，它包括静态输入非线性模块f，动态线性模块G和静态输出非线性模块h。

图1 Hammerstein-Wiener模型

本文采用基于最小二乘支持向量机（least squares support vector machine，LS-SVM）的Hammerstein-Wiener模型对非线性系统进行建模，通过模型辨识得到非线性系统的预测模型［17］。

对于单输入单输出系统，基于SVM的Hammerstein-Wiener模型辨识步骤如下：

步骤1 采集足够的输入输出数据样本。

步骤2 根据Hammerstein-Wiener各个模块模型结构，得到Hammerstein-Wiener非线性模型的输出数学表达式

步骤3 采用SWM表示Hammerstein-Wiener模型表达式中的非线性函数。

步骤4 定义相应的优化问题和约束条件，利用最小二乘对上述非线性函数参数进行寻优，得到待定参数d1、d0，再进行奇异值分解得到回归参数ai、bj，从而得到非线性部分f（x）和g（y）。

步骤5 选择高斯径向基函数（radial basis function，RBF）为核函数，得到Hammerstein-Wiener模型的最终表达式

系统静态非线性模块的输出即是整个模型的输出，对系统未来时刻输出进行多步预测，得到多步预测输出，将其反馈到输入端，从而对非线性系统进行提前预测控制。

3 基于混合优化算法的预测控制滚动优化

3.1 拟牛顿信赖域CPSO混合优化算法

信赖域法是一种求解非线性约束优化问题的算法，它可以同时确定搜索方向和步长，简化搜索过程。将拟牛顿法引入到信赖域中［16］，采用BFGS拟牛顿公式修正信赖域模型子问题中的Bk，同时，在每个迭代步，优先使用信赖域方法，当试探步不成功时，采用拟牛顿步继续迭代，避免重新求解子问题及信赖域的盲目缩小，加快寻优速度而且算法具有二次终止性。

该算法基本思想是迭代求解信赖域模型子问题：

式中，f（xk）为f（x）在xk处的函数值；s为尝试迭代步；Bk为近似于Hessian矩阵Δ2f（xk）的对称矩阵；gk＝Δf（xk）为f（x）在xk处的梯度；Δk为信赖域半径。

CPSO是一种基于种群的随机优化技术，对待优化目标能实现快速全局搜索。将CPSO作为全局搜索器［14，18］，用拟牛顿信赖域［16］加快局部搜索，能提高收敛速度，取得较好的寻优效果。

基于拟牛顿信赖域的混沌粒子群混合优化算法步骤如下：

步骤1 确定群体规模M、最大函数评价次数Mm、确定算法权重w，粒子速度范围［－vmax，vmax］初始化粒子群随机初始位置和速度。

步骤2 将每个粒子的个体极值Pi设置为当前位置，根据适应度函数（n是训练样本数，分别为实际值和预测值）计算每个粒子的适应度值，取适应度值最好的粒子的个体极值为CPSO最初的全局极值Pg，函数评价次数k＝M。

步骤3 若k≥Mm，则Pbest＝min｛Pg，Pl｝，转步骤10，否则继续。

步骤4 用PSO速度和位置更新公式

对粒子速度和位置进行更新。

步骤5 根据各个粒子的适应度值更新Pi和Pg，并记录全局最优粒子下标gbest，更新k。

步骤7 以全局最优粒子位置xgbest为初始点，运行拟牛顿信赖域算法，更新xlbest，Pl和k。

步骤8 若k≥Mm，则Pbest＝min｛Pg，Pl｝，转步骤10，否则继续。

步骤9 对当前粒子群Pg进行混沌扰动。首先将Pg映射为定义域［0 1］之间的混沌变量r，若r＜Pm（Pm为混沌变换概率），则用Logistic［19］映射式Z′k＝μZk（1－Zk）（其中μ＝4，0＜Zk＜1）进行迭代，得到n个混沌变量，这些变量通过逆映射获得n个粒子，对粒子适应度值进行计算和排序，从而获到最优解P′g，令Pg＝P′g，更新k，转步骤3。

步骤10 输出粒子群最优值，算法结束。

其中，拟牛顿信赖域算法达到收敛即停止迭代，不需要满足最大函数评价次数的终止条件。

3.2 寻优性能比较

为了验证拟牛顿信赖域CPSO混合优化算法的优越性，本文将其与遗传算法（genetic algorithm，GA）、CPSO算法分别对4个给定初始点的无约束优化问题测试函数［20］进行寻优，并对各自的函数最优值、仿真时间和收敛精度进行比较。

表1 3种算法寻优性能比较

对每个函数进行寻优，参数设置如下：信赖域半径Δ0＝‖g（x0）‖，最小信赖域半径Δmin＝0.2，μ＝0.75，正定矩阵B0＝I，精确度ε＝10－5，c1＝0.25，c2＝2，δ＝0.1。粒子群规模M＝30，加速常数c′1＝c′2＝2，惯性权重w取文献［19］中的线性惯性权重，wmax＝0.9，wmin＝0.4，vmax＝（xmax－xmin）／2，vmax＝（xmax－xmin）／2，vmin＝－vmax，ε＝0.01，收敛阈值ε′＝10－3，最大函数评价次数3 000，混沌概率pm＝0.2，混沌搜索区间［－0.5pg，2.5pg］。CPSO算法的参数取值同上，GA算法取种群规模30，交叉概率0.8，变异概率0.002。

对每个函数测试50次，测试函数维数均为2，寻优结果取50次搜索得到最优值的平均值，仿真时间取50次仿真的平均时间。得到最终结果如表1所示，其中：QN-TRCPSO为拟牛顿信赖域混沌粒子群混合优化算法，CPSO为混沌粒子群优化算法，GA为遗传算法，f＊为标准的函数最优值，f为搜索得到的函数最优值，T为仿真时间，CR为收敛率。

3.3 基于混合优化算法的滚动优化预测控制

鉴于拟牛顿信赖域CPSO混合优化算法函数寻优的准确性和快速性，本文将其应用于广义预测控制的滚动优化。以Hammerstein-Wiener模型为预测模型，以拟牛顿信赖域CPSO混合优化算法为滚动优化策略的广义预测控制结构如图2所示。其预测控制步骤如下：

步骤1 设置预测控制初始参数、拟牛顿信赖域初始参数、CPSO初始参数，基于SVM的Hammerstein-Wiener模型辨识初始参数。

步骤2 通过Hammerstein-Wiener模型辨识得到系统多步预测输出。

步骤3 由系统多步预测输出和参考轨迹输出计算性能指标式（5）。

步骤4 采用拟牛顿信赖域CPSO混合优化算法进行寻优，输出使性能指标达到最小的控制量u（k）。

步骤5 转步骤2，将得到的u（k）重新作用于系统，得到下一步多步预测输出，进行下一周期的预测控制。

图2 基于Hammerstein-Wiener模型的广义预测控制

3.4 滚动优化效果比较

取以下非线性系统为研究对象：

式中，E（k）为［－0.2，0.2］范围内的均匀白噪声，来验证拟牛顿信赖域CPSO混合优化算法寻优的优越性。分别采用LS-SVM、CPSO及本文的拟牛顿信赖域混合优化算法（QN-TR-CPSO）作为滚动优化策略，得到仿真结果如图3所示。

图3 3种优化算法仿真比较

其中，Yr为参考输入，LS、CPSO和QN-TR-CPSO分别为采用最小二乘、混沌粒子群及拟牛顿信赖域混沌粒子群优化算法作为滚动优化策略得到的预测控制结果。可以看出，拟牛顿信赖域混沌粒子群优化算法使输出快速稳定跟随给定输入，超调较小，调节速度最快，取得了良好的控制效果。

4 仿真实例研究

4.1 非线性系统Hammerstein-Wiener模型辨识

以如下非线性Hammerstein-Wiener系统为被控对象，验证所给出的基于Hammerstein-Wiener模型和混合优化算法的广义预测控制策略的实时控制效果。

被控系统包括两个静态非线性模块和一个动态线性模块，其中，静态输入部分为一非线性死区函数：

动态线性部分的传递函数：

非线性输出部分的逆：

采集200组系统实时输入输出数据，利用前面所述的最小二乘法辨识基于SVM的Hammerstein-Wiener模型线性部分递归参数和非线性部分参数［17］，得到系统的预测模型，其中各个预测参数如下：

图4为基于SVM的Hammerstein-Wiener预测模型与系统实际输出的对比。可以看出，Hammerstein-Wiener模型对复杂非线性对象有较好的辨识能力，该方法简单且辨识精度高。作为全局的N-L-N模型，辨识参数会随非线性动态系统变更而更新，能随时反映系统动态过程。

图4 实际输出与预测模型输出对比

4.2 基于Hammerstein-Wiener模型的广义预测控制

将上述非线性Hammerstein-Wiener系统模型作为预测模型，采用拟牛顿信赖域混沌粒子群混合优化算法作为滚动优化策略，构成了基于Hammerstein-Wiener模型的非线性系统的广义预测控制系统。控制系统结构如图2所示，取预测长度P＝6，控制长度M＝2，控制加权系数λj＝0.6，柔化系数α＝0.4，最大函数评价次数Mm＝3000，得到非线性系统实时跟踪控制仿真曲线如图5所示。

仿真结果表明，在系统输入发生变化时，非线性系统能够快速稳定跟随系统给定输入，输出波动小，控制精度高，取得了良好的控制效果。

图5 非线性系统的广义预测控制

5 结 论

针对非线性系统建模困难和广义预测控制滚动优化性能存在的问题，提出了一种基于Hammerstein-Wiener辨识模型和新型滚动优化算法的广义预测控制方法。采用基于支持向量机的非线性Hammerstein-Wiener模型作为预测模型，有效提高了预测模型对研究对象的逼近精度；同时给出了一种拟牛顿信赖域混沌粒子群混合优化算法，并将其作为广义预测控制的滚动优化策略，仿真测试表明该算法的优越性。最后，将该控制策略应用于非线性Hammerstein-Wiener系统的预测控制，结果表明系统输出能快速稳定跟随给定输入，取得了良好的控制效果。

［1］Clarke D W，Mohtadi C，Tuffs P S.Generalized predictive control［J］.Automatica，1987，23（2）：137－148.

［2］Yuan Z D.Generalized predictive control based on ARMAX model［J］.Control Theory ＆Applications，1988，5（1）：12－17.（袁震东.基于ARMAX模型的广义预侧控制［J］.控制理论与应用，1988，5（1）：12－17.）

［3］Akpan V A，Hassapis G D.Nonlinear model identification and adaptive model predictive control using neural networks［J］.Trans.on ISA，2011，50（2）：177－194.

［4］Wang S X，Dong C，Liu H R.Generalized nonlinear predictive controller based on T-S fuzzy model and small-world optimization algorithm［J］.Control and Decision，2011，26（5）：673－678.（王爽心，董旸，刘海瑞.基于T-S模型和小世界优化算法的广义非线性预测控制［J］.控制与决策，2011，26（5）：673－678.）

［5］Zhang S T，Bai S Z.Controller design of uncertain nonlinear systems based on T-S fuzzy model［J］.Control Theory ＆Applications，2009，7（2）：139－143.

［6］Zhu Y C.Estimation of an N-L-N Hammerstein-Wiener model［J］.Automatica，2002，38：1067－1614.

［7］Wills A，Schön T B，Ljung L，et al.Identification of Hammerstein-Wiener models［J］.Automatica，2013，49（1）：70－81.

［8］Ding B，Huang B.Output feedback model predictive control for nonlinear systems represented by Hammerstein-Wiener model［J］.Control Theory ＆Applications，2007，1（5）：1302－1310.

［9］Bloemen H H J，Van Den Boom T J J，Verbruggen H B.Modelbased predictive control for Hammerstein-Wiener systems［J］.International Journal of Control，2001，74（5）：482－495.

［10］Tong Z N，Xiao L，Peng K X，et，al.Constrained generalized predictive control of mould level based on genetic algorithm［J］.Control and Decision，2009，24（11）：1735－1739.（童朝南，肖磊，彭开香，等.一基于遗传算法的结晶器液位约束广义预测控制［J］.控制与决策，2009，24（11）：1735－1739.）

［11］Xiao B X，Zhu Z G，Liu Y F.Research of hybrid optimized generalized predictive controller based on particle swarm optimization［J］.Journal of System Simulation，2007，19（4）：820－824.（肖本贤，朱志国，刘一福.基于粒子群算法混合优化的广义预测控制器研究［J］.系统仿真学报，2007，19（4）：820－824.）

［12］Wang Z H，Sun Y X.Generalized predictive control based on particle swarm optimization for linear／nonlinear process with constraints［C］∥Proc.of the Second International Conference on Computational Intelligence and Natural Computing，2010：303－306.

［13］Zou Q Y，Ji J W，Ma L L.Greenhouse air temperature predictive control using the particle swarm optimisation algorithm［C］∥Proc.of the 4th International Symposium on Intelligent Information Technology in Agriculture，2007：495－499.

［14］Jiang H M，Kwong C K，Chen Z Q，et al.Chaos particle swarm optimization and T-S fuzzy modeling approaches to constrained predictive control［J］.Knowledge-Based Systems，2012，39（1）：194－201.

［15］Huynh D C，Dunnigan M W，Finney S J.Energy efficient control of an induction machine using a chaos particle swarm optimization algorithm［C］∥Proc.of the IEEE International Conference on Power and Energy，2010：450－455.

［16］Gertz E M.A quasi-Newton trust-region method［J］.Mathematical Programming，2004，100（3）：447－470.

［17］Gui W H，Song H Y，Yang C H.Hammerstein-Wiener model identified by least-square-support-vector machine and its application［J］.Control Theory ＆Applications，2008，25（3）：393－397.（桂卫华，宋海鹰，阳春华.Hammerstein-Wiener模型最小二乘向量机辨识及其应用［J］.控制理论与应用，2008，25（3）：393－397.）

［18］Liu J M，Gao Y L.Chaos particle swarm optimization algorithm［J］.Computer Applications，2008，28（2）：322－325.（刘军民，高岳林.混沌粒子群优化算法［J］.计算机应用，2008，28（2）：322－325.）

［19］Liu L，Qian W M，Qian F.An Improved Chaos-Particle Swarm Optimization Algorithm［J］.Journal of East China University of Science and Technology（Natural Science Edition），2010，36（2）：267－272.（刘玲，钟伟民，钱锋.改进的混沌粒子群优化算法［J］.华东理工大学学报（自然科学版），2010，36（2）：267－272.）

［20］Sun W Y，Yuan Y X.Optimization theory and methods：non-linear programming［M］.New York：Springer-Verlag，2006：637－647.

Generalized predictive control based on Hammerstein-Wiener model

LI Tai1，2，HOU Xiao-yan2，LIN He-yun1
（1.School of Electrical Engineering，Southeast University，Nanjing 210096，China；2.School of Electronic Information，Jiangsu University of Science and Technology，Zhenjiang 212003，China）

A novel generalized predictive control（GPC）strategy based on the Hammerstein-Wiener model is proposed.The dynamic characteristics of the nonlinear system are described by the Hammerstein-Wiener model based on the support vector machine，so a prediction model of the controlled object is obtained.Furthermore，an optimization algorithm of chaotic particle swarm combined with quasi-Newton trust region（QN-TR）is proposed in order to avoid the deficiency of slow convergence speed and low accuracy of the genetic algorithm and the chaotic particle swarm optimization（CPSO）algorithm，so a rolling optimization strategy of the predictive control is obtained.Function tests and rolling optimization of the GPC to the nonlinear object reflect the superiority of the algorithm.Finally，the results of the simulation example for the generalized predictive controller show that it can meet the demand of real-time and stable operation of the system，and a good control effect is obtained.

generalized predictive control（GPC）；Hammerstein-Wiener model；quasi-Newton trust region（QN-TR）；chaotic particle swarm optimization（CPSO）

TP 319

A

10.3969／j.issn.1001-506X.2015.08.24

李 泰（1982－），男，讲师，博士后，主要研究方向为复杂非线性系统的智能控制、建模与优化。

E-mail：teilytl＠126.com

侯小燕（1990－），女，硕士研究生，主要研究方向为混沌优化、非线性系统的预测控制。

E-mail：1181501030＠qq.com

林鹤云（1965－），男，教授，博士生导师，主要研究方向为复杂电力传动系统智能控制。

E-mail：185441638＠qq.com

1001－506X201508－1874－06

网址：www.sys－ele.com

2014－09－03；

2014－11－13；网络优先出版日期：2015－01－04。

网络优先出版地址：http：／／www.cnki.net／kcms／detail／11.2422.TN.20150104.1720.011.html

国家自然科学基金（51307074）；江苏省博士后基金（1301005B）资助课题