牛顿法潮流计算

1、潮流计算实例潮流计算实例 #oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);oQ(i)=oQ(i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j);endoP(i)=oP(i)+P(i);oQ(i)=oQ(i)+Q(i);endfx=oP,oQ;%求雅克比矩阵%当1=：时候求H,N,M,L如下：fori=1:4forj=1:4ifi=jH(i,j)=-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j);

2、N(i,j)=-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);L(i,j)=H(i,j);M(i,j)=-N(i,j);endendendH,N,M,L%当1=：时H,N,M,L如下：fori=1:4forj=1:5ifi=jH(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j);N(i,i)=N(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);M(i,i)=M(i,i)-U(i)*U(j)

3、*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);L(i,i)=L(i,i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j);endendN(i,i)=N(i,i)-2*(U(i)J2*G(i,i);L(i,i)=L(i,i)+2*(U(i)J2*B(i,i);endJ=H,N;M,L%J为雅克比矩阵ox=-(inv(J)*fx);fori=1:4oe(i)=ox(i);oU(i)=ox(i+4)*U(i);endfori=1:4e(i)=e(i)+oe(i);U(i)=U(i)+oU(i);end

4、count=count+1;endox,U,e,count%求节点注入的净功率i=5;forj=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j)+P(i);Q(i)=U(i)*U(j)*(G(irj)*sin(e(i)-e(j)-B(irj)*cos(e(i)-e(j)+Q(i);endS(5)=P(5)+Q(5)*sqrt(-1);S%求节点注入电流I=Y*U运行结果Y值：CommandWindow-+i?(pNew忙MATLAB?WatchthisVidec.seeCemcs.orreadCettinqstarted.

5、10.8333-32.500ji-1.6667+E.OOOOi-1.6667+5.OOOOi-2.5000-+7.EQijOi-5.0000+15.00uui-1.6667+5.uuLiiji12.9167-38.75001-10.ijijuu+30.0000i0-1.2500+3.75001-1.6667+5.uuuiji-10.uuuu4-3flL0000i12.9167-38.75001-1.2SO0-+3.750uiu-2.5000十7.5fl00iU-1.2500十3.750Di3.750D-11.2500iu-5.Liijijij+15.00001-1.25004-3.75001u

6、uG.2SO0-18.75001迭代过程:05.00005.00007.50005.0000030,0000a5.000030.000003.75007.500003.7500a-33.10005.0000E.uOOj7.500j-10.53331.66671.66612.500u5.0000-3S.975030.00001.6667-12.341710.COOO5.000030.uuOO-38.75003.75001.6S6710.CuDO-12.91671.25007.5Q0003.7500-11.250?2.5000i.2500-5.750011.1333-1.6667-1.6667-2

7、.5000-31.60uii5.ijijUU5.ijUUO7.50uU-1.666712.9917-10.0000j5.UUULi-38.525030.ijUUUU-1.6667-10.Duuu12.9167-1.250j5l00003u.ijuuu-38.75003.750u-2.50000-1.25003.750j7.5000U3.75OU-11.250005.290815.25127.924S5.1555O30.54260E133730-4267a3.80187.3217a3.744aN-01.54361.51122.13691.9393O9-9S7601.9QES1O.33ES1.工曰

8、口3.0153O131710Til=0-1.5439-1.6112-2.13681.03S3a-S.OS73a-la365S一iQm33550-1.2050-3.0453O1,31710潮流计算实例潮流计算实例潮流计算实例潮流计算实例-34.71365.28995.29077.9237-L1.83701.54361.51112.13685.1580-39.581230.539501.9391-12.69379.986705.139230.4236-39.36423,S0131.965510,3343-12.7U471.20497.G2LOU3.7640-11.38603.044801.3169-3.161711.4379-1.5436-1.5111-2.1368-35.use5.28995.29077.9237-1.339113,5937-9.986705115sl-