%gpc1（AplantBplantqpaifa）.m
clear?all;
close?all;
Amodel=[1?-1.4?0.49?-0.343];
%A（q^-1）=1+a1*q^-1+a2*q^-2+......+ana*q^-na
Bmodel=[-0.2?1];
%B（q^-1）=bo+b1*q^-1+b2*q^-2+......+bnb*q^-nb
%Amodel=[1-0.9355];
%Bmodel=[0.06449];
na=length（Amodel）-1;
nb=length（Bmodel）-1;
k=4;
p=14;
SPt=10;
aifa=0.99;
Ts=1;
T_final=500;
lambda=0;
M=6;
Aplant=Amodel;
Bplant=Bmodel;
F（11:na）=Amodel（11:na）-Amodel（12:na+1）;
F（1na+1）=Amodel（1na+1）;
%F1（q^-1）=（1-a1）+（a1-a2）*q^-1+...+ana*q^-na
E（11）=1;
%F1=poly2sym（F1）
G_ef（1:）=Bmodel;
for?j=2:k+p-1
E（1j）=F（j-11）;
F（j1:na）=F（j-12:na+1）-E（1j）*（Amodel（12:na+1）-Amodel（11:na））;
F（jna+1）=F（j-11）*Amodel（1na+1）;
G_ef（j1:nb+j）=conv（EBmodel）;
end
for?i=1:p
????for?j=1:i
????????G（ij）=G_ef（ii-j+1）;
????end
end
for?i=1:p
G_poly（i1:nb+k-1）=G_ef（k+i-1i+1:nb+k+i-1）;
end
G1=G（1:p1:M）;
ypast（:1）=zeros（na+11）;
Upast（:1）=zeros（nb+k+11）;
outU=[];
outY=[];
for?I=0:Ts:T_final
????delt_ypast=ypast（1:na1）-ypast（2:na+11）;
deltU=Upast（1:nb+k1）-Upast（2:nb+k+11）;
yt=ypast（11）-Aplant（12:na+1）*delt_ypast（1:na1）+Bplant（1:）*deltU（k:nb+k1）;
outY=[outY;yt];
ypast=[yt;ypast];
ypast=ypast（1:na+11）;
for?i=1:p
ypre1（i1）=F（k+i-1:）*ypast（:1）+G_poly（i:）*deltU（1:nb+k-11）;
end
if?k==1
W（11）=yt;
else
W（11）=F（k-1:）*ypast（:1）+G_ef（k-11:nb+k-1）*deltU（1:nb+k-11）;
end
for?i=2:p+1
W（i1）=aifa*W（i-11）+（1-aifa）*SPt;
end
deltU1=inv（G1‘*G1+lambda*eye（M））*G1‘*（W（2:p+11）-ypre1（:1））;
ut=Upast（11）+deltU1（11）;
Upast=[ut;Upast];
Upast=Upast（1:nb+k+11）;
outU=[outU;ut];
end
t=[0:Ts:T_final];
plot（toutY）hold?onplot（toutU‘r‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1645??2009-07-22?22:25??0579950311111\廣義預測控制（GPC）仿真.m

?????文件???????2754??2010-01-24?14:08??539347816GPC\GPC\GPC.m

?????文件???????2792??2010-01-24?14:07??539347816GPC\GPC\GPC_zijiaozheng.m

?????文件?????142848??2019-03-18?12:26??539347816GPC\GPC\GPC說明文檔.doc

?????文件????????315??2019-03-18?08:19??63750275test\mohu.asv

?????文件????????343??2019-03-18?08:24??63750275test\mohu.m

?????文件???????2045??2019-03-18?15:57??63750275test\test.m

?????文件???????1037??2011-02-19?11:52??63535277GPC-matlab\my_diophantine.m

?????文件???????2095??2011-03-02?10:55??63535277GPC-matlab\my_main.m

?????文件???????1100??2010-03-29?08:17??47651480GPC\GPC.M

?????文件????????752??2019-01-19?21:38??GPC\diophantine.m

?????文件????????453??2019-01-20?21:47??GPC\getG.m

?????文件????????178??2019-01-19?21:39??GPC\getGij.m

?????文件???????2034??2019-03-17?18:40??GPC\gpc.m

?????文件???????2792??2010-01-24?14:07??GPC\GPC_zijiaozheng.m

?????文件?????145408??2010-01-24?14:09??GPC\GPC說明文檔.doc

?????文件?????????97??2019-01-23?22:57??GPC\rankn.m

?????文件????????593??2019-01-19?15:55??GPC\tf2AB.m

?????文件?????136704??2019-03-14?20:47??GPC\多變量廣義預測控制.doc

?????文件????????259??2009-08-24?10:33??chap4\chap4_01_sindiophantine.m

?????文件????????225??2009-05-03?13:20??chap4\chap4_02_multidiophantine.m

?????文件???????1252??2009-07-03?21:35??chap4\chap4_03_MVC.m

?????文件???????2317??2009-08-24?11:48??chap4\chap4_04_MVSTC_indirect.m

?????文件???????2425??2009-08-24?15:04??chap4\chap4_05_MVSTC_direct.m

?????文件???????1522??2009-07-03?21:31??chap4\chap4_06_GMVC.m

?????文件???????2384??2009-08-23?16:18??chap4\chap4_07_GMVSTC_indirect.m

?????文件???????2694??2009-08-24?15:03??chap4\chap4_08_GMVSTC_direct.m

?????文件???????1728??2009-07-03?22:16??chap4\chap4_09_GPC_NoIden.m

?????文件???????2486??2009-07-03?21:25??chap4\chap4_10_GPC_Iden.m

?????文件???????2240??2009-08-23?20:00??chap4\chap4_11_GPC_NoIden_C.m

............此處省略18個文件信息