CN109885807A - Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法 - Google Patents
Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012937 correction Methods 0.000 claims abstract description 18
- 238000012850 discrimination method Methods 0.000 claims description 9
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 description 4
- 108010076504 Protein Sorting Signals Proteins 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000002085 persistent effect Effects 0.000 description 2
- 238000001311 chemical methods and process Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供了Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法。当干扰噪声为白噪声时,在最小二乘辨识方法的基础上引入权重系数对上一时刻和当前时刻的修正项取权重和。进而引入最新估计,用最新参数估计代替之前的参数估计。本发明的有益效果是:Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法的辨识精度高,收敛速度快,抗干扰能力强。
Description
技术领域
本发明涉及系统辨识领域,尤其涉及一种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所述的方法,其特征在于:对加权后的修正项用最新估计代替之前的估计具体为:根据最新估计思想,在求解的时候作为最新估计是已知,在一般情况下比更接近真实值,因此将辨识方法改进为:
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111427266A (zh) * | 2020-03-20 | 2020-07-17 | 北华航天工业学院 | 一种针对带扰动的非线性系统辨识方法 |
CN111478316A (zh) * | 2020-04-01 | 2020-07-31 | 华南理工大学 | 用于电力系统灵敏度在线辨识的噪声自助集成回归方法 |
-
2019
- 2019-01-22 CN CN201910061995.5A patent/CN109885807A/zh active Pending
Cited By (2)
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 |