CN113625555B - 一种基于递推子空间辨识的自适应逆控制agv转速控制方法 - Google Patents

Abstract

本发明提供一种基于递推子空间辨识的自适应逆控制AGV转速控制方法，利用递推子空间辨识算法获得AGV运动学***的等效逆模型的基本参数；在此基础上利用自适应逆控制法进行AGV转速自适应逆控制器的设计，并将AGV运动学***的等效逆模型作为AGV转速自适应逆控制器的传递函数，实现对AGV转速的自适应跟踪控制。本发明基于递推子空间辨识的自适应逆控制AGV转速控制方法，该AGV转速控制方法能很好适应AGV参数变化的调速***设计，有效克服控制对象的非线性，并有效抑制外部扰动，有利于改善对AGV路径跟踪控制的品质和精度。

Description

一种基于递推子空间辨识的自适应逆控制AGV转速控制方法

技术领域

本发明涉及AGV调速技术领域，更具体地说，涉及一种基于递推子空间辨识的自适应逆控制AGV转速控制方法。

背景技术

全向AGV的导向轮多采用麦克纳姆轮或球形轮，在轨迹跟踪过程中，自动导引小车的控制模型通常不会考虑万向轮从动对AGV主体的作用，由于忽略了这一影响因素，控制***可能出现不稳定状态。多数研究者建模时将AGV假设为刚体，其运动过程未考虑因车体负载不均或磨损导致的轮子半径变化问题，而AGV实际运动过程中，由于负载不均导致左右轮形变和大小不一都会对轨迹跟踪控制的精度造成较大的影响。

采用机理建模时难以确定模型结构和模型参数，因此目前更多的研究集中于数据驱动的黑箱模型或半机理半数据的灰箱模型。AGV调速***是对AGV进行控制的主要***，其动态性能直接影响到运输***的经济性和安全性。为了合理地整定AGV控制器参数、准确地分析跟踪***的稳定特性，均需要研究AGV运动学***的详细模型。实时辨识是获取模型的一种有效手段。AGV调速***是具有非最小相位、非线性的控制***，应根据其自身特点、不同运行方式和工况条件等对其进行实时动态的建模、辨识。子空间辨识法是一类直接估计线性时不变状态变量模型的黑盒建模方法，其具有无需非线性迭代优化、无需参数初始化、解算速度快、实现对***矩阵参数的一致估计的优点。自适应控制是一种能够自动修正控制规律及控制器特性、以适应对象和扰动的动态特性变化的控制方法。自适应逆控制是典型的自适应控制方式之一。为此在AGV路径跟踪***建模时考虑车体参数的变化，实现自适应参数的控制是AGV动态路径跟踪控制的关键。

发明内容

本发明的目的在于克服现有技术中的缺点与不足，提供一种基于递推子空间辨识的自适应逆控制AGV转速控制方法，该AGV转速控制方法能很好适应AGV参数变化的调速***设计，有效克服控制对象的非线性，并有效抑制外部扰动，有利于改善对AGV路径跟踪控制的品质和精度。

为了达到上述目的，本发明通过下述技术方案予以实现：一种基于递推子空间辨识的自适应逆控制AGV转速控制方法，其特征在于：利用递推子空间辨识算法获得AGV运动学***的等效逆模型的基本参数；在此基础上利用自适应逆控制法进行AGV转速自适应逆控制器的设计，并将AGV运动学***的等效逆模型作为AGV转速自适应逆控制器的传递函数，实现对AGV转速的自适应跟踪控制。

包括以下步骤：

步骤1，采集待辨识AGV运动学***的输入数据和输出数据；

步骤2，利用递推子空间辨识算法，得到子空间预估器模型，构造当前时刻输入数据和输出数据的Hankel矩阵，求解当前时刻的AGV运动学***的预测输出数据；

