【正文】 )。應(yīng)用計算例題驗證程序的計算效果。程序要求根據(jù)已知的電力網(wǎng)的數(shù)學模型（節(jié)點導納矩陣）及各節(jié)點參數(shù)，完成該電力系統(tǒng)的潮流計算，要求計算出節(jié)點電壓、功率等參數(shù)。五、實驗流程六、實驗結(jié)果參數(shù)輸入：運行結(jié)果：七、實驗體會通過這次實驗，讓我第一次接觸到了MATLAB,并深切體會到了它的強大之處；潮流計算的計算機算法的實現(xiàn)不僅鞏固了我的學過的知識，還讓我學到一些MATLAB的編程，雖然在實驗的過程中出現(xiàn)了很多的錯誤，但在老師的細心指導下，問題都解決啦；計算機為我們省去了大量的人工計算，希望在以后的學習中能接觸到更多的軟件，學習到更多的知識。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。%修正節(jié)點電壓 j=0。endendendend end Jacbi。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。endendendend end k=0。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。k=0。I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。end end DetaS。B2(i,6)==3h=h+1。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。endend end OrgS。endend end for i=1:nif i~=isbamp。for i=1:nif i~=isbamp。end end B2。for i=1:n%對PQ節(jié)點處理if B2(i,6)==2j=j+1。%九：求解修正方程，獲取節(jié)點電壓的不平衡量DetaU=zeros(2*n2,1)。Jacbi(2*h,2*k1)=0。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。for i=1:n%對PV節(jié)點的處理if B2(i,6)==3h=h+1。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。for i=1:n%對PQ節(jié)點的處理if B2(i,6)==2h=h+1。end end I。%七：創(chuàng)建I，用于存儲節(jié)點電流參數(shù)i=zeros(n1,1)。t=t+1。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。%五：創(chuàng)建DetaS,用于存儲有功功率、無功功率和電壓幅值的不平衡量h=0。t=0。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 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)))。%二：創(chuàng)建OrgS，用于存儲初始功率參數(shù)h=0。Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3))。q=B1(i,2)。Y(q,p)=Y(p,q)。Times=1。X=input(39。)。請輸入誤差精度:pr=39。isb=input(39。)。三、實驗要求每人一組，在實驗課時內(nèi)，用MATLAB調(diào)試和修改運行程序，用算例計算輸出潮流結(jié)果。3150=180。8760+DPT180。1102Mvar=2等值電路中流過線路等值阻抗的功率為S1=S+DST+DS0+jQB2=40+j30++++=++180。248。31500kVA 247。I%=2231。31500247。247。100232。SDST=DPT+jDQT=2231。線路和變壓器的參數(shù)如下：線路（每回）：r=, x=, b=180。=ST39。T182+=8++180。10248。231。10248。231。ST1=S1LD+Sc=8+j6+=8+ST2=S2LDSc=8+=8+（4）計算高壓側(cè)電壓。k1248。247。若取其假定正方向為順時針方向，則可得k2 DE187。已知低壓母線電壓10kV，負荷功率為16+j12MVA，試求變壓器的功率分布和高壓側(cè)電壓。為了計算線路Ⅰ的電壓損耗，要用A點的電壓和功率SA1。取SI續(xù)進行計算。104180。40W=18+=180。104S2DQB213。50W=+B213。104S2DQBI=180。線路Ⅰ： ZI=(+)180。+180。進一步可求得E點電壓 +=(+)MVA=+10239。10+2200180。7+2200180。(+)=+ 求C點和D點的運算負荷，為 DSCE +=(+)=+ 210=2600+j1600+300+j160++=+SD=600+j200+1600+j1000=2200+j1200kVA循環(huán)功率Sc(V=*A*VBVN()180。（線路參數(shù)LJ185：z=+；LJ70：z=+）解 線路等值阻抗ZL1=10180。kV、Vamp。=，Vc= 與計及電壓降落橫分量的計算結(jié)果相比，誤差很小。= (VbDVT)+(dVT)=22()2+=1111Vc=Vc180。Pc39。RT+Qc39。XLQ139。RL+Q139。39。+=(+)MVA=+39。V247。+Sb=+++=+線路中的功率損耗為230。=+Sc39。(+)MVA2231。變壓器繞組中的功率損耗為230。180。W= 2221000SN21000180。104180。180。= 21 XL=180。當變壓器取主軸時，求母線c的電壓。 如圖所示簡單系統(tǒng)，額定電壓為110KV 雙回輸電線路，長度為80km，采用LGJ150導線，其單位長度的參數(shù)為：r=，x=,b=180。V2DV24=(P23R23+Q23X23)DV23==187。(P12R12+Q12X12)180。=S23+S24+S2=+S12又39。運行結(jié)果：運行結(jié)果復(fù)制如下：第1次迭代后的節(jié)點電壓分別為:第1次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:00第2次迭代后的節(jié)點電壓分別為:第2次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第3次迭代后的節(jié)點電壓分別為:第3次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第4次迭代后的節(jié)點電壓分別為:第4次迭代后功率偏差△P △Q電壓偏差△V的平方分別為: *收斂后的節(jié)點電壓用極坐標表示為: 1號節(jié)點電壓的幅值為:2號節(jié)點電壓的幅值為:3號節(jié)點電壓的幅值為:4號節(jié)點電壓的幅值為:平衡節(jié)點的功率 第二篇：電力系統(tǒng)分析潮流計算例題3．1 電網(wǎng)結(jié)構(gòu)如圖3—11所示，其額定電壓為10KV。end fprintf(39。%d號節(jié)點電壓的相角度數(shù)為39。for i=1:n1 fprintf(39。enddisp(39。fprintf(39。fprintf(39。endendend end X=inv(F1)*(FX)。F1(n+i1,w)=2*E(i)。F1(n+i1,w)=0。F1(n+i1,w)=F1(n+i1,w)+G(i,k)*F(k)+B(i,k)*E(k)。F1(n+i1,w)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,w)=B(i,w)*E(i)G(i,w)*F(i)。%△P FX(n+i1)=V(im)^2E(i)^2F(i)^2。%△P FX(n+i1)=Q(i)。%F(x)矩陣F1=zeros(n1,n1)。Y=G+j*B。 15 0。 0 0。%V=PV節(jié)點的U值E=[1 1 1 ]39。%v=PV節(jié)點個數(shù)P=[ 0]。%m=PQ節(jié)點個數(shù) v=1。%PP=PV節(jié)點的P值 V=[]。 0 0。B=[ 00。000]。%總的節(jié)點數(shù)FX=ones(2*n2,1)。%最后的節(jié)點電壓矩陣 while max(FX)=10^(5)for i=1:m %PQ節(jié)點FX(i)=P(i)。%△Qend end for i=m+1:n1 %PV節(jié)點 FX(i)=PP(im)。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。end a=a+1。disp(E+j*F)。disp(FX)。EF=E+j*F。fprintf(39。for i=1:nPPH=PPH+EF(n)*conj(Y(n,i))*conj(EF(i))。disp(PPH)。P32++=(R23+jX23)=(1+j2)=+++=(R24+jX24)=(+j3)=+則： S23=S3+DS23=+S24=S4+DS24=+39。故： S12=S12+DS12=+（2）再用已知的線路始端電壓V1=12，求出線路各點電壓。V1DV12=(P24R24+Q24X24)DV24==187。=