Hi every one
I need help for solving DC power flow using WLS function. I have tried
following and need guidance for developing it further.
****************************************************************************************************
cd = loadcase('case4gs');
casesystem = size(cd.bus, 1);
cd.bus(:,6) = zeros(length(cd.bus(:,6)), 1)
cd.branch((cd.branch(:,9) ~= 0),9) = zeros(size(cd.branch((cd.branch(:,9)
~= 0),9)));
pfres = runpf(cd);
%% Importing line Data
linedata = cd.branch(:, 1:6); % Calling "linedata6.m" for Line Data...
linedata(:,6) = ones(size(linedata, 1),1);
fb = linedata(:,1); % From bus number...
tb = linedata(:,2); % To bus number...
r = linedata(:,3); % Resistance, R...
x = linedata(:,4); % Reactance, X...
b = linedata(:,5)./2; % Ground Admittance, B/2...
a = linedata(:,6); % Tap setting value.. %Assumed 1
y = 1./(r + 1j*x);
nbus = max(max(fb),max(tb)); % no. of buses...
nbranch = length(fb); % no. of branches...ยจ
%% Formation of Y-Bus(Admittance) Matrix
ybus = zeros(nbus,nbus); % Initialise YBus...
% Formation of the Off-Diagonal Elements...
for k=1:nbranch
ybus(fb(k),tb(k)) = ybus(fb(k),tb(k))-y(k)/a(k);
ybus(tb(k),fb(k)) = ybus(fb(k),tb(k));
end
% Formation of Diagonal Elements...
for m =1:nbus
for n =1:nbranch
if fb(n) == m
ybus(m,m) = ybus(m,m) + y(n)/(a(n)^2) + 1j*b(n);
elseif tb(n) == m
ybus(m,m) = ybus(m,m) + y(n) + 1j*b(n);
end
end
end
ybus
%% Formation of B-Bus
Bp=-imag(ybus-diag(sum(ybus)));
Bpp=imag(ybus)
****************************************************************************************************
Regards
Saeed Ahmed
