广义最小方差自校正控制直接算法_第1页
广义最小方差自校正控制直接算法_第2页
广义最小方差自校正控制直接算法_第3页
广义最小方差自校正控制直接算法_第4页
广义最小方差自校正控制直接算法_第5页
已阅读5页,还剩2页未读 继续免费阅读

下载本文档

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

文档简介

1、%广义最小方差自校正控制(直接算法)考虑如下系统:y(k)-1.7y(k-1)0.7y(k-2)=u(k-4)0.5u(k-5)(k)0.2(k-1)式中Hk)为方差为0.1的白噪声。取P(z,)=1,R(z,)=1,Q(z,)=2,期望输出yr(k)为幅值为10的方波信号。clearall;closeall;a=1-1.70.7;b=12;c=10.2;d=4;na=length(a)-1;nb=length(b)-1;nc=length(c)-1;nf=nb+d-1;ng=na-1;Pw=1;R=1;Q=2;%加权多项式np=length(Pw)-1;nr=length(R)-1;nq=l

2、ength(Q)-1;L=400;uk=zeros(d+nf,1);yk=zeros(d+ng,1);yek=zeros(nc,1);yrk=zeros(nc,1);xik=zeros(nc,1);%xiek=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);%RELS初值设定thetaek=zeros(na+nb+d+nc,d);P=10八6*eye(na+nb+d+nc);fork=1:Ltime(k)=k;y(k)=-a(2:na+1)*yk(1:na

3、)+b*uk(d:d+nb)+c*xi(k);xik;%递推最小二乘法phie=yk(d:d+ng);uk(d:d+nf);-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);%最优预测输出估计值%xie=y(k)-phie'*thetae(:,k);%白噪声估计值%提取辨识参数ge=thetae(1:ng+1,k)'f

4、e=thetae(ng+2:ng+nf+2,k)'ce=1thetae(ng+nf+3:ng+nf+2+nc,k)'ifabs(ce(2)>0.9ce(2)=sign(ce(2)*0.9;endiffe(1)<0.1fe(1)=0.1;end%e,f,g=singlediophantine(ae,be,ce,d);CQ=conv(ce,Q);FP=conv(fe,Pw);CR=conv(ce,R);GP=conv(ge,Pw);u1=-Q(1)*CQ(2:nc+nq+1)*uk(1:nc+nq)/fe(1)-FP(2:np+nf+1)*uk(1:np+nf);u2=

5、CR*yr(k+d:-1:k+d-min(d,nr+nc);yrk(1:nr+nc-d);u(k)=(u1+u2-GP*y(k);yk(1:np+ng)/(Q(1)*Q(1)/fe(1)+fe(1);%更新数据fori=d:-1:2thetaek(:,i)=thetaek(:,i-1);endthetaek(:,1)=thetae(:,k);fori=d+nf:-1:2uk(i)=uk(i-1);enduk(1)=u(k);fori=d+ng:-1:2yk(i)=yk(i-1);endyk(1)=y(k);fori=nc:-1:2yek(i)=yek(i-1);yrk(i)=yrk(i-1);

6、endifnc>0yek(1)=ye;yrk(1)=yr(k);xik(1)=xi(k);endendfigure(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(0L-2020);subplot(2,1,2);plot(time,u);xlabel('k');ylabel('u(k)');axis(0L

7、-1010);figure(2);subplot(211);plot(1:L,thetae(1:ng+1,:),1:L,thetae(ng+nf+3:ng+2+nf+nc,:);xlabel('k');ylabel('参数估计g,c');legend('g_0','g_1','c_1');axis(0L-34);subplot(212);plot(1:L,thetae(ng+2:ng+2+nf,:);xlabel('k');ylabel('参数估计f);klegend('f_0','f_1','f_2','f_3','f_4');axis(0L08);106S0n-5-10050100150200250300

温馨提示

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

评论

0/150

提交评论