廣義最小方差控制_第1頁
廣義最小方差控制_第2頁
廣義最小方差控制_第3頁
廣義最小方差控制_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、精選優(yōu)質(zhì)文檔-傾情為你奉上%廣義最小方差控制(顯示控制)考慮如下系統(tǒng)：式中x(k)為方差為0.1的白噪聲。取，期望輸出yr(k)為幅值為10的方波信號。clear all;close all;a=1 -1.7 0.7;b=1 2;c=1 0.2;d=4;na=length(a)-1;nb=length(b)-1;nc=length(c)-1;nf=nb+d-1;ng=na-1;P=1;R=1;Q=2; %加權(quán)多項式np=length(P)-1;nr=length(R)-1;nq=length(Q)-1;L=400;uk=zeros(d+nb,1);yk=zeros(na,1);yrk=zero

2、s(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);e,f,g=singlediophantine(a,b,c,d);CQ=conv(c,Q);FP=conv(f,P);CR=conv(c,R);GP=conv(g,P);for k=1:L time(k)=k; y(k)=-a(2:na+1)*yk+b*uk(d:d+nb)+c*xi(k);xik; u1=-Q(1)*CQ(2:nc+nq+1)*uk(1:nc+nq)/b(1)-FP(

3、2:np+nf+1)*uk(1:np+nf); u2=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)*CQ(1)/b(1)+FP(1); %更新數(shù)據(jù) for i=d+nb:-1:2 uk(i)=uk(i-1); end uk(1)=u(k); for i=na:-1:2 yk(i)=yk(i-1); end yk(1)=y(k); for i=nc:-1:2 yrk(i)=yrk(i-1); xik(i)=xik(i-1); end if nc>0 yrk(1)=yr(k); xik(1)=xi(k); endendsubplot(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)');subplot(2,1,2);plot(tim

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論