sub1.asv

可用传输容计算程序

function Tp=transfer1(A1,P1,D)%A1为节点导纳矩阵的逆矩阵，P1为节点注入功率列向量(注入设为0)，D为支路数据矩阵；
A1=[A1(1:12,:)]
for i1=1:length(D)
    if D(i1,2)~=13&D(i1,3)~=13
       Tp(i1,:)=(1/D(i1,4))*(A1(D(i1,2),:)-A1(D(i1,3),:))*P1;
    elseif D(i1,3)==13
            Tp(i1,:)=(1/D(i1,4))*A1(D(i1,2),:)*P1;
        
    elseif D(i1,2)==13
            Tp(i1,:)=(1/D(i1,4))*(-A1(D(i1,3),:))*P1;
       
    end
end