1、.%線路數(shù)據(jù) Line=1-節(jié)點(diǎn)編號(hào) 2-線路首端節(jié)點(diǎn)號(hào) 3-線路末端節(jié)點(diǎn)號(hào) 4-支路電阻 5-支路電抗 6-支路電納注意：此處取的是B/2Line=1 1 2 0.04 0.25 0.25 2 1 3 0.1 0.35 0 3 2 3 0.08 0.30 0.25 ; %變壓器數(shù)據(jù) transform=1-支路編號(hào) 2-支路首節(jié)點(diǎn)編號(hào) 3-支路末節(jié)點(diǎn)編號(hào) 4-支路電阻p.u. 5-支路電抗p.u. 6-變壓器變比p.u.transform= 1 2 4 0 0.015 1.052 3 5 0 0.03 1.05;%數(shù)據(jù)預(yù)處理 Nbus=5; % 節(jié)點(diǎn)數(shù) nline=sizeLine,1; %

2、 線路個(gè)數(shù) ntrans=sizetransform,1; % 變壓器個(gè)數(shù) slack=5; % 平衡節(jié)點(diǎn)號(hào) Npq=4; % PQ節(jié)點(diǎn)的個(gè)數(shù) Ss=0;%計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Y=zerosNbus;if nline>=1%判斷是否存在線路 for k=1:nline %以下處理線路 t1=Linek,2; t2=Linek,3; b2=Linek,6;%分別取出線路的首端節(jié)點(diǎn)編號(hào)t1、末端節(jié)點(diǎn)編號(hào)t2和對(duì)地電納b2 Yl=1/Linek,4+j*Linek,5;%計(jì)算線路的支路電導(dǎo)Yl Yt1,t1=Yt1,t1+Yl+j*b2;%修正第k條線路首端節(jié)點(diǎn)的自導(dǎo)納 Yt1,t2=Yt1,t2

3、-Yl;%修正第k條線路首端節(jié)點(diǎn)與末端節(jié)點(diǎn)之間的互導(dǎo)納 Yt2,t1=Yt2,t1-Yl;%修正第k條線路末端節(jié)點(diǎn)與首端節(jié)點(diǎn)之間的互導(dǎo)納 Yt2,t2=Yt2,t2+Yl+j*b2;%修正第k條線路末端節(jié)點(diǎn)的自導(dǎo)納 endendif ntrans>=1%判斷是否存在變壓器 for k=1:ntrans %以下處理變壓器 t1=transformk,2;t2=transformk,3;t3=transformk,6;%分別取出變壓器的首端節(jié)點(diǎn)編號(hào)t1、末端節(jié)點(diǎn)編號(hào)t2和變比t3 Yt=1/transformk,4+j*transformk,5; Yt1=Yt/t3; Yt2=Yt*1-t3

4、/t3*t3; Yt3=Yt*t3-1/t3; Yt1,t1=Yt1,t1+Yt1+Yt2; Yt1,t2=Yt1,t2-Yt1; Yt2,t1=Yt2,t1-Yt1; Yt2,t2=Yt2,t2+Yt1+Yt3; endendG=realY;B=imagY; %區(qū)分節(jié)點(diǎn)導(dǎo)納矩陣的實(shí)部和虛部GB%賦初值delt1=0;delt2=0;delt3=0;delt4=0;u1=1;u2=1;u3=1;u4=1;p1=-0.30;q1=-0.18;p2=-0.55;q2=-0.13;p3=0;q3=0;p4=0.5;q4=1.10;p5=0.8;q5=0.50;k=0;precision=1;Npq=

5、4; %Npq分別是網(wǎng)絡(luò)中的PQ節(jié)點(diǎn)數(shù)%Unbalance=-JacobiCorrectionwhile precision>0.00001 %設(shè)定誤差上限，判斷是否繼續(xù)迭代 u5=1.06; delt5=0; %設(shè)定平衡節(jié)點(diǎn)電壓相角與幅值 k; u; delt; for m=1:Npq for n=1:Nbus ptn=um*un*Gm,n*cosdeltm-deltn+Bm,n*sindeltm-deltn; %由節(jié)點(diǎn)電壓求得的PQ節(jié)點(diǎn)注入有功功率 qtn=um*un*Gm,n*sindeltm-deltn-Bm,n*cosdeltm-deltn; %由節(jié)點(diǎn)電壓求得的PQ節(jié)點(diǎn)注入無功

6、功率 end Unbalance2*m-1=pm-sumpt; %計(jì)算PQ節(jié)點(diǎn)有功功率不平衡量 Unbalance2*m =qm-sumqt; %計(jì)算PQ節(jié)點(diǎn)無功功率不平衡量end %Unbalance是節(jié)點(diǎn)不平衡量矩陣 for m=1:Npq for n=1:Nbus h0n= um*un*Gm,n*sindeltm-deltn-Bm,n*cosdeltm-deltn; n0n=-um*un*Gm,n*cosdeltm-deltn+Bm,n*sindeltm-deltn; j0n=-um*un*Gm,n*cosdeltm-deltn+Bm,n*sindeltm-deltn; l0n=-um*

7、un*Gm,n*sindeltm-deltn-Bm,n*cosdeltm-deltn; end Hm,m=sumh0-um2*Gm,m*sindeltm-deltm-Bm,m*cosdeltm-deltm; Nm,m=sumn0+um2*Gm,m*cosdeltm-deltm+Bm,m*sindeltm-deltm-2*um2*Gm,m; Jm,m=sumj0+um2*Gm,m*cosdeltm-deltm+Bm,m*sindeltm-deltm; Lm,m=suml0+um2*Gm,m*sindeltm-deltm-Bm,m*cosdeltm-deltm+2*um2*Bm,m; Jacobi

8、2*m-1,2*m-1=Hm,m; Jacobi2*m-1,2*m =Nm,m; Jacobi2*m ,2*m-1=Jm,m; Jacobi2*m ,2*m =Lm,m; end %計(jì)算m=n情況下的Jacobi矩陣中的子矩陣元素 for m=1:Npq for n=1:Npq if m=n else Hm,n=-um*un*Gm,n*sindeltm-deltn-Bm,n*cosdeltm-deltn; Jm,n= um*un*Gm,n*cosdeltm-deltn+Bm,n*sindeltm-deltn; Nm,n=-Jm,n; Lm,n= Hm,n; Jacobi2*m-1,2*n-1=

9、Hm,n; Jacobi2*m-1,2*n =Nm,n; Jacobi2*m ,2*n-1=Jm,n; Jacobi2*m ,2*n =Lm,n; end end end %計(jì)算mn情況下的Jacobi矩陣中的子矩陣元素 Correction=-JacobiUnbalance' %計(jì)算電壓相角和幅值的修正量 precision = maxabsCorrection; %取誤差最大值 for m=1:Npq deltm=deltm+Correction2*m-1; %修正PQ節(jié)點(diǎn)電壓相角 um= um+Correction2*m ; %修正PQ節(jié)點(diǎn)電壓幅值 end k=k+1; %迭代輪數(shù)+1endk, u, delt, Jacobi, precisionfor m=1:Nbus Um=um*cosdeltm+j*sindeltm; %采用直角坐標(biāo)系表示電壓 Im=YNbus,m*Um; %計(jì)算注入平衡節(jié)點(diǎn)的電流endSslack=UNbus*sumconjI %計(jì)算注入平衡節(jié)點(diǎn)的功率for m=1:Nbus for n=1:Nbus Sm,n=Um*conjUm-conjUn*conj-Ym,n