预测某市未来几工业总产值.docVIP

预测某市工业总产值 摘要 工业生产总值是一个城市经济发展的重要衡量标准，对城市的评估也离不开观测工业总产值，以往年工业总产值为基础，对未来几年工业总产值进行预测。 模型：Logistic阻滞增长模型，假设城市对工业总产值存在最大值，建立微分方程，先利用Matlab编程求解出各个参数，再对未来的产值进行预测。 得到接下来几年的工业产值（单位：元） 30622.465 34167.359 37899.63 41782.286 45771.265 然后比较拟合数值与实际数值，得到拟合相对误差和绝对误差。 模型仅使用了几年的工业总产值，实际一个城市的工业总产值受多方面影响，随时间后移，科技发展，工业总产值的最大值应有升高。 关键词 工业总产值预测，Logistic模型，最小二乘拟合 问题描述 下表是某市历年工业生产总值，是建立数学模型，分析并预测未来几年工业总产值的变化 某市工业总产值表 年份 工业生产总值 年份 工业生产总值 年份 工业生产总值 1980 598.75 1990 1642.75 2000 7022.98 1981 620.12 1991 1947.18 2001 7806.18 1982 634.65 1992 2429.96 2002 8730.00 1983 663.53 1993 3327.04 2003 11708.49 1984 728.12 1994 4255.19 2004 14595.29 1985 862.73 1995 5349.53 2005 16876.78 1986 952.21 1996 5126.22 2006 19631.23 1987 1073.84 1997 5649.93 2007 23108.63 1988 1304.66 1998 5763.67 2008 25968.38 1989 1524.67 1999 6213.24 2009 24888.08 模型假设 假设表中所给数据能反映该市工业总产值的基本情况； 假设该市经济状况平稳，没有出现政府作出重大决策，科技发生重大突破，经济危机的影响； 假设该市的工业生产模式没有发生改变； 自然资源和环境所容许的最大工业总产值为常数，并工业总产值的净相对增长率是工业总产值的线性递减函数，设为，表示工业总产值相对增产率随的增加而减少，其中r为固有增长率。当时，工业总产值净增长率趋于零。 模型建立与求解 模型建立求解 使用阻滞增长模型： 方程右端因子体现工业总产值的增长趋势，因子则体现了自然资源和环境对工业总产值的阻滞作用。显然越大，前一因子越大，后一因子越小，表明工业总产值的增长是两个因子共同作用的结果。 用分离变量法可解得方程的解为 模型的参数估计与校验 用Logisitc模型做预测，需要估计方程解中的三个参数和，其中可以计算初始时刻的工业总产值，也可以用数据拟合出来。采用1980~2009年工业总产值作为拟合数据，用MATLAB中最小二乘曲线拟合命令lsqcurvefit来估计参数。然后用拟合的参数带入方程的解，用你来预测2010~2020年的工业总产值。利用程序算得参数, 故第t年工业总产值为： 预测结果 年份 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 预测产值 30622 34167 37899 41782 45771 49817 53868 57871 61776 65538 拟合绝对误差： 拟合相对误差： 附MATLAB程序： %logistic模型模拟，预测2010~2020年工业总产值clc,clear %logistic模型模拟，预测2010~2020年工业总产值 clc,clear tdata=1980:2009; ndata=[598.75 620.12 634.65 663.53 728.12 862.73 952.21 1073.84 1304.66 1524.67... 1642.75 1957.18 2429.96 3327.04 4255.19 5349.53 5126.22 5649.93 5763.67 6213.24... 7022.98 7806.18 8730.00 11708.49 14595.29 16876.78 19631.23 23108.63 25968.38 24888.08 ]; a0=[26000 0.1 598.75]; a=lsqcurvefit(@logistic,a0,tdata,ndata) tfit=1980:2019; nfit=a(1)./(1+(a(1)