CN108153151B - Mimo全格式无模型控制器基于系统误差的参数自整定方法 - Google Patents

Abstract

本发明公开了一种MIMO全格式无模型控制器基于系统误差的参数自整定方法，利用系统误差集作为BP神经网络的输入，BP神经网络进行前向计算并通过输出层输出惩罚因子、步长因子等MIMO全格式无模型控制器待整定参数，采用MIMO全格式无模型控制器的控制算法计算得到针对被控对象的控制输入向量，以系统误差函数的值最小化为目标，采用梯度下降法，并结合控制输入分别针对各个待整定参数的梯度信息集，进行系统误差反向传播计算，在线实时更新BP神经网络的隐含层权系数、输出层权系数，实现控制器基于系统误差的参数自整定。本发明提出的MIMO全格式无模型控制器基于系统误差的参数自整定方法，能有效克服控制器参数的在线整定难题，对MIMO系统具有良好的控制效果。

Description

MIMO全格式无模型控制器基于系统误差的参数自整定方法

技术领域

本发明属于自动化控制领域，尤其是涉及一种MIMO全格式无模型控制器基于系统误差 的参数自整定方法。

背景技术

MIMO(Multiple Input and Multiple Output，多输入多输出)系统的控制问题，一直以来 都是自动化控制领域所面临的重大挑战之一。

MIMO控制器的现有实现方法中包括MIMO全格式无模型控制器。MIMO全格式无模型控制器是一种新型的数据驱动控制方法，不依赖被控对象的任何数学模型信息，仅依赖于MIMO被控对象实时测量的输入输出数据进行控制器的分析和设计，并且实现简明、计算负担小及鲁棒性强，对未知非线性时变MIMO系统也能够进行很好的控制，具有良好的应用前景。MIMO全格式无模型控制器的理论基础，由侯忠生与金尚泰在其合著的《无模型自适应控制—理论与应用》(科学出版社，2013年，第116页)中提出，其控制算法如下：

其中，u(k)为k时刻的控制输入向量，u(k)＝[u₁(k),…,u_mu(k)]^T，mu为控制输入个数， Δu(k)＝u(k)-u(k-1)；e(k)为k时刻的误差向量，e(k)＝[e₁(k),…,e_my(k)]^T，my为输出个 数；Δy(k)＝y(k)-y(k-1)，y(k)为k时刻的MIMO系统输出实际值向量，y(k)＝ [y₁(k),…,y_my(k)]^T；

为k时刻的MIMO系统伪分块雅克比矩阵估计值，

为

的 第i块(i＝1,…,Ly+Lu)，

为矩阵

的2范数；λ为惩罚因子，ρ₁,…,ρ_Ly+Lu为步长因子，Ly为控制输出线性化长度常数，Lu为控制输入线性化长度常数。

然而，MIMO全格式无模型控制器在实际投用前需要依赖经验知识来事先设定惩罚因子 λ和步长因子ρ₁,…,ρ_Ly+Lu等参数的数值，在实际投用过程中也尚未实现惩罚因子λ和步长因 子ρ₁,…,ρ_Ly+Lu等参数的在线自整定。参数有效整定手段的缺乏，不仅使MIMO全格式无模 型控制器的使用调试过程费时费力，而且有时还会严重影响MIMO全格式无模型控制器的控 制效果，制约了MIMO全格式无模型控制器的推广应用。也就是说：MIMO全格式无模型控 制器在实际投用过程中还需要解决在线自整定参数的难题。

为此，为了打破制约MIMO全格式无模型控制器推广应用的瓶颈，本发明提出了一种 MIMO全格式无模型控制器基于系统误差的参数自整定方法。

发明内容

为了解决背景技术中存在的问题，本发明的目的在于，提供一种MIMO全格式无模型控 制器基于系统误差的参数自整定方法。

为此，本发明的上述目的通过以下技术方案来实现，包括以下步骤：

