第一篇：潮流音樂

潮流音樂：

某人-汪蘇瀧、浪漫冬季-小峰峰、風度-汪蘇瀧、單色冰激凌-單色凌、雨停曲-阿悄、停止跳動-汪蘇瀧、失她-萬軒嚴、顛覆-星弟、oh baby-孫斌、當你聽我說-胡夏、默戀-Evin。柒、那首歌，很好聽-孫斌、江湖有事-汪蘇瀧、友情演出-徐良，小凌、三國殺-汪蘇瀧、自怨自愛-貫施欽、不像情歌-小賤、試探-孫斌、為你脫掉-單色凌、最心疼的玩笑-KingStar、不止心傷-萬軒嚴、累不累-汪蘇瀧、任性愛-阿悄、洪湖岸邊-喬洋、當賣火柴的小女孩遇上放煙花的小男孩-阿悄,小賤、我也不知道-汪蘇瀧、后會無期-徐良、傷以成曲-五華、沒有這首歌-后弦、分手做朋友-星弟、壞女孩-徐良、美瞳-徐良、不寫情歌-徐良、飛機場-徐良、天龍八部之宿敵-許嵩、.請答應我的獨白-小峰峰、不良-徐良、情話毀滅-五華、我的女孩你別碰-KingStar、逆戰(zhàn)-張杰、分手才快樂-弦小杰、雨天是放聲哭泣的時間-秀愛組合、光年-Xun、繁華過后-貫詩欽、每天想我一點點-KingStar、壞人-小米、分手的距離-本兮、你聽了嗎，那首歌-孫斌、雪花謠-許嵩、把煩惱丟下去-單色凌、請安靜的忘記我-孫子涵、依然愛你-王力宏、致命的甜蜜-孫子涵、對待壞女孩-孫子涵、.全世界宣布我愛你-孫子涵、老男孩-劉忻、危險走紅-單色凌、天空物語-KingStar、電話里的秘密-徐良

蝴蝶迷咒-漁圈、試探你的溫柔-單色凌、蝶變-單色凌、本來這個時候-小賤、如果不能在一起-曉曉文、小背叛-KingStar、有沒有-by2、分不清的傷-單色凌、觸碰純白-單色凌、告訴我-阿悄、歸零-萬軒嚴、自愛-阿悄、考試什么的都去死吧-徐良、哩哩哩-小賤、舍子花-阿悄、les不哭-莫小娘、哭泣的格桑花-莫小娘、說好了不見面-小賤、當我唱起這首歌-小賤、未成年-本兮、怎么辦我愛你-本兮、慣性想象-（本兮&單色凌）、奇怪，我不懂得愛-本兮、你在看孤獨的風景-（本兮&單小緣）、當我們一起走過-蘇打綠、我想我不夠好-單色凌、泡泡泡泡-本兮、不想和你做朋友-（星弟&小賤）、殘留的回憶-小

5、再見了笑容-小

5、這個冬天我還在-小峰峰、對不起，沒關系-小峰峰、兩個傻瓜-徐良、假裝開心-小峰峰、傷之戀-小峰峰、很有愛-后弦、信條-后弦、喜歡寂寞-蘇打綠、幸福額度-蘇打綠、南泥灣-單小源、黑色眼淚-灰原窮、雨夜-灰原窮、小愛-貫詩欽、你不是我的，我不是你的-貫詩欽、一樣愛著你-by2 拜拜，愛過-SARA、能不能-mxd小寶、關于我她不懂-灰原窮、藤-漁圈...好聽歌曲：

《因為了解》《甩蔥歌》《巴赫舊約》《不配做你的男朋友》《莫斯科的天空》《幸福是被你需要》《為你唱首歌》《離開以后》《放手也是幸?！贰段覑畚易约骸贰稄涂袒貞洝贰督釿ao》《苦笑》《你還愛著他》《不同》《離開以后》《獨自旅行》《素顏》《never bereplaced remix》《你的爛借口》《天真》《曾經(jīng)、我們的愛》《你不像她》《蘇州城外的微笑》《浪漫式乖乖》《假裝淡定》《全世界只有你和我》

《相愛和分開都是罪》《銘記我們的愛情》《JOE》《停在昨天》《最后是我開了口》《欠你的太多》《配角》《始終不夠》《《專屬味道》《專屬電影》《淚落在琴上》《沒有忘掉他》《犯賤》《一點點的藍》《小星星》《別讓我們的愛變成回憶》《大女孩》《牽牽》《不必了》《我太任性》《過去式》《好安靜》《放晴》《最后的風度》《你還不懂》《我的心

好冷》《然后沒有然后》《娃娃臉》《朋友關系》《默寫》《你還裝什么》《慢慢習慣》《寂寞作祟》《假動作》

《 幸福的記號》《我們一起唱》《像笑話般的愛你》《草戒指》《犯賤》《未見面就說再見》《愛情劇終》《一千遍我愛你》《習慣一個人》《習慣》 《你讓我懂》《月光》《最后，是我開了口》《咬耳朵》《微博控》《我總是一個人在練習》《愛很美》《認輸》《故意忘記你》《專屬感動》《死一樣的痛過》《假如愛是一種錯覺》《我不是沒臉的男孩》《不分手的戀愛》

《唯你懂我心》《放不下》《萬一奧特曼打不贏小怪獸》《放棄》《相約在雨季》《我會很誠實》《我們的回憶》《止步愛情》《不舍》《你還欠我一個擁抱》《最美的風景》《失眠》《好想對你說》》《卑微的承諾》《回不去的甜》《hi小姐》《劇終》《若相惜》《回憶唱給你聽》《不要走太遠》《如果等待只能成為等待》《要命的煩惱》《始終不夠》《分開后的你》《傷之戀》《潛意識失憶》《迷失終點》《苦笑》《結局后才明白》 《卑戀》《因為你喜歡聽》《寂寞在唱什么歌》

《難以啟齒的柔弱》《我們的愛在贖罪》《說一句我不走了》《最后一刻才明白》《多情的人總被無情的傷》《我們一起長大》《戀上你的無賴》《超級自戀男》《掩飾寂寞》《偶爾是悲傷》《距離》《浪漫式愛情》《我的眼淚是你看不見的傷心》《寂夜》《無淚的遺憾》《放棄后的心疼》《放手》《不要愛情了》《末日》《主動》《笑著說分手》《等不到你》《被遺忘的愛》《心酸的承諾》《不分開好不好》《說什么我愛你》《偏偏懷念》《如果你還不明白》《唯你懂我心》《和平分手》《18歲的紀念》《蒼 白》《你真的走了》

第二篇：潮流報告

matlab潮流計算報告

王振 電氣3班

引言

