《经济预测与决策技术及MATLAB实现》课件 第11章 景气预测法

经济预测与决策技术及MATLAB实现第11章景气预测法11.1景气预测基本理论11.2案例分析练习与提高（十一）

11.2.1国房景气指数11.2.2上海房地产景气指数第11章景气预测法11.1景气预测基本理论11.1.1景气指标体系基本概念景气预测是通过正确地选择景气指标体系来实现的经济的景气状态，是通过一系列经济指标来描述的，称为景气指标。景气指标是从众多的经济指标中挑选出来的，分为先行指标、同步指标和滞后指标三类。景气循环又称经济波动，也称经济周期。一个标准的经济周期，通常包括扩张和收缩两个时期，分为四个阶段：复苏、高涨、衰退和萧条。11.1.2景气循环法预测过程景气循环法预测步骤（1）确定时差关系的参照系——基准循环，这是关键的一步；（2）选择景气指标；（3）划分先行、同步、滞后指标；（4）对先行、同步、滞后指标分别编制扩散指数和合成指数；（5）对计算结果进行分析，了解当前经济状况，预测未来经济波动。1．基准循环的确定方法（1）以重要的经济指标（GNP、GDP、工业总产值等）的周期为基准循环；（2）专家意见及专家评分；（3）经济大事记和经济循环年表；（4）初选几项重要指标计算历史扩散指数；（5）以一致合成指数转折点为基础。2．景气指标的筛选（1）聚类分析：是将要分析的m个指标变量各自看作一类，然后计算各类之间的关系密切程度（相关系数或距离），并将关系最密切的两类归为一类，把另外一些彼此之间相似程度大的指标聚合为另一类。（2）非参数检验：检验同类指标之间是否具有显著性差异。（3）指标筛选流程

对搜集到得所有指标，进行分类，分别在每一类中进行聚类分析，并将聚类结果中单独一类的指标直接纳入预警指标体系。对于指标个数的大于两个的子类，则进行显著性差异检验。如果有显著性差异，则对此类中的指标继续进行聚类分析。直到同类的指标不再有显著性差异为止如果同类指标没有显著性差异，对于指标个数为2的子类，通过主观判断，从二者选择一个纳入预警指标体系；对指标个数大于2的子类，对指标进行相关分析(求指标间的相关系数)，并对每一指标的相关系数求和，数值最大者，认为它在此子类中与同类其它指标具有较高的相关性，因此代表性较强，纳入预警指标体系。3．景气指标的划分（1）时差相关分析.是利用相关系数验证时间序列先行、一致或滞后关系的一种方法。计算方法是以一个重要的能够敏感反应当前市场活动的指标作为基准指标，然后使被选择指标超前或滞后若干期，计算它们的相关系数。l:为时差或延迟数，取负数表示超前，取正数表示滞后，L为最大延迟数，是数据取齐后的数据个数。在选择景气指标时，取最大的时差相关系数，作为反应被选指标与基准指标超前、滞后或同步的时差相关关系，相应l

的表示超前或滞后期。（2）K-L信息量K-L信息量，用以判定两个概率分布的接近程度。设基准指标做归一化处理得到序列为p，被选择指标归一化处理得到序列为q，则K-L信息量计算公式：l:为时差或延迟数，取负数表示超前，取正数表示滞后，L为最大延迟数，是数据取齐后的数据个数。在选择景气指标时，取最小的K-L信息量，作为反应被选指标与基准指标超前、滞后或同步的时差相关关系，相应l

