【正文】 1)在電網(wǎng)規(guī)劃階段,通過潮流計算,合理規(guī)劃電源容量及接入點,合理規(guī)劃網(wǎng)架,選擇無功補償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。 (2)在編制年運行方式時,在預(yù)計負荷增長及新設(shè)備投運基礎(chǔ)上,選擇典型方式進行潮流計算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對規(guī)劃、基建部門提出改進網(wǎng)架結(jié)構(gòu),加快基建進度的建議。 (3)正常檢修及特殊運行方式下的潮流計算,用于日運行方式的編制,指導(dǎo)發(fā)電廠開機方式,有功、無功調(diào)整方案及負荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。 　　(4)預(yù)想事故、設(shè)備退出運行對靜態(tài)安全的影響分析及作出預(yù)想的運行方式調(diào)整方案。 總結(jié)為在電力系統(tǒng)運行方式和規(guī)劃方案的研究中，都需要進行潮流計算以比較運行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟性。同時，為了實時監(jiān)控電力系統(tǒng)的運行狀態(tài)，也需要進行大量而快速的潮流計算。因此，潮流計算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運算。在系統(tǒng)規(guī)劃設(shè)計和安排系統(tǒng)的運行方式時，采用離線潮流計算；在電力系統(tǒng)運行狀態(tài)的實時監(jiān)控中，則采用在線潮流計算。 —拉夫遜法潮流計算基本步驟牛頓—拉夫遜法潮流計算的基本步驟不外乎有如下幾步。 形成節(jié)點導(dǎo)納矩陣。 設(shè)各節(jié)點電壓的初始值 和 求出功率偏差根據(jù)公式算出雅可比矩陣。解如下線性方程組，得到依據(jù)如下公式計算修正量，并形成迭代公式：這樣反復(fù)迭代計算，直至所有節(jié)點|ΔU|ε和|Δδ|ε為止。計算線路功率分布，PV節(jié)點無功功率和平衡節(jié)點注入功率。第四章程序流程圖及程序代碼根據(jù)牛頓—拉夫遜法潮流計算的基本步驟可以得到如下程序流程圖%*******************主程序開始********************clc。clear allformat long。%********************定義全局變量*****************global nSW nPQ nPV。 %平衡節(jié)點，PQ節(jié)點，PV節(jié)點個數(shù)global nb。 %節(jié)點個數(shù)global nl。 %支路數(shù) global bus。 %bus:節(jié)點數(shù)據(jù)global line。 %line:支路數(shù)據(jù)global Y。 %Y:節(jié)點導(dǎo)納矩陣global nodenum。 %節(jié)點編號對應(yīng)表global lPQ。 %線路損耗矩陣global myf。 %輸出文件句柄%*******************子程序調(diào)用***************openfile。 %打開bus line數(shù)據(jù)文件 change。 %節(jié)點重新編號 ybus。 %建立節(jié)點導(dǎo)納矩陣 NR。 %牛頓拉夫遜法解非線性方程序組 PQ。 %計算節(jié)點注入功率 flow。 %計算線路功率及損耗 ret。 %恢復(fù)原節(jié)點編號 output。 %計算結(jié)果輸出%*************************主程序結(jié)束*****************%******************打開bus line的數(shù)據(jù)文件*****************function openfile global nb nl。global line。global bus。global myf。 [dfile,pathname]=uigetfile(39。*.m39。,39。Select Data File39。)。 %打開數(shù)據(jù)文件if pathname == 0 error(39。 you must select a valid data file39。)else lfile =length(dfile)。 % strip off .m eval(dfile(1:lfile2))。 %執(zhí)行數(shù)據(jù)文件中語句end [nl,ml]=size(line)。 %求線路數(shù)nl[nb,mb]=size(bus)。 %求節(jié)點數(shù)nb outfile=39。39。for I=1:length(dfile) if dfile(I)==39。.39。 break。 end outfile=strcat(outfile,dfile(I))。 endoutfile=strcat(outfile,39。.txt39。)。 %輸出文件名形成myf=fopen(outfile,39。w39。)。 %打開輸出文件%**********************給節(jié)點編號*********************function changeglobal nb。global nl。global nPQ。global bus。global line。global nodenum。 nSW = 0。 nPV = 0。 nPQ = 0。 for I = 1:nb, % nb為總節(jié)點數(shù) type= bus(I,6)。 if type == 3, %求平衡節(jié)點數(shù) nSW = nSW + 1。 SW(nSW,:)=bus(I,:)。 elseif type == 2, %求PV節(jié)點數(shù) nPV = nPV +1。 PV(nPV,:)=bus(I,:)。else %求PQ節(jié)點數(shù) nPQ = nPQ + 1。 PQ(nPQ,:)=bus(I,:)。 endend bus=[PQ。PV。SW]。 %按PQ PV SW順序排列的節(jié)點數(shù)據(jù)矩陣newbus=[1:nb]39。nodenum=[newbus bus(:,1)]。 %形成節(jié)點號對應(yīng)表bus(:,1)=newbus。 for I=1:nl %根據(jù)節(jié)點號對應(yīng)表重新排列線路矩陣節(jié)點編號 for J=1:2 for k=1:nb if line(I,J)==nodenum(k,2) line(I,J)=nodenum(k,1)。 break end end endend %******************建立節(jié)點導(dǎo)納矩陣*****************function Y = ybus(bus,line)global nl。global nb。global bus。global line。global Y。global myf。 Y=zeros(nb,nb)。 %對導(dǎo)納矩陣賦初值0 for k=1:nl I=line(k,1)。 %讀入線路參數(shù) J=line(k,2)。 Zt=line(k,3)+j*line(k,4)。 if J~=0 %接地支路分母為零 Yt=1/Zt。 %非接地支路為阻抗的倒數(shù) end Ym=line(k,5)+j*line(k,6)。 %求Ym K=line(k,7)。 %求K 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為