潮流计算

高斯-赛德尔法：
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-B 2*f)+eu20*(G2*f+B2*eu);P3-e u30*(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)-U f20*(G2*f1+B2*eu1);
F2=Q2-Uf20*(G2*eu1-B2*f1)+U e20*(G2*f1+B2*eu1);
F3=P3-Ue30*(G3*eu1-B3*f1)-U f30*(G3*f1+B3*eu1);
F4=Q3-Uf30*(G3*eu1-B3*f1)+U e30*(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。