%价格的预测
for i = 1 : 30
ahat = AR3(price);%preice改为最后一天最后两百个数组(在最后数据会不变,忽略这些数据)
price(end+1) = ahat(end);%把预测的数加到现有的数组
end
function [ahat] = AR3(array)
%AR3 AR(n)预测后一项值
% n为AR(n)的参数设计
array = array';array = array(:);array= array';
array1 = diff(array);
array2 = diff(array1);
d = ar(array2,81,'ls'); %,用最小二乘法估计模型的参数,得到多项式d,n=81
bhat = predict(d,[array2';0;0]);%多项式的d进行预测,样本数据必须为列向量,预测结果为price31-brece34504
ahat = [array(1),array+bhat'] ;
delta = abs((ahat(1:end - 1) -array)./array);
end
%
for i = 1 : 30
ahat = AR4(volume10);%preice改为最后一天最后两百个组成数组(在最后数据会变为0,忽略这些数据)
volume10(end+1) = ahat(end);%把预测一个数加到现有的数组
end
function [ahat,delta] = AR4(array)
%AR4 AR(n)预测后一项值
% n为AR(n)的参数设计
array = array';array = array(:);array= array';
d = ar(array,20,'ls'); %,用最小二乘法估计模型的参数,得到多项式d,n=81
bhat = predict(d,[array']);%多项式的d进行预测,样本数据必须为列向量,预测结果为price31-brece34504
ahat = [array(1),array+bhat'] ;
delta = abs((ahat(1:end - 1) -array)./array);
end