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