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
- estimation
- correction term
- weighting
- latest
- white noise
- 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 28
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 3
- 108010076504 Protein Sorting Signals Proteins 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 238000006467 substitution reaction 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 |
|---|---|---|
| Billings et al. | Piecewise linear identification of non-linear systems | |
| CN105929340B (zh) | 一种基于arima估算电池soc的方法 | |
| CN109885807A (zh) | Hammerstein系统在白噪声干扰下的加权最新估计最小二乘辨识方法 | |
| CN109753762A (zh) | 基于类别修正的配电网两阶段网络拓扑识别方法及装置 | |
| CN108512528B (zh) | 一种cim函数下的比例控制和归一化lmp滤波方法 | |
| CN103684349B (zh) | 一种基于递推协方差阵估计的卡尔曼滤波方法 | |
| Sentoni et al. | State‐space nonlinear process modeling: Identification and universality | |
| Andrieu et al. | Continuous discrete observer with updated sampling period | |
| WO2004090782A1 (en) | Accurate linear parameter estimation with noisy inputs | |
| CN109918617A (zh) | Hammerstein系统在有色噪声下的加权最新估计最小二乘辨识方法 | |
| CN117630678A (zh) | 电动汽车电池电压异常波动故障在线诊断方法及系统 | |
| CN109274107B (zh) | 一种计及奇异值的低频振荡信号参数辨识方法 | |
| CN116055340A (zh) | 分布式网络未知参数估计方法、装置及电子设备 | |
| CN116227674A (zh) | 一种基于在线更新模型的碳排放预测方法 | |
| WO2019149376A1 (en) | Method and system for processing input data using a neural network and normalizations | |
| Savkin et al. | Weak robust controllability and observability of uncertain linear systems | |
| CN111274752B (zh) | 一种基于两级开环和二元函数的功率放大器行为建模方法 | |
| Voros | Recursive identification of Hammerstein systems with polynomial nonlinearities | |
| CN110556927A (zh) | 一种感应式无线电能传输系统的状态观测的方法及系统 | |
| CN1485798A (zh) | 一种神经网络均衡器的优化训练方法 | |
| Liu et al. | Modeling of mems gyroscope random error based on Kalman filter | |
| CN114566227B (zh) | 一种基于牛顿迭代算法的ph中和过程模型辨识方法 | |
| Jing | On parameter separation of the Wiener system identification | |
| Pandey et al. | Variational approach to joint linear model and state estimation | |
| CN113359429B (zh) | 一种在线识别间隙环节特征参数的方法及装置 |
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 |
|
| RJ01 | Rejection of invention patent application after publication |