【文章內(nèi)容簡(jiǎn)介】 次迭代，直至達(dá)到精度退出循環(huán)； 7)計(jì)算平衡節(jié)點(diǎn)的功率和線路功率，輸出最后計(jì)算結(jié)果； ① 公式推導(dǎo)② 流程圖三matlab編程代碼clear。% 如圖所示1，2，3，4為PQ節(jié)點(diǎn),5為平衡節(jié)點(diǎn)y=0。% 輸入原始數(shù)據(jù)，求節(jié)點(diǎn)導(dǎo)納矩陣y(1,2)=1/(+)。y(4,5)=0。y(1,3)=1/(+)。y(1,4)=1/(+)。y(1,5)=1/(+)。y(2,3)=1/(+)。y(2,5)=1/(+)。y(3,4)=1/(+)。for i=1:5for j=i:5y(j,i)=y(i,j)。endendY=0。% 求節(jié)點(diǎn)導(dǎo)納矩陣中互導(dǎo)納for i=1:5for j=1:5if i~=jY(i,j)=y(i,j)。endendend% 求節(jié)點(diǎn)導(dǎo)納矩陣中自導(dǎo)納for i=1:5Y(i,i)=sum(y(i,:))。endY% Y為導(dǎo)納矩陣G=real(Y)。B=imag(Y)。% 輸入原始節(jié)點(diǎn)的給定注入功率S(1)=+。S(2)=。S(3)=。S(4)=。S(5)=0。P=real(S)。Q=imag(S)。% 賦初值，U為節(jié)點(diǎn)電壓的幅值，a為節(jié)點(diǎn)電壓的相位角U=ones(1,5)。U(5)=。a=zeros(1,5)。x1=ones(8,1)。x2=ones(8,1)。k=0。while max(x2)1e6for i=1:4for j=1:4H(i,j)=0。N(i,j)=0。M(i,j)=0。L(i,j)=0。oP(i)=0。oQ(i)=0。endend% 求有功、無(wú)功功率不平衡量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)))。oQ(i)=oQ(i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。endoP(i)=oP(i)+P(i)。oQ(i)=oQ(i)+Q(i)。endx2=[oP,oQ]39。% 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)))。N(i,j)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。L(i,j)=H(i,j)。M(i,j)=N(i,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)))。N(i,i)=N(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。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)))。L(i,i)=L(i,i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))endendN(i,i)=N(i,i)2*(U(i))^2*G(i,i)。L(i,i)=L(i,i)+2*(U(i))^2*B(i,i)。endJ=[H,N。M,L]% J為雅克比矩陣x1=((inv(J))*x2)。% x1為所求△x的列向量% 求節(jié)點(diǎn)電壓新值，準(zhǔn)備下一次迭代for i=1:4oa(i)=x1(i)。oU(i)=x1(i+4)*U(i)。endfor i=1:4a(i)=a(i)+oa(i)。U(i)=U(i)+oU(i)。endk=k+1。endk,U,a% 求節(jié)點(diǎn)注入功率i=5。for j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))+P(i)。Q(i)=U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))+Q(i)。endS(5)=P(5)+Q(5)*sqrt(1)。S% 求節(jié)點(diǎn)注入電流I=Y*U39。四運(yùn)行結(jié)果節(jié)點(diǎn)導(dǎo)納矩陣經(jīng)過(guò)五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點(diǎn)電壓的幅值和相角（弧度數(shù)）節(jié)點(diǎn)注入功率和電流五 結(jié)果分析在這次學(xué)習(xí)和實(shí)際操作過(guò)程里：首先，對(duì)電力系統(tǒng)分析中潮流計(jì)算的部分特別是潮流計(jì)算的計(jì)算機(jī)算法中的牛頓拉夫遜法進(jìn)行深入的研讀，弄明白了其原理、計(jì)算過(guò)程、公式推導(dǎo)以及設(shè)計(jì)流程。牛頓拉夫遜法是求解非線性方程的迭代過(guò)程，其計(jì)算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。利用x(*)=x(k+1)+DX(k+1)進(jìn)行多次迭代，通過(guò)迭代判據(jù)得到所需要的精度值即準(zhǔn)確值x(*)。六 結(jié)論通過(guò)這個(gè)任務(wù)，自己在matlab編程，潮流計(jì)算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問(wèn)題：理論知識(shí)儲(chǔ)備不足，對(duì)matlab的性能和特點(diǎn)還不能有一個(gè)全面的把握，對(duì)word軟件也不是很熟練，相信通過(guò)以后的學(xué)習(xí)能彌補(bǔ)這些不足，達(dá)到一個(gè)新的層次。第三篇：電力系統(tǒng)潮流計(jì)算程序設(shè)計(jì)電力系統(tǒng)潮流計(jì)算程序設(shè)計(jì)姓名：韋應(yīng)順學(xué)號(hào)：2011021052 電力工程學(xué)院牛頓—拉夫遜潮流計(jì)算方法具有能夠?qū)⒎蔷€性方程線性化的特點(diǎn)，而使用MATLAB語(yǔ)言是由于MATLAB語(yǔ)言的數(shù)學(xué)邏輯強(qiáng)，易編譯?！尽俊尽縁unction tisco %這是一個(gè)電力系統(tǒng)潮流計(jì)算的程序 n=input（‘n請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=’）； m=input(‘請(qǐng)輸入支路數(shù)：m=’)。ph=input(‘n請(qǐng)輸入平衡母線的節(jié)點(diǎn)號(hào)：ph=’)； B1=input（‘n請(qǐng)輸入支路信號(hào)：B1=’）。%它以矩陣形式存貯支路的情況，每行存貯一條支路 %第一列存貯支路的一個(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)的類(lèi)型：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為修正值)。Y=zeros（n）；Y(p,q)=Y(p,q)1./(B1(i3)*B1(i5)。e=zeros（1，n）；Y(p,q)=Y(p,q)。f=zeros（1，n）；no=2*ph=1； Y(q,q)=Y(q,q)+1./B1(i3)+B1(i4)/2；End for i=1：nG=real(Y)。if A（i2）=0B=imag(Y)。p=A(i1)；Y(p p)=1./A(i2)。for i=1：n End e(i)=real(B2(i3))。End f(i)=imag(B2(i3))。For i=1:m S(i)=B2(i1)B2(i2)。p=B1(i1)。V(i)=B2(i3)。p=B1(i2)。end Y(p,p)=Y(p,p)+1./(B1(i3)*B1(i5)^2+B1(i4)./2P=real(S)。Q=imag(S)。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)。J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)。[De,Di]=hxf(J,D,F,ph,n,no)。t=0。whilemax(abs(De))ipamp。max(abs(Dfi)ipt=t+1。e=e+De；f=f+Df。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)；J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)；[De,Df]=hxf(J,Df,ph,n,no)。end v=e+f*j。for i=1:n hh(i)=conj(Y(ph,i)*v(i))。end S(ph)=sum(hh)*v(ph)。B2(ph,1)=S(ph)。V=abs(v)。jd=angle(v)*180/p。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))]。for i=1:ma(i)=conj((v(B1(i1))/B1(i5)v(B1(i2))/B1(i3))。b(i)=v(B1(i1))*a(i)j*B1(i4)*v(B1(i))^2/2。c(i)=v(B1(i2))*a(i)j*B1(i4)*v(B1(i2))^2/2。end result2=[B1(:,6),B1(:,1),B1(:,2),real(b’),imag(b’),real(c’),imag(c’), real(b’+c’),imag(b’+c’)]。printcut(result1,S,b,c,result2)。type resultm function [C,D,Df]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)%該子程序是用來(lái)求取Df for i=1:nIfi=phC(i)=0。D(i)=0。For j=i:nC(i)=C(i)+G(i,j)*e(j)B(i,j)*f(j)。D(i)=D(i)+G(i,j)*f(j)+B(i,j)*e(j)。endP1=C(i)*e(i)+D(i)*f(i)。Q1=C(i)*f(i)D(i)*e(i)。V1=e(i)^2+f(i)^2。IfB2(i4)=2 p=2*i1。Df(p)=P(i)P1。p=p+1。else p=2*i1。Df(p)=P(i)P1。p=p+1。Df(p)=Q(i)Q1。end end end Df=Df’。If ph=n Df(noJ=[]。endfunction [De,Df]=hxf(J,Df,ph,n,no)%該子函數(shù)是為求取De Df DX=JDf。DX1=DX。x1=length(DX1)。if ph=n DX(no)=0。DX(no+1)=0。For i=(no+2):(x1+2)DX(i)=DX1(i2)。End ElseDX=[DX1,0,0]。End k=0。[x,y]=size(DX)。For i=1:2:x K=k+1。Df(k)=DX(i)。De(k)=DX(i+1)。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)。X2=G(i,j)*e(i)+B(i,j)*f(i)。%該子程序是用來(lái)求取jacci矩陣for i=1:n X3=0。switch B2(i4)X4=0。case 3 P=2*i1。continue q=2*j1。case 1 J(p,q)=X1。for j=1:n m=p+1。ifJ=amp。J=ph J(m,q)=X3。X1=G(i)*f(i)B(i,j)*e(i)。q=q+1。X2=G(i,j)*e(i)+B(i,j)*f(i)。J(p,q)=X2。X3=X2。J(m,q)=X4。X4=X1。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。p=2*i1。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。q=2*j1。X3=0。J(p,q)=X1。X4=0。m=p+1。P=2*i1。J(p,q)=X2。q=2*j1。J(m,q)=X4。J(p,q)=X1。Else if j=amp。j=jph m=p+1。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。J(m,q)=X3。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。q=q+1。X3= C(i)+G(i,j)*e(i)B(i,j)*f(i)。J(p,q)=X2。X4= C(i)+G(i,j)*f(i)B(i,j)*e(i)。J(m,q)=X4。P=2*i1。end q=2*j1。end J(p,q)=X1。end m=p+1。end J(m,q)=X3。if ph=n q=q+1。J(no:)=[]。J(p,q)=X2。J(no:)=[]。J(m,q)=X4。J(:,no)=[]。End J(:,no)=[]。End2實(shí)例驗(yàn)證 【例題】設(shè)有一系統(tǒng)網(wǎng)絡(luò)結(jié)線見(jiàn)圖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)信