步骤3，根据AGV的控制动作，确定下一时刻的实际输出数据，计算***实际输出数据和预测输出数据的误差，若超过容许误差，则利用增加变遗忘因子的递推子空间辨识算法更新当前时刻的子空间预估器模型，得到下一时刻的子空间预估器模型，否则不更新当前时刻的子空间预估器模型；

步骤4，重复步骤2和3，实时更新子空间预估器模型，确定AGV运动学***的等效逆模型；

步骤5，采用自适应逆控制法设计AGV转速自适应逆控制器，将AGV转速作为目标信号，将步骤4得到的AGV运动学***的等效逆模型作为控制器的传递函数；将目标信号输入AGV转速自适应逆控制器得到控制信号，再用控制信号激励AGV，实现对AGV转速的自适应跟踪控制。

步骤1中，采用AGV的位置作为AGV运动学***的输入数据，采用AGV轮子的速度作为AGV运动学***的输出数据。

步骤2中，所述利用递推子空间辨识算法，得到子空间预估器模型是指：采用递推子空间辨识算法在线获得AGV运动学***的输入数据u和输出数据y，实现对AGV运动学***矩阵[A,B,C,D]的在线递推估计，AGV运动学***矩阵[A,B,C,D]为子空间预估器模型：

线性离散时不变***形式为：

x_k+1＝Ax_k+Bu_k+Ke_k (1)

y_k＝Cx_k+Du_k+e_k (2)

其中u_k∈R^l为***的输入测量值，y_k∈R^m为***的输出测量值，x_k∈Rⁿ为***的过程状态，e_k∈R^m为平稳、零均值的白噪声新息序列，K为稳定的状态Kalman增益，A，B，C，D为相对应维数的AGV运动学***矩阵，A为n×n维***矩阵，B为n×r维控制矩阵，C为m×n维输出矩阵，D为m×r维前馈矩阵。

步骤2中，所述构构造当前时刻输入数据和输出数据的Hankel矩阵，求解当前时刻的AGV运动学***的预测输出数据是指：

设定采样时刻k∈{1,2,…,2i+j-1}，构造***输入u_k的i行j列Hankel矩阵：

其中，下标p和f分别代表“过去”和“未来”，也表示Hankel矩阵的维数；β为指数加权遗忘因子，满足0<β<1；同理，输出y_k的Hankel矩阵Y_p、Y_f由(1)-(2)迭代可以得到递推的预测输出：

其中，Γ_i∈R^im*n为广义可观测矩阵，和/>为下三角矩阵，分别表示：

Y_f的最优预测输出可由/>的行空间到/>的行空间的正交投影得到：

其中W_p为过去的输入输出数据矩阵，即L_w、L_u分别表示状态和确定输入子空间预估器矩阵，可通过如下的QR分解和SVD分解求得；

利用QR分解：

其中R为下三角阵，Q为正交矩阵。实施正交投影(9)，可得：

其中上角标表示Moore-Penrose广义逆，将得到的L_w和L_u用于预测控制器的设计；

考虑如下的预测控制目标代价函数：

其中N₂为预测时域，N_u为控制时域，λ为控制作用加权系数，r_t+k为未来时刻t+k的设定值；

的第1列用于预测未来输出值，采用增量型形式表示为：

作为AGV运动学***的预测输出数据。

步骤3中，所述根据AGV的控制动作，确定下一时刻的实际输出数据，计算***实际输出数据和预测输出数据的误差，若超过容许误差，则利用增加变遗忘因子的递推子空间辨识算法更新当前时刻的子空间预估器模型，得到下一时刻的子空间预估器模型，否则不更新当前时刻的子空间预估器模型是指：

根据AGV加入的控制动作u_k+1，确定k+1时刻的实际输出数据y_k+1，计算实际输出数据y_k+1和式(13)中预测输出数据的误差，若该误差超过容许误差，则利用变遗忘因子的递推子空间辨识算法更新子空间预估器模型，得到AGV运动学***矩阵[A,B,C,D]作为下一时刻的子空间预估器模型：