隨著我國電力工業(yè)的迅猛發(fā)展，集中抄表系統(tǒng)得到了廣泛的應用。集中抄表系統(tǒng)是由主站通過傳輸媒體(無線、有線、電力線載波等信道或IC 卡等介質(zhì))集中抄讀多個電能表電能量記錄值的自動化系統(tǒng)。它主要由采集用戶電能表電能量信息的采集終端(或采集模塊)、集中器、信道和主站等設備組成。它可以提供配電網(wǎng)終端實時的功率、電壓數(shù)據(jù)。充分利用這些測量數(shù)據(jù)，可以把傳統(tǒng)潮流計算的非線性方程轉(zhuǎn)化為線性方程求解。提出應用PQ 法、牛頓-拉夫遜法進行潮流計算，但由于配電網(wǎng)線路R/X 較大，潮流的收斂性難于保證。結合配電網(wǎng)輻射狀的特點，以支路電流或母線電壓為研究對象建立運算模型。這些算法通過求解功率、電壓的偏移量進行迭代求解，具有很高的執(zhí)行效率，但是當量測量增多時，需要進行多次迭代計算，計算速度不高。本文應用集中抄表系統(tǒng)采集的數(shù)據(jù)，當所有采集點電壓、功率可準確量測時，通過矩陣變換消去中間變量，給出一種快速線性潮流計算方法；當少量節(jié)點數(shù)據(jù)未知時，也給出一種線性計算方法。并利用30節(jié)點配電網(wǎng)算例驗證了所提算法的正確性及有效性。

潮流計算是電力系統(tǒng)規(guī)劃、運行的基本研究方法。隨著現(xiàn)代電力系統(tǒng)大系統(tǒng)、強非線性與多元件的特點Et益突出，其計算量與計算復雜度急劇增加。舊的計算機軟件在處理潮流計算時，其速度已無法滿足大電網(wǎng)模擬和實時控制的仿真要求，而高效的潮流問題的相關軟件的研究已成為大規(guī)模電力系統(tǒng)仿真計算的關鍵。隨著計算機技術的不斷發(fā)展和成熟，基于MATLAB潮流計算研究近年來得到了長足的發(fā)展，為真正解決大電網(wǎng)快速、詳細的仿真技術開辟了新思路。針對這一現(xiàn)狀，以某電力網(wǎng)絡為例，分析了BASIC、FORTRAN 和MATLAB高級語言潮流計算的異同，指出它們的優(yōu)缺點，并針對潮流計算模型結構的特點，提出了基于MAT—LAB的潮流算法.問題分析

l 潮流計算的數(shù)學模型

電力系統(tǒng)潮流計算的基本方程: ?YijUj?j?1n.Pi?jQiUi*Pi?jQi.,(1,2,3,...,n)

?ZijIj?j?1n.Ui*.,(1,2,3,...,n)(1)式(1)中。各有，z個非線性復數(shù)方程，對它作不同的應用和處理，就形成了不同的潮流計算方法。

Newton I aphson法，收斂性好，是非線性方程數(shù)值求解的有效方法。它把非線性方程線性化，線性方程的系數(shù)矩陣在結構上是稀疏、非對稱矩陣，結合稀疏矩陣技術，計算機內(nèi)存占用量大大減少，計算速度也比以往的方法優(yōu)越；P—Q分解法是在Newton—Laphson法的基礎上將有功功率P和無功功率Q分開交替迭代的潮流計算方法。方法計算過程簡單，所占計算機內(nèi)存約為Newton—Laphson法的1／2～1／4，計算速度也顯著加快，是目前常用的潮流計算方法。由于近代電

力系統(tǒng)網(wǎng)絡節(jié)點數(shù)量極大，節(jié)點導納陣y 的稀疏度也極高，稀疏技術的形成是必不可少的。為了減少潮流計算程序中求逆、求積運算所耗費的大量機時，用三角分解法形成因子表解算修正方程已成為現(xiàn)代電力系統(tǒng)潮流計算的主要方法。2 高級語言的類型

計算機潮流計算的基本要求是計算方法具有一定的可靠性和收斂性；盡量選用占計算機內(nèi)存較少的存儲方式；在可靠收斂的前提下，選用計算速度較快的方法；人機交互環(huán)境便于數(shù)據(jù)輸入、校核和修改，且具有一定的靈活性。BASIC是一種解決數(shù)學問題的語言，其取消編譯工作及連接過程，語言的執(zhí)行時間相應較慢。BA—SIC語句的標識符可為常數(shù)和變量，矩陣也可規(guī)定為變量進行計算。

FORTRAN是1957年發(fā)明的用于科學計算的語言，提供實型、雙精度型和復數(shù)型等數(shù)據(jù)結構，其雙精度型數(shù)據(jù)結構使所求得的解保留更多的有效位數(shù)，保證更高的精度；復數(shù)的引人為復數(shù)計算提供了方便。由于其在數(shù)學運算上的廣泛性，也曾是潮流計算的一種主要計算機語言，但在矩陣的計算上其程序設計較為繁瑣。MATI AB是集數(shù)值計算、符號運算及圖形處理等強大功能于一體的科學計算語言?，F(xiàn)已成工程計算中普遍采用的工具。MATLAB擁有600多個工程數(shù)學運算函數(shù)，可實現(xiàn)潮流計算中的矩陣求積、求逆、稀疏矩陣形成、復數(shù)運算以及初等數(shù)學運算等。MATI AB語言允許用戶以數(shù)學形式的語言編寫程序，比BASIC語言和FORTRAN語言等更為接近書寫的數(shù)學表達格式，其程序易調(diào)試。在計算要求相同的情況下，使用MATIAB編程，工作量將會大為減少。

圖(1)面向?qū)ο笮蛄谢夹g

3.配電網(wǎng)絡潮流牛頓計算方法：

在采用直角坐標系的前提條件F。

基于午頓法的潮流計算可由下列方程描述：

?J?V??W

式中，．，為雅可比矩陣，AW 為節(jié)點注入功率殘差向量。有關表達式如下：

則雅克比矩陣可以寫成：

上式中：Y為節(jié)點導納矩陣。在配電網(wǎng)中，由于支路較短，各支路兩端電壓落差很小，而且一般情況下不考慮對地電納，因此有

因此，在這種情況下，可以利用矩陣求逆運算的松馳方法近似計算

這樣，我們可以得出

圖2 多環(huán)配電網(wǎng)的重構流程

Fig.2 Reconfiguration flow chart of multi-loop distribution network 4.基于MATLAB電力系統(tǒng)潮流計算

試用P—Q分解法計算電力系統(tǒng)的潮流分布，分析BASIC、FORTRAN和MATLAB高級語言潮流計算特點。

(1)形成系數(shù)矩陣B、B(設該矩陣為4×4矩陣)。

?a11?a?21?.B??B?????.?.???a41

a12a22...a42....a14?....a24??....??

....?....??...a44??

在程序設計上，用MATI AB語言編寫程序形成系數(shù)矩陣B、B”遠比BASIC、FORTRAN簡單，矩陣輸入、輸出與數(shù)學書寫格式相似。在命令窗口輸入：

?a11?a?21?a31> B=??.?.???a41a12a22a32..a42...a14?...a24??....a34??

.....?.....??....a44??再在命令窗口輸入：>>B一

窗口將顯示出B、B”矩陣。