的表示超前或滞后期。K-L信息量越小，越接近于零，说明指标x与基准指标y越接近。一般情况下，当K-L信息显小于50时，就可以认为两个指标之间的概率分布较为接近。4．景气指数计算方法（1）扩散指数（DI）:是（先行、一致或滞后）指标组内在某个时期扩张（上升）指标个数占组内所采用指标个数的比率。DI(t)为t时刻的扩散指数；Xi(t)为第i个变量指数在t时刻的波动测定值；Wi为第i个变量指标分配的权数；n为变量指标总数；s为两比较指标值的时间差。从扩散指数的定义可知，DI(t)始终处于0到100％之间当0<DI(t)<50%：经济向扩张方向运动，经济系统运行于不景气空间后期；当50%<DI(t)<100%：经济运行于景气空间，经济状况发生重大转折，随着DI(t)向100%不断逼近，经济越来越热当100%>DI(t)>50%：经济处于景气空间后期，经济系统处于降温阶段；当50%>DI(t)>0：经济运行发生重大转折，经济系统处于全面收缩阶段，进入一个新的不景气空间前期。利用扩散指数可预测和监控经济运行，能有效分析和预测经济波动的转折点，但不能说明经济波动的程度。（2）合成指数又称景气综合指数，是由一类特征指标以各自的变化幅度为权数的加权综合平均数，即以多个指标的加权平均。计算方法1）单指标的对称变化率和标准化对称变化率就是本期比上期的增长量与两期的平均值比较得出的数值。Cit中表示第i个指标第t年的对称变化率，用百分数表示。和分别表示第i个指标第t年和第t-1年的原始数据。用Ai表示第i个指标对称变化率时间序列的序时平均数，N表示标准化的期数，则有：用Sit表示第i个指标t期对称变化率的标准化数值，则有2）多指标对称变化率标准化后的加权平均数Rt:多指标综合的平均对称变化率，Sit表示第i个指标第t年的对称变化率的标准化数值，Wi表示第i个指标的权数，i=1,2,...,k表示指标项目数。3）平均变化率以同步指数标准化标准化因子F的计算公式是Rt表示超前指标或滞后指标的综合平均对称变化率第t期的数值，Pt表示一致指标的综合平均对称变化率时间序列的t期数据，N表示期数。平均变化率以同步指数标准化的计算公式是：t=2,3,...,m表示期数。4）计算合成指数首先计算初始指数的时间序列，即环比原始指数，t=2,3,...,m表示期数。用

表示所选基准年份的平均值，由此得到合成指数为：11.1.3景气综合评分——预警系统1．预警系统原理是选择一组反映经济发展状况的敏感指标，运用有关的数据处理方法，将多个指标合并为一个综合性指标，通过一组类似于交通管制信号红、黄、绿灯的标识，利用这组指标和综合指标对当时的经济状况发出不同的信号，通过观察信号的变动情况，来判断未来经济增长的趋势。2．预警指标状态区域划分和临界点确定（1）区域划分：将判断区域分为过热、偏热、正常、偏冷和过冷五个域，分别以红灯、黄灯、绿灯、浅蓝灯和蓝灯表示。绿灯区居中，代表常态区或稳定区，概率一般定为50%；红灯区和蓝灯区属于极端区，经济含义为过热和过冷，概率各占10%；黄灯区和浅蓝灯区为相对稳各占15%。（2）临界值的确定：先给出各种灯区的分界线检查值，再求出综合指标的各时期的评分数，最后将这些综合评分数与各灯区的分界线对应。如假定参与综合的指标数为M，每个指标的满分为5分，则总满分为5M。按惯例取满分的85％为红灯区与红黄灯区的分界线。其总分数为5M×85％=4.25M分，满分的75％和50％为绿灯区的上、下分界线，它对应的分数为5×M×75％=3.75M和5×M×50％＝2.5M，满分的40％为浅蓝灯与蓝灯的分界线，对应的分数为5×M×40％＝2M。3．计算景气综合评分将每类指标所包含的各指标所处灯区对应的分值加总，得到景气综合评分，并判断综合评分落在哪个区域，显示哪种灯色。根据综合指数的信号判断宏观经济的运行状况，根据各预警指标信号分析导致综合指数信号处于当前状态的具体原因，为宏观调控提供参考和建议。11.2案例分析11.2.1国房景气指数一、我国国房景气指数体系

国房景气指数是全国房地产开发景气指数的简称，其过程为：1．确定指标体系计算国房景气指数的指标包括：（1）土地出让收入指数；（2）本年完成开发土地面积指数（3）房地产开发投资指数；（4）本年资金来源指数；（5）商品房销售价格指数；（6）新开工面积指数；（7）房屋竣工面积指数；（8）空置面积指数。分析报告用的指数指标包括：（1）新开工面积指数；（2）房屋施工面积指数；（3）国有单位投资指数；（4）国内贷款指数；（5）利用外资指数；（6）自筹资金指数；（7）住宅销售价格指数；（8）办公楼销售价格指数；（9）商业营业用房销售价格指数；（10）个人商品房销售额指数；（11）竣工房屋价值指数；（12）住宅空置面积指数；（13）办公楼空置面积指数；（14）商业营业用房空置面积指数。2．建立指标数据库及数据调整对指标体系确定的各项指标，建立历史资料数据库。有些指标如房地产销售价格指标必须先计算后才能得到。有些指标如资金来源、空置面积等必须先进行推算，以保证有必要的数据长度。3．消除量纲影响的方法由于各指标的计量单位不同，不能直接相加综合平均，必须采用消除量纲影响的方法来对数据进行同度量处理，同时，也满足了合成指数编制中对指标进行标准化处理的要求。国房景气指数采用功效系数法，消除量纲的影响比较符合国房景气指数的特点。4．权数的确定

