【正文】 mplex(bus[statemp1].V,dd)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。statemp=line[i].start。=0。i{int statemp,endtemp。goto LOOP。//修正節(jié)點(diǎn)電壓// Bus::PrintfNodeV(bus,N)。gauss::gauss_slove(JacAug,x,NN)。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//計(jì)算節(jié)點(diǎn)功率差值Bus::PrintfNodeScal(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。x=new double[NN]。for(i=0。double **JacAug。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。for(i=0。//輸出結(jié)點(diǎn)參數(shù)Complex **X。}Bus::ScanfBusData(bus)。ibus[i].=0。//輸出支路參數(shù)Bus *bus=new Bus[N]。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。N)==0){return 0。//輸入支路個(gè)數(shù)if((Lamp。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。謝謝include “” include include include include“” include“” include“” using namespace std。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)：71End2實(shí)例驗(yàn)證 【例題】設(shè)有一系統(tǒng)網(wǎng)絡(luò)結(jié)線見圖1，各支路阻抗和各節(jié)點(diǎn)功率均已以標(biāo)幺值標(biāo)示于圖1中，其中節(jié)點(diǎn)2連接的是發(fā)電廠，設(shè)節(jié)點(diǎn)1電壓保持U1＝，試計(jì)算其中的潮流分布，請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=5 請(qǐng)輸入支路數(shù)：m=7 請(qǐng)輸入平衡母線的節(jié)點(diǎn)號(hào)：ph=l 請(qǐng)輸入支路信息：BI=[ l 2 0．02+0．06i O l 1；1 3 0．08+0．24i 0 1 2；2 3 0．06+0．18i 0 l 3： 2 4 0．06+0．18i O l 4： 2 5 0．04+0．12i 0 l 5： 3 4 0．01+0．03i 0 l 6： 4 5 0．08+0．24i O 1 7] 請(qǐng)輸入節(jié)點(diǎn)信息：B2=[ 0 0 1．06 3；0．2+0．20i 0 1 1；一O．45一O．15i 0 l l；一0．40．05i 0 l 1；一0．6—0．1i 0 1 l] 請(qǐng)輸入節(jié)點(diǎn)號(hào)及對(duì)地阻抗： A=[l 0；2 0；3 0；4 0；5 O ] 請(qǐng)輸入修正值：ip=0．000 0l參考文獻(xiàn)[1]陳珩．電力系統(tǒng)穩(wěn)定分析[M]．北京：中國(guó)電力出版社，2002：139—187．[2]鄭阿奇．MATLAB實(shí)用教程[M]．北京：電子工業(yè)出版社，2005：1243．[3] 束洪春，孫士云，等．云電送粵交商流混聯(lián)系統(tǒng)全過 程動(dòng)態(tài)電壓研究[J】．中國(guó)電力，2008，4l(10)：l4． SHU Hong—ch吼，SUN Shiyun，et a1．Research on fun prc39。J(:,no)=[]。J(no:)=[]。J(no:)=[]。end J(m,q)=X3。end J(p,q)=X1。P=2*i1。X4= C(i)+G(i,j)*f(i)B(i,j)*e(i)。X3= C(i)+G(i,j)*e(i)B(i,j)*f(i)。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。Else if j=amp。J(m,q)=X4。J(p,q)=X2。m=p+1。J(p,q)=X1。q=2*j1。p=2*i1。X4=X1。X3=X2。X2=G(i,j)*e(i)+B(i,j)*f(i)。X1=G(i)*f(i)B(i,j)*e(i)。ifJ=amp。case 1 J(p,q)=X1。case 3 P=2*i1。%該子程序是用來求取jacci矩陣for i=1:n X3=0。End End case 2 Function for j=1:n J=jacci(Y,G,B,PQ,e,f,V,C,D,B2,n,ph,no)X1=G(i,j)*f(i)B(i,j)*e(i)。Df(k)=DX(i)。[x,y]=size(DX)。End ElseDX=[DX1,0,0]。DX(no+1)=0。x1=length(DX1)。endfunction [De,Df]=hxf(J,Df,ph,n,no)%該子函數(shù)是為求取De Df DX=JDf。end end end Df=Df’。p=p+1。else p=2*i1。Df(p)=P(i)P1。V1=e(i)^2+f(i)^2。endP1=C(i)*e(i)+D(i)*f(i)。For j=i:nC(i)=C(i)+G(i,j)*e(j)B(i,j)*f(j)。type resultm function [C,D,Df]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)%該子程序是用來求取Df for i=1:nIfi=phC(i)=0。end result2=[B1(:,6),B1(:,1),B1(:,2),real(b’),imag(b’),real(c’),imag(c’), real(b’+c’),imag(b’+c’)]。b(i)=v(B1(i1))*a(i)j*B1(i4)*v(B1(i))^2/2。resulte1=[A(:,1),real(v),imag(v),V,jd,real(S’)，imag(S’),real(B2(:1)),imag(B2(:1)),real(B2(:2)),imag(B2(:,2))]。V=abs(v)。end S(ph)=sum(hh)*v(ph)。end v=e+f*j。e=e+De；f=f+Df。whilemax(abs(De))ipamp。[De,Di]=hxf(J,D,F,ph,n,no)。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)。end Y(p,p)=Y(p,p)+1./(B1(i3)*B1(i5)^2+B1(i4)./2P=real(S)。V(i)=B2(i3)。For i=1:m S(i)=B2(i1)B2(i2)。for i=1：n End e(i)=real(B2(i3))。if A（i2）=0B=imag(Y)。e=zeros（1，n）；Y(p,q)=Y(p,q)。%它以矩陣形式存貯支路的情況，每行存貯一條支路 %第一列存貯支路的一個(gè)端點(diǎn) %第二列存貯支路的另一個(gè)端點(diǎn) %第三列存貯支路阻抗%第四列存貯支路的對(duì)地導(dǎo)納%第五列存貯變壓器的變比，注意支路為1 %第六列存貯支路的序號(hào)B2=input(‘n請(qǐng)輸入節(jié)點(diǎn)信息：B2=’)； %第一列為電源側(cè)的功率 %第二列為負(fù)荷側(cè)的功率 %第三列為該點(diǎn)的電壓值%第四列為該點(diǎn)的類型：1為PQ，2為PV節(jié)點(diǎn)，3為平衡節(jié)點(diǎn) A=input（‘n請(qǐng)輸入節(jié)點(diǎn)號(hào)及對(duì)地阻抗：A=’）； ip=input（‘n請(qǐng)輸入修正值：ip=’）； %ip為修正值)。【】【】Function tisco %這是一個(gè)電力系統(tǒng)潮流計(jì)算的程序 n=input（‘n請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=’）； m=input(‘請(qǐng)輸入支路數(shù)：m=’)。六 結(jié)論通過這個(gè)任務(wù)，自己在matlab編程，潮流計(jì)算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識(shí)儲(chǔ)備不足，對(duì)matlab的性能和特點(diǎn)還不能有一個(gè)全面的把握，對(duì)word軟件也不是很熟練，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，達(dá)到一個(gè)新的層次。牛頓拉夫遜法是求解非線性方程的迭代過程，其計(jì)算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。S% 求節(jié)點(diǎn)注入電流I=Y*U39。Q(i)=U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))+Q(i)。endk,U,a% 求節(jié)點(diǎn)注入功率i=5。U(i)=U(i)+oU(i)。oU(i)=x1(i+4)*U(i)。M,L]% J為雅克比矩陣x1=((inv(J))*x2)。L(i,i)=L(i,i)+2*(U(i))^2*B(i,i)。M(i,i)=M(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。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)))。L(i,j)=H(i,j)。% x2為不平衡量列向量% 求雅克比矩陣% 當(dāng)i~=j時(shí)，求H,N,M,Lfor i=1:4for j=1:4if i~=jH(i,j)=U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。oQ(i)=oQ(i)+Q(i)。oQ(i)=oQ(i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。oQ(i)=0。L(i,j)=0。N(i,j)=0。k=0。x1=ones(8,1)。U(5)=。Q=imag(S)。S(5)=0。S(3)=。% 輸入原始節(jié)點(diǎn)的給定注入功率S(1)=