Matlab计算完全消耗系数

数据来源wiod

1.计算

clear,clc
eta=cell(14,2); 
year=2000:2013
for i=1:14
    eta{i,2}=1999+i
end
[Type,Sheet,Fromat]=xlsfinfo('F:\yida_project\IND_NIOT_nov16.xlsx');
for i=1:length(Sheet)
    [Num,Text,RawData] =xlsread('F:\yida_project\IND_NIOT_nov16.xlsx',Sheet{i});
end
sheet=cell(14,1)
for i=1:14
sheet{i,1}=Sheet{1,i+3}
end
for e=1:14
P=xlsread('F:\yida_project\IND_NIOT_nov16.xlsx',sheet{e,1},'B2:BE57');
Q=xlsread('F:\yida_project\IND_NIOT_nov16.xlsx',sheet{e,1},'B58:BE113');
GO=xlsread('F:\yida_project\IND_NIOT_nov16.xlsx',sheet{e,1},'B121:BE121');
O=(P+Q)./GO
s=size(O);
for i=1:s
    for j =1:s
     K(i,j)=isnan(O(i,j));
     end 
end
  for i=1:s
    for j=1:s
    if K(i,j)==1
        O(i,j)=0;
    end      
   end  
end
I=eye(s);
eta{e,1}=(I-O)^(-1) -I
end
for i=1:14
filetitle=['F:\yida_project\' 'data' num2str(i) '.xlsx'];
dataTmp=eta{i,1};
xlswrite(filetitle,dataTmp);
end

2.合并

clear,clc

QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data1.xlsx','Sheet1','E24:V56');
 wq1=[sum(QWE(:,:))];
 
  QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data2.xlsx','Sheet1','E24:V56');
 wq2=[sum(QWE(:,:))];
 
  QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data3.xlsx','Sheet1','E24:V56');
 wq3=[sum(QWE(:,:))];
 
  QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data4.xlsx','Sheet1','E24:V56');
 wq4=[sum(QWE(:,:))];
 
  QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data5.xlsx','Sheet1','E24:V56');
 wq5=[sum(QWE(:,:))];
 
  QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data6.xlsx','Sheet1','E24:V56');
 wq6=[sum(QWE(:,:))];
 
   QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data7.xlsx','Sheet1','E24:V56');
 wq7=[sum(QWE(:,:))];
 
 QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data8.xlsx','Sheet1','E24:V56');
 wq8=[sum(QWE(:,:))];
 
 QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data9.xlsx','Sheet1','E24:V56');
 wq9=[sum(QWE(:,:))];
 
 QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data10.xlsx','Sheet1','E24:V56');
 wq10=[sum(QWE(:,:))];
 
 QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data11.xlsx','Sheet1','E24:V56');
 wq11=[sum(QWE(:,:))];
 
  QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data12.xlsx','Sheet1','E24:V56');
 wq12=[sum(QWE(:,:))];
 
  QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data13.xlsx','Sheet1','E24:V56');
 wq13=[sum(QWE(:,:))];
 
 QWE=xlsread('C:\Users\Martin Song\Desktop\yida_project\data14.xlsx','Sheet1','E24:V56');
 wq14=[sum(QWE(:,:))];
 
wq=[wq1;wq2;wq3;wq4;wq5;wq6;wq7;wq8;wq9;wq10;wq11;wq12;wq13;wq14];