没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
- 1 -
实验一 在线广义与增广最小二乘辨识器的设计与实现
一、实验原理
1、M 序列
伪随机二进制序列(Pseudo Random Binary Sequence, PRBS)是广泛应用的
一种序列,所谓“二进制”是批序列中每个随机变量只有 0 或 1 两种逻辑状态。
伪随机序列可由多级线性反馈移位寄存器组成的随机信号发生器产生,其中具有
最长循环周期的线性移位寄存器序列是伪随机二进制序列最常见的一种形式,简
称 M 序列。 M 序列由于具有近似白噪声的性质,能保证较好的系统辨识精度,
是普遍采用的一种新辨识用输入信号。具体产生方式如图 1 所示。
x
1
x
2
x
3
x
4
移位脉冲
M序列
图 1 M 序列产生的原理图
2、在线的增广最小二乘法
当噪声均值为 0 时,最小二乘参数估计计算法为无偏估计,当噪声的均值不
为 0 时,最小二乘参数估计计算法为有偏估计。为了解决最小二乘参数估计的有
偏性,将噪声模型的辨识同时考虑进去,因此称为增广最小二乘法。
增广最小二乘法(Extended Least Squares, ELS)是最小二乘法的一种简单推
广,它只是扩充了参数向量 θ 和数据向量 φ 的维度,在辨识过程中同时考虑了噪
声模型的参数。
设 SISO 系统采用的数学模型为:
( )
( )
( )
( )
( )
( )
1 1 1
A z y k B z u k D z k
− − −
=+
其中
( )
1 1 2
12
1
n
n
A z a z a z a z
− − − −
= + + + +
( )
1 1 2
0 1 2
n
n
B z b b z b z b z
− − − −
= + + + +
( )
1 1 2
0 1 2
m
m
D z d d z d z d z
− − − −
= + + + +
参数矩阵
1 2 0 1 0 1
T
n n m
a a a b b b d d d
=
输出方程
( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
1 2 0
10
12
1
n
nm
y k a y k a y k a y k n b u k
bu k b u k n d k d k m
= − − − − − − − + +
− + + − + + + −
令
- 2 -
( ) ( ) ( ) ( ) ( ) ( ) ( )
1
T
n N y n N y N u n N u N n N n N m
+ = − + − − + + + −
且有
( )
( )
( )
( )
( )
( )
( )
( )
1
1
2
2
,
T
T
T
n
yn
yn
n
N Y N
y n N
nN
+
+
+
+
==
+
+
沿用最小二乘的思路
( ) ( ) ( ) ( ) ( )
1
ˆ
TT
N N N N Y N
−
=
令
( )
( )
( )
( )
( )
( )
1 , 1
11
T
n N Y N
n N Y n N
n N y n N
+
+ + = + + =
+ + + +
,于是有:
( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
1
1
ˆ
1 1 1 1 1
1 1 1 1
TT
T T T T
N N N N Y N
N N n N n N N Y N n N y n N
−
−
+ = + + + +
= + + + + + + + + + +
定义协方差阵为
( ) ( ) ( )
1
T
P N N N
−
=
结合矩阵求逆引理
( )
1
1
A BC A
−
−
+ = −
( )
1
1 1 1
A B I CA B CA
−
− − −
+
有
( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( )
( ) ( ) ( )
( ) ( )
1
1
1
1
1 1 1
1 1 1
11
1
1
1 1 1
TT
T
T
T
T
P N N N n N n N
P N P N n N n N
I P N n N n N P N
P N n N
I n N P N
n N P N n N
−
−
−
−
+ = + + + + +
= + + + + +
= + + + + +
++
= − + +
+ + + + +
记增益矩阵
( )
( ) ( )
( ) ( ) ( )
1
1
1 1 1
T
P N n N
KN
n N P N n N
++
+=
+ + + + +
有
( ) ( ) ( ) ( )
ˆ
1 1 1 1N P N N Y N
+ = + + +
整理可得
( ) ( ) ( ) ( ) ( ) ( )
ˆ ˆ ˆ
1 1 1 1
T
N N K N y n N n N N
+ = + + + + − + +
( )
( ) ( )
( ) ( ) ( )
1
1
1 1 1
T
P N n N
KN
n N P N n N
++
+=
+ + + + +
( ) ( ) ( ) ( )
1 1 1
T
P N I K N n N P N
+ = − + + +
3、在线的广义最小二乘法
剩余7页未读,继续阅读
资源评论
hit-lu
- 粉丝: 2
- 资源: 4
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功