GPC广义预测控制matlab例程

广义预测控制，根据当前输入和预测未来输出，控制系统的输入，使输出达到预期的性能指标。

function [EFGH]=Diophantine（NaNLPMLAB）
%Solving diophantine equation
%Na----Order of demoninator polynoMial;
%NL----Order of numerator polynoMial;
%P-----Predictive horizon;
%M-----Control horizon;
%L-----Step of tiMe delay;
%E-----P-dimension vector composed of noise coeficients;
%F-----P×（Na+1） matrix relevant to past controlled variable;
%G-----P×M matrix relevant to future manipulated variable;
%H-----P×NL matrix relevant to future manipulated variable;

a=[A0];
aa=[0A];

for i=1:L
   b（i）=0;
end

b=[bB];
aa=a-aa;
E（1）=1;
F（11:Na+1）=-aa（2:Na+2）;
g（1）=b（1）;
H（11:NL）=b（2:NL+1）;

for i=1:P-1
   E（i+1）=F（i1）;
   for j=1:Na;
      F（i+1j）=F（ij+1）-aa（j+1）*F（i1）;
   end
   F（i+1Na+1）=-aa（Na+2）*F（i1）;
end

for i=1:P-1
   g（i+1）=E（i+1）*b（1）+H（i1）;
   for j=1:NL-1
      H（i+1j）=E（i+1）*b（j+1）+H（ij+1）;
   end
   H（i+1NL）=E（i+1）*b（NL+1）;
end

for j=1:M
   for i=j:P
      G（ij）=g（i-j+1）;
   end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2269  2009-04-17 18:21  预测控制\GPCMatlab.rar

     文件       1270  2009-06-08 23:07  预测控制\GPCPrg1\CONPARA.cpp

     文件       1276  2008-09-30 21:19  预测控制\GPCPrg1\CONPARA.h

     文件       3917  2009-06-10 11:19  预测控制\GPCPrg1\GPCPrg.clw

     文件       4209  2008-09-30 19:24  预测控制\GPCPrg1\GPCPrg.cpp

     文件       4930  2008-09-30 21:05  预测控制\GPCPrg1\GPCPrg.dsp

     文件        535  2008-09-30 19:24  预测控制\GPCPrg1\GPCPrg.dsw

     文件       1356  2008-09-30 19:24  预测控制\GPCPrg1\GPCPrg.h

     文件      53760  2009-06-09 10:21  预测控制\GPCPrg1\GPCPrg.opt

     文件       1250  2009-06-10 06:25  预测控制\GPCPrg1\GPCPrg.plg

     文件       1742  2008-09-30 19:24  预测控制\GPCPrg1\GPCPrgDoc.cpp

     文件       1475  2008-09-30 19:24  预测控制\GPCPrg1\GPCPrgDoc.h

     文件      10297  2009-06-10 06:25  预测控制\GPCPrg1\GPCPrgView.cpp

     文件       2280  2008-10-26 16:50  预测控制\GPCPrg1\GPCPrgView.h

     文件       2371  2008-10-03 15:52  预测控制\GPCPrg1\MainFrm.cpp

     文件       1440  2008-10-03 15:51  预测控制\GPCPrg1\MainFrm.h

     文件       1197  2009-06-08 23:07  预测控制\GPCPrg1\MODELPARA.cpp

     文件       1262  2008-09-30 21:20  预测控制\GPCPrg1\MODELPARA.h

     文件       4311  2008-09-30 19:24  预测控制\GPCPrg1\ReadMe.txt

     文件       1099  2008-09-30 20:48  预测控制\GPCPrg1\Resource.h

     文件       1023  2008-09-30 20:00  预测控制\GPCPrg1\Sheet.cpp

     文件       1370  2008-09-30 21:03  预测控制\GPCPrg1\Sheet.h

     文件        208  2008-09-30 19:24  预测控制\GPCPrg1\StdAfx.cpp

     文件       1054  2008-09-30 19:24  预测控制\GPCPrg1\StdAfx.h

     文件      24064  2008-10-07 10:43  预测控制\GPCPrg1\参数说明.doc

     文件     189440  2009-06-10 19:33  预测控制\GPCPrg1\GPCPrg.ncb

     文件       1078  2008-09-30 19:24  预测控制\GPCPrg1\res\GPCPrg.ico

     文件        398  2008-09-30 19:24  预测控制\GPCPrg1\res\GPCPrg.rc2

     文件       1078  2008-09-30 19:24  预测控制\GPCPrg1\res\GPCPrgDoc.ico

     文件       1078  2009-10-09 16:08  预测控制\GPCPrg1\res\Toolbar.bmp

............此处省略35个文件信息