【正文】 j1,1)*sqrt(1)。%創(chuàng)建DetaS h=0。end end I。for i=1:nif B2(i,6)==3h=h+1。for i=1:nif B2(i,6)==2j=j+1。三、實驗程序：function [e,f,p,q]=flow_out(g,b,kind,e,f)%計算潮流后efpq的終值 s=flow(g,b,kind,e,f)。endfunction s=flow(g,b,kind,e,f)%計算當前ef與規(guī)定的pqv的差值 l=length(e)。for j=1:(2*l2)J(2*i1,j)=s(1,j)。elsepq(1,2*j1)=(g(i,j)*e(i)+b(i,j)*f(i))。end end%數(shù)據(jù)輸入g=[ 0 0 0 0 0 ]。3．了解并掌握PSDBPA電力系統(tǒng)分析程序單線圖和地理接線圖的使用。常用的控制語句主要包括：(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é)果二進制文件名的控制語句“/NEW_BASE, FILE = 文件名”；(7)指定潮流圖和地理接線圖使用的結(jié)果文件控制語句“/PF_MAP，F(xiàn)ILE=文件名”；(8)指定網(wǎng)絡數(shù)據(jù)的控制語句“/NETWORK_DATA”；(9)指定潮流數(shù)據(jù)文件結(jié)束的控制語句“(END)”； （PFO文件）潮流計算結(jié)果文件內(nèi)容主要分下述幾個方面： 1）程序控制語句列表。無功旋轉(zhuǎn)備用只統(tǒng)計有功出力大于0并且基準電壓小于30kV的發(fā)電機。牛頓拉夫遜方法是解非線性代數(shù)方程組的一種基本方法，在潮流計算中也得到應用。例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é)點算例： ：，并進行計算，結(jié)果如下： 1）系統(tǒng)數(shù)據(jù) 2）計算過程迭代信息及詳細的輸出列表：小結(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ā)電/ /負荷//無功補償/ 類型 擁有者 分區(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é)說明：（如 抽水蓄能電機）。不正常時，將留在硬盤上，可隨時刪除） 儲存一個潮流作業(yè)計算時屏幕顯示的信息。我這用了牛頓法直角坐標系來編寫程序的。pv(2,2*i1)=2*e(i)。endpq(2,2*i1)=s+b(i,i)*e(i)g(i,i)*f(i)。if kind(1,i)==1s=PQ_out(g,b,e,f,i)。for j=1:ls1=s1+g(i,j)*e(j)b(i,j)*f(j)。在各種潮流計算的算法中選擇一種，按照計算方法編制程序。DetaU=zeros(2*n2,1)。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。%創(chuàng)建I i=zeros(n1,1)。B2(i,6)==3h=h+1。DetaU=inv(Jacbi)*DetaS。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。h=0。for i=1:nif B2(i,6)==3t=t+1。end end Y。節(jié)點號和對地參數(shù):X=39。n1=input(39。180。231。106S/km變壓器（每臺）:DP0=86kW,DPs=200kW,I0%=,Vs%=解 最大負荷時變壓器的繞組功率損耗為V%230。10+247。故循環(huán)功率為 Sc187。SA1=SI++(+)MVA=+=++1102DVI=PA1RI+180。=180。=VE=VCDVCE== 圖所示110kV閉式電網(wǎng)，A點為某發(fā)電廠的高壓母線，其運行電壓為117kV。10=(+)kVA=580+j129kVA =*17180。=180。180。S1246。15DPo+jDQo=2180。80180。+180。,i)disp(angle(EF(i))*180/pi)。第%d次迭代后的節(jié)點電壓分別為:n39。F1(n+i1,n+w1)=F1(n+i1,n+w1)G(i,k)*E(k)+B(i,k)*F(k)。%△Qfor w=1:n FX(i)= FX(i)E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)F(i)*G(i,w)*F(w)F(i)*B(i,w)*E(w)。0000 0。%Q=PQ節(jié)點的Q值 PP=[]。%X=△X n=m+v+1。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。)。解：（1）先假設(shè)各節(jié)點電壓均為額定電壓，求線路始端功率。變電所中裝有兩臺三相110/11kV的變壓器，每臺的容量為15MVA,其參數(shù)為：DP0=，DPs=128kW,Vs%=,Io%=。兩臺變壓器并聯(lián)運行時，它們的等值電阻、電抗及勵磁功率分別為1DPsVN21128180。110232。2 于是可得 S139。+180。kV時的初始如圖所示。10+j1200180。1102Mvar=線路Ⅱ： Z213。DSTb242+182(+)MVA=+ =2110Sb=SLDb+DSTb+DSob+jDQBI+jDQBIII=24+j18++++=++92(+)MVA=+ =1102Sc=SLDc+DSTc+DSoc+jDQBIII+jDQBII=12+j9++++=+（3）計算閉式網(wǎng)絡的功率分布。=101247。k1=(10+)180。j1MVA=8+ 21039。232。248。熟悉MATLAB，并掌握MATLAB程序的基本調(diào)試方法。B1=input(39。Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。endend end %三：對PV節(jié)點的處理，注意這時不可再將h初始化為0 for i=1:nif i~=isbamp。for i=1:n%六：對PV節(jié)點的處理，注意這時不可再將h初始化為0if i~=isbamp。if i==j%對角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。else%非對角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3)