上证综指时间段为2014年4月1号到2015年3月31号的收盘价灰色预测.docVIP

上证综指时间段为2014年4月1号到2015年3月31号的收盘价灰色预测 代码如下： function []=greymodel(y) % 本程序主要用来计算根据灰色理论建立的模型的预测值。 % 应用的数学模型是 GM(1,1)。 % 原始数据的处理方法是一次累加法。 y=input(请输入数据 ); n=length(y); yy=ones(n,1); yy(1)=y(1); for i=2:n yy(i)=yy(i-1)+y(i); end B=ones(n-1,2); for i=1:(n-1) B(i,1)=-(yy(i)+yy(i+1))/2; B(i,2)=1; end BT=B; for j=1:n-1 YN(j)=y(j+1); end YN=YN; A=inv(BT*B)*BT*YN; a=A(1); u=A(2); t=u/a; i=1:n+2; yys(i+1)=(y(1)-t).*exp(-a.*i)+t; yys(1)=y(1); for j=n+2:-1:2 ys(j)=yys(j)-yys(j-1); end x=1:n; xs=2:n+2; yn=ys(2:n+2); plot(x,y,^r,xs,yn,*-b); title(预测结果和实际值比较图); det=0; sum1=0; sumpe=0; for i=1:n sumpe=sumpe+y(i); end pe=sumpe/n; for i=1:n; sum1=sum1+(y(i)-pe).^2; end s1=sqrt(sum1/n); sumce=0; for i=2:n sumce=sumce+(y(i)-yn(i)); end ce=sumce/(n-1); sum2=0; for i=2:n; sum2=sum2+(y(i)-yn(i)-ce).^2; end s2=sqrt(sum2/(n-1)); c=(s2)/(s1); disp([后验差比值为：,num2str(c)]); if c0.35 disp(系统预测精度好) else if c0.5 disp(系统预测精度合格) else if c0.65 disp(系统预测精度勉强) else disp(系统预测精度不合格) end end end disp([下个拟合值为 ,num2str(ys(n+1))]); disp([再下个拟合值为,num2str(ys(n+2))]); 模拟案例(1) greymodel(y) 请输入数据 [2.874,3.278,3.337,3.390,3.679] 后验差比值为：0.23876 系统预测精度好 下个拟合值为 3.7507 再下个拟合值为3.8928 模拟案例(2) greymodel(y) 请输入数据 [724.57, 746.62, 778.27, 800.8, 827.75,871.1, 912.37, 954.28, 995.01, 1037.2] 后验差比值为：0.067811 系统预测精度好 下个拟合值为 1079.3804 再下个拟合值为1125.6546 实际案例分析 短期预测 上证综指(000001)，时间段为2015年3月10号到3月31号的收盘价，网址/trade/lsjysj_zhishu_000001.html?year=2015season=1 greymodel(y) 请输入数据 [3302.41 3286.07 3290.90 3349.32 3372.91 3449.31 3502.85 3577.30 3582.27 3617.32 3687.73 3691.41 3660.73 3682.10 3691.10 3786.57 ] 后验差比值为：0.26072 系统预测精度好 下个拟合值为 3835.8636 再下个拟合值为3873.8142 实际值：2015/4/1 为3810.29； 预测精度：3835.8636/3810.29=100.6712% ，偏差+0.6712%。 实际值：2015/4/2 为3825.78； 预测精度：3873.8142/3825.78=101.2555%，偏差+1.2555%。 长期预测 上证综指(000001)，时间段为2014年4月1号到2015年3月31号的收盘价，网址/trade/lsjysj_zhishu_000001.html?year=2015season=1。 greymodel(y) 请输入