【正文】 程，并對(duì)功率平衡方程求偏導(dǎo)，得出對(duì)應(yīng)的雅可比矩陣，給未知節(jié)點(diǎn)賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點(diǎn)電壓不平衡量（未知的）構(gòu)成了誤差方程，解誤差方程，得到節(jié)點(diǎn)電壓不平衡量，節(jié)點(diǎn)電壓加上節(jié)點(diǎn)電壓不平衡量構(gòu)成節(jié)點(diǎn)電壓新的初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計(jì)算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。2 基本步驟和設(shè)計(jì)流程圖形成了雅克比矩陣并建立了修正方程式，運(yùn)用牛頓拉夫遜法計(jì)算潮流的核心問題已經(jīng)解決，已有可能列出基本計(jì)算步驟并編制流程圖。由課本總結(jié)基本步驟如下：1)形成節(jié)點(diǎn)導(dǎo)納矩陣Y；2)設(shè)各節(jié)點(diǎn)電壓的初值，如果是直角坐標(biāo)的話設(shè)電壓的實(shí)部e和虛部f；如果是極坐標(biāo)的話則設(shè)電壓的幅值U和相角a；3)將各個(gè)節(jié)點(diǎn)電壓的初值代入公式求修正方程中的不平衡量以及修正方程的系數(shù)矩陣的雅克比矩陣；4)解修正方程式，求各節(jié)點(diǎn)電壓的變化量，即修正量； 5)計(jì)算各個(gè)節(jié)點(diǎn)電壓的新值，即修正后的值；6)利用新值從第（3）步開始進(jì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% 求有功、無功功率不平衡量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)過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點(diǎn)電壓的幅值和相角（弧度數(shù)）節(jié)點(diǎn)注入功率和電流五 結(jié)果分析在這次學(xué)習(xí)和實(shí)際操作過程里：首先，對(duì)電力系統(tǒng)分析中潮流計(jì)算的部分特別是潮流計(jì)算的計(jì)算機(jī)算法中的牛頓拉夫遜法進(jìn)行深入的研讀，弄明白了其原理、計(jì)算過程、公式推導(dǎo)以及設(shè)計(jì)流程。牛頓拉夫遜法是求解非線性方程的迭代過程，其計(jì)算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。利用x(*)=x(k+1)+DX(k+1)進(jìn)行多次迭代，通過迭代判據(jù)得到所需要的精度值即準(zhǔn)確值x(*)。六 結(jié)論通過這個(gè)任務(wù)，自己在matlab編程，潮流計(jì)算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識(shí)儲(chǔ)備不足，對(duì)matlab的性能和特點(diǎn)還不能有一個(gè)全面的把握，對(duì)word軟件也不是很熟練，相信通過以后的學(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語言是由于MATLAB語言的數(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)的類型：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)%該子程序是用來求取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)P