sub1.m

可用传输容计算程序

function [Tp,A]=sub1(A1,P1,D)%A1为节点导纳矩阵的逆矩阵，P1为节点注入功率列向量(注入设为0)，D为支路数据矩阵；
A1=[A1(1:12,:);zeros(1,29);A1(13:29,:)];

for i1=1:length(D)
   
      Tp(i1,:)=(1/D(i1,4))*(A1(D(i1,2),:)-A1(D(i1,3),:))*P1;
      A(i1,:)=(1/D(i1,4))*(A1(D(i1,2),:)-A1(D(i1,3),:));
   
end