chaoliu.m

可用传输容计算程序

function [O,PL]=chaoliu(D,P,Xl)
B(30,30)=0;%定义节点导纳矩阵为一个30X30的零矩阵；
for i1=1:length(D)
    B(D(i1,2),D(i1,3))=D(i1,6); B(D(i1,3),D(i1,2))=D(i1,6);%定义导纳矩阵非对角元素即互导纳；
end
for i2=1:30%定义导纳矩阵对角元素；
    for i3=1:length(D)
        if i2==D(i3,2)|i2==D(i3,3)
            B(i2,i2)=B(i2,i2)-D(i3,6);
        end
    end
end
B(13,:)=[];B(:,13)=[];%消去导纳矩阵平衡节点所在的行和列；
P(13,:)=[];%消去平衡节点注入功率；
P1=P(:,2)-P(:,3);%节点注入功率减去负荷功率；
A1=inv(B);
O=A1*P1;%直流潮流模型算出各节点相角
O=[O(1:12);0;O(13:29)];
Xm=diag([Xl]);%Xm为线路电抗对角矩阵；
for i4=1:length(D)%M为支路节点关系矩阵起始节点为1，终止节点为-1；
    M(i4,D(i4,2))=1;M(i4,D(i4,3))=-1;
end
PL=inv(Xm)*M*O;%PL为各支路的线路功率；