本程序参考EI论文《电力系统碳排放流的计算方法初探》，文中根据潮流来推算电力系统碳排放流方向，根据电流和电势类比出碳流和碳势，程序较为简单和基础，具有较强的可扩展性和适用性，注释清晰，干货满满，下面对文章和程序作简要介绍。
问题
电力系统碳排放流分析的理念为低碳电力带来了新的研究方向。根据潮流计算的结果准确而系统地求解电力系统中碳排放流的分布成为亟待解决的问题。
创新点（解决办法）
在电力系统碳排放流分析理论的基础上，进一步分析了碳排放流和电力系统潮流计算之间的异同，以及电力系统碳排放流的影响因素、计算体系和计算思路；根据碳排放流求解的需要，定义了一些关键矩阵和向量；结合电力系统潮流计算方法与高等电力网络分析方法，在忽略网损的情况下，建立了电力系统碳排放流的基本计算方法，并通过算例系统验证了该方法的正确性。
文中结果
程序结果：复现结果稍有微小偏差！
部分程序
%% 求系统潮流power_flow;%% 论文复现K=5;%发电机数N=14;%系统节点数M=11;%负荷数%% 支路潮流分布矩阵PBPB=Pij*baseMVA;for k=1:length(branch(:,1))branch(k,14)=PB(branch(k,1),branch(k,2));branch(k,15)=PB(branch(k,2),branch(k,1));%% 机组注入分布矩阵PG=zeros(K,N);for k=1:length(gen(:,1))PG(k,gen(k,1))=gen(k,2);%% 辅助变量PZxigama=ones(1,K+N);PZ=[PB;PG];%% 节点有功通量矩阵PN=diag(xigama*PZ);disp('节点有功通量：')disp((xigama*PZ)')%% 发电机组碳排放强度向量EG=[875;525;0;520;0];%% 节点碳势分布向量EN=(PN-PB')^(-1)*PG'*EG;disp('节点碳势分布向量：')disp(EN)%% 负荷分布矩阵PL=zeros(M,N);bus_copy=bus;if bus_copy(kk,3)>0PL(k,kk)=bus_copy(kk,3);bus_copy(kk,3)=0;
部分内容源自网络，侵权联系删除！
欢迎感兴趣的小伙伴点击文末阅读原文获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！
更多创新性预测及优化程序请点击公众号首页：《预测优化》|《综合能源》链接！
