matlab潮流计算.pdf

《matlab潮流计算.pdf》由会员分享，可在线阅读，更多相关《matlab潮流计算.pdf（14页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选文档2 V(i)=B2(i,4); %PV 节点电压给定模值end for i=1:n %给定各节点注入功率S(i)=B2(i,1)-B2(i,2); %i 节点注入功率SG-SL B(i,i)=B(i,i)+B2(i,5); %i 节点无功补偿量end %- P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数ICT1 、a；不满足收敛要求的节点数IT2 while IT2=0 % N0=2*n 雅可比矩阵的阶数；N=N0+1 扩展列IT2=0;a=a+1; for i=1:n if i=i

2、sb %非平衡节点C(i)=0;D(i)=0; for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);% (Gij*ej-Bij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);% (Gij*fj+Bij*ej)end P1=C(i)*e(i)+f(i)*D(i);%节点功率 P计算 ei (Gij*ej-Bij*fj)+fi(Gij*fj+Bij*ej) Q1=C(i)*f(i)-e(i)*D(i);%节点功率 Q 计算 fi(Gij*ej-Bij*fj)-ei(Gij*fj+Bij*ej) %求 i 节点有功和

3、无功功率P,Q的计算值V2=e(i)2+f(i)2; %电压模平方%以下针对非PV 节点来求取功率差及Jacobi 矩阵元素 - if B2(i,6)=3 %非 PV 节点DP=P(i)-P1; %节点有功功率差DQ=Q(i)-Q1; %节点无功功率差%以上为除平衡节点外其它节点的功率计算- %求取 Jacobi 矩阵 - for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de=-dQ/df X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df=dQ/de X3=X2; %

4、X2=dp/df X3=dQ/de X4=-X1; % X1=dP/de X4=dQ/df p=2*i-1;q=2*j1-1; J(p,q)=X3;J(p,N)=DQ;m=p+1; % X3=dQ/de J(p,N)=DQ节点无功功率差J(m,q)=X1;J(m,N)=DP;q=q+1; % X1=dP/de J(m,N)=DP节点有功功率差J(p,q)=X4;J(m,q)=X2; % X4=dQ/df X2=dp/df elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*

5、e(i)-G(i,i)*f(i);% dP/df 精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV

6、 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B

7、(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开

8、始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+

9、1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J

10、(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q

11、=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行

12、规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j

13、1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=

14、-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1

15、=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p

16、,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m

17、=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节

18、点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i

19、)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角

20、元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,

21、i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1

22、 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列

23、P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)

24、=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end en

25、d end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); %

26、dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,

27、j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*

28、e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k

29、1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end e

30、lse %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; %

31、 PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为

32、求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*

33、e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B

34、(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*

35、i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi

36、 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jaco

37、bi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; e

38、lseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素

39、及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/

40、df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1; % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(

41、i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5

42、;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2

43、)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化精选文档3 X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列 Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列 P J(m,q)=X2; end end else %下面是针对PV 节点来求取Jacobi 矩阵的元素 - DP=P(i)-P1;

44、 % PV 节点有功误差DV=V(i)2-V2; % PV 节点电压误差for j1=1:n if j1=isb&j1=i %非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; elseif j1=i&j1=isb %非平衡

45、节点 &对角元X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节点电压误差m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节点有功误差J(m,q)=X2; end end end end end % 以 上 为 求 雅 可 比 矩 阵 的 各 个 元 素 及 扩 展 列 的 功 率 差 或 电 压 差- for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列 P、 Q 或 U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的P、 Q 或U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除K 行 K 列后的非对角元素进行规格化