用matlab建立经济模型,进行数值模拟。
财政收入预测:首先要知道影响财政收入的主要因素有哪些。我们可以把它们分为:国民收入、工业总产值、总人口、就业人口、固定资产投资等。我们只需要列出主要因素。我们对它们做如下回归分析:(从左到右)
年国民收入(亿元)、工业总产值(亿元)、农业总产值(亿元)
总人口(万人)就业人口(万人)固定资产投资(亿元)财政收入(亿元)
1952 598 349 461 57482 20729 44 184
1953 586 455 175 58796 21364 89 216
1954 707 520 491 60266 21832 97 248
1955 737 558 529 61465 22328 98 254
1956 825 715 556 62828 23018 150 268
837 798 575 64653 23711 139 286
1028 1235 598 65994 26600 256 357
1114 1681 509 67207 26173 338 444
1079 1870 444 66207 25880 380 506
757 1156 434 65859 25590 138 271
677 964 461 67295 25110 66 230
779 1046 514 69172 26640 85 266
943 1250 584 70499 27736 129 323
1152 1581 632 72538 28670 175 393
1322 1911 687 74542 29805 212 466
1249 1647 697 76368 30814 156 352
1187 1565 680 78534 31915 127 303
1372 2101 688 80671 33225 207 447
1638 1747 676 82992 34432 312 564
1780 3156 790 85229 35620 355 638
1833 3365 789 87177 35854 354 658
1978 3684 855 89211 36652 374 691
1993 3696 891 90859 37369 393 655
2121 4254 932 92421 38168 462 692
2052 4309 955 93717 38834 443 657
2189 4925 971 94974 39377 454 732
2475 5590 1058 96259 39856 550 922
2702 6065 1150 97542 40581 564 890
2791 6592 1194 98705 41896 568 826
2927 6862 1273 100072 43280 496 810
自己输入年份。朋友记住:我不输入年份到1956。要想了解数据,1956后面只有七组,前面还有八组。你不能在这里插入表格。抱歉,退后一格。现在按顺序输入即可。我已经瘦了很多,眼睛都花了.....
假设国民收入、工业总产值、农业总产值、总人口、就业人口、固定资产投资分别为x1,x2,x3,x4,x5,x6,财政收入为y,假设它们之间的关系为:
y = a * x 1+b * x2+c * x3+d * x4+e * X5+f * X6
用非线性回归方法求解:
(1)为回归模型建立m文件model.m,如下所示:
函数yy =模型(beta0,x)
a =β0(1);
b =β0(2);
c =β0(3);
d =β0(4);
e =β0(5);
f =β0(6);
x1=x(:,1);
x2=x(:,1);
x3=x(:,1);
x4=x(:,1);
x5=x(:,1);
x6=x(:,1);
YY = a * x 1+b * x2+c * x3+d * x4+e * X5+f * X6;
主要程序如下:shuju.m
x =[598 349 461 57482 20729 44 586 455 175 58796 21364 89 707 520 491 60266 21832 97
737 558 529 61465 22328 98
825 715 556 62828 23018 150
837 798 575 64653 23711 139
1028 1235 598 65994 26600 256
1114 1681 509 67207 26173 338
1079 1870 444 66207 25880 380
757 1156 434 65859 25590 138
677 964 461 67295 25110 66
779 1046 514 69172 26640 85
943 1250 584 70499 27736 129
1152 1581 632 72538 28670 175
1322 1911 687 74542 29805 212
1249 1647 697 76368 30814 156
1187 1565 680 78534 31915 127
1372 2101 688 80671 33225 207
1638 1747 676 82992 34432 312
1780 3156 790 85229 35620 355
1833 3365 789 87177 35854 354
1978 3684 855 89211 36652 374
1993 3696 891 90859 37369 393
2121 4254 932 92421 38168 462
2052 4309 955 93717 38834 443
2189 4925 971 94974 39377 454
2475 5590 1058 96259 39856 550
2702 6065 1150 97542 40581 564
2791 6592 1194 98705 41896 568
2927 6862 1273 100072 43280 496];
y =[184 216 248 254 268 286 357 444 506 271 230 266 323 393 466 352 303 447 564 638 658 691 655 692 657 732 922 890 826 810]';
beta 0 =[0.5-0.03-0.60 0.01-0.02 0.35];
betafit=nlinfit(x,y,'模型',beta0)
结果是:betafit =
0.5243
-0.0294
-0.6304
0.0112
-0.0230
0.3658
即:y = 0.5243 x 1-0.0294 x2-0.6304 x3+0.0112x 4-0.0230 X5+0.3658 X6。
做了很久,希望能帮到你。希望大家文明上网,不要随便骂人,呵呵!