设定在t+1时刻，新的输入输出对{u(t+1),y(t+1)}已知，构造下述数据向量：

则当t+1时刻获取新的数据后，更新的数据Hankel矩阵为：

对更新的数据矩阵做QR分解，并通过一系列Givens旋转变换，对分解得到的下三角阵R进行SVD分解即可从求解获得***矩阵；

由(6)可知，***矩阵C的估计值可由下式直接获得：

***矩阵A则可通过求解下述线性方程获得：

实际上，的估计值则可描述为下述优化问题的解：

其中因此***矩阵A可通过递推最小二乘法求解，其递推算法如下：

其中γ(t)为随时间变化的遗忘因子，且0<γ(t)<1，常取接近1的常数；P(t)为逆相关矩阵；

***矩阵B和D包含在矩阵中，为得到其估计值，须从式(5)中消去/>项；

定义矩阵为子空间/>的正交补投影因子：

式(5)两边同时左乘右乘/>得：

其中

M_l∈R；即有：

在每一时刻用最近的估计值代替则/>和/>的估计值可以类似于式(26)～(28)递推求解。

步骤4中，AGV运动学***的等效逆模型则为实时更新的AGV运动学***矩阵[A,B,C,D]。

与现有技术相比，本发明具有如下优点与有益效果：

1、本发明基于递推子空间辨识的自适应逆控制AGV转速控制方法，该AGV转速控制方法能很好适应AGV参数变化的调速***设计，有效克服控制对象的非线性，并有效抑制外部扰动，有利于改善对AGV路径跟踪控制的品质和精度。

2、本发明基于递推子空间辨识的自适应逆控制AGV转速控制方法能够为建立高精度的AGV运动学***模型打下坚实基础，为AGV运动学***的控制***设计提供依据和参考，工程实用性较高。

附图说明

图1是本发明基于递推子空间辨识的自适应逆控制AGV转速控制方法的流程图；

图2是本发明AGV转速自适应逆控制器的示意图；

图3是本发明AGV转速自适应逆控制器的控制原理图；

图4(1)和图4(2)是本发明利用递推子空间辨识算法辨识AGV运动学***的参数a₁、a₂仿真效果图；

图5(1)和图5(2)是本发明利用递推子空间辨识算法辨识AGV运动学***的参数b₀、b₁仿真效果图；

图6是本实施例AGV采用的模型图。

具体实施方式

下面结合附图与具体实施方式对本发明作进一步详细的描述。

实施例

如图1至图6所示，本发明基于递推子空间辨识的自适应逆控制AGV转速控制方法是这样的：利用递推子空间辨识算法获得AGV运动学***的等效逆模型的基本参数；在此基础上利用自适应逆控制法进行AGV转速自适应逆控制器的设计，并将AGV运动学***的等效逆模型作为AGV转速自适应逆控制器的传递函数，实现对AGV转速的自适应跟踪控制。

具体包括以下步骤：

步骤1，采集待辨识AGV运动学***的输入数据和输出数据，在AGV监控平台采集至少100组实验数据，其中，采用AGV的位置作为AGV运动学***的输入数据，采用AGV轮子的速度作为AGV运动学***的输出数据。

步骤2，利用递推子空间辨识算法，得到子空间预估器模型，构造当前时刻输入数据和输出数据的Hankel矩阵，求解当前时刻的AGV运动学***的预测输出数据。

采用递推子空间辨识算法在线获得AGV运动学***的输入数据u和输出数据y，实现对AGV运动学***矩阵[A,B,C,D]的在线递推估计，AGV运动学***矩阵[A,B,C,D]为子空间预估器模型：

线性离散时不变***形式为：

x_k+1＝Ax_k+Bu_k+Ke_k (1)

y_k＝Cx_k+Du_k+e_k (2)

