【文章內(nèi)容簡(jiǎn)介】 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為不平衡量列向量% 求雅克比矩陣% 當(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è)新的層次。第四篇：MATPOWER在電力系統(tǒng)潮流計(jì)算中的應(yīng)用,MATPOWERMATPOWER在電力系統(tǒng)潮流計(jì)算中的應(yīng)用MATPOWER是一個(gè)用MATLAB的M文件編寫,用來解決電力系統(tǒng)潮流計(jì)算和優(yōu)化潮流計(jì)算問題的軟件包。MATPOWER特點(diǎn)是簡(jiǎn)單、易懂且程序代碼公開，這為電力系統(tǒng)專業(yè)學(xué)生深入學(xué)習(xí)和理解掌握潮流計(jì)算中的難點(diǎn)（如節(jié)點(diǎn)導(dǎo)納矩陣、算法及迭代過程等）提供了一個(gè)開放、便捷的平臺(tái)。，僅供學(xué)習(xí)電力系統(tǒng)專業(yè)的讀者參考。讀者利用MATPOWER進(jìn)行潮流計(jì)算時(shí)。圖1 2機(jī)5節(jié)點(diǎn)系統(tǒng)對(duì)圖1的2機(jī)5節(jié)點(diǎn)系統(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當(dāng)采用牛頓—，在MATLAB的名令窗口輸入以下命令即可 runpf(39。case2_539。)計(jì)算輸出結(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