Modern Power System Analysis with MATLAB® Applications by Pearson

Modern Power System Analysis with MATLAB Applications

196

if Bustype(ii)==2

qvec(ii) = 0;

end

end

%

% Formation of Bprime and Bdouble prime matrices

%

Bprime=-Bmat;

for i=1:Nbus

if Bustype(i)==1

ijk=i;

for ii=1:Nbus

Bprime(i,ii)=0;

Bprime(ii,i)=0;

end

end

end

Bprime(ijk,ijk)=1;

Bprime

Bdouble=Bprime;

for i=1:Nbus

if qvec(i)==0

ijk=i;

for ii=1:Nbus

Bdouble(ijk,ii)=0;

Bdouble(ii,ijk)=0;

end

Bdouble(ijk,ijk)=1;

end

end

Bdouble

%

% General settings

%

flag = 0;

iter = 1;

itmax = 4;

maxerror=0;

tol = 1e-6;

itmax

tol

for Count=1:4

M04 Modern Power System Analysis 01 XXXXX.indd 196 1/6/2020 12:07:28 PM

Power Flow Analysis

197

for i=1:Nbus

Pcal(i)=0;

Qcal(i)=0;

end

Count

for i=1:Nbus ...