【正文】 J(p,q)=X2。j=jph m=p+1。J(m,q)=X4。De(k)=DX(i+1)。Df(p)=P(i)P1。for i=1:ma(i)=conj((v(B1(i1))/B1(i5)v(B1(i2))/B1(i3))。Q=imag(S)。結(jié)束語第五篇：電力系統(tǒng)潮流計(jì)算程序設(shè)計(jì)電力系統(tǒng)潮流計(jì)算程序設(shè)計(jì)姓名：韋應(yīng)順學(xué)號(hào)：2011021052 電力工程學(xué)院牛頓—拉夫遜潮流計(jì)算方法具有能夠?qū)⒎蔷€性方程線性化的特點(diǎn)，而使用MATLAB語言是由于MATLAB語言的數(shù)學(xué)邏輯強(qiáng)，易編譯。④利用matlab軟件，進(jìn)行上述各種情況潮流的計(jì)算及仿真。實(shí)驗(yàn)要求：。0j206。直流潮流模型是把非線性電力潮流問題簡(jiǎn)化為線形電路問題，從而使分析計(jì)算非常方便，直流潮流專門用于研究電網(wǎng)中有功潮流的分布。S% 求節(jié)點(diǎn)注入電流I=Y*U39。endendend% 當(dāng)i=j時(shí)，求H,N,M,Lfor i=1:4for j=1:5if i~=j H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。k=0。% 求節(jié)點(diǎn)導(dǎo)納矩陣中互導(dǎo)納for i=1:5for j=1:5if i~=jY(i,j)=y(i,j)。牛頓迭代法其最大優(yōu)點(diǎn)是在方程f(x)=0的單根附近時(shí)誤差將呈平方減少，而且該法還可以用來求方程的重根、復(fù)根。JS[k]=j。if(fp!=NULL)fclose(fp)。sL[i].NumI=i2。float YG[NBUS][NBUS],YB[NBUS][NBUS]。int Type。5． 同樣的符號(hào)不能既定義成單獨(dú)變量，又定義成數(shù)組，這樣程序會(huì)出錯(cuò)。3． 掌握子函數(shù)的聲明，定義，調(diào)用的方法。下面以形成導(dǎo)納矩陣為例，說明數(shù)據(jù)的讀寫。FILE *fpout。sB[i].GenP=d3。l/* Check the Y matrix */ if((fp=fopen(“”,“w”))==NULL){ printf(“Can not open the file named 39。} void Gauss(float a[NBUS][NBUS], float b[NBUS], int n)/*定義高斯法 */ { int JS[NBUS]。j用追趕法，解方程組，求將結(jié)果寫到TXT第二篇：電力系統(tǒng)潮流計(jì)算南 京 理 工 大 學(xué)《電力系統(tǒng)穩(wěn)態(tài)分析》課程報(bào)告姓名XX學(xué) 號(hào)： 5*** 自動(dòng)化學(xué)院 電氣工程基于牛頓拉夫遜法的潮流計(jì)算例題編程報(bào)學(xué)院(系): 專業(yè): 題目: 任課教師 碩士導(dǎo)師 告楊偉 XX2015年6月10號(hào)基于牛頓拉夫遜法的潮流計(jì)算例題編程報(bào)告摘要：電力系統(tǒng)潮流計(jì)算的目的在于：確定電力系統(tǒng)的運(yùn)行方式、檢查系統(tǒng)中各元件是否過壓或者過載、為電力系統(tǒng)繼電保護(hù)的整定提供依據(jù)、為電力系統(tǒng)的穩(wěn)定計(jì)算提供初值、為電力系統(tǒng)規(guī)劃和經(jīng)濟(jì)運(yùn)行提供分析的基礎(chǔ)。y(1,4)=1/(+)。P=real(S)。endoP(i)=oP(i)+P(i)。endfor i=1:4a(i)=a(i)+oa(i)。為了符合這些要求，有時(shí)可以降低計(jì)算精度。iPij=229。用式（4）不需要迭代就可以求出節(jié)點(diǎn)電壓相角，再用式（2）計(jì)算各支路的有功潮流，這就是直流潮流的解算過程。（圖1）②發(fā)電廠資料：③變電所資料：1）變電所4低壓母線的電壓等級(jí)分別為：10kv，35kv，10kv，35kv。表5為調(diào)節(jié)前后各節(jié)點(diǎn)的電壓標(biāo)幺值。for i=1：n End e(i)=real(B2(i3))。end v=e+f*j。endP1=C(i)*e(i)+D(i)*f(i)。DX(no+1)=0。ifJ=amp。m=p+1。P=2*i1。2007，35(1)：8184．[11] 方思立，劉增煌，孟慶和．大型汽輪發(fā)電機(jī)自并勵(lì)勵(lì) 磁系統(tǒng)的應(yīng)用條件【J]．中國(guó)電力，1994，27(12)：61．63． FANG Si．Ii，LIU Zenghu鋤g，MENG Qin爭(zhēng)hc．m application conditions of large turbine generator selfexcitation system【J】．Electric Powef，1994，27(12)： 61．63．[12]梁小冰，黃方能．利用EMTDC進(jìn)行長(zhǎng)持續(xù)時(shí)間過程 的仿真研究【J】．電網(wǎng)技術(shù)，2002，26(9)：55．57． LIANG Xiaobing，HUANG Fan爭(zhēng)眥ng．How to cany out simulalion of long dul‘a(chǎn)tion processes by use of EMTDC【J】．Power System 11echnology，2002，26(9)： 5557．[13]王卉，陳楷，彭哲，等．?dāng)?shù)字仿真技術(shù)在電力系統(tǒng)中 的應(yīng)用及常用的幾種數(shù)字仿真工具【J】．繼電器，2004，32(21)：7l一75．wANG Hui，CHEN Kai，PENG zhe，et a1．Application of digital simulation眥hniques棚d severaJ simulation tools in power system[J】．Relay，2004，32(21)：71X4= C(i)+G(i,j)*f(i)B(i,j)*e(i)。J(p,q)=X1。case 1 J(p,q)=X1。x1=length(DX1)。For j=i:nC(i)=C(i)+G(i,j)*e(j)B(i,j)*f(j)。e=e+De；f=f+Df。if A（i2）=0B=imag(Y)。由matlab編程調(diào)節(jié)后，可得到表4的發(fā)電廠電壓和變壓器分接頭電壓得標(biāo)幺值。潮流計(jì)算可以選擇電力系統(tǒng)的運(yùn)行方式，便于定期對(duì)電力系統(tǒng)中的元件進(jìn)行檢修。式（4）為直流潮流方程，因?yàn)楹雎粤私拥刂罚瑫r(shí)忽略了支路電阻，所以沒有有功功率損耗。j206。電力系統(tǒng)潮流的計(jì)算和分析是電力系統(tǒng)運(yùn)行和規(guī)劃工作的基礎(chǔ)。% x1為所求△x的列向量% 求節(jié)點(diǎn)電壓新值，準(zhǔn)備下一次迭代for i=1:4oa(i)=x1(i)。endend% 求有功、無功功率不平衡量for i=1:4for j=1:5oP(i)=oP(i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。S(4)=。y(4,5)=0。i){ t=。int NN=3。i電力系統(tǒng)潮流上機(jī)課程設(shè)計(jì)教學(xué)指導(dǎo)書/*Make Y Matrix*/ for(i=1。sB[i].Volt=d1。int L。（我提供的四個(gè)文件，；導(dǎo)納陣形成文件需要的頭文件,，潮流計(jì)算的原始數(shù)據(jù)；,用高斯消去法解線方程組。掌握“*.txt”文件的打開，讀寫，關(guān)閉的函數(shù)。7． 文件的讀寫目錄用“”表示。struct Shunt { int Num,NumI。d1=d2=d3=d4=d5=d6=ps=qs=。sL[i].R=d1。重點(diǎn)學(xué)習(xí)如何聲明和調(diào)用子程序。} } if(L==0)break。為了便于用迭代法解方程組，需要將上述功率方程改寫成功率平衡方程，并對(duì)功率平衡方程求偏導(dǎo)，得出對(duì)應(yīng)的雅可比矩陣，給未知節(jié)點(diǎn)賦電壓初值，將初值帶入功率