确定参加计算的8个指标的权数。根据我国房地产开发业的发展状况，征求部分专家意见初步确定了国家级各分类指数权数值。分类指数权重代号权重取值性质房地产开发投资指数W120-30同步资金来源指数W2<10同步土地出让收入指数W3<10超前土地开发面积W410-20超前新开工面积指数W5<10同步竣工面积指数W6<10滞后空置面积指数W710-20滞后商品房销售价格指数W810-20同步5．确定基准对比时期分析我国房地产特点，确定1995年3月份为国房指数的基准对比期。6．季节和价格因素调整季节因素调整：用X-11程序进行季节因素调整。价格因素调整：房地产开发投资、到位资金、土地转让收入和销售价格4个指标涉及到价格调整问题。7．指数计算模型（1）分类指数的计算分类指数（K1-K8）的计算采用了环比指数的计算方法：Z(t)表示报告期t期各项指标数据，Z(t-1)表示上一年同期各项指标数据。（2）国房指数的计算在计算8个分类指标后，采用固定权重的算术加权平均法，得到综合指数

将综合指数与基准综合指数对比，得到国房指数G：对空置面积分类指数K7，由于空置面积属于逆指标，其指数值越小，表明景气状况越好，因此，在计算综合指数时，K7是以倒数形式出现。8．计算结果和景气评价根据上述流程和确定的权数，计算国房景气指数值说明房地产业的发展状况。国房指数以100为景气分界点，100以上为景气空间，100以下为不景气空间。为更好发挥国房指数的景气监测作用，应设计合理预警界限，当国房指数突破某预警界限时，系统能发出“报警”信号。1．景气区域的划分过热区：表示房地产发展已经过热，必须降温，该区为极端区趋热区：也称偏热区，表示房地产发展有过热倾向，需引起关注，该区为警告区；适度区：表示房地产发展适度，该区为稳定区；趋冷区：也称偏冷区，表示房地产萎缩倾向，需引起关注，该区为警告区；过冷区：表示房地产处于萧条状态，需要扶持，该区为极端区二、国房指数预警系统根据我国房地产发展趋势的判断，我们认为“适度”区的落地点概率为80%，“偏热”区和“偏冷”区的落地点概率为7.5%，“过热”区和“过冷”区的落地点概率为2.5%的比例比较合理。根据

原理，即假定国房指数服从正态分布：2．分界限的确定首先确定指数波动的中心线，以此作为“适度“区地中心。选取国房指数样本，用其平均值来估计均值和均方差“适度”区地上、下限为：“偏热”区与“过热”区的分界限为“偏冷”区与“过冷”区的分界限为3．对各序列在各时期所处的景气状态标上各色信号。分别以红灯、黄灯、绿灯、浅蓝灯和蓝灯表示。【例11-1】给出2000年至2015年的我国国房景气指数数据如表11-2所示，其中2008年（含）以前每年国家公布12个月数据，从2009年开始，每年1月份数据不公布，即每年公布11个月数据，试划分出我国国房景气指数在这期间的“过热”、“偏热”、“适度”、“偏冷”和“过冷”五个区域，以及景气和不景气空间，来说明我国的房地产发展状况。年\月1234567891011122000102.2102.6101.8101.9102.2102.7102.9103.3103.5103.4103.7104.02001105.2105.0106.2106.9106.5106.0106.4106.4105.9105.6105.1105.32002103.9103.5103.9104.4104.6104.7104.3104.5105.0105.2104.9104.22003108.9109.1107.3106.9106.3107.0107.3106.9106.7106.7106.5106.22004106.6106.6106.1105.5105.0104.8104.8104.9105.0104.7104.2103.72005103.1102.7102.4102.3102.2102.1102.0101.8101.4101.0100.7100.62006100.8101.1101.5101.6101.9102.9103.5103.3103.1103.4103.9103.02007102.4101.8101.2102.7103.3103.6104.0104.5105.0105.7106.6106.52008106.1105.6104.7104.1103.3103.1102.4101.8101.299.798.596.52009

94.994.794.895.996.698.0100.1101.1102.0102.8103.72010

105.5105.9105.7105.1105.1104.7104.1103.5103.6103.2101.82011

