电力系统三相潮流matlab编程实现.7z

function [delPQ] =deltaPQ（Ybus）         %计算功率偏差
G=real（Y）;
B=imag（Y）;
[nymy]=size（Y）;
[nbmb]=size（bus）;
nSW=0;                        %计数：平衡节点个数
nPV=0;                        %计数：PV节点个数
nPQ=0;                        %计数：PQ节点个数
for I=1:nb                    % nb为总节点数
    type=bus（I6）;
    if type==3
       nSW=nSW+1;             % increment swing bus counter
    elseif type==2
       nPV=nPV +1;            % increment PV bus counter
    else
       nPQ=nPQ+1;             % increment PQ bus counter
    end
end
ang=bus（:3）;                 %把bus第3列相角存放到矩阵ang中
V=bus（:2）;                   %把bus第2列相角存放到矩阵V中
A=zeros（nymy）;               %矩阵A存放节点相角差
P=bus（:4）;                   %把bus第4列输入有功功率存放到矩阵P中
Q=bus（:5）;                   %把bus第5列输入无功功率存放到矩阵Q中
for k=1:ny                    %求各节点相角差存放到矩阵A中
    A（k:）=ang（k）-ang（:）;
end
for j=1:（nPQ+nPV）             %计算PQ、PV节点的有功功率偏差dP存放到矩阵delP中
    delP（j）=P（j）;
    for k=1:ny
        delP（j）=delP（j）-V（j）*V（k）*（G（jk）*cos（A（jk））+B（jk）*sin（A（jk）））;
    end
end
for j=1:nPQ                  %计算PQ节点的无功功率偏差dQ存放到矩阵delQ中
    delQ（j）=Q（j）;
    for k=1:ny
        delQ（j）=delQ（j）-V（j）*V（k）*（G（jk）*sin（A（jk））-B（jk）*cos（A（jk）））;
    end
end
delPQ=[delPdelQ]‘;          %按dP、dQ节点的顺序形成一个单列矩阵存放到矩阵delPQ中