系统建模与仿真-哈尔滨工业大学

《系统辨识》

实验手册

哈尔滨工业大学控制与仿真中心

2018年5月

目录

实验1 白噪声和M序列的产生---------------------------------------------------------- 2 实验2 脉冲响应法的实现---------------------------------------------------------------- 5 实验3 递推最小二乘法的实现---------------------------------------------------------- 9 附录实验报告模板---------------------------------------------------------------------- 13

实验1 白噪声、M 序列的产生

一、实验目的

1、熟悉并掌握产生均匀分布随机序列方法以及进而产生高斯白噪声方法

2、熟悉并掌握M 序列生成原理及仿真生成方法

二、实验原理

1、混合同余法

混合同余法是加同余法和乘同余法的混合形式，其迭代式如下：

11

1(*)mod /n n n n x a x b M

R x M +++=+⎧⎨

=⎩ 式中a 为乘子，0x 为种子，b 为常数，M 为模。混合同余法是一种递归算法，即先提供一个种子0x ，逐次递归即得到一个不超过模M 的整数数列。

2、正态分布随机数产生方法

由独立同分布中心极限定理有：设随机变量12,,....,,...n X X X 相互独立，服从同一分布，且具有数学期望和方差：

2(),()0,(1,2,...)k k E X D X k μσ==>=

则随机变量之和1

n

k i X =∑的标准化变量

:

()

n

n

n

k

k k

X

E X X

n Y μ

--=

=

∑∑∑近似服从(0,1)N 分布。

如果n X 服从[0, 1]均匀分布，则上式中0.5μ=，2

1

12

σ=

。即

0.5n

k

X

n Y -=

∑近似服从(0,1)N 分布。

3、M 序列生成原理

用移位寄存器产生M 序列的简化框图如下图所示。该图表示一个由4个双稳态触发器顺序连接而成的4级移位寄存器，它带有一个反馈通道。当移位脉冲来到时，每级触发器的状态移到下一级触发器中，而反馈通道按模2加法规则反馈到第一级的输入端。

三、实验内容

1、生成均匀分布随机序列

（1）利用混合同余法生成[0, 1]区间上符合均匀分布的随机序列，并计算该序列的均值和方差，与理论值进行对比分析。要求序列长度为1200，推荐参数为a=65539，M=2147483647，0

（2）将[0, 1]区间分为不重叠的等长的10个子区间，绘制该随机序列落在每个子区间的频率曲线图，辅助验证该序列的均匀性。

（3）对上述随机序列进行独立性检验。（该部分为选作内容）

2、生成高斯白噪声

利用上一步产生的均匀分布随机序列，令n=12，生成服从N(0,1)的白噪声，序列长度为100，并绘制曲线。

3、生成M 序列

M 序列的循环周期取为63126=-=P N ，时钟节拍Sec 1=∆t ，幅度1=a ，逻辑“0”为a ，逻辑“1”为-a ，特征多项式65()F s s s =⊕。

生成M 序列的结构图如下所示。

要求编写Matlab程序生成该M序列，绘制该信号曲线，并分析验证M序列的性质。

四、实验步骤

1．分别画出三部分实验内容的程序框图（流程图）；

2．编制MATLAB的M文件；

3．运行编制的M文件；

4．查看程序运行结果并进行分析；

5．填写实验报告。

五、实验报告

格式参见附录一。

实验2 相关分析法辨识脉冲响应

一、实验目的

通过仿真实验掌握利用相关分析法辨识脉冲响应的原理和方法。

二、实验原理

一个单入单出线性定常系统的动态特性可用它的脉冲响应函数g(σ)来描述。

这样，只要记录x(t)、y(t)的值，并计算它们的互相关函数，即可求得脉冲响应函数g(τ)。

而在系统有正常输入的情形下，辨识脉冲响应的原理图如下图所示。

0 ()()()y t g x t d σσσ

∞

=-⎰则000

()11lim ()()(){lim ()()}T T

T T x t y t x t dt g x t x t dt d T T ττσστσ

∞→∞→∞--=--⎰⎰⎰上式两端同乘，进而取时间均值，有

0 ()()()xy x R g R d τστσσ

∞

=--⎰则这就是著名的维纳霍夫积分方程。0

() ()(), ()() ()()()()

()()x x xy x xy x t R k R k R g R d kg R g k

τδττσδτστστσστττ∞=-=--=-==⎰如果输入是，这时的自相关函数为则根据维纳霍夫积分方程可得

或者 白噪声

三、实验内容

下图为本实验的原理框图。系统的传递函数为)(s G ,其中

S e c 26T S e c ,3812021..,===T K ；)()(k z k u 和分别为系统的输入和输出变量；)(k v 为测量白噪声，服从正态分布，均值为零，方差为2v σ，记作),(~)(20v N k v σ；

)(k g 0为系统的脉冲响应理论值，)(ˆk g 为系统脉冲响应估计值，)(~k g 为系统脉冲响应估计误差。

系统的输入采用M 序列（采用实验1中的M 序列即可），输出受到白噪声

)(k v 的污染。根据过程的输入和输出数据{})(),(k z k u ，利用相关分析法计算出系

统的脉冲响应值)(ˆ

k g ，并与系统的脉冲响应理论值)(k g 0比较，得到系统脉冲响

应估计误差值)(~k g

，当∞→k 时，应该有0→)(~k g 。

1、模拟过程传递函数)(s G ，获得过程的输入和输出数据

})(),(k z k u （采样

时间取1秒）。

(1) 惯性环节

其中，T 为惯性环节的时间常数，K 为惯性环节的静态放大倍数。若采样时间记作0T ，则惯性环节的输出可写成：

k g =

)(ˆ

]

2T t k /∆