在BASIC、FORTRAN 語言編程時，必須進行矩陣說明，BASIC用DIM B(4，4)語句說明矩陣，F(xiàn)ORTRAN 用DOUBLE PRECISION B(1：4，1：4)語句同時說明矩陣和矩陣精度；在建立矩陣時，BASIC、FORTRAN語言用READ語句和DA—TA語句輸入矩陣；FORTRAN語言用WRITE語句輸出矩陣，而BASIC語言只能用二重循環(huán)語句輸出二維矩陣，相當繁瑣。在數(shù)據(jù)精度方面，MAT—LAB所有數(shù)據(jù)以雙精度類型進行數(shù)據(jù)的存儲和運算，其數(shù)值有效位數(shù)可達16～17位，范圍為10 ∞～ 10 ；FORTRAN 也提供雙精度類型，但編程中必須用DOUBI E語句說明，而其數(shù)據(jù)有效范圍10?38?1038.BASIC不提供雙精度類型，其數(shù)值只保留7位有效位數(shù)，有效范圍為一1．70141E一38～ 1．70141E+38 引。要達到題目所要求的計算精度，且兼顧矩陣程序設計的難易程度，MATLAB則成為首選的潮流計算的計算機語言。對于大矩陣的潮流計算，MATLAB提供的M 文件可將輸入矩陣按格式先寫入一個文本文件中，在編程時，按文本文件名在命令窗調(diào)用，大矩陣就被輸入到內(nèi)存中了。MATLAB也直接能創(chuàng)建復數(shù)矩陣，這兩點也是BASIC、FORTRAN語言不可比擬的。(2)求系數(shù)矩陣B、B 的逆陣，計算各節(jié)點電壓的相位角△。相位角△ 計算公式：

?(B?)?1(?P(K)/U(K))?(U(K)??i)

式中：?p(k).U(K)為已知量，相位角??i計算主要進行矩陣的求積運算。MATI AB和BASIC語言都用命令inv(B)或INV(B)來實現(xiàn)矩陣的求逆功能，MATI AB和BA—SIC語言求積運算均與數(shù)學書寫格式一致，但BA—SIC語言必須對矩陣進行MAT說明，而且執(zhí)行時間相應較慢，計算精度低。FORTRAN語言雖然在參考資料上提供矩陣求逆、求積程序，但它按線性代數(shù)上矩陣求逆、求積步驟編程l=1，編寫程序至少需用一個三重循環(huán)語句，所以與MATI AB相比，工作量很大，效率很低。

(3)計算平衡節(jié)點功率和線路功率。平衡節(jié)點功率為

S?U?YUt*?P1?Q1

*tit?1~.n.(k)線路功率為

Sij?UiIij?Pij?Qij

平衡節(jié)點功率和線路功率的計算屬復數(shù)運算。MATI AB和FORTRAN語言都提供復數(shù)功能句,MATLAB以數(shù)學上的復數(shù)書寫格式編寫程序方便，而且不容易出錯；而FORTRAN 必須對復數(shù)變量進行CMPI X語句的復數(shù)類型說明，較為繁瑣。BASIC不提供復數(shù)功能語句，一般將實數(shù)、虛數(shù)分開計算，最后用輸出語句寫在一起，所以一般不用BASIC語言編寫復數(shù)運算程序。

另外，MATI AB提供潮流計算稀疏矩陣的建立命令spconvert，可將外部數(shù)據(jù)轉(zhuǎn)化為稀疏矩陣，而且MATI AB 函數(shù)lu也可直接實現(xiàn)LR三角分解l1，從而為現(xiàn)代電力系統(tǒng)潮流計算提供了方便。

~..模型求解

利用matlab軟件，根據(jù)潮流計算先進行賦值，利用下面的語句：

%本程序的功能是用牛頓-拉夫遜法進行潮流計算 n=input('請輸入節(jié)點數(shù):n=');nl=input('請輸入支路數(shù):nl=');isb=input('請輸入平衡母線節(jié)電號:isb=');pr=input('請輸入誤差精度:pr=');B1=input('請輸入由支路參數(shù)形成的矩陣:B1=');%變壓器側為1，否則為0 B2=input('請輸入各節(jié)點參數(shù)形成的矩陣:B2=');X=input('請輸入由節(jié)點號及其對地阻抗形成的矩陣:X=');求出?P與?Q利用一下程序得到： P=real(S);Q=imag(S);ICT1=0;IT2=1;while IT2~=0 IT2=0;t1=1;t2=1;for i=1:n if i~=isb C(i)=0;D(i)=0;for j1=1:n

C(i)=C(i)+U(i)*U(j1)*(G(i,j1)*cos(cta(i)-cta(j1))+B(i,j1)*sin(cta(i)-cta(j1)));

D(i)=D(i)+U(i)*U(j1)*(G(i,j1)*sin(cta(i)-cta(j1))-B(i,j1)*cos(cta(i)-cta(j1)));end DP(t1)=P(i)-C(i);t1=t1+1;if B2(i,6)==2 DQ(t2)=Q(i)-D(i);t2=t2+1;end

end end t1=t1-1;t2=t2-1;DPQ=[DP';DQ'];得到相應的值：

心得體會

把潮流計算的任務完成后，有一些體會，主要體現(xiàn)在以下三個方面：

（1）、對Matlab這個工具本身有了一些了解，包括它的集成開發(fā)環(huán)境，調(diào)試運行方式，程序的基本語句以及控制方式。其中感受最深刻的是Matlab的數(shù)據(jù)結構，很多書上介紹說Matlab的數(shù)據(jù)單元是矩陣，但自己看了以后并沒有什么體會，在編寫調(diào)試程序時，出于C語言的習慣，我們把節(jié)點的初始電壓，節(jié)點的輸出功率都設置成單個數(shù)據(jù)的形式，但在寫功率平衡方程用到For語句時，單個數(shù)據(jù)就很不方便，于是又改成了矩陣形式，在后續(xù)的程序中，也均使用矩陣做為數(shù)據(jù)單元。通過編寫Matlab程序，對“Matlab的基本數(shù)據(jù)單元是矩陣”這句話有了一定的了解。

（2）、對潮流計算也比以前有了進一步的體會，在學習潮流計算時，雖然依次學習了節(jié)點導納矩陣，功率方程、雅可比矩陣，但不能將它們聯(lián)系起來，更不知道其中的原委，通過程序的編寫，知道了其中的聯(lián)系，也知道了每個方程、矩陣在計算中的作用。

（3）、在接到這個課外課題時，由于以前沒有用過Matlab，剛開始感覺有些茫然，不知道從何做起，但后來通過看書，對Matlab有了一定了解，又覺得和C語言差不多，沒什么難的，真正開始寫程序時，發(fā)現(xiàn)既不是我們想的那么簡單，也不是難得無法下手，說它不是那么簡單，是因為Matlab和C語言有一定區(qū)別，它的數(shù)據(jù)結構主要是矩陣，另外功率方程里的求和部分也不像WORD里那樣，可以用公式編輯器寫出，得用循環(huán)實現(xiàn)。說它不是難得無法下手，是因為通過我們查閱資料和自己調(diào)試，最后完成了潮流計算的程序并得到了收斂的結果。

通過這個任務，自己在Matlab編程，潮流計算，WORD文檔的編輯等方面均有提高，但也暴露出了一些問題：理論知識儲備不足，對Matlab的性能和特點還不能有一個全面的把握，對WORD軟件也不是很熟練，相信通過以后的學習能彌補這些不足，達到一個新的層次。

參考文獻

[1]

黃平，孫德寶．基于低壓電力載波遠程集中抄表系統(tǒng)的研制與實現(xiàn)

[J]．交通與計算機，2004，22(1)：99-101．

Huang Ping，Sun Debao．Study of remote collect meter reading system based on power carrier wave of low voltage[J]．Traffic and Computer，2004，22(1)：99-101(in Chinese)．

