系统辨识广义最小二乘PPT课件_第1页
系统辨识广义最小二乘PPT课件_第2页
系统辨识广义最小二乘PPT课件_第3页
系统辨识广义最小二乘PPT课件_第4页
系统辨识广义最小二乘PPT课件_第5页
已阅读5页，还剩9页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

最小二乘与广义最小二乘法系统辨识,组别：第2组组员：史玄玄薛竹韵岳壮壮,.,目录,1.最小二乘法简介2.广义最小二乘法简介3.用最小二乘法与广义最小二乘法辨识系统4.总结,.,1.最小二乘法介绍,.,最小二乘估值算法,.,2.广义最小二乘法,.,广义最小二乘算法具体步骤,.,3.系统辨识实例：,.,.,最小二乘辨识结果：,.,（2）广义最小二乘对系统的参数估计,求出系统最小二乘估计uk=1.1470.201-0.787-1.159-1.0520.8661.1521.5730.6260.433-0.9580.810-0.0440.947-1.474-0.719-0.086-1.0991.4501.1510.4851.6330.0431.3261.706-0.3400.8901.1441.177-0.390;%输入mm=0.01;%白噪声的方差nn=normrnd(0,sqrt(mm),1,31);%产生个均值为0方差为mm的1x31的白噪声序列yk(1)=0;%得出测量输出ykyk(2)=0;fori=3:31yk(i)=-1.6*yk(i-1)-0.8*yk(i-2)+0.4*uk(i-1)+0.5*uk(i-2)+nn(i)+1.2*nn(i-1)+0.9*nn(i-2);endfori=1:29%求矩阵用A表示A(i,:)=-yk(i+1)-yk(i)uk(i+1)uk(i);endsita=inv(A*A)*A*yk(3:31);,.,迭代过程forj=1:60fori=3:31yk(i)=-1.6*yk(i-1)-0.8*yk(i-2)+0.4*uk(i-1)+0.5*uk(i-2)+nn(i)+1.6*nn(i-1)+0.8*nn(i-2);endfori=1:29A(i,:)=-yk(i+1)-yk(i)uk(i+1)uk(i);endsita=(inv(A*A)*A*yk(3:31);e(1)=yk(1);%计算残差e(2)=yk(2)+sita(1)*yk(1)-sita(3)*uk(1);fori=3:31e(i)=yk(i)+sita(1)*yk(i-1)+sita(2)*yk(i-2)-sita(3)*uk(i-1)-sita(4)*uk(i-2);end,.,广义最小二乘仿真结果,.,4.总结,通过编程计算，发现在噪声方差比较小的情况下，各种方法所获得的估值比较理想.但随着噪声方差的增大，用
人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论