電力系統(tǒng)三相短路計(jì)算地MATLAB代碼

1、實(shí)用文檔電力系統(tǒng)三相短路計(jì)算a main.mclear tim%丁開文件dfile,pathname=uigetfile('*m', 'Select Data File');if pathname = 0 error( ' you must select a valid data file')else lfile =length(dfile); eval(dfile(1:lfile-2); end%定義輸出文件output_file=fopen('output.dat' , 'w');%開始計(jì)時(shí)tic;%求解節(jié)點(diǎn)導(dǎo)

2、納£!陣，其中 Ymatrix1 是考慮了變比，且支路未近似的導(dǎo)納矩陣；Ymatrix2 是近似變比為1 ,但是支路未近似計(jì)算的節(jié)點(diǎn)導(dǎo)納矩陣；Ymatrix3 是近似變比為1 ,采取近似支路參數(shù)1的導(dǎo)納矩陣；Ymatrix4 是近似變比為1,采取近似支路參數(shù) 2的導(dǎo)納矩陣。Y = Ymatrix2(bus,line);%寸故障點(diǎn)進(jìn)行導(dǎo)納修正fixY = FixY(Y,bus,fault);%求注入電流Iinj = Inode(bus,calcSettings);U = fixYIinj;%得到故障支路與其他支路電流Bcurrent = Ibranch( line,U,fault,Y

3、);%如果發(fā)生支路三相短路，那么對(duì)應(yīng)該支路的電流修正為-999999-j999999Ib = ReviseBcurrent( fault,Bcurrent ); %結(jié)束計(jì)時(shí) tim=toc;fprintf( 'n程序運(yùn)行結(jié)果)； fprintf( 'n計(jì)算完成，共用時(shí) %4.4fs,相關(guān)結(jié)果已保存在 output.datn' ,tim);%輸出結(jié)果fprintf_result(output_file, Ib);fprintf_result1(Ib);b FixY.mfunction fixY = FixY( Y,bus,fault ) %寸形成的導(dǎo)納矩陣進(jìn)行故障點(diǎn)的修正

4、 nb,mb=size(bus);nf,mf= size(fault);文案大全用matlab實(shí)現(xiàn)電力系統(tǒng)潮流計(jì)算晏鳴宇fixY = Y;%寸發(fā)電機(jī)節(jié)點(diǎn)導(dǎo)納修正for k=1:nbbusType=bus(k,7);if (busType=1)fixY(bus(k,1),bus(k,1) = fixY(bus(k,1),bus(k,1) + 1/1i/bus(k,8); end end %寸節(jié)點(diǎn)短路和支路短路的導(dǎo)納矩陣進(jìn)行修正 for k=1:nf nodeI=fault(k,1);nodeJ=fault(k,2);dis=fault(k,3);if (nodeI=0)fixY(nodeJ,no

5、deJ) = 999999+1i*999999; continue ;endif (nodeJ=0)fixY(nodeI,nodeI) = 999999+1i*999999;continue ;endif (dis=0)&&(nodeI*nodeJ=0)fixY(nodeI,nodeI) = 999999+1i*999999;continue ;endif (dis=1)&&(nodeI*nodeJ=0)fixY(nodeJ,nodeJ) = 999999+1i*999999;continue ;endif (dis=1)&&(dis=0)&am

6、p;&(nodeI*nodeJ=0)fixY(nodeI,nodeI) = fixY(nodeI,nodeI) - fixY(nodeI,nodeJ)/dis;fixY(nodeJ,nodeJ) = fixY(nodeJ,nodeJ) - fixY(nodeI,nodeJ)/(1-dis);fixY(nodeI,nodeJ)=0;fixY(nodeJ,nodeI)=0;endend endc fprintf_result.mfunction output_args = fprintf_result( output_file, Ib )%a得到的短路電流輸入到輸出文件中n,m=size(

7、Ib);fprintf( output_file,' No. No. vector of I value of In');用matlab實(shí)現(xiàn)電力系統(tǒng)潮流計(jì)算晏鳴宇for k=1:nI=Ib(k,1)；J=Ib(k,2);I01=real(Ib(k,3);I02=imag(Ib(k,3);I1=Ib(k,4);if (I02>=0) fprintf( output_file, endif (I02<0) I02=abs(I02); fprintf( output_file, endfprintf( output_file, endend'%3d %3d %10

8、.6f+j%10.6f %10.6f''%3d %3d %10.6f-j%10.6f %10.6f''n');,I,J,I01,I02,I1)；,I,J,I01,I02,I1)；fprintf_result1.mfunction output_args = fprintf_result1( Ib )%UNTITLED ?u ' ?' | e ? e ?o e y? o a n,m=size(Ib);fprintf( ' No. No. vector of I value of In');for k=1:nI=Ib(k,1)

9、；J=Ib(k,2);I01=real(Ib(k,3);I02=imag(Ib(k,3);I1=Ib(k,4);if (I02>=0)fprintf('%3d %3d %10.6f+j%10.6f %10.6f',I,J,I01,I02,I1);end if (I02<0) I02=abs(I02);fprintf('%3d %3d %10.6f-j%10.6f %10.6f',I,J,I01,I02,I1);end fprintf('n');end endfunctionIbranch.mBcurrent = Ibranch( li

10、ne,U,fault,Y )用matlab實(shí)現(xiàn)電力系統(tǒng)潮流計(jì)算晏鳴宇% 計(jì)算短路電流%記錄短路故障參數(shù)，如短路節(jié)點(diǎn)，如為支路短路，記錄距離節(jié)點(diǎn)的距離 %t匕段計(jì)算采用的支路參數(shù)未近似，如果計(jì)算近似的時(shí)候需要修改 nl,ml=size(line);Bcurrent=zeros(nl+1,4);faultI=fault(1,1);faultJ=fault(1,2);dis=fault(1,3);faultNode = 0;if (faultI=0)faultNode = faultJ;endif (faultJ=0)faultNode = faultI;endif (dis=1)&&

11、;(faultI*faultJ=0) faultNode = faultJ;endif (dis=0)&&(faultI*faultJ=0) faultNode = faultI;endif (faultNode=0)Bcurrent(nl+1,1) = faultNode;Bcurrent(nl+1,2) = faultNode;Iij = 0;Iij1=0;end%計(jì)算非故障支路的短路電流for k=1:nli=line(k,1);j=line(k,2);Ui=U(i);if j=0Uj=U(j);elseUj=0;endif line(k,2)=0Ym=line(k,5)

12、+1i*line(k,6);Iij=Ui*Ym;Iij1=abs(Iij);endif line(k,2)=0Zt=line(k,3)+1i*line(k,4);用matlab實(shí)現(xiàn)電力系統(tǒng)潮流計(jì)算晏鳴宇Yt=1/Zt;Ym=line(k,5)+1i*line(k,6);Iij=(Ui-Uj)*Yt+Ui*Ym;Iij1=abs(Iij);endBcurrent(k,1)=i;Bcurrent(k,2)=j;Bcurrent(k,3)=Iij;Bcurrent(k,4)=Iij1;end%如果為節(jié)點(diǎn)短路，修正短路點(diǎn)的電流大小if (faultNode=0)Bcurrent(nl+1,1) = f

13、aultNode;Bcurrent(nl+1,2) = faultNode;Ifault = 0;branchCurrent=0;for k=1:nlI=line(k,1);J=line(k,2);if (I*J=0) continue ;endbranchCurrent = (U(I)-U(J)/(line(k,3)+1i*line(k,4); if (I=faultNode)Ifault = Ifault - branchCurrent ;else if (J=faultNode) Ifault = Ifault + branchCurrent ; endendendBcurrent(nl

14、+1,3) = Ifault;Bcurrent(nl+1,4) = abs(Bcurrent(nl+1,3);end%如果為支路短路，修正短路支路的短路電流大小 if (dis=0)&&(dis=1)&&(faultI*faultJ=0) Bcurrent(nl+1,1) = faultI;Bcurrent(nl+1,2) = faultJ;Bcurrent(nl+1,3) = U(faultI)*Y(faultI,faultJ)/dis + U(faultJ)*Y(faultI,faultJ)/(1-dis);Bcurrent(nl+1,4) = abs(Bc

15、urrent(nl+1,3);endend用matlab實(shí)現(xiàn)電力系統(tǒng)潮流計(jì)算晏鳴宇f Inode.mfunction Iinj = Inode( bus,calcSettings ) % 計(jì)算節(jié)點(diǎn)注入電流 nb,mb=size(bus);Iinj = zeros(nb,1);for k=1:nbbusType=bus(k,7);if (calcSettings=1)v = 1; elsev = bus(k,2); end%對(duì)發(fā)電機(jī)節(jié)點(diǎn)電流進(jìn)行修正if (busType=1)Iinj(bus(k,1),1) = Iinj(bus(k,1),1) + v/1i/bus(k,8); end end

16、endg ReviseBcurrent.mfunction Ib = ReviseBcurrent( fault,Bcurrent )%如果發(fā)生支路短路，對(duì)原來的計(jì)算電流進(jìn)行修正，使該支路短路電流輸出為 clear faultl faultJ dis nt,mt=size(Bcurrent);Ib=zeros(nt,mt);faultI=fault(1,1);faultJ=fault(1,2); dis=fault(1,3);for k=1:nt-1 i=Bcurrent(k,1); j=Bcurrent(k,2); Ib(k,:)=Bcurrent(k,:);if (faultI*fault

17、J=0)&&(dis=1)&&(dis=0)&&(i=faultI)&&(j=faultJ) Ib(k,1)=i； Ib(k,2)=j;Ib(k,3)=-999999-1i*999999;Ib(k,4)=-999999;endif (faultI*faultJ=0)&&(dis=1)&&(dis=0)&&(i=faultJ)&&(j=faultI) Ib(k,1)=i;Ib(k,2)=j;Ib(k,3)=-999999-1i*999999;-999999-j99999

18、9Ib(k,4)=-999999;用matlab實(shí)現(xiàn)電力系統(tǒng)潮流計(jì)算晏鳴宇endIb(nt,:尸Bcurrent(nt,:); endh Ymatrixl.mfunction Y = Ymatrix1( bus,line )%考慮變壓器，并且支路參數(shù)不近似的節(jié)點(diǎn)導(dǎo)納矩陣nb,mb=size(bus);nl,ml=size(line);Y=zeros(nb,nb);for k=1:nlI=line(k,1);J=line(k,2);Zt=line(k,3)+1i*line(k,4);Yt=1/Zt;Ym=line(k,5)+1i*line(k,6);K=line(k,7);if (K=0)&am

19、p;&(J=0)Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt+Ym;Y(I,J尸Y(I,J)-Yt;Y(J,I)=Y(I,J);endif (K=0)&&(J=0)Y(I,I)=Y(I,I)+Ym;endif K>0Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt/(K*K);Y(I,J)=Y(I,J)-Yt/K;Y(J,I)=Y(I,J);endif K<0Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J尸丫(J,J)+K*K*Yt;Y(I,J)=Y(I,J)+K*Yt;Y(J,I)=Y(I,J);en

20、dend endi Ymatrix2.mfunctionY = Ymatrix2( bus,line )用matlab實(shí)現(xiàn)電力系統(tǒng)潮流計(jì)算晏鳴宇%考慮變壓器變比近似為1,支路參數(shù)不等效nb,mb=size(bus);nl,ml=size(line);Y=zeros(nb,nb);for k=1:nlI=line(k,1);J=line(k,2);Zt=line(k,3)+1i*line(k,4);Yt=1/Zt;Ym=line(k,5)+1i*line(k,6);if J=0Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J尸Y(J,J)+Yt+Ym;Y(I,J尸Y(I,J)-Yt;Y(J,I)=Y(I,J);endif J=0Y(I,I)=