精选优质文档-----倾情为你奉上用追加支路法计算电力网节点阻抗矩阵佘名寰 编写电力网节点阻抗矩阵在电力系统短路电流计算中获得广泛运用本文通过例题介绍用追加支路法计算电力网节点阻抗矩阵的方法和程序，程序计及零序互感的影响程序用MATLAB语言编写，线路参数均采用标么值本文可供电气专业人员计算相关问题时参考例2.1】 图2-1所示为一个3节点网络，两台发电机，四条线路发电机一台中性点接地，另一台不接地发电机次暂态电抗和线路阻抗原始数据在表2-1、表2-2中给出，线路3、4间有零序互感电抗试计算该网络图的正序和零序节点阻抗矩阵①①②③23451G1图2-1例2.1网络图G④6表2-1 正序网络数据首端节点编号末端节点编号回路编号自感标么阻抗R (pu.)自感标么阻抗X (pu.) 4160.00.20004122132333123450.00.00.00.00.00.16000.08000.06000.06000.1300表2-2 零序网络数据首端节点编号末端节点编号回路编号自感标么阻抗R 自感标么阻抗X 互感标么阻抗RM互感标么阻抗XM4122132333123450.00.00.00.00.00.02000.14000.10000.12000.17000.00.00.05000.05002.3 用追加支路法形成节点阻抗矩阵在网络改变如增加或断开一条支路需重新计算短路电流时，若仍用节点导纳矩阵求逆矩阵的方法求新的节点阻抗矩阵，计算工作量比较大。
这时采用追加支路法形成节点阻抗矩阵比较简便追加支路法不需计算逆矩阵，对小型网络求阻抗矩阵尤为方便2.3.1 追加支路法形成节点阻抗矩阵的基本公式⑴ 对参考节点追加辐射支路：从参考节点到节点q引入一条阻抗为z的支路，q为新节点，该支与其它支路无耦合，则节点阻抗矩阵的元素Z=z, Zqi=Ziq=0 （2-22） ⑵ 追加一条辐射支路到一个新节点：从k节点到q节点追加一条阻抗为z支路，该支与其它支路无耦合，k是原有结点，q是新结点，k不是参考结点，p个结点原来已确定，则： Z=Zkk+z Ziq=Zik, i=1,2,…,p Zqi=Zki, i=1,2,…,p （2-23）⑶ 从k结点到参考结点追加一条链枝：该支路阻抗为z，与其它支路无耦合，k结点为原来已确定的p个结点中的一个，则：置q=p+1, Ziq=Zik, Zqi=Zki, i=1,2,…,p Z=Zkk+z （2-24）用克朗降阶法消去矩阵第q行和第q列；⑷ 追加一条链枝：从i结点到k结点追加一条支路，I,k都是原已确定结点，p是已经确定的结点总数，则： 置q=p+1, Zjq=Zjk-Zji, Zqj=Zkj-Zij, j=1,2,…,p Z=Zii+Zkk-Zik-Zki+z, （2-25）用克朗降阶法消去矩阵第q行和第q列；⑸ 追加一条有互感的支路： 零序网络节点阻抗矩阵计及线路间的互感时，其计算方法见参考文献⑧P333公式（12-38），（12-42）。
假定m,n和p是结点总数为p的网络中已确定的结点，从结点p到结点q追加一条阻抗为zpq支路，其与阻抗为zmn的支路相耦合，互感阻抗为zm,新形成的节点阻抗矩阵第q行元素为 Zqk=ZPK-(zm/zmn)(Zmk-Znk) (k=1,2,….p) Z=Zpq+zpq-(zm/zmn)( zm+Zmq-Znq) （2-26） 式中ZPK，Zmk，Znk 为原p阶阻抗矩阵元素，Zpq，Zmq，Znq为第q列元素 当追加一条有互感的链支pk时可分两步走假设一个虚拟节点q,先追加阻抗为zpk的树支pq,再追加阻抗为零无互感的链支qk,可套用相关公式2.3.2 追加支路法形成节点阻抗矩阵的M函数⑴ 追加无互感支路形成节点阻抗的M函数 faddbra.m function[zz]=faddbra(zb,nstart,nend) global n m p=1; for i=1:n+1 for j=1:n+1 zzm(i,j)=0.0; end end for i=1:n for j=1:n zz(i,j)=0.0; end end for i=1:n+1 nnod(i)=0; end for i=1:m if (nstart(i)>n) nstart(i)=0;end end for i=1:m k=nstart(i); L=nend(i); if (L==0) L=n+1;end if ((k==0)&(nnod(L)==0))==1 if (L>p) p=L;end nnod(L)=1; zzm(L,L)=zb(i); elseif ((k==0)&(nnod(L)==1))==1 q=p+1; for i1=1:p for j=1:p zzm(i1,q)=zzm(i1,L); zzm(q,j)=zzm(L,j); end end zzm(q,q)=zzm(L,L)+zb(i); for i1=1:p for j=1:p zzm(i1,j)=zzm(i1,j)-(zzm(i1,q)*zzm(q,j)/zzm(q,q)); end end for i1=1:p for j=1:p zzm(i1,q)=0.0 ; zzm(q,j)=0.0 ; end end zzm(q,q)=0.0 ; elseif ((nnod(k)==1)&(nnod(L)==0)&(L~=n+1))==1 if (L>p) p=L;end nnod(L)=1; for i1=1:p for j=1:p zzm(i1,L)=zzm(i1,k); zzm(L,j)=zzm(k,j); end end zzm(L,L)=(zzm(k,k)+zb(i)); elseif ((nnod(k)==1)&(nnod(L)==1))==1 q=p+1; for i1=1:p for j=1:p zzm(i1,q)=zzm(i1,k)-zzm(i1,L); zzm(q,j)=zzm(k,j)-zzm(L,j); end end zzm(q,q)=(zzm(k,k)+zb(i))+zzm(L,L)-2*zzm(k,L); for i1=1:p for j=1:p zzm(i1,j)=zzm(i1,j)-(zzm(i1,q)*zzm(q,j)/zzm(q,q)); end end for i1=1:p for j=1:p zzm(i1,q)=0.0; zzm(q,j)=0.0; end end zzm(q,q)=0.0; elseif ((nnod(k)==1)&(L==n+1))==1 q=p+1; for i1=1:p for j=1:p zzm(i1,q)=zzm(i1,k); zzm(q,j)=zzm(k,j); end end zzm(q,q)=zzm(k,k)+zb(i); for i1=1:p for j=1:p zzm(i1,j)=zzm(i1,j)-(zzm(i1,q)*zzm(q,j)/zzm(q,q)); end end for i1=1:p for j=1:p zzm(i1,q)=0.0 ; zzm(q,j)=0.0 ; end end zz。