【文章內(nèi)容簡介】 end% 求有功、無功功率不平衡量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為不平衡量列向量% 求雅克比矩陣% 當i~=j時，求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% 當i=j時，求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é)點電壓新值，準備下一次迭代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é)點注入功率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é)點注入電流I=Y*U39。四運行結(jié)果節(jié)點導(dǎo)納矩陣經(jīng)過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點電壓的幅值和相角（弧度數(shù)）節(jié)點注入功率和電流五 結(jié)果分析在這次學(xué)習(xí)和實際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機算法中的牛頓拉夫遜法進行深入的研讀，弄明白了其原理、計算過程、公式推導(dǎo)以及設(shè)計流程。牛頓拉夫遜法是求解非線性方程的迭代過程，其計算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。利用x(*)=x(k+1)+DX(k+1)進行多次迭代，通過迭代判據(jù)得到所需要的精度值即準確值x(*)。六 結(jié)論通過這個任務(wù)，自己在matlab編程，潮流計算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識儲備不足，對matlab的性能和特點還不能有一個全面的把握，對word軟件也不是很熟練，相信通過以后的學(xué)習(xí)能彌補這些不足，達到一個新的層次。第四篇：MATPOWER在電力系統(tǒng)潮流計算中的應(yīng)用,MATPOWERMATPOWER在電力系統(tǒng)潮流計算中的應(yīng)用MATPOWER是一個用MATLAB的M文件編寫,用來解決電力系統(tǒng)潮流計算和優(yōu)化潮流計算問題的軟件包。MATPOWER特點是簡單、易懂且程序代碼公開，這為電力系統(tǒng)專業(yè)學(xué)生深入學(xué)習(xí)和理解掌握潮流計算中的難點（如節(jié)點導(dǎo)納矩陣、算法及迭代過程等）提供了一個開放、便捷的平臺。，僅供學(xué)習(xí)電力系統(tǒng)專業(yè)的讀者參考。讀者利用MATPOWER進行潮流計算時。圖1 2機5節(jié)點系統(tǒng)對圖1的2機5節(jié)點系統(tǒng)（具體參數(shù)參見：現(xiàn)代電力系統(tǒng)分析王錫凡），按“version 2”：functionmpc=case2_5 % MATPOWER Case Format:Version 2 =39。239。%%Power Flow Data%% %% system MVA base =100。%% Bus data% bus_i type PdQdGsBs area VmVabaseKV zone VmaxVmin =[ 1 1 160 80 0 0 1 1 0 100 1 。2 1 200 100 0 0 1 1 0 100 1 。3 1 370 130 0 0 1 1 0 100 1 。4 2 0 0 0 0 1 0 100 1 。5 3 0 0 0 0 1 0 100 1 。]。%% generator data% bus PgQgQmaxQmin Vg mbase status PmaxPmin =[ 4 500 0 999909999 100 1 600 0。5 0 0 999909999 100 1 600 0。]。%% branch data% fbustbus r x b rataArataBrataC ratio angle status angminangmax =[ 2 1 0 0 0 0 0 1360 360。3 1 0 0 0 0 0 0 1360 360。3 2 0 0 0 0 0 1360 360。3 5 0 0 0 0 0 0 1360 360。2 4 0 0 0 0 0 0 1360 360。]。Return當采用牛頓—，在MATLAB的名令窗口輸入以下命令即可 runpf(39。case2_539。)計算輸出結(jié)果如下：MATPOWER Version , 07Feb2011AC Power Flow(Newton)Newton39。s method power flow converged in 5 in seconds ======================= System Summary ======================= How many? How much? P(MW)Q(MVAr)Buses 5 Total Gen Capacity to Generators 2 Online Capacity to Committed Gens2Generation Loads 3 Load Fixed 3 Fixed Dispatchable0 of Shunts 0 Shunt(inj) Branches 5 Losses(I^2 * Z) Transformers 2 Branch Charging(inj)—— Interties 0 Total Intertie Flow Areas 1 Minimum Maximu