分塊矩陣的三角分解與局部重新分解

分塊矩陣的三角分解與局部重新分解

ID:36285556

大?。?73.00 KB

頁數(shù):14頁

時(shí)間:2019-05-08

分塊矩陣的三角分解與局部重新分解_第1頁
分塊矩陣的三角分解與局部重新分解_第2頁
分塊矩陣的三角分解與局部重新分解_第3頁
分塊矩陣的三角分解與局部重新分解_第4頁
分塊矩陣的三角分解與局部重新分解_第5頁
資源描述:

《分塊矩陣的三角分解與局部重新分解》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。

1、分塊矩陣的三角分解與局部重新分解及其在電力系統(tǒng)中的應(yīng)用一、引言眾所周知,矩陣的三角分解法是一種非常有效的手段,已廣泛應(yīng)用于電力系統(tǒng)應(yīng)用程序,諸如:潮流(牛頓-拉夫遜法)、快速解藕潮流、狀態(tài)估計(jì)、安全分析、最優(yōu)潮流、暫態(tài)穩(wěn)定等。本文對(duì)分塊矩陣的三角分解與局部重新分解列出了幾種算法,作了論證和分析比較,并給出了在電力系統(tǒng)中的應(yīng)用二、分塊矩陣的分解方法方法1、按行列逐一形成法對(duì)于非奇異n階方陣Α,可分割成分塊矩陣:式中,Α12,,和Α22分別為n1階和n2階方陣,n1+n2=n,Α11又為非奇異陣,L1和U1分別為下三角陣和上三角陣。分塊矩陣的劃分由一組正整數(shù)n1,n2,…,nq確定,且n1+n2+

2、…+nq=n。將Α寫作?。?),并根據(jù)式(1)記:則又可將?。?)分割成分塊矩陣:式中,和分別為n1階和n1-n2階方陣,為非奇異陣。再記逐次作類似上述處理,矩陣Α便得(推導(dǎo)略):式中,子矩陣都是n1階方陣,q=n/n1是正整數(shù)。若n1,n2,…,nq不盡相等,也可作類似處理。若Α為對(duì)稱矩陣,則U=。若Α為關(guān)聯(lián)對(duì)稱矩陣,則U≠,但稀疏結(jié)構(gòu)相同。方法2、直接計(jì)算法將Α、L、D和U寫成展開形式,就不難證明可由下列諸式來表述式(5)(如果,k=1,2,…,q都非奇異),因而可藉以直接計(jì)算。式中,k=1,2,…,q,分塊矩陣的劃分由一組正整數(shù)n1,n2,…,nq確定,且n1+n2+…+nq=n。在電力

3、系統(tǒng)應(yīng)用程序中,矩陣Α是稀疏的,子矩陣、有很多是零陣,作LDU分解后所得的三角矩陣也有高度的稀疏性。在很多程序中,n1=n2=…=nq=2。方法3、簡(jiǎn)單局部重新分解法只重新分解矩陣元素變化了的子矩陣,如圖1的陰影部分。即從矩陣元素開始有變化的k值起用式(9)-(12)計(jì)算一遍。因此,要更新的子矩陣的大小以及重新分解的工作量取決于變化了的矩陣元素所處的位置,如接近矩陣頂部,則重新分解的工作量節(jié)約甚微以至沒有。方法4、特殊排序的簡(jiǎn)單局部重新分解法通過特殊排序,將變化了的矩陣元素集中到矩陣右下部,從而可減少重新分解的工作量。但有可能嚴(yán)重地降低了矩陣的稀疏性,從而又增加計(jì)算量。方法5、重新分解有關(guān)行列

4、圖2示電力系統(tǒng)牛頓一拉夫遜法潮流計(jì)算的有關(guān)情況。圖2a為單線結(jié)線圖。圖2b為導(dǎo)納矩陣結(jié)構(gòu),其中黑點(diǎn)表示非零元素。圖2c為下三角矩陣結(jié)構(gòu),其中黑點(diǎn)表示與圖2b對(duì)應(yīng)的非零元素,圓圈表示三角分解時(shí)產(chǎn)生的非零“注入”元素。由電力系統(tǒng)理論可知,雅可比矩陣的稀疏性與系統(tǒng)導(dǎo)納矩陣是一樣的,雅可比矩陣常以2階方陣作為子矩陣(元素)以構(gòu)成分塊矩陣。由式(9)可見,k=7,i=5,雖然=0,但≠0,因而≠0。設(shè)矩陣Α已作LDU分解,現(xiàn)在、變化了(設(shè)電力系統(tǒng)節(jié)點(diǎn)3與5之間有雙回線,一回?cái)嚅_),則:按方法3,重新分解第3一8行列;按方法4,將節(jié)點(diǎn)號(hào)3、5改成大的節(jié)點(diǎn)號(hào)(例如7、8),然后重新分解;按本法(方法5),

