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];
  • 2
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

mengke25

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值