2021-11-07

%价格的预测
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