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

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