【正文】 差值 l=length(e)。for i=1:(l1)s1=0。for j=1:ls1=s1+g(i,j)*e(j)b(i,j)*f(j)。ends(2*i1)=kind(2,i)e(i)*s1f(i)*s2。elses(2*i)=kind(3,i)^2f(i)^2e(i)^2。J=zeros(2*l2,2*l2)。if kind(1,i)==1s=PQ_out(g,b,e,f,i)。J(2*i,j)=s(2,j)。for j=1:(2*l2)J(2*i1,j)=s(1,j)。endend endfunction pq=PQ_out(g,b,e,f,i)%計算pq節(jié)點的雅克比矩陣 l=length(e)。for j=1:(l1)if j==i s=0。endpq(1,2*i1)=sg(i,i)*e(i)b(i,i)*f(i)。for k=1:ls=s(g(i,k)*f(k)+b(i,k)*e(k))。s=0。endpq(2,2*i1)=s+b(i,i)*e(i)g(i,i)*f(i)。for k=1:ls=s(g(i,k)*e(k)b(i,k)*f(k))。elsepq(1,2*j1)=(g(i,j)*e(i)+b(i,j)*f(i))。pq(2,2*j)=pq(1,2*j1)。end endfunction pv=PV_out(g,b,e,f,i)%計算pv節(jié)點的雅克比矩陣 l=length(e)。for j=1:(l1)if j==i s=0。endpv(1,2*i1)=sg(i,i)*e(i)b(i,i)*f(i)。for k=1:ls=s(g(i,k)*f(k)+b(i,k)*e(k))。pv(2,2*i1)=2*e(i)。elsepv(1,2*j1)=(g(i,j)*e(i)+b(i,j)*f(i))。end end%數(shù)據(jù)輸入g=[ 0 0 0 0 0 ]。e=[1 1 ]。kind=[1 1 2 0]。e f四、例題及運行結(jié)果在上圖所示的簡單電力系統(tǒng)中，系統(tǒng)中節(jié)點2為PQ節(jié)點，節(jié)點3為PV節(jié)點，節(jié)點4為平衡節(jié)點，已給定 P1s+jQ1s= P2s+jQ2s= P3s= V3s= V4s=∠0176。答：高斯迭代法（高斯塞德爾法），牛頓拉夫遜法以及PQ分解法。？為什么？答:不一樣,因為每次迭代的電壓、有功、無功都是與前一次不同的新值，所以每次迭代過程中，雅可比矩陣都是變化的。我這用了牛頓法直角坐標(biāo)系來編寫程序的。第五篇：電力系統(tǒng)分析計算機算法PSDBPA潮流計算實驗報告電力系統(tǒng)分析的計算機算法實驗報告學(xué)生姓名 課 程 電力系統(tǒng)分析的計算機算法 學(xué) 號專 業(yè) 電氣工程及其自動化 指導(dǎo)教師 邱曉燕二Ο一四 年 六 月 二日實驗一潮流計算一、實驗?zāi)康?．了解并掌握電力系統(tǒng)計算機算法的相關(guān)原理。3．了解并掌握PSDBPA電力系統(tǒng)分析程序單線圖和地理接線圖的使用。無論是電力系統(tǒng)規(guī)劃、設(shè)計還是運行，對其安全穩(wěn)定進(jìn)行分析都是極其重要的。本實驗課程基于PSDBPA平臺，結(jié)合《電力系統(tǒng)分析計算機算法》課程，旨在引導(dǎo)學(xué)生將理論知識和實際工程相結(jié)合，掌握電力系統(tǒng)穩(wěn)態(tài)、暫態(tài)分析的原理、分析步驟以及結(jié)論分析。三、原理和說明PSDBPA電力系統(tǒng)分析程序穩(wěn)態(tài)分析主要是潮流計算，軟件中潮流程序的計算方法有P_Q分解法，牛頓_拉夫遜法，改進(jìn)的牛頓－拉夫遜算法。注：采用PQ分解法和牛頓－拉夫遜法相結(jié)合，以提高潮流計算的收斂性能，程序通常先采用PQ分解法進(jìn)行初始迭代，然后再轉(zhuǎn)入牛頓－拉夫遜法求解潮流。除了潮流計算功能外，該軟件還具有自動電壓控制、聯(lián)絡(luò)線功率控制、系統(tǒng)事故分析（N1開斷模擬）、網(wǎng)絡(luò)等值、靈敏度分析、節(jié)點P－V、Q－V和P－Q曲線、確定系統(tǒng)極限輸送水平、負(fù)荷靜特性模型、靈活多樣的分析報告、詳細(xì)的檢錯功能等功能。不正常時，將留在硬盤上，可隨時刪除） 儲存一個潮流作業(yè)計算時屏幕顯示的信息。該文件安裝時放在與潮流程序相同的目錄中。常用的控制語句主要包括：(1)指定潮流文件開始的一級控制語句“(POWERFLOW, CASEID=方式名, PROJECT=工程名)”(2)指定計算方法和最大迭代次數(shù)的控制語句“/SOL_ITER, DECOUPLED=PQ法次數(shù), NEWTON=牛拉法次數(shù)”；(3)指定計算結(jié)果輸出的控制語句“/P_OUTPUT_LIST, ?”；(4)指定計算結(jié)果輸出順序的控制語句“/RPT_SORT= ?”；(5)指定計算結(jié)果分析列表的控制語句“/P_ANALYSIS, LEVEL= ？”；(6)指定潮流結(jié)果二進(jìn)制文件名的控制語句“/NEW_BASE, FILE = 文件名”；(7)指定潮流圖和地理接線圖使用的結(jié)果文件控制語句“/PF_MAP，F(xiàn)ILE=文件名”；(8)指定網(wǎng)絡(luò)數(shù)據(jù)的控制語句“/NETWORK_DATA”；(9)指定潮流數(shù)據(jù)文件結(jié)束的控制語句“(END)”； （PFO文件）潮流計算結(jié)果文件內(nèi)容主要分下述幾個方面： 1）程序控制語句列表。3）錯誤信息。4）誤差控制參數(shù)列表。6）計算結(jié)果輸出：詳細(xì)計算結(jié)果列表：按節(jié)點、與該節(jié)點相聯(lián)接支路順序，并根據(jù)用戶的要求（通過控制語句控制）可按照字母、分區(qū)或區(qū)域排序輸出潮流計算結(jié)果。7）錯誤信息統(tǒng)計。表1 IEEE 9節(jié)點算例節(jié)點參數(shù)表2 IEEE 9節(jié)點算例線路參數(shù)表3 IEEE 9節(jié)點算例變壓器參數(shù)注：表1表3中功率基準(zhǔn)值為100MVA；電阻、電感值為標(biāo)幺值。例1：(POWERFLOW,CASEID=IEEE9,PROJECT=IEEE_9BUS_TEST_SYSTEM)/SOL_ITER,DECOUPLED=2,NEWTON=15,OPITM=0./P_INPUT_LIST,ZONES=ALL /P_OUTPUT_LIST,ZONES=ALL /RPT_SORT=ZONE /NEW_BASE,FILE= /PF_MAP,FILE = /NETWORK_DATA BS GEN1 BGEN1BSTATIONA STATIONB STATIONC BGEN2BE GEN2 10 25 BGEN3 BE GEN3 .Ltransmission linesLGEN1 .. LGEN1 L.. L.. L.. L.. L..T GEN1230..0576 GEN2230..0625 GEN3230..0586 230.(END)四、實驗過程及結(jié)果（一）IEEE9節(jié)點算例： ：，并進(jìn)行計算，結(jié)果如下： 1）系統(tǒng)數(shù)據(jù) 2）計算過程迭代信息及詳細(xì)的輸出列表：小結(jié)總結(jié) * 并聯(lián)無功補償數(shù)據(jù)列表/電容器（Mvar）//電抗器（Mvar）/區(qū)域/分區(qū)最大容量使用容量備用未安排容量最大容量使用容量備用未安排容量01總結(jié)TRANSMISSION LINES CONTAINING COMPENSATIONOWN ZONE BUS1BASE1 ZONE BUS2BASE2ID PERCENTCASE CONTAINS NO TRANSMISSION LINES WITH SERIES COMPENSATION* 節(jié)點相關(guān)數(shù)據(jù)列表節(jié)點電壓/發(fā)電/ /負(fù)荷//無功補償/ 類型 擁有者 分區(qū)電壓/角度kVMWMVAR 功率因數(shù)MWMVAR使用的存在的未安排PU/度發(fā)電機1 S01發(fā)電機2 E01發(fā)電機3 E01母線101母線201母線301母線A01母線B01母線C01整個系統(tǒng)電容器總和電抗器總和 * 旋轉(zhuǎn)備用數(shù)據(jù)列表有功功率無功功率區(qū)域/分區(qū)最大值實際出力備用最大值最小值已發(fā)無功吸收無功備用(MW)(MW)(MW)(MVAR)(MVAR)(MVAR)(MVAR)(MVAR)01總結(jié)說明：（如 抽水蓄能電機）。當(dāng)最大出力值小于實際出力時，統(tǒng)計時最大出力值用實際出力值代替。無功旋轉(zhuǎn)備用只統(tǒng)計有功出力大于0并且基準(zhǔn)電壓小于30kV的發(fā)電機。牛頓拉夫遜法迭代次數(shù)為 5次。T母線1母線3 (END)五、總結(jié)及思考題實驗中遇到的問題及解決方法：路徑錯誤——————重設(shè)各個參數(shù)路徑 卡片無法識別—————將參數(shù)規(guī)范化本次實驗使我初步掌握了PSDBPA軟件在電力系統(tǒng)潮流計算中的使用方法，收獲良多，為今后的工作打下了基礎(chǔ)。電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種基本計算。電力系統(tǒng)潮流計算的結(jié)果是電力系統(tǒng)穩(wěn)定計算和故障分析的基礎(chǔ)。此后，潮流計算曾采用了各種不同的方法，這些方法的發(fā)展主要是圍繞著對潮流計算的一些基本要求進(jìn)行的。牛頓拉夫遜方法是解非線性代數(shù)方程組的一種基本方法，在潮流計算中也得到