資源描述：
《短路故障計(jì)算》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、一、三相短路故障計(jì)算MATLAB程序n=4;%獨(dú)立節(jié)點(diǎn)數(shù)nl=3;%支路數(shù)B1=[130.51i010;230.59i010;341.43i010];%線路參數(shù)形成的矩陣X=[10.2i;24i;30;40];%對(duì)地阻抗形成的矩陣V0=[1;1;1;1];%各節(jié)點(diǎn)的初電壓標(biāo)幺值形成的列矩陣D=[40];%短路點(diǎn)阻抗組成的矩陣NF=1;%短路點(diǎn)數(shù)目B=[0;0;0;1];%常數(shù)項(xiàng)矩陣(短路點(diǎn)對(duì)應(yīng)為1,其余節(jié)點(diǎn)對(duì)應(yīng)為0)Y=zeros(n);%初始化節(jié)點(diǎn)導(dǎo)納矩陣fori=1:n%初始化節(jié)點(diǎn)是否非直接接地ifX(i,2)~=0;%判斷i節(jié)點(diǎn)是否非直接接地p=
2、X(i,1);Y(p,p)=1./X(i,2);%算自導(dǎo)中接地阻抗的導(dǎo)納endend%形成節(jié)點(diǎn)導(dǎo)納矩陣fori=1:nlifB1(i,6)==0%判斷變壓器的變比是否在低壓側(cè)(it=0說明i側(cè)為低壓側(cè)；it=1說明i側(cè)為高壓側(cè))p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./B1(i,3)*B1(i,5);%節(jié)點(diǎn)互導(dǎo)納Y(q,p)=Y(p,q);%節(jié)點(diǎn)互導(dǎo)納Y(p,p)=Y(p,p)+1./B1(i,3)*B1(i,5)^2+B1(i,4)./2;%節(jié)點(diǎn)導(dǎo)納=自導(dǎo)納+X修正
3、量Y(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./2;enddisp('導(dǎo)納矩陣');disp(Y);%輸出導(dǎo)納矩陣A=Y;[n,m]=size(A);%解線性方程組，形成矩陣因子表Afori=1:nA(i,i)=1./A(i,i);forj=i+1:nA(i,j)=A(i,j)*A(i,i);%矩陣規(guī)格化endfork=i+1:nforj=i+1:nA(k,j)=A(k,j)-A(k,i)*A(i,j);endendenddisp('矩陣A的因子表為');disp(A);B=input('請(qǐng)輸入常數(shù)項(xiàng)矩陣(短路點(diǎn)對(duì)應(yīng)為1,其余節(jié)點(diǎn)
4、對(duì)應(yīng)為0),:B=');%利用因子表對(duì)常數(shù)項(xiàng)進(jìn)行的前代過程（修正矩陣B）fori=1:nB(i)=B(i)*A(i,i);forj=i+1:nB(j)=B(j)-A(j,i)*B(i);endenddisp('利用因子表對(duì)常數(shù)項(xiàng)進(jìn)行會(huì)帶的結(jié)果為:B=');disp(B)%利用因子表的上三角回帶過程fori=n-1:-1:1forj=i+1:-1:2B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1);endenddisp('在因子表的基礎(chǔ)上求解線性方程組的解為：X=');disp(B);%解出各節(jié)點(diǎn)電壓即得到Zik，此時(shí)B為列向量V0=inp
5、ut('請(qǐng)輸入由各節(jié)點(diǎn)的初始電壓標(biāo)幺值形成的列矩陣：V0=');D=input('請(qǐng)輸入由短路號(hào)、短路點(diǎn)阻抗組成的矩陣：D=');NF=input('請(qǐng)輸入短路點(diǎn)的數(shù)目：NF=');Z=zeros(n);V=zeros(n);I=zeros(nl);fork=1:NF%求各短路點(diǎn)的電流標(biāo)幺值disp('短路點(diǎn)與其他各節(jié)點(diǎn)的互阻抗');fori=1:nZ(i,D(k,1))=B(i,1);disp(Z(i,D(k,1)));endI(D(k,1),D(k,1))=V0(D(k,1),1)./(Z(D(k,1),D(k,1))+D(k,2));ft=num
6、2str(D(k,1));ts1=('點(diǎn)短路時(shí)');ts2=('電流的標(biāo)幺值If=');dn=strcat(ft,ts1,ts2);disp(dn);disp(I(D(k,1),D(k,1)));fori=1:nV(i,i)=V0(i,1)-I(D(k,1),D(k,1))*Z(i,D(k,1));%求各節(jié)點(diǎn)的電壓標(biāo)幺值電壓end%求非接地支路的短路電流標(biāo)幺值fori=1:nlifB1(i,6)==0%判斷該支路是否為含有變壓器支路且變比在低壓側(cè)k=B1(i,5);elsek=1./B1(i,5);endp=B1(i,1);q=B1(i,2);I(i,
7、i)=(V(p,p)-V(q,q)./k)./B1(i,3);%各支路短路電流=支路殘壓/支路阻抗(標(biāo)幺值)enddisp('各節(jié)點(diǎn)的電壓標(biāo)幺值U為(節(jié)點(diǎn)號(hào)從小到大排)：');fori=1:ndisp(V(i,i));enddisp('各接地支路短路電流的標(biāo)幺值I為(順序同您輸入B時(shí)一樣');fori=1:nifX(i,2)~=0;%判斷i節(jié)點(diǎn)是否為直接接地e=0;b=X(i,2);Ii0=(e-V(i,i))./b;disp(Ii0);endenddisp('各非接地支路短路電流的標(biāo)幺值I為(順序同您輸入B時(shí)一樣)：');fori=1:nldisp(
8、I(i,i));endend二、程序輸出結(jié)果導(dǎo)納矩陣0-6.9608i00+1.9608i00