5、只需重新分解第3和第5一8行列.因?yàn)楦鶕?jù)式(9)-(11)可知,隨著和的改變,、、、、、、、、、、、、將相應(yīng)變化。重新分解的行列號(hào)可由補(bǔ)路徑圖(圖3)得到。圖28節(jié)點(diǎn)的網(wǎng)絡(luò)(a)實(shí)際系統(tǒng)圖;(b)導(dǎo)納矩陣結(jié)構(gòu);(c)下三角矩陣結(jié)圖38節(jié)點(diǎn)網(wǎng)絡(luò)的路徑圖一、在電力系統(tǒng)中的應(yīng)用舉例1、牛頓一拉夫遜法潮流計(jì)算潮流計(jì)算是應(yīng)用極其廣泛的,它是在向量F和雅可比矩陣J已知的情況下,用三角分解法求解下列方程中的向量ΔX。如果J只有部分以至少量子塊發(fā)生變化,則可作局部重新分解。本文建議,對(duì)#α節(jié)點(diǎn)和#3節(jié)點(diǎn)各考慮兩個(gè)方程,它們的雅可比矩陣的子塊分別為:式中,i、j是節(jié)點(diǎn)號(hào);P、Q分別為有功功率和無功功率;e、f分

6、別是電壓的實(shí)部和虛部;V是電壓的模。這樣,在節(jié)點(diǎn)類別變換時(shí)不致發(fā)生Α的結(jié)構(gòu)性變更。1、穩(wěn)態(tài)狀態(tài)估計(jì)加權(quán)最小二乘法穩(wěn)態(tài)狀態(tài)估計(jì)的基本關(guān)系式:信息陣Α和向量B已知,求向量ΔX。A劃分為2?2階的子矩陣:式中,i、j、e、f的含義與式(44)中相同,Zk為第k個(gè)量測(cè),m為測(cè)量總數(shù)。為Zk的方差。也可應(yīng)用分塊矩陣的三角分解或局部重新分解。二、結(jié)論很多問題應(yīng)用分塊矩陣的思想來處理是有益的。本文對(duì)分塊矩陣的三角分解與局部重新分解作了研究,導(dǎo)出了幾種算法,作了證明,并進(jìn)行了分析比較。在電力系統(tǒng)應(yīng)用程序中,分塊矩陣的三角分解與局部重新分解是很有效的。總黃酮生物總黃酮是指黃酮類化合物,是一大類天然產(chǎn)物,廣泛存在

7、于植物界,是許多中草藥的有效成分。在自然界中最常見的是黃酮和黃酮醇,其它包括雙氫黃(醇)、異黃酮、雙黃酮、黃烷醇、查爾酮、橙酮、花色苷及新黃酮類等。簡(jiǎn)介  近年來,由于自由基生命科學(xué)的進(jìn)展,使具有很強(qiáng)的抗氧化和消除自由基作用的類黃酮受到空前的重視。類黃酮參與了磷酸與花生四烯酸的代謝、蛋白質(zhì)的磷酸化、鈣離子的轉(zhuǎn)移、自由基的清除、抗氧化活力的增強(qiáng)、氧化還原作用、螯合作用和基因的表達(dá)。它們對(duì)健康的好處有

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動(dòng)畫的文件,查看預(yù)覽時(shí)可能會(huì)顯示錯(cuò)亂或異常,文件下載后無此問題,請(qǐng)放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫負(fù)責(zé)整理代發(fā)布。如果您對(duì)本文檔版權(quán)有爭(zhēng)議請(qǐng)及時(shí)聯(lián)系客服。
3. 下載前請(qǐng)仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時(shí)可能由于網(wǎng)絡(luò)波動(dòng)等原因無法下載或下載錯(cuò)誤,付費(fèi)完成后未能成功下載的用戶請(qǐng)聯(lián)系客服處理。