102.9103.0103.2103.2101.8101.5101.1100.4100.399.998.92012

97.996.995.694.994.794.694.694.494.695.795.62013

97.997.697.497.397.397.497.397.396.996.497.22014

96.996.495.895.094.894.894.894.794.894.393.92015

93.893.192.692.492.693.093.593.493.393.493.3X=[…];%表数据，共185个，按时间序列排序plot(X)m=mean(X)%均值为适度中心线s=std(X)%标准差Y1=m-1.96*s%偏冷区与过冷区分界线Y2=m-1.28*s%适度区下限Y3=m%适度区中心线Y4=m+1.28*s%适度区上限Y5=m+1.96*s%偏热区与过热区分界线holdon%绘制五个区域图x=1:length(X);Z1=Y1.*ones(size(x));Z2=Y2.*ones(size(x));Z3=Y3.*ones(size(x));Z4=Y4.*ones(size(x));Z5=Y5.*ones(size(x));plot(x,Z1,'-b',x,Z2,'-c',x,Z3,':k',x,Z4,'-k',x,Z5,'-r')axis([0,200,88,115])text(60,112,'过热区')text(60,108,'偏热区')text(60,105,'适度区')text(60,98,'适度区')text(60,95,'偏冷区')text(60,91,'过冷区')text(50,102,'适度中心线’)%绘制景气不景气空间图Z6=100.*ones(size(x));plot(x,Z6,'-.')text(160,100,'景气分界线')text(160,103,'景气空间')text(160,98,'不景气空间')t1=[0:8]*12;%2008年前每年12个月数据t2=[10:17]*11;%2009年开始每年11个月数据set(gca,'XTick',[t1,t2])set(gca,'XTickLabel',{'2000';'2001';'2002';'2003';'2004';'2005';'2006';'2007';...'2008';'2009';'2010';'2011';'2012';'2013';'2014';'2015';'2016'})set(gca,'YTick',[93.296.1100101.5106.9109.8])set(gca,'YTickLabel',{'93.2';'96.1';'100';'101.5';'106.9';'109.8’})xlabel('年份')ylabel('景气指数')holdoff11.2.2上海房地产景气指数【例11-2】现已搜集上海房地产市场有关的指标数据（见表11-3），试求：（1）使用聚类分析和假设检验等方法筛选景气指数；（2）利用时差相关性分析判断超前指标、同步指标和滞后指标；（3）计算超前指标、同步指标的合成指数，说明上海房地产市场的景气发展状况。

