duichen.m

可以用来进行电力系统短路计算中对称短路计算。

%Z是网络的阻抗矩阵
Z=[j*0.15   j*0.1   j*0.1
   j*0.1    j*0.075 j*0.1
   j*0.1    j*0.1    0    ]
N1=3; %是网络的节点数
k=3 ; %k是短路的节点号
 Y=zeros(3,3);                  %形成3*3的零矩阵
 for i=1:N1
     for p=1:N1
         if Z(i,p)~=0
         Y(i,i)=Y(i,i)+1/Z(i,p);%自导纳等于所有连接在该节点上的导纳之和
            if  i~=p
        Y(i,p)=-1/Z(i,p);       %互导纳等于两节点间导纳的负值
            end
         end
 end
end
 ZZ=eye(N1)/Y  ;       %对节点导纳矩阵求逆得到节点阻抗矩阵
 for p=1:N1
     if p==k;
         I(p)=1;
     else 
         I(p)=0;
     end 
 end
 ZZf=ZZ*I';
 If=1/ZZf(k)          %求得短路电流
 %计算各点的故障电压分量
 for p=1:N1
     if p==k
         II(p)=-If;
     else 
         II(p)=0;   
     end
 end
 vv=ZZ*II'
 
 %假设故障前空载，各节点电压为1
 for p=1:N1
     vvv(p)=1
 end
 %计算各点故障后的电压值：正常运行的电压分量+故障分量
 for p=1:N1
     if p==k
         U(p)=0;   %短路点的电压为零
     else 
         U(p)=vvv(p)-vv(p)
     end
 end
 %计算各支路电流
 for m=1:N1
     for n=1:N1
         if m<n
             I(m,n)=(U(m)-U(n))/Z(m,n)  %支路电流等于节点的电压差除以支路阻抗
         end
     end
 end