系统辨识实验手册.doc

《系统辨识》实验手册哈尔滨工业大学 控制与仿真中心2012年8月 / 目 录实验1 白噪声和M序列的产生---------------------------------------------------------- 2实验2 脉冲响应法的实现---------------------------------------------------------------- 5实验3 最小二乘法的实现--------------------------------------------------------------- 9实验4 递推最小二乘法的实现---------------------------------------------------------- 12附录 实验报告模板---------------------------------------------------------------------- 16实验1 白噪声和M序列的产生一、实验目的1、熟悉并掌握产生均匀分布随机序列方法以及进而产生高斯白噪声方法2、熟悉并掌握M序列生成原理及仿真生成方法二、实验原理1、混合同余法混合同余法是加同余法和乘同余法的混合形式，其迭代式如下：式中a为乘子，为种子，b为常数，M为模。

混合同余法是一种递归算法，即先提供一个种子，逐次递归即得到一个不超过模M的整数数列2、正态分布随机数产生方法由独立同分布中心极限定理有：设随机变量相互独立，服从同一分布，且具有数学期望和方差：则随机变量之和的标准化变量:近似服从分布如果服从[0, 1]均匀分布，则上式中，即近似服从分布3、M序列生成原理用移位寄存器产生M序列的简化框图如下图所示该图表示一个由4个双稳态触发器顺序连接而成的4级移位寄存器，它带有一个反馈通道当移位脉冲来到时，每级触发器的状态移到下一级触发器中，而反馈通道按模2加法规则反馈到第一级的输入端三、实验内容1、生成均匀分布随机序列（1）利用混合同余法生成[0, 1]区间上符合均匀分布的随机序列，并计算该序列的均值和方差，与理论值进行对比分析要求序列长度为1200，推荐参数为a=65539，M=2147483647，0

3、生成M序列M序列的循环周期取为，时钟节拍，幅度，逻辑“0”为a，逻辑“1”为-a，特征多项式生成M序列的结构图如下所示C1C2C3C4C5C6CPM(6)M(5)+M(4)M(3)M(2)M(1)M(0)要求编写Matlab程序生成该M序列，绘制该信号曲线，并分析验证M序列的性质四、实验步骤1．分别画出三部分实验内容的程序框图（流程图）；2．编制MATLAB的M文件；3．运行编制的M文件；4．查看程序运行结果并进行分析；5．填写实验报告五、实验报告格式参见附录一实验2 脉冲响应法一、实验目的通过仿真实验掌握利用相关分析法辨识脉冲响应的原理和方法二、实验原理一个单入单出线性定常系统的动态特性可用它的脉冲响应函数g(σ)来描述这样，只要记录x(t)、y(t)的值，并计算它们的互相关函数，即可求得脉冲响应函数g(τ)而在系统有正常输入的情形下，辨识脉冲响应的原理图如下图所示三、实验内容下图为本实验的原理框图系统的传递函数为,其中；分别为系统的输入和输出变量；为测量白噪声，服从正态分布，均值为零，方差为，记作；为系统的脉冲响应理论值，为系统脉冲响应估计值，为系统脉冲响应估计误差相关分析法v(k)u(k)z(k)系统的输入采用M序列（采用实验1中的M序列即可），输出受到白噪声的污染。

根据过程的输入和输出数据，利用相关分析法计算出系统的脉冲响应值，并与系统的脉冲响应理论值比较，得到系统脉冲响应估计误差值，当时，应该有1、模拟过程传递函数，获得过程的输入和输出数据（采样时间取1秒）u(k)y(k)(1) 惯性环节其中，T为惯性环节的时间常数，K为惯性环节的静态放大倍数若采样时间记作，则惯性环节的输出可写成：(2) 传递函数仿真（串联） u(k)x(k)y(k)令，则的表达框图为：2、互相关函数的计算其中，r为周期数，表示计算互相关函数所用的数据是从第二个周期开始的，目的是等过程仿真数据进入平稳状态可分别令r =1、3，对比仿真结果）3、c的补偿补偿量c应取，不能取因为是周期函数，则有，故不能取4、计算脉冲响应估计值● 脉冲响应估计值 ● 脉冲响应估计误差 四、实验步骤(1) 掌握相关分析辨识方法的基本原理；(2) 设计实验方案，画出程序框图；(3) 编制实验程序；(4) 调试并运行程序，记录数据；(5) 分析实验结果，完成实验报告五、实验报告格式参见附录一实验3 最小二乘法的实现一、实验目的理解并掌握系统辨识中的最小二乘法原理二、实验原理给定系统 （1）其中，为待辨识的未知参数，是不相关随机序列。

