AR模型预测程序

function [Porder d1]=ARorder（testsigp1 p2 ）
%求解AR模型的阶数及模型系数
%输入tests为时间序列为行向量，p1为阶数起点，p2为阶数终点 
%输出Porder是模型阶数，d1为模型参数

%%%确定AR模型阶数
k=1;
for p=p1:p2
m=armax（testsig[p 0]）;