财政收入预测问题.doc

财政收入预测问题摘要财政收入是衡量一个地区和国家经济实力的重要标准，控制着国民经济的命脉, 对财政收入进行定量分析并对其做出比较准确的预测可以为相关部门或者企业制定发展规则，实施相关措施提供可靠的理论预测参考。科学、合理地预测财政收入，对于克服年度预算收支确定的随意性和盲目性，正确处理财政与经济的相互关系具有十分重要的意义。本文为预测未来五年的财政收入，运用灰色预测法，建立了相应的预测模型。首先，根据已知信息可知，财政收入与国民收入、工业总产值、农业总产值、总人口、固定资产投资等五种因素有关，结合附表中给出的数据，对财政收入有关的因素进行回归拟合，并结合残差图进行分析，得出财政收入与国民收入、工业总产值、农业总产值、固定资产投资等四种因素的线性相关关系较为显著，表明这四种因素对财政收入的影响较大，而总人口对其影响相对较小一些。其次，对回归分析的得到的方程进行整合，得出财政收入与五种因素的一个关系式：然后，运用灰色预测法，通过软件对数据处理分析，确立了其预测模型方程，即：结合上面的关系式，可以预测出后五年的财政收入分别为：年份（年）财政收入(亿元) 最后对建立的模型进行检验，并对模型作进一步的评价与推广。关键词：软件 财政收入 回归分析 灰色预测法 线性相关1 问题重述国家的财政收入与国民收入、工业总产值、农业总产值、总人口、固定资产投资等因素有关，根据附表中所给出的年的原始数据，对其进行处理分析，构造预测模型，并利用该模型对年至年的财政收入作出预测。2 问题假设与符号说明2.1问题假设1 财政收入只与本文给出的五种因素有关，其它因素不考虑；2 未来五年内无经济危机及重大自然灾害发生；3 本文收集的数据均真实、可靠。2.2 符号说明 ：年份时间 ：回归参数:财政收入与各因素的参数 ：判定系数：预测值（）； ：发展灰数：内生控制灰数 ：数据矩阵：常数矩阵 ：原始数据列：累加生成列 ：相应的模拟误差序列：相应的累减生成序列 ： 序列的均方差：残差的均方差 ：小残差概率：残差平方和 ：各种因素残差概率：绝对误差 ：相对误差序列：关联系数 ：关联度：分辨率 ：平均模拟相对误差：点相当对误差 ：方差比3 问题分析首先，由本文给出的关系可知，财政收入与国民收入、工业总产值、农业总产值、总人口、固定资产投资等五种因素有关。其次，为便于预测，附表中给出了1978-2009年的原始数据，根据这些数据，运用软件分别绘制出财政收入与国民收入、工业总产值、农业总产值、总人口、固定资产投资之间的散点图，从散点图可以看出，它们之间是否存在线性相关关系以及线性关系是否显著。然后对存在线性相关的进行回归分析，对于不存在的需要转化为线性相关后，对其进行回归分析，得出财政收入与各种因素之间的回归方程和判定系数，判定系数越大，回归直线对所测数据的拟合程度就越好，间接反映出了每种因素对财政收入的影响程度。最后，把五种因素作为自变量，财政收入作为因变量，整合为一个方程。运用灰色预测法，建立预测模型并由此预测出后五年的财政收入4 模型的建立与求解4.1 回归分析由题目中的已知信息可知，财政收入与国民收入、工业总产值、农业总产值、总人口、固定资产投资等因素有关，首先画出散点图见下图财政收入与五大因素的散点图由以上的散点图，财政收入与各因素大致上为线性关系。设回归模型为根据题目中给出的数据，利用软件对财政收入和各个因素进行回归分析，得出财政收入和各因素的判定系数和回归方程如下表：影响因素回归系数判定系数回归方程国民收入（亿元）工业总产值（亿元）农业总产值（亿元）总人口（万人）固定资产投资（亿元） 通过以上图中的判定系数可看出财政收入和国民收入、工业总产值、农业总产值、固定资产投资线性关系较强且呈正线性相关，和总人口线性关系中等。说明得出的回归方程较可靠。4.2 建立预测模型由上面的计算可知，财政收入和国民收入、工业总产值、农业总产值、总人口、固定资产投资呈线性关系，由此建立以下模型用编程可求解出观察下面财政收入和各因素的残差图可知除2001，2002，2007，2008，2009这五年数据外其余的残差的置信区间均包含零点，这五年的数据应视为异常年，将其剔除后重新计算，可得4.3 财政收入预测4.3.1 灰色预测法介绍1、灰色预测法是一种对含有不确定因素的系统进行预测的方法。灰色系统是介于白色系统和黑色系统之间的一种系统。2、白色系统是指一个系统的内部特征是完全已知的，即系统的信息是完全充分的。而黑色系统是指一个系统的内部信息对外界来说是一无所知的，只能通过它与外界的联系来加以观测研究。灰色系统内的一部分信息是已知的，另一部分信息时未知的，系统内各因素间具有不确定的关系。灰色预测方法的特点表现在：首先是它把离散数据视为连续变量在其变化过程中所取的离散值，从而可利用微分方程式处理数据；而不直接使用原始数据而是由它产生累加生成数，对生成数列使用微分方程模型。这样，可以抵消大部分随机误差，显示出规律性。灰色预测的类型灰色时间序列预测；即用观察到的反映预测对象特征的时间序列来构造灰色预测模型，预测未来某一时刻的特征量，或达到某一特征量的时间。畸变预测；即通过灰色模型预测异常值出现的时刻，预测异常值什么时候出现在特定时区内。 系统预测；通过对系统行为特征指标建立一组相互关联的灰色预测模型，预测系统中众多变量间的相互协调关系的变化。拓扑预测；将原始数据作曲线，在曲线上按定值寻找该定值发生的所有时点，并以该定值为框架构成时点数列，然后建立模型预测该定值所发生的时点。在本文中建立灰色预测模型之前，需先对原始时间序列进行数据处理，经过数据处理后的时间序列即称为生成列。灰色系统常用的数据处理方式有累加和累减两种。4.3.2 预测求解为了减小误差，分别选用10、20、30个数据进行预测=98001 108068 119096 135174 159587 184089 213132 259259 302853 343464=18718 21826 26937 35260 48109 59811 70143 78061 83024 88479 98001108068 119096 135174 159587 184089 213132 259259 302853 343464=3645 4063 4546 4890 5331 5986 7244 9041 10274 12051 15037 17001 18718 21826 26937 35260 48109 59811 70143 78061 83024 88479 98001 108068 119096 135174 159587 184089 213132 259259 302853 343464第1步 构造累加生成序列=98001 206069 325165 460339 619926 804015 1017147 1276406 1579259 1922723=18718 40544 67481 102741 150850 210661 280804 358865 441889 530368 628369 736437 855533 990707 1150294 1334383 1547515 1806774 2109627 2453091=3645 7708 12254 17144 22475 28461 35705 44746 5502067071 82108 99109 117827 139653 166590 201850 249959309770 379913 457974 540998 629477 727478 835546 954642 1089816 1249403 1433492 1646624 1905883 2208736 2552200第2步 构造数据矩阵和数据向量 第3步 计算=第4步 得出预测模型 第5步 残差检验计算绝对残差和相对残差序列绝对残差序列：=0, 6947, 1258, 2147, 436, 2389, 4176, 6024, 7752, 424相对残差序列：=0, 6.43%, 1.06%, 1.59%, 0.27%, 1.30%, 1.96%, 2.32%, 2.56%, 0.12%绝对残差序列：=2，5528，4403，645，6973，12682，16148，16200，12151，7281，4974，1489，3009，4720，687，465，2758，18237，26719，27102相对残差序列：=0.01%，25.33%，16.35%，1.83%，14.49%，21.20%，23.02%，20.75%，14.64%，8.23%，5.08%，1.38%，2.53%，3.49%，0.43%，0.25%，1.29%，7.03%，8.82%，7.89%绝对残差序列：=0，7447，8681，10310，12136，14085，15820，17463，20182，22947，25180，29213，34388，39200，43188，45323，44491，46598，52135，62451，78442，97065，115212，136941，162449，188357，212190，243129，277794，304876，345408，401468相对残差序列：=0，183.29%，190.96%，210.84%，227.65%，235.30%，218.39%，193.15%，196.44%，190.42%，167.45%，171.83%，183.72%，179.60%，160.33%，128.54%，92.48%，77.91%，74.33%，80.00%，94.48%，109.70%，117.56%，126.72%，136.40%，139.34%，132.96%，132.07%，130.34%，117.60%，114.05%，116.89%经过以上对10，20，32个数据的对比，得出10个数据的残差最小，因此选用10个数据预测出来的模型更可靠。第6步 进行关联度检验第一种相关因素：国民收入(1) 计算序列与的绝对残差序列(k)0，6932，1227，2195，508，2493，4318，5833，7500，752min(k) = min0，6932，1227，2195，508，2493，4318，5833，7500752=0max(k) = max0，6932，1227，2195，508，2493，4318，5833，7500，752=7500(2) 计算关联系数由于只有两个序列（即一个参考序列，一个被比较序列）故不再寻求第二级最小差和最大差。求得1，0.351，0.7535，0.631，0.881，0.601，0.4648，0.3913，0.33，0.833(3) 计算关联度=0.6237r=0.6237是满足P=0.5时的检验准则r0.6的。第7步 后验差检验（1）计算：（2）计算 序列的均方差：=85464.15（3）计算残差的均值：=3175.8（4）计算残差的均方差：=2778.4（5）计算C：=0.0325（6）计算小残差概率：=85464.15*0.6237=533040，3756.2，1948.8，980.8，2667.8，682.8，1142.2，2657.2，4324.2，2423.8所有都小于，故小残差概率=1，而同时C=0.03250.35，故模型=611315.3312536296513314.3312536296合格。第8步、预测：第一种相关因素据以上的方法，可预测出后五年其他各种相关因素工业总产值、农业总产值、总人口、固定资产（具体求解过程见附录D）序号年份国民收入（亿元）工业总产值（亿元）农业总产值（亿元）总人口（亿元）固定资产投资（亿元）1120104011651900733982313433025018812201146753422165244498135080307881132012544884258478497221358353788781420136350303014225555913659346624615201474009035150062081137356573762据以上回归出来的财政收入和各个因素的关系。可预测后五年的财政收入2010年的财政收入 2011年的财政收入 2012年的财政收入 2013年的财政收入 2014年的财政收入 5 模型检验与误差分析通过对整合得到的关系式分析求解，预测出财政收入费用与实际费用相比较，得出相对误差如下表：表5-1 财政收入的相对误差百分比年份国民收入（亿元）工业总产值(亿元)农业总产值（亿元）总人口（万人）固定资产投资（亿元）实际财政收入(亿元)预测财政收入（亿元）相对误差（%）1978364517451026962595501132125711.04197940631914127097542564114611944.19198045462192137298705568116011912.671981489022561560100072961117612697.9119825331238317771016541200121213047.591983598626461978103008136913671344-1.681984724431062316104357183316431487-10.491985904138672564105851254320051798-10.3219861027444932789107507312121221959-8.3219871205152523233109300379221992135-3.00198815037658738651110264754235724233.0819891700172784266112704441026652496-6.7719901871877175062114333451729372444-16.7919912182691025342115823559531492997-4.8319922693711700586711717180803483386811.05199335260164546964118517130724349509817.22199448109224459573119850170425218617918.20199559811286801213612112120019624265605.0919967014333835140151223892291474087355-0.72199778061375431444212362624941865186740.271998830243900414818124761284069876103044.3319998847941034147701257862985511444116752.022000980014555614945126743329181339511578-13.562001108068.4951215781127627372141638615776-3.7220021190965389716537128453435001890418610-1.56200313517462436173821292275556721715225904.03200415958773904214131299887047726396267981.52200518408987365224201307568877431649328153.6920062131321031622404013144810999838760396432.282007259259124799286271321291373245132249326-3.892008302853146183340001328021728286133059548-2.8620093434641569583522613347419413868477729116.84由上表可以看出，在1984、1985、1990、1992、1993、1994、2000这几年财政收入的预测误差较大，但整体误差相对较小。在建立模型过程中，为了便于计算，对数据进行了一些近似处理，必然使数据存在一定的误差，并且为了建模方便，题目中作了很多假设，所以有一些误差存在，但也在误差允许范围内。因此，我们预测出后五年的财政收入可信度较高。6 模型的评价与推广模型的评价优点 1 建立的模型联系实际，解决了各个因素之间的关系，算法简便，易于推广。 2 运用软件求解，得出的数据更加精确，具有科学性。 3 对于其它预测类的问题具有较高的参考价值。缺点 1 模型的建立是在假定某些条件下完成的，但相差不大 2 只对本文给出的五种因素进行考虑，因此建立的模型与实际存在一定的差距，但影响不大。模型的推广针对本文所建立的模型，不仅适用于财政收入问题方面的预测，也适用于粮食产量、人口增长、天气预报等方面的预测，作为参考使用的价值较高。7 参考文献1 钱小军主编数量方法高等教育出版设1999.82 吴建国等编,数学建摸案例精编中国水利水电出版社，北京，2005.3 孙祥 徐流美 吴清编著 MATLAB7.0基础教程 清华大学出版社 2005.54 姜启源 谢金星 叶俊编著数学模型（第三版）高等教育出版社 2003.88 附录附表 表8-1 1978-2009年的原始数据 年份国民收入（亿元）工业总产值(亿元)农业总产值（亿元）总人口（万人）固定资产投资（亿元）财政收入(亿元)1978364517451026962595501132197940631914127097542564114619804546219213729870556811601981489022561560100072961117619825331238317771016541200121219835986264619781030081369136719847244310623161043571833164319859041386725641058512543200519861027444932789107507312121221987120515252323310930037922199198815037658738651110264754235719891700172784266112704441026651990187187717506211433345172937199121826910253421158235595314919922693711700586711717180803483199335260164546964118517130724349199448109224459573119850170425218199559811286801213612112120019624219967014333835140151223892291474081997780613754314442123626249418651199883024390041481812476128406987619998847941034147701257862985511444200098001455561494512674332918133952001108068.4951215781127627372141638620021190965389716537128453435001890420031351746243617382129227555672171520041595877390421413129988704772639620051840898736522420130756887743164920062131321031622404013144810999838760200725925912479928627132129137324513222008302853146183340001328021728286133020093434641569583522613347419413868477附录function yy=model(beta0,X)a=beta0(1);b=beta0(2);c=beta0(3);d=beta0(4);e=beta0(5);%相关系数x1=X(:,1);x2=X(:,2);x3=X(:,3);x4=X(:,4);x5=X(:,5);yy=a*x1+b*x2+c*x3+d*x4+e*x5;%财政收入与五个因素的关系式附录Cfunction f=huise1a=3645 1745 1026 96259 550 11324063 1914 1270 97542 564 11464546 2192 1372 98705 568 11604890 2256 1560 100072 961 11765331 2383 1777 101654 1200 12125986 2646 1978 103008 1369 13677244 3106 2316 104357 1833 16439041 3867 2564 105851 2543 200510274 4493 2789 107507 3121 212212051 5252 3233 109300 3792 219915037 6587 3865 111026 4754 235717001 7278 4266 112704 4410 266518718 7717 5062 114333 4517 293721826 9102 5342 115823 5595 314926937 11700 5867 117171 8080 348335260 16454 6964 118517 13072 434948109 22445 9573 119850 17042 521859811 28680 12136 121121 20019 624270143 33835 14015 122389 22914 740878061 37543 14442 123626 24941 865183024 39004 14818 124761 28406 987688479 41034 14770 125786 29855 1144498001 45556 14945 126743 32918 49512 15781 127627 37214 16386119096 53897 16537 128453 4350062436 17382 129227 55567 21715159587 73904 21413 129988 70477 26396184089 87365 22420 130756 88774 31649213132 103162 24040 131448 109998 38760259259 124799 28627 132129 137324 51322302853 146183 34000 132802 172828 61330343464 156958 35226 133474 194138 68477;R=;for i=1:5 subplot(3,2,i) y=a(:,6);x=a(:,i); plot(x,y,*); %画出散点图 hold on b,brint,r,rint,stats=regress(y,ones(32,1),x);%对数据进行线性回归 x1=x(1):x(end); s=b(1)+b(2)*x1; plot(x1,s,k) %画出进行回归后的直线 R(i)=sqrt(stats(1);%求相关系数rendR附录Dfunction lpx=;s=0;yp=;t=14945 15781 16537 17382 21413 22420 24040 28627 34000 35226;n=length(t);for i=1:n s=s+t(i);%对原始值进行累加 x=x,s;endfprintf( %d,x);z(1)=x(1);for i=2:n z(i)=0.5*(x(i)+x(i-1);%求紧邻均值endformat long gz;z(1)=;B=-z,ones(n-1,1);Y=15781 16537 17382 21413 22420 24040 28627 34000 35226;format long ga=inv(B*B)*B*Y;%求发展系数和灰色作用量aa=input(shuru:)b=input(shuru:)for k=1:(n+5)y=(t(1)-b/a)*exp(-a*(k-1)+b/a;yp=yp,y;endyp;%通过预测模型预测z(1)=yp(1);for i=2:(n+5)z(i)=yp(i)-yp(i-1);%还原值endzfor i=1:nw(i)=abs(t(i)-z(i)/t(i);%误差百分比endformat long gw附录Efunction f=ping2()x=1978 3645 1745 1026 96259 550 11321979 4063 1914 1270 97542 564 11461980 4546 2192 1372 98705 568 11601981 4890 2256 1560 100072 961 11761982 5331 2383 1777 101654 1200 12121983 5986 2646 1978 103008 1369 13671984 7244 3106 2316 104357 1833 16431985 9041 3867 2564 105851 2543 20051986 10274 4493 2789 107507 3121 21221987 12051 5252 3233 109300 3792 21991988 15037 658