💥💥💞💞欢迎来到本博客❤️❤️💥💥
🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。
⛳️座右铭:行百里者,半于九十。
📋📋📋本文目录如下:🎁🎁🎁
目录
💥1 概述
📚2 运行结果
🎉3 参考文献
🌈4 Matlab代码、数据、文章讲解
💥1 概述
文献来源:
本文提出了三相配电系统的线性潮流。考虑平衡和非平衡运行以及负载的ZIP模型。该方法不需要任何与R/X率。尽管它很简单,但与传统的后向扫描算法相比,它非常准确。
直流功率流是电力系统分析和运行研究最多的方法之一[1]。然而,这种线性近似不适用于配电系统,因为它们高R/X比率和不平衡操作[2]。这封信通过在复平面上使用线性近似来解决这个问题。不需要径向拓扑。不考虑光伏节点,但在电网规范强制这些发电机单位功率因数运行的情况下,可以包括配电发电机。
📚2 运行结果
部分代码:
function ShowResults(Res,Feeder);
fprintf('---------------%s---------------\n',Feeder.Options.Name);
fprintf('Total Losses : %5.4f kW\n',Res.Perd);
%fprintf('N Iterations : %i \n',Res.iter);
if Feeder.Options.DeltaLoadFlow
%% Voltajes de linea
disp('------- LINE TO LINE VOLTAGES ----------------------------------------------------------');
format short
fprintf('NODE\tVAB(pu)\t VAB(deg)\t VBC(pu)\tVBC(deg)\tVCA(pu)\tVCA(deg)\n');
for k = 1:Feeder.NumN
V = abs(Res.Vpu_line(k,:));
A = angle(Res.Vpu_line(k,:))*180/pi;
n = Feeder.Nodes_ID(k);
if (n<10)
fprintf('N%i \t',Feeder.Nodes_ID(k));
else
if (n<100)
fprintf('N%i \t',Feeder.Nodes_ID(k));
else
fprintf('N%i\t',Feeder.Nodes_ID(k));
end
end
for m = 1:3
if (V(m)>0.1)
fprintf('%5.4f < ',V(m));
if (A(m)>0)
fprintf(' %5.4f \t ',A(m));
else
fprintf('%5.4f \t ',A(m));
end
else
fprintf('. \t ')
end
end
fprintf('\n');
end
else
% Voltajes de fase
disp('--------PHASE VOLTAGES ---------------------------------------------------------');
format short
fprintf('NODE\tVAn(pu)\t VAn(deg)\t VBn(pu)\tVBn(deg)\tVCn(pu)\tVCn(deg)\n');
for k = 1:Feeder.NumN
V = abs(Res.Vpu_phase(k,:));
A = angle(Res.Vpu_phase(k,:))*180/pi;
n = Feeder.Nodes_ID(k);
if (n<10)
fprintf('N%i \t',Feeder.Nodes_ID(k));
else
if (n<100)
fprintf('N%i \t',Feeder.Nodes_ID(k));
else
fprintf('N%i\t',Feeder.Nodes_ID(k));
end
end
for m = 1:3
if (V(m)>0)
fprintf('%5.4f < ',V(m));
if (A(m)>0)
fprintf(' %5.4f \t ',A(m));
else
fprintf('%5.4f \t ',A(m));
end
else
fprintf('. \t ')
end
end
fprintf('\n');
end
end
% %disp([Feeder.Nodes_ID,abs(Res.Vpu_line(:,1)),angle(Res.Vpu_line(:,1))*180/pi,abs(Res.Vpu_line(:,2)),angle(Res.Vpu_line(:,2))*180/pi,abs(Res.Vpu_line(:,3)),angle(Res.Vpu_line(:,3))*180/pi]);
%disp('Reguladores')
%lin = Feeder.Regulators(:,1);
%N1 = Feeder.Nodes_ID(Feeder.Topology(lin,1));
%N2 = Feeder.Nodes_ID(Feeder.Topology(lin,2));
%disp([N1 N2]);
%Res.Vpu_phase = conj(Vs')/Vbase_fase;
% Res.Vpu_line = Vlin/Vbase_lin;
% Res.Inodes = -conj(Inodal'); % entrando al nodo
% Res.Ilines = conj(Is');
% Res.iter = iter;
% Res.err = err;
% Res.Perd = Sperdidas;
% Res.Sort = Sort;
🎉3 参考文献
部分理论来源于网络,如有侵权请联系删除。