【正文】 (39。請輸入支路數:n1=39。)。B1=input(39。請輸入節(jié)點參數:B2=39。)。q=B1(i,2)。Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。Y(q,p)=Y(p,q)。%將OrgS、DetaS初始化OrgS=zeros(2*n2,1)。for i=1:n%對PQ節(jié)點的處理if i~=isbamp。endend end %三：對PV節(jié)點的處理，注意這時不可再將h初始化為0 for i=1:nif i~=isbamp。endend end OrgS。PVU(t,1)=B2(i,3)。B2(i,6)==2h=h+1。for i=1:n%六：對PV節(jié)點的處理，注意這時不可再將h初始化為0if i~=isbamp。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。for i=1:nif i~=isbh=h+1。h=0。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))。endif k==(n1)%將用于內循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0。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))。endif k==(n1)%將用于內循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0。DetaU。end end for i=1:n%對PV節(jié)點的處理if B2(i,6)==3j=j+1。%！!初始功率參數在迭代過程中是不累加的，所以在這里必須將其初始化為零矩陣h=0。for j=1:nOrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。for j=1:nOrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。for i=1:nif i~=isbamp。end end t=0。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。h=0。%創(chuàng)建Jacbi Jacbi=zeros(2*n2)。for j=1:nif j~=isbk=k+1。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。for j=1:nif j~=isbk=k+1。Jacbi(2*h,2*k)=2*real(B2(i,3))。Jacbi(2*h,2*k)=0。DetaU=inv(Jacbi)*DetaS。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。Times=Times+1。通過實驗教學加深學生對復雜電力系統(tǒng)潮流計算計算方法的理解，學會運用電力系統(tǒng)的數學模型，掌握潮流計算的過程及其特點，熟悉各種常用應用軟件，熟悉硬件設備的使用方法，加強編制調試計算機程序的能力，提高工程計算的能力，學習如何將理論知識和實際工程問題結合起來。將事先編制好的電力系統(tǒng)潮流計算的計算程序原代碼由自備移動存儲設備導入計算機。k=0。l=length(dv)/2。endl=length(e)。s2=s2+g(i,j)*f(j)+b(i,j)*e(j)。s=zeros(2*l2,1)。s2=s2+g(i,j)*f(j)+b(i,j)*e(j)。end endfunction J=J_out(g,b,kind,e,f)%計算節(jié)點的雅克比矩陣 l=length(e)。for j=1:(2*l2)J(2*i1,j)=s(1,j)。J(2*i,j)=s(2,j)。for k=1:ls=s(g(i,k)*e(k)b(i,k)*f(k))。endpq(1,2*i)=s+b(i,i)*e(i)g(i,i)*f(i)。s=0。pq(1,2*j)=b(i,j)*e(i)g(i,j)*f(i)。pv=zeros(2,2*l2)。s=0。pv(2,2*i)=2*f(i)。b=[ 0 0 0 ]。[e,f,p,q]=flow_out(g,b,kind,e,f)。高斯迭代法是直接迭代，對初值要求比較低，程序簡單，內存小，但收斂性差，速度慢，多用于配電網或輻射式網絡中；牛頓拉夫遜法是將非線性方程線性化之后再迭代的，對初值要求比較高，收斂性好，速度快，迭代次數少，運行時間短，被廣泛使用；PQ分解法是在極坐標牛頓法的基礎上進行三個簡化所得，有功、無功分開迭代，迭代次數比牛頓多一倍但運算量小，整體速度更快，運行時間更短，多用于110KV以上的高壓電網中，電力網絡的節(jié)點分幾類？各類節(jié)點的已知量和待求量是什么？答: PQ節(jié)點：P、Q為已知量，V、d為待求量；PV節(jié)點：給定P、V，求Q、d；平衡節(jié)點：給定V、d，求P、Q。通過編寫這次程序可以更深一步的理解潮流計算的步驟，也明白了在潮流計算中要注意的一些細節(jié)。二、實驗背景隨著科學技術的飛速發(fā)展，電力系統(tǒng)也在不斷地發(fā)展，電網通過互聯(lián)變得越來越復雜，同時也使系統(tǒng)穩(wěn)定問題越來越突出。清晰認知電力系統(tǒng)分析的意義?？蛇M行交流系統(tǒng)潮流計算，也可進行包括雙端和多端直流系統(tǒng)的交直流混合潮流計算。 定義潮流程序卡片格式文件，用戶可更改及調整該文件。2）輸入、輸出文件及輸出的內容列表。5）迭代過程。IEEE 9節(jié)點例題：圖1 IEEE9節(jié)點系統(tǒng)接線圖節(jié)點參數、線路參數及變壓器參數分別見表1～表3。有功出力為負值的發(fā)電機（包括電動機）作為負荷處理，不統(tǒng)計在內。* 潮流計算迭代過程和平衡節(jié)點相關信息數據計算結果收斂。獲益匪淺。潮流計算的方法有最基本的手算迭代方法，而利用電子計算機進行潮流計算從20世紀50年代中期就已經開始。PSDBPA仿真軟件中潮流計算模型建模的注意事項?。從數學上說，潮流計算是求解一組由潮流方程描述的非線性代數方程組。它的任務是根據給定的運行條件和網路結構確定整個系統(tǒng)的運行狀態(tài)，如各母線上的電壓（幅值及相角）、網絡中的功率分布以及功率損耗等。各區(qū)域平衡機出力數據列表區(qū)域平衡機電壓額定有功有功出力無功出力有功負荷無功負荷所屬分區(qū)SYSTEM發(fā)電機1 01* 沒有遇到錯誤信息 23:03:48 3）單線圖：（二）課本習題：E25 ：：(POWERFLOW,CASEID=IEEE9,PROJECT=IEEE_9BUS_TEST_SYSTEM)/SOL_ITER,DECOUPLED=2,NEWTON=15,OPITM=0 /P_OUTPUT_LIST,ZONES=ALL /RPT_SORT=ZONE /NEW_BASE,FILE= /PF_MAP,FILE = /母線4999999B母線1 B母線2 BE母線3 999.L支路數據L母線1母線2 L母線2母線4 L母線4母線1 .T變壓器數據，包括普通變壓器、移相器、帶調節(jié)的變壓器等。對應于上述系統(tǒng)及數據的潮流計算數據（）見例1。分析報告列表：并根據用戶的要求（通過控制語句控制），輸出各種潮流分析報告。如為致命性錯誤，則中斷計算。打開TextEdit應用程序時先讀入該文件。、輸出相關文件 *.dat潮流計算數據文件*.bse潮流計算二進制結果文件（可用于潮流計算的輸入或穩(wěn)定計算）*.pfo潮流計算結果文件*.map 供單線圖格式潮流圖及地理接線圖格式潮流圖程序使用的二進制結果文件*.pff，*.pfd 中間文件（正常計算結束后將自動刪除。采用什么算法以及迭代的最大步數可以由用戶指定。PSDBPA軟件包主要由潮流和暫穩(wěn)程序構成，具有計算規(guī)模大、計算速度快、數值穩(wěn)定性好、功能強等特點，已在我國電力系統(tǒng)規(guī)劃、調度、生產運行及科研部門得到了廣泛應用。2．了解和掌握PSDBPA電力系統(tǒng)分析程序穩(wěn)態(tài)分析方法（即潮流計算）。六、實驗心得這次實驗是通過matlab編寫出一個潮流計算的程序。 容許誤差ε=105節(jié)點導納矩陣：各節(jié)點電壓：節(jié)點efvζ 各節(jié)點功率：節(jié)點PQ 2 – 結果：五、思考討論題？簡述各種算法的優(yōu)缺點。f=[0 0 0 0]。pv(1,2*j)=b(i,j)*e(i)g(i,j)*f(i)。endpv(1,2*i)=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))。pq(2,2*j1)=pq(1,2*j)。endpq(2,2*i)=s+g(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。pq=zeros(2,2*l2)。endelses=PV_out(g,b,e,f,i)。for i=1:(l1)。if kind(1,i)==1s(2*i)=kind(3,i)f(i)*s1+e(i)*s2。s2=0。q(i)=f(i)*s1e(i)*s2。s2=0。f(i)=f(i)dv(2*i)。J_ni=inv(J