没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
《系统建模与仿真》
实验手册(学生用)
2006 年 6 月
目 录
实验 1 最小二乘法的实现------------------------------------------------------------------------------- 2
实验 2 龙格-库塔法的实现------------------------------------------------------------------------------ 5
实验 3 随动控制系统仿真------------------------------------------------------------------------------- 6
1
实验 1 最小二乘法的实现
一、实验目的
熟悉并掌握递推最小二乘法的算法原理
二、实验原理——递推最小二乘法
给定系统
(1)
其中 , 为待辨识的未知参数, 是不相关随机序列。 为
系 统 的 输 出 , 为 系 统 的 输 入 。 分 别 测 出 个 输 出 、 输 入 值
,则可写出 个方程,具体写成矩阵
形式,有
(2)
设
,
则式(2)可写为
(3)
式中:y 为 N 维输出向量; 为 N 维噪声向量; 为 维参数向量; 为
测量矩阵。为了尽量减小噪声 对 估值的影响,应取 ,即方程数目大于未知
2
数数目。
的最小二乘估计为
(4)
为了实现实时控制,必须采用递推算法,这种辨识方法主要用于在线辨识。设已获得
的观测数据长度为 N,将式(3)中的 、 和 分别用 来代替,即
(5)
用 表示 的最小二乘估计,则
(6)
令 ,则
(7)
如果再获得一组新的观测值 和 ,则又增加一个方程
(8)
式中
将式(5)和式(8)合并,并写成分块矩阵形式,可得
(9)
于是,类似地可得到新的参数估值
(10)
式中
3
(11)
应用矩阵求逆引理,从求得 与 的递推关系式出发,经过一系列的推导,最终可求
得递推最小二乘法辨识公式:
(12)
(13)
(14)
为了进行递推计算,需要给出 和 的初值 和 。
推荐取值方法为:假定 ,c 是充分大的常数,I 为 单
位矩阵,则经过若干次递推之后能得到较好的参数估计。
三、实验内容
给定系统
(15)
即 。假设实际系统的参数为 , , , , ,但
是不已知,即不可测。取 的零均值白噪声。输入信号取为
(16)
要求编制 MATLAB 程序,运用递推最小二乘法对这一系统的参数进行在线辨识,并将辨
识结果与实际参数进行对比。
四、实验步骤
1.写出系统结构、实际参数、噪声源及输入信号等内容;
2.画出程序框图;
3.编制 MATLAB 的 M 文件;
4.运行上一步编制的 M 文件;
5.将辨识结果与实际参数对比分析;
6.修改程序、重新运行,直至结果符合精度要求;
7.填写实验报告。
五、实验报告
见附录一。
4
剩余20页未读,继续阅读
资源评论
- zsn2212016-11-10不错,通过实践验证
- ZP19002012-10-11系统建模与仿真实验
- jinyanan19892014-02-17不错,已实践通过
lanbeilei
- 粉丝: 3
- 资源: 10
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功