|
对WIOD去除行列名称,假设文件为WIOT_2000.csv等。
WIOT中,2464*2464为中间使用额,第2465列开始,每5列一循环为各国实际最终需求,最后一列为总产出。
求直接消耗系数与完全消耗系数,MATLAB 代码如下,仅供参考:
clc,clear;
dataframe = cell(15,1);
for y = 0:14
tempfile = sprintf('WIOT_%d.csv', y+2000);
dataframe{y+1} = readtable(tempfile);
end
for y = 0:14
X = dataframe{y+1};
X = table2array(X);
O = X(1:2464,end);
%direct consumption coefficient matrix.
A = X(1:2464,1:2464)./O;
%complete consumption coefficient matrix.
I = eye(2464,2464);
B = pinv(I-A) - I;
tempresults1 = sprintf('DCC%d.csv', y+2000);
tempresults2 = sprintf('CCC%d.csv', y+2000);
writematrix(A, tempresults1);
writematrix(B, tempresults2);
end
|