预测某市未几年工业总产值

预测某市工业总产值一、摘要工业生产总值是一个城市经济发展的重要衡量标准，对城市的评估也离不开观测工业总产值，以往年工业总产值为基础，对未来几年工业总产值进行预测。模型：Logistic阻滞增长模型，假设城市对工业总产值存在最大值，建立微分方程，先利用Matlab编程求解出各个参数，再对未来的产值进行预测。得到接下来几年的工业产值（单位：元）30622.46534167.35937899.6341782.28645771.265然后比较拟合数值与实际数值，得到拟合相对误差和绝对误差。模型仅使用了几年的工业总产值，实际一个城市的工业总产值受多方面影响，随时间后移，科技发展，工业总产值的最大值应有升高。二、关键词工业总产值预测，Logistic模型，最小二乘拟合三、问题描述下表是某市历年工业生产总值，是建立数学模型，分析并预测未来几年工业总产值的变化某市工业总产值表年份工业生产总值年份工业生产总值年份工业生产总值1980598.7519901642.7520007022.981981620.1219911947.1820017806.181982634.6519922429.9620028730.001983663.5319933327.04200311708.491984728.1219944255.19200414595.291985862.7319955349.53200516876.781986952.2119965126.22200619631.2319871073.8419975649.93200723108.6319881304.6619985763.67200825968.3819891524.6719996213.24200924888.08四、模型假设（1）假设表中所给数据能反映该市工业总产值的基本情况；（2）假设该市经济状况平稳，没有出现政府作出重大决策，科技发生重大突破，经济危机的影响；（3）假设该市的工业生产模式没有发生改变；（4）自然资源和环境所容许的最大工业总产值为常数mN，并工业总产值的净相对增长率是工业总产值的线性递减函数，设为()(1)mNrNrN，表示工业总产值相对增产率随()Nt的增加而减少，其中r为固有增长率。当()mNtN时，工业总产值净增长率()rN趋于零。五、模型建立与求解1，模型建立求解使用阻滞增长模型：00()()1()()mdNtNtrNtdtNNtN方程右端因子()Nt体现工业总产值的增长趋势，因子()1mNtN则体现了自然资源和环境对工业总产值的阻滞作用。显然()Nt越大，前一因子越大，后一因子越小，表明工业总产值的增长是两个因子共同作用的结果。用分离变量法可解得方程的解为0()0()1(1)mrttmNNtNeN2，模型的参数估计与校验用Logisitc模型做预测，需要估计方程解中的三个参数,mNr和0N，其中0N可以计算初始时刻的工业总产值，也可以用数据拟合出来。采用1980~2009年工业总产值作为拟合数据，用MATLAB中最小二乘曲线拟合命令lsqcurvefit来估计参数。然后用拟合的参数带入方程的解，用你来预测2010~2020年的工业总产值。利用程序算得参数0100369,r=0.1617,473.98mNN,故第t年工业总产值为：0.1617(1980)1003691210.758ttNe3，预测结果年份2010201120122013201420152016201720182019预测产值30622341673789941782457714981753868578716177665538拟合绝对误差：拟合相对误差：附MATLAB程序：%logistic模型模拟，预测2010~2020年工业总产值clc,cleartdata=1980:2009;ndata=[598.75620.12634.65663.53728.12862.73952.211073.841304.661524.67...1642.751957.182429.963327.044255.195349.535126.225649.935763.676213.24...7022.987806.188730.0011708.4914595.2916876.7819631.2323108.6325968.3824888.08];a0=[260000.1598.75];a=lsqcurvefit(@logistic,a0,tdata,ndata)tfit=1980:2019;nfit=a(1)./(1+(a(1)/a(3)-1)*exp(-a(2)*(tfit-1982)));figure(1);plot(tdata,ndata,'k*',tfit,nfit,'-k');legend('实际产值','预测产值','location','NorthWest');xlabel('年份');ylabel('产值');title('logistic模型对工业总产值拟合预测');es=abs(ndata-nfit(1:30));%绝对误差esc=abs((ndata-nfit(1:30))./ndata);%相对误差figure(2);六、模型分析与评价（1）观察拟合绝对误差在，在1992年以前差别极小，但是考虑到产值的总大小在增大，观测相对误差图，发现2003年以后的相对误差均小于10%，故可以认为该模型预测的值在短期内是可信的；（2）有拟合预测曲线可以看出，图形比较接近一个指数函数，可以认为工业总产值距离饱和值还比较遥远，由拟合的数据，在2019年才达到2mN即增长的拐点，一直处于较快增长的情况；（3）Logistic模型仅考虑了工业产值饱和值对增长速率的影响，相比于实际情况还是过于简洁，尤其是饱和值应当随着时间的增长而变大，是本模型的改进方向；（4）Logistic模型虽然本身是为了预测人口而建立的，但在本问题中对工业总产值的预测也是较为可信的，它可以运用在实际多领域上，除了对种群增长的使用，也可以适用于冰块融化，传染病蔓延，经济增长等。参考书目：数学建模教程，梅正阳、韩志斌，北京：科学出版社，2012.8