步骤(1)：针对具有mu个输入(mu为大于或等于2的整数)与my个输出(my为大 于或等于2的整数)的MIMO(Multiple Input and Multiple Output，多输入多输出)系统，采 用MIMO全格式无模型控制器进行控制；确定MIMO全格式无模型控制器的控制输出线性 化长度常数Ly，Ly为大于或等于1的整数；确定MIMO全格式无模型控制器的控制输入线 性化长度常数Lu，Lu为大于或等于1的整数；所述MIMO全格式无模型控制器参数包含惩 罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu；确定MIMO全格式无模型控制器待整定参数，所述MIMO 全格式无模型控制器待整定参数，为所述MIMO全格式无模型控制器参数的部分或全部，包 含惩罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu的任意之一或任意种组合；确定BP神经网络的输入层 节点数、隐含层节点数、输出层节点数，所述输出层节点数不少于所述MIMO全格式无模型 控制器待整定参数个数；初始化所述BP神经网络的隐含层权系数、输出层权系数；

步骤(2)：将当前时刻记为k时刻；

步骤(3)：基于MIMO系统的第jy个输出期望值与第jy个输出实际值(1≤jy≤my)，采用第jy个误差计算函数计算得到k时刻的第jy个误差，记为e_jy(k)；将第jy个误差及其函数组、第jy个输出期望值、第jy个输出实际值的任意之一或任意种组合，记为集合{误差组jy}，并放入集合{系统误差集}；针对MIMO系统的其他my-1个输出，重复执行本步骤，直 至所述集合{系统误差集}包含全部{{误差组1},…,{误差组my}}的集合，并同时得到由my个误差所构成的误差向量e(k)＝[e₁(k),…,e_my(k)]^T，然后进入步骤(4)；

步骤(4)：将步骤(3)得到的所述集合{系统误差集}作为BP神经网络的输入，所述BP神经网络进行前向计算，计算结果通过所述BP神经网络的输出层输出，得到所述MIMO 全格式无模型控制器待整定参数的值；

步骤(5)：基于步骤(3)得到的所述误差向量e(k)、步骤(4)得到的所述MIMO全 格式无模型控制器待整定参数的值，采用MIMO全格式无模型控制器的控制算法，计算得到MIMO全格式无模型控制器针对被控对象在k时刻的控制输入向量u(k)＝[u₁(k),…, u_mu(k)]^T；

步骤(6)：针对步骤(5)得到的所述控制输入向量u(k)中的第ju个控制输入u_ju(k)(1≤ju≤mu)，计算所述第ju个控制输入u_ju(k)分别针对各个所述MIMO全格式无模型控制器待整定参数在k时刻的梯度信息，具体计算公式如下：

当所述MIMO全格式无模型控制器待整定参数中包含惩罚因子λ且Lu＝1时，所述第ju 个控制输入u_ju(k)针对所述惩罚因子λ在k时刻的梯度信息为：

当所述MIMO全格式无模型控制器待整定参数中包含惩罚因子λ且Lu>1时，所述第ju 个控制输入u_ju(k)针对所述惩罚因子λ在k时刻的梯度信息为：

当所述MIMO全格式无模型控制器待整定参数中包含步长因子ρ_i且1≤i≤Ly时，所述第 ju个控制输入u_ju(k)针对所述步长因子ρ_i在k时刻的梯度信息为：

当所述MIMO全格式无模型控制器待整定参数中包含步长因子ρ_Ly+1时，所述第ju个控 制输入u_ju(k)针对所述步长因子ρ_Ly+1在k时刻的梯度信息为：

当所述MIMO全格式无模型控制器待整定参数中包含步长因子ρ_i且Ly+2≤i≤Ly+Lu且 Lu>1时，所述第ju个控制输入u_ju(k)针对所述步长因子ρ_i在k时刻的梯度信息为：

其中，Δu_j(k)＝u_j(k)-u_j(k-1)，Δy_j(k)＝y_j(k)-y_j(k-1)，

为k时刻的MIMO系统伪 分块雅克比矩阵估计值，

为

的第i块(i＝1,…,Ly+Lu)，

为矩阵

的 第jy行第ju列元素，

为矩阵

的2范数；

上述全部所述梯度信息的集合记为{梯度信息ju}，放入集合{梯度信息集}；

针对步骤(5)得到的所述控制输入向量u(k)中的其他mu-1个控制输入，重复执行本 步骤，直至所述集合{梯度信息集}包含全部{{梯度信息1},…,{梯度信息mu}}的集合，然后 进入步骤(7)；

步骤(7)：以系统误差函数的值最小化为目标，采用梯度下降法，结合步骤(6)得到的所述集合{梯度信息集}，进行系统误差反向传播计算，更新BP神经网络的隐含层权系数、输出层权系数，作为后一时刻BP神经网络进行前向计算时的隐含层权系数、输出层权系数；

