潮流计算matlab程序

这个压缩包包含了从主程序到Y矩阵到PQ节点再到最后的算法，使用了牛顿拉夫逊计算法，理论上可运行4，14，22，30，118，300，1047节点（因为1047节点可能要算半天也算不出来），可运行。

% 功能：计算不平衡量
% 日期：2017年11月17日
% 名字：卢煜
function [deltaPdeltaQdeltaPQ]=CalculatePQ（YfaVbVanodeNumbalanceNodepvnodeNumPQ）
deltaP=zeros（nodeNum1）;
deltaQ=zeros（nodeNum1）;
deltaPQ=zeros（nodeNum1）;
for i=1:nodeNum
    if i==balanceNode
    continue
    end
    for j=1:nodeNum
    deltaP（i）=deltaP（i）+Yf（ij）*Vb（j）*cos（Va（i）-Va（j）-a（ij））;
    deltaQ（i）=deltaQ（i）+Yf（ij）*Vb（j）*sin（Va（i）-Va（j）-a（ij））;
    end
    deltaP（i）=deltaP（i）*（-Vb（i））;
    deltaP（i）=deltaP（i）+P（i1）;
    deltaQ（i）=deltaQ（i）*（-Vb（i））;
    deltaQ（i）=deltaQ（i）+Q（i1）;
end
%% 处理PV节点
for i=1length（pvnodeNum）
deltaQ（pvnodeNum1）=0;
end
%% 处理平衡节点
deltaP（balanceNode1）=0;
deltaQ（balanceNode1）=0;
%% 形成不平衡量矩阵
for i=1:nodeNum
    deltaPQ（i）=deltaP（i）;
    deltaPQ（nodeNum+i）=deltaQ（i）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-01-11 22:55  潮流计算程序\
     文件         844  2018-01-12 00:48  潮流计算程序\CalculatePQ.m
     文件        1093  2018-01-12 00:48  潮流计算程序\CalculateS.m
     文件        1760  2018-01-12 00:48  潮流计算程序\FormJacobi.m
     文件        1406  2018-01-12 00:47  潮流计算程序\FormY.m
     文件      185785  2011-12-19 10:03  潮流计算程序\IEEE1047.dat
     文件       16634  2013-12-18 12:18  潮流计算程序\IEEE118.dat
     文件        1110  2013-11-21 09:41  潮流计算程序\IEEE14.dat
     文件        1729  2017-12-30 11:00  潮流计算程序\IEEE22.dat
     文件        4017  2013-11-21 09:41  潮流计算程序\IEEE30.dat
     文件       61718  2011-12-19 10:03  潮流计算程序\IEEE300.dat
     文件         321  2011-12-19 10:02  潮流计算程序\IEEE4.dat
     文件        1266  2018-01-12 00:57  潮流计算程序\MainPragram.m
     文件        2387  2018-01-12 00:47  潮流计算程序\ReadData.m