《经济预测与决策技术及MATLAB实现》课件 第10、11章 灰色测方法、景气预测法.ppt

4）计算合成指数首先计算初始指数的时间序列，即环比原始指数，t=2,3,...,m表示期数。用表示所选基准年份的平均值，由此得到合成指数为：11.1.3景气综合评分——预警系统1．预警系统原理是选择一组反映经济发展状况的敏感指标，运用有关的数据处理方法，将多个指标合并为一个综合性指标，通过一组类似于交通管制信号红、黄、绿灯的标识，利用这组指标和综合指标对当时的经济状况发出不同的信号，通过观察信号的变动情况，来判断未来经济增长的趋势。2．预警指标状态区域划分和临界点确定（1）区域划分：将判断区域分为过热、偏热、正常、偏冷和过冷五个域，分别以红灯、黄灯、绿灯、浅蓝灯和蓝灯表示。绿灯区居中，代表常态区或稳定区，概率一般定为50%；红灯区和蓝灯区属于极端区，经济含义为过热和过冷，概率各占10%；黄灯区和浅蓝灯区为相对稳各占15%。clearX0=[7.188.609.509.6010.26 10.98 12.02 13.22 14.1414.76];%产生累加生成列X1X1=cumsum(X0)%构造数据矩阵B和数据向量Yfork=2:length(X0)z(k)=(1/2)*(X1(k)+X1(k-1));endzB=[(-z(2:end))ones(length(z)-1,1)]Y=(X0(2:end))%计算模型参数bata=inv(B*B)*B*Y;a=bata(1)b=bata(2)c=-a*(X0(1)-b/a)%根据模型求出原始数据预测序列值X2X2(1)=X0(1)fork=1:length(X0)-1X2(k+1)=(-a)*(X0(1)-b/a)*exp(-a*k);endX2%计算残差序列和相对残差序列E=X0-X2F=E./X0E=0-0.20900.0518-0.5338-0.6092-0.6779-0.4838-0.1912-0.2443-0.6681F=0-0.02430.0055-0.0556-0.0594-0.0617-0.0403-0.0145-0.0173-0.0453残差序列和相对残差序列的后面7项符号一致，都为负值，且误差相对较大，因此可以使用残差修正模型提高其精度。（2）下面利用残差对上述模型进行修正E0=-E(4:end)%产生累加生成列E1E1=cumsum(E0);%构造数据矩阵B2和数据向量Y2fork=2:length(E0)z2(k)=(1/2)*(E1(k)+E1(k-1));endB2=[(-z2(2:end))ones(length(z2)-1,1)];Y2=(E0(2:end));%计算模型参数a2,b2bata2=inv(B2*B2)*B2*Y2;a2=bata2(1)b2=bata2(2)c2=-a2*(E0(1)-b2/a2)%根据修正模型预测X3X3(1)=X0(1);fork=1:length(X0)-1ifk3X3(k+1)=-a*(X0(1)-b/a)*exp(-a*k);%修正之前模型elseX3(k+1)=-a*(X0(1)-b/a)*exp(-a*k)-(-a2)*(E0(1)-b2/a2)*exp(-a2*k);%修正之后模型endendX3%计算修正之后残差序列和相对残差序列E2=X0-X3F2=E2./X0E2=0-0.20900.0518-0.0631-0.1810-0.2885-0.12970.13100.0487-0.4016F2=0-0.02430.0055-0.0066-0.0176-0.0263-0.01080.00990.0034-0.0272（3）对2022至2024年进行预测。k=length(X0):length(X0)+2X3(k+1)=-a*(X0(1)-b/a)*exp(-a*k)-(-a2)*(E0(1)-b2/a2)*exp(-a2*k)X3=7.18008.80909.44829.663110.441011.268512.