步骤(8)：所述控制输入向量u(k)作用于被控对象后，得到被控对象在后一时刻的my 个输出实际值，返回到步骤(2)，重复步骤(2)到步骤(8)。

在采用上述技术方案的同时，本发明还可以采用或者组合采用以下进一步的技术方案：

所述步骤(3)中的所述第jy个误差计算函数的自变量包含第jy个输出期望值与第jy个 输出实际值。

所述步骤(3)中的所述第jy个误差计算函数采用

其中

为 k时刻设定的第jy个输出期望值，y_jy(k)为k时刻采样得到的第jy个输出实际值；或者采用

其中

为k+1时刻的第jy个输出期望值，y_jy(k)为k时刻 采样得到的第jy个输出实际值。

所述步骤(3)中的所述第jy个误差及其函数组，包含k时刻的第jy个误差e_jy(k)、k时 刻及之前所有时刻的第jy个误差的累积即

k时刻第jy个误差e_jy(k)的一阶后向差分 e_jy(k)-e_jy(k-1)、k时刻第jy个误差e_jy(k)的二阶后向差分e_jy(k)-2e_jy(k-1)+e_jy(k-2)、 k时刻第jy个误差e_jy(k)的高阶后向差分的任意之一或任意种组合。

所述步骤(7)中的所述系统误差函数的自变量包含my个误差、my个输出期望值、my个输出实际值的任意之一或任意种组合。

所述步骤(7)中的所述系统误差函数为

其中，e_jy(k)为第 jy个误差，Δu_ju(k)＝u_ju(k)-u_ju(k-1)，a_jy与b_ju为大于或等于0的常数，1≤jy≤my， 1≤ju≤mu。

本发明提供的MIMO全格式无模型控制器基于系统误差的参数自整定方法，能够实现良 好的控制效果，并有效克服惩罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu需要费时费力进行整定的难 题。

附图说明

图1为本发明的原理框图；

图2为本发明采用的BP神经网络结构示意图；

图3为两输入两输出MIMO系统在惩罚因子λ和步长因子ρ₁,ρ₂,ρ₃,ρ₄同时自整定时第 1个输出的控制效果图；

图4为两输入两输出MIMO系统在惩罚因子λ和步长因子ρ₁,ρ₂,ρ₃,ρ₄同时自整定时第 2个输出的控制效果图；

图5为两输入两输出MIMO系统在惩罚因子λ和步长因子ρ₁,ρ₂,ρ₃,ρ₄同时自整定时的 控制输入图；

图6为两输入两输出MIMO系统在惩罚因子λ和步长因子ρ₁,ρ₂,ρ₃,ρ₄同时自整定时的 惩罚因子λ变化曲线；

图7为两输入两输出MIMO系统在惩罚因子λ和步长因子ρ₁,ρ₂,ρ₃,ρ₄同时自整定时的 步长因子ρ₁,ρ₂,ρ₃,ρ₄变化曲线；

图8为两输入两输出MIMO系统在惩罚因子λ固定而步长因子ρ₁,ρ₂,ρ₃,ρ₄自整定时第 1个输出的控制效果图；

图9为两输入两输出MIMO系统在惩罚因子λ固定而步长因子ρ₁,ρ₂,ρ₃,ρ₄自整定时第 2个输出的控制效果图；

图10为两输入两输出MIMO系统在惩罚因子λ固定而步长因子ρ₁,ρ₂,ρ₃,ρ₄自整定时的 控制输入图；

图11为两输入两输出MIMO系统在惩罚因子λ固定而步长因子ρ₁,ρ₂,ρ₃,ρ₄自整定时的 步长因子ρ₁,ρ₂,ρ₃,ρ₄变化曲线。

具体实施方式

下面结合附图和具体实施例对本发明进一步说明。

图1给出了本发明的原理框图。针对具有mu个输入(mu为大于或等于2的整数)与my个输出(my为大于或等于2的整数)的MIMO系统，采用MIMO全格式无模型控制器进 行控制；确定MIMO全格式无模型控制器的控制输出线性化长度常数Ly，Ly为大于或等于1 的整数；确定MIMO全格式无模型控制器的控制输入线性化长度常数Lu，Lu为大于或等于 1的整数；MIMO全格式无模型控制器参数包含惩罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu；确定 MIMO全格式无模型控制器待整定参数，其为MIMO全格式无模型控制器参数的部分或全部， 包含惩罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu的任意之一或任意种组合；在图1中，MIMO全格 式无模型控制器待整定参数为惩罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu；确定BP神经网络的输入 层节点数、隐含层节点数、输出层节点数，其中输出层节点数不少于MIMO全格式无模型控 制器待整定参数个数；初始化BP神经网络的隐含层权系数、输出层权系数。

