第1个回答 2015-06-19
用matlab确定ARMA模型的阶数,要对差分以后的序列进行拟合和预测,求出最好的阶数。代码如下供参考:
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]);