利用MATLAB进行基础的投入产出表分析实验报告

利用MATLAB求直接消耗系数 完全消耗系数等基本投入产出分析

学 生 实 验 报 告

学 院： XX学院 课程名称：宏观经济数量分析方法与模型

专业班级：

姓 名：

学 号：

利用MATLAB求直接消耗系数 完全消耗系数等基本投入产出分析

学生实验报告

一、实验综述

1、实验目的及要求

学会投入产出表的数据分析

掌握基础的MATLAB矩阵运算

2、实验仪器、设备或软件

MATLAB7.0

二、实验过程（实验步骤、记录、数据、分析）

1.2007年投入产出表

利用MATLAB求直接消耗系数 完全消耗系数等基本投入产出分析

2.直接消耗系数矩阵a

>> x=[6877.2 24657.5 259.3 379.7 1784.7 385.6

10248.3 326099.9 37925.5 10912.9 8849.6 27220.2

11.3 159.8 598 122.9 127.6 980.3

797.5 11705.8 4737.4 2264.7 2817.8 2798.1

851.7 12336.3 2022.7 789.6 1407.4 5811.2

1447.9 19918 2665.3 2978.2 5761.8 14103.6] %输入x矩阵

x =

1.0e+005 *

0.0688 0.2466 0.0026 0.0038 0.0178 0.0039

0.1025 3.2610 0.3793 0.1091 0.0885 0.2722

0.0001 0.0016 0.0060 0.0012 0.0013 0.0098

0.0080 0.1171 0.0474 0.0226 0.0282 0.0280

0.0085 0.1234 0.0202 0.0079 0.0141 0.0581

0.0145 0.1992 0.0267 0.0298 0.0576 0.1410

>> X=[48893 514859.1 62721.7 32430.9 43648 116306.3] %输入X矩阵

X =

1.0e+005 *

0.4889 5.1486 0.6272 0.3243 0.4365 1.1631

>> for i=1:6,a(:,i)=x(:,i)/X(i);end %计算直接消耗系数矩阵

>> a

a =

0.1407 0.0479 0.0041 0.0117 0.0409 0.0033

0.2096 0.6334 0.6047 0.3365 0.2027 0.2340

0.0002 0.0003 0.0095 0.0038 0.0029 0.0084

0.0163 0.0227 0.0755 0.0698 0.0646 0.0241

0.0174 0.0240 0.0322 0.0243 0.0322 0.0500

0.0296 0.0387 0.0425 0.0918 0.1320 0.1213

利用MATLAB求直接消耗系数 完全消耗系数等基本投入产出分析

3.完全需求系数矩阵（列昂惕夫逆矩阵）b1

>> I=eye(6); %定义数据全为1的6阶矩阵

>> b1=inv(I-a); %对（I-a）矩阵就逆

>> b1

b1 =

1.2128 0.1772 0.1260 0.0886 0.1030 0.0613

0.8230 3.0700 2.0412 1.2440 0.8933 0.9252

0.0016 0.0030 1.0126 0.0065 0.0057 0.0110

0.0472 0.0889 0.1457 1.1172 0.1040 0.0617

0.0481 0.0900 0.0992 0.0703 1.0712 0.0879

0.0893 0.1641 0.1732 0.1854 0.2149 1.2010

4.完全消耗系数矩阵b=(I-a)^-1-I

>> b=b1-I

b =

0.2128 0.1772 0.1260 0.0886 0.1030 0.0613

0.8230 2.0700 2.0412 1.2440 0.8933 0.9252

0.0016 0.0030 0.0126 0.0065 0.0057 0.0110

0.0472 0.0889 0.1457 0.1172 0.1040 0.0617

0.0481 0.0900 0.0992 0.0703 0.0712 0.0879

0.0893 0.1641 0.1732 0.1854 0.2149 0.2010

5.分配系数矩阵h

>> for j=1:6,h(j,:)=x(j,:)/X(j);end

>> h

h =

0.1407 0.5043 0.0053 0.0078 0.0365 0.0079

0.0199 0.6334 0.0737 0.0212 0.0172 0.0529

0.0002 0.0025 0.0095 0.0020 0.0020 0.0156

0.0246 0.3609 0.1461 0.0698 0.0869 0.0863

0.0195 0.2826 0.0463 0.0181 0.0322 0.1331

0.0124 0.1713 0.0229 0.0256 0.0495 0.1213

6.生产依存度（L）和生产依存度系数（l）

>> L=sum(a) %对直接消耗系数矩阵求列和

L =

0.4138 0.7670 0.7686 0.5380 0.4754 0.4411

>> for i=1:6,l(i)=6*L(i)/sum(L);end

利用MATLAB求直接消耗系数 完全消耗系数等基本投入产出分析

>> l

l =

0.7295 1.3519 1.3548 0.9483 0.8379 0.7775

7.生产拉动力（U）和生产拉动力系数（u）

>> U=sum(h)

U =

0.2173 1.9551 0.3038 0.1444 0.2244 0.4171

>> for i=1:6,u(i)=6*U(i)/sum(U);end

>> u

u =

0.3997 3.5960 0.5588 0.2657 0.4127 0.7671

8.生产支撑度（F）和生产支撑度系数（f）

>> for i=1:6,F(i)=sum(h(i,:));end 对h矩阵求行和

>> F

F =

0.7024 0.8182 0.0319 0.7746 0.5320 0.4030

>> for i=1:6,f(i)=6*F(i)/sum(F);end

>> f

f =

1.2920 1.5049 0.0586 1.4247 0.9784 0.7413

9.产业影响力（R）和影响力系数（r）

>> R=sum(b1)

R =

2.2220 3.5932 3.5979 2.7120 2.3921 2.3481

>> for i=1:6,r(i)=6*R(i)/sum(R);end

>> r

r =

0.7905 1.2783 1.2800 0.9648 0.8510 0.8354

利用MATLAB求直接消耗系数 完全消耗系数等基本投入产出分析

10.感应度（S）和感应度系数（s）

>> for i=1:6,S(i)=sum(b1(i,:));end

>> S

S =

1.7690 8.9967 1.0404 1.5647 1.4666 2.0278

>> for i=1:6,s(i)=6*S(i)/sum(S);end

>> s

s =

0.6293 3.2007 0.3701 0.5567 0.5218 0.7214

11.各部门产出分别增加10% 15% 15% 10% 10% 10%时总产出变化量sum(DX)

增加后各部门产出xx1 各部门总产出XX1=xx1./L 则各部门总产出变化DX=XX1-X 总产出变化量sum(DX)

>> xx1(1)=x1(1)*(1+0.1);

>> xx1(2)=x1(2)*(1+0.15);

>> xx1(3)=x1(3)*(1+0.15);

>> xx1(4)=x1(4)*(1+0.1);

>> xx1(5)=x1(5)*(1+0.1);

>> xx1(6)=x1(6)*(1+0.1);

>> xx1

xx1 …… 此处隐藏：1926字，全部文档内容请下载后查看。喜欢就下载吧 ……