資源描述：
《《matlab潮流計算》word版》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在應用文檔-天天文庫。
1、%線路數(shù)據(jù)Line=[1-節(jié)點編號2-線路首端節(jié)點號3-線路末端節(jié)點號4-支路電阻5-支路電抗6-支路電納(注意：此處取的是B/2)]Line=[1120.040.250.252130.10.3503230.080.300.25];%變壓器數(shù)據(jù)transform=[1-支路編號2-支路首節(jié)點編號3-支路末節(jié)點編號4-支路電阻(p.u.)5-支路電抗(p.u.)6-變壓器變比(p.u.)]transform=[12400.0151.0523500.031.05];%數(shù)據(jù)預處理Nbus=5;%節(jié)點數(shù)nline=size(Line,1);%線路個數(shù)ntrans=size(transform,1);%
2、變壓器個數(shù)slack=5;%平衡節(jié)點號Npq=4;%PQ節(jié)點的個數(shù)Ss=0;%計算節(jié)點導納矩陣Y=zeros(Nbus);ifnline>=1%判斷是否存在線路fork=1:nline%以下處理線路t1=Line(k,2);t2=Line(k,3);b2=Line(k,6);%分別取出線路的首端節(jié)點編號t1、末端節(jié)點編號t2和對地電納b2Yl=1/(Line(k,4)+j*Line(k,5));%計算線路的支路電導YlY(t1,t1)=Y(t1,t1)+Yl+j*b2;%修正第k條線路首端節(jié)點的自導納Y(t1,t2)=Y(t1,t2)-Yl;%修正第k條線路首端節(jié)點與末端節(jié)點之間的互導納Y(t
3、2,t1)=Y(t2,t1)-Yl;%修正第k條線路末端節(jié)點與首端節(jié)點之間的互導納Y(t2,t2)=Y(t2,t2)+Yl+j*b2;%修正第k條線路末端節(jié)點的自導納endendifntrans>=1%判斷是否存在變壓器fork=1:ntrans%以下處理變壓器t1=transform(k,2);t2=transform(k,3);t3=transform(k,6);%分別取出變壓器的首端節(jié)點編號t1、末端節(jié)點編號t2和變比t3Yt=1/(transform(k,4)+j*transform(k,5));Yt1=Yt/t3;Yt2=Yt*(1-t3)/(t3*t3);Yt3=Yt*(t3-1)
4、/t3;Y(t1,t1)=Y(t1,t1)+Yt1+Yt2;Y(t1,t2)=Y(t1,t2)-Yt1;Y(t2,t1)=Y(t2,t1)-Yt1;Y(t2,t2)=Y(t2,t2)+Yt1+Yt3;endendG=real(Y);B=imag(Y);%區(qū)分節(jié)點導納矩陣的實部和虛部GB%賦初值delt(1)=0;delt(2)=0;delt(3)=0;delt(4)=0;u(1)=1;u(2)=1;u(3)=1;u(4)=1;p(1)=-0.30;q(1)=-0.18;p(2)=-0.55;q(2)=-0.13;p(3)=0;q(3)=0;p(4)=0.5;q(4)=1.10;p(5)=0.8
5、;q(5)=0.50;k=0;precision=1;Npq=4;%Npq分別是網(wǎng)絡中的PQ節(jié)點數(shù)%[Unbalance]=-[Jacobi][Correction]whileprecision>0.00001%設定誤差上限，判斷是否繼續(xù)迭代u(5)=1.06;delt(5)=0;%設定平衡節(jié)點電壓相角與幅值k;u;delt;form=1:Npqforn=1:Nbuspt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));%由節(jié)點電壓求得的PQ節(jié)點注入有功功率qt(n)=u(m)*u(n)*(G(m,n
6、)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));%由節(jié)點電壓求得的PQ節(jié)點注入無功功率endUnbalance(2*m-1)=p(m)-sum(pt);%計算PQ節(jié)點有功功率不平衡量Unbalance(2*m)=q(m)-sum(qt);%計算PQ節(jié)點無功功率不平衡量end%[Unbalance]是節(jié)點不平衡量矩陣form=1:Npqforn=1:Nbush0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));n0(n)=-u(m)*u(n)*(G
7、(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));l0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt
