【正文】 牛頓拉夫遜方法是解非線(xiàn)性代數(shù)方程組的一種基本方法，在潮流計(jì)算中也得到應(yīng)用。電力系統(tǒng)潮流計(jì)算的結(jié)果是電力系統(tǒng)穩(wěn)定計(jì)算和故障分析的基礎(chǔ)。T母線(xiàn)1母線(xiàn)3 (END)五、總結(jié)及思考題實(shí)驗(yàn)中遇到的問(wèn)題及解決方法：路徑錯(cuò)誤——————重設(shè)各個(gè)參數(shù)路徑 卡片無(wú)法識(shí)別—————將參數(shù)規(guī)范化本次實(shí)驗(yàn)使我初步掌握了PSDBPA軟件在電力系統(tǒng)潮流計(jì)算中的使用方法，收獲良多，為今后的工作打下了基礎(chǔ)。無(wú)功旋轉(zhuǎn)備用只統(tǒng)計(jì)有功出力大于0并且基準(zhǔn)電壓小于30kV的發(fā)電機(jī)。例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)四、實(shí)驗(yàn)過(guò)程及結(jié)果（一）IEEE9節(jié)點(diǎn)算例： ：，并進(jìn)行計(jì)算，結(jié)果如下： 1）系統(tǒng)數(shù)據(jù) 2）計(jì)算過(guò)程迭代信息及詳細(xì)的輸出列表：小結(jié)總結(jié) * 并聯(lián)無(wú)功補(bǔ)償數(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é)點(diǎn)相關(guān)數(shù)據(jù)列表節(jié)點(diǎn)電壓/發(fā)電/ /負(fù)荷//無(wú)功補(bǔ)償/ 類(lèi)型 擁有者 分區(qū)電壓/角度kVMWMVAR 功率因數(shù)MWMVAR使用的存在的未安排PU/度發(fā)電機(jī)1 S01發(fā)電機(jī)2 E01發(fā)電機(jī)3 E01母線(xiàn)101母線(xiàn)201母線(xiàn)301母線(xiàn)A01母線(xiàn)B01母線(xiàn)C01整個(gè)系統(tǒng)電容器總和電抗器總和 * 旋轉(zhuǎn)備用數(shù)據(jù)列表有功功率無(wú)功功率區(qū)域/分區(qū)最大值實(shí)際出力備用最大值最小值已發(fā)無(wú)功吸收無(wú)功備用(MW)(MW)(MW)(MVAR)(MVAR)(MVAR)(MVAR)(MVAR)01總結(jié)說(shuō)明：（如 抽水蓄能電機(jī)）。7）錯(cuò)誤信息統(tǒng)計(jì)。4）誤差控制參數(shù)列表。常用的控制語(yǔ)句主要包括：(1)指定潮流文件開(kāi)始的一級(jí)控制語(yǔ)句“(POWERFLOW, CASEID=方式名, PROJECT=工程名)”(2)指定計(jì)算方法和最大迭代次數(shù)的控制語(yǔ)句“/SOL_ITER, DECOUPLED=PQ法次數(shù), NEWTON=牛拉法次數(shù)”；(3)指定計(jì)算結(jié)果輸出的控制語(yǔ)句“/P_OUTPUT_LIST, ?”；(4)指定計(jì)算結(jié)果輸出順序的控制語(yǔ)句“/RPT_SORT= ?”；(5)指定計(jì)算結(jié)果分析列表的控制語(yǔ)句“/P_ANALYSIS, LEVEL= ？”；(6)指定潮流結(jié)果二進(jìn)制文件名的控制語(yǔ)句“/NEW_BASE, FILE = 文件名”；(7)指定潮流圖和地理接線(xiàn)圖使用的結(jié)果文件控制語(yǔ)句“/PF_MAP，F(xiàn)ILE=文件名”；(8)指定網(wǎng)絡(luò)數(shù)據(jù)的控制語(yǔ)句“/NETWORK_DATA”；(9)指定潮流數(shù)據(jù)文件結(jié)束的控制語(yǔ)句“(END)”； （PFO文件）潮流計(jì)算結(jié)果文件內(nèi)容主要分下述幾個(gè)方面： 1）程序控制語(yǔ)句列表。不正常時(shí)，將留在硬盤(pán)上，可隨時(shí)刪除） 儲(chǔ)存一個(gè)潮流作業(yè)計(jì)算時(shí)屏幕顯示的信息。注：采用PQ分解法和牛頓－拉夫遜法相結(jié)合，以提高潮流計(jì)算的收斂性能，程序通常先采用PQ分解法進(jìn)行初始迭代，然后再轉(zhuǎn)入牛頓－拉夫遜法求解潮流。本實(shí)驗(yàn)課程基于PSDBPA平臺(tái)，結(jié)合《電力系統(tǒng)分析計(jì)算機(jī)算法》課程，旨在引導(dǎo)學(xué)生將理論知識(shí)和實(shí)際工程相結(jié)合，掌握電力系統(tǒng)穩(wěn)態(tài)、暫態(tài)分析的原理、分析步驟以及結(jié)論分析。3．了解并掌握PSDBPA電力系統(tǒng)分析程序單線(xiàn)圖和地理接線(xiàn)圖的使用。我這用了牛頓法直角坐標(biāo)系來(lái)編寫(xiě)程序的。答：高斯迭代法（高斯塞德?tīng)柗ǎ?，牛頓拉夫遜法以及PQ分解法。kind=[1 1 2 0]。end end%數(shù)據(jù)輸入g=[ 0 0 0 0 0 ]。pv(2,2*i1)=2*e(i)。endpv(1,2*i1)=sg(i,i)*e(i)b(i,i)*f(i)。end endfunction pv=PV_out(g,b,e,f,i)%計(jì)算pv節(jié)點(diǎn)的雅克比矩陣 l=length(e)。elsepq(1,2*j1)=(g(i,j)*e(i)+b(i,j)*f(i))。endpq(2,2*i1)=s+b(i,i)*e(i)g(i,i)*f(i)。for k=1:ls=s(g(i,k)*f(k)+b(i,k)*e(k))。for j=1:(l1)if j==i s=0。for j=1:(2*l2)J(2*i1,j)=s(1,j)。if kind(1,i)==1s=PQ_out(g,b,e,f,i)。elses(2*i)=kind(3,i)^2f(i)^2e(i)^2。for j=1:ls1=s1+g(i,j)*e(j)b(i,j)*f(j)。endfunction s=flow(g,b,kind,e,f)%計(jì)算當(dāng)前ef與規(guī)定的pqv的差值 l=length(e)。for j=1:ls1=s1+g(i,j)*e(j)b(i,j)*f(j)。ends=flow(g,b,kind,e,f)。dv=J_ni*s。三、實(shí)驗(yàn)程序：function [e,f,p,q]=flow_out(g,b,kind,e,f)%計(jì)算潮流后efpq的終值 s=flow(g,b,kind,e,f)。在各種潮流計(jì)算的算法中選擇一種，按照計(jì)算方法編制程序。第四篇：電力系統(tǒng)分析潮流實(shí)驗(yàn)報(bào)告南昌大學(xué)實(shí)驗(yàn)報(bào)告電力系統(tǒng)潮流計(jì)算實(shí)驗(yàn) 學(xué)生姓名： 學(xué) 號(hào)： 專(zhuān)業(yè)班級(jí)： 實(shí)驗(yàn)類(lèi)型：□ 驗(yàn)證 □ 綜合 ■ 設(shè)計(jì) □ 創(chuàng)新 實(shí)驗(yàn)日期： 實(shí)驗(yàn)成績(jī)：一、實(shí)驗(yàn)?zāi)康模罕緦?shí)驗(yàn)通過(guò)對(duì)電力系統(tǒng)潮流計(jì)算的計(jì)算機(jī)程序的編制與調(diào)試，獲得對(duì)復(fù)雜電力系統(tǒng)進(jìn)行潮流計(jì)算的計(jì)算機(jī)程序，使系統(tǒng)潮流計(jì)算能夠由計(jì)算機(jī)自行完成，即根據(jù)已知的電力網(wǎng)的數(shù)學(xué)模型（節(jié)點(diǎn)導(dǎo)納矩陣）及各節(jié)點(diǎn)參數(shù)，由計(jì)算程序運(yùn)行完成該電力系統(tǒng)的潮流計(jì)算。end end B2。for i=1:nif B2(i,6)==2j=j+1。DetaU=zeros(2*n2,1)。Jacbi(2*h,2*k1)=0。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。for i=1:nif 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:nif B2(i,6)==2h=h+1。end end I。%創(chuàng)建I i=zeros(n1,1)。t=t+1。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。%創(chuàng)建DetaS h=0。B2(i,6)==3h=h+1。B2(i,6)==2h=h+1。%十：開(kāi)始循環(huán)********************************************************************** while abs(max(DetaU))pr OrgS=zeros(2*n2,1)。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。DetaU=inv(Jacbi)*DetaS。Jacbi(2*h,2*k)=0。Jacbi(2*h,2*k)=2*real(B2(i,3))。for j=1:nif j~=isbk=k+1。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。for j=1:nif j~=isbk=k+1。%八：創(chuàng)建Jacbi(雅可比矩陣)