[2]

謝開貴，周家啟．樹狀網(wǎng)絡潮流計算的新算法[J]．中國電機工程 學報，2001，21(9)：116-120．

Xie Kaigui，Zhou Jiaqi．A new load flow algorithm for radial distribution networks[J]．Proceeding of the CSEE，2001，21(9)： 116-120(in Chinese)．

[3]

張堯，王琴，宋文南，等．樹狀網(wǎng)的潮流算法[J]．中國電機工程 學報，1998，18(3)：217-220．

Zhang Yao，Wang Qin，SongWennan，et al．A load flow algorithm for radial distribution power networks[J]．Proceeding of the CSEE，1998，18(3)：217-220(in Chinese)．

[4]

車仁飛，李仁俊． 一種少環(huán)配電網(wǎng)三相潮流計算新方法[J]．中國 電機工程學報，2003，23(1)，74-79．

Che Renfei，Li Renjun．A new three-phase power flow method for weakly meshed distribution systems[J]．Proceeding of the CSEE，2003，23(1)，74-79(in Chinese)．

[5]

Cheng Carol S，Shirmohammadi Dariush．A three-phase power flow method for real-time distribution system analysis[J]．IEEE Transaction on Power Systems，1995，10(2)：671-679．

[6]

Zimmerman Ray D，Chiang Hsiaodong．Fast decoupled power flow for unbalanced radial distribution systems[J]．IEEE Transaction on Power Systems，1995，10(4)：2045-2051．

附錄

計算潮流的程序：

%本程序的功能是用牛頓-拉夫遜法進行潮流計算 n=input('請輸入節(jié)點數(shù):n=');nl=input('請輸入支路數(shù):nl=');isb=input('請輸入平衡母線節(jié)電號:isb=');pr=input('請輸入誤差精度:pr=');B1=input('請輸入由支路參數(shù)形成的矩陣:B1=');%變壓器側為1，否則為0 B2=input('請輸入各節(jié)點參數(shù)形成的矩陣:B2=');X=input('請輸入由節(jié)點號及其對地阻抗形成的矩陣:X=');Y=zeros(n);U=zeros(1,n);cta=zeros(1,n);V=zeros(1,n);O=zeros(1,n);S1=zeros(nl);for i=1:n if X(i,2)~=0;p=X(i,1);Y(p,p)=X(i,2);end end for i=1:nl

if B1(i,6)==0 p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;end %求導納矩陣

G=real(Y);B=imag(Y);for i=1:n cta(i)=angle(B2(i,3));U(i)=abs(B2(i,3));%V(i)=B2(i,4);end for i=1:n S(i)=B2(i,1)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);end P=real(S);Q=imag(S);ICT1=0;IT2=1;while IT2~=0 IT2=0;t1=1;t2=1;for i=1:n if i~=isb C(i)=0;D(i)=0;for j1=1:n

C(i)=C(i)+U(i)*U(j1)*(G(i,j1)*cos(cta(i)-cta(j1))+B(i,j1)*sin(cta(i)-cta(j1)));

D(i)=D(i)+U(i)*U(j1)*(G(i,j1)*sin(cta(i)-cta(j1))-B(i,j1)*cos(cta(i)-cta(j1)));end DP(t1)=P(i)-C(i);t1=t1+1;if B2(i,6)==2 DQ(t2)=Q(i)-D(i);t2=t2+1;end end end t1=t1-1;t2=t2-1;

DPQ=[DP';DQ'];%求DP,DQ for i=1:t1+t2 if abs(DPQ(i))>pr IT2=IT2+1;end end H=zeros(t1,t1);N=zeros(t1,t2);K=zeros(t2,t1);L=zeros(t2,t2);for i=1:t1 for j1=1:t1 if j1~=isb&j1~=i

H(i,j1)=0-U(i)*U(j1)*(G(i,j1)*sin(cta(i)-cta(j1))-B(i,j1)*cos(cta(i)-cta(j1)));elseif j1~=isb&j1==i H(i,j1)=U(i)^2*B(i,j1)+D(i);end end end for i=1:t1 for j1=1:t2 if j1~=isb&j1~=i

N(i,j1)=0-U(i)*U(j1)*(G(i,j1)*cos(cta(i)-cta(j1))+B(i,j1)*sin(cta(i)-cta(j1)));elseif j1~=isb&j1==i N(i,j1)=0-U(i)^2*G(i,j1)-C(i);end end end for i=1:t2 for j1=1:t1 if j1~=isb&j1~=i K(i,j1)= U(i)*U(j1)*(G(i,j1)*cos(cta(i)-cta(j1))+B(i,j1)*sin(cta(i)-cta(j1)));elseif j1~=isb&j1==i K(i,j1)=U(i)^2*G(i,j1)-C(i);end end end for i=1:t2 for j1=1:t2 if j1~=isb&j1~=i

L(i,j1)=0-U(i)*U(j1)*(G(i,j1)*sin(cta(i)-cta(j1))-B(i,j1)*cos(cta(i)-

cta(j1)));elseif j1~=isb&j1==i L(i,j1)=U(i)^2*B(i,j1)-D(i);end end end J=[H,N;K,L];%求雅可比矩陣 modify=-JDPQ;Dcta=modify([1:t1],:);t3=U(:,[1:t2]);DU=diag(t3,0)*modify([t1+1:t1+t2],:);t4=1;for i=1:t1 if B2(i,6)~=1 cta(1,i)=cta(1,i)+Dcta(t4,1);t4=t4+1;end end t5=1;for i=1:t2 if B2(i,6)==2 U(1,i)=U(1,i)+DU(t5,1);t5=t5+1;end end ICT1=ICT1+1;end %修正原值 for i=1:n UU(i)=U(i)*cos(cta(i))+1i*U(i)*sin(cta(i));end for p=1:n c(p)=0;for q=1:n c(p)=c(p)+conj(Y(p,q))*conj(UU(q));end s(p)=UU(p)*c(p);end disp('------------------');disp('各節(jié)點電壓U為（節(jié)點從小到大排列）:');disp(UU);disp('------------------');disp('各節(jié)點電壓相角為（節(jié)點從小到大排列）:');

disp(180*angle(UU)/pi);disp('------------------');disp('按公式計算全部線路功率，結果如下:');for i=1:nl if B1(i,6)==0 p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);end

Si(p,q)=UU(p)*(conj(UU(p))*conj(B1(i,4)./2)+(conj(UU(p)*B1(i,5))-conj(UU(q)))*conj(1./(B1(i,3)*B1(i,5))));%各條支路首端功率Si f=[p,q,Si(p,q)];disp(f);end for i=1:nl if B1(i,6)==0 p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);end

Sj(q,p)=UU(q)*(conj(UU(q))*conj(B1(i,4)./2)+(conj(UU(q)./B1(i,5))-conj(UU(p)))*conj(1./(B1(i,3)*B1(i,5))));%各條支路末端功率Sj f=[q,p,Sj(q,p)];disp(f);end disp('------------------');disp('各條支路的功率損耗DS為(順序同您輸入B1時一樣):');for i=1:nl if B1(i,6)==0 p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);end DS(i)=Si(p,q)+Sj(q,p);%各條支路功率損耗DS disp(DS(i));end Sp=0;for i=1:n Sp=Sp+UU(isb)*conj(Y(isb,i))*conj(UU(i));end disp('平衡節(jié)點的功率:');disp(Sp);13

