版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、用最小方差自校正控制 算法对以下系统进行闭环控制:y(k) 1.7y(k 1) 0.7y(k2) u(k 4)0.5u(k5)(k)0.2 (k 1)式中(k)为方差为0.1的白噪声,取期望输出yr(k)为幅值为10的方波 信号。解:上式可以化为:(1 1.7z 10.7z2)y(k) z4(1 0.5z 1)u(k) (1 0.2z 1) (k)则有11A(z )11.7z11B(z )10.5zC(z 1)10.2z1d 4Diopha ntine 方程为:又有6取初值 P(0)=10 , (0)0递推公式为:?(k)办K(k)P(k)0.7zC(z1)H(z1)A(z1)F(z1) zd
2、G(z1)B(z1)F(z1)1) K(k)y(k)?T(k d)?(k 1)P(k 1)?(k d)1?r(k d)P(k 1)?(k d)I K(k)?T(k d)P(k 1)程序清单如下:clear all; close all;na=length(a)-1; nb=length(b)-1; nc=length(c)-1;%多项式 A 、B、C 的阶次nh=nb+d-1; ng=na-1;%nh、ng 为多项式 H、G 的阶次L=400;uk=zeros(d+nh,1); %输入初值:yk=zeros(d+ng,1); %输出初值yek=zeros(nc,1); %最优输出预测估计初值
3、yrk=zeros(nc,1); %期望输出初值 xik=zeros(nc,1); %白噪声初值 yr=10*ones(L/4,1);-ones(L/4,1);ones(L/4,1);-ones(L/4+d,1); % 期望输出 xi=sqrt(0.1)*randn(L,1); % 白噪声序列%递推估计初值 thetaek=zeros(na+nb+d+nc,d);P=10A6*eye( na+n b+d+nc); %P(k )的初始值for k=1:Ltime(k)=k;y(k)=-a(2:na+1)*yk(1:na)+b*uk(d:d+nb)+c*xi(k);xik; % 采集输出数据%递推
4、增广最小二乘法公式估计参数 phie=yk(d:d+ng);uk(d:d+nh);-yek(1:nc);K=P*phie/(1+phie*P*phie); thetae(:,k)=thetaek(:,1)+K*(y(k)-phie*thetaek(:,1);P=(eye(na+nb+d+nc)-K*phie)*P;ye=phie*thetaek(:,d); %预测输出的估计值%提取辨识参数ge=thetae(1:ng+1,k); he=thetae(ng+2:ng+nh+2,k);ce=1 thetae(ng+nh+3:ng+nh+2+nc,k);if abs(ce(2)0.9 ce(2)=s
5、ign(ce(2)*0.9;endif he(1)0yek(1)=ye;yrk(1)=yr(k);xik(1)=xi(k);endend figure(1); subplot(2,1,1); plot(time,yr(1:L),r:,time,y); xlabel(k); ylabel(y_r(k) 、 y(k); legend(y_r(k),y(k); axis(0 L -20 20); subplot(2,1,2);plot(time,u); xlabel(k); ylabel(u(k); axis(0 L -40 40);figure(2) subplot(211) plot(1:L,thetae(1:ng+1,:),1:L,thetae(ng+nh+3:ng+2+nh+nc,:);xlabel(k); ylabel(参数估计 g、c); legend(g_0,g_1,c_1); axis(0 L -3 4);subplot(212) plot(1:L,thetae(ng+2:ng+2+nh,:);xlabel(k); ylabel(参数估计 h); legend(h_0,h_1,h_2,h_3
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 华师大版初中科学1.2 水的三态变化(30课件)
- 20XX年1月华懋达集团年会庆典概念方案
- 2024年烟台货运资格证模拟考试题
- 算法设计与分析 课件 5.9-动态规划应用-最优二叉搜索树
- 2024年宣城客运资格证考试答题
- 2024年贵州客运从业资格证的考试题目是什么题
- 吉首大学《结构试验》2021-2022学年第一学期期末试卷
- 吉首大学《当代中国电影》2021-2022学年期末试卷
- 《机床夹具设计》试题4
- 吉林艺术学院《音乐文论写作Ⅱ》2021-2022学年第一学期期末试卷
- 6.2共筑生命家园 课件-2024-2025学年统编版道德与法治九年级上册
- 2024年秋新沪教牛津版英语三年级上册课件 Unit 5 第2课时
- 中图版七年级下册信息技术 2.1规划影片任务 教学设计
- 2024中科院心理咨询师考试复习题库(官方版)-上单选题汇
- 小学未成年人思想道德建设工作实施方案
- 化工公司安全知识竞赛题库(共1000题)
- GB/T 44421-2024矫形器配置服务规范
- 福建省福州市(2024年-2025年小学二年级语文)统编版期中考试试卷(含答案)
- 2024-2024部编版九年级语文上册期末考试测试卷(附答案)
- 争做“四有好老师”-当好“四个引路人”
- 2024-2025学年八年级生物上册第一学期 期末综合模拟测试卷( 人教版)
评论
0/150
提交评论