请教如何根据BIC图来判断ARMA模型的p，q值
展开全部
%下面要对差分以后的序列进行拟合和预测，求出最好的阶数
z=[DX;zeros(12,1)];
z=iddata(z);
test=[];
for p=1:12
for q=1:12
m=armax(z(1:200),[p q]);
AIC=aic(m);
test=[test;p q AIC];
end
end
for k=1:size(test,1)
if test(k,3)==min(test(:,3))
p_test=test(k,1);
q_test=test(k,2);
break;
end
end
%拟合
m1=armax(z(1:200),[p_test q_test]);
figure(5);
e = resid(m1,z);
plot(e);
set(gca,'Xlim',[0 ls]);
figure(6);
subplot(2,1,1)
autocorr(e.outputdata)
subplot(2,1,2)
parcorr(e.outputdata)
set(gca,'Xlim',[0 ls]);
%预测过程
pr=predict(m1,z,12);
po=pr.outputdata;
figure(7)
plot(po,'r')
hold on
plot(y,'b');
set(gca,'Xlim',[0 ls]);
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
