【正文】 Matlab 中的變量 ① 變量命名規(guī)則： 1)變量名是不包含空格的單個詞 2)變量名區(qū)分大小寫 3)變量名必須以字母開頭的字母、數(shù)字、下劃線的組合，最多 19 個字符。他可以高效率的解決工業(yè)計算問題，特別是關(guān)于矩陣和矢量的計算。 Matlab 是一種交換式、面向?qū)ο蟮脑O(shè)計語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運算，同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強大的 功能。 ④ 適應(yīng)性好，亦能處理變壓器變比調(diào)整、系統(tǒng)元件的不同描述和與其他程序配合的能力強。 ② 內(nèi)存需要少?，F(xiàn)在潮流計算有很多方法。但其最終的計算精確度是不受影響的，因為計算的精度取決于最終的誤差要求 1? 和 2? ，如果誤差要求和牛頓－拉夫遜法是一樣的，那么 PQ 分解法最終的計算結(jié)果和牛頓－拉夫遜法的計算結(jié)果的精度就是一樣的。 ③ 以對稱的系數(shù)矩陣 ‘B 和 B 代替變化的系 數(shù)矩陣 J，是逆運算量所需的儲存容量都大為減少。 PQ 分解法的特點 ① 以一個（ n1）階和一個（ m1）階系數(shù)矩陣 ‘B 和 B 代替原來的（ n+m2）階系數(shù)矩陣 J，提高了計算速度，降低了對儲存容量的要求。 河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 復(fù)雜電力系統(tǒng)潮流計算的計算機算法 21 ⑨ 運用各節(jié)點電壓的新值自第三步開始進(jìn)入下一次迭代 。 ⑦ 解修正方程，求節(jié)點電壓大小的變量 )0(iV? 。 ⑤ 求節(jié)點電壓相位角的新值 )1(i?? = )0(i? + )0(i?? 。 ② 設(shè)節(jié)點電壓的初值 )0(i? 和 )0(iV ③ 計算有功功率的不平衡量 )0(iP? ，從而求出 )0(iP? / )(oiV 。B 和 39。39。 2）由于約束條件m a xjiji ???? ???， ijji ??? ?? ，不宜過大，可認(rèn)為0?ij? 通過上述的兩個化簡可以將雅克比矩陣化為對稱的矩陣 39。因此，為了減少計算量，根據(jù)基于極坐標(biāo)的牛頓 －拉夫遜法的特點，建立了 PQ分解 法的潮流計算方法。顯然，雖修正方程式有兩種不同表示方式，但牛頓一拉夫遜法潮流計算的基本步驟卻總不外乎如下幾步： ① 形成節(jié)點導(dǎo)納矩陣； ② 設(shè)各節(jié)點電壓的初值； ③ 把節(jié)點電壓的初值代入功率方程，求修正方程式中的不平衡量； ④ 將各節(jié)點電壓的初值代入雅可比矩陣系數(shù)求解公式求修正方程式的系數(shù)矩陣 —— 雅可比矩陣 H、 N、 J、 L 塊的各個元素； ⑤ 解修正方程式，求各節(jié)點電壓的變化量，即修正量； ⑥ 計算各節(jié)點電壓的新值，即修正后的值； 河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 復(fù)雜電力系統(tǒng)潮流計算的計算機算法 20 ⑦ 運用各節(jié)點電壓的新值自第三步開始進(jìn)入下一次迭代； ⑧ 計算平衡節(jié)點功率和線路功率 PQ 分解法潮流計算 PQ 分解法潮流計算概述 由于該設(shè)計是運用牛頓拉夫遜直角坐標(biāo)法進(jìn)行的，對 PQ 分解法只做簡單的介紹。 ③ 雅可比矩陣不是對稱矩陣 。22 式（ ） 雅可比矩陣的特點 : ① 矩陣中各元素是節(jié)點電壓的函數(shù) ,在迭代過程中 ,這些元素隨著節(jié)點電壓的變化而變化 。2。 式 () ij? 時，雅克比矩陣的元素為： ????????????????????????????????????????????0。 運用牛頓拉夫遜法時，可直接用以求解功率方程。 由此可得 ?),( )0()0(2)0(11 nxxxf ? 10120211011 yxxfxxfxxf nn????????????? ? 式（ ） 改寫為 ?1y ?),( )0()0(2)0(11 nxxxf ?（011xf?? +021xf?? +nxf??1 1)x? 式（ ） 化為矩陣形式為： 河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 復(fù)雜電力系統(tǒng)潮流計算的計算機算法 17 ????????????????????????????????????????????????????????????????????????????nnnnnnxxxxfxfxfxnfxfxfxnfxfxfxxxfyxxxfyxxxfy????????2101021011010210110101011)0()0(2)0(12)0()0(2)0(122)0()0(2)0(111......2),(),(),( 式 （ ） 簡寫 為 xJf ??? 式 （ ） ① 牛頓 — 拉夫遜法的特點 1)牛頓 拉夫遜法是迭代法，逐漸逼近的方法； 2)修正方程是線性化方程，它的線性化過程體現(xiàn)在把非線性方程在 (0)x 按泰勒級數(shù)展開，并略去高階小量； 3)用牛頓 — 拉夫遜法解題時，其初始值要求嚴(yán)格 (較接近真解 )，否則迭代不收斂。以任何一式為例， 110120211011)0()0(2)0(11)0(2)0(21)0(11 ),(),(yxxfxxfxxfxxxfxxxxxxfnnnnn?????????????????????????? 式中 :01021011 ,nxfxfxf ?????? ? 分別 表示以 )0()0(2)0(1 , nxxx ? 代入這些偏導(dǎo)數(shù)時的所計算得 1? 則是包含nxxx ??? , 21 ? 的高次方與 1f 的高階偏導(dǎo)數(shù)乘積的函數(shù)。 河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 復(fù)雜電力系統(tǒng)潮流計算的計算機算法 16 設(shè)有非線性方程組 ??????????nnnnnyxxxfyxxxfyxxxf),(),(),(2122121211????? 式 （ ） 其近似解 )0(1x ， )0(2x ， ? ， )0(nx 。如果不能滿足要求，則應(yīng)修改某些變量的給定值，甚至修改系統(tǒng)的 運行方式，重新進(jìn)行計算。 因此，潮流計 算可以歸結(jié)為求解一組非線性方程組，并使其解答滿足一定的約束條件。節(jié)點之間電壓的相位差應(yīng)滿足 ： m a x| | | | | |ij i j i j? ? ? ? ?? ? ? ? 式 （ ） 為了保證系統(tǒng)運行的穩(wěn)定性，要求某些輸電線路兩端的電壓相位 差 不超過一定 的數(shù)值。因此，這一約束條件對 PQ 節(jié)點而言。這些要求構(gòu) 成了潮流問題中某些變量的約束條件，常用的約束條件如下： 節(jié)點電壓應(yīng)滿足 : m axm in iii VVV ?? （ ni ,...2,1? ） 式 （ ） 從保證電能質(zhì)量和供電安全的要求來看，電力系統(tǒng)的所有電氣設(shè)備都必須運行在額定電壓附近。 河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 復(fù)雜電力系統(tǒng)潮流計算的計算機算法 15 以上三種節(jié)點四個運行參數(shù) P、 Q、 V、θ中已知量都是兩個，待求量也是兩個，只是類型不同而已。 關(guān)于平衡節(jié)點的選擇，一般選擇系統(tǒng)中擔(dān)任調(diào)頻調(diào)壓的某一發(fā)電廠（或發(fā)電機）。對該節(jié)點，給定其電壓值，并在計算中取該節(jié)點電壓向量的方向作為參考軸，相當(dāng)于給定該點電壓向量的角度為零。 PV 節(jié)點上的發(fā)電機稱之為 PV 機（或PV 給定型發(fā)電機）。 這種節(jié)點在運行中往往要有一定可調(diào)節(jié)的無功電源，用以維 持給定的電壓值。在潮流計算中，系統(tǒng)大部分節(jié)點屬于 PQ 節(jié)點）。當(dāng)某些發(fā)電機的出力 P、 Q給定時，也作為 PQ節(jié)點。所以叫“ PQ 節(jié)點”。所以，根據(jù)電力系統(tǒng)中各節(jié)點性質(zhì)的不同，很 自然地把節(jié)點分成三種類型。然而在電力系統(tǒng)中，給出發(fā)電機或負(fù)荷連接母線上電壓或電流（都是向量）的情況是很少的，一般是給出發(fā)電機母線上發(fā)電機的有功功率（ P）和母線電壓的幅值（ V），給出負(fù)荷母線上負(fù)荷消耗的有功功率（ P）和無功功率（ Q）。?k 的變 壓器的計算公式。 11； 0??jj ；Tjiij YkkYY ???????? ?????? ? *39。?k 。 ； ijijjiij yyYY ????? 39。ijy ，如下圖 圖 這種情況節(jié)點 i、 j有關(guān)的元素應(yīng)作如下修改 ： ijijii yyY ??? 39。 新增的對角元 ijY ，由于在節(jié)點 j 只有一個支路 ijy ，將為 ijY = ijy ； 新增的 非 對角 元 ijY = jiY ，則為 ijY = jiY = ijy ；原有矩陣中的對角元 iiY 將增加△ iiY = ijy 。 1）從原有網(wǎng)絡(luò)引出一支路，同時增加一節(jié)點，如下圖 所示 圖 設(shè) i 為原有網(wǎng)絡(luò)中節(jié)點， j 為新增加節(jié)點，新增加支路導(dǎo)納為 ijy 。由于改變一個支路的參數(shù)或它的投入、退出狀態(tài)只影響該支路兩端節(jié)點的自導(dǎo)納和它們之間的互導(dǎo)納，可不必重新形成與新運行狀況相對應(yīng)的節(jié)點導(dǎo)納矩陣，僅需對原有的矩陣作某些修改。 4）節(jié)點導(dǎo)納矩陣是方陣，其階數(shù)就等于網(wǎng)絡(luò)中出去參考節(jié)點外的節(jié)點數(shù) n。因此，一般情況下，節(jié)點導(dǎo)納矩陣的對角元素往往大于非對角元素的負(fù)值； 3）節(jié)點導(dǎo)納矩陣一般是對稱矩陣，這是網(wǎng)絡(luò)的互易特性所決定的。 ② 節(jié)點導(dǎo)納矩陣的形成 1）節(jié)點導(dǎo)納矩陣的對角元素就等于各該節(jié)點所連接導(dǎo)納的總和。而且，由于每個節(jié)點所連接的支路數(shù)總有一定限度，隨著網(wǎng)絡(luò)中節(jié)點數(shù)的增加，非零元素數(shù)相對越來越少，節(jié)點導(dǎo)納矩陣的稀疏度，即零元素數(shù)與總元素數(shù)的比值也越來越高。而且，如節(jié)點j、 i之間沒有直接聯(lián)系，也不計兩支路之間，例如兩相鄰電力線路之間有互感時，jiY = ijY =0。由此可見 ,節(jié)點 j、 i 之間的互導(dǎo)納 jiY數(shù)值上就等于連接節(jié)點 j、 i支路導(dǎo)納的負(fù)值。而由式（ ）可見，互導(dǎo)納 jiY 在數(shù)值上就等于在節(jié)點 i 施加單位電壓 ，其它節(jié)點全部接地時，經(jīng)節(jié)點 j 注入網(wǎng)絡(luò)的電流。由此可見，節(jié)點 i 的自導(dǎo)納 iiY 數(shù)值上就等于與該節(jié)點直接相連接的所有支路導(dǎo)納的總和。由式 ()可見 ,自導(dǎo)納 iiY 數(shù)值上就等于節(jié)點 i 施加單位電壓 ,其它節(jié)點全部接地時 ,經(jīng)節(jié)點 i 注入網(wǎng)絡(luò)的電流。網(wǎng)絡(luò)中沒有接地支路時 ,各節(jié)點電壓可指各該節(jié)點與某一個被選定作參考節(jié)點之間的對地電壓差。以下僅介紹本設(shè)計所用到的節(jié)點電壓方程。在電路中，我們曾運用過節(jié)電電壓方程、回路電流方程等。 設(shè)節(jié)點 4的電壓差為 ： 用簡化的回路電流法解簡化等值電路 通過近似方法，從功率中求取相應(yīng)的電流，電壓近似認(rèn) 為是額定電壓： 流經(jīng)阻 抗 12Z 功率為： S a S c S b U 1 2 3 U 4 1 Z 12 Z 23 Z 34 4 S 2