带模型辨识的广义预测控制matlab程序

虽然是单输入单输出的GPC程序，但是我作了很多注释，大家可以先学着

function [E  F  G H] = my_diophantine（ABPM）%P为预测时域  E:P*1  F:（na+1）*P  G:P*M H:nb*P
%A=[a1a2...ana]‘
%B=[b0b1b2...bnb]‘

na = length（A）;
nb = length（B）-1;

%--------递推矩阵FT（na+1）*（na+1）的输入--------------
FT = zeros（na+1）;
for i = 1:na
    FT（ii+1） = 1;
end
for i = 2:na
    FT（i1） = A（i-1）-A（i）;
end
FT（11） = 1-A（1）;
FT（na+11） = A（na）;
%---------FT是用于计算序列F时的递推矩阵--------------

E=zeros（P1）;
F=zeros（na+1P）;
E（11）=1;
F（:1）=FT（:1）;

for i=2:P
    E（i1）=F（1i-1）;
    F（:i）=FT*F（:i-1）;
end

%-------------递推矩阵HTnb*nb的输入----------------- 
HT = zeros（nbnb）;
for i=1:（nb-1）
    HT（ii+1） = 1;
end
%-------HT是用于计算序列H时的递推矩阵---------------
g = zeros（P1）;
H = zeros（nbP）;
g（1） = E（1）*B（1）;
H（:1） = B（2:nb+1）;

for i=2:P
    g（i） = E（i）*B（1） + H（1i-1）;
    H（:i） = HT*H（:i-1）+E（i）*H（:1）;
end

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

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

     文件       2095  2011-03-02 10:55  my_main.m

     文件       1037  2011-02-19 11:52  my_diophantine.m

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

                 3132                    2