广义预测控制matlab仿真程序

广义预测控制matlab仿真程序！希望对大家学预测控制有所帮助！

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）;
% note that nAD = n~a + 1
nAD = size（AD）;
nAD = nAD（2）;
% compute horizons
 N1 = d +1;
 N2 = d + N;
% Compute F（z^-1）
% inilialization vector f
f（1:）= [1 zeros（1nAD-2）];
% 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:）;

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

     文件        740  2009-09-25 11:05  gpc\diophantine.m

     文件        450  2009-09-24 09:58  gpc\getG.m

     文件        174  2009-09-25 10:47  gpc\getGij.m

     文件       1959  2009-09-26 07:59  gpc\gpc.m

     文件         85  2009-09-25 11:10  gpc\rankn.m

     文件        836  2009-09-24 09:57  gpc\staten.m

     文件        577  2009-09-24 09:57  gpc\tf2AB.m

     目录          0  2009-09-26 07:58  gpc

----------- ---------  ---------- -----  ----

                 4821                    8