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