【正文】 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 dP/df X5=0。 p=2*i1。J(p,q)=X5。 m=p+1。J(m,N)=DP。J(p,q)=X6。 elseif j1==iamp。對(duì)角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 dP/df X5=2*e(i)。 p=2*i1。J(p,q)=X5。 m=p+1。J(m,N)=DP。J(p,q)=X6。 end end end end end========= 以上為求雅可比矩陣的各個(gè)元素 ===================== for k=3:N0 N0=2*n （從第三行開始，第一、二行是平衡節(jié)點(diǎn)） k1=k+1。 N=N0+1 即 N=2*n+1擴(kuò)展列△P、△Q for k2=k1:N1 擴(kuò)展列△P、△Q J(k,k2)=J(k,k2)./J(k,k)。 對(duì)角元規(guī)格化 if k~=3 不是第三行 ============================================================ k4=k1。%消去運(yùn)算 end J(k3,k)=0。 end========================================== for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 end else for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 end end end ====上面是用線性變換方式將Jacobi矩陣化成單位矩陣===== for k=3:2:N01 L=(k+1)./2。 修改節(jié)點(diǎn)電壓實(shí)部 k1=k+1。 修改節(jié)點(diǎn)電壓虛部 end 修改節(jié)點(diǎn)電壓 for k=3:N0 DET=abs(J(k,N))。 不滿足要求的節(jié)點(diǎn)數(shù)加1 end end ICT2(a)=IT2。end用高斯消去法解w=J*Vdisp(39。)。disp(39。)。for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 E(k)=e(k)+f(k)*j。各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為(節(jié)點(diǎn)號(hào)從小到大排列)：39。disp(E)。disp(EE)。39。disp(39。)。VV=V*UB。disp(39。)。各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：39。disp(sida)。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。enddisp(39。)。disp(39。)。disp(SS)。39。disp(39。)。q=B1(i,2)。 Siz(i)=Si(p,q)。 Siz(i)=Si(p,q)。 SSi(p,q)=Si(p,q)*SB。S(39。,39。)=39。 disp(ZF)。 disp(39。)。各條支路的末端功率Sj為(順序同您輸入B1時(shí)一致)：39。for i=1:nl p=B1(i,1)。 if B1(i,6)==0 Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。 else Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。 end disp(Sj(q,p))。 ZF=[39。,num2str(q),39。,num2str(p),39。,num2str(SSj(q,p))]。 %disp(SSj(q,p))。39。enddisp(39。)。q=B1(i,2)。 disp(DS(i))。 ZF=[39。,num2str(p),39。,num2str(q),39。,num2str(DDS(i))]。 %disp(DDS(i))。39。endfigure(1)。plot(V)。節(jié)點(diǎn)號(hào)39。ylabel(39。)。subplot(2,2,2)。xlabel(39。)。電壓角度39。grid on。bar(S)。節(jié)點(diǎn)注入有功39。ylabel(39。)。subplot(2,2,4)。xlabel(39。)。支路首端無功39。grid on。變電所4變電所3變電所2變電所135kV母線35kV母線10kV母線10kV母線一次側(cè)電壓220kV一次側(cè)電壓220kV線路長(zhǎng)為90km線路長(zhǎng)為70km線路長(zhǎng)為90km線路長(zhǎng)為80km線路長(zhǎng)為80km線路長(zhǎng)為100km母線2母線1。母線32*QFQ5022*TQN1002TQN10023*QFs502電廠二電廠一發(fā)電廠資料：母線1和2為發(fā)電廠高壓母線，發(fā)電廠一總裝機(jī)容量為350MW，母線3為機(jī)壓母線，機(jī)壓母線上裝機(jī)容量為100MMW，最大負(fù)荷和最小負(fù)荷分別為40MW和20MW；發(fā)電廠二總裝機(jī)容量為200WM。二、 課程設(shè)計(jì)基本內(nèi)容：1. 對(duì)給定的網(wǎng)絡(luò)查找潮流計(jì)算所需的各元件等值參數(shù)，畫出等值電路圖2. 輸入各支路數(shù)據(jù)，各節(jié)點(diǎn)數(shù)據(jù)利用給定的程序進(jìn)行在變電所在某一負(fù)荷情況下的潮流計(jì)算，并對(duì)計(jì)算結(jié)果進(jìn)行分析。1) 4個(gè)變電所的負(fù)荷同時(shí)以2%的比例增大；2) 4個(gè)變電所的負(fù)荷同時(shí)以2%的比例下降3) 1和4號(hào)變電所的負(fù)荷同時(shí)以2%的比例下降，而2和3號(hào)變電所的負(fù)荷同時(shí)以2%的比例上升；4. 在不同的負(fù)荷情況下，分析潮流計(jì)算的結(jié)果，如果各母線電壓不滿足要求，進(jìn)行電壓的調(diào)整。6. 利用DDRTS軟件，進(jìn)行繪制系統(tǒng)圖進(jìn)行上述各種情況潮流的分析，并進(jìn)行結(jié)果的比較。2. 等值電路的計(jì)算過程電壓是衡量電力系統(tǒng)電能質(zhì)量的標(biāo)準(zhǔn)之一。保證用戶的電壓接近額定值是電力系統(tǒng)調(diào)度的基本任務(wù)之一。本文經(jīng)過手算形成了等值電路圖，并編寫好了程序得出節(jié)點(diǎn)電壓標(biāo)幺值，使其滿足所要求的調(diào)整范圍。⑵為了保證整個(gè)系統(tǒng)潮流計(jì)算的完整性，我們把凡具有母線處均選作節(jié)點(diǎn)，這樣，共有10條母線即選定10個(gè)節(jié)點(diǎn)，我們確定發(fā)電廠一母線為平衡節(jié)點(diǎn)，節(jié)點(diǎn)號(hào)為①，發(fā)電廠二母線為PV節(jié)點(diǎn)，節(jié)點(diǎn)號(hào)為②，其余節(jié)點(diǎn)均為PQ節(jié)點(diǎn)，節(jié)點(diǎn)號(hào)見等值電路圖。根據(jù)所給實(shí)際電路圖和題中的已知條件，有以下公式計(jì)算各輸電線路的阻抗和對(duì)地支路電容的標(biāo)幺值和變壓器的阻抗標(biāo)幺值。 ⑷ 計(jì)算各線路參數(shù)：支路阻抗 支路對(duì)地電容 1 3支路80km輸電線路： 1 4支路80km輸電線路：3 4支路90km輸電線路：4 5支路70km同桿架設(shè)雙回線：2 5支路90km同桿架設(shè)雙回線：2 6支路100km銅桿架設(shè)雙回線： ⑸ 計(jì)算各變器參數(shù)： 變壓器1，2： 兩臺(tái)變壓器并聯(lián)，得阻抗變壓器3，4： 電阻的標(biāo)么值為： 電抗的標(biāo)么值為：兩臺(tái)變壓器并聯(lián)，得阻抗為：+以上計(jì)算結(jié)果列表：首端號(hào)末端號(hào)阻抗有名值阻抗標(biāo)幺值電納有名值電納標(biāo)幺值13++14++25++26++34++45+37++0048++0059++00610++00⑹ 假設(shè)我們所采用的變壓器有5個(gè)抽頭，電壓調(diào)節(jié)范圍為2*%， 對(duì)應(yīng)的分接頭開始時(shí)設(shè)變壓器高壓側(cè)主分接頭, 形成節(jié)點(diǎn)參數(shù)矩陣時(shí)，還需要我們計(jì)算出各節(jié)點(diǎn)所接發(fā)電機(jī)功率，節(jié)點(diǎn)負(fù)荷功率（可根據(jù)負(fù)荷有功功率，發(fā)電機(jī)有功功率及已知的功率因數(shù)計(jì)算）變電所1 變電所2 變電所3 變電所4 平衡節(jié)點(diǎn)為1節(jié)點(diǎn)，所設(shè)節(jié)點(diǎn)電壓的初始值為，PV節(jié)點(diǎn)為2節(jié)點(diǎn)，四. 用牛頓——拉夫遜法進(jìn)行潮流計(jì)算分析與DDRTS潮流計(jì)算比較(一) 變電所負(fù)荷為題目所給數(shù)據(jù)進(jìn)行求解 變電所的負(fù)荷分別為70MW 、60MW、50MW、45MW，按照3中的算法框圖，通過潮流計(jì)算程序[1]進(jìn)行參數(shù)的輸入并求解。輸入節(jié)點(diǎn)參數(shù)形成的矩陣B2； 1）節(jié)點(diǎn)所接發(fā)電機(jī)的功率Sg； 2）節(jié)點(diǎn)負(fù)荷的功率Sl； 3）節(jié)點(diǎn)電壓的初始值； 4）PV節(jié)點(diǎn)電壓V的給定值； 5）節(jié)點(diǎn)所接的無功補(bǔ)償設(shè)備的容量； 6）節(jié)點(diǎn)分類號(hào)igl（1—平衡節(jié)點(diǎn)、2—PQ節(jié)點(diǎn)、3—PV節(jié)點(diǎn)）、B2矩陣：根據(jù)所求參數(shù)，以及B1矩陣的含義，列寫B(tài)1矩陣如下：發(fā)電廠一、二的總裝機(jī)容量分別為100 MW和200 MW。并且，根據(jù)前面敘述的節(jié)點(diǎn)分類，形成B2矩陣如下： 由于各節(jié)點(diǎn)電壓不應(yīng)超過各自的允許范圍，由題目要求知，變壓器高壓側(cè)電壓標(biāo)幺值變化范圍如下表：節(jié)點(diǎn)3456電壓下限11電壓上限同時(shí)保證發(fā)電廠發(fā)出的有功功率要大于變電所有功功率，并且兩者差值越小表明線路有功功率損耗越小，越滿足實(shí)際要求。由附錄2的結(jié)果可以觀察到各個(gè)節(jié)點(diǎn)電壓標(biāo)幺值分別為: 節(jié)點(diǎn)12345678910電壓節(jié)點(diǎn)12345678910電壓實(shí)際值230230222224226227節(jié)點(diǎn)12345678910電壓相角0結(jié)果分析：變電所1，~，變電所3的低壓母線電壓與要求的相比，略微偏低，極端電壓已經(jīng)達(dá)到最大，所以，進(jìn)行無功補(bǔ)償，將節(jié)點(diǎn)9的無功功率減少4var.Q=調(diào)節(jié)后的參數(shù)為：節(jié)點(diǎn)12345678910電壓節(jié)點(diǎn)12345678910電壓實(shí)際值2302302242