【正文】 行后各行上三角元素J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。e(L)=e(L)J(k,N)。迭代次數(shù)：39。%計算各節(jié)點電壓的模值sida(k)=atan(f(k)./e(k))*180./pi。)。disp(39。各節(jié)點的功率S為(節(jié)點號從小到大排列)：39。)。SSi(p,q)=Si(p,q)。disp(ZF)。q=B1(i,2)。S(39。39。disp(DS(i))。,num2str(DDS(i))]。請輸入節(jié)點數(shù)：n=39。請輸入平衡母線節(jié)點號：isb=39。+0。+0。00002。Y=zeros(n)。else%左節(jié)點處于K側(cè)p=B1(i,2)。)。%i節(jié)點無功補償量end %===================== P=real(S)。a=a+1。,Q39。% dP/df=dQ/deX3=X2。J(m,N)=DP。% dP/dfX3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。對角元X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。m=p+1。非對角元X1=G(i,j1)*e(i)B(i,j1)*f(i)。%節(jié)點功率P計算eiΣ(Gij*ejBij*fj)+fiΣ(Gij*fj+Bij*ej)Q1=C(i)*f(i)e(i)*D(i)。a=0。%PV節(jié)點電壓給定模值end for i=1:n%給定各節(jié)點注入功率S(i)=B2(i,1)B2(i,2)。%對角元1側(cè)end %求導(dǎo)納矩陣disp(39。% % %for i=1:nl%支路數(shù)if B1(i,6)==0%左節(jié)點處于1側(cè)p=B1(i,1)。請輸入各節(jié)點參數(shù)形成的矩陣： B2=39。00002。+0。B1=[1+0。isb=1。n=12。,num2str(q),39。q=B1(i,2)。disp(ZF)。SSj(q,p)=Sj(q,p)。)。)=39。Siz(i)=Si(p,q)。disp(39。%計算各節(jié)點的注入電流的共軛值endS(p)=E(p)*C(p)。39。%顯示各節(jié)點的實際電壓標(biāo)幺值E用復(fù)數(shù)表示 disp(39。disp(ICT2)。%不滿足要求的節(jié)點數(shù)ICT1=ICT1+1。%消去運算（當(dāng)前行3列元素消為0）end%用當(dāng)前行K2列元素減去當(dāng)前行3列元素乘以第三行K2列元素J(k3,k)=0。% 用K行K列對角元素去除K行K列后的非對角元素進行規(guī)格化endJ(k,k)=1。J(m,q)=X1。% dP/deX2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。J(m,q)=X1。% dP/deX2=B(i,j1)*e(i)G(i,j1)*f(i)。J(p,q)=X4。% dQ/deX4=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。q=q+1。% X2=dp/df X3=dQ/deX4=X1。的計算值V2=e(i)^2+f(i)^2。for i=1:nif i~=isb%非平衡節(jié)點C(i)=0。Q=imag(S)。disp(Y)%G=real(Y)。q=B1(i,1)。e=zeros(1,n)。0+002。+0。)。)。end%本程序的功能是用牛頓——拉夫遜法進行10節(jié)點潮流計算 %本程序的功能是用牛頓——拉夫遜法進行潮流計算 clear。,39。for i=1:nlp=B1(i,1)。,num2str(SSj(q,p))]。enddisp(Sj(q,p))。各條支路的末端功率Sj為(順序同您輸入B1時一致)：39。,num2str(q),39。elseSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))...conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。)。for q=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q))。%顯示各節(jié)點的電壓大小V的模值 disp(39。disp(E)。)。%不滿足要求的節(jié)點數(shù)加1endendICT2(a)=IT2。%當(dāng)前行第k列元素已消為0endelse%是第三行k=3%====================== 第三行k=3的前代運算 ========================for k3=k1:N0%從第四行到2n行（最后一行）for k2=k1:N1%從第四列到2n+1列（即擴展列）J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。% N=N0+1 即 N=2*n+1擴展列△P、△Q 或 △Ufor k2=k1:N1% 從k+1列的Jacobi元素到擴展列的△P、△Q 或 △UJ(k,k2)=J(k,k2)./(J(k,k)+eps)。% PV節(jié)點電壓誤差m=p+1。對角元X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。% PV節(jié)點電壓誤差m=p+1。非對角元X1=G(i,j1)*e(i)B(i,j1)*f(i)。q=q+1。% dP/dfX3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。J(m,N)=DP。% dP/df=dQ/deX3=X2。,Q39。a=a+1。%i節(jié)點無功補償量end %===================== P=real(S)。)。else%左節(jié)點處于K側(cè)p=B1(i,2)。Y=zeros(n)。0+002。+0。B1=[1+0。isb=1。n=11。+180。(+)+()180。14==VdDVL3==功率分布：FLR4：SL4=ZZ+Z*T5*L4*SfT5(VG3Vd)Z+ZL4****180。(+)=()MVA 電壓分布：Sc1=SL2+DSL2=()+(+)=()180。(+)=()MVA352Sg=SfDST61+SG5DST62DST63Si35=(+)MVA計算每一個FLR的功率分布和電壓分布計算如下： FLR1：DVT2=PR+QX40180。Sd=(+)MVADSL339。b0180。r0180。l4=180。180。130=180。40=180。根據(jù)教材查得r0=180。352=180。VN150180。VN=180。10= 322(10180。10=180。180。10=180。180。180。25=50MWQ5=50180。15=60MWQ1=60180。運行的具體結(jié)果和分析也在詳細(xì)設(shè)計中詳述。12MW1180。xd======89kW237。I%=%=,d1116MWVA4180。它既是對電力系統(tǒng)規(guī)劃設(shè)計和運行方式的合理性、可靠性及經(jīng)濟性進行定量分析的依據(jù)，又是電力系統(tǒng)靜態(tài)和暫態(tài)穩(wěn)定計算的基礎(chǔ)?？梢哉f,它是電力系統(tǒng)分析中最基本、最重要的計算,是系統(tǒng)安全、經(jīng)濟分析和實時控制與調(diào)度的基礎(chǔ)。ComputerRomofHighScholinHami, Ham i, Xinjiang83900, China)Acartsb:Withthedevolpmentoftheputerlangseuinrectnyears,thereshcraofM AT LABpower wolfalgohmtirsgotsubstanlidevolpment.Thepaper, usingPQmethod, analyzedthediferenscandimsilarsetiofBASIC, FO RT RANandMAT LABadvaelangseu, andpointedoutt heiradvatnagesand M AT LABpowerflowalgoirthmwords:electirpowersystem。電力系統(tǒng)卷 [ M].第 :中國電力出版社 , 2001.31.社 199.6778 [ 10]趙晉泉 ,侯志儉 , [J ].電力系統(tǒng)自動化 , 1999 , 23(3):32[ 8]陳 [ M].北京 :高等教育出版社 ,20.139193.34.[ 11][英]勞頓 M A,塞 M [M].北京 :機械工業(yè)出版社 , 1992.826827.[ 12]樓順天 ,陳生潭 ,T LAB5.X程序設(shè)計語言 [ M].西安 :西安電子科技大學(xué)出版社 , 2001.5778。＊126西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版)第 32卷 的建立命令 spconvetr,從而將外部數(shù)據(jù)轉(zhuǎn)化為稀疏 矩陣,而且 MATLA B函數(shù) lu也可直接實現(xiàn) LR三 角分解[了方便。U s isi =P 1 ＊YUi =1n∑￣Sij= UiI=Pjiji+Qij平衡節(jié)點功率和線路功率的計算屬復(fù)數(shù)運算。相位角 iΔiδ計算公式為:PQ分解法計算電力系統(tǒng)的潮流分布 ,其試用(k)i步驟為:(B ′)(ΔP(k)/ U(k))=(U式中 , ΔP(k), U(k)為已知量,相位角 Δδ的計算主(1)形成系數(shù)矩陣 B ′, B ″(設(shè)該矩陣為44矩i要(k)Δδ)是矩陣的求積運陣)。其中, N ewton個工程數(shù)學(xué)運算函數(shù),+38[ 15 ]。Uii為節(jié)點消了編譯工作及連接過程,語言的執(zhí)行時間相應(yīng)較i和無功功率。T M744[文獻標(biāo)識碼] A[文章編號] 16719387(204)12012403潮流計算是電力系統(tǒng)規(guī)劃、運行的基本研究方法,隨著現(xiàn)代電力系統(tǒng)大系統(tǒng)、強非線性與多元件特 點的日益突出,其計算量與計算復(fù)雜度急劇增[ 1]加。2哈密市高級中學(xué)計算機室 ,新疆哈密 83900)[摘 要] 隨著計算機語言技術(shù)的不斷發(fā)展和成熟 ,基于 MATLAB的潮流計算研究近年來得到了長足的發(fā) 展。針對這一現(xiàn)狀,本文以某電力網(wǎng)絡(luò)為例,分 析了 BASIC、FORT RAN和 M ATLAB 高級語言潮 流計算的異同,并提出了基于 MAT LAB 的潮流算 法,以期為電力系統(tǒng)潮流計算提供參考。＊Ui一定的可靠性和收斂性。I為節(jié)點 Y為節(jié) 慢。對于大矩陣的潮流計 為接近書寫的數(shù)學(xué)表達格式 ,且程序易調(diào)試。FO RTRAN 語言雖成的系數(shù)矩陣 B′和B″,遠(yuǎn)比 BASIC,FORT RAN簡然在參考資料上提供矩陣求逆、求積程序[ 1618] ,但單,矩陣輸入、輸出與數(shù)學(xué)書寫格式相似。在建立矩陣時, 得較為繁瑣。以雙精度類型進行數(shù)據(jù)的存儲和運算 ,數(shù)據(jù)精確。98.[ 15 ]劉文遠(yuǎn) ,李興成 ,趙SIC語言與農(nóng)業(yè)實用程序 [ 041.[ 16 ]潭浩強 ,崔武子 , RTRAN程序設(shè)計(二級)教程[社 , 123.[ 17 ]鮑有文 ,周海燕 ,趙重敏 , RTRAN77程序設(shè)計試題匯編 [ M].北京 :清華大學(xué)出版[ 20]鄧建中 ,葛仁杰 ,[ M].西安 :西安交通大學(xué)出版社 , 1998.21社 ,.68。本文主要運用的事潮流計算，潮流計算是電力網(wǎng)絡(luò)設(shè)計與運行中最基本的運算，對電力網(wǎng)絡(luò)的各種設(shè)計方案及各種運行方式進行潮流計算，可以得到各種電網(wǎng)各節(jié)點的電壓，并求得網(wǎng)絡(luò)的潮流及網(wǎng)絡(luò)中的各元件的電力損耗，進而求得電能損耗。在電力系統(tǒng)運行方式和規(guī)劃方案的研究中，都需要進行潮流計算以比較運行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟性?？傮w設(shè)計SimpowerSystems使用Simulink環(huán)境，可以將該系統(tǒng)中的發(fā)電機、變壓器，線路等模型聯(lián)結(jié)起來，形成電力系統(tǒng)仿真模擬圖。Dps=50kW237。16MWVA63MWVA236。25MW”236。Dps=121kW237。tan()=156MVar發(fā)電機G3：P3=3180。VN73180。10=180。103=(20180。10DS02=2180。103=(63180。103I0%180。[Vs(12)%+Vs(23)%Vs(13)%]=21Vs3%=180。VN=180。10=180。180。40==x0180。130==x0180。3=2180。l3=2180。60=24WBL4=b0180。20==180。106180。80DST5DS05jQL4DSL4=(+)MVAFLR5: 152DST7=2180。+180。VN=()180。VN=()180。35==(+)：ST631=ST63+DST63=(+)180。對于含有變壓器的支路，第一列為低壓側(cè)節(jié)點編號，第二列為高壓側(cè)節(jié)點編號，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理，第三列為支路的串列阻抗參數(shù)，第四列為支路的對地導(dǎo)納參數(shù)，第五烈為含變壓器支路的變壓器的變比，第六列為變壓器是否是否含有變壓器的參數(shù)，其中“1”為含有變壓器，“0”為不含有變壓器。)。)。+01。00002。0+002。V=zeros(1,n)。%非對角元Y(q,p)=Y(p,q)。%分解出導(dǎo)納陣的實部和虛部for i=1:n%給定各節(jié)點初始電壓的實部和虛部e(i)=real(B2(i,3))。IT2=1。for j1=1:nC(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。%節(jié)