【導(dǎo)讀】或集體已經(jīng)發(fā)表或撰寫的成果作品。對本文的研究作出重要貢獻的個人和集體，均已在文中以明確方式表明。本人完全意識到本申明的法律后果由本人承擔(dān)。本人授權(quán)南昌大學(xué)可以將本論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進行檢。索，可以采用影印、縮印或掃描等復(fù)制手段保存和匯編本學(xué)位論文。保密□，在年解密后適用本授權(quán)書。它根據(jù)給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀態(tài)。來定量地分析比較供電方案或運行方式的合理性，可靠性和經(jīng)濟性。使用方便，有著其他高級語言無法比擬的強大的矩陣處理功能。MATLAB成為電力系統(tǒng)潮流計算的首選計算機語言。法，最后介紹了利用MATLAB制作潮流計算軟件的過程。

　　

【正文】 。 % 計算電壓修正量 bus(1:nPQ,2) = bus(1:nPQ,2) dU。 % 修正電壓 bus(1:nb1,3) = bus(1:nb1,3) dAng。 % 修正相角 if (max(abs(dU))EPS)amp。(max(abs(dAng))EPS) break end % 判斷是否滿足 精度誤差，如滿足則跳出，否則返回繼續(xù)迭代 end bus = PQ_(bus,Y,nPQ,nPV)。 % 計算每個節(jié)點的有功和無功注入的子程序 [bus,line] = ReNum_(bus,line,nodenum)。 % 對節(jié)點恢復(fù)編號的子程序 YtYm = YtYm_(line)。 % 計算線路的等效 Yt 和 Ym 的子程序，以計算線路潮流 bus_res = bus_res_(bus)。 % 計算節(jié)點數(shù)據(jù)結(jié)果的子程序 S_res = S_res_(bus,line,YtYm)。 % 計算線路潮流的子程序 myf = fopen(39。39。,39。a39。)。 fprintf(myf,39。牛頓－拉夫遜法潮流計算結(jié)果 \n\n 節(jié)點計算結(jié)果： \n 節(jié) 25 點 節(jié)點電壓 節(jié)點相角（角度） 節(jié)點注入功率 \n39。)。 for i = 1:nb fprintf(myf,39。% 39。,bus_res(i,1))。 fprintf(myf,39。.6f 39。,bus_res(i,2))。 fprintf(myf,39。.6f 39。,bus_res(i,3))。 fprintf(myf,39。.6f + j .6f\n39。,real(bus_res(i,4)),imag(bus_res(i,4)))。 end fprintf(myf,39。\n 線路計算結(jié)果： \n 節(jié)點 I 節(jié)點 J 線路功率 S(I,J) 線路功率S(J,I) 線路損耗 dS(I,J)\n39。)。 for i = 1:nl fprintf(myf,39。% 39。,S_res(i,1))。 fprintf(myf,39。% 39。,S_res(i,2))。 fprintf(myf,39。.6f + j .6f 39。,real(S_res(i,3)),imag(S_res(i,3)))。 fprintf(myf,39。.6f + j .6f 39。,real(S_res(i,4)),imag(S_res(i,4)))。 fprintf(myf,39。.6f + j .6f\n39。,real(S_res(i,5)),imag(S_res(i,5)))。 end fclose(myf)。 % 迭代結(jié)束后繼續(xù)在 “”寫入節(jié)點計算結(jié)果和線路計算結(jié)果 程序結(jié)束 子程序 1 作 用為打開數(shù)據(jù)文件 function [bus,line] = OpDF_ [dfile,pathname]=uigetfile(39。*.m39。,39。Select Data File39。)。 % 數(shù)據(jù)文件類型為 m文件，窗口標(biāo)題為“Select Data File” if pathname == 0 error(39。 you must select a valid data file39。) % 如果沒有選擇有效文件，則出現(xiàn)錯誤提示 else lfile =length(dfile)。 % 讀取文件字符串長度 eval_r(dfile(1:1file2))。 % 去除后綴，打開文件 ！注意！新浪博客中 eval函數(shù)會自動加上 _r后綴，此處的函數(shù)名是 eval而不是 eval_r，拷貝后請去除 _r后綴 子程序 2 作用為對節(jié)點重排序，并修改相應(yīng)的線路數(shù)據(jù) function [bus,line,nPQ,nPV,nodenum] = Num_(bus,line) [nb,mb]=size(bus)。 [nl,ml]=size(line)。 nSW = 0。 % nSW 為平衡節(jié)點個數(shù) 26 nPV = 0。 % nPV為 PV節(jié)點個數(shù) nPQ = 0。 % nPQ 為 PQ節(jié)點個數(shù) for i = 1:nb, % nb 為總節(jié)點數(shù) type= bus(i,6)。 if type == 3, nSW = nSW + 1。 SW(nSW,:)=bus(i,:)。% 計算并儲存平衡節(jié)點 elseif type == 2, nPV = nPV +1。 PV(nPV,:)=bus(i,:)。% 計算并儲存 PV節(jié)點 else nPQ = nPQ + 1。 PQ(nPQ,:)=bus(i,:)。% 計算并儲 存 PQ 節(jié)點 end end bus=[PQ。PV。SW]。 % 對 bus 矩陣按 PQ、 PV、平衡節(jié)點的順序重新排序 nodenum=[[1:nb]39。 bus(:,1)]。% 生成新舊節(jié)點對照表 for i=1:nl for j=1:2 for k=1:nb if line(i,j)==nodenum(k,2) line(i,j)=nodenum(k,1)。 break end end end end % 按排序以后的節(jié)點順序?qū)?line 矩陣重新編號 子程序 3 作用為計算節(jié)點導(dǎo)納矩陣 function Y = y_(bus,line) [nb,mb]=size(bus)。 [nl,ml]=size(line)。 Y=zeros(nb,nb)。 % 對導(dǎo)納矩陣賦初值 0 for k=1:nl I=line(k,1)。 27 J=line(k,2)。 Zt=line(k,3)+j*line(k,4)。 % 讀入線路參數(shù) if Zt~=0 Yt=1/Zt。 % 接地支路不計算 Yt end Ym=line(k,5)+j*line(k,6)。 K=line(k,7)。 if (K==0)amp。(J~=0) % 普通線路 : K=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)。 end if (K==0)amp。(J==0) % 對地支路 : K=0,J=0,R=X=0 Y(I,I)=Y(I,I)+Ym。 end if K0 % 變壓器線路 : Zt 和 Ym 為折算到 i側(cè)的值 ,K 在 j 側(cè) Y(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)。 end if K0 % 變壓器線路 : Zt 和 Ym 為折算到 K 側(cè)的值 ,K 在 i側(cè) Y(I,I)=Y(I,I)+Yt+Ym。 Y(J,J)=Y(J,J)+K*K*Yt。 Y(I,J)=Y(I,J)+K*Yt。 Y(J,I)=Y(I,J)。 end end 子程序 4 作用為計算功率偏差 28 function [dP,dQ] =dPQ_(Y,bus,nPQ,nPV) % nPQ、 nPV為相應(yīng)節(jié)點個數(shù) n = nPQ + nPV +1。 % 總節(jié)點個數(shù) dP = bus(1:n1,4)。 dQ = bus(1:nPQ,5)。 % 對 dP 和 dQ賦初值 PV節(jié)點不需計算 dQ 平衡節(jié)點不參與計算 for i = 1:n1 for j = 1:n dP(i,1) = dP(i,1)bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)bus(j,3))+imag(Y(i,j))*sin(bus(i,3)bus(j,3)))。 if inPQ+1 dQ(i,1) = dQ(i,1)bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)bus(j,3))imag(Y(i,j))*cos(bus(i,3)bus(j,3)))。 end end end % 利用循環(huán)計算求取 dP 和 dQ 子程序 5 作用為計算雅克比矩陣 function J = Jac_(bus,Y,nPQ) [nb,mb]=size(bus)。 H = zeros(nb1,nb1)。 N = zeros(nb1,nPQ)。 K = zeros(nPQ,nb1)。 L = zeros(nPQ,nPQ)。 % 將雅克比矩陣分塊，即： J = [H N。 K L]，并初始化 Qi = zeros(nb1,1)。 Pi = zeros(nb1,1)。 for i = 1:nb1 for j = 1:nb Pi(i,1)=Pi(i,1)+bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)bus(j,3))+imag(Y(i,j))*sin(bus(i,3)bus(j,3)))。 Qi(i,1)=Qi(i,1)+bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)bus(j,3))imag(Y(i,j))*cos(bus(i,3)bus(j,3)))。 end end % 初始化并計算 Qi和 Pi 29 for i = 1:nb1 for j = 1:nb1 if i~=j H(i,j)=bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)bus(j,3))imag(Y(i,j))*cos(bus(i,3)bus(j,3)))。 else H(i,j)=Qi(i,1)+imag(Y(i,j))