【正文】 （ 220） i i iiiiS P jQIUU??????? （ 221） 1 ( 1, 2 , .. ., )ni ij jjI Y U i n????? （ 222） 由式（ 221）、（ 222）可得 1( 1 , 2 , ..., )nii ij jjiP jQ Y U i nU?? ?? ??? （ 223） 將式（ 219）帶入式（ 223）的右端，展開并分出實(shí)部和虛部，便得 ???? ????nj jijjijinj jijjijii eBfGffBeGeP 11 )()( ? ?? ? ????njnj jijjijijijjijii eBfGefBeGfQ 1 1 )()( （ 224） 按照分類， PQ節(jié)點(diǎn)的有功功率和無功功率是給定的，第 i個(gè)節(jié)點(diǎn)的給定功率設(shè)為isP 和 isQ 。同時(shí)可以看到，方程式（ 225）（ 226）具備方程組（ 28）的形式： UJW ???? （ 227） 式中 ????????????????????????????????????????????1n21n1m21mmm11..UPUPQPQPW ??????????????????????????????????????????????111111..nnmmmmfefefefeU 上式方程中雅可比矩陣 J的各元素，可以對式（ 225）（ 226）求偏導(dǎo)獲得 畢業(yè)設(shè)計(jì)（論文） 15 1 1 1 1 1 1 1 11 1 1 1 1 11 1 1 1 1 1 1 11 1 1 1 1 11 1 1 1m m m m n nm m m m n nm m m m m m mm m m m nP P P P P P P Pe f e f e f e fQ Q Q Q Q Q Q Qe f e f e f e fP P P P P P Pe f e f e f eJ? ? ? ?? ? ? ?? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ??111 1 1 1 1 11 1 1 1 1 1 1 11 1 1 1 1 12 2 2 2 21 1 1 1 111mnm m m m m m m mm m m m n nm m m m m m m mm m m m n nm m m m mmmPfQ Q Q Q Q Q Q Qe f e f e f e fP P P P P P P Pe f e f e f e fU U U U Ue f e f e?? ? ? ?? ? ? ? ? ? ? ?? ? ? ?? ? ? ? ????? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ? ?? ? ? ? ?2 2 21 1 11 1 1 11 1 1 1 1 1 1 11 1 1 1 1 12 2 2 2 2 21 1 1 1 1 11 1 1 1m m mm m n nn n n n n n n nm m m m n nn n n n n nm m m mU U Uf e fP P P P P P P Pe f e f e f e fU U U U U Ue f e f e f? ? ?? ? ? ?? ? ? ? ? ? ? ?? ? ? ?? ? ? ? ? ???? ? ? ? ? ?? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ?221111nnnnUUef?????????????????????????? 當(dāng) ji? 時(shí) , 雅可比矩陣中非 對角元素為 22()0iiij i ij ijjiiij i ij ijjjjPQG e B fefPQB e G ffeUUef??? ??? ? ? ? ? ??? ???? ???? ? ? ??? ???? ??????? ?? （ 228） 當(dāng) ji? 時(shí) ,雅可比矩陣中對角元素為 : 111122()()()()22niij j ij j ii i ii ijiniij j ij j ii i ii ijjniij j ij j ii i ii ijiniij j ij j ii i ii ijjiijiiiPG e B f G e B fePG f B e G f B efQG f B e G f B eeQG e B f G e B ffUeeUff?????? ?? ? ? ? ???????? ? ? ? ? ???????? ? ? ?? ??????? ? ? ? ? ?????? ??????????? ???????? （ 229） 畢業(yè)設(shè)計(jì)（論文） 16 根據(jù)上述原理，可以把潮流計(jì)算的 求解過程大致可以分為以下步驟： (1) 形成節(jié)點(diǎn)導(dǎo)納矩陣； (2)將各節(jié)點(diǎn)電壓設(shè)初值 U (3)將節(jié)點(diǎn) 初值代入相關(guān)求式，求出修正方程式的常數(shù)項(xiàng)向量； (4)將節(jié)點(diǎn)電壓初值代入求式，求出雅可比矩陣元素； (5)求解修正方程，求修正向量； (6)求取節(jié)點(diǎn)電壓的新值； (7)檢查是否收斂，如不收斂，則以各節(jié)點(diǎn)電壓的新值作為初值自第 3 步重新開始進(jìn)行狹義次迭代，否則轉(zhuǎn)入下一步； (8)計(jì)算支路功率分布， PV 節(jié)點(diǎn)無功功率和平衡節(jié)點(diǎn)注入功率。因此，這一約束條件對 PQ節(jié)點(diǎn)而言。 因此，潮流計(jì)算可以歸結(jié)為求解一組非線性方程組，并使其解答滿足一定的約束條件?，F(xiàn)有很多潮流計(jì)算方法。它可以高效率地解決工業(yè)計(jì)算問題，特別是關(guān)于矩陣和矢量的計(jì)算。不同領(lǐng)域、不同層次的用戶通過相應(yīng)工具的學(xué)習(xí)和應(yīng)用，可以方便地進(jìn)行計(jì)算、分析及設(shè)計(jì)工作。矩陣的基本數(shù)學(xué)運(yùn)算包括矩陣的四則運(yùn)算、與常數(shù)的運(yùn)算、逆運(yùn)算、行列式運(yùn)算、秩運(yùn)算、特征值運(yùn)算等基本函數(shù)運(yùn)算，這里進(jìn)行簡單介紹。在 MATLAB7 中兩者的區(qū)別不太大。這種運(yùn)算不同于前面講的數(shù)學(xué)運(yùn)算，為有所區(qū)別，我們稱之為數(shù)組運(yùn)算。在數(shù)組運(yùn)算中 有了“對應(yīng)關(guān)系”的規(guī)定，數(shù)組與常數(shù)之間的除法運(yùn)算沒有任何限制。 邏輯關(guān)系運(yùn)算 邏輯運(yùn)算是 MATLAB 中數(shù)組運(yùn)算所特有的一種運(yùn)算形式，也是幾乎所有的高級語言普遍適用的一種運(yùn)算。 (4)對地阻抗矩陣 X 矩陣 X每一行由下列參數(shù)構(gòu)成： ① 節(jié)點(diǎn)號(hào)； ② 對地阻抗。 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。 Y(p,p)=1./x(i,2)。 q=B1(i,1)。 %對角高壓側(cè) Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。N0=2*n。a=a+1。J(p,N)=DQ。q=q+1。q=2*j11。 J(m,p)=X1。J(m,q)=X2。 X3=X2。J(p,q)=X3。J(m,N)=DP。 elseif j1==i amp。 X4=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。J(p,N)=DQ。q=q+1。 DV=V(i)^2V2。 X5=0。J(p,q)=X5。J(m,N)=DP。 elseif j1==i amp。 。 X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。J(p,q)=X6。m=p+1。 p=2*i1。 j1~=i %PU 節(jié)點(diǎn)非對角 X1=G(i,j1)*e(i)B(i,j1)*f(i)。J(m,q)=X2。 J(m,p)=X1。q=2*j11。 X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。J(p,q)=X4。m=p+1。 p=2*i1。 j1~=i %PQ 節(jié)點(diǎn)非對角 X1=G(i,j1)*e(i)B(i,j1)*f(i)。q=q+1。J(p,N)=DV。J(m,q)=X2。 J(m,p)=X1。q=2*j11。a=0。以下則是此子程序重要部分的編寫： 雅可比矩陣形成的流程圖如圖 32所示 畢業(yè)設(shè)計(jì)（論文） 23 圖 32 雅可比矩陣框圖 雅可比矩陣的程序代碼： ICT1=0。 %非對角 Y(q,p)=Y(p,q)。 q=B1(i,2)。 for i=1:n if x(i,2)~=0。Y(q,p)=Y(p,q)。 1 程序中需要輸入的數(shù)據(jù)： （ 1） n 為節(jié)點(diǎn)數(shù)、 n1為支路數(shù)、 isb 為平衡母線節(jié)點(diǎn)號(hào)、 pr 為誤差精度 畢業(yè)設(shè)計(jì)（論文） 20 （ 2） 輸入由支路參數(shù)形成的矩陣 B1 矩陣 B1的每一行由下列參數(shù)構(gòu)成的： ① 某支路的首段號(hào) P； ② 末端號(hào) Q，且 PQ； ③ 支路的阻抗（ R+jX）； ④ 支路的對地阻抗； ⑤ 支路的變比 K； ⑥ 折算到哪一側(cè)的標(biāo)志（如果支路的首段 P 處于高壓側(cè)則輸入“ 1”，否則請輸入“ 0”）。有了“對應(yīng)元素”的規(guī)定，數(shù)組的運(yùn)算實(shí)質(zhì)上就是針對數(shù)組內(nèi)部的每個(gè)元素進(jìn)行的。而乘除法運(yùn)算有相當(dāng)大的區(qū)別，數(shù)組的乘除法是指兩同維數(shù)組對應(yīng)元素之間的乘除法，它們的運(yùn)算符為“ .*”和“ ./”或“ .\”。但需注意進(jìn)行數(shù)除時(shí)，常數(shù)通常只能做除數(shù)。在傳統(tǒng)的 MATLAB 算法中，右除是先計(jì)算矩陣的逆再相乘，而左除則不需要計(jì)算逆矩陣直接進(jìn)行除運(yùn)算。 原始數(shù)據(jù)輸入格式的設(shè)計(jì)，主要應(yīng)從使用的角度出發(fā)，原則是簡單明了，便于修改。通過 M 語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時(shí)間，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。 MATLAB 是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計(jì)語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運(yùn)算，同時(shí)在數(shù)值分析、自動(dòng)控制模擬、數(shù)字信號(hào)處理、動(dòng)態(tài)分析、繪圖等方面也具有強(qiáng)大的功能。如果不能滿足要求，則應(yīng)修改某些變量的給定值，甚至修改系統(tǒng)的運(yùn)行方式，重新進(jìn)行計(jì)算。 ③