广义预测控制MATLAB算法

广义预测控制的MATLAB实现。文件夹共8个子m文件。其中gpc文件是主程序，其他为调用函数。说明很仔细，希望可以帮到各位学习预测控制的同学。

function[EF] = diophantine（ANd）
% Computes polynomials E（z^-1） e F（z^-1） 
% delta = 1-z^（-1）
delta = [1 -1];
% A = 1 + a1 z^（-1） + ... + a2 z（-na）
AD = conv（Adelta）;%所求AD既是A~的参数 AD=（1，a1~a2~...ana~）
% note that nAD = n~a + 1
nAD = size（AD）;
nAD = nAD（2）;%nAD是A~参数的个数
% compute horizons
 N1 = d +1;
 N2 = d + N;%在这里N2相当于P?
% Compute F（z^-1）
% inilialization vector f  ！！！！！f是na次多项式
f（1:）= [1 zeros（1nAD-2）];%nAD-2的意思是nAD 已经是阶数+1，第一行第一列已有，是1，所以减2f的第一行是（1 0 0 0）
% i = 0 ... n?-1
for j = 1: N2;
% Note that for i = 1 corresponds to f（j0）
for i = 1:nAD-2
   f（j+1i） = f（ji+1）-f（j1）*AD（i+1）;
end
   f（j+1nAD-1） = -f（j1）*AD（nAD）;
end
F = f（1+N1:1+N2:）;%这个什么意思？
% Computes E（z^-1）
E = zeros（N2）;
e（1） = 1;               % for the special case  1/~A
E（11） = e（1）;
for i = 2: N2
    e（i） = f（i1）;
    E（i1:i）=e;
end
E = E（N1:N2:）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         450  2009-09-24 09:58  getG.m
     文件         107  2017-11-28 15:15  getGii.m
     文件        2084  2017-12-05 15:25  gpc.m
     文件          85  2009-09-25 11:10  rankn.m
     文件         836  2009-09-24 09:57  staten.m
     文件         577  2009-09-24 09:57  tf2AB.m
     文件         945  2017-11-27 22:41  diophantine.m