其中u_k∈R^l为***的输入测量值，y_k∈R^m为***的输出测量值，x_k∈Rⁿ为***的过程状态，e_k∈R^m为平稳、零均值的白噪声新息序列，K为稳定的状态Kalman增益，A，B，C，D为相对应维数的AGV运动学***矩阵，A为n×n维***矩阵，B为n×r维控制矩阵，C为m×n维输出矩阵，D为m×r维前馈矩阵。

为保证***的可辨识性以及收敛性分析要求，假定***满足下述条件：

1)(A，B)能控，(A，C)能观测，***为最小实现。

2)外部输入u与噪声e不相关，且u满足充分激励条件。

设定采样时刻k∈{1,2,…,2i+j-1}，构造***输入u_k的i行j列Hankel矩阵：

其中，下标p和f分别代表“过去”和“未来”，也表示Hankel矩阵的维数；β为指数加权遗忘因子，满足0<β<1；同理，输出y_k、噪声e_k的Hankel矩阵Y_p、Y_f、E_p、E_f也有类似定义。输出y_k的Hankel矩阵Y_p、Y_f由(1)-(2)迭代可以得到递推的预测输出：

其中，Γ_i∈R^im*n为广义可观测矩阵，和/>为下三角矩阵，分别表示：

Y_f的最优预测输出可由/>的行空间到/>的行空间的正交投影得到：

其中W_p为过去的输入输出数据矩阵，即L_w、L_u分别表示状态和确定输入子空间预估器矩阵，可通过如下的QR分解和SVD分解求得；

利用QR分解：

其中R为下三角阵，Q为正交矩阵。实施正交投影(9)，可得：

其中上角标表示Moore-Penrose广义逆，将得到的L_w和L_u用于预测控制器的设计；

考虑如下的预测控制目标代价函数：

其中N₂为预测时域，N_u为控制时域，λ为控制作用加权系数，r_t+k为未来时刻t+k的设定值；

的第1列用于预测未来输出值，采用增量型形式表示为：

作为AGV运动学***的预测输出数据。

其中

F＝[I_m … I_m]^T (15)

由L_w构造：

其中1≤k≤N₂。将(13)代入(12)，求取控制律，得到

只需将Δu_f的第一个值Δu_t作为控制输入，下一时刻再重新计算控制输入，因此t时刻的控制输入为：

u_t＝u_t-1+Δu_t (19)。

步骤3，根据AGV的控制动作，确定下一时刻的实际输出数据，计算***实际输出数据和预测输出数据的误差，若超过容许误差，则利用增加变遗忘因子的递推子空间辨识算法更新当前时刻的子空间预估器模型，得到下一时刻的子空间预估器模型，否则不更新当前时刻的子空间预估器模型：

根据AGV加入的控制动作u_k+1，确定k+1时刻的实际输出数据y_k+1，计算实际输出数据y_k+1和式(13)中预测输出数据的误差，若该误差超过容许误差，则利用变遗忘因子的递推子空间辨识算法更新子空间预估器模型，得到AGV运动学***矩阵[A,B,C,D]作为下一时刻的子空间预估器模型：

设定在t+1时刻，新的输入输出对{u(t+1),y(t+1)}已知，构造下述数据向量：

则当t+1时刻获取新的数据后，更新的数据Hankel矩阵为：

对更新的数据矩阵做QR分解，并通过一系列Givens旋转变换，对分解得到的下三角阵R进行SVD分解即可从求解获得***矩阵；

由(6)可知，***矩阵C的估计值可由下式直接获得：

***矩阵A则可通过求解下述线性方程获得：

实际上，的估计值则可描述为下述优化问题的解：

其中因此***矩阵A可通过递推最小二乘法求解，其递推算法如下：

其中γ(t)为随时间变化的遗忘因子，且0<γ(t)<1，常取接近1的常数；P(t)为逆相关矩阵；

***矩阵B和D包含在矩阵中，为得到其估计值，须从式(5)中消去/>项；定义矩阵/>为子空间/>的正交补投影因子：

式(5)两边同时左乘右乘/>得：

其中

M_l∈R；即有：

在每一时刻用最近的估计值代替则/>和/>的估计值可以类似于式(26)～(28)递推求解。

步骤4，重复步骤2和3，实时更新子空间预估器模型，确定AGV运动学***的等效逆模型；AGV运动学***的等效逆模型则为实时更新的AGV运动学***矩阵[A,B,C,D]。

步骤5，采用自适应逆控制法设计AGV转速自适应逆控制器，将AGV转速作为目标信号，将步骤4得到的AGV运动学***的等效逆模型作为控制器的传递函数；将目标信号输入AGV转速自适应逆控制器得到控制信号，再用控制信号激励AGV，实现对AGV转速的自适应跟踪控制。

根据现代控制理论，***的传递函数阵和状态空间表达式之间的转换公式G)s)＝C(sI-A)^-1B+D，可通过AGV运动学***状态空间表达式得到AGV运动学***的传递函数：

本实施例对本发明方法得到的AGV运动学***进行***辨识仿真实验，原AGV运动学***的参数a₁＝-1.971，a₂＝0.971，b₀＝0.6007，b₁＝0.6126。在k＝30时刻，AGV运动学***参数变化为a₁＝-0.97，a₂＝0.6702，b₀＝0.195，b₁＝2.0187，验证递推子空间辨识算法效果佳。结果表明递推子空间辨识算法在计算的50步内能够准确辨识出***参数，证实了辨识算法的有效性。

本实施例采用的AGV的模型如图6所示。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (2)

1.一种基于递推子空间辨识的自适应逆控制AGV转速控制方法，其特征在于：利用递推子空间辨识算法获得AGV运动学***的等效逆模型的基本参数；在此基础上利用自适应逆控制法进行AGV转速自适应逆控制器的设计，并将AGV运动学***的等效逆模型作为AGV转速自适应逆控制器的传递函数，实现对AGV转速的自适应跟踪控制；

包括以下步骤：

步骤1，采集待辨识AGV运动学***的输入数据和输出数据；

步骤2，利用递推子空间辨识算法，得到子空间预估器模型，构造当前时刻输入数据和输出数据的Hankel矩阵，求解当前时刻的AGV运动学***的预测输出数据；

步骤3，根据AGV的控制动作，确定下一时刻的实际输出数据，计算***实际输出数据和预测输出数据的误差，若超过容许误差，则利用增加变遗忘因子的递推子空间辨识算法更新当前时刻的子空间预估器模型，得到下一时刻的子空间预估器模型，否则不更新当前时刻的子空间预估器模型；

步骤4，重复步骤2和3，实时更新子空间预估器模型，确定AGV运动学***的等效逆模型；

步骤5，采用自适应逆控制法设计AGV转速自适应逆控制器，将AGV转速作为目标信号，将步骤4得到的AGV运动学***的等效逆模型作为控制器的传递函数；将目标信号输入AGV转速自适应逆控制器得到控制信号，再用控制信号激励AGV，实现对AGV转速的自适应跟踪控制；

步骤2中，所述利用递推子空间辨识算法，得到子空间预估器模型是指：采用递推子空间辨识算法在线获得AGV运动学***的输入数据u和输出数据y，实现对AGV运动学***矩阵[A,B,C,D]的在线递推估计，AGV运动学***矩阵[A,B,C,D]为子空间预估器模型：

线性离散时不变***形式为：

x_k+1＝Ax_k+Bu_k+Ke_k(1)

y_k＝Cx_k+Du_k+e_k(2)

