【正文】 修正、迭代 給定 PQ 節(jié)點(diǎn)初值和各 節(jié)點(diǎn)電壓相角初值 V1= ， V2(0)=V3(0)=， V4= δ 2(0)=δ 3(0)=0， δ 4(0)=0 1 作第一次有功迭代，按公式計(jì)算節(jié)點(diǎn)有功功率不平衡量 迭代中的△ P： 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 12 2 做第一次無功迭代，按公式計(jì)算無功功率不平衡量，計(jì)算時(shí)電壓相角最新的修正值。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語言（如 C、Fortran）的編輯模式，代表了當(dāng)今國際科學(xué)計(jì)算軟件的先進(jìn)水平 MATLAB 和 Mathematica、 Maple 并稱為三大數(shù)學(xué)軟件。在新的版本中也加入了對 C， FORTRAN， C++ ， JAVA 的支持。 常數(shù)與矩陣的運(yùn)算即是同該矩陣的每一元素進(jìn)行計(jì)算。 數(shù)組的加、減與矩陣的加、減運(yùn)算完全相同。 另外，矩陣的數(shù)組運(yùn)算中還有冪運(yùn)算（運(yùn)算符 .^）、指數(shù)運(yùn)算（ exp）、對數(shù)運(yùn)算 (log)、和開方運(yùn)算（ sqrt）、等，有了“對應(yīng)元素”的規(guī)定，數(shù)組的運(yùn)算實(shí)質(zhì)上就是針對數(shù)組內(nèi)部的每個(gè)元素進(jìn)行的。 2） 可用 ↑ 、 ↓ 鍵來重現(xiàn)已輸入的數(shù)據(jù)或命令。 4） 常用數(shù)有特定的名字，如 pi（ =） 、 Inf（ =∞ ）、 NaN 則表示不定型求得的結(jié)果（如 0/0）。7 8 9] a = 1 2 3 4 5 6 7 8 9 7） 需要顯示命令的計(jì)算結(jié)果時(shí)，則語句后面不加 “ ； ” 號，否則要加 “ ； ” 號。4 5 6。 但這個(gè)命令不是存儲數(shù)據(jù)，而是存儲輸入與屏幕上輸出的內(nèi)容。 10） save 命令 該命令存儲定義的變量或演算結(jié)果，也可以用來存儲指定的變量。 1) 在 MATLAB 中表示一個(gè)矢量要用方括號， 而列矢量的輸入只需在行矢量輸入格式基礎(chǔ)上加轉(zhuǎn)置符 (‘) 即可。 (加轉(zhuǎn)置符 ) x = 1 4 2 5 3 6 注意上面兩式的區(qū)別。 第三條命令產(chǎn)生 X至 y步長為 n 的行矢量。 3) 矩陣的加 , 減 , 乘 , 除等 , 和其它 語言書寫一樣。 命令格式： [ n， d ] = paralltl（ n1， d1， n2， d2） [a， b， c， d] = paralltl（ a1， b1， c1， d1， a2， b2， c2， d2） 2） 系統(tǒng)的串聯(lián) series 命令實(shí)現(xiàn)兩個(gè)系統(tǒng)的串聯(lián)，示意圖如下： u1 y u2 y2 串聯(lián)后系統(tǒng)的傳遞函數(shù)為 121212( ) ( ) ( ) nng s g s g s dd?? ? ? ? 命令格式： [ n， d ] = series（ n1， d1， n2， d2） [ a， b， c， d ] = series（ a1， b1， c1， d1， a2， b2， c2， d2） 系統(tǒng) 1 系統(tǒng) 2 g1（ s） g2（ s） 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 18 3） 系統(tǒng)的反饋 feedback 命令實(shí)現(xiàn)兩個(gè)系統(tǒng)的反饋連接，示意圖如下： u1 + y1 177。 第六章 潮流計(jì)算流程圖及源程序 本次課程設(shè)計(jì)采用極坐標(biāo)下的牛頓 拉夫遜計(jì) 算網(wǎng)絡(luò)的潮流計(jì)算。)。 n1=input(39。請輸入 PQ 結(jié)點(diǎn)數(shù)： n2=39。)。 K=input(39。請輸入支路阻抗矩陣： C=39。)。 S=input(39。N=zeros(n1+n2,n2)。 if K(m,R)~=0 C(m,m)=C(m,m)+1/(C(m,R) /( K(m,R) * (K(m,R)1))) 。 end end end for m=1:n for R=1:n if m~=R Z(m)=Z(m)+1/C(m,R)。結(jié)點(diǎn)導(dǎo)納矩陣 :39。迭代中的雅克比矩陣 :39。 O=angle(U)。 P=real(S)。 end PT1(m)=sum(PT)。 end QT1(m)=sum(QT)。 PR2=max(abs())。 end end end for m=1:n1+n2 for R=1:n2 if m==R N(m,m)=U1(m)^2*G(m,m)PT1(m)。 end end end for m=1:n2 for R=1:n2 if m==R L(m,m)=U1(m)^2*B(m,m)QT1(m)。 disp(JJ)。 DA=inv(JJ)*PQ。 end UD2=diag(UD)。 end for m=1:n1+n2 o(k+1,m)=180/pi*O(m)。 c(m)=U1(m)*sin(O(m))。 for m=1:n for R=1:n if m~=R C1(m,R)=1/C(m,R)。 end end end disp(39。 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 27 disp(39。 disp(39。 disp(39。 disp(39。 disp(39。 程序結(jié)果 請輸入結(jié)點(diǎn)數(shù)： n=4 請輸入 PV結(jié)點(diǎn)數(shù)： n1=1 請輸入 PQ結(jié)點(diǎn)數(shù)： n2=2 請輸入平衡結(jié)點(diǎn)： isb=4 請輸入精確度： pr= 請輸入變比矩陣看 :K=[0 0 0 0。 + 0 + +。 0 0 0。同時(shí)也看到了研究性學(xué)習(xí)的效果，從研究中去學(xué)習(xí)，理論結(jié)合實(shí)際，將理論運(yùn)用到實(shí)際，同時(shí)在實(shí)踐中發(fā)現(xiàn)問題，然后解決問題， 通過這次課程設(shè)計(jì)我在 MATLAB 編程，潮流計(jì)算， WORD 文檔的編輯方面均有所提高，但也暴露出了一些問題；我發(fā)現(xiàn)了自己的基礎(chǔ)知識有很多的不足，對 MATLAB的性能和特點(diǎn)還不能有一個(gè)全面的把握，這些基礎(chǔ)的缺乏給我的設(shè)計(jì)計(jì)劃造成了不小的障礙，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，從 而達(dá)到一個(gè)新的層次。在進(jìn)行課題設(shè)計(jì)的過程中，加深了我對潮流計(jì)算的認(rèn)識，尤其是對牛頓拉夫遜潮流計(jì)算的求解思路有 了比較透徹的理解。inf + + 0] 請輸入支路導(dǎo)納矩陣： y=[0 0。0 0 0 0。)。)。)。)。)。)。 end end end for m=1:n for R=1:n if (C(m,R)~=inf)amp。 for R=1:n PH1(R)=U(isb)*conj(Y(isb,R))*conj(U(R))。 end k=k+1。 for m=1:n1+n2 O(m)=O(m)+OO(m)。 for m=1:n1+n2 OO(m)=DA1(m)。39。 end end end JJ=[H N。 end end end for m=1:n2 for R=1:n1+n2 if m==R J(m,m)=U1(m)^2*G(m,m)PT1(m)。 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 23 for m=1:n1+n2 for R=1:n1+n2 if m==R H(m,m)=U1(m)^2*B(m,m)+QT1(m)。 1(k+1,m)=(m)。 PP1(k+1,m)=PP(m)。 while PRpr 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 22 for m=1:n2 UD(m)=U1(m)。 k=0。 G=real(Y)。 disp(Y)。 else Y(m,R)=1/C(m,R)。 C(m,R)=C(m,R)/K(m,R)。QT1=zer