【正文】 功（ MW） 無功（ Mvar） BUS 3 BUS 4 BUS 7 BUS 8 BUS 12 BUS 15 BUS 16 BUS 18 BUS 20 BUS 21 BUS 23 BUS 24 BUS 25 BUS 26 BUS 27 BUS 28 BUS 29 BUS 39 3）發(fā)電機(jī)數(shù)據(jù)見表 D2 表 D2 發(fā)電機(jī)數(shù)據(jù) 發(fā)電機(jī) 節(jié)點(diǎn)號(hào) 有功（ MW） 無功（ Mvar） 電壓（ .） PV PQ PQ PV PV PQ PV PV PV BUS30 BUS32 BUS33 BUS34 BUS35 BUS36 BUS37 BUS38 BUS39 平衡節(jié)點(diǎn) BUS31 (電壓幅角 ) （幅值） 4）線路參數(shù)見表 D3 表 D3 線路參數(shù) 序號(hào) 節(jié)點(diǎn) I 節(jié)點(diǎn) J R(.) X(.) B/2(.) LN1 LN2 LN3 LN4 LN5 LN6 LN7 LN8 LN9 BUS2 BUS39 BUS3 BUS25 BUS4 BUS18 BUS5 BUS14 BUS6 BUS1 BUS1 BUS2 BUS2 BUS3 BUS3 BUS4 BUS4 BUS5 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 LN10 LN11 LN12 LN13 LN14 LN15 LN16 LN17 LN18 LN19 LN20 LN21 LN22 LN23 LN24 LN25 LN26 LN27 LN28 LN29 LN30 LN31 LN32 LN33 LN34 BUS8 BUS7 BUS11 BUS8 BUS9 BUS39 BUS11 BUS13 BUS14 BUS15 BUS16 BUS17 BUS19 BUS21 BUS24 BUS18 BUS27 BUS22 BUS23 BUS24 BUS26 BUS27 BUS28 BUS29 BUS29 BUS5 BUS6 BUS6 BUS7 BUS8 BUS9 BUS10 BUS10 BUS13 BUS14 BUS15 BUS16 BUS16 BUS16 BUS16 BUS17 BUS17 BUS21 BUS22 BUS23 BUS25 BUS26 BUS26 BUS26 BUS28 LN35: BUS4 接有并聯(lián)電容器， B4＝ LN36: BUS5 接有并聯(lián)電容器， B4＝ 5)變壓器參數(shù)見表 D4 表 D4 變壓器參數(shù) 序號(hào) 節(jié)點(diǎn) I 節(jié)點(diǎn) J R(.) X(.) 變比 (%) TR37 TR38 TR39 TR40 TR41 TR42 TR43 TR44 TR45 TR46 TR47 TR48 BUS11 BUS13 BUS30 BUS31 BUS32 BUS34 BUS33 BUS35 BUS36 BUS37 BUS38 BUS20 BUS12 BUS12 BUS2 BUS6 BUS10 BUS20 BUS19 BUS22 BUS23 BUS25 BUS29 BUS19 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 %IEEE美國新英格蘭 10 機(jī) 39 節(jié)點(diǎn)測試系統(tǒng) % 1 2 3 4 5 6 % bus volt angle p q type bus=[ 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 2 31 3 32 1 33 1 34 2 35 2 36 1 37 2 38 2 39 2]。 （ 2）編程語言為 MATLAB。 4 3 0 0。 8 7 0 0。 17 16 0 0。 24 23 0 0。 11 12 0 0 。 36 23 0 0 。,39。 [nl,ml]=size(line)。計(jì)算結(jié)果 39。 % 生成電 壓對(duì)角矩陣 end end dAngU = J\[dP。 % 計(jì)算每個(gè)節(jié)點(diǎn)的有功和無功注入的子程序 [bus,line] = ReNum(bus,line,nodenum)。)。% 39。 end fprintf(myf,39。,S_res(i,2))。% + j%\n39。 % number of PQ bus counter for i = 1:nb, % nb為總節(jié)點(diǎn)數(shù) type= bus(i,6)。 bus=[PQ。 [nl,ml]=size(line)。 if (K==0)amp。 Y(I,I)=Y(I,I)+Ym。 Y(J,I)=Y(I,J)。 N = zeros(nb1,nPQ)。 end end % 初始化并計(jì)算 Qi和 Pi for i = 1:nb1 for j = 1:nb1 if i~=j H(i,j)=bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)bus(j,3))imag(Y(i,j))*cos(bus(i,3)bus(j,3)))。 end end % 分別計(jì)算 L矩陣的對(duì)角及非對(duì)角元素 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 end end end J = [H N。 bus_temp = zeros(nb,mb)。 % 對(duì) YtYm矩陣賦初值 0 YtYm(:,1:2) = line(:,1:2)。 if (K==0)amp。 YtYm(k,5) = Yt*(1K)/K/K。 % 相角采用角度制 bus_res(:,4) = bus(:,4) + (sqrt(1))*bus(:,5)。 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))。 toc。 eps1=。 sprintf(39。 y(3)=exp(x(2)*x(3))+20*x(3)+(10*pi3)/3。 for i=1:max1, A=Fd(x)。 y(2)=3*x(1)*x(2)+x(1)*x(3)11*x(1)。 end 。 y=[y(1)。 dx=A\b。y(2)。,i) [x dx F(x)] %在屏幕上輸出每次的 x(i),dx(i),F(x(i)) if (max(abs(dx))eps1)amp。 x0=[ ]39。 for i=1:n y=y+(1)^(i1)*1/(2*i1)。 % 利用公式計(jì)算非接地支路的潮流