一、选取房地产统计指标并收集资料数据1．经过上海房地产市场的数据收集，建立上海房地产市场21指标数据库国内生产总值(x1)，房地产业增加值(x2)，房地产开发投资额(x3)，建筑业总产值(x4)，社会消费品零售总额(x5)，人均可支配收入(x6)，房地产业各项贷款余额(x7)，房地产开发企业资金(x8)，供应土地占地面积(x9)，房地产从业人数(x10)，房屋竣工面积(x11)，房屋新开工面积(x12)，房屋施工面积(x13)，房屋出租面积(x14)，房地产开发企业总收入(x15)，商品房销售面积(x16)，住宅商品房销售额(x17)，新建商品住宅销售价格指数，固定资产投资价格指数(x19)，居民消费价格指数(x20)，住宅房屋销售价格(x21)2．房地产统计指标分类总量指标（x1-x4），购买力指标（x5~x7），生产类指标（x8-x13）、交易类指标（x14-x19）和价格类指标（x19~x21）二、筛选景气指标以生产类指标：房地产开发企业资金(x8)，供应土地占地面积(x9)，房地产从业人数(x10)，房屋竣工面积(x11)，房屋新开工面积(x12)，房屋施工面积(x13)共6个指标为例进行分析。年份房地产开发企业资金(x8)/百亿元供应土地占地面积(x9)/百万m2房地产从业人数(x10)/万人房屋竣工面积(x11)/百万m2房屋新开工面积(x12)/百万m2房屋施工面积(x13)/百万m2201032.2924.859.1419.4130.31112.95201134.4431.487.6522.4136.44129.83201239.6922.317.4223.0527.24132.50201350.9322.247.1022.5427.06135.17201452.7018.626.8623.1327.82146.90201555.3210.006.4626.4726.05150.95201664.098.916.3625.5128.41151.11201753.8510.256.3833.8826.18153.62201853.3013.086.4631.1626.87146.72201954.4513.136.3826.7030.63148.03202055.1622.645.8928.7834.41157.40在进行聚类以前，应首先对这6个指标进行无量纲处理，为了保持原有指标的变异特性，选择归一化处理：A=[32.29 24.85 9.14 19.41 30.31 112.9534.44 31.48 7.65 22.41 36.44 129.8339.69 22.31 7.42 23.05 27.24 132.5050.93 22.24 7.10 22.54 27.06 135.1752.70 18.62 6.86 23.13 27.82 146.9055.32 10.00 6.46 26.47 26.05 150.9564.09 8.91 6.36 25.51 28.41 151.1153.85 10.25 6.38 33.88 26.18 153.6253.30 13.08 6.46 31.16 26.87 146.7254.45 13.13 6.38 26.70 30.63 148.0355.16 22.64 5.89 28.78 34.41 157.40];%归一化处理B=sum(A);X=[A(:,1)/B(1),A(:,2)/B(2),A(:,3)/B(3),A(:,4)/B(4),A(:,5)/B(5),A(:,6)/B(6)];%将X进行转置，进行聚类分析Y=pdist(X')%计算指标之间的距离Z=linkage(Y)%形成系统聚类树dendrogram(Z)%显示系统聚类树T=cluster(Z,4)%不同閥值的分类结果对这6个指标，采用聚类法进行聚类结果：分成四类：x11，x13；x8；x10，xl2；x9；由于x8和x9只有一个指标的子类，因此将这两个指标直接纳入预警指标体系。对于x11和x13，x10和x12分别进行显著性检验--秩和ranksum检验，判断指标之间是否具有显著性差异。对x11和x13的秩和ranksum检验：[p1,h1,stats1]=ranksum(X(:,4),X(:,6))运行结果如下：p1=0.5545h1=0检验结果h=0表示x11和x13的总体差别不显著，一定程度上商品房竣工面积比施工面积更能够衡量一定时间点上的商品房供应量．故将商品房竣工面积x11纳入预警指标体系。对x10和x12的秩和ranksum检验：[p2,h2,stats2]=ranksum(X(:,3),X(:,5))运行结果如下：p2=1h2=0检验结果h=0表示x10和x12的总体差别不显著，房地产开发企业从业人数比房屋新开工面积更能说明市场景气现象，因此将房地产开发企业从业人数x10纳入预警指标体系。按照此算法，将总量指标、购买力指标、交易类指标、价格类指标各进行分类指标筛选，最终得出预警指标体系：类型原始指标预警指标总量指标x1~x4国内生产总值x1、房地产业增加值x2、房地产开发投资额x3购买力指标x5~x7全体居民人均可支配收入x6、房地产贷款余额x7生产类指标x8~x13房地产开发企业本年资金x8、供应土地占地面积x9、房地产开发企业平均从业人数x10、房屋竣工面积x11交易类指标x14~x18房屋出租面积x14，商品房销售面积x16、住宅商品房销售额x17价格类指标X19~x21居民消费价格指数x20、住宅房屋平均销售价格x21四、利用时差相关分析划分超前、同步、滞后指标选定基准指标与比较指标，可以直接用原始数据计算超前或滞后若干期的相关系数；也可以先求原始数据的变化率，利用变化率数据再计算超前或滞后若干期的相关系数。下面计算商品房销售面积与其它景气指标的变化率数据，并以最大的时差相关系数为两年来分析。三、确定基准指标计算房地产景气指数的首要工作是确定基准指标，商品房销售面积，是房地产市场供给和需求共同作用的结果，它同房地产市场的兴衰具有直接同步性。%其它景气指标矩阵D=[171.6612.4219.813.0445.4632.2924.859.1419.4112.6220.6123.95104.51.43191.9612.8922.543.4747.6334.4431.487.6522.41 13.0517.9119.82104.51.36201.8215.0123.813.8651.6339.6922.317.4223.0513.5518.9822.09102.211.39218.1819.1728.204.2255.7950.9322.247.122.5412.0623.8232.64102.41.62235.6821.2132.064.6065.6452.7018.626.8623.1311.4220.8529.23102.61.64251.2324.0434.694.9968.4755.3210 6.4626.4712.0224.3143.2102.32.15281.7931.0337.095.4362.3964.098.916.3625.5113.2327.0652.33103.272.59306.3328.7338.575.9072.1953.8510.256.3833.8816.4616.9233.36101.552.49360.1231.8540.336.4282.5053.3013.086.4631.1618.6817.6738.64101.182.9381.5533.0142.316.9493.5054.4513.136.3826.720.3616.9644.57103.783.29387.0133.9346.997.22103.1755.1622.645.8928.7820.5617.8952.69100.13.67];MATLAB程序如下：%商品房销售面积为基准指标y=[20.6117.9118.9823.8220.8524.3127.0616.9217.6716.9617.89]';z=price2ret(y)%基准指标变动率A=price2ret(D)%其它指标变动率[m,n]=size(A);R=[];B=[];%存放相关系数矩阵R和最大相关系数所对应的位置Bfork=1:nx=A(:,k);r1=corrcoef(x(1:end-2),z(3:end));%先行2期r2=corrcoef(x(1:end-1),z(2:end));%先行1期r3=corrcoef(x,z);%同期0r4=corrcoef(x(2:end),z(1:end-1));%滞后1期r5=corrcoef(x(3:end),z(1:end-2));%滞后2期r=[r1(1,2),r2(1,2),r3(1,2),r4(1,2),r5(1,2)];%从先行到滞后共5期相关系数[a,b]=max(abs(r));%找出最大相关系数a和所对应的位置b%其中b=1，表先行2期；b=2，表先行1期；b=3，表同步；b=4，表滞后1期；b=5，表滞后2期R=[R;r];B=[B;b];end指标名称先行2期先行1期同步0期滞后1期滞后2期r最大值位置指标类型国内生产总值0.4559-0.4790-0.0766-0.37150.13282先行房地产行业增加值0.0271-0.54040.81270.0172-0.20063同步房地产开发投资额0.4656-0.04390.31810.3793-0.01961先行人均可支配收入0.40650.2364-0.1510-0.0408-0.05211先行房地产贷款余额0.08790.7178-0.4604-0.14580.12512先行房地产开发企业资金0.1651-0.31910.81420.0764-0.38053同步供应土地占地面积0.8071-0.2325-0.2818-0.0992-0.12461先行房地产从业人数-0.2357-0.1824-0.0412-0.1806-0.15431先行房屋竣工面积-0.26890.2796-0.58350.19800.72485滞后房屋出租面积-0.2552-0.0811-0.5825-0.23890.16523同步商品房销售面积-0.1949-0.30471.0000-0.3047-0.19493同步住宅商品房销售额-0.3478-0.27000.9664-0.3186-0.21993同步居民消费价格指数0.3457-0.39560.11740.2336-0.67095滞后住宅房屋销售价格-0.5227-0.41930.7631-0.3228-0.06343同步基于指标变动率的时差相关分析结果，得到上海房地产市场的景气指标划分如下：（1）先行指标：国内生产总值、房地产开发投资额、人均可支配收入、房地产贷款余额、供应土地占地面积和房地产从业人数，共6个指标。（2）同步指标：房地产行业增加值、房地产开发企业经营资金、房屋出租面积、商品房销售面积、住宅商品房销售额和住宅房屋销售价格，共6个指标。（3）滞后指标：房屋竣工面积和居民消费价格指数2个指标。五、计算各类指标的扩散指数DI扩散指数包括领先指标的扩散指数(LDI)、同步指标的扩散指数(CDI)和滞后指标的扩散指数(DDI)。用扩散指数可判断房地产景气波动的转折点和景气状态%先行指标数据矩阵X1=[…];%同步指标数据矩阵X2=[…];%滞后指标X3=[…];[m1,n1]=size(X1);[m2,n2]=size(X2);[m3,n3]=size(X3);d1=[];d2=[];d3=[];d4=[];d5=[];d6=[];%先行指标扩散指数fori=2:m1;D1=X1(i,:)-X1(i-1,:);d11=length(find(D1>0));d12=length(find(D1==0))d1=[d1,d11];d2=[d2,d12];endK1=100*d1/n1+50*d2/n1;

%同步指标扩散指数fori=2:m2;D2=X2(i,:)-X2(i-1,:);d21=length(find(D2>0));d22=length(find(D2==0))d3=[d3,d21];d4=[d4,d22];endK2=100*d3/n2+50*d4/n2;%滞后指标扩散指数fori=2:m3;D3=X3(i,:)-X3(i-1,:);d31=length(find(D3>0));d32=length(find(D3==0))d5=[d5,d31];d6=[d6,d32];endK3=100*d5/n3+50*d6/n3;t=2011:2020;plot(t,K1,'-ob',t,K2,'-^k',t,K3,'-+r')xlabel('年份’);ylabel('DI指标')legend('先行指标','同步指标','滞后指标')K1=83.333366.666766.666766.666766.66675010010083.333383.3333K2=50