第三篇：潮流等于時尚？范文

潮流等于時尚？

姓名：蔡瑩瑩

班級：服設21102

學號：1114031059 中文摘要：

同人類歷史一樣，服裝的產(chǎn)生于發(fā)展經(jīng)歷了一個漫長的時期，它伴隨著原始人從遠古走來，又緊跟著現(xiàn)代人向未來走去。無論是東方還是西方、古代還是當代，服裝的演變直接的反應了這個時代所代表的時尚。時尚就是在一定“時間”里或長期的“時間”里“崇尚”某些事物;也就是“為時尚早；長時崇尚”，包含這兩種意思。能夠使人倡導或示范的事物?！盀闀r尚早”就是流行時尚。流行時尚就是在特定時段內(nèi)率先由少數(shù)人實驗、而后來為社會大眾所崇尚和仿效的生活樣式。

關鍵詞：潮流 時尚 風格 流行

論文正文

現(xiàn)如今，人們看見一個人的打扮很有個性，也不見得是很時尚的就會說：嘿，你好潮！那么我會經(jīng)常想潮流是否就是等于時尚的呢，答案肯定是否定的，但是它們之間是肯定存在著千絲一縷的關系。

香奈兒女士曾經(jīng)說過：潮流稍縱即逝，只有風格永存。我們就不說這句話是否是對或錯，但是這句話卻是體現(xiàn)了時尚它的精髓，但什么也不是一成不變的。時尚(fashion)不僅僅就是風格(style)，它也要融入源源不斷地新鮮“血液”。雖然潮流是會在一段時期內(nèi)漸漸消失在人們的視線之中，但是潮流（trend）會比我們所謂的風格更快的走入人們的眼中，影響更大，范圍更廣，更有賣點，更有話題。

當然，在我眼中，時尚是必須要有風格的，沒有風格就不成就所謂的時尚，沒有風格就會讓人感覺到非?；靵y。在有了風格之后，加入一些當下的一些潮流元素，才不會顯得呆板，沒有生氣。不僅社會在進步，時尚必然更是要進步，甚至要走在時代最前沿才能永遠保持活躍。

此時，有人會問，潮流究竟是什么個意思？到底什么是潮流？

潮流，它是設計師們靠創(chuàng)造力和自身所迸發(fā)出的有理性的激情。它像空中的風箏，飛得搖擺不定；它像縷縷輕煙，去得形影無綜；它像天上的流星，一現(xiàn)即逝；它喜歡即興表演，因為它沒有過去與未來，你根本就來不及琢磨。如果稍有半點遲疑，它就會“春風忽來花如海了”。它總是新的，總是能閃亮，因為那些流傳下來的經(jīng)典正是源自它；你或許可以說它是源于激情的，因為它們總是不厭倦地做著最前沿的嘗試。

潮流，一方面反映了具有相當數(shù)量人群的意愿和行動需求，同時還體現(xiàn)了一個時代的精神風貌。就服裝的本身來講，幾乎每一個人部在不斷地尋求一種新的式樣、一種新的視覺形象，以滿足人們心理上的需要。但是，另一方面人們又會受到社會的習俗、觀念、生活過程等影響而產(chǎn)生種種制約。這樣，某些服裝式樣、色彩等往往在人們的個性表現(xiàn)和社會規(guī)范(制約)之間起著平衡協(xié)調(diào)的作用，而流行則正是在不同時期、不同條件下對這一特征的充分反映。因此當人們談到潮流的時候，總是會將它與一些潮流概念相仿的俗語聯(lián)系起來，如風尚、時髦、時興、時尚等。這里可以體現(xiàn)出時尚與潮流是有著密不可分的關系，缺一不可。

在課本中，書上是這樣解釋潮流，也可以說是流行： 流行，是指某一事物在某一時期、某一地區(qū)為廣大群眾所接受、所昌愛，并帶有傾向性色彩的杜會現(xiàn)象。服裝的流行定義同樣呂有這些特點，并有非常明顯的時間性和地域性，它必然是隨著時代潮流利社會發(fā)展而產(chǎn)生的，它設計和體現(xiàn)看人們心理上的滿足感、刺激感、新鮮感和愉悅感。流行已經(jīng)成為當代人類生活的一個基本特征，這一點在服裝樣式造型(包括色彩)上面的表現(xiàn)尤為突出。流行服裝是人們穿著的一種新潮流，是時新的時裝樣式，是時裝作品中被人們普遍接受的服裝，是流行服裝的高潮。如歐洲每年的高級時裝發(fā)布會，其中成衣制造高從發(fā)布會上發(fā)現(xiàn)和選擇認為茍合時尚潮流的、能夠引起流行的服裝信息，目口風格特點和造型特征、材質(zhì)的選澤、流行色的運用、配件的搭配等，經(jīng)過提煉、概括和簡化，進行再設計，制作出不同批呈的適合市場需求成衣，作為新的流行款式投放市場，從而構成一定規(guī)模的流行。

那什么又是時尚呢？

各種不同顏色的頭發(fā)可能是我們說的所謂的時尚，但是時尚就是那些么？答案肯定是否定的。漂染彩發(fā)是為了追求與眾不同，所以這里的時尚是“個性化”——時尚是時代帶來的，工業(yè)和科技發(fā)展的時代使人們越來越崇尚個性化。計算機用50年的時間改變了這個世界，網(wǎng)絡又用了不到10年的時間再次改變了這個世界。時尚是對一個時代的傳達，所以時尚不是靠個別天才就能在奇思妙想中誕生的，而是一個時代人們心里吹出的風，它不是形式，而是內(nèi)容。

時尚就是在一定“時間”里或長期的“時間”里“崇尚”某些事物;也就是“為時尚早；長時崇尚”，包含這兩種意思。

而時尚，是走在最前者，是先驅(qū)。而且，在如今的社會中，名人炒作，社會名流等也總會在時尚圈走一遭，讓人們更加深刻的記住他們。不管他們是否是真的時尚，或是趕趕潮流，現(xiàn)實是越來越多的人走進這五彩繽紛的花花世界中。正所謂：世間萬物，皆雙面性。時尚不是時髦。時尚其實是一種鮮時的個性。時髦靚麗是一種時尚，樸素大方也是一種時尚；粉墨登場是一種時尚，素面朝天是一種時尚；肌膚白皙，嬌美可人是時尚，小麥色的皮膚和爽朗的笑聲也是時尚；看一本流行的書是時尚，沉醉于“小橋流水人家”的古韻中也是時尚；物欲橫流的社會忙著賺錢是時尚，不為所動，一杯清茶，一身書香氣也是時尚。時尚其實伴隨著每個人。認真生活，保持個性的人就是時尚的人。

事實上，如果是作為一名服裝設計師，我們首先要搞清楚——在我們眼中，時尚和潮流各代表著什么，我們要如何來看待它們之間的關系，還要平衡他們之間的關系。

