function?[Porder?d1]=ARorder（testsigp1?p2?）
%求解AR模型的階數(shù)及模型系數(shù)
%輸入tests為時(shí)間序列為行向量，p1為階數(shù)起點(diǎn)，p2為階數(shù)終點(diǎn)?
%輸出Porder是模型階數(shù)，d1為模型參數(shù)

%%%確定AR模型階數(shù)
k=1;
for?p=p1:p2
m=armax（testsig[p?0]）;