交通经济学代码补充

 交通经济学作业——投入产出表分析，代码实现直接消耗矩阵，完全消耗矩阵，逆阵系数，感应度及影响力系数，中间投入率，中间需求率以及产业部门间联系深度。

%直接消耗系数矩阵的计算方法
%方法一：X为总产出，Q为中间投入矩阵，a为直接消耗系数矩阵
load('X.mat')
load('Q.mat')
n=149;
a=Q*inv(diag(X));
%方法二：X为总产出，Q为中间投入矩阵
% Yr=repmat(X',n,1);
% A=Q./Yr;

%完全消耗系数矩阵计算方法
%b=inv(I-a)-I
b=inv(eye(n)-a)-eye(n);

%间接消耗系数矩阵的计算方法;
%有两种计算方法，1.c=b*a;2.b=c+a
c=b*a;
% c=b-a;

%计算中间需求率和中间投入率
Qh=sum(Q,2);
Qz=sum(Q);
Xh=X;
Xz=X';
ratio_xuqiu=Qh./(X);
ratio_touru=Qz./(X');

%产业部门间联系广度与深度
r=Q./Qz;

%直接折旧系数
% aDj=D./X';

%逆阵系数
%使用inv(I-a)来计算,a为我们的直接消耗系数矩阵
a_inv=inv(eye(n)-a);

%计算感应度系数及影响力系数
W=a_inv;
%计算分子部分
meanh=mean(W,2);
meanz=mean(W);
%计算分母部分
meanhall=mean(meanh(:));
meanzall=mean(meanz(:));
% meanall=mean(W(:));
%感应度系数
S=meanh/meanhall;
%影响力系数
T=meanz/meanzall;