第9章 马尔可夫预测方法.pptVIP

【例9-3】 已知2008年3月27日―5月30日上证A股收盘价指数数据如表9-2，试用加权马氏链法预测下一个交易日的指数范围。 3 3602.66 44 5 3816.5 33 3 3645.58 22 3 3664.84 11 2 3568.88 43 5 3837.74 32 4 3732.85 21 3 3642.65 10 3 3629.33 42 4 3735.7 31 4 3759.61 20 2 3581.82 9 2 3541.75 41 5 3805.75 30 1 3439.85 19 4 3790.69 8 2 3530.58 40 4 3791.52 29 1 3303.09 18 4 3777.25 7 3 3644.39 39 5 3837.08 28 1 3270.77 17 3 3616.38 6 3 3657.49 38 4 3755.65 27 1 3247.44 16 2 3512.97 5 4 3718.98 37 5 3917.7 26 1 3381.68 15 1 3492.78 4 3 3612.64 36 5 3946.45 25 1 3453.86 14 3 3643.24 3 4 3782.37 35 5 3875.22 24 2 3513.24 13 4 3756.12 2 5 3802.85 34 3 3696.74 23 1 3459 12 2 3578.92 1 状态 收盘价 序号 状态 收盘价 序号 状态 收盘价 序号 状态 收盘价 序号 （1）先将价格指数按从小到大划分为五个区间：低于3500的为状态1，3500(含)至3600的为状态2，3600(含)至3700的为状态3，3700(含)至3800的为状态4，3800(含)以上的股价为状态5；相应地找出每个数据所在的状态， （2）计算各阶自相关系数 rk=autocorr(X,8) %rk= 1.0000 0.7844 0.5760 0.4516 0.3584 0.2347 0.0511 -0.0977 取滞后5期 （3）将各阶自相关系数规范化可以得到各阶滞时的马尔可夫链权重 wk=rk(2:6)./sum(rk(2:6)) （4）不同滞时的马尔可夫链的转移概率矩阵： 1）滞后一步转移矩阵 程序类似例9-2中的程序，只是将区间范围值替换即可 2）滞后二步转移矩阵 程序同1），只是将程序中的“＋1”改为“＋2”，注意一下最后两个数据不转换 3）滞后三步转移矩阵 程序同1），只是将程序中的“＋1”改为“＋3”，注意一下最后三个数据不转换 4）滞后四步转移矩阵 程序同1），只是将程序中的“＋1”改为“＋4”，注意一下最后四个数据不转换 5）滞后五步转移矩阵 程序同1），只是将程序中的“＋1”改为“＋5”，注意一下最后五个数据不转换 （5）选定最后五天为初始日（5月26日至30日），其所在的状态结合其相应的转移概率矩阵即可预测出该时段证券价格指数的状态概率 0.0743 0.1051 0.3588 0.2571 0.2046 （加权和） 0.1000 0.2000 0.3000 0.2000 0.2000 0.3261 1 3 5-30 0 0.1667 0.5000 0.1667 0.1667 0.2395 2 2 5-29 0.2222 0 0.3333 0.3333 0.1111 0.1878 3 3 5-28 0 0 0.2000 0.6000 0.2000 0.1490 4 2 5-27 0 0 0.5000 0 0.5000 0.0976 5 2 5-26 5 4 3 2 1 状态空间概率 权 重 滞时/天 状 态 初始 日 （6）将同一状态的各预测概率加权和作为证券价格指数处于该状态的预测概率 将表中“状态空间概率”用矩阵Q表示： Q=[P5(2,:);P4(2,:);P3(3,:);P2(2,:);P1(3,:)] wk=wk p1=sum(wk(5:-1:1).*Q(:,1)) p2=sum(wk(5:-1:1).*Q(:,2)) p3=sum(wk(5:-1:1).*Q(:,3)) p4=sum(wk(5:-1:1).*Q(:,4)) p5=sum(wk(5:-1:1).*Q(:,5)) 结果为： p1 = 0.2046；p2 = 0.2571；p3 =0.3588；p4 =0.1051；p5 =0.0743 9.2.4 期望利润预测 在企业当中，除需要摸清销路的变化情况外，还要对利润的变化进行预测。如：某商品的销售状态有畅销和滞销两种状态，则通过