兩機五節(jié)點網(wǎng)絡潮流計算牛拉法設計說明書(編輯修改稿)
2024-08-16 08:07 本頁面
【文章內容簡介】 ，然后轉入牛頓法迭代。 直角牛頓拉夫遜法潮流計算求解過程以下討論的是用直角坐標形式的牛頓—拉夫遜法潮流的求解過程。當采用直角坐標時，潮流問題的待求量為各節(jié)點電壓的實部和虛部兩個分量由于平衡節(jié)點的電壓向量是給定的，因此待求兩共需要2(n1)個方程式。事實上，除了平衡節(jié)點的功率方程式在迭代過程中沒有約束作用以外，其余每個節(jié)點都可以列出兩個方程式。對PQ節(jié)點來說，是給定的，因而可以寫出 （221）對PV節(jié)點來說，給定量是，因此可以列出 （222）求解過程大致可以分為以下步驟：（1）形成節(jié)點導納矩陣Y（2）將各節(jié)點電壓設初值U，（3）將節(jié)點初值代入式(221)或式(222)，求出修正方程式的常數(shù)項向量（4）將節(jié)點電壓初值代入求式，求出雅可比矩陣元素（5）求解修正方程，求修正向量（6）求取節(jié)點電壓的新值（7）檢查是否收斂，如不收斂，則以各節(jié)點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉入下一步（8）計算支路功率分布，PV節(jié)點無功功率和平衡節(jié)點柱入功率。以直角坐標系形式表示①. 迭代推算式 采用直角坐標時,節(jié)點電壓相量及復數(shù)導納可表示為: (223)將以上二關系式代入上式中,展開并分開實部和虛部。假定系統(tǒng)中的第1,2,m號為P—Q節(jié)點,第m+1,m+2,n1為P—V節(jié)點,根據(jù)節(jié)點性質的不同,得到如下迭代推算式: ⑴對于PQ節(jié)點 (224)⑵對于PV節(jié)點 (225)⑶對于平衡節(jié)點 平衡節(jié)點只設一個,電壓為已知,不參見迭代,其電壓為: (226)②. 修正方程式迭代式共包括2(n1),代入方程并按泰勒級數(shù)展開,略去二次方程及以后各項,得到一組線性方程組或線性化了的方程組，常稱修正方程組: (227) (228)③.雅可比矩陣各元素的算式式(228)中, 雅可比矩陣中的各元素可通過對式(224)和(225), 雅可比矩陣中非對角元素為 (229)當時,雅可比矩陣中對角元素為: (2210)由式(229和(2210)看出,雅可比矩陣的特點:⒈矩陣中各元素是節(jié)點電壓的函數(shù),在迭代過程中,這些元素隨著節(jié)點電壓的變化而變化。⒉導納矩陣中的某些非對角元素為零時,則必有。⒊雅可比矩陣不是對稱矩陣。雅可比矩陣各元素的表示如下式: (2211) (2212) (2213) (2214) (2215) (2216)第三章 設計過程題目：設計原始資料網(wǎng)絡接線如圖所示，各支路導納均以標幺值標于圖中。其中節(jié)點1為平衡節(jié)點，保持=+j0為定值：節(jié)點5為PQ節(jié)點，給定的注入功率如圖。計算精度要求各節(jié)點電壓修正量不大于。圖31 系統(tǒng)接線圖圖32 等值阻抗圖解：節(jié)點1為平衡節(jié)點,U1=+J0為一值,其它四個節(jié)點都是PQ節(jié)點給定的注入功率分別為:2 =+,S3=,S4=,S5=.Y =[ + + 0 0+ + + + + + 0 0 + + 0 + 0 + ]。:取 U1=+J0U2=+j0U3=+j0U4=+j0U5=+j0 3，計算雅可比矩陣中各元素當時, 雅可比矩陣中非對角元素為 當時,雅可比矩陣中對角元素為: 把數(shù)據(jù)代入上邊公式，可得：雅克比矩陣各元素H22= J22= N22= L22=H23=5 J23= N23= L23=5H24=5 J24= N24= L24=5H25= J25= N25= L25=相似可得雅可比矩陣中其它元素采用矩陣求逆委員會駝算求各節(jié)點電壓的修正量求得雅可比矩陣的逆陣節(jié)點功率不平衡量節(jié)點電壓修正量,從而求得電壓新值修正方程的解DY: 節(jié)點電壓的第1次近似值: 各點的電壓實部ei(單位:V)為(節(jié)點號從小到大排列): 各點的電壓虛部fi(單位:V)為(節(jié)點號從小到大排列): 0 平衡節(jié)點的功率：S1 = + 參數(shù)矩陣：[1,2,+*j。1,3,+*j。2,3,+*j。2,4,+*j。2,5,+*j。3,4,+*j。4,5,+*j。]源程序：clcclearN=input(39。請輸入節(jié)點總數(shù)為:N=39。)。nl=input(39。請輸入支路數(shù)為：nl=39。)。isb=input(39。請輸入平衡節(jié)點為:isb=39。)。pr=input(39。請輸入迭代精度為:pr=39。)。B1=input(39。請輸入由支路參數(shù)形成的矩陣為:B1=39。)。Y=zeros(N)。for i=1:nl p=B1(i,1)。q=B1(i,2)。 Y(p,q)=Y(p,q)1./B1(i,3)。 Y(q,p)=Y(p,q)。 Y(q,q)=Y(q,q)+1./B1(i,3)。 Y(p,p)=Y(p,p)+1./B1(i,3)。end disp(39。導納矩陣為:Y=
化學相關推薦
鄂ICP備17016276號-1