CN109885807A - Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法 - Google Patents

Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法 Download PDF

Info

Publication number
CN109885807A
CN109885807A CN201910061995.5A CN201910061995A CN109885807A CN 109885807 A CN109885807 A CN 109885807A CN 201910061995 A CN201910061995 A CN 201910061995A CN 109885807 A CN109885807 A CN 109885807A
Authority
CN
China
Prior art keywords
latest estimated
correction term
weighting
white noise
estimation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201910061995.5A
Other languages
English (en)
Inventor
吴爱国
康瑞环
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Graduate School Harbin Institute of Technology
Original Assignee
Shenzhen Graduate School Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Graduate School Harbin Institute of Technology filed Critical Shenzhen Graduate School Harbin Institute of Technology
Priority to CN201910061995.5A priority Critical patent/CN109885807A/zh
Publication of CN109885807A publication Critical patent/CN109885807A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明提供了Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法。当干扰噪声为白噪声时,在最小二乘辨识方法的基础上引入权重系数对上一时刻和当前时刻的修正项取权重和。进而引入最新估计,用最新参数估计代替之前的参数估计。本发明的有益效果是:Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法的辨识精度高,收敛速度快,抗干扰能力强。

Description

Hammerstein系统在白噪声干扰下的加权最新估计最小二乘 辨识方法
技术领域
本发明涉及系统辨识领域,尤其涉及一种Hammerstein系统在白噪声干扰下的在线递推辨识方法。
背景技术
在实际系统中经常存在这样或是那样的干扰作用,比如测量误差,系统误差。干扰往往具有随机性,一些情况下可以为白噪声。Hammerstein系统是非线性环节与线性环节串联的形式。因为具有非线性特性,所以可以很好地用来描述放大功率器、挖掘臂、化工过程等。Hammerstein系统在白噪声干扰下的最小二乘辨识方法的精度还可进一步提高。
发明内容
为了解决现有技术中的问题,本发明提供了Hammerstein系统在白噪声干扰下的加权最小二乘辨识方法,针对白噪声干扰的情况,在最小二乘辨识方法的基础上引入权重系数,对上一时刻和当前时刻的修正项取权重和,进而引入最新估计,用最新参数估计代替之前的参数估计。
本发明通过以下技术方案来实现:
Hammerstein系统在白噪声干扰下的加权最小二乘辨识方法,根据系统的输入输出关系,利用最小二乘原理定义准则函数J(θ,t),并对其极小化,得到最小二乘辨识方法;然后将当前时刻参数估计定义为上一时刻的估计加上当前时刻修正项的和,采用加权思想,修正项改写为当前时刻修正项与上一时刻修正项的权重和:最后,采用最新估计的思想,对加权后的修正项用最新估计代替之前的估计。
定义准则函数
对准则函数极小化可以得到最小二乘算法:
定义修正项:
根据加权思想,对算法进行提升,递推形式改进为:
其中,户(t)为信息向量相关的变量,L(t)为中间变量,为待辨识参数向量的估计,α是权重参数,为信息向量,其定义为
作为本发明的进一步改进,在加权的基础上引入最新估计:
最终可以得到辨识方法
本发明的有益效果是:与现有技术相比,Hammerstein系统在白噪声干扰下,本发明有更好的辨识精度,收敛速度;在噪信比较高的情况下能表现出更好的抗干扰能力。
附图说明
图1是本发明的加权最新估计最小二乘辨识方法原理图。
图2是本发明的仿真实验图,其噪声方差为σ2=1.02
图3是本发明的仿真实验图,其噪声方差为σ2=1.22
具体实施方式
结合附图说明及具体实施方式对本发明作进一步说明。
如图1所示,本发明的Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法,包括:根据系统的输入输出关系,利用最小二乘原理定义准则函数J(θ,t):
并对J(θ,t)极小化得出最小二乘辨识方法:
其中,y(t)系统输出,户(t)为信息向量相关的变量,其引入的目的是避免递推公式中取逆,L(t)为中间变量无具体明确含义,为待辨识参数向量的估计。
当干扰噪声为白噪声,最小二乘辨识方法可以通过修正项加权和引入最新估计的方法提高辨识精度。
Step 1引入权重参数
定义修正项:
引入加权参数α,对修正项加权
Step 2引入最新估计
的精度优于因此引入最新估计
最终辨识算法为:
下面本发明通过实验仿真针对Hammerstein方程误差模型在白噪声干扰下加权最新估计辨识方法的有效性。
Hammertein方程误差的输入输出关系为:
其中,{u(t)}是持续激励信号序列,{y(t)}是可测的系统输出序列,{v(t)}是不可测的白噪声序列。f(·)是阶次已确定的非线性基的的线性函数,其系数为z-1是单位后移算子,且
系统输入输出关系:
定义信息向量:
其中,
其中对应参数
A(z)=1-1.60z-1+0.80z-2
B(z)=0.85z-1+0.65z-2
组合参数θ真实值
θ=[-1.60 0.80 0.85 0.65 0.425 0.325 0.2125 0.2]
在实验中用相对误差百分比δ来衡量模型的辨识精度,||·||2表示向量的2范数
此外,用噪信比NSR来衡量系统的抗干扰能力
{u(t)}是持续激励信号序列,其均值为0,方差为σ2=1.02。{v(t)}为白噪声序列,为验证抗干扰能力采用两种噪声强度分别为方差为σ2=1.02和σ2=1.22。其中RLS表示最小二乘算法,WRLS表示加权最新估计最小二乘辨识算法。
从图3中可以看出本发明提出的算法与原算法相比在上述模型中有更好的辨识精度和收敛速度,并且与图2对比可以看出当噪信比增大的时候,本发明的算法的精度下降更少。即本发明的提算法具有更好的抗干扰能力。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (3)

