2011年5月13日更新。【已另外写了篇博文,见以下链接:
最近在写一篇关于中国部门间温室气体排放关系的论文,已快完成。
原文是用Excel计算的,自动化程度不高,最关键的是Excel矩阵计算功能薄弱。所以下一步打算将这个数据库和计算过程完全移植到MATLAB上,一来方便自己以后使用,二来也方便和同行交流(不用一遍又一遍地解释怎么用Excel计算)。
===MATLAB Code,非常简单的代码===
%% China EIO-LCA 2007 Model for Greenhouse Gas Emissions
% Written by JI Junping(or Jacky JI, my English name)
% Last updated inShenzhen, China on May 13th, 2011
%% Load all necessaryvariables
% R is a matrix withdiagonal elements representing the greenhouse gas
% emissions per RMB yuanof output for 2007.
% I is the identitymatrix.
% A is the input-outputdirect requirements matrix.
% y is the vector ofdesired output.
% For more information onEIO-LCA, please visit: http://www.eiolca.net/
% For more information onmy work on China EIO-LCA model, please visit:
%http://goo.gl/lHaPKorhttp://goo.gl/VyxfS
% If you have anyquestion or comment, contact me by Email: jackyjjp(at)gmail.com
loadEIOLCA2007.mat
%% Calculate sectoralproduction-based greenhouse gas emissions induced by sectoral final demand
% B is greenhouse gasemissions matrix.
% B_column_sum is thevector of greenhouse gas emissions induced by
% sectoral final demand.
% B_row_sum is the vectorof greenhouse gas emission produced in the
% production process ofeach sector.
B = R*inv(I-A)*inv(y)
B_column_sum = sum(B)
B_row_sum = sum(B,2)
转载本文请联系原作者获取授权,同时请注明本文来自计军平科学网博客。
链接地址:http://blog.sciencenet.cn/blog-241374-347142.html
上一篇:需注意的统计问题
下一篇:论文:基于投入产出分析的中国部门温室气体排放研究