用matlab实现最小二乘递推算法辨识系统参数_第1页
用matlab实现最小二乘递推算法辨识系统参数_第2页
用matlab实现最小二乘递推算法辨识系统参数_第3页
全文预览已结束

付费下载

下载本文档

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

文档简介

1、自动化专业综合设计报告自动化专业综合设计报告设计题目：最小二乘递推算法辨识系统参数所在实验室: 自动化系统仿真实验室指导教师:学生姓名 班级计082-2班学号撰写时间： 2012-3-1 成绩评定： 一. 设计目的1、学会用Matlab实现最小二乘法辨识系统参数。2、进一步熟悉Matlab的界面及基本操作；3、了解并掌握Matlab中一些函数的作用与使用；二. 设计要求最小二乘递推算法辨识系统参数，利用matlab编程实现，设初始参数为零z(k)-1.5*z(k-1)+0.7*z(k-2)=1*u(k-1)+0.5*u(k-2)+v(k);选择如下形式的辨识模型：z(k)+a1*z(k-1)+

2、a2*z(k-2)=b1*u(k-1)+b2*u(k-2)+v(k);三. 实验程序m= 3;N=100;uk=rand(1,N);for i=1:Nuk(i)=uk(i)*(-1)A(i-1);endyk=zeros(1,N);for k=3:Nyk(k)=1.5*yk(k-1)-0.7*yk(k-2)+uk(k-1)+0.5*uk(k-2);end%j=100;kn=0;%y=yk(m:j)'%psi=yk(m-1:j-1);yk(m-2:j-2);uk(m-1:j-1);uk(m-2:j-2)'%pn=inv(psi'*psi);%theta=(inv(psi&#

3、39;*psi)*psi'*y);theta=0;0;0;0;pn=10A6*eye(4);for t=3:Nps=(yk(t-1);yk(t-2);uk(t-1);uk(t-2);pn=pn-pn*ps*ps'*pn*(inv(1+ps'*pn*ps);theta=theta+pn*ps*(yk(t)-ps'*theta);thet=theta'a1=thet(1);a2=thet (2);b1=thet(3);b2=thet(4);a1t(t)=a1;a2t(t)=a2;b1t(t)=b1;b2t(t)=b2;endt=1:N;plot(t,a1t(

4、t),t,a2t(t),t,b1t(t),t,b2t(t);text(20,1.47,'a1');text(20,-0.67,'a2');text(20,0.97,'b1');text(20,0.47,'b2');四. 设计实验结果及分析实验结果图:仿真结果表明，大约递推到第十步时，参数辨识的结果基本到稳态状态，即 a1 = 1.5999, b1=1, c1=0.5, d1=-0.7。五. 设计感受这周的课程设计告一段落了，时间短暂，意义重大。通过这次次练习的机会， 重新把matlab课本看了一遍，另外学习了系统辨识的有关内容，收获颇丰。对 matlab的使用更加纯熟，也锻炼了自己在课本中搜索信息和知识的能力。在设 计过程中虽然遇到了一些问题，但经过一次乂一次的思考，一遍乂一遍的检查 终丁找出了原因所在，也暴

人人文库> 全部分类> 行业资料 > 信息产业

温馨提示

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

最新文档

评论

0/150

提交评论