预测控制多种算法matlab实现

预测控制课程详细建议，包含DMC，MAC，Gpc和MPC等各类算法介绍

%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个文件信息