【正文】 速恒頻雙饋式發(fā)電機的風(fēng)電場接入系統(tǒng)后對電壓質(zhì)量的影響,在雙饋式發(fā)電機簡化等值電路的基礎(chǔ)上建立了風(fēng)電場的確定性潮流模型,建立了風(fēng)力發(fā)電機的隨機分析模型,并在這二者的基礎(chǔ)上運用基于半不變量法的隨機潮流進行計算。結(jié) 論本設(shè)計要求運用MATLAB進行某電網(wǎng)的潮流計算，經(jīng)過幾個月的努力，初步完成了設(shè)計要求，現(xiàn)總結(jié)如下：（1）經(jīng)過查閱相關(guān)文獻和書籍，基本掌握了潮流計算的手工算法。經(jīng)過對幾種教材仔細的閱讀和揣摩，對幾種手工算法的特點和原理有了比較清晰的認識和了解，能夠運用手工算法對簡單的電力系統(tǒng)網(wǎng)絡(luò)進行分析和計算。（2）對目前流行的潮流計算計算機算法進行了總結(jié)和歸納。 通過對幾種計算機算法的比較，決定采用牛頓拉夫遜法（直角坐標）進行潮流計算。牛頓迭代法是求方程根的重要方法之一，是常用的解非線性方程組的方法，也是當前廣泛采用的計算潮流的方法。其最大優(yōu)點是在方程的單根附近具有平方收斂，收斂性較好，而且該法還可以用來求方程的重根、復(fù)根。但是對初始值的要求比較嚴格。（3）用MATLAB編程并仿真MATLAB語言允許用戶以數(shù)學(xué)形式的語言編寫程序, 其比BASIC 語言和FORTRAN 等更為接近書寫的數(shù)學(xué)表達格式, 且程序易于調(diào)試。在計算要求相同的情況下, 使用MATLAB 編程, 工作量將會大為減少。其程序的編寫也因MATLAB提供了許多功能函數(shù)而變得簡單易行。設(shè)計已經(jīng)做完，能夠順利的運行，但是還有很多不完善的地方，MATLAB的用戶界面由于時間有限，比較復(fù)雜，所以沒能做出來，十分遺憾，還有就是所編寫的程序通用性不是很強，只能針對特定的網(wǎng)絡(luò)，其適用性還有待加強。在畢業(yè)設(shè)計中遇到了很多困難和挫折，但是通過我自己的努力，加上老師和同學(xué)的幫助，是我順利渡過難關(guān)。通過克服這些困難，我個人的能力有了很大提升，我的知識也得到了運用和升華，我想這也就是本次畢業(yè)設(shè)計最核心的目的，使我進入社會后能更快地適應(yīng)工作崗位，出色地完成工作任務(wù)。這一階段的設(shè)計使我的大學(xué)生活過的無比充實，也為我的大學(xué)生涯畫上了一個圓滿的句號，在今后的工作和學(xué)習(xí)中，我應(yīng)該繼續(xù)嚴格要求自己，充滿自信的對待每一件事，充分發(fā)揮自己的價值，成為一個對社會有用的人。致 謝在大學(xué)的學(xué)習(xí)和生活即將結(jié)束，我衷心的感謝所有曾經(jīng)指導(dǎo)過我的老師和幫助我的同學(xué)們。他們在我的學(xué)習(xí)生涯中給了我很大的幫助，這次的設(shè)計能夠成功的完成，與他們的幫助是分不開的，在此，特別感謝我的指導(dǎo)老師對我的關(guān)懷和教導(dǎo)。本次的設(shè)計和論文是在**老師的直接指導(dǎo)下完成的，在論文的選題及設(shè)計思路上得到了老師多次的指導(dǎo)與幫助，并且多次給予了我寶貴的意見，不厭其煩的幫我修改論文，接了我無數(shù)次的電話，但老師始終耐心的為我解答各種在設(shè)計中遇到的問題。，在本次畢業(yè)設(shè)計過程中，能夠得到老師的指導(dǎo)使我感到十分的榮幸，老師嚴謹?shù)闹螌W(xué)作風(fēng)、熱情的待人態(tài)度給我留下了深刻的印象。在論文完成之際，再次表示對老師衷心的感謝，同時也對在設(shè)計過程中幫助過我的其他老師和同學(xué)，在此也表示誠摯的謝意。對給予我?guī)椭^的老師和同學(xué)們再次表示萬分的感謝！參 考 文 獻1 [M].北京：中國電力出版社，2007. 2 李維波. MATLAB 在電氣工程中的應(yīng)用［M］. 北京: 中國電力出版社，2007.3 武曉朦，張飛廷. 電力系統(tǒng)的P Q 分解法潮流計算[J].現(xiàn)代電子技術(shù),2002 ,142(11): 105106 4 謝威，彭志煒，張朝綱，—拉夫遜的潮流就算方法[J].許昌學(xué)院學(xué)報. ,25(2):27305 [J].山西冶金，2007，106（2）：42446 劉軍,劉學(xué)軍. MATLAB 在電力系統(tǒng)分析中的應(yīng)用[J]. 電力系統(tǒng)及其自動化學(xué)報, ,12(2): 23257 李家坤,劉姣姣. MATLAB仿真技術(shù)在電力電子教學(xué)中的應(yīng)用[J]. ,26(1): 75778 連小洲，焦?jié)? MATLAB在電力系統(tǒng)計算中的應(yīng)用[J]. 江西電力職業(yè)技術(shù)學(xué)院學(xué)報，,18(4): 10119 張寧,江紅梅,張渭. 基于MATLAB的電力系統(tǒng)潮流計算[J].西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版)，,32(12): 12412610 周衛(wèi)星，[J].科技咨詢導(dǎo)報，2007，10：707111 徐一哲，[J].中外企業(yè)家，2009,5(下)：20620812 趙芳,陳士方,[J].電力情報，2000,03：222513 劉浩,戴居豐. 基于系統(tǒng)分割的保留非線性快速 PQ 分解狀態(tài)估計[J].2005,06：737614 [J].[碩士學(xué)位論文]，長沙理工大學(xué)，200515 [J].[碩士學(xué)位論文],2006.16 羅緯. 電力市場環(huán)境下基于最優(yōu)潮流的輸電容量充裕度研究[J].[碩士學(xué)位論文],上海交通大學(xué),200817 劉方. 電力系統(tǒng)動態(tài)最優(yōu)潮流的模型與算法研究[J].[博士學(xué)位論文],重慶大學(xué),2007.18 曾超，彭建春，金靈滿，周元清. 基于改進布羅伊登法的交直流潮流計算[J].繼電器，2008,06(01)6819 閏麗梅，陳娟，徐建軍. 基于直流潮流和分布因子三母線系統(tǒng)脆性源辨識技術(shù)[J].電氣技術(shù)，2007，06：555820 劉蘇琴，[J],江蘇電機工程，2007，04(26)：1316附錄A 基于MATLAB的牛頓拉夫遜法（直角坐標）潮流計算程序清單clear。clcy=0。%輸入原始數(shù)據(jù)，求節(jié)點導(dǎo)納矩陣y (1,2)=1/(+)。 y (1,3)=1/(+)。 y (1,4)=1/(+)。y(1,5)=1/(+)。y(2,3)=1/(+)。y(2,5)=1/(+)。y(3,4)=1/(+)。y(4,5)=0。for i=1:5for j=i:5y(j,i)=y(i,j)。endendY=0。%求互導(dǎo)納for i=1:5for j=1:5if i~=jY(i,j)=y(i,j)。endendend%求自導(dǎo)納for i=1:5Y(i,i)=sum(y(i,:))。endY %Y 為導(dǎo)納矩陣G=real(Y)。B=imag(Y)。%原始節(jié)點功率S(1)=+。S(2)=。S(3)=。S(4)=。S(5)=0。P=real(S)。Q=imag(S)。%賦初值U=ones(1,5)。U(5)=。e=zeros(1,5)。ox=ones(8,1)。fx=ones(8,1)。count=0 %計算迭代次數(shù)while max(fx)1e5for 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。 endendfor i=1:4for j=1:5 oP(i)=oP(i)U(i)*U(j)*(G(i,j)*cos(e(i)e(j))+B(i,j)*sin(e(i)e(j)))。oQ(i)=oQ(i)U(i)*U(j)*(G(i,j)*sin(e(i)e(j))B(i,j)*cos(e(i)e(j)))。endoP(i)=oP(i)+P(i)。 oQ(i)=oQ(i)+Q(i)。endfx=[oP,oQ]39。%求雅克比矩陣%當i~=j時候求H,N,M,L 如下：for i=1:4for j=1:4if i~=j H(i,j)=U(i)*U(j)*(G(i,j)*sin(e(i)e(j))B(i,j)*cos(e(i)e(j)))。N(i,j)=U(i)*U(j)*(G(i,j)*cos(e(i)e(j))+B(i,j)*sin(e(i)e(j)))。L(i,j)=H(i,j)。M(i,j)=N(i,j)。endendendH,N,M,L%當i=j 時H,N,M,L如下：for i=1:4for j=1:5if i~=jH(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(e(i)e(j))B(i, j)*cos (e(i)e(j)))。 N(i,i)=N(i,i)U(i)*U(j)*(G(i, j)*cos(e(i)e(j))+B(i,j)*sin(e(i)e(j)))。M(i,i)=M(i,i)U(i)*U(j)*(G(i,j)*cos(e(i)e(j))+B(i,j)*sin(e(i)e(j)))。 L(i,i)=L(i,i)U(i)*U(j)*(G(i,j)*sin(e(i)e(j))B(i,j)*cos(e(i)e(j))。endendN(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)。endJ=[H,N。M,L] %J 為雅克比矩陣ox=((inv(J))*fx)。for i=1:4oe(i)=ox(i)。 oU(i)=ox(i+4)*U(i)。endfor i=1:4e(i)=e(i)+oe(i)。 U(i)=U(i)+oU(i)。endcount=count+1。endox,U,e,count%求節(jié)點注入的凈功率i=5。for j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)e(j))+B(i,j)*sin(e(i)e(j)))+P(i)。Q(i)=U(i)*U(j)*(G(i,j)*sin(e(i)e(j))B(i,j)*cos(e(i)e(j)))+Q(i)。endS(5)=P(5)+Q(5)*sqrt(1)。S%求節(jié)點注入電流I=Y*U39。