在我看來，時尚是各種各樣的風格，我們加入自己對當下的一些帶理性的創(chuàng)意就成就了時尚。但是這樣說來容易，實則不易。難道時尚的就是這么輕易的被創(chuàng)造？——其實它一定是跟隨著社會走的，絕不能與社會脫節(jié)。比如第二次世界大戰(zhàn)結束后，人們沉浸在勝利的喜悅中，對衣著要求有了新的向往。法匡著名的時裝大師迪奧正是順應了這一思潮，設計出“新風貌”（new look）的裝，一舉成名，風靡全球，以至于流行整個20世紀40—60年代。“新風貌”時裝幾乎司以稱為第二次世界大戰(zhàn)后時裝史上的劃時代作品，可見影響之大、之深。

這些有個性的時裝逐漸地大范圍變成人們口中所說的潮流，因為畢竟正真走在時尚前沿的人數(shù)比例是很小的，其他的人大多數(shù)是跟風，這就是所謂的潮流了。

流行與時尚它們也是互相較量的，究竟是時尚引導著潮流，還是潮流引導著時尚？眾說紛紜。

時尚與潮流的相互關系：

（1）時尚的流行范圍比時髦狹小，持續(xù)時間短，所受到的社會支持也不強烈。所以，時尚必須不斷地服從新奇原則，才能引起人們的向往，形成流行。（2）時尚與群眾行為相類似，兩者都具有高度的暗示性、無結構性和反抗性。（3）時尚和流行行為具有抒情性。（4）時尚常常流行于社會階層相類似的人群中。時尚的變遷，親手遵循著一定的規(guī)律：①時尚的產(chǎn)生、流行與社會文明成正比。②時尚是循著極端而變。③社會上參與時尚流行的成員，符合常態(tài)分布。④時尚流行過程中，經(jīng)常改變原形。⑤時尚的變遷受社會文化背景的制約

時尚和潮流都是不斷變化不斷反復，他們是有著密切的聯(lián)系的。但這絕不意味著潮流等于時尚，它們不是化等號的。我認為：時尚是一種對生活積極向上的精神，是一種熱忱的態(tài)度，是經(jīng)典的耐人尋味的風格。潮流是從時尚里以一種跳躍的思想分支出來的。他是時尚的助推器，讓我們與時尚走的更近。

有一些人認為時尚就是名牌兒（brand)，如LV，CHANEL,GUCCI,DIOR等等這些奢侈品牌(luxury brand)，也沒人會把這些奢侈品叫潮牌；而潮流，就是炒作，明星效應，跟風，赤橙黃綠青藍紫鉆，各色閃電半閃四分之三閃，潮牌的各種合作，各種骷髏，各種天價，也沒人會把潮牌說成是時尚品牌。

以上所述，又會讓我們想到這樣一個問題——潮牌就象征著潮流，奢侈品牌象征著時尚，那么是否穿著潮牌就會顯得低俗，穿奢侈品會顯得高貴？其實答案是沒有覺得的。它們兩個是有類似的地方，但也都是一種自我的生活態(tài)度，風格，不受影響的，不是跟風的，潮，更注重生活的隨性，舒適，興趣愛好或許會改變的比較快，不會太久。時尚追逐的更多是logo以外的理念，設計元素，個人風格，相對可能持久些。

還有人認為時尚是高雅文化，屬于精英文化 潮流屬于大眾文化，是經(jīng)大眾流行后的統(tǒng)稱。但不管怎樣，我認為時尚與潮流既矛盾但是又是相互吸引。時尚名牌的保質(zhì)期可以有很長一段的時間都不會過期，但是潮流品牌就是應季的服裝了，這也很大的表現(xiàn)了時尚引領潮流，但是這現(xiàn)象是不可避免的，畢竟社會中是沒有那么多的所謂的有錢人或是名流，或是高干子弟等他們有他們的著裝，而普通人有自己對時尚的理解，可能就是穿著仿制或是類似與那些奢侈品的服裝。

在美國紐約，有分著上城區(qū)和下城區(qū)，在電視劇中會說著上城區(qū)的人們是如何穿著的，下城區(qū)的人們又是如何穿著的。這就是說時尚還會分階級嗎，上城區(qū)就是所謂的時尚，下城區(qū)代表著普通人就是所謂的潮流？但是如今，越來越受人們關注的一樣事物就是——街拍（street snap),街拍就是現(xiàn)在主流的潮流穿著，但那些街拍里面也有一些穿著高雅很有風格的男男女女們，說來說去，不管怎樣，我們看待時尚或是潮流都是我們對生活的一種積極的態(tài)度，可能對于時尚的話年齡歲數(shù)大的占比例會多一些，而潮牌，年輕人的比例會多些，也有一些“例外”的人。但這些相同的不同的都是互相所影響的，沒有相互影響，也就沒有現(xiàn)在的時尚與潮流。

第四篇：潮流計算畢業(yè)論文

科學技術學院

畢業(yè)設計（論文）開題報告

題

目：

電力系統(tǒng)潮流分析計算機輔助設計

學 科 部：

信息學科部

專

業(yè)：

電氣工程及其自動化

班

級：

電氣082班

學

號：

7022808070

姓

名：

黃義軍

指導教師：

劉愛國

填表日期：

2011 年 月 日

一、選題的依據(jù)及意義：

電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種基本電氣計算。它的任務是根據(jù)給定的運行條件和網(wǎng)路結構確定整個系統(tǒng)的運行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡中的功率分布以及功率損耗等。電力系統(tǒng)潮流計算的結果是電力系統(tǒng)穩(wěn)定計算和故障分析的基礎。

潮流計算經(jīng)歷了一個由手工, 利用交、直流計算臺到應用數(shù)字電子計算機的發(fā)展過程?，F(xiàn)在的潮流算法都以計算機的應用為前提。

利用電子計算機進行潮流計算從20世紀50年代中期就已經(jīng)開始。此后，潮流計算曾采用了各種不同的方法，這些方法的發(fā)展主要是圍繞著對潮流計算的一些基本要求進行的。一般要滿足四個基本要求： a)可靠收斂 b)計算速度快 c)使用方便靈活 d)內(nèi)存占用量少

它們也是對潮流算法進行評價的主要依據(jù)。

在電力系統(tǒng)運行方式和規(guī)劃方案的研究中，都需要進行潮流計算以比較運行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟性。同時，為了實時監(jiān)控電力系統(tǒng)的運行狀態(tài)，也需要進行大量而快速的潮流計算。因此，潮流計算是電力系統(tǒng)中應用最廣泛、最基本和最重要的一種電氣運算。在系統(tǒng)規(guī)劃設計和安排系統(tǒng)的運行方式時，采用離線潮流計算；在電力系統(tǒng)運行狀態(tài)的實時監(jiān)控中，則采用在線潮流計算。

二、國內(nèi)外研究現(xiàn)狀及發(fā)展趨勢（含文獻綜述）：

在用數(shù)字計算機求解電力系統(tǒng)潮流問題的開始階段，人們普遍采用以節(jié)點導納矩陣為基礎的高斯-賽德爾迭代法（一下簡稱導納法）。這個方法的原理比較簡單，要求的數(shù)字計算機的內(nèi)存量也比較小，適應當時的電子數(shù)字計算機制作水平和電力系統(tǒng)理論水平，于是電力系統(tǒng)計算人員轉(zhuǎn)向以阻抗矩陣為主的逐次代入法（以下簡稱阻抗法）。

