电力系统分析_夏道止_第四章_例4-3_潮流计算_MATLAB源码

%电压
clc
clear
theta=[0 0 0 0]
U=[1 1 1.05 1.05]
%功率
G1=0+0j;
G2=-0.5-0.3j;
G3=0.2;
G4=0-0j
%自导纳
Y=zeros(4);  %设置导纳矩阵初值为0
Y(1,1)=2.9049 -11.5040i;
Y(2,2)=0 - 4.6647i;
Y(3,3)=2.4900 - 4.7040i;
Y(4,4)=2.0750 - 3.9092i;

%互导纳
Y(1,2)=0 + 5.3191i;
Y(2,1)=Y(1,2);
Y(1,3)=-1.6600 + 3.1619i;
Y(3,1)=Y(1,3);
Y(1,4)=-1.2450 + 2.3714i;
Y(4,1)=Y(1,4);
Y(3,4)=-0.8300 + 1.5809i;
Y(4,3)=Y(3,4);

for tatal_loop=1:4

%功率平衡