用相关法辨识系统的脉冲响应 matlab,利用相关分析法辨识脉冲响应

利用相关分析法辨识脉冲响应

自1205 刘彬 41251141

1 实验方案设计

1.1 生成输入数据和噪声

用M 序列作为辨识的输入信号，噪声采用标准正态分布的白噪声。 生成白噪声时，首先利用乘同余法生成U[0,1]均匀分布的随机数，再利用U[0,1]均匀分布的随机数生成标准正态分布的白噪声。 1.2 过程仿真

模拟过程传递函数)(s G ，获得输出数据y(k)。)(s G 采取串联传递函数仿真，

2

12111

11)(T s T s T T K s G ++=

，用M 序列作为辨识的输入信号。

1.3 计算互相关函数

∑++=-=

p

p N r N i p

Mz i z k i u rN k R )1(1

)()(1

)(

其中r 为周期数，1+=p N i 表示计算互相关函数所用的数据是从第二个周期开始的，目的是等过程仿真数据进入平稳状态。

1.4 计算脉冲响应估计值、脉冲响应理论值、脉冲响应估计误差

脉冲响应估计值[]

)1()()1()(?2

--?+=p Mz Mz

p

p

N R k R

t

a N N k g

脉冲响应理论值[]

21//2

10)(T t k T t k e e T T K

k g ?-?---=