Scribd
Upload
14 views

Loadflowfinal Kad Matlab Code

Uploaded by

sdmrnmist
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
14 views

Loadflowfinal Kad Matlab Code

Uploaded by

sdmrnmist
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

6/24/23, 7:51 AM loadflowfinal_kad

Contents
Formation of Y matrix from the given data

Solution using Gauss-Seidel Method


Calculation of voltage and phase angle of slack, generator and load buses
Calculation of real and reactive power of Slack Bus

%{
EECE 305 Open Ended Lab
Group Members:
202116029 : K. M. Taskin Zaman Kafi
202116034 : Mahmuda Akter
202116029 : MD Abrar Shahriar Kabir
%}

clc
clear all
close all

% Given Bus Data


bus_data = [ 1 1.06 0.0 0 0 0 0
2 1.0 0.0 20 20 40 30
3 1.0 0.0 -45 -15 0 0
4 1.0 0.0 -40 -5 0 0
5 1.0 0.0 -60 -10 0 0 ];

% Given Transmission Line Data


tldata = [ 1 2 0.02 0.06 0.03
1 3 0.08 0.24 0.025
2 3 0.06 0.25 0.02
2 4 0.06 0.18 0.02
2 5 0.04 0.12 0.015
3 4 0.01 0.03 0.01
4 5 0.08 0.24 0.025];

Formation of Y matrix from the given data

b1 = tldata(:,1); %starting bus


b2 = tldata(:,2); %ending bus
R = tldata(:,3); %resistance between lines
X = tldata(:,4); %impedance between lines
j = sqrt(-1);
B = j*tldata(:,5); %susceptance
n = length(b1);
bus = max(max(b1),max(b2)); %number of bus
Z = R + j*X; %Impedance Matrix
y = ones(n,1)./Z ; %admittance

Ymat = zeros(bus,bus);
%Finding OFF DIAGONAl elements
for m = 1:n
Ymat(b1(m),b2(m)) = Ymat(b1(m),b2(m)) - y(m);
Ymat(b2(m),b1(m)) = Ymat(b1(m),b2(m));
end
file:///G:/Anonno/MIST Documents/L3, T1/EECE 306/loadflowfinal_kad.html 1/5
6/24/23, 7:51 AM loadflowfinal_kad
%Finding the diagonal elements
for m=1:bus
for k=1:n
if b1(k)== m
Ymat(m,m) = Ymat(m,m)+y(k)+ B(k);
elseif b2(k)== m
Ymat(m,m)= Ymat(m,m)+y(k)+ B(k);
else
end
end
end
Ymat

Ymat =

Columns 1 through 4

6.2500 -18.6950i -5.0000 +15.0000i -1.2500 + 3.7500i 0.0000 + 0.0000i


-5.0000 +15.0000i 10.0744 -31.1971i -0.9077 + 3.7821i -1.6667 + 5.0000i
-1.2500 + 3.7500i -0.9077 + 3.7821i 12.1577 -37.4771i -10.0000 +30.0000i
0.0000 + 0.0000i -1.6667 + 5.0000i -10.0000 +30.0000i 12.9167 -38.6950i
0.0000 + 0.0000i -2.5000 + 7.5000i 0.0000 + 0.0000i -1.2500 + 3.7500i

Column 5

0.0000 + 0.0000i
-2.5000 + 7.5000i
0.0000 + 0.0000i
-1.2500 + 3.7500i
3.7500 -11.2100i

Solution using Gauss-Seidel Method


Given Base = 100MVA

base=100;
VB = bus_data(:,2); %Bus Voltage in p.u.
PB = bus_data(:,3); %Bus Phase
P = bus_data(:,4); %Real Power of Bus
Q = bus_data(:,5); %Reactive Power of Bus
bus_n = length(bus_data(:,1));
VR = zeros(bus_n,1); %Resultant Voltage
VR = VB;

V_NEW= VR;
for jk= 1:200
for i=2:bus_n
SUM=0;
for j=1:bus_n
if i~=j
SUM= SUM+Ymat(i,j)*V_NEW(j);
end
end
VR(i) = (((((P(i)-1j*Q(i)))/base)/conj(V_NEW(i)))-SUM)/Ymat(i,i);

file:///G:/Anonno/MIST Documents/L3, T1/EECE 306/loadflowfinal_kad.html 2/5


6/24/23, 7:51 AM loadflowfinal_kad
V_NEW(i) = VR(i);
end
end
V_NEW
VR = V_NEW

imaginary_part = zeros(bus_n,1);
for ph=1:bus_n
imaginary_part(ph,1)= imag(VR(ph,1));
end
real_part = zeros(bus_n,1);
for re = 1:bus_n
real_part(re,1) = real(VR(re,1));
end

V_NEW =

1.0600 + 0.0000i
1.0464 - 0.0503i
1.0187 - 0.0929i
1.0179 - 0.0978i
1.0118 - 0.1093i

VR =

1.0600 + 0.0000i
1.0464 - 0.0503i
1.0187 - 0.0929i
1.0179 - 0.0978i
1.0118 - 0.1093i

Calculation of voltage and phase angle of slack, generator and load buses

theta= zeros(bus_n,1);
for th =1:bus_n

theta(th,1) = atand(imaginary_part(th,1)/real_part(th,1));

end
fprintf('The column matrix containg Voltage Magnitude:')
magg = abs(VR);
magg
fprintf('The column matrix containg Phase:')
theta

for print=1:bus_n
fprintf('The voltage of BUS %d is %f \n ',print,magg(print,1))
end
fprintf('\n \n \n')
for print=1:bus_n
fprintf('The phase of BUS %d is %f \n ',print,theta(print,1))
end
fprintf('\n \n \n')

file:///G:/Anonno/MIST Documents/L3, T1/EECE 306/loadflowfinal_kad.html 3/5


6/24/23, 7:51 AM loadflowfinal_kad

The column matrix containg Voltage Magnitude:


magg =

1.0600
1.0476
1.0229
1.0225
1.0176

The column matrix containg Phase:


theta =

0
-2.7529
-5.2086
-5.4879
-6.1672

The voltage of BUS 1 is 1.060000


The voltage of BUS 2 is 1.047578
The voltage of BUS 3 is 1.022892
The voltage of BUS 4 is 1.022541
The voltage of BUS 5 is 1.017641

The phase of BUS 1 is 0.000000


The phase of BUS 2 is -2.752888
The phase of BUS 3 is -5.208649
The phase of BUS 4 is -5.487914
The phase of BUS 5 is -6.167155

Calculation of real and reactive power of Slack Bus

S=0;
slack=1;
SUM=0;
for j=2:bus_n
SUM=SUM+Ymat(slack,j)*VR(j,1);
end
S=VR(slack,1)*(VR(slack,1)*Ymat(slack,slack)+SUM);
S;
P_slack= real(S)*100;
fprintf('Real power of slack bus is P= %d \n' ,P_slack)
Q_slack = imag(S)*100;
fprintf('Reactive power of slack bus is Q= %d \n' ,Q_slack)

Real power of slack bus is P= 1.296117e+02


Reactive power of slack bus is Q= 7.047667e+00

file:///G:/Anonno/MIST Documents/L3, T1/EECE 306/loadflowfinal_kad.html 4/5

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy