版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、<p><b> 1 潮流計(jì)算概述</b></p><p> 在電力系統(tǒng)運(yùn)行和規(guī)劃中,都需要研究電力系統(tǒng)穩(wěn)定運(yùn)行情況,確定電力系統(tǒng)的穩(wěn)態(tài)運(yùn)行狀態(tài)。給定電力系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)、參數(shù)和決定電力系統(tǒng)運(yùn)行狀況的邊界條件,,確定電力系統(tǒng)運(yùn)行的方法之一是朝流計(jì)算。從數(shù)學(xué)上說(shuō),朝流計(jì)算是要求解一組有潮流方程描述的非線性方程組。電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)分析中最重要最基本的計(jì)算,是電力運(yùn)行、規(guī)劃以
2、及安全性、可靠性分析和優(yōu)化的基礎(chǔ),也是各種電磁暫態(tài)和機(jī)電暫態(tài)分析的基礎(chǔ)和出發(fā)點(diǎn)。 潮流計(jì)算方法的發(fā)展是與人們所使用的計(jì)算工具的發(fā)展相聯(lián)系的。早期,除了手動(dòng)潮流外,人們用交流計(jì)算臺(tái)通過(guò)物理模擬的方法來(lái)分析電力系統(tǒng)穩(wěn)態(tài)運(yùn)行狀態(tài)。這種方法雖然直觀,物理概念清楚,但受到系統(tǒng)規(guī)模等因素的影響,分析大網(wǎng)絡(luò)的超流會(huì)遇到困難。 作為研究電力系統(tǒng)穩(wěn)定運(yùn)行情況的一種基本電氣計(jì)算 ,電力系統(tǒng)常規(guī)潮流計(jì)算的任務(wù)是根據(jù)給定的網(wǎng)絡(luò)結(jié)構(gòu)及運(yùn)行條件,求
3、出整個(gè)網(wǎng)絡(luò)的運(yùn)行狀態(tài),其中包括各母線的電壓、網(wǎng)絡(luò)中的功率分布以及功率損耗等。 潮流計(jì)算的結(jié)果,無(wú)論是對(duì)于現(xiàn)有系統(tǒng)運(yùn)行方式的分析研究,還是對(duì)規(guī)劃中供電方案的分析比較,都是必不可少的。它為判別這些運(yùn)行方式計(jì)規(guī)劃設(shè)計(jì)方案的合理性、安全可靠性</p><p><b> 2 題目要求</b></p><p> 在圖1所示的簡(jiǎn)單電力系統(tǒng)中,系統(tǒng)中節(jié)點(diǎn)1、2、3為節(jié)點(diǎn),
4、節(jié)點(diǎn)4為節(jié)點(diǎn),節(jié)點(diǎn)5為平衡節(jié)點(diǎn),已給定,,,,,,網(wǎng)絡(luò)各元件參數(shù)的標(biāo)幺值如表1所示,給定電壓的初始值如表2所示,收斂系數(shù)。試求: 采用極坐標(biāo)下的分解法計(jì)算圖1網(wǎng)絡(luò)的潮流分布。</p><p><b> 圖1 簡(jiǎn)單電力系統(tǒng)</b></p><p> 表2.1 網(wǎng)絡(luò)各元件參數(shù)的標(biāo)幺值</p><p> 表2.2 各節(jié)點(diǎn)電壓(初值)標(biāo)幺值參
5、數(shù) </p><p> 3 P-Q分解計(jì)算方法</p><p> 3.1 形成節(jié)點(diǎn)導(dǎo)納矩陣</p><p> 自導(dǎo)納的形成。對(duì)節(jié)點(diǎn)i其自導(dǎo)納Yii是節(jié)點(diǎn)i以外的所有節(jié)點(diǎn)都接地時(shí)節(jié)點(diǎn)i對(duì)地的總導(dǎo)納。顯然,Yii應(yīng)等于與節(jié)點(diǎn)i相接的各支路導(dǎo)納之和,即</p><p> 式中,yi0為節(jié)點(diǎn)i與零電位節(jié)點(diǎn)之間的支路導(dǎo)納;yij為節(jié)點(diǎn)i與節(jié)點(diǎn)
6、j之間的支路導(dǎo)納。</p><p> 互導(dǎo)納的形成。對(duì)節(jié)點(diǎn)i與節(jié)點(diǎn)k之間的互導(dǎo)納是節(jié)點(diǎn)i、k之間的支路導(dǎo)納的負(fù)值,即不難理解。若節(jié)點(diǎn)i和k沒(méi)有支路直接相連時(shí),便有Yik=0。</p><p> 含變壓器支路的處理。若節(jié)點(diǎn)p、q間接有變壓器,如下圖所示,則可作出其∏型等值電路為:</p><p> 圖1 變壓器∏型等值電路</p><p>
7、; 則p、q的自導(dǎo)納和節(jié)點(diǎn)間的互導(dǎo)納分別為</p><p> 3.2 計(jì)算不平衡功率△P、△Q并形成修正方程式</p><p> 對(duì)每一個(gè)PQ節(jié)點(diǎn)或每一個(gè)PV節(jié)點(diǎn)都根據(jù)下列公式計(jì)算出有功功率增量△P</p><p> 而對(duì)于每一個(gè)PQ節(jié)點(diǎn)還可以根據(jù)下面的公式計(jì)算出無(wú)功功率增量△Q</p><p> 在有功功率增量和無(wú)功功率增量不滿足
8、如下約束條件時(shí)</p><p> 利用PQ分解法則可以形成如下修正方程</p><p> 3.3 利用因子表法求解修正方程</p><p> 在電網(wǎng)計(jì)算中經(jīng)常遇到這樣的問(wèn)題,對(duì)方程組需要反復(fù)多次求解,而每次求解僅改變常數(shù)項(xiàng)F,系數(shù)矩陣保持不變。按照一般的高斯消去法,對(duì)每一改變的常數(shù)項(xiàng),形成包括常數(shù)項(xiàng)及系數(shù)矩陣在內(nèi)的增廣矩陣,然后消去回代求出其解??梢钥闯?,每次
9、對(duì)增廣矩陣中A矩陣元素的消元都是重復(fù)的,為了避免這種重復(fù),我們把對(duì)相同的系數(shù)矩陣重復(fù)進(jìn)行的消去與對(duì)不同的常數(shù)項(xiàng)進(jìn)行的消去分開(kāi)進(jìn)行,因此對(duì)系數(shù)矩陣的消去只需進(jìn)行一次,并在消去的過(guò)程中將對(duì)常數(shù)項(xiàng)進(jìn)行消去運(yùn)算的運(yùn)算因子保存下來(lái),形成所謂因子表,這就是因子表法。因?yàn)橐蜃颖碛涗浟烁咚瓜シ▽?duì)常數(shù)項(xiàng)進(jìn)行消去的全部信息,利用它便可對(duì)不同常數(shù)項(xiàng)進(jìn)行消去,形成上三角矩陣,最后求出全部未知數(shù)。</p><p> 在使用PQ分解法時(shí)
10、,其系數(shù)矩陣是在迭代過(guò)程中保持不變的,所以為了節(jié)省內(nèi)存和縮短運(yùn)算時(shí)間我們采取了因子表法。同時(shí)由于電網(wǎng)的節(jié)點(diǎn)導(dǎo)納矩陣矩陣是稀疏陣和對(duì)稱陣,于是我們可以采取只保存系數(shù)矩陣的上三角陣來(lái)使運(yùn)算更為簡(jiǎn)化。</p><p> 若線性方程組一般形式如下:</p><p> 其中稱為系數(shù)矩陣,稱為未知數(shù)向量, 稱為常數(shù)項(xiàng)向量。將矩陣A的元素進(jìn)行如下處理:</p><p><
11、;b> 得到因子表</b></p><p><b> 其中 ;</b></p><p> 再利用因子表進(jìn)行前代過(guò)程,求出每次迭代后的常數(shù)項(xiàng)。其前代公式是:</p><p><b> 求得向量;</b></p><p> 再由因子表與前代得到的向量F,得到方程組</
12、p><p> 求解出此方程即可得到線性方程組的解向量。</p><p> 3.4 多次迭代最終求得V和以及全線路功率</p><p> 利用上面所介紹的方法求解修正方程組</p><p><b> 可以求得和。</b></p><p> 再利用求得每次迭代后的結(jié)果。多次迭代當(dāng)其滿足約束條件和
13、時(shí),迭代結(jié)束。迭代結(jié)束后即可得到各節(jié)點(diǎn)的V和,再根據(jù)V、來(lái)計(jì)算PV節(jié)點(diǎn)的無(wú)功功率Q和平衡節(jié)點(diǎn)的功率以及網(wǎng)絡(luò)中的功率分布。</p><p> PV節(jié)點(diǎn)及平衡節(jié)點(diǎn)無(wú)功功率計(jì)算公式為:</p><p> 平衡節(jié)點(diǎn)有功功率計(jì)算公式為:</p><p> 以下圖所標(biāo)示的正方向,輸電線路功率的計(jì)算公式如下:</p><p> 圖2 支路功率計(jì)算
14、</p><p> 對(duì)其進(jìn)行實(shí)部虛部進(jìn)行分解可得P、Q計(jì)算公式為:</p><p> 4 Matlab簡(jiǎn)介</p><p> 目前電子計(jì)算機(jī)已廣泛應(yīng)用于電力系統(tǒng)的分析計(jì)算,潮流計(jì)算是其基本應(yīng)用軟件之一?,F(xiàn)有很多潮流計(jì)算方法。對(duì)潮流計(jì)算方法的要求:計(jì)算速度快;內(nèi)存需要少;計(jì)算結(jié)果有良好的可靠性和可信性;適應(yīng)性好,亦即能處理變壓器變比調(diào)整;系統(tǒng)元件的不同描述和與其
15、它程序配合的能力強(qiáng)、簡(jiǎn)單。</p><p> MATLAB是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計(jì)語(yǔ)言,廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界,主要用于矩陣運(yùn)算,同時(shí)在數(shù)值分析、自動(dòng)控制模擬、數(shù)字信號(hào)處理、動(dòng)態(tài)分析、繪圖等方面也具有強(qiáng)大的功能。</p><p> MATLAB程序設(shè)計(jì)語(yǔ)言結(jié)構(gòu)完整,且具有優(yōu)良的移植性,它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。它可以高效率地解決工業(yè)計(jì)算問(wèn)題,特別是關(guān)于矩陣和矢量的計(jì)
16、算。MATLAB與C語(yǔ)言和FORTRAN語(yǔ)言相比更容易被掌握。通過(guò)M語(yǔ)言,可以用類(lèi)似數(shù)學(xué)公式的方式來(lái)編寫(xiě)算法,大大降低了程序所需的難度并節(jié)省了時(shí)間,從而可把主要的精力集中在算法的構(gòu)思而不是編程上。</p><p> 另外,MATLAB提供了一種特殊的工具:工具箱(TOOLBOXES).這些工具箱主要包括:信號(hào)處理(SIGNAL PROCESSING)、控制系統(tǒng)(CONTROL SYSTEMS)、神經(jīng)網(wǎng)絡(luò)(NEU
17、RAL NETWORKS)、模糊邏輯(FUZZY LOGIC)、小波(WAVELETS)和模擬(SIMULATION)等等。不同領(lǐng)域、不同層次的用戶通過(guò)相應(yīng)工具的學(xué)習(xí)和應(yīng)用,可以方便地進(jìn)行計(jì)算、分析及設(shè)計(jì)工作。</p><p> MATLAB設(shè)計(jì)中,原始數(shù)據(jù)的填寫(xiě)格式是很關(guān)鍵的一個(gè)環(huán)節(jié),它與程序使用的方便性和靈活性有著直接的關(guān)系。原始數(shù)據(jù)輸入格式的設(shè)計(jì),主要應(yīng)從使用的角度出發(fā),原則是簡(jiǎn)單明了,便于修改。<
18、/p><p><b> 矩陣的運(yùn)算</b></p><p> 矩陣是MATLAB數(shù)據(jù)存儲(chǔ)的基本單元,而矩陣的運(yùn)算是MATLAB語(yǔ)言的核心,在MATLAB語(yǔ)言系統(tǒng)中幾乎一切運(yùn)算均是以對(duì)矩陣的操作為基礎(chǔ)的。矩陣的基本數(shù)學(xué)運(yùn)算包括矩陣的四則運(yùn)算、與常數(shù)的運(yùn)算、逆運(yùn)算、行列式運(yùn)算、秩運(yùn)算、特征值運(yùn)算等基本函數(shù)運(yùn)算,這里進(jìn)行簡(jiǎn)單介紹。</p><p>&
19、lt;b> 四則運(yùn)算</b></p><p> 矩陣的加、減、乘運(yùn)算符分別為“+,—,*” ,用法與數(shù)字運(yùn)算幾乎相同,但計(jì)算時(shí)要滿足其數(shù)學(xué)要求 在MATLAB中矩陣的除法有兩種形式:左除“\”和右除“/”。在傳統(tǒng)的MATLAB算法中,右除是先計(jì)算矩陣的逆再相乘,而左除則不需要計(jì)算逆矩陣直接進(jìn)行除運(yùn)算。通常右除要快一點(diǎn),但左除可避免被除矩陣的奇異性所帶來(lái)的麻煩。</p><
20、p> 與常數(shù)的運(yùn)算 </p><p> 常數(shù)與矩陣的運(yùn)算即是同該矩陣的每一元素進(jìn)行運(yùn)算。但需注意進(jìn)行數(shù)除時(shí),常數(shù)通常只能做除數(shù)。</p><p><b> 基本函數(shù)運(yùn)算</b></p><p> 矩陣的函數(shù)運(yùn)算是矩陣運(yùn)算中最實(shí)用的部分,常用的主要有以下幾個(gè):</p><p> det(a)
21、 求矩陣a的行列式</p><p> eig(a) 求矩陣a的特征值</p><p> inv(a)或a ^ (-1) 求矩陣a的逆矩陣</p><p> rank(a) 求矩陣a的秩</p>
22、<p> trace(a) 求矩陣a的跡(對(duì)角線元素之和)</p><p> 我們?cè)谶M(jìn)行工程計(jì)算時(shí)常常遇到矩陣對(duì)應(yīng)元素之間的運(yùn)算。這種運(yùn)算不同于前面講的數(shù)學(xué)運(yùn)算,為有所區(qū)別,我們稱之為數(shù)組運(yùn)算。</p><p><b> 基本數(shù)學(xué)運(yùn)算</b></p><p> 數(shù)組的加、減與矩陣的
23、加、減運(yùn)算完全相同。而乘除法運(yùn)算有相當(dāng)大的區(qū)別,數(shù)組的乘除法是指兩同維數(shù)組對(duì)應(yīng)元素之間的乘除法,它們的運(yùn)算符為“.*”和“./”或“.\”。前面講過(guò)常數(shù)與矩陣的除法運(yùn)算中常數(shù)只能做除數(shù)。在數(shù)組運(yùn)算中有了“對(duì)應(yīng)關(guān)系”的規(guī)定,數(shù)組與常數(shù)之間的除法運(yùn)算沒(méi)有任何限制。</p><p> 另外,矩陣的數(shù)組運(yùn)算中還有冪運(yùn)算(運(yùn)算符為 .^ )、指數(shù)運(yùn)算(exp)、對(duì)數(shù)運(yùn)算(log)、和開(kāi)方運(yùn)算(sqrt)等。有了“對(duì)應(yīng)元素
24、”的規(guī)定,數(shù)組的運(yùn)算實(shí)質(zhì)上就是針對(duì)數(shù)組內(nèi)部的每個(gè)元素進(jìn)行的。矩陣的冪運(yùn)算與數(shù)組的冪運(yùn)算有很大的區(qū)別。</p><p> 邏輯關(guān)系運(yùn)算 </p><p> 邏輯運(yùn)算是MATLAB中數(shù)組運(yùn)算所特有的一種運(yùn)算形式,也是幾乎所有的高級(jí)語(yǔ)言普遍適用的一種運(yùn)算。</p><p><b> 5 程序及說(shuō)明</b></p><
25、p> 本程序采用MATLAB程序設(shè)計(jì)語(yǔ)言進(jìn)行程序的設(shè)計(jì)。MATLAB程序設(shè)計(jì)語(yǔ)言結(jié)構(gòu)完整,且具有優(yōu)良的移植性,它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。它可以高效率地解決工業(yè)計(jì)算問(wèn)題,特別是關(guān)于矩陣和矢量的計(jì)算。MATLAB與C語(yǔ)言和FORTRAN語(yǔ)言相比更容易被掌握。通過(guò)M語(yǔ)言,可以用類(lèi)似數(shù)學(xué)公式的方式來(lái)編寫(xiě)算法,大大降低了程序所需的難度并節(jié)省了時(shí)間,從而可把主要的精力集中在算法的構(gòu)思而不是編程上。MATLAB設(shè)計(jì)中,原始數(shù)據(jù)的填寫(xiě)
26、格式是很關(guān)鍵的一個(gè)環(huán)節(jié),它與程序使用的方便性和靈活性有著直接的關(guān)系。原始數(shù)據(jù)輸入格式的設(shè)計(jì),主要應(yīng)從使用的角度出發(fā),原則是簡(jiǎn)單明了,便于修改。</p><p><b> 5.1程序流程圖</b></p><p> 圖5.1 程序主流程圖</p><p> 圖5.2 迭代部分流程圖</p><p><b>
27、 5.2具體程序</b></p><p> %電力系統(tǒng)PQ分解法潮流計(jì)算</p><p> disp('電力系統(tǒng)極坐標(biāo)下的PQ分解法潮流計(jì)算:');</p><p><b> clear</b></p><p> n=input('請(qǐng)輸入結(jié)點(diǎn)數(shù):n=');</p
28、><p> n1=input('請(qǐng)輸入PV結(jié)點(diǎn)數(shù):n1=');</p><p> n2=input('請(qǐng)輸入PQ結(jié)點(diǎn)數(shù):n2=');</p><p> isb=input('請(qǐng)輸入平衡結(jié)點(diǎn):isb=');</p><p> pr=input('請(qǐng)輸入精確度:pr=');<
29、;/p><p> K=input('請(qǐng)輸入變比矩陣看:K=');</p><p> C=input('請(qǐng)輸入支路阻抗矩陣:C=');</p><p> y=input('請(qǐng)輸入支路導(dǎo)納矩陣:y=');</p><p> U=input('請(qǐng)輸入結(jié)點(diǎn)電壓矩陣:U=');<
30、/p><p> S=input('請(qǐng)輸入各結(jié)點(diǎn)的功率:S=');</p><p> Z=zeros(1,n);N=zeros(n2,n2+n1);L=zeros(n1+n2,n2);QT1=zeros(1,n1+n2);</p><p><b> for m=1:n</b></p><p><b&
31、gt; for R=1:n</b></p><p> C(m,m)=C(m,m)+y(m,R);</p><p> if K(m,R)~=0</p><p> C(m,m)=C(m,m)+1/((K(m,R)*C(m,R))/(K(m,R)-1));</p><p> C(R,R)=C(R,R)+1/((K(m,R)^2
32、*C(m,R))/(1-K(m,R)));</p><p> C(m,R)=C(m,R)*K(m,R);</p><p> C(R,m)=C(m,R);</p><p><b> end</b></p><p><b> end</b></p><p><b&g
33、t; end</b></p><p><b> for m=1:n</b></p><p><b> for R=1:n</b></p><p><b> if m~=R</b></p><p> Z(m)=Z(m)+1/C(m,R);</p>
34、<p><b> end </b></p><p><b> end</b></p><p><b> end</b></p><p><b> for m=1:n</b></p><p><b> for R=1:n&l
35、t;/b></p><p><b> if m==R</b></p><p> Y(m,m)=C(m,m)+Z(m);</p><p><b> else</b></p><p> Y(m,R)=-1/C(m,R);</p><p><b> end
36、</b></p><p><b> end</b></p><p><b> end</b></p><p> disp('結(jié)點(diǎn)導(dǎo)納矩陣:');</p><p><b> disp(Y);</b></p><p>
37、disp('迭代中關(guān)于B的矩陣:');</p><p> G=real(Y);</p><p> B=imag(Y);</p><p> O=angle(U);</p><p> U1=abs(U);</p><p><b> k=0;</b></p>&l
38、t;p><b> PR=1;</b></p><p> P=real(S);</p><p> Q=imag(S);</p><p> while PR>pr</p><p> for m=1:n2</p><p> UD(m)=U1(m);</p><
39、p><b> end</b></p><p> for m=1:n1+n2</p><p><b> for R=1:n</b></p><p> PT(R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));</p><p
40、><b> end</b></p><p> PT1(m)=sum(PT);</p><p> PP(m)=P(m)-PT1(m);</p><p> PP1(k+1,m)=PP(m);</p><p><b> end</b></p><p> for m
41、=1:n2</p><p><b> for R=1:n</b></p><p> QT(R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(R)));</p><p> end </p><p> QT1(m)=sum(QT);
42、 </p><p> QQ(m)=Q(m)-QT1(m); </p><p> QQ1(k+1,m)=QQ(m);</p><p><b> end</b></p><p> PR1=max(abs(PP));</p><p> PR2=max(abs(QQ));<
43、/p><p> PR=max(PR1,PR2);</p><p> for m=1:n1+n2</p><p> for R=1:n1+n2</p><p> B1(m,R)=B(m,R);</p><p><b> end</b></p><p><b>
44、 end</b></p><p> for m=1:n2</p><p> for R=1:n2</p><p> B2(m,R)=B(m,R);</p><p><b> end</b></p><p><b> end</b></p>
45、<p> JJ=[B1 L;N B2];</p><p><b> disp(JJ);</b></p><p> for m=1:n1+n2</p><p> PP2(m)=PP(m)/U(m);</p><p><b> end</b></p><p>
46、; for m=1:n2</p><p> QQ2(m)=QQ(m)/U(m);</p><p><b> end</b></p><p> PQ=[PP2';QQ2'];</p><p> DA=-inv(JJ)*PQ; </p><p><b> DA1
47、=DA';</b></p><p> for m=1:n1+n2</p><p> OO(m)=DA1(m)/U(m);</p><p><b> end</b></p><p> for m=n:n1+n2+n2</p><p> UU(m-n1-n2)=DA1(m
48、);</p><p><b> end</b></p><p> for m=1:n1+n2</p><p> O(m)=O(m)+OO(m);</p><p><b> end</b></p><p> for m=1:n2</p><p>
49、; U1(m)=U1(m)+UU(m);</p><p><b> end</b></p><p> for m=1:n1+n2</p><p> o(k+1,m)=180/pi*O(m);</p><p><b> end</b></p><p> for m=
50、1:n2</p><p> u(k+1,m)=U1(m);</p><p><b> end</b></p><p><b> k=k+1;</b></p><p><b> end</b></p><p><b> for m=1:
51、n</b></p><p> b(m)=U1(m)*cos(O(m));</p><p> c(m)=U1(m)*sin(O(m));</p><p><b> end</b></p><p><b> U=b+i*c;</b></p><p><b
52、> for R=1:n</b></p><p> PH1(R)=U(isb)*conj(Y(isb,R))*conj(U(R));</p><p><b> end</b></p><p> PH=sum(PH1);</p><p><b> for m=1:n</b>&
53、lt;/p><p><b> for R=1:n</b></p><p><b> if m~=R</b></p><p> C1(m,R)=1/C(m,R);</p><p><b> else</b></p><p> C1(m,m)=C(m,
54、m);</p><p><b> end</b></p><p><b> end</b></p><p><b> end</b></p><p><b> for m=1:n</b></p><p><b>
55、 for R=1:n</b></p><p> if (C(m,R)~=inf)&(m~=R)</p><p> SS(m,R)=U1(m)^2*conj(C1(m,m))+U(m)*(conj(U(m))-conj(U(R)))*conj(C1(m,R));</p><p><b> end</b></p>
56、;<p><b> end</b></p><p><b> end</b></p><p> disp('迭代中的△P:');disp(PP1);</p><p> disp('迭代中的△Q:');disp(QQ1);</p><p> di
57、sp('迭代中相角:');disp(o);</p><p> disp('迭代中電壓的模:');disp(u);</p><p> disp('平衡結(jié)點(diǎn)的功率:');disp(PH);</p><p> disp('全部線路功率分布:');disp(SS);</p><p>
58、; 5.3 程序計(jì)算結(jié)果</p><p><b> 結(jié)點(diǎn)導(dǎo)納矩陣:</b></p><p> 0.8381 - 3.8179i -0.4044 + 1.6203i 0 0 -0.4337 + 2.2586i</p><p> -0.4044 + 1.6203i 0.7769
59、 - 3.4370i -0.3726 + 1.8557i 0 0 </p><p> 0 -0.3726 + 1.8557i 1.1608 - 7.0930i -0.5224 + 4.1792i -0.2739 + 1.2670i</p><p> 0 0
60、 -0.5224 + 4.1792i 0.5499 - 4.3991i 0 </p><p> -0.4337 + 2.2586i 0 -0.2739 + 1.2670i 0 0.7077 - 3.5257i</p><p> 迭代中關(guān)于B的矩陣:</p><p>
61、; -3.8179 1.6203 0 0 0 0 0</p><p> 1.6203 -3.4370 1.8557 0 0 0 0</p><p> 0 1.8557 -7.0930 4.1792
62、 0 0 0</p><p> 0 0 4.1792 -4.3991 0 0 0</p><p> 0 0 0 0 -3.8179 1.6203 0</p><p> 0
63、 0 0 0 1.6203 -3.4370 1.8557</p><p> 0 0 0 0 0 1.8557 -7.0930</p><p> -3.8179 1.6203 0 0 0
64、 0 0</p><p> 1.6203 -3.4370 1.8557 0 0 0 0</p><p> 0 1.8557 -7.0930 4.1792 0 0 0</p><p> 0
65、 0 4.1792 -4.3991 0 0 0</p><p> 0 0 0 0 -3.8179 1.6203 0</p><p> 0 0 0 0 1.6203 -3.4370
66、 1.8557</p><p> 0 0 0 0 0 1.8557 -7.0930</p><p> -3.8179 1.6203 0 0 0 0 0</p><p> 1.6203 -
67、3.4370 1.8557 0 0 0 0</p><p> 0 1.8557 -7.0930 4.1792 0 0 0</p><p> 0 0 4.1792 -4.3991 0 0
68、 0</p><p> 0 0 0 0 -3.8179 1.6203 0</p><p> 0 0 0 0 1.6203 -3.4370 1.8557</p><p> 0 0
69、 0 0 0 1.8557 -7.0930</p><p> -3.8179 1.6203 0 0 0 0 0</p><p> 1.6203 -3.4370 1.8557 0 0 0
70、 0</p><p> 0 1.8557 -7.0930 4.1792 0 0 0</p><p> 0 0 4.1792 -4.3991 0 0 0</p><p> 0 0 0
71、 0 -3.8179 1.6203 0</p><p> 0 0 0 0 1.6203 -3.4370 1.8557</p><p> 0 0 0 0 0 1.8557 -7.0930</p><p&
72、gt; -3.8179 1.6203 0 0 0 0 0</p><p> 1.6203 -3.4370 1.8557 0 0 0 0</p><p> 0 1.8557 -7.0930 4.1792
73、 0 0 0</p><p> 0 0 4.1792 -4.3991 0 0 0</p><p> 0 0 0 0 -3.8179 1.6203 0</p><p> 0
74、 0 0 0 1.6203 -3.4370 1.8557</p><p> 0 0 0 0 0 1.8557 -7.0930</p><p> -3.8179 1.6203 0 0 0
75、 0 0</p><p> 1.6203 -3.4370 1.8557 0 0 0 0</p><p> 0 1.8557 -7.0930 4.1792 0 0 0</p><p> 0 0
76、 4.1792 -4.3991 0 0 0</p><p> 0 0 0 0 -3.8179 1.6203 0</p><p> 0 0 0 0 1.6203 -3.4370 1.8557<
77、;/p><p> 0 0 0 0 0 1.8557 -7.0930</p><p> -3.8179 1.6203 0 0 0 0 0</p><p> 1.6203 -3.4370 1.8557
78、 0 0 0 0</p><p> 0 1.8557 -7.0930 4.1792 0 0 0</p><p> 0 0 4.1792 -4.3991 0 0 0</p>
79、<p> 0 0 0 0 -3.8179 1.6203 0</p><p> 0 0 0 0 1.6203 -3.4370 1.8557</p><p> 0 0 0 0 0
80、 1.8557 -7.0930</p><p><b> 迭代中的△P:</b></p><p> -0.1983 -0.1800 -0.2482 0.3225</p><p> 0.0104 0.0183 -0.0040 0.0055</p><p> 0.0037
81、-0.0004 0.0039 -0.0012</p><p> -0.0000 -0.0004 -0.0004 0.0004</p><p> -0.0001 0.0001 -0.0001 0.0000</p><p> -0.0000 0.0000 0.0000 -0.0000</p>&
82、lt;p> 0.0000 -0.0000 0.0000 0.0000</p><p><b> 迭代中的△Q:</b></p><p> 0.0339 -0.0510 0.1423</p><p> -0.0536 -0.0461 -0.0490</p><p> 0.0
83、002 0.0026 -0.0006</p><p> 0.0012 0.0001 0.0011</p><p> 0.0001 -0.0001 -0.0001</p><p> -0.0000 0.0000 -0.0000</p><p> -0.0000 0.0000 0.000
84、0</p><p><b> 迭代中相角:</b></p><p> -6.3664 -7.9886 -3.6795 0.7049</p><p> -5.8635 -7.1725 -3.1710 1.2602</p><p> -5.7670 -7.0765 -3.0640
85、 1.3463</p><p> -5.7734 -7.0905 -3.0725 1.3432</p><p> -5.7753 -7.0903 -3.0735 1.3427</p><p> -5.7754 -7.0898 -3.0731 1.3430</p><p> -5.7753
86、 -7.0899 -3.0731 1.3430</p><p><b> 迭代中電壓的模:</b></p><p> 1.0090 1.0003 1.0201</p><p> 0.9797 0.9642 1.0038</p><p> 0.9802 0.9653
87、1.0040</p><p> 0.9807 0.9657 1.0043</p><p> 0.9807 0.9657 1.0042</p><p> 0.9807 0.9657 1.0042</p><p> 0.9807 0.9657 1.0042</p><p&
88、gt;<b> 平衡結(jié)點(diǎn)的功率:</b></p><p> 0.3531 + 0.1786i</p><p><b> 全部線路功率分布</b></p><p> 0 0.0413 - 0.0432i 0 0 -0.2613 -
89、 0.1555i</p><p> -0.0410 - 0.0506i 0 -0.1390 - 0.0757i 0 0 </p><p> 0 0.1332 - 0.1598i 0 -0.3276 - 0.1401i -0.09
90、20 - 0.2516i</p><p> 0 0 0.3500 + 0.1743i 0 0 </p><p> 0.2679 + 0.1313i 0 0.0852 + 0.0473i 0
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 眾賞文庫(kù)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)
- 電力系統(tǒng)課程設(shè)計(jì)-- 潮流計(jì)算
- 電力系統(tǒng)課程設(shè)計(jì)-----潮流計(jì)算
- 電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)
- 電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)
- 電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)
- 電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)
- 電力系統(tǒng)課程設(shè)計(jì)--潮流計(jì)算
- 課程設(shè)計(jì)--電力系統(tǒng)設(shè)計(jì)潮流計(jì)算
- 電力系統(tǒng)穩(wěn)態(tài)潮流計(jì)算課程設(shè)計(jì)
- 電力系統(tǒng)穩(wěn)態(tài)課程設(shè)計(jì)--潮流計(jì)算
- 電力系統(tǒng)潮流短路計(jì)算課程設(shè)計(jì)
- 電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)報(bào)告
- 電力系統(tǒng)分析課程設(shè)計(jì)-- 電力系統(tǒng)潮流計(jì)算
- 電力系統(tǒng)分析課程設(shè)計(jì)--電力系統(tǒng)潮流計(jì)算
- 潮流計(jì)算課程設(shè)計(jì)--基于matlab的電力系統(tǒng)潮流計(jì)算
- 用matlab電力系統(tǒng)潮流計(jì)算(課程設(shè)計(jì))
- 電力系統(tǒng)綜合課程設(shè)計(jì)--電力系統(tǒng)的潮流計(jì)算和故障分析
- 電力系統(tǒng)分析課程設(shè)計(jì)---潮流計(jì)算
- 基于matlab的電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)
評(píng)論
0/150
提交評(píng)論