【正文】 J(no:)=[]。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。X4=X1。End End case 2 Function for j=1:n J=jacci(Y,G,B,PQ,e,f,V,C,D,B2,n,ph,no)X1=G(i,j)*f(i)B(i,j)*e(i)。p=p+1。b(i)=v(B1(i1))*a(i)j*B1(i4)*v(B1(i))^2/2。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)?！尽俊尽縁unction tisco %這是一個電力系統(tǒng)潮流計算的程序 n=input（‘n請輸入節(jié)點(diǎn)數(shù)：n=’）； m=input(‘請輸入支路數(shù)：m=’)。3）錯誤信息。除了潮流計算功能外，該軟件還具有自動電壓控制、聯(lián)絡(luò)線功率控制、系統(tǒng)事故分析（N1開斷模擬）、網(wǎng)絡(luò)等值、靈敏度分析、節(jié)點(diǎn)P－V、Q－V和P－Q曲線、確定系統(tǒng)極限輸送水平、負(fù)荷靜特性模型、靈活多樣的分析報告、詳細(xì)的檢錯功能等功能。實(shí)驗(yàn)要求：。0j206。直流潮流模型是把非線性電力潮流問題簡化為線形電路問題，從而使分析計算非常方便，直流潮流專門用于研究電網(wǎng)中有功潮流的分布。S% 求節(jié)點(diǎn)注入電流I=Y*U39。endendend% 當(dāng)i=j時，求H,N,M,Lfor i=1:4for j=1:5if i~=j H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。k=0。% 求節(jié)點(diǎn)導(dǎo)納矩陣中互導(dǎo)納for i=1:5for j=1:5if i~=jY(i,j)=y(i,j)。牛頓迭代法其最大優(yōu)點(diǎn)是在方程f(x)=0的單根附近時誤差將呈平方減少，而且該法還可以用來求方程的重根、復(fù)根。高斯迭代法是直接迭代，對初值要求比較低，程序簡單，內(nèi)存小，但收斂性差，速度慢，多用于配電網(wǎng)或輻射式網(wǎng)絡(luò)中；牛頓拉夫遜法是將非線性方程線性化之后再迭代的，對初值要求比較高，收斂性好，速度快，迭代次數(shù)少，運(yùn)行時間短，被廣泛使用；PQ分解法是在極坐標(biāo)牛頓法的基礎(chǔ)上進(jìn)行三個簡化所得，有功、無功分開迭代、將一個變系數(shù)的2n階J陣轉(zhuǎn)化成兩個常系數(shù)且對稱的n階子陣，迭代次數(shù)比牛頓多一倍但運(yùn)算量小，整體速度更快，運(yùn)行時間更短，多用于110KV以上的高壓電網(wǎng)中。if max(max(abs(op)),max(abs(oq)))x=A(3,4)*A(3,j)*(G(4,j)*cos(2*pi/360*(A(4,4)A(4,j)))+B(4,j)*sin(2*pi/360*(A(4,4)A(4,j))))。sum=sum+h。y(j,m)=y(m,j)。yd=zeros(n,n)。DVn248。230。231。(13)247。231。232。247。DVn248。230。231。230。247。231。V2231。231。230。2121222(9)H=L=231。 235。第一篇：實(shí)驗(yàn)二 電力系統(tǒng)潮流計算實(shí)驗(yàn)電力系統(tǒng)分析實(shí)驗(yàn)報告學(xué)生姓名： 學(xué) 號： 專業(yè)班級：實(shí)驗(yàn)類型：□ 驗(yàn)證 □ 綜合 ■ 設(shè)計 □ 創(chuàng)新 實(shí)驗(yàn)日期： 2012528 實(shí)驗(yàn)成績：一、實(shí)驗(yàn)?zāi)康模罕緦?shí)驗(yàn)通過對電力系統(tǒng)潮流計算的計算機(jī)程序的編制與調(diào)試，獲得對復(fù)雜電力系統(tǒng)進(jìn)行潮流計算的計算機(jī)程序，使系統(tǒng)潮流計算能夠由計算機(jī)自行完成，即根據(jù)已知的電力網(wǎng)的數(shù)學(xué)模型（節(jié)點(diǎn)導(dǎo)納矩陣）及各節(jié)點(diǎn)參數(shù)，由計算程序運(yùn)行完成該電力系統(tǒng)的潮流計算。231。V100231。247。247。247。B2n247。B11B12L230。247。DV1246。230。 OO247。就可得到P1230。Pn247。231。B21B22L231。DV1246。以上兩式就是PQ分解法達(dá)到修正方程式，其中系數(shù)矩陣只不過是系統(tǒng)導(dǎo)納矩陣的虛部，因而是對稱矩陣，而且在迭代過程中維持不變。y=zeros(n,n)。yd(m,j)=(k(m,j)1)*k(m,j)/z(m,j)。end op(1,m)=A(1,m)sum。sum=sum+x。2．在潮流計算中，電力網(wǎng)絡(luò)的節(jié)點(diǎn)分幾類？各類節(jié)點(diǎn)的已知量和待求量是什么？ 答:根據(jù)給定的控制變量和狀態(tài)變量的電力網(wǎng)絡(luò)的節(jié)點(diǎn)可分為以下幾類:PQ節(jié)點(diǎn)（負(fù)荷節(jié)點(diǎn)）：Pi、Qi為已知量，Vi、di為待求量；（該類節(jié)點(diǎn)數(shù)量最多）如：負(fù)荷節(jié)點(diǎn)、變電站節(jié)點(diǎn)（聯(lián)絡(luò)節(jié)點(diǎn)、浮游節(jié)點(diǎn)）、給定P、Q的發(fā)電機(jī)節(jié)點(diǎn)和給定QGi的無功電源節(jié)點(diǎn)。電力系統(tǒng)潮流計算，一般來說，各個母線所供負(fù)荷的功率是已知的，各個節(jié)點(diǎn)的電壓是未知的（平衡節(jié)點(diǎn)外）可以根據(jù)網(wǎng)絡(luò)結(jié)構(gòu)形成節(jié)點(diǎn)導(dǎo)納矩陣，然后由節(jié)點(diǎn)導(dǎo)納矩陣列寫功率方程，由于功率方程里功率是已知的，電壓的幅值和相角是未知的，這樣潮流計算的問題就轉(zhuǎn)化為求解非線性方程組的問題了。endendend% 求節(jié)點(diǎn)導(dǎo)納矩陣中自導(dǎo)納for i=1:5Y(i,i)=sum(y(i,:))。while max(x2)1e6for i=1:4for j=1:4H(i,j)=0。N(i,i)=N(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。四運(yùn)行結(jié)果節(jié)點(diǎn)導(dǎo)納矩陣經(jīng)過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點(diǎn)電壓的幅值和相角（弧度數(shù)）節(jié)點(diǎn)注入功率和電流五 結(jié)果分析在這次學(xué)習(xí)和實(shí)際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機(jī)算法中的牛頓拉夫遜法進(jìn)行深入的研讀，弄明白了其原理、計算過程、公式推導(dǎo)以及設(shè)計流程。二、直流潮流算法的形成過程對下圖所示等值電路圖，對于之路（i,j），如果忽略其并聯(lián)支路，例如忽略線路的充電電容。i,j185。、輸出相關(guān)文件 *.dat潮流計算數(shù)據(jù)文件*.bse潮流計算二進(jìn)制結(jié)果文件（可用于潮流計算的輸入或穩(wěn)定計算）*.pfo潮流計算結(jié)果文件*.map 供單線圖格式潮流圖及地理接線圖格式潮流圖程序使用的二進(jìn)制結(jié)果文件*.pff，*.pfd 中間文件（正常計算結(jié)束后將自動刪除。如為致命性錯誤，則中斷計算。ph=input(‘n請輸入平衡母線的節(jié)點(diǎn)號：ph=’)； B1=input（‘n請輸入支路信號：B1=’）。J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)。c(i)=v(B1(i2))*a(i)j*B1(i4)*v(B1(i2))^2/2。Df(p)=Q(i)Q1。X2=G(i,j)*e(i)+B(i,j)*f(i)。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。J(m,q)=X3。J(m,q)=X4。J(p,q)=X2。j=jph m=p+1。J(m,q)=X4。De(k)=DX(i+1)。Df(p)=P(i)P1。for i=1:ma(i)=conj((v(B1(i1))/B1(i5)v(B1(i2))/B1(i3))。Q=imag(S)。五、思考題1． 電力系統(tǒng)潮流計算方法有哪些？各種方法的原理？ 2． 電力系統(tǒng)潮流計算的作用？3． PSDBPA仿真軟件中潮流計算模型中不同控制語句的作用？ 4． PSDBPA仿真軟件中潮流計算模型建模的注意事項(xiàng)?第五篇：電力系統(tǒng)潮流計算程序設(shè)計電力系統(tǒng)潮流計算程序設(shè)計姓名：韋應(yīng)順學(xué)號：2011021052 電力工程學(xué)院牛頓—拉夫遜潮流計算方法具有能夠?qū)⒎蔷€性方程線性化的特點(diǎn)，而使用MATLAB語言是由于MATLAB語言的數(shù)學(xué)邏輯強(qiáng)，易編譯。2）輸入、輸出文件及輸出的內(nèi)容列表?？蛇M(jìn)行交流系統(tǒng)潮流計算，也可進(jìn)行包括雙端和多端直流系統(tǒng)的交直流混合潮流計算。15）通過式（2）的潮流計算公式求出各條支路的功率情況。x239。雖然這些方法都具有一定的精度，但計算量較大，這顯然不適應(yīng)形成電網(wǎng)規(guī)劃方案時多次而反復(fù)的潮流計算要求。endS(5)=P(5)+Q(5)*sqrt(1)。M(i,j)=N(i,j)。x2=ones(8,1)。endendY=0。二 牛頓拉夫遜法潮流計算 1 基本原理牛頓法是取近似解x(k)之后，在這個基礎(chǔ)上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根。答:潮流計算目前比較主要的方法有三種:高斯迭代法（高斯塞德爾法），牛頓拉夫遜法以及PQ分解法。A([3],[1,2])=A([3],[1,2])+b。for j=1:nh=A(3,m)*A(3,j)*(G(m,j)*cos(2*pi/360*(A(4,m)A(4,j)))+B(m,j)*sin(2*pi/360*(A(4,m)A(4,j))))。end end end for m=1:n for j=1:n if k(m,j)~=0 y(m,j)=k(m,j)/z(m,j)。Y=zeros(n,n)。232。B1n246。DVQ2247。V2Dq2247。M231。248。=231。232。將以上兩式的左右兩側(cè)用以下矩陣左乘1B1n246。M247。0=Q247。231。MM231。=P1246。231。247。2VVBVVBVBL2n2n247。DV/V通過實(shí)驗(yàn)教學(xué)加深學(xué)生對復(fù)雜電力系統(tǒng)潮流計算計算方法的理解，學(xué)會運(yùn)用電力系統(tǒng)的數(shù)學(xué)模型，掌握潮流計算的過程及其特點(diǎn)，熟悉各種常用應(yīng)用軟件，熟悉硬件設(shè)備的使用方法，加強(qiáng)編制調(diào)試計算機(jī)程序的能力，提高工程計算的能力，學(xué)習(xí)如何將理論知識和實(shí)際工程問題結(jié)合起來。235。M247。247。247。231。231。231。1246。231。247。V1246。231。D246。231。247。2247。247。它們與功率誤差方程式j(luò)=nj=1DPi=PisVi229。z(1,2)=+*i。yd(j,m)=(1k(