【正文】 點(diǎn)處理if B2(i,6)==2j=j+1。DetaU。%九：求解修正方程，獲取節(jié)點(diǎn)電壓的不平衡量DetaU=zeros(2*n2,1)。endif k==(n1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0。Jacbi(2*h,2*k1)=0。else%非對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。if i==j%對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。for i=1:n%對(duì)PV節(jié)點(diǎn)的處理if B2(i,6)==3h=h+1。endif k==(n1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。else%非對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。if i==j%對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。for i=1:n%對(duì)PQ節(jié)點(diǎn)的處理if B2(i,6)==2h=h+1。h=0。end end I。for i=1:nif i~=isbh=h+1。%七：創(chuàng)建I，用于存儲(chǔ)節(jié)點(diǎn)電流參數(shù)i=zeros(n1,1)。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。t=t+1。for i=1:n%六：對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0if i~=isbamp。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。B2(i,6)==2h=h+1。%五：創(chuàng)建DetaS,用于存儲(chǔ)有功功率、無(wú)功功率和電壓幅值的不平衡量h=0。PVU(t,1)=B2(i,3)。t=0。endend end OrgS。for j=1:nOrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。endend end %三：對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0 for i=1:nif i~=isbamp。for j=1:nOrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*Imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。for i=1:n%對(duì)PQ節(jié)點(diǎn)的處理if i~=isbamp。%二：創(chuàng)建OrgS，用于存儲(chǔ)初始功率參數(shù)h=0。%將OrgS、DetaS初始化OrgS=zeros(2*n2,1)。Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3))。Y(q,p)=Y(p,q)。q=B1(i,2)。Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。Y(q,p)=Y(p,q)。q=B1(i,2)。Times=1。)。X=input(39。請(qǐng)輸入節(jié)點(diǎn)參數(shù):B2=39。)。B1=input(39。請(qǐng)輸入誤差精度:pr=39。)。isb=input(39。請(qǐng)輸入支路數(shù):n1=39。)。%清空內(nèi)存n=input(39。三、實(shí)驗(yàn)要求每人一組，在實(shí)驗(yàn)課時(shí)內(nèi)，用MATLAB調(diào)試和修改運(yùn)行程序，用算例計(jì)算輸出潮流結(jié)果。熟悉MATLAB，并掌握MATLAB程序的基本調(diào)試方法。disp(PPH)。平衡節(jié)點(diǎn)的功率39。for i=1:nPPH=PPH+EF(n)*conj(Y(n,i))*conj(EF(i))。,i)disp(angle(EF(i))*180/pi)。fprintf(39。%d號(hào)節(jié)點(diǎn)電壓的幅值為:39。EF=E+j*F。收斂后的節(jié)點(diǎn)電壓用極坐標(biāo)表示為:39。disp(FX)。第%d次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:n39。disp(E+j*F)。第%d次迭代后的節(jié)點(diǎn)電壓分別為:n39。end a=a+1。for i=1:n1E(i)=E(i)+X(i)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(n+i1,n+w1)=2*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,n+w1)=0。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,n+w1)=F1(n+i1,n+w1)G(i,k)*E(k)+B(i,k)*F(k)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(n+i1,n+w1)=G(i,i)*E(i)+B(i,i)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,n+w1)=G(i,w)*E(i)+B(i,w)*F(i)。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。%△Qfor w=1:n FX(i)= FX(i)E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)F(i)*G(i,w)*F(w)F(i)*B(i,w)*E(w)。%△Qend end for i=m+1:n1 %PV節(jié)點(diǎn) FX(i)=PP(im)。%△Qfor w=1:n FX(i)= FX(i)E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)F(i)*G(i,w)*F(w)F(i)*B(i,w)*E(w)。%最后的節(jié)點(diǎn)電壓矩陣 while max(FX)=10^(5)for i=1:m %PQ節(jié)點(diǎn)FX(i)=P(i)。%F(x)導(dǎo)數(shù)矩陣 a=0。%總的節(jié)點(diǎn)數(shù)FX=ones(2*n2,1)。X=[]。000]。 0。B=[ 00。0000 0。 0 0。%E=PQ,PV,Vθ節(jié)點(diǎn)e的初值F=[0 0 0 0 0]39。%PP=PV節(jié)點(diǎn)的P值 V=[]。%P=PQ節(jié)點(diǎn)的P值 Q=[ 0]。%m=PQ節(jié)點(diǎn)個(gè)數(shù) v=1。3150=180。8760+252180。8760+DPT180。180。1102Mvar=2等值電路中流過(guò)線路等值阻抗的功率為S1=S+DST+DS0+jQB2=40+j30++++=++180。106180。248。248。31500kVA 247。=2231。I%=2231。kVA247。31500247。248。247。232。100232。231。SDST=DPT+jDQT=2231。246。線路和變壓器的參數(shù)如下：線路（每回）：r=, x=, b=180。3．6 如圖所示網(wǎng)絡(luò)，變電所低壓母線上的最大負(fù)荷為40MW，cosj=，Tmax=4500h。=ST39。j1MVA=8+ 21039。T182+=8++180。VA=，VA=109kV（5）計(jì)及從高壓母線輸入變壓器T1和T2的功率。10248。k2=(10+)180。231。1246。10248。k1=(10+)180。231。1246。ST1=S1LD+Sc=8+j6+=8+ST2=S2LDSc=8+=8+（4）計(jì)算高壓側(cè)電壓。VBDE10180。k1248。10248。247。=101247。若取其假定正方向?yàn)轫槙r(shí)針?lè)较颍瑒t可得k2 DE187。因兩臺(tái)變壓器電抗相等，故S1LD=S2LD=11SLD=(16+j12)MVA=8+j6MVA 22（2）求循環(huán)功率。已知低壓母線電壓10kV，負(fù)荷功率為16+j12MVA，試求變壓器的功率分布和高壓側(cè)電壓。+180。為了計(jì)算線路Ⅰ的電壓損耗，要用A點(diǎn)的電壓和功率SA1。（4）計(jì)算電壓損耗。取SI續(xù)進(jìn)行計(jì)算。DSTb242+182(+)MVA=+ =2110Sb=SLDb+DSTb+DSob+jDQBI+jDQBIII=24+j18++++=++92(+)MVA=+ =1102Sc=SLDc+DSTc+DSoc+jDQBIII+jDQBII=12+j9++++=+（3）計(jì)算閉式網(wǎng)絡(luò)的功率分布。104180。40S=180。40W=18+=180。104180。104S2DQB213。106180。50W=+B213。1102Mvar=線路Ⅱ： Z213。104S2DQBI=180。106180。線路Ⅰ： ZI=(+)180。網(wǎng)絡(luò)各組件參數(shù)為：變電所b SNXT==20MVA，DS0=+，RT=，變電所c SNXT=127W=10MVA，DS0=+，RT=，負(fù)荷功率 SLDb=24+j18MVA，SLDc=12+j9MVA試求電力網(wǎng)絡(luò)的功率分布及最大電壓損耗。+180。+180。進(jìn)一步可求得E點(diǎn)電壓 +=(+)MVA=+10239。10+j1200180。10+2200180。7+j1200180。7+2200180。()Z229。(+)=+ 求C點(diǎn)和D點(diǎn)的運(yùn)算負(fù)荷，為 DSCE +=(+)=+ 210=2600+j1600+300+j160++=+SD=600+j200+1600+j1000=2200+j1200kVA循環(huán)功率Sc(V=*A*VBVN()180。(+)=+ZL3=3180。（線路參數(shù)LJ185：z=+；LJ70：z=+）解 線路等值阻抗ZL1=10180。kV時(shí)的初始如圖所示。kV、Vamp。線路LL2和L3導(dǎo)線型號(hào)均為L(zhǎng)J185，線路長(zhǎng)度分別為10km，4km和3km，線路L4為2km長(zhǎng)的LJ70導(dǎo)線；各負(fù)荷點(diǎn)負(fù)荷amp。=，Vc= 與計(jì)及電壓降落橫分量的計(jì)算結(jié)果相比，