其中u_k∈R^l为***的输入测量值，y_k∈R^m为***的输出测量值，x_k∈Rⁿ为***的过程状态，e_k∈R^m为平稳、零均值的白噪声新息序列，K为稳定的状态Kalman增益，A，B，C，D为相对应维数的AGV运动学***矩阵，A为n×n维***矩阵，B为n×r维控制矩阵，C为m×n维输出矩阵，D为m×r维前馈矩阵；

步骤2中，所述构造当前时刻输入数据和输出数据的Hankel矩阵，求解当前时刻的AGV运动学***的预测输出数据是指：

设定采样时刻k∈{1,2,…,2i+j-1}，构造***输入u_k的i行j列Hankel矩阵：

其中，下标p和f分别代表“过去”和“未来”，也表示Hankel矩阵的维数；β为指数加权遗忘因子，满足0<β<1；同理，输出y_k的Hankel矩阵Y_p、Y_f由(1)-(2)迭代可以得到递推的预测输出：

其中，Γ_i∈R^im*n为广义可观测矩阵，和/>为下三角矩阵，分别表示：

Y_f的最优预测输出可由/>的行空间到/>的行空间的正交投影得到：

其中W_p为过去的输入输出数据矩阵，即L_w、L_u分别表示状态和确定输入子空间预估器矩阵，可通过如下的QR分解和SVD分解求得；

利用QR分解：

其中R为下三角阵，Q为正交矩阵；实施正交投影(9)，可得：

其中上角标表示Moore-Penrose广义逆，将得到的L_w和L_u用于预测控制器的设计；

考虑如下的预测控制目标代价函数：

其中N₂为预测时域，N_u为控制时域，λ为控制作用加权系数，r_t+k为未来时刻t+k的设定值；

的第1列用于预测未来输出值，采用增量型形式表示为：

作为AGV运动学***的预测输出数据；

步骤3中，所述根据AGV的控制动作，确定下一时刻的实际输出数据，计算***实际输出数据和预测输出数据的误差，若超过容许误差，则利用增加变遗忘因子的递推子空间辨识算法更新当前时刻的子空间预估器模型，得到下一时刻的子空间预估器模型，否则不更新当前时刻的子空间预估器模型是指：

根据AGV加入的控制动作u_k+1，确定k+1时刻的实际输出数据y_k+1，计算实际输出数据y_k+1和式(13)中预测输出数据的误差，若该误差超过容许误差，则利用变遗忘因子的递推子空间辨识算法更新子空间预估器模型，得到AGV运动学***矩阵[A,B,C,D]作为下一时刻的子空间预估器模型：

设定在t+1时刻，新的输入输出对{u(t+1),y(t+1)}已知，构造下述数据向量：

则当t+1时刻获取新的数据后，更新的数据Hankel矩阵为：

对更新的数据矩阵做QR分解，并通过一系列Givens旋转变换，对分解得到的下三角阵R进行SVD分解即可从求解获得***矩阵；

由(6)可知，***矩阵C的估计值可由下式直接获得：

***矩阵A则可通过求解下述线性方程获得：

实际上，的估计值则可描述为下述优化问题的解：

其中因此***矩阵A可通过递推最小二乘法求解，其递推算法如下：

其中γ(t)为随时间变化的遗忘因子，且0<γ(t)<1，常取接近1的常数；P(t)为逆相关矩阵；

***矩阵B和D包含在矩阵中，为得到其估计值，须从式(5)中消去/>项；

定义矩阵为子空间/>的正交补投影因子：

式(5)两边同时左乘右乘/>得：

其中

M_l∈R；即有：

在每一时刻用最近的估计值代替则/>和/>的估计值可以类似于式(26)～(28)递推求解；

步骤4中，AGV运动学***的等效逆模型则为实时更新的AGV运动学***矩阵[A,B,C,D]。

2.根据权利要求1所述的基于递推子空间辨识的自适应逆控制AGV转速控制方法，其特征在于：步骤1中，采用AGV的位置作为AGV运动学***的输入数据，采用AGV轮子的速度作为AGV运动学***的输出数据。