【系统辨识】随机系统的数学模型


前言

数学模型有确定性和随机性之分,模型中变量受到随机干扰的系统称为随机系统,反之为确定系统。本文会介绍随机系统的各个数学模型


一、随机系统模型

随机系统模型是指在确定系统模型的基础上加上一个随机噪声,因此数学模型中带有噪声项的都可以认为是随机系统模型。
首先将随机系统分为了“时间序列模型”、“方程误差类模型”和“输出误差类模型”三大类,接下来从这三个类别依次讨论。以下讨论的系统均为单输入单输出系统。

为了方便起见,设 u ( t ) {u(t)} u(t)为系统输入序列, y ( t ) {y(t)} y(t)为系统观测输出序列, v ( t ) {v(t)} v(t)为零均值随机白噪声序列,然后定义:
A ( z ) = 1 + a 1 z − 1 + a 2 z − 2 + … + a n a z − n a , a i ∈ R , A(z)=1+a_1z^{-1}+a_2z^{-2}+…+a_{n_a}z^{-n_a},a_i\in R, A(z)=1+a1z1+a2z2++anazna,aiR,
B ( z ) = b 1 z − 1 + b 2 z − 2 + … + b n b z − n b , b i ∈ R , B(z)=b_1z^{-1}+b_2z^{-2}+…+b_{n_b}z^{-n_b},b_i\in R, B(z)=b1z1+b2z2++bnbznb,biR,
C ( z ) = 1 + c 1 z − 1 + c 2 z − 2 + … + c n c z − n c , c i ∈ R , C(z)=1+c_1z^{-1}+c_2z^{-2}+…+c_{n_c}z^{-n_c},c_i\in R, C(z)=1+c1z1+c2z2++cncznc,ciR,
D ( z ) = 1 + d 1 z − 1 + d 2 z − 2 + … + d n d z − n d , d i ∈ R , D(z)=1+d_1z^{-1}+d_2z^{-2}+…+d_{n_d}z^{-n_d},d_i\in R, D(z)=1+d1z1+d2z2++dndznd,diR,
F ( z ) = 1 + f 1 z − 1 + f 2 z − 2 + … + f n f z − n f , f i ∈ R , F(z)=1+f_1z^{-1}+f_2z^{-2}+…+f_{n_f}z^{-n_f},f_i\in R, F(z)=1+f1z1+f2z2++fnfznf,fiR,

1.时间序列模型

时间序列模型的特征是系统中有两个变量,一个是观测 y ( t ) y(t) y(t),一个是随机白噪声 v ( t ) v(t) v(t).

1.1 自回归模型(AR模型)

A ( z ) y ( t ) = v ( t ) A(z)y(t)=v(t) A(z)y(t)=v(t)

1.2 滑动平均模型(MA模型)

y ( t ) = D ( z ) v ( t ) y(t)=D(z)v(t) y(t)=D(z)v(t)

1.3 自回归滑动平均模型(ARMA模型)

A ( z ) y ( t ) = D ( z ) v ( t ) A(z)y(t)=D(z)v(t) A(z)y(t)=D(z)v(t)

1.4 确定性ARMA模型

A ( z ) y ( t ) = D ( z ) u ( t ) A(z)y(t)=D(z)u(t) A(z)y(t)=D(z)u(t)
这里是一个确定系统模型,因为随机的 v ( t ) v(t) v(t)被替换为了确定的 u ( t ) u(t) u(t).

1.5 带积分ARMA模型

A ( z ) ( 1 − z − 1 ) d y ( t ) = D ( z ) v ( t ) A(z)(1-z^{-1})^dy(t)=D(z)v(t) A(z)(1z1)dy(t)=D(z)v(t)

2.方程误差类模型

A ( z ) y ( t ) = B ( z ) u ( t ) + ω ( t ) A(z)y(t)=B(z)u(t)+\omega(t) A(z)y(t)=B(z)u(t)+ω(t)
其中 ω ( t ) \omega(t) ω(t)为白噪声或有色噪声。
对于方程误差类模型的分析,需要把重点放在 ω ( t ) \omega(t) ω(t)上。

2.1 受控自回归模型(CAR模型)

A ( z ) y ( t ) = B ( z ) u ( t ) + v ( t ) A(z)y(t)=B(z)u(t)+v(t) A(z)y(t)=B(z)u(t)+v(t)

2.2 受控自回归滑动平均模型(CARMA模型)

A ( z ) y ( t ) = B ( z ) u ( t ) + D ( z ) v ( t ) A(z)y(t)=B(z)u(t)+D(z)v(t) A(z)y(t)=B(z)u(t)+D(z)v(t)

2.3 受控自回归自回归模型(CARAR模型)

A ( z ) y ( t ) = B ( z ) u ( t ) + 1 C ( z ) v ( t ) A(z)y(t)=B(z)u(t)+\frac{1}{C(z)}v(t) A(z)y(t)=B(z)u(t)+C(z)1v(t)

2.4 受控自回归自回归滑动平均模型(CARARMA模型)

A ( z ) y ( t ) = B ( z ) u ( t ) + D ( z ) C ( z ) v ( t ) A(z)y(t)=B(z)u(t)+\frac{D(z)}{C(z)}v(t) A(z)y(t)=B(z)u(t)+C(z)D(z)v(t)

3.输出误差类模型

y ( t ) = B ( z ) A ( z ) u ( t ) + ω ( t ) y(t)=\frac{B(z)}{A(z)}u(t)+\omega(t) y(t)=A(z)B(z)u(t)+ω(t)
模型特征是包含了有理方式项 B ( z ) A ( z ) u ( t ) \frac{B(z)}{A(z)}u(t) A(z)B(z)u(t)

3.1 输出误差模型(OE模型)

y ( t ) = B ( z ) A ( z ) u ( t ) + v ( t ) y(t)=\frac{B(z)}{A(z)}u(t)+v(t) y(t)=A(z)B(z)u(t)+v(t)

3.2 输出误差滑动平均模型(OMEA模型)

y ( t ) = B ( z ) A ( z ) u ( t ) + D ( z ) v ( t ) y(t)=\frac{B(z)}{A(z)}u(t)+D(z)v(t) y(t)=A(z)B(z)u(t)+D(z)v(t)

3.3 输出误差自回归模型

y ( t ) = B ( z ) A ( z ) u ( t ) + 1 C ( z ) v ( t ) y(t)=\frac{B(z)}{A(z)}u(t)+\frac{1}{C(z)}v(t) y(t)=A(z)B(z)u(t)+C(z)1v(t)

3.4 输出误差自回归滑动平均模型(OEARMA模型或Box-Jenkins模型)

y ( t ) = B ( z ) A ( z ) u ( t ) + D ( z ) C ( z ) v ( t ) y(t)=\frac{B(z)}{A(z)}u(t)+\frac{D(z)}{C(z)}v(t) y(t)=A(z)B(z)u(t)+C(z)D(z)v(t)

3.5 一般输出误差模型

F ( z ) y ( t ) = B ( z ) A ( z ) u ( t ) + D ( z ) C ( z ) v ( t ) F(z)y(t)=\frac{B(z)}{A(z)}u(t)+\frac{D(z)}{C(z)}v(t) F(z)y(t)=A(z)B(z)u(t)+C(z)D(z)v(t)


参考文献

丁锋.系统辨识新论[M].北京:科学出版社,2013:45.

  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

麦斯威尔逊

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值