用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。

做了很久,希望能帮到你。希望大家文明上网,不要随便骂人,呵呵!