多变量差分方程模型的极大似然递推最小二乘辨识算法
【专利摘要】本发明公开了一种多变量差分方程模型的极大似然递推最小二乘辨识算法，包括根据现有极大似然原理，构建出一种多变量差分方程模型中子系统极大似然准则函数；以多变量差分方程模型中子系统最大似然准则函数为基础，构建出多变量差分方程模型的极大似然递推最小二乘辨识算法的实现流程；构建一套多变量差分方程模型的极大似然递推最小二乘辨识算法。本发明采用极大似然原理和递推辨识方法，应用于线性多变量系统的参数估计。
【专利说明】
多变量差分方程模型的极大似然递推最小二乘辨识算法
技术领域
[0001] 本发明涉及一种多变量差分方程模型的极大似然递推最小二乘辨识算法。
【背景技术】
[0002] 数学模型在控制领域及其他工程领域具有非常重要的地位，它用于描述系统变量 之间的相互关系。许多工业中的实际系统可以建模为多变量系统，它的特点在于系统中具 有多个输入输出变量。与单变量系统相比，由于多变量系统变量多，维数高，结构复杂，故多 变量系统的建模和辨识比单变量系统复杂。多变量系统可以用不同的数学模型来描述，比 如状态空间模型、传递函数模型等。本发明适用于多变量差分方程模型的参数辨识。
[0003] 在系统辨识和参数估计领域，极大似然辨识方法是一种非常有用的参数辨识方 法，它最初由英国统计学家Fi sher发展起来，是一种基于概率论的辨识方法。极大似然估计 方法的基本思想是以观测数据和待辨识的未知参数为自变量构建一个似然函数或对数似 然函数，通过极大化这个似然函数获得参数估计值。由于极大似然估计具有一致性、有效性 和渐近正态性等统计性质，受到了国内外专家和学者的广泛关注，在许多领域都有广泛的 应用。
[0004] 按照辨识算法的执行方式来划分，递推辨识和迭代辨识是系统辨识中非常重要的 两类辨识技术，与迭代辨识相比，递推辨识能够在线估计系统的参数，具有占用内存空间 小、计算量小、可以在线估计系统参数的优点。
【发明内容】
[0005] 本发明的目的在于提供一种利用极大似然原理和递推辨识方法，应用于线性多变 量系统参数估计的多变量差分方程模型的极大似然递推最小二乘辨识算法。
[0006] 本发明的技术解决方案是：
[0007] -种多变量差分方程模型的极大似然递推最小二乘辨识算法，其特征是:包括下 列步骤：
[0008] (1)根据现有极大似然原理，构建出一种多变量差分方程模型中子系统极大似然 准则函数：
[0010]上式符号说明也作为递推时刻信参数向量，作为信息向量；
[0011] (2)以多变量差分方程模型中子系统最大似然准则函数为基础，构建出多变量差 分方程模型的极大似然递推最小二乘辨识算法的实现流程：
[0012] 第一步：启动算法；
[0013] 第二步:对递推时刻t进行初始化，初始值为1;
[0014]第三步:采集输入-输出数据u(t)和y(t)，构造出信息向量 今
[0016] 第五步:计算出U(t)和Pi(t);
[0017] 第六步:计算奋#);
[0018] 第七步:实时刷新氣{|). ?
[0019] 第八步:递推时刻t加1，重复上述步骤；
[0020] 上述各符号的含义：
[0033] (3)结合步骤(2)流程构建一套多变量差分方程模型的极大似然递推最小二乘辨 识算法，如下：
[0046]对上述算法中符号的说明：
[0056]上述算法的具体步骤：
[0058] 13)采集输入-输出数据1!(〇和7&)，分别通过式（25)、（26)、（27)构造 和1(0通过式(23)构造參rW . .、、 s v <> .， ：，
[0059] c)分别通过式（28)、（29)、（30)计算
接着通过式 (24)构造私，(，?）；
[0060] d)分别通过式(20)、（21)、（22)计算1^(〇?4〇以及匕">;
[0061] e)通过式(19)刷新所估参数. '9
[0062] f )t值增加1，重复上述步骤。
[0063]本发明采用极大似然原理和递推辨识方法，应用于线性多变量系统的参数估计。
【附图说明】
[0064]下面结合附图和实施例对本发明作进一步说明。
[0065] 图1是本发明极大似然递推最小二乘辨识算法流程图。
【具体实施方式】
[0066] -种多变量差分方程模型的极大似然递推最小二乘辨识算法，包括下列步骤：
[0067] (1)根据现有极大似然原理，构建出一种多变量差分方程模型中子系统极大似然 准则函数：
[0069] 上式符号说明…作为递推时刻信参数向量，#<|:M乍为信息向量；
[0070] (2)以多变量差分方程模型中子系统最大似然准则函数为基础，构建出多变量差 分方程模型的极大似然递推最小二乘辨识算法的实现流程：
[0071] 第一步:启动算法；
[0072] 第二步:对递推时刻t进行初始化，初始值为1;
[0073] 第三步:采集输入-输出数据u⑴和y(t)，构造出信息向量?Μ?
[0074] 第四步:计算滤波信息向量参⑷以及构造出％./!.〇 ^
[0075] 第五步:计算出Li(t)和Pi(t);
[0076] 第六步:计算6");
[0077] 第七步:实时刷新表a :1 .
[0078]第八步:递推时刻t加1，重复上述步骤；
[0079]上述各符号的含义：
[0080] 输入变量:銶_飪滅# ..
[0092] (3)结合步骤(2)流程构建一套多变量差分方程模型的极大似然递推最小二乘辨 识算法，如下：
[0117] b)采集输入-输出数据u(t)和y(t)，分别通过式(25)、(26)、(27)构造 ··* .、 和K/ )通过式(23)构造麵貌》. .， 9
[one]幻分别通过式⑶以四以如计算人/热和接着通过式 (24)构造^/'(〇. <·、 9
[0119] d)分别通过式(20)、（21)、（22)计算U(t)、Pi(t)以及. 9.
[0120] e)通过式(19)刷新所估参数? s < Π . 9
[0121] f )t值增加1，重复上述步骤。
【主权项】
1. 一种多变量差分方程模型的极大似然递推最小二乘辨识算法，其特征是：包括下列 步骤： (1) 根据现有极大似然原理，构建出一种多变量差分方程模型中子系统极大似然准则 函数：上式符号说明：0,作为递推时刻t的参数向Μ乍为信息向量； (2) 以多变量差分方程模型中子系统最大似然准则函数为基础，构建出多变量差分方 程模型的极大似然递推最小二乘辨识算法的实现流程： 第一步:启动算法； 第二步:对递推时刻t进行初始化，初始值为1; 第三步:采集输入-输出数据u (t)和y (t)，构造出信息向量. 第四步:计算滤波信息向量，叙以及备私/_)，构造出參?/以)； 第五步:计算出Li(t)和Pi(t); 第六步:计算％ ; 第七步:实时刷新#邊:）. 9 第八步:递推时刻t加1，重复上述步骤； 上述各符号的含义： 输入变量：妓以）€捷、（， 输出变量:. ，： 子系统参数向量：子系统参数向量義€ ν ， 子系统息向量：子系统信息向量：子系统滤波信息向量：子系统滤波信息向量、、' ?· 协方差矩阵： ' ， 增益向量(3)结合步骤(2)流程构建一套多变量差分方程模型的极大似然递推最小二乘辨识算 法，如下： s·· ^ S ^ 、 、、 ?· f >< \<>\ -· ' i > < ··；·· V > ? > \ Λ. < f K t ·· 、'.：，. '·； ·· 对上述算法中符号的说明： 定义输入变量为(0 € 输出变量为，. ? 9 定义卿丨€臟.^力人% €.1ΙΚ?@作为子系统参数向量； 定义#| € Ρ:(2内为子系统参数向量； 定义€'$/(0 € I1 为子系统信息向量； 定义作为子系统信息向量； 定义llh f(i) € Slx巧八 ill f _ € llx% € E1 为子系统滤 波信息向量； 定义辦f(/) € 为子系统滤波信息向量； _卩为协方差矩阵； 徽_-作为增益向量； 上述算法的具体步骤： a) 令t = l，设置初始b) 采集输入-输出数据u(t)和y(t)，分别通过式(25)、（26)、（27)构造、七;(?) 和^⑴，通过式(23)构造参|的； c) 分别通过式(28)、（29)、（30)计算^接着通过式(24)构 造 d) 分别通过式(20)、（21)、（22)计算 UUhPKt)以及 e) 通过式(19)刷新所估参数# 5 0七值增加1，重复上述步骤。
【文档编号】G06F17/17GK105868163SQ201610173971
【公开日】2016年8月17日
【申请日】2016年3月24日
【发明人】李俊红, 杨奕, 朱建红, 李晨, 杨赛, 张晴, 李建国
【申请人】南通大学