能运行很好的GPC程序。

GPC 课程设计随便用,没问题。所有中间过程都有详细注释。

function[EF] = diophantine（ANd）                                          %A 既输出Y系数 N预测步数 d纯迟延
% 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）;                                                        %conv是计算卷积的公式 。A1的第二项至最后一项是书中2.14公式A（z^-1）*△。
% note that nAD = n~a + 1
nAD = size（AD）;
nAD = nAD（2）;                                                              %获取AD矩阵的列数 nAD（2）是列数
% compute horizons
 N1 = d +1;
 N2 = d + N;
% Compute F（z^-1）
% inilialization vector f
f（1:）= [1 zeros（1nAD-2）];                                                %  zeros（1nAD-2）创建一个一行A矩阵列数减一的全零阵   ，加上前边的1以后即f（1:）= [1 0 0.....（na-1项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）;                                     %  f（j+1i） = f（ji+1）-f（j1）*AD（i+1）;书中F矩阵求取的递推公式
end                                                                        % f（j+1nAD-1） = -f（j1）*AD（nAD）最后一列。
   f（j+1nAD-1） = -f（j1）*AD（nAD）;                                         %虽然f（1:）= [1 zeros（1nAD-2）]; 初始化了但是在该部计算中没有对结果产生影响，因为系统固有的计算周期原因d>=1，
                                                                           %所以矩阵f（2，：）可计算出书中17页F1（z^-1）的初始值，之后由递推公式可以求出整个f矩阵值。
end
F = f（1+N1:1+N2:）;                                                        %得到了去除纯迟延后的F矩阵。也就说了F和纯迟延无关。                                           
% Computes E（z^-1）
E = zeros（N2）;                                                             %行向量E
e（1） = 1;               % for the special case  1/~A                       %E（z^-1）的初值
E（11） = e（1）;
for i = 2: N2
    e（i） = f（i1）;                                                         %e是f的第一列
    E（i1:i）=e;
end
E = E（N1:N2:）;

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

     文件       2006  2015-08-03 15:35  gpc\diophantine.asv

     文件       2079  2015-08-04 19:56  gpc\diophantine.m

     文件        953  2015-08-07 10:43  gpc\getG.m

     文件       1210  2015-08-06 22:44  gpc\getGij.m

     文件       5106  2015-08-17 12:15  gpc\gpc.asv

     文件       5461  2015-08-18 16:42  gpc\gpc.m

     文件         85  2015-08-04 19:56  gpc\rankn.m

     文件        836  2015-08-04 19:56  gpc\staten.m

     文件       1034  2015-08-03 14:30  gpc\tf2AB.asv

     文件       1229  2015-08-06 21:23  gpc\tf2AB.m

     目录          0  2015-09-03 12:59  gpc

     文件        832  2015-08-06 22:46  gpc\getG.asv

     文件        742  2015-08-06 22:41  gpc\getGij.asv

     文件       3760  2015-08-06 22:15  gpc\gpctest.m

     文件       2982  2015-08-31 13:13  gpc\RLS.m

     文件       3052  2015-08-31 13:46  gpc\simullink\RLS.asv

     文件       3052  2015-08-31 13:51  gpc\simullink\RLS.m

     文件      21111  2015-08-31 13:51  gpc\simullink\testgpc.mdl

     文件      19902  2015-09-01 13:28  gpc\zzz.mdl

     目录          0  2015-09-03 12:59  gpc\simullink

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

                75432                    20