本程序参考中文核心论文《基于电力系统碳排放流理论的碳排放分摊模型研究》，文中运用潮流追踪对系统网损引起的碳排放进行分配定量，最后得到系统实时碳排放的产权分配结果，程序较为简单和基础，具有较强的可扩展性和适用性，注释清晰，干货满满，下面对文章和程序作简要介绍。
创新点
提出一种电力系统碳排放产权界定分配的模型。分摊模型以电力系统碳排放结构分解为切入点，借鉴负荷分析理论将电力系统碳排放以其所对应的电能用途分成了综合用电负荷碳排放、网损碳排放、厂用电碳排放等３部分。
确定结构成分并与电力系统碳排放流的计算指标对应后，计算得到系统中各碳排放流成分分布；运用潮流追踪对系统网损引起的碳排放进行 分配定量，最后得到系统实时碳排放的产权分配结果。在对IEEE14标准节点进行算例分析后，其结果证明了本文所述模型的可行性与正确性。
文中结果
程序结果
部分程序
%% 论文复现K=5;%发电机数N=14;%系统节点数M=11;%负荷数%% 式(2)的计算PE=zeros(N);%支路潮流分布矩阵for k=1:length(branch(:,1))if branch(k,14)>0PE(branch(k,1),branch(k,2))=branch(k,14);PE(branch(k,2),branch(k,1))=-branch(k,14);PG=zeros(K,N);%机组注入分布矩阵for k=1:length(gen(:,1))PG(k,gen(k,1))=gen(k,2);xigama=ones(1,K+N);PZ=[PE;PG];PN=diag(xigama*PZ);%节点有功通量矩阵EG=[875;525;0;520;0];%碳排放强度EN=(PN-PE')^(-1)*PG'*EG;%节点碳势分布向量%% 式(3)的计算PL=zeros(M,N);%负荷分布矩阵bus_copy=bus;RL=PL*EN;%负荷碳流率向量RL=RL/1000;%转换单位%% 式(4)的计算Pl=zeros(N);%支路有功损耗分布矩阵for k=1:length(branch(:,1))if sum(branch(k,[14 16]))>0Pl(branch(k,1),branch(k,2))=sum(branch(k,[14 16]));Pl(branch(k,2),branch(k,1))=-sum(branch(k,[14 16]));Rl=diag(EN)*Pl;%支路功率损耗引起的碳排放分布Rl=Rl/1000;%转换单位%% 式(5)的计算Pa=[124.5;40;60;19;20];%厂用电功率？？？？Ra=diag(EG)*Pa;%厂用电碳流率向量Ra=Ra/1000;%转换单位%% 式(6-10)的计算Ad=eye(N)-PE*PN^(-1);%上游分布矩阵PG_l=zeros(K,1);%机组对系统有功网损的贡献RG_l=zeros(K,1);%各个机组应承担的碳排放J=[1 1;2 2;3 3;4 6;5 8];%各个发电机与其连接的节点entaN=zeros(K,N);%第j个元素为1，其余分量为0的N维向量entaN(k,J(k,2))=1;PG_l(k)=((PG(k,J(k,2))/PN(J(k,2),J(k,2))*entaN(k,:)*Ad^(-1))*Pl)*ones(N,1);RG_l(k)=PG_l(k)*EG(k);
部分内容源自网络，侵权联系删除！
欢迎感兴趣的小伙伴点击文末阅读原文获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！
更多创新性预测及优化程序请点击公众号首页：《预测优化》|《综合能源》链接！
