matlab和powerworld潮流计算

MATLAB文件以计算示例3总线电源系统中的潮流。

请使用Newton-Raphson和Gauss-Seidel方法的文件夹一个。执行Newton_Raphson_task_1a.m文件并获得电源和无功功率

功率流经所有传输线（两端）并由使用这两种方法在1和2母线上的发电机。画出三个总线系统，并添加P和Q功率值和方向，执行Gauss_seidel_task_1a程序。是否有相同的结果？为什么需要这么多迭代？这个程序实现了上述提到的潮流计算方法

​想了解详情请下载附件http://www.apollocode.net/a/442.html

 

% Example 5-4 Power Flow in a 3-bus Test System
clear
j = sqrt(-1);
V = zeros(3,1);
S = zeros(3,1);
Mismatch = zeros(3,1);
Copyright © 2007 by Ned Mohan. 12
% ---------- Input line impedances ------------- %
Z = [0 0.0047 + 0.0474i 0.0062 + 0.0632i
0.0047 + 0.0474i 0 0.0047 + 0.0474i
0.0062 + 0.0632i 0.0047 + 0.0474i 0];
%---------Base Values ---------%
kVLL=345;
MVA3Ph=100;
Zbase=kVLL^2/MVA3Ph;
XL_km=0.376; % ohm/km at 60 Hz
RL_km= 0.037; B_km=4.5; % B in micro-mho/km
%---------Line Susceptances--------%
B13_Micro_Mho=4.5*200; %200 km long
B12_Micro_Mho=4.5*150; %150 km long
B23_Micro_Mho=4.5*150; %150 km long
%---------Line impedances------------%
Z13_ohm=(RL_km+j*XL_km)*200; %200 km long
Z12_ohm=(RL_km+j*XL_km)*150; %150 km long
Z23_ohm=(RL_km+j*XL_km)*150; %150 km long
%------- line impedances in per unit--------%
Z13=Z13_ohm/Zbase;
Z12=Z12_ohm/Zbase;
Z23=Z23_ohm/Zbase;
%-------- susceptances in per unit----------%
B13=B13_Micro_Mho*Zbase*10^-6;
B12=B12_Micro_Mho*Zbase*10^-6;
B23=B23_Micro_Mho*Zbase*10^-6;
%---------- YBUS Creation-------------%
Y(1,1)=1/Z12 + 1/Z13;
Y(1,2)=-1/Z12;
Y(1,3)=-1/Z13;
Y(2,1)=-1/Z12;
Y(2,2)=1/Z12 + 1/Z23;
Y(2,3)=-1/Z23;
Y(3,1)=-1/Z13;
Y(3,2)=-1/Z23;
Y(3,3)=1/Z13 + 1/Z23;
Y % Print Y=G+jB Admittance Matrix