将当前时刻记为k时刻；将第jy个输出期望值

与第jy个输出实际值y_jy(k)之差作 为k时刻的第jy个误差e_jy(k)；将第jy个误差e_jy(k)、k时刻及之前所有时刻的第jy个误差 的累积即

k时刻第jy个误差e_jy(k)的一阶后向差分e_jy(k)-e_jy(k-1)的组合，记为 集合{误差组jy}，并放入集合{系统误差集}；针对MIMO系统的其他my-1个输出，重复执 行直至集合{系统误差集}包含全部{{误差组1},…,{误差组my}}的集合，并同时得到由my个 误差所构成的误差向量e(k)＝[e₁(k),…,e_my(k)]^T；然后将集合{系统误差集}作为BP神经网络 的输入，BP神经网络进行前向计算，计算结果通过BP神经网络的输出层输出，得到MIMO 全格式无模型控制器待整定参数的值；基于误差向量e(k)、MIMO全格式无模型控制器待整 定参数的值，采用MIMO全格式无模型控制器的控制算法，计算得到MIMO全格式无模型 控制器针对被控对象在k时刻的控制输入向量u(k)＝[u₁(k),…,u_mu(k)]^T；针对控制输入向量 u(k)中的第ju个控制输入u_ju(k)，计算所述第ju个控制输入u_ju(k)分别针对各个所述MIMO 全格式无模型控制器待整定参数在k时刻的梯度信息，并将全部所述梯度信息的集合记为{梯 度信息ju}，放入集合{梯度信息集}；针对控制输入向量u(k)中的其他mu-1个控制输入， 重复执行直至集合{梯度信息集}包含全部{{梯度信息1},…,{梯度信息mu}}的集合；随后， 结合所述集合{梯度信息集}，以系统误差函数的值最小化为目标，图1中以综合考虑全部my 个误差贡献的系统误差函数

的值最小化为目标，采用梯度下降法，进行系统误差反 向传播计算，更新BP神经网络的隐含层权系数、输出层权系数，作为后一时刻BP神经网络 进行前向计算时的隐含层权系数、输出层权系数；控制输入向量u(k)作用于被控对象后，得 到被控对象在后一时刻的my个输出实际值，然后重复执行本段落所述的工作，进行后一时 刻的MIMO全格式无模型控制器基于系统误差的参数自整定过程。

图2给出了本发明采用的BP神经网络结构示意图。BP神经网络可以采用隐含层为单层 的结构，也可以采用隐含层为多层的结构。在图2的示意图中，为简明起见，BP神经网络采 用了隐含层为单层的结构，即采用由输入层、单层隐含层、输出层组成的三层网络结构，输 入层节点数设为3×my个，隐含层节点数设为8个，输出层节点数设为待整定参数个数(图 2中待整定参数个数为Ly+Lu+1个)。输入层的节点数分成my组，每组包含3个节点，其中 第jy组的3个节点与集合{误差组jy}中的误差e_jy(k)、误差的累积

误差e_jy(k)的一 阶后向差分e_jy(k)-e_jy(k-1)分别对应。输出层的节点，与惩罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu分别对应。BP神经网络的隐含层权系数、输出层权系数的更新过程具体为：以系统误差函数的值最小化为目标，图2中以综合考虑全部my个误差贡献的系统误差函数

的值最小化为目标，采用梯度下降法，结合所述集合{梯度信息集}，进行系统误差反向传播计算，从而更新BP神经网络的隐含层权系数、输出层权系数。

以下是本发明的一个具体实施例。

被控对象为典型非线性的两输入两输出MIMO系统：

y₁(k)＝x₁₁(k)

y₂(k)＝x₂₁(k)

其中，α(k)＝1+0.1sin(2kπ/1500)，β(k)＝1+0.1cos(2kπ/1500)。

系统输出期望值y^*(k)如下：

在本具体实施例中，mu＝my＝2。

