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）;%所求AD既是A~的參數?AD=（1，a1~a2~...ana~）
%?note?that?nAD?=?n~a?+?1
nAD?=?size（AD）;
nAD?=?nAD（2）;%nAD是A~參數的個數
%?compute?horizons
?N1?=?d?+1;
?N2?=?d?+?N;%在這里N2相當于P?
%?Compute?F（z^-1）
%?inilialization?vector?f??！！！！！f是na次多項式
f（1:）=?[1?zeros（1nAD-2）];%nAD-2的意思是nAD?已經是階數+1，第一行第一列已有，是1，所以減2f的第一行是（1?0?0?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）;
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:）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????450??2009-09-24?09:58??getG.m
?????文件?????????107??2017-11-28?15:15??getGii.m
?????文件????????2084??2017-12-05?15:25??gpc.m
?????文件??????????85??2009-09-25?11:10??rankn.m
?????文件?????????836??2009-09-24?09:57??staten.m
?????文件?????????577??2009-09-24?09:57??tf2AB.m
?????文件?????????945??2017-11-27?22:41??diophantine.m