【正文】 因此，除了要求計算方法盡可能適應(yīng)各種修改 .調(diào)整以外，還要注意輸入和輸出的方便性和靈活性，加強人機聯(lián)系，以便使計算人員能及時監(jiān)視計算過程并適當(dāng)?shù)乜刂朴嬎愕倪M行。為了得到一個合理的運行方式，往往需要不斷根據(jù)計算結(jié)果，修改原始數(shù)據(jù)。 為了彌補這個缺點，潮流程序的編制必須盡可能使計算人員在計算機計算的過程中加強對計算機過程的監(jiān)視和控制，并便于作各種修改和調(diào)整。這樣，計算的過程就相當(dāng)于運算人員丟系統(tǒng)進行操作 .調(diào)整的過程，非常直觀，物理概念也很清楚。過去在很長時間內(nèi)，電力系統(tǒng)潮流計算是借助于交流臺 進行的。潮流是確定電力網(wǎng)絡(luò)運行狀態(tài)的基本因素，潮流問題是研究電力系統(tǒng)穩(wěn)態(tài)問題的基礎(chǔ)和前提。)。),title(39。迭代次數(shù)’）， ylabel(39。 end disp(39。 disp(DS(i))。q=B1(i,1)。q=B1(i,2)。 end disp(39。 end Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(xonj(E(q)./B1(i,5))conj(E(p)))*xonj(1./(B1(i,3)*B1(i,5))))。 else p=B1(i,2)。各條支路的末端功率 Sj為（順序同您的輸入 B1時一樣）：‘）； 30 for i=1:n1 if B1(i,6)==0 p=B1(i,1)。 disp(Si())。q=B1(i,1)。q=B1(i,2)。 disp(39。 end disp(39。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。 disp(O)。各節(jié)點的電壓相角 O為（節(jié)點號從小到大的排列）： 39。 disp(V)。各節(jié)點的電壓大小 V為（節(jié)點號從小到大的排列）： 39。 disp(E)。各節(jié)點的實際電壓標(biāo)么值 E為（節(jié)點號從小到大的排列）： 39。 end E=e+f*j。 for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。)。 disp(39。)。 end end %用高斯消去法解“ w=J*V” disp(39。 end for i=1:n Dy(k)=sqrt(e(k)^2+f(k)^2)。 if DET=pr IT2=IT2+1 end end ICT2(a)=IT2 ICT1=ICT1+1。 f(L)=f(L)J(k1,N)。 e(L)=e(L)J(k,N)。 end J(k3,k)=0。 end J(k3,k)=0。 if k~=3 k4=k1。 for k2=k1:N1 J(k,k2)=J(k,k2)./J(k,k)。 end end end end end %求雅可比矩陣 for k=3:N0 k1=k+1。J(p,q)=X6。J(m,N)=DP。m=p+1。J(p,q)=X5。 p=2*i1。 X5=2*e(i)。j1~=isb X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。J(m,q)=X2。q=q+1。 J(m,q)=X1。J(p,N)=DV。q=2*j11。 X6=0。 X2=B(i,j1)*e(i)G(i,j1)*f(i)。 for j1=1:n if j1~=isbamp。 end end else DP=P(i)P1。J(p,q)=X4。J(m,N)=DP。m=p+1。J(p,q)=X3。 p=2*i1。 X3=X2。j1~=i X1=G(i,j1)*e(i)B(i,j1)*f(i)。 DQ=Q(i)Q1。 V2=e(i)^2+f(i)^2。 %求 39。 end P1=C(i)*e(i)+f(i)*D(i)。 for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。 if i~=isb C(i)=0。a=a+1。a=0。NO=2*n。 27 ICT1=0。 end P=rea(S)。 end for i=1:n S(i)=B2(i,1)B2(i,2)。 f(i)=imag(B2(i,3))。B=imag(Y)。 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。 Y(p,q)=Y(p,q)。q=B1(i,1)。q=B1(i,2)。 Y(p,p)=1./X(i,2)。 for i=1:n if X(i,2)~=0。O=zeros(1,n)。f=zeros(1,n)。 Y=zeros(n)。請輸入由節(jié)點參數(shù)形成的矩陣： X=39。)。 B2=input(39。請輸入由支路參數(shù)形成的矩陣： B1=39。)。 pr=input(39。請輸入平衡母線節(jié)點號： isb=39。)。 n1=input(39。請輸入節(jié)點數(shù)： n=39。 邏輯關(guān)系運算 邏輯運算是 MATLAB 中數(shù)組運算所特有的一種運算形式，也是幾乎所有的高級語言普遍適用的一種運算。有了“對應(yīng)元素”的規(guī)定，數(shù)組的運算實質(zhì)上就是針對數(shù)組內(nèi)部的每個元素進行的。在數(shù)組運算中有了“對應(yīng)關(guān)系”的規(guī)定，數(shù)組與常數(shù)之間的除法運算沒有任何限制。而乘除法運算有相當(dāng)大的區(qū)別，數(shù)組的乘除法是指兩同維數(shù)組對應(yīng)元素之間的乘除法，它們的運算符為“ .*”和“ ./”或“ .\”。 這種運算不同于前面講的數(shù)學(xué)運算，為有所區(qū)別，我們稱之為數(shù)組運算。但需注意進行數(shù)除時，常數(shù)通常只能做除數(shù)。在 MATLAB6中兩者的區(qū)別不太大。在傳統(tǒng)的 MATLAB 算法中，右除是先計算矩陣的逆再相乘，而左除則不需要計算逆矩陣直接進行除運算。 矩陣的基本數(shù)學(xué)運算包括矩陣的四則運算、與常數(shù)的運算、逆運算、行列式運算、秩運算、特征值運算等基本函數(shù)運算，這里進行簡單介紹。 原始數(shù)據(jù)輸入格式的設(shè)計，主要應(yīng)從使用的角度出發(fā)，原則是簡單明了，便于修改。不同領(lǐng)域、不同層次的用戶通過相應(yīng)工具的學(xué)習(xí)和應(yīng)用，可以方便地進行計算、分析及設(shè)計工作。通過 M語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時間，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。它可以高效率地解決工業(yè)計算問題，特別是關(guān)于矩陣和矢量的計算。 MATLAB 是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運算，同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強大的功能。現(xiàn)有很多潮流計算方法。( 1, 2, , 。 ⒉導(dǎo)納矩陣中的某些非對角元素為零時 ,雅可比矩陣中對應(yīng)的元素也是為零 .若 22 0ijY? ,則必有 0ijJ ? 。 19 以直角坐標(biāo)系形式表示 ① . 迭代推算式 采用直角坐標(biāo)時 ,節(jié)點電壓相量及復(fù)數(shù)導(dǎo)納可表示為 : i i iij ij ijV e jfY G jB???? (323) 將以上二關(guān)系式代入上式中 ,展開并分開實部和虛部 。 事實上，除了平衡節(jié)點的功率方程式在迭代過程中沒有約束作用以外，其余每個節(jié)點都可以列出兩個方程式。 18 牛頓 拉夫遜法潮流求解過程 以下討論的是用直角坐標(biāo)形式的牛頓 — 拉夫遜法潮流的求解過程。解決這個問題的辦法可以用高斯法迭代 1~2 次，以此迭代結(jié)果作為牛頓法的初值。 )i q n i s?? (316) 這樣一般能得到滿意的結(jié)果。如果初值選擇不當(dāng)，算法有可能根本不收斂或收斂到一個無法運行的節(jié)點上。牛頓法所需的內(nèi)存量及每次迭代所需時間均較高斯法多。而且其迭代次數(shù)與所計算網(wǎng)絡(luò)的規(guī)?；緹o關(guān)。牛頓法當(dāng)初始估計值 (0)x 和方程的精確解足夠接近時，收斂速度非常快，具有平方收斂特性。k 為迭代次數(shù)。 ( ) ( ) ( )( ( )k k kf x x f x? ? ? (314) ( 1 ) ( ) ( )k k kx x x? ? ? ? (315) 上兩式中： 39。接著就從 (1)x 出發(fā)，重復(fù)上述計算過程。由此可以求得第一次迭代的修正量 ( 0 ) 39。 對于非線性代數(shù)方程組： ( ) 0fx? 即 12( , , , ) 0inf x x x ? ( 1,2, , )in? (311) 在待求量 x 的某一個初始估計值 (0)x 附近，將上式展開成泰勒級數(shù)并略去二階及以上的高階項，得到如下的經(jīng)線性化的方程組： ( 0 ) 39。其要點是把非線性方程式的求解過程變成反復(fù)地對相應(yīng)的線性方程式進行求解的過程。實際電力系統(tǒng)的潮流技術(shù)那主要采用牛頓 拉 夫遜法。對現(xiàn)有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設(shè)計以及對電力系統(tǒng)進行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計算為基礎(chǔ)。即節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷。 第三章 牛頓－拉夫遜法概述 牛頓 拉夫遜法基本原理 電力系統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。常用的方法是迭代法和牛頓法，在計算過程中，或得出結(jié)果之后用約束條件進行檢驗。這一約束的主要意義就在于此。 3. 節(jié)點的有功功率和無功功率應(yīng)滿足 4. m in m a xm in m a xG i G i G iG i G i G iP P PQ Q Q?? ???? ? (219) PQ節(jié)點的有功功率和無功功率，以及 PU 節(jié)點的有功功率，在給定是就必須滿足上述條件，因此，對平衡節(jié)點的 P 和 Q 以及 PU 節(jié)點的 Q 應(yīng)按上述條件進行檢驗。 PU 節(jié)點電壓幅值必須按上述條件給定。 電力系統(tǒng)運行必須滿足一定技術(shù)和經(jīng)濟上的要求。 U。 以上三類節(jié)點 4 個運行參數(shù) P。 關(guān)于平衡節(jié)點的選擇，一般選擇系統(tǒng)中擔(dān)任調(diào)頻調(diào)壓的某一發(fā)電廠 (或發(fā)電機 )，有時也可能按其他原則選擇，例如，為提高計算的收斂性。也就是說，對平衡節(jié)點給定的運行參數(shù)是 U 和 ? ，因此有城為 U? 節(jié)點，而待求量是該節(jié)點的 P。 PU 節(jié)點上的發(fā)電機稱為 PU 機(或 PU 給定型發(fā)電機 ) ③ 平衡節(jié)點 在潮流計算中，這類節(jié)點一般只 設(shè)一個。用以維持給定的電壓值。 ② PU 節(jié)點 這類節(jié)點給出的參數(shù)是該節(jié)點的有功功率 P 及電壓幅值 U，待求量為該節(jié)點的無功功率 Q 及電壓向量的相角 ? 。 PQ 節(jié)點上的發(fā)電機稱之為 PQ 機 (或 PQ 給定型發(fā)電機 )。通常變電所母線都是 PQ 節(jié)點，當(dāng)某些發(fā)電機的輸出功率 P。主要目的是由這些已知量去求電力系統(tǒng)內(nèi)的各種電氣量。 用一般的電路理論求解網(wǎng)絡(luò)方程，目的是給出電壓源 (或電流源 )研究網(wǎng)絡(luò)內(nèi)的電流 (或電壓 )分布，作為基礎(chǔ)的方程式，一般用線性代數(shù)方程式表示。例如：利用