MIMO全格式无模型控制器的控制输出线性化长度常数Ly的数值通常根据被控对象的 复杂程度和实际的控制效果进行设定，一般在1到5之间，过大会导致计算量大，所以一般 常取1或3，在本具体实施例中Ly取为1；MIMO全格式无模型控制器的控制输入线性化长 度常数Lu的数值也通常根据被控对象的复杂程度和实际的控制效果进行设定，一般在1到 10之间，过小会影响控制效果，过大会导致计算量大，所以一般常取3或5，在本具体实施例中Lu取为3。

BP神经网络采用由输入层、单层隐含层、输出层组成的三层网络结构，输入层节点数设 为6个，隐含层节点数设为8个，输出层节点数设为待整定参数个数。

针对上述具体实施例，共进行了两组试验验证。

第一组试验验证时，图2中BP神经网络的输出层节点数预设为5个，对惩罚因子λ和步长因子ρ₁,ρ₂,ρ₃,ρ₄进行同时自整定，图3为第1个输出的控制效果图，图4为第2个输 出的控制效果图，图5为控制输入图，图6为惩罚因子λ变化曲线，图7为步长因子ρ₁,ρ₂, ρ₃,ρ₄变化曲线。结果表明，本发明的方法通过对惩罚因子λ和步长因子ρ₁,ρ₂,ρ₃,ρ₄进行同 时自整定，能够实现良好的控制效果，并且可以有效克服惩罚因子λ和步长因子ρ₁,ρ₂,ρ₃,ρ₄需要费时费力进行整定的难题。

第二组试验验证时，图2中BP神经网络的输出层节点数为4个，首先将惩罚因子λ固定取值为第一组试验验证时惩罚因子λ的平均值，然后对步长因子ρ₁,ρ₂,ρ₃,ρ₄进行自整定，图8为第1个输出的控制效果图，图9为第2个输出的控制效果图，图10为控制输入图，图 11为步长因子ρ₁,ρ₂,ρ₃,ρ₄变化曲线。结果同样表明，本发明的方法在惩罚因子λ固定时通 过对步长因子ρ₁,ρ₂,ρ₃,ρ₄进行自整定，能够实现良好的控制效果，并且可以有效克服步长 因子ρ₁,ρ₂,ρ₃,ρ₄需要费时费力进行整定的难题。

应该特别指出的是，在上述具体实施例中，将第jy个输出期望值

与第jy个输出实 际值y_jy(k)之差作为k时刻的第jy个误差e_jy(k)，也就是

仅为所述 第jy个误差计算函数中的一种方法；也可以将k+1时刻第jy个输出期望值

与k时刻 第jy个输出y_jy(k)之差作为第jy个误差e_jy(k)，也就是

所述第jy 个误差计算函数还可以采用自变量包含第jy个输出期望值与第jy个输出实际值的其他计算方 法，举例来说，

对上述具体实施例的被控对象而言，采用 上述不同的系统误差计算函数，都能够实现良好的控制效果。

还应该特别指出的是，在上述具体实施例中，集合{误差组jy}包含了误差e_jy(k)、误差的 累积

误差e_jy(k)的一阶后向差分e_jy(k)-e_jy(k-1)的组合，仅为其中一种类型的组 合；集合{误差组jy}还可以采用其他类型的组合，举例来说，包含误差e_jy(k)、误差的累积

误差e_jy(k)的一阶后向差分e_jy(k)-e_jy(k-1)、误差e_jy(k)的二阶后向差分 e_jy(k)-2e_jy(k-1)+e_jy(k-2)、误差e_jy(k)的三阶或四阶或更高阶的后向差分等函数的任意 之一或任意种组合。对上述具体实施例的被控对象而言，采用上述不同的集合{误差组jy}， 都能够实现良好的控制效果。

更应该特别指出的是，在上述具体实施例中，在以系统误差函数的值最小化为目标来更 新BP神经网络的隐含层权系数、输出层权系数时，所述系统误差函数采用了综合考虑全部 my个误差贡献的系统误差函数

仅为所述系统误差函数中的一种函数；所述系统误 差函数还可以采用自变量包含my个误差、my个输出期望值、my个输出实际值的任意之一或 任意种组合的其他函数，举例来说，所述系统误差函数采用

或

也就是采用

的另一种函数形式；再举例来说，所述系统误 差函数采用

其中，e_jy(k)为第jy个误差，Δu_ju(k)＝u_ju(k)- u_ju(k-1)，a_jy与b_ju为大于或等于0的常数，1≤jy≤my，1≤ju≤mu；显然，当b_ju等于0时，所述系统误差函数仅考虑了

