![系统辨识实验报告_第1页](http://file3.renrendoc.com/fileroot_temp3/2022-1/17/9650b444-b25c-4cea-90fb-718b3c93f84b/9650b444-b25c-4cea-90fb-718b3c93f84b1.gif)
![系统辨识实验报告_第2页](http://file3.renrendoc.com/fileroot_temp3/2022-1/17/9650b444-b25c-4cea-90fb-718b3c93f84b/9650b444-b25c-4cea-90fb-718b3c93f84b2.gif)
![系统辨识实验报告_第3页](http://file3.renrendoc.com/fileroot_temp3/2022-1/17/9650b444-b25c-4cea-90fb-718b3c93f84b/9650b444-b25c-4cea-90fb-718b3c93f84b3.gif)
![系统辨识实验报告_第4页](http://file3.renrendoc.com/fileroot_temp3/2022-1/17/9650b444-b25c-4cea-90fb-718b3c93f84b/9650b444-b25c-4cea-90fb-718b3c93f84b4.gif)
![系统辨识实验报告_第5页](http://file3.renrendoc.com/fileroot_temp3/2022-1/17/9650b444-b25c-4cea-90fb-718b3c93f84b/9650b444-b25c-4cea-90fb-718b3c93f84b5.gif)
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、系统辨识实验报告一、相关分析法(1)实验原理图1 实验原理图本实验的原理图如图1。过程传递函数中;输入变量,输出变量,噪声服从,为过程的脉冲响应理论值,为过程脉冲响应估计值,为过程脉冲响应估计误差。过程输入采用M序列,其输出数据加白噪声得到输出数据。利用相关分析法估计出过程的脉冲响应值,并与过程脉冲响应理论值比较,得到过程脉冲响应估计误差值。M序列阶次选择说明:首先粗略估计系统的过渡过程时间TS (通过简单阶跃响应)、截止频率fM (给系统施加不同周期的正弦信号或方波信号,观察输出)。本次为验证试验,已知系统模型,经计算,。根据式及式,则取值为1,此时,由于与N选择时要求完全覆盖,则选择六阶M
2、移位寄存器,即N=63。(2) 编程说明人机对话噪声标准差:sigma;生成数据周期数:r生成数据生成M序列;生成白噪声序列过程仿真得到理论输出数据计算脉冲响应估计值计算互相关函数,得到脉冲响应估计值计算脉冲响应估计误差计算脉冲响应理论值,得到脉冲响应估计误差图2 程序流程图(3) 分步说明 生成M序列:M序列的循环周期,时钟节拍,幅度,移位寄存器中第5、6位的内容按“模二相加”,反馈到第一位作为输入。其中初始数据设为1,0,1,0,0,0。程序如下: 生成白噪声序列:程序如下: 过程仿真得到输出数据:如图2所示的过程传递函数串联,可以写成形如,其中。图2 过程仿真方框图程序如下: 计算脉冲响
3、应估计值:互相关函数采用公式,互相关函数所用的数据是从第二个周期开始的,其中为周期数,取1-3之间。则脉冲响应估计值为:,。补偿量。程序如下: 计算脉冲响应估计值:脉冲响应的理论值由式可计算得到。这时可得到过程脉冲相应估计误差。脉冲响应估计误差为:程序如下:(4) 数据记录当噪声标准差sigma=0.1,生成数据周期r为2时:脉冲响应估计误差为0.0121。脉冲响应估计曲线为图3所示。当噪声标准差sigma=0.5,生成数据周期r为1时:脉冲响应估计误差为 0.0347。脉冲响应估计曲线为图4所示。当噪声标准差sigma=0.5,生成数据周期r为3时:脉冲响应估计误差为 0.0258。脉冲响应
4、估计曲线为图5所示。当噪声标准差sigma=1,生成数据周期r为3时:脉冲响应估计误差为 0.0279。脉冲响应估计曲线为图6所示。 图3 sigma=0.1,r=2时脉冲响应估计曲线 图4 sigma=0.5,r=1时脉冲响应估计曲线 图5 sigma=0.5,r=3时脉冲响应估计曲线 图6 sigma=1,r=3时脉冲响应估计曲线(5) 结果分析实验中可以看到脉冲响应估计的曲线与理论曲线的重合度还是比较高的,脉冲响应估计误差也比较小,实验证明相关分析法的估计效果还是不错的。同时,经过实验可以得出结论:固定数据周期r,给定不同的噪声标准差sigma可以发现,噪声的方差越大,也就是信噪比越大,
5、估计的效果越不好;固定噪声标准差sigma,选择不同的数据生成周期r可以发现,数据周期越大,估计的周期越多,估计的效果越好。二、最小二乘法1. 基本最小二乘(离线辨识)残差为: 最小二乘目标:残差平方和最小(一阶导为0,二阶导>0)。从上式看出,逆存在才有解,满足条件的u(k):(1) 伪随机;(2) 白噪声;(3)有色随机信号。程序如下: 结果如下:result1 = -0.8287; 0.1275; -0.0024; 1.9884; -1.27232. 递推最小二乘(在线辨识) RLS 为了启动RLS,需给初值:。计算框图见书P66。程序如下: 结果如下:result2 = -0.8
6、284; 0.1274; -0.0024; 1.9883; -1.2717图7 递推最小二乘法参数过渡过程数据饱和:(1)原因:,不再起修正作用,引起误差变大。 (2)为了克服数据饱和现象,可以用降低老的数据影响的方法:渐消记忆法(遗忘因子法) 限定记忆法(固定窗法)当为不相关序列,最小二乘有一致性与无偏性,但往往为相关序列,为克服最小二乘有偏估计的缺点,引入辅助变量法和广义最小二乘法,增广最小二乘法等。3. 渐消记忆法(遗忘因子法) 一般程序如下:结果如下:result3 = -0.6862; 0.1114; 0.0640; 2.0356; -0.9429图8 渐消记忆法参数过渡过程4. 限
7、定记忆法(固定窗法)程序如下:结果如下:(较前三种方法偏差较大)result4 = -1.0190; 0.1725; -0.4531; 2.1092; -1.37245. 辅助变量法(IV)(1).辅助变量Z (2) .计算步骤: 先根据实测数据最小二乘求粗略(为有偏估计) (3) .递推:RIV 初始条件: 缺点: P0的选择非常敏感,一个改进方法是,用递推最小二乘辨识算法作为启动方法,然后转换到辅助变量法。程序如下: 结果如下:result5 = -0.9369; 0.1207; -0.0254; 1.9781; -1.5331 图9 辅助变量法参数过渡过程6. 广义最小二乘法(GLS)(
8、1)广义最小二乘法的基本思想:由于在n+k个采样周期的时差范围内具有自相关性,从而使的最小二乘估计为有偏的,所以引入一个所谓成形滤波器(白化滤波器),把相关噪声转化成白噪声。如果知道有色噪声序列的相关性:令,有广义最小二乘法(GLS)是建立在最小二乘法(LS)的基础上的。基本最小二乘法只是广义最小二乘法在时的特例。(2)广义最小二乘法计算步骤:广义最小二乘法的关键问题是如何用比较简单的方法找到成形滤波器的系数。其计算是逐次逼近法。应用输入输出数据按最初模型求出的最小二乘估计。这个估计值是不精确的,它只是被估参数的一次近似。计算残差e(k),并拟合成形滤波器的模型:得到 其中 应用所得的成形滤波
9、器,对输入输出数据滤波:其中,m为噪声模型的阶,一般事先不知道,实际经验表明指定m为2或3可以得到比较满意的输出。按新的输入、输出模型求出参数的第二次估计值。结果如下:Result6 = -0.6538; 0.2926; 0.0454; 1.8776; -1.60427. 广义递推最小二乘法(GLS) 广义最小二乘法的递推计算过程可分成两个部分:(1) 按递推最小二乘法(RLS),随着N的增大,不断计算(逐步接近于无偏)和(逐步使噪声白化);(2) 在递推过程中,和是时变的,则过滤信号及残差是由时变系统产生,要不断计算。因而,递推广义最小二乘法由两组普通的递推最小二乘法组成,它们是通过滤波算法联系起来的: 结果如下:result7 = -1.0164; 0.1754; -0.0159; 2.0056; -1.6438噪声传递系数的估计结果:-0.0307; 0.0900图10 广义递推最小二乘参数过渡过程8. 增广矩阵法(ELS/RELS)(增广最小二乘法)增广矩阵法是把观测矩阵适当增大,使得有偏估计的程度得到一定改善。 这一方程结构适用于递推最小二乘法,但向量中是未知的。解决这个矛
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 零售业中的顾客安全保障措施
- DB3715T 69-2025研学旅游指导师服务规范
- 专业技术人才海外培训服务合同(版)
- 上海股权转让合同文本
- 二手房转让合同定金协议书范本
- 中外合资企业劳动合同样本
- 个人保证担保融资合同协议
- NBA赛事中国区电视转播合同
- 互利共赢投资合作合同
- 个人物流配送服务合同模板
- 关于防范遏制矿山领域重特大生产安全事故的硬措施课件
- 2025年中国成都餐饮业市场运营态势分析及投资前景预测报告
- 2024年xx县第三小学安全工作管理制度汇编
- 项目合作备忘录范文
- 妇产科医生个人年终述职报告课件
- 《费曼学习法》读后感
- 趣味成语课程设计
- 钢筋焊接工艺性试验方案
- 2024年四川省凉山州中考物理适应性试卷(附答案解析)
- 幼儿园师资培训课件【区域进阶】科学区各年龄段目标制定与投放材料:区域材料玩出新高度课件
- 人教版初中数学同步讲义八年级上册专题提升02 解分式方程与分式方程的实际应用(30题)(原卷版)
评论
0/150
提交评论