20世紀60年代初，數(shù)字計算機已經(jīng)發(fā)展到第二代，計算機的內(nèi)存和計算速度發(fā)生了很大的飛躍，從而為阻抗法的采用創(chuàng)造了條件。阻抗矩陣是滿矩陣，阻抗法要求計算機儲存表征系統(tǒng)接線和參數(shù)的阻抗矩陣。這就需要較大的內(nèi)存量。而且阻抗法每迭代一次都要求順次取阻抗矩陣中的每一個元素進行計算，因此，每次迭代的計算量很大。

阻抗法改善了電力系統(tǒng)潮流計算問題的收斂性，解決了導納法無法解決的一些系統(tǒng)的潮流計算，在當時獲得了廣泛的應用，曾為我國電力系統(tǒng)設計、運行和研究作出了很大的貢獻。但是，阻抗法的主要缺點就是占用計算機的內(nèi)存很大，每次迭代的計算量很大。當系統(tǒng)不斷擴大時，這些缺點就更加突出。為了克服阻抗法在內(nèi)存和速度方面的缺點，后來發(fā)展了以阻抗矩陣為基礎的分塊阻抗法。這個方法把一個大系統(tǒng)分割為幾個小的地區(qū)系統(tǒng)，在計算機內(nèi)只需存儲各個地區(qū)系統(tǒng)的阻抗矩陣及它們之間的聯(lián)絡線的阻抗，這樣不僅大幅度的節(jié)省了內(nèi)存容量，同時也提高了計算速度。

克服阻抗法缺點的另一途徑是采用牛頓-拉夫遜法（以下簡稱牛頓法）。牛頓法是數(shù)學中求解非線性方程式的典型方法，有較好的收斂性。解決電力系統(tǒng)潮流計算問題是以導納矩陣為基礎的，因此，只要在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓潮流程序的計算效率。自從20世紀60年代中期采用了最佳順序消去法以后，牛頓法在收斂性、內(nèi)存要求、計算速度方面都超過了阻抗法，成為直到目前仍被廣泛采用的方法。

在牛頓法的基礎上，根據(jù)電力系統(tǒng)的特點，抓住主要矛盾，對純數(shù)學的牛頓法進行了改造，得到了P-Q分解法。P-Q分解法在計算速度方面有顯著的提高，迅速得到了推廣。

牛頓法的特點是將非線性方程線性化。20世紀70年代后期，有人提出采用更精確的模型，即將泰勒級數(shù)的高階項也包括進來，希望以此提高算法的性能，這便產(chǎn)生了保留非線性的潮流算法。另外，為了解決病態(tài)潮流計算，出現(xiàn)了將潮流計算表示為一個無約束非線性規(guī)劃問題的模型，即非線性規(guī)劃潮流算法。

近20多年來，潮流算法的研究仍然非常活躍，但是大多數(shù)研究都是圍繞改進牛頓法和P-Q分解法進行的。此外，隨著人工智能理論的發(fā)展，遺傳算法、人工神經(jīng)網(wǎng)絡、模糊算法也逐漸被引入潮流計算。但是，到目前為止這些新的模型和算法還不能取代牛頓法和P-Q分解法的地位。由于電力系統(tǒng)規(guī)模的不斷擴大，對計算速度的要求不斷提高，計算機的并行計算技術也將在潮流計算中得到廣泛的應用，成為重要的研究領域。

三、本課題研究內(nèi)容

1.熟悉電力系統(tǒng)潮流計算的相關理論。

2.在綜合分析各種電力系統(tǒng)特點的基礎上，運用所學專業(yè)知識，提出一種合理高效的潮流計算算法。

3.熟練運用程序設計語言如C語言。

4.通過軟件編程實現(xiàn)所提出的算法，并通過典型系統(tǒng)進行驗證。

四、本課題研究方案

1、確定一種計算方法，如牛頓-拉夫遜法。

2、結合C語言，編寫一套適用的程序完成潮流計算。

3、選取一典型模型進行驗證，試驗程序是否可靠。

五、研究目標、主要特色及工作進度：

研究目標：提出一種合理高效的潮流計算算法，在保證電力系統(tǒng)供電可靠性和電能質(zhì)量的前提下，盡可能提高潮流計算的效率，降低人力資源消耗。從而提高電力系統(tǒng)運行的經(jīng)濟性。進度安排：

第1周: 收集相關參考資料和相關文獻。

第2周: 總結整理資料，熟習課題。

第3周: 提出初步設計方案。

第4周: 熟悉電力系統(tǒng)潮流計算的相關理論及計算機語言。

第5周: 實習

第6周: 寫實習報告

第7周: 確定一種計算方法。

第8周: 提出一種合理的程序設計方法。

第9周： 畫出設計程序整體流程圖。

第10周: 將整體程序模塊化，并定義出每個模塊的功能。

六、參考文獻：

[1] Tankut Yalcinoz, Onur Ko¨ ksoy.A multiobjective optimization

method to environmental economic diaspatch.2007,29（1）:42-50 [2] X.S.Han，H.B.Gooi.Effective economic dispatch model and algorithm.Electrical Power and Energy Systems.2007, 29(1):113-120 [3] 何仰贊，溫增銀.電力系統(tǒng)分析.武漢：華中科技大學出版社，2002 [4] 王錫凡，方萬良，杜正春.現(xiàn)代電力系統(tǒng)分析.北京：科學出版社，2003 [5] 宋文南，李樹鴻，張堯.電力系統(tǒng)潮流計算.天津：天津大學出版社，1990 [6] 王晶，翁國慶，張有冰.電力系統(tǒng)的MATLAB6/SIMULINK仿真與應用.西安：西安電子科技大學出版社，2008.[7] 王祖佑.電力系統(tǒng)穩(wěn)態(tài)運行計算機分析.北京:水利電力出版社，1987.[8] 周全仁，張清益.電網(wǎng)計算與程序設計.長沙:湖南科學技術出版社，1983.[9] 許主平，周少武，鄒軍安。電力系統(tǒng)計算機輔助設計。北京：中國電力出版社，2001。

第五篇：電力系統(tǒng)潮流計算

南 京 理 工 大 學

《電力系統(tǒng)穩(wěn)態(tài)分析》

課程報告

姓名

XX

學 號： 5*** 自動化學院 電氣工程

基于牛頓-拉夫遜法的潮流計算例題編程報學院(系): 專

業(yè): 題

目: 任課教師 碩士導師 告

楊偉 XX

2015年6月10號

基于牛頓-拉夫遜法的潮流計算例題編程報告

摘要：電力系統(tǒng)潮流計算的目的在于：確定電力系統(tǒng)的運行方式、檢查系統(tǒng)中各元件是否過壓或者過載、為電力系統(tǒng)繼電保護的整定提供依據(jù)、為電力系統(tǒng)的穩(wěn)定計算提供初值、為電力系統(tǒng)規(guī)劃和經(jīng)濟運行提供分析的基礎。潮流計算的計算機算法包含高斯—賽德爾迭代法、牛頓-拉夫遜法和P—Q分解法等，其中牛拉法計算原理較簡單、計算過程也不復雜，而且由于人們引入泰勒級數(shù)和非線性代數(shù)方程等在算法里從而進一步提高了算法的收斂性和計算速度。同時基于MATLAB的計算機算法以雙精度類型進行數(shù)據(jù)的存儲和運算, 數(shù)據(jù)精確度高，能進行潮流計算中的各種矩陣運算，使得傳統(tǒng)潮流計算方法更加優(yōu)化。

