sub2.asv

可用传输容计算程序

function [A1duan,Aduan,Tpduan]=sub2(D,A1,A,P1,n)
%A1duan为断线后的节点导纳矩阵的逆，即节点阻抗阵；
%Aduan为断线后的线路有功与节点注入有功的关系矩阵；
%Tpduan为断线后的各线路功率；
%n为断线支路；
S(30,41)=0;%S为节点和支路关系矩阵，线路初始节点为1，终止节点为-1，不相关节点值为0；
for i1=1:length(D)
    S(D(i1,2),i1)=1;S(D(i1,3),i1)=-1;
end
S(13,:)=[];%去除参考节点；
a=D(n,4)-S(:,n)'*A1*S(:,n);
A1duan=A1+A1*S(:,n)*(1/a)*S(:,n)'*A1;
[k1,k2]=size(A);
for m=1:k1
    p=(1/a)*S(:,m)'*A1*S(:,n);
    Aduan(m,:)=A(m,:)+p*A(n,:)*(D(n,4)/D(m,4));
end
Tpduan=Aduan*P1;