1.Hammerstein系统在白噪声干扰下的加权最小二乘辨识方法,其特征在于:所述方法根据系统的输入输出关系,利用最小二乘原理定义准则函数J(θ,t),并对其极小化,得到最小二乘辨识方法;然后将当前时刻参数估计定义为上一时刻的估计加上当前时刻修正项的和,采用加权思想,修正项改写为当前时刻修正项与上一时刻修正项的权重和;最后,采用最新估计的思想,对加权后的修正项用最新估计代替之前的估计。
2.根据权利要求1所述的方法,其特征在于:定义准则函数
对准则函数极小化可以得到最小二乘算法:
定义修正项:
根据加权思想,对算法进行提升,递推形式改进为:
其中,P(t)为信息向量相关的变量,L(t)为中间变量,为待辨识参数向量的估计,α是权重参数,为信息向量,其定义为
3.根据权利要求1所述的方法,其特征在于:对加权后的修正项用最新估计代替之前的估计具体为:根据最新估计思想,在求解的时候作为最新估计是已知,在一般情况下更接近真实值,因此将辨识方法改进为:
CN201910061995.5A 2019-01-22 2019-01-22 Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法 Pending CN109885807A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910061995.5A CN109885807A (zh) 2019-01-22 2019-01-22 Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910061995.5A CN109885807A (zh) 2019-01-22 2019-01-22 Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法

Publications (1)

Publication Number Publication Date
CN109885807A true CN109885807A (zh) 2019-06-14

Family

ID=66926615

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910061995.5A Pending CN109885807A (zh) 2019-01-22 2019-01-22 Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法

Country Status (1)

Country Link
CN (1) CN109885807A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111427266A (zh) * 2020-03-20 2020-07-17 北华航天工业学院 一种针对带扰动的非线性系统辨识方法
CN111478316A (zh) * 2020-04-01 2020-07-31 华南理工大学 用于电力系统灵敏度在线辨识的噪声自助集成回归方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111427266A (zh) * 2020-03-20 2020-07-17 北华航天工业学院 一种针对带扰动的非线性系统辨识方法
CN111478316A (zh) * 2020-04-01 2020-07-31 华南理工大学 用于电力系统灵敏度在线辨识的噪声自助集成回归方法

Similar Documents

Publication Publication Date Title
Singh et al. Robust estimation of power system harmonics using a hybrid firefly based recursive least square algorithm
CN110579638B (zh) 基于卡尔曼滤波的交叉供电系统动态电压跌落检测方法
CN109885807A (zh) Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法
CN113765880B (zh) 一种基于时空关联性的电力系统网络攻击检测方法
Wang et al. A class of improved least sum of exponentials algorithms
CN109753762A (zh) 基于类别修正的配电网两阶段网络拓扑识别方法及装置
CN110059294A (zh) 基于核递归最大互相关熵算法的非线性格兰杰因果性检测方法
CN104979829A (zh) 一种基于改进的FastICA算法的谐波电流估计方法
CN108281961B (zh) 一种自适应鲁棒扩展卡尔曼的参数辨识方法
Kumar et al. Efficient detection of false data injection attacks on AC state estimation in smart grids
CN103699650A (zh) 消息传播预测方法及装置
CN109918617A (zh) Hammerstein系统在有色噪声下的加权最新估计最小二乘辨识方法
Diversi Bias‐eliminating least‐squares identification of errors‐in‐variables models with mutually correlated noises
CN106452621B (zh) 一种峰值功率、峰均值功率比的确定方法及装置
Moravej et al. A hybrid least squares–clonal selection based algorithm for harmonics estimation
Zheng On estimation of autoregressive signals in the presence of noise
Liang et al. Adaptive filtering for stochastic systems with generalized disturbance inputs
Terespolsky et al. Developing an approximation to the heidler function-with an analytical transformation into the frequency domain
Kulikov et al. The accurate continuous-discrete extended Kalman filter for continuous-time stochastic systems
Cui et al. An adaptive improved unsaturated bistable stochastic resonance method based on weak signal detection
CN115409245A (zh) 电力系统的预测辅助状态估计方法、装置、设备及介质
CN115047400A (zh) 三相电能表准确度校验方法、系统、终端设备及介质
KR100926749B1 (ko) 2차 차분 전력분석 공격 방법 및 2차 차분 전자기분석 공격방법
Savkin et al. Weak robust controllability and observability of uncertain linear systems
Djurović et al. Estimation of single-tone signal frequency by using the L-DFT

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20190614