一 研究內(nèi)容

通過一道例題來認真分析牛頓-拉夫遜法的原理和方法（采用極坐標形式的牛拉法），同時掌握潮流計算計算機算法的相關知識，能看懂并初步使用MATLAB軟件進行編程，培養(yǎng)自己電力系統(tǒng)潮流計算機算法編程能力。

例題如下：用牛頓-拉夫遜法計算下圖所示系統(tǒng)的潮流分布，其中系統(tǒng)中5為平衡節(jié)點，節(jié)點5電壓保持U=1.05為定值，其他四個節(jié)點分別為PQ節(jié)點，給定的注入功率如圖所示。計算精度要求各節(jié)點電壓修正量不大于10-6。

二 牛頓-拉夫遜法潮流計算 1 基本原理

牛頓法是取近似解x(k)之后，在這個基礎上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根。牛頓迭代法其最大優(yōu)點是在方程f(x)=0的單根附近時誤差將呈平方減少，而且該法還可以用來求方程的重根、復根。電力系統(tǒng)潮流計算，一般來說，各個母線所供負荷的功率是已知的，各個節(jié)點的電壓是未知的（平衡節(jié)點外）可以根據(jù)網(wǎng)絡結構形成節(jié)點導納矩陣，然后由節(jié)點導納矩陣列寫功率方程，由于功率方程里功率是已知的，電壓的幅值和相角是未知的，這樣潮流計算的問題就轉(zhuǎn)化為求解非線性方程組的問題了。為了便于用迭代法解方程組，需要將上述功率方程改寫成功率平衡方程，并對功率平衡方程求偏導，得出對應的雅可比矩陣，給未知節(jié)點賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點電壓不平衡量（未知的）構成了誤差方程，解誤差方程，得到節(jié)點電壓不平衡量，節(jié)點電壓加上節(jié)點電壓不平衡量構成節(jié)點電壓新的初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。2 基本步驟和設計流程圖

形成了雅克比矩陣并建立了修正方程式，運用牛頓-拉夫遜法計算潮流的核心問題已經(jīng)解決，已有可能列出基本計算步驟并編制流程圖。由課本總結基本步驟如下：

1)形成節(jié)點導納矩陣Y；

2)設各節(jié)點電壓的初值，如果是直角坐標的話設電壓的實部e和虛部f；如果是極坐標的話則設電壓的幅值U和相角a；

3)將各個節(jié)點電壓的初值代入公式求修正方程中的不平衡量以及修正方程的系數(shù)矩陣的雅克比矩陣；

4)解修正方程式，求各節(jié)點電壓的變化量，即修正量； 5)計算各個節(jié)點電壓的新值，即修正后的值；

6)利用新值從第（3）步開始進入下一次迭代，直至達到精度退出循環(huán)； 7)計算平衡節(jié)點的功率和線路功率，輸出最后計算結果； ① 公式推導

② 流程圖

三

matlab編程代碼

clear;

% 如圖所示1，2，3，4為PQ節(jié)點,5為平衡節(jié)點

y=0;

% 輸入原始數(shù)據(jù)，求節(jié)點導納矩陣

y(1,2)=1/(0.07+0.21j);

y(4,5)=0;y(1,3)=1/(0.06+0.18j);

y(1,4)=1/(0.05+0.10j);

y(1,5)=1/(0.04+0.12j);

y(2,3)=1/(0.05+0.10j);

y(2,5)=1/(0.08+0.24j);

y(3,4)=1/(0.06+0.18j);

for i=1:5

for j=i:5

y(j,i)=y(i,j);

end

end

Y=0;

% 求節(jié)點導納矩陣中互導納

for i=1:5

for j=1:5

if i~=j

Y(i,j)=-y(i,j);

end

end

end

% 求節(jié)點導納矩陣中自導納

for i=1:5

Y(i,i)=sum(y(i,:));

end

Y

% Y為導納矩陣

G=real(Y);

B=imag(Y);% 輸入原始節(jié)點的給定注入功率

S(1)=0.3+0.3j;

S(2)=-0.5-0.15j;

S(3)=-0.6-0.25j;

S(4)=-0.7-0.2j;

S(5)=0;

P=real(S);

Q=imag(S);

% 賦初值，U為節(jié)點電壓的幅值，a為節(jié)點電壓的相位角

U=ones(1,5);

U(5)=1.05;

a=zeros(1,5);

x1=ones(8,1);

x2=ones(8,1);

k=0;

while max(x2)>1e-6

for i=1:4

for j=1:4

H(i,j)=0;

N(i,j)=0;

M(i,j)=0;

L(i,j)=0;

oP(i)=0;

oQ(i)=0;

end

end

% 求有功、無功功率不平衡量

for i=1:4

for j=1:5

oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)));

oQ(i)=oQ(i)-U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)));

end

oP(i)=oP(i)+P(i);

oQ(i)=oQ(i)+Q(i);

end

x2=[oP,oQ]';

% x2為不平衡量列向量

% 求雅克比矩陣

% 當i~=j時，求H,N,M,L

for i=1:4

for j=1:4

if i~=j

H(i,j)=-U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)));

N(i,j)=-U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)));

L(i,j)=H(i,j);

M(i,j)=-N(i,j);

end

end

end

% 當i=j時，求H,N,M,L

for i=1:4

for j=1:5

if 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)));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)));

M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)));

L(i,i)=L(i,i)-U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)))

end

end

N(i,i)=N(i,i)-2*(U(i))^2*G(i,i);

L(i,i)=L(i,i)+2*(U(i))^2*B(i,i);

end

J=[H,N;M,L]

% J為雅克比矩陣

x1=-((inv(J))*x2);

% x1為所求△x的列向量

% 求節(jié)點電壓新值，準備下一次迭代

for i=1:4

oa(i)=x1(i);

oU(i)=x1(i+4)*U(i);

end

for i=1:4

a(i)=a(i)+oa(i);

U(i)=U(i)+oU(i);

end

k=k+1;

end

k,U,a

% 求節(jié)點注入功率

i=5;

for j=1:5

P(i)=U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)))+P(i);

Q(i)=U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)))+Q(i);

end

S(5)=P(5)+Q(5)*sqrt(-1);

S

% 求節(jié)點注入電流

I=Y*U'

四

運行結果

節(jié)點導納矩陣

經(jīng)過五次迭代后的雅克比矩陣

迭代次數(shù)以及節(jié)點電壓的幅值和相角（弧度數(shù)）

節(jié)點注入功率和電流

五 結果分析

在這次學習和實際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機算法中的牛頓-拉夫遜法進行深入的研讀，弄明白了其原理、計算過程、公式推導以及設計流程。牛頓-拉夫遜法是求解非線性方程的迭代過程，其計算公式為?F?J?X，式中J為所求函數(shù)的雅可比矩陣；?X為需要求的修正值；?F為不平衡的列向量。利用x(*)=x(k+1)+?X(k+1)進行多次迭代，通過迭代判據(jù)得到所需要的精度值即準確值x(*)。六 結論

通過這個任務，自己在matlab編程，潮流計算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識儲備不足，對matlab的性能和特點還不能有一個全面的把握，對word軟件也不是很熟練，相信通過以后的學習能彌補這些不足，達到一個新的層次。