潮流计算

  • 格式:docx
  • 大小:15.33 KB
  • 文档页数:3

高斯-赛德尔法:

Matlab代码:

P2=0;

Q2=0;

P3=-1.0567;

Q3=-0.5447;

U10=1;

U20=1;

U30=1;

y10=0.00145-0.00744i;

y20=0.102i;

y30=0.102i;

y11=0.4368-11.4325i;

y12=-0.4353+11.4251i;

y13=0;

y21=y12;

y22=1.5681-17.0625i;

y23=-1.1328+5.7394i;

y31=y13;

y32=y23;

y33=1.1328-5.6374i;

Y1=[y11 y12 0];

Y2=[y21 0 y23];

Y3=[y31 y32 0];

m=50;

while m>0

A2=(P2-Q2*i)/U20;

A3=(P3-Q3*i)/U30;

U21=(A2/U20-Y2*[U10;U20;U30])/y22;

U20=U21;

U31=(A3/U30-Y3*[U10;U20;U30])/y33;

U30=U31;

m=m-1;

end

S1=U10*(Y1*[U10;U20;U30])';

S12=(U10^2)*y10'+U10*(U10'-U20')*y12';

S21=(U20^2)*y20'+U20*(U20'-U10')*y21';

dS12=S12+S21;

S23=(U20^2)*y20'+U20*(U20'-

U30')*y23';

S32=(U30^2)*y30'+U30*(U30'-U20')*y32';

dS23=S23+S32;

fprintf('高斯法计算结果:\n');

fprintf('U2的结果是:');

disp(U20);

fprintf('U3的结果是:');

disp(U30);

fprintf('节点1的功率为:');

disp(S1);

fprintf('节点12之间的功率损耗为:');

disp(dS12);

fprintf('节点23之间的功率损耗为:');

disp(dS23);

仿真计算结果:

高斯法计算结果:

U2的结果是:1.0214 - 0.0974i

U3的结果是:1.0159 - 0.2868i

节点1的功率为:1.1047 - 0.2796i

节点12之间的功率损耗为: -0.0232

- 0.2116i

节点23之间的功率损耗为: -0.1204

- 0.4084i 牛拉法:

Matlab代码:

P2=0;

Q2=0;

P3=-1.0567;

Q3=-0.5447;

Ue10=1;

Uf10=0;

Ue20=0.9;

Uf20=0;

Ue30=0.8;

Uf30=0;

y10=0.00145-0.00744i;

y20=0.102i;

y30=0.102i;

y11=0.4368-11.4325i;

y12=-0.4353+11.4251i;

y13=0;

y21=y12;

y22=1.5681-17.0625i;

y23=-1.1328+5.7394i;

y31=y13;

y32=y23;

y33=1.1328-5.6374i;

G11=0.4368;

B11=-11.4325;

G12=-0.4353;

B12=11.4251;

G13=0;

B13=0;

G21=G12;

B21=B12;

G22=1.5681;

B22=-17.0625;

G23=-1.1328;

B23=5.7394;

G31=G13;

B31=B13;

G32=G23;

B32=B23;

G33=1.1328;

B33=-5.6374;

G1=[G11 G12 G13];

G2=[G21 G22 G23];

G3=[G31 G32 G33];

B1=[B11 B12 B13];

B2=[B21 B22 B23];

B3=[B31 B32 B33];

m=50;

while m>0

syms eu20 f20 eu30 f30

syms f p v

eu=[1;eu20;eu30];

f=[0;f20;f30];

p =

[P2-eu20*(G2*eu-B2*f)-f20*(G2*f+B2*eu);Q2-f20*(G2*eu-B2*f)+eu20*(G2*f+B2*eu);P3-eu30*(G3*eu-B3*f)-f30*(G3*f+B3*eu);Q3-f30*(G3*eu-B3*f)+eu30*(G3*f+B3*eu)];

v = [eu20, f20, eu30,

f30];

R = jacobian(p,v);

J=subs(R,{eu20,f20,eu30,f30},[Ue20,Uf20,Ue30,Uf30]);

disp(J);

eu1=[1;Ue20;Ue30];

f1=[0;Uf20;Uf30];

F1=P2-Ue20*(G2*eu1-B2*f1)-Uf20*(G2*f1+B2*eu1);

F2=Q2-Uf20*(G2*eu1-B2*f1)+Ue20*(G2*f1+B2*eu1);

F3=P3-Ue30*(G3*eu1-B3*f1)-Uf30*(G3*f1+B3*eu1);

F4=Q3-Uf30*(G3*eu1-B3*f1)+Ue30*(G3*f1+B3*eu1);

J0=J\[F1;F2;F3;F4];

H=[Ue20;Uf20;Ue30;Uf30]-J0;

Ue20=H(1,1);

Uf20=H(2,1);

Ue30=H(3,1);

Uf30=H(4,1);

m=m-1;

end

U10=1;

U20=Ue20+i*Uf20;

U30=Ue30+i*Uf30;

S1=U10*([y11 y12

y13]*[U10;U20;U30])';

S12=(U10^2)*y10'+U10*(U10'-U20')*y12';

S21=(U20^2)*y20'+U20*(U20'-U10')*y21';

dS12=S12+S21;

S23=(U20^2)*y20'+U20*(U20'-U30')*y23';

S32=(U30^2)*y30'+U30*(U30'-U20')*y32';

dS23=S23+S32;

fprintf('牛拉法计算结果:\n');

fprintf('U2的结果是:');

disp(U20);

fprintf('U3的结果是:');

disp(U30);

fprintf('节点1的功率为:');

disp(S1);

fprintf('节点12之间的功率损耗为:');

disp(dS12);

fprintf('节点23之间的功率损耗为:');

disp(dS23); 仿真计算结果:

牛拉法计算结果:

U2的结果是:0.9026 - 0.0972i

U3的结果是:0.6697 - 0.2504i节点1的功率为:1.1544 + 1.0781i

节点12之间的功率损耗为: -0.0247

- 0.2911i

节点23之间的功率损耗为: -0.1401

- 0.5673i