【導(dǎo)讀】基于牛頓拉夫遜法潮流計(jì)算的matlab實(shí)驗(yàn)報(bào)告。應(yīng)用MATLAB語言編寫具有一定通用性的牛。潮流計(jì)算方法為牛頓-拉夫遜法。程序具有較強(qiáng)通用性。形成節(jié)點(diǎn)導(dǎo)納矩陣。設(shè)節(jié)點(diǎn)電壓iief，i=1,2?按式（3-3），（3-4）計(jì)算雅克比矩陣元素。求解修正方程式，得()kie?計(jì)算平衡節(jié)點(diǎn)及PV節(jié)點(diǎn)功率。IEEE-美國新英格蘭10機(jī)39節(jié)點(diǎn)測試系統(tǒng)。1）系統(tǒng)容量基準(zhǔn)值為100MVA。2)負(fù)荷數(shù)據(jù)見表D-1. 節(jié)點(diǎn)號(hào)有功無功節(jié)點(diǎn)號(hào)有功無功。平衡節(jié)點(diǎn)BUS-31（幅值）

　　

【正文】 Zt~=0 Yt=1/Zt。 end Ym=line(k,5)+j*line(k,6)。 K=line(k,7)。 if (K==0)amp。(J~=0) % 普通線路 : K=0 YtYm(k,3) = Yt。 YtYm(k,4) = Ym。 YtYm(k,5) = Ym。 end if (K==0)amp。(J==0) % 對(duì)地支路 : K=0,J=0,R=X=0 YtYm(k,4) = Ym。 end if K0 % 變壓器線路 : Zt和 Ym為折算到 i側(cè)的值 ,K在 j側(cè) YtYm(k,3) = Yt/K。 YtYm(k,4) = Ym+Yt*(K1)/K。 YtYm(k,5) = Yt*(1K)/K/K。 end if K0 % 變壓器線路 : Zt和 Ym為折算到 K側(cè)的值 ,K在 i側(cè) YtYm(k,3) = Yt*K。 YtYm(k,4) = Ym+Yt*(1+K)。 YtYm(k,5) = Yt*(K^2+K)。 end end 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 計(jì)算并返回節(jié)點(diǎn)數(shù)據(jù)結(jié)果 function bus_res = bus_res(bus) [nb,mb]=size(bus)。 bus_res = zeros(nb,4)。 % bus_res矩陣儲(chǔ)存著節(jié)點(diǎn)計(jì)算結(jié)果 bus_res(:,1:2) = bus(:,1:2)。 bus_res(:,3) = bus(:,3) *180 / pi。 % 相角采用角度制 bus_res(:,4) = bus(:,4) + (sqrt(1))*bus(:,5)。 % 注入功率 end 計(jì)算并返回線路潮流 function S_res = S_res(bus,line,YtYm) [nl,ml]=size(line)。 S_res = zeros(nl,5)。 % S_res矩陣儲(chǔ)存著線路潮流計(jì)算結(jié)果 S_res(:,1:2) = line(:,1:2)。 % 前兩列為節(jié)點(diǎn)編號(hào) for k=1:nl I = S_res(k,1)。 J = S_res(k,2)。 if (J~=0)amp。(I~=0) S_res(k,3)=bus(I,2)^2*(conj(YtYm(k,3))+conj(YtYm(k,4)))bus(I,2)*bus(J,2)*(cos(bus(I,3))+j*sin(bus(I,3)))*(conj(cos(bus(J,3))+j*sin(bus(J,3))))*conj(YtYm(k,3))。 S_res(k,4)=bus(J,2)^2*(conj(YtYm(k,3))+conj(YtYm(k,5)))bus(I,2)*bus(J,2)*(conj(cos(bus(I,3))+j*sin(bus(I,3))))*(cos(bus(J,3))+j*sin(bus(J,3)))*conj(YtYm(k,3))。 S_res(k,5)=S_res(k,3) + S_res(k,4)。 % 利用公式計(jì)算非接地支路的潮流 else if(J==0 ) S_res(k,3)=bus(I,2)^2*conj(YtYm(k,4))。 S_res(k,5)=S_res(k,3)+S_res(k,4)。 else S_res(k,4)=bus(J,2)^2*conj(YtYm(k,5))。 S_res(k,5)=S_res(k,3)+S_res(k,4)。 % 利用公式計(jì)算接地支路的潮流 end end end end 作業(yè) y=pifun(10000) 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 y = clear all tic。 y=pifun(100000)。 toc。 Elapsed time is seconds. function y= pifun(n) y=0。 for i=1:n y=y+(1)^(i1)*1/(2*i1)。 end y=y*4。 end a=0:pi/12:2*pi。 x=5*cos(a)。 y=3*sin(a)。 plot(x,y) 5 4 3 2 1 0 1 2 3 4 53210123 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 clear format long max1=100。 eps1=。 eps2=。 x0=[ ]39。 x=x0。 for i=1:max1, A=Fd(x)。 b=F(x)。 dx=A\b。 %矩陣左除，即 b 除以 A x=xdx。 sprintf(39。第 %d 次迭代結(jié)果： \n39。,i) [x dx F(x)] %在屏幕上輸出每次的 x(i),dx(i),F(x(i)) if (max(abs(dx))eps1)amp。(max(abs(F(x)))eps2) break end end ans = 第 1 次迭代結(jié)果： ans = ans = 第 2 次迭代結(jié)果： ans = 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 ans = 第 3 次迭代結(jié)果： ans = ans = 第 4 次迭代結(jié)果： ans = ans = 第 5 次迭代結(jié)果： ans = ans = 第 6 次迭代結(jié)果： 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 ans = ans = 第 7 次迭代結(jié)果： ans = ans = 第 8 次迭代結(jié)果： ans = 0 0 x = function J=Fd(x) J=[3,x(3)*sin(x(2)*x(3)),x(2)*sin(x(2)*x(3))。 2*x(1), 162*(x(2)+),cos(x(3))。 0,x(3)*exp(x(2)*x(3)),x(2)*exp(x(2)*x(3))+20]。 end 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 function y=F(x) y(1)=3*x(1)cos(x(2)*x(3))。 y(2)=x(1)^281*(x(2)+)^2+sin(x(3))+。 y(3)=exp(x(2)*x(3))+20*x(3)+(10*pi3)/3。 y=[y(1)。y(2)。y(3)]。 end format long max1=100。 eps=10^10。 x0=[ ]39。 x=x0。 for i=1:max1, A=Fd(x)。 b=F(x)。 dx=A\b。 %矩陣左除，即 b 除以 A x=xdx。 sprintf(39。第 %d 次迭代結(jié)果： \n39。,i) [x dx F(x)] %在屏幕上輸出每次的 x(i),dx(i),F(x(i)) if (max(abs(dx))eps)amp。(max(abs(F(x)))eps) break end end ans = 第 1 次迭代結(jié)果： ans = +002 * ans = 第 2 次迭代結(jié)果： 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 ans = ans = 第 3 次迭代結(jié)果： ans = ans = 第 4 次迭代結(jié)果： ans = ans = 第 5 次迭代結(jié)果： ans = 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 ans = 第 6 次迭代結(jié)果： ans = 0 0 0 x = 1 5 4 function y=F(x) y(1)=x(1)5*x(2)^2+7*x(3)^2+12。 y(2)=3*x(1)*x(2)+x(1)*x(3)11*x(1)。 y(3)=2*x(2)*x(3)+40*x(1)。 y=[y(1)。y(2)。y(3)]。 end function J=Fd(x) J=[1,10*x(2),14*x(3)。 3*x(2)+x(3)11, 3*x(1),x(1)。 40,2*x(3),2*x(2)]。 end