为系统的输出，为系统的输入分别测出个输出、输入值，则可写出个方程，具体写成矩阵形式，有 （2）设，则式（2）可写为 （3）式中：y为N维输出向量；为N维噪声向量；为维参数向量；为测量矩阵为了尽量减小噪声对估值的影响，应取，即方程数目大于未知数数目的最小二乘估计为 （4）三、实验内容对象的数学模型如下：其中，是服从正态分布的白噪声N输入信号采用4阶M序列，幅度为1选择如下形式的辨识模型：设输入信号的取值是从k =1到k =16的M序列，则待辨识参数为=其中，被辨识参数、观测矩阵z L、H L的表达式为 , ， 要求编制仿真程序，获取系统输入输出数据，并运用最小二乘法对这一系统的参数进行辨识，并将辨识结果与实际参数进行对比四、实验步骤1．写出系统结构、实际参数、噪声源及输入信号等内容；2．画出程序框图（流程图）；3．编制MATLAB的M文件；4．运行编制的M文件；5．查看程序运行结果并进行分析；6．填写实验报告五、实验报告格式参见附录一实验4 递推最小二乘法的实现一、实验目的熟悉并掌握递推最小二乘法的算法原理。

二、实验原理 给定系统 （1）其中，为待辨识的未知参数，是不相关随机序列为系统的输出，为系统的输入分别测出个输出、输入值，则可写出个方程，具体写成矩阵形式，有 （2）设，则式（2）可写为 （3）式中：y为N维输出向量；为N维噪声向量；为维参数向量；为测量矩阵为了尽量减小噪声对估值的影响，应取，即方程数目大于未知数数目的最小二乘估计为 （4）为了实现实时控制，必须采用递推算法，这种辨识方法主要用于辨识设已获得的观测数据长度为N，将式（3）中的、和分别用来代替，即 （5）用表示的最小二乘估计，则 （6）令，则 （7）如果再获得一组新的观测值和，则又增加一个方程 （8）式中将式（5）和式（8）合并，并写成分块矩阵形式，可得 （9） 于是，类似地可得到新的参数估值 （10）式中 （11） 应用矩阵求逆引理，从求得与的递推关系式出发，经过一系列的推导，最终可求得递推最小二乘法辨识公式： （12） （13） （14）为了进行递推计算，需要给出和的初值和。

推荐取值方法为：假定，c是充分大0的常数，I为单位矩阵，则经过若干次递推之后能得到较好的参数估计三、实验内容给定系统 （15）即假设实际系统的参数为，，，，，但是不已知，即不可测取的零均值随机序列输入信号取为 （16） 要求编制MATLAB程序，运用递推最小二乘法对这一系统的参数进行辨识，并将辨识结果与实际参数进行对比四、实验步骤1．写出系统结构、实际参数、噪声源及输入信号等内容；2．画出程序框图；3．编制MATLAB的M文件；4．运行上一步编制的M文件；5．将辨识结果与实际参数对比分析；6．修改程序、重新运行，直至结果符合精度要求；7．填写实验报告五、实验报告 格式参见附录一附录实验1 ××××××实验报告哈尔滨工业大学航天学院控制科学与工程系专业： 班级： 姓名： 日期： 年 月 日1．实验题目： 2．实验目的3．实验主要原理4．实验对象或参数5．程序框图6．程序代码7．实验结果及分析8．结论 友情提示：方案范本是经验性极强的领域，本范文无法思考和涵盖全面，供参考!最。