的贡献，表明最小化的目标是系统误差最小，也就是追求精度高；而当b_ju大于0时，所述系统误差函数同时考虑

的贡献和

的贡献，表明最小化的目标在追求系统误差小的同时，还追求控制输入变化小，也就是既追求精度高又追 求操纵稳。对上述具体实施例的被控对象而言，采用上述不同的系统误差函数，都能够实现 良好的控制效果；与系统误差函数仅考虑

贡献时的控制效果相比，在系统误差函数同 时考虑

的贡献和

的贡献时其控制精度略有降低而其操纵平稳性则有提高。

最后应该特别指出的是，所述MIMO全格式无模型控制器待整定参数，包含惩罚因子λ 和步长因子ρ₁,…,ρ_Ly+Lu的任意之一或任意种组合；在上述具体实施例中，第一组试验验证时 惩罚因子λ和步长因子ρ₁,ρ₂,ρ₃,ρ₄实现了同时自整定，第二组试验验证时惩罚因子λ固定 而步长因子ρ₁,ρ₂,ρ₃,ρ₄实现了自整定；在实际应用时，还可以根据具体情况，选择待整定 参数的任意种组合，举例来说，步长因子ρ₁,ρ₂固定而惩罚因子λ、步长因子ρ₃,ρ₄实现自整 定；此外，MIMO全格式无模型控制器待整定参数，包括但不限于惩罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu，举例来说，根据具体情况，还可以包括MIMO系统伪分块雅克比矩阵估计值

等参数。

上述具体实施方式用来解释说明本发明，仅为本发明的优选实施例，而不是对本发明进 行限制，在本发明的精神和权利要求的保护范围内，对本发明作出的任何修改、等同替换、 改进等，都落入本发明的保护范围。

Claims (4)

1.MIMO全格式无模型控制器基于系统误差的参数自整定方法，其特征在于包括以下步骤：

步骤(1)：针对具有mu个输入与my个输出的MIMO(Multiple Input and MultipleOutput，多输入多输出)系统，其中mu为大于或等于2的整数，my为大于或等于2的整数，采用MIMO全格式无模型控制器进行控制；确定MIMO全格式无模型控制器的控制输出线性化长度常数Ly，Ly为大于或等于1的整数；确定MIMO全格式无模型控制器的控制输入线性化长度常数Lu，Lu为大于或等于1的整数；所述MIMO全格式无模型控制器参数包含惩罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu；确定MIMO全格式无模型控制器待整定参数，所述MIMO全格式无模型控制器待整定参数，为所述MIMO全格式无模型控制器参数的部分或全部，包含惩罚因子λ和步长因子ρ₁,…,ρ_Ly+Lu的任意之一或任意种组合；确定BP神经网络的输入层节点数、隐含层节点数、输出层节点数，所述输出层节点数不少于所述MIMO全格式无模型控制器待整定参数个数；初始化所述BP神经网络的隐含层权系数、输出层权系数；

步骤(2)：将当前时刻记为k时刻；

步骤(3)：基于MIMO系统的第jy个输出期望值与第jy个输出实际值，其中1≤jy≤my，采用第jy个误差计算函数计算得到k时刻的第jy个误差，记为e_jy(k)；所述第jy个误差计算函数的自变量包含第jy个输出期望值与第jy个输出实际值；将第jy个误差及其函数组、第jy个输出期望值、第jy个输出实际值的任意之一或任意种组合，记为集合{误差组jy}，并放入集合{系统误差集}；针对MIMO系统的其他my-1个输出，重复执行本步骤，直至所述集合{系统误差集}包含全部{{误差组1},…,{误差组my}}的集合，并同时得到由my个误差所构成的误差向量e(k)＝[e₁(k),…,e_my(k)]^T，然后进入步骤(4)；

步骤(4)：将步骤(3)得到的所述集合{系统误差集}作为BP神经网络的输入，所述BP神经网络进行前向计算，计算结果通过所述BP神经网络的输出层输出，得到所述MIMO全格式无模型控制器待整定参数的值；

步骤(5)：基于步骤(3)得到的所述误差向量e(k)、步骤(4)得到的所述MIMO全格式无模型控制器待整定参数的值，采用MIMO全格式无模型控制器的控制算法，计算得到MIMO全格式无模型控制器针对被控对象在k时刻的控制输入向量u(k)＝[u₁(k),…,u_mu(k)]^T；

步骤(6)：针对步骤(5)得到的所述控制输入向量u(k)中的第ju个控制输入u_ju(k)，其中1≤ju≤mu，计算所述第ju个控制输入u_ju(k)分别针对各个所述MIMO全格式无模型控制器待整定参数在k时刻的梯度信息，具体计算公式如下：

当所述MIMO全格式无模型控制器待整定参数中包含惩罚因子λ且Lu＝1时，所述第ju个控制输入u_ju(k)针对所述惩罚因子λ在k时刻的梯度信息为：

当所述MIMO全格式无模型控制器待整定参数中包含惩罚因子λ且Lu>1时，所述第ju个控制输入u_ju(k)针对所述惩罚因子λ在k时刻的梯度信息为：

当所述MIMO全格式无模型控制器待整定参数中包含步长因子ρ_i且1≤i≤Ly时，所述第ju个控制输入u_ju(k)针对所述步长因子ρ_i在k时刻的梯度信息为：

当所述MIMO全格式无模型控制器待整定参数中包含步长因子ρ_Ly+1时，所述第ju个控制输入u_ju(k)针对所述步长因子ρ_Ly+1在k时刻的梯度信息为：

当所述MIMO全格式无模型控制器待整定参数中包含步长因子ρ_i且Ly+2≤i≤Ly+Lu且Lu>1时，所述第ju个控制输入u_ju(k)针对所述步长因子ρ_i在k时刻的梯度信息为：

其中，Δu_j(k)＝u_j(k)-u_j(k-1)，Δy_j(k)＝y_j(k)-y_j(k-1)，

为k时刻的MIMO系统伪分块雅克比矩阵估计值，

为

的第i块，其中i＝(1,…,Ly+Lu)，

为矩阵

的第jy行第ju列元素，

为矩阵

的2范数；

上述全部所述梯度信息的集合记为{梯度信息ju}，放入集合{梯度信息集}；

针对步骤(5)得到的所述控制输入向量u(k)中的其他mu-1个控制输入，重复执行本步骤，直至所述集合{梯度信息集}包含全部{{梯度信息1},…,{梯度信息mu}}的集合，然后进入步骤(7)；

步骤(7)：以系统误差函数的值最小化为目标，采用梯度下降法，结合步骤(6)得到的所述集合{梯度信息集}，进行系统误差反向传播计算，更新BP神经网络的隐含层权系数、输出层权系数，作为后一时刻BP神经网络进行前向计算时的隐含层权系数、输出层权系数；所述系统误差函数的自变量包含my个误差、my个输出期望值、my个输出实际值的任意之一或任意种组合；

步骤(8)：所述控制输入向量u(k)作用于被控对象后，得到被控对象在后一时刻的my个输出实际值，返回到步骤(2)，重复步骤(2)到步骤(8)。

2.根据权利要求1所述的MIMO全格式无模型控制器基于系统误差的参数自整定方法，其特征在于，所述步骤(3)中的所述第jy个误差计算函数采用

其中

为k时刻设定的第jy个输出期望值，y_jy(k)为k时刻采样得到的第jy个输出实际值；或者采用

其中

为k+1时刻的第jy个输出期望值，y_jy(k)为k时刻采样得到的第jy个输出实际值。

3.根据权利要求1所述的MIMO全格式无模型控制器基于系统误差的参数自整定方法，其特征在于，所述步骤(3)中的所述第jy个误差及其函数组，包含k时刻的第jy个误差e_jy(k)、k时刻及之前所有时刻的第jy个误差的累积即

k时刻第jy个误差e_jy(k)的一阶后向差分e_jy(k)-e_jy(k-1)、k时刻第jy个误差e_jy(k)的二阶后向差分e_jy(k)-2e_jy(k-1)+e_jy(k-2)、k时刻第jy个误差e_jy(k)的高阶后向差分的任意之一或任意种组合。

4.根据权利要求1所述的MIMO全格式无模型控制器基于系统误差的参数自整定方法，其特征在于，所述步骤(7)中的所述系统误差函数为

其中，e_jy(k)为第jy个误差，Δu_ju(k)＝u_ju(k)-u_ju(k-1)，a_jy与b_ju为大于或等于0的常数，1≤jy≤my，1≤ju≤mu。

Images