EI论文复现：电力系统碳排放流的计算方法初探程序代码！

本程序参考EI论文《电力系统碳排放流的计算方法初探》，文中根据潮流来推算电力系统碳排放流方向，根据电流和电势类比出碳流和碳势，程序较为简单和基础，具有较强的可扩展性和适用性，注释清晰，干货满满，下面对文章和程序作简要介绍。

问题：

电力系统碳排放流分析的理念为低碳电力带来了新的研究方向。根据潮流计算的结果准确而系统地求解电力系统中碳排放流的分布成为亟待解决的问题。

创新点（解决办法）：

在电力系统碳排放流分析理论的基础上，进一步分析了碳排放流和电力系统潮流计算之间的异同，以及电力系统碳排放流的影响因素、计算体系和计算思路；根据碳排放流求解的需要，定义了一些关键矩阵和向量；结合电力系统潮流计算方法与高等电力网络分析方法，在忽略网损的情况下，建立了电力系统碳排放流的基本计算方法，并通过算例系统验证了该方法的正确性。

文中结果：

程序结果：复现结果稍有微小偏差，但是无任何影响！

部分程序：

%% 求系统潮流
power_flow;
%% 论文复现
K=5;%发电机数
N=14;%系统节点数
M=11;%负荷数
%% 支路潮流分布矩阵PB
PB=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));
end
for k=1:14
    for kk=1:14
        if PB(k,kk)<0
           PB(k,kk)=0;
        end
    end
end
%% 机组注入分布矩阵
PG=zeros(K,N);
for k=1:length(gen(:,1))
    PG(k,gen(k,1))=gen(k,2);
end
%% 辅助变量PZ
xigama=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;
for k=1:M
    for kk=1:N
        if bus_copy(kk,3)>0
            PL(k,kk)=bus_copy(kk,3);
            bus_copy(kk,3)=0;
            break
        end
    end
end

欢迎感兴趣的小伙伴关注下方公众号获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！