Lab Report
Lab Report
Lab Report
Aim: Bus admittance matrix (Ybus) formation for power systems using inspection method.
Theory:
Bus admittance matrix or Ybus is matrix which gives the information about the admittances of
lines connected to the node as well as the admittance between the nodes. Principal diagonal
elements are called self-admittances of node and is equal to the algebraic sum of all the
admittances terminating at the node. Off diagonal elements are called mutual admittances and are
equal to the admittances between the nodes. The size of ybus is n*n where n is the number of
buses in the system and m= n+1( the total number of buses including the reference buses).
I bus = Y bus * V bus where I bus = vector of impressed bus currents
Inspection method makes use of KVL at all the nodes to get the current equations.
From these equations, Ybus can be directly written. It is the simplest and direct
method of obtaining all the diagonal elements as well as off diagonal elements in the matrix
of any power system. Bus admittance matrix is a sparse matrix. It is often used in solving
load flow problems. Sparsity is one of its greatest advantages as it heavily reduces computer
clc
clear all
n=input(‘enter no. of buses:’);
nl= input(‘enter no. of lines:’);
sb= input(‘enter starting bus of each line:’);
eb= input(‘enter ending bus of each line’);
zser= input(‘enter resistance and reactance of each line:’);
yshty= input(‘enter shunt admittance of the bus:’);
i=0
k=1
while i<nl
zser1(i+1)=zser(k)+j*zser(k+1)
i=i+1
k=k+2
end
zser2=reshape(zser1,nl,1);
yser=ones(nl,1)./zser2;
ybus=zeros(n,n);
for i=1:nl
ybus(sb(i),sb(i))=ybus(sb(i),sb(i))+(j*yshty(i))+yser(i);
ybus(eb(i),eb(i))=ybus(eb(i),eb(i))+(j*yshty(i))+yser(i);
ybus(sb(i),eb(i))=-yser(i);
ybus(eb(i),sb(i))=-yser(i);
end
Input:
n=5
nl=7
sb=[1 1 2 2 2 3 4]
eb=[2 3 3 4 5 4 5]
zser=[0.02 0.06 0.08 0.24 0.06 0.18 0.06 0.18 0.04 0.12 0.01 0.03 0.08 0.24]
yshty=[0.03 0.025 0.02 0.02 0.015 0.01 0.025]
Output:
bus admittance matrix ybus
6.25+(-18.70)j -5.00+( 15.00)j -1.25+( 3.75)j 0.00+( 0.00)j 0.00+( 0.00)j
-5.00+( 15.00)j 10.83+(-32.42)j -1.67+( 5.00)j -1.67+( 5.00)j -2.50+( 7.50)j
-1.25+( 3.75)j -1.67+( 5.00)j 12.92+(-38.70)j-10.00+( 30.00)j 0.00+( 0.00)j
0.00+( 0.00)j -1.67+( 5.00)j-10.00+( 30.00)j 12.92+(-38.70)j -1.25+( 3.75)j
0.00+( 0.00)j -2.50+( 7.50)j 0.00+( 0.00)j -1.25+( 3.75)j 3.75+(-11.21)j
Conclusion:
Ybus generation using inspection method are very useful to calculate on diagonal element and
off diagonal elements of admittance matrix. As we can see from outputs it gets to admittance
matrix and this matrix is sparse matrix. This sparse matrix is used to solve power flow problems
such as, currents of the buses, voltage and power angle of the buses.
Expt.No-1(ii)
Aim: Bus admittance matrix (Y Bus) formation for power systems without mutual coupling
using singular transformation method.
Theory: The Y Matrix is designated by Ybus and called the bus admittance matrix. Y matrix is a
symmetric and square matrix that completely describes the configuration of power transmission
lines. In realistic systems which are quite large containing thousands of buses, the Y matrix is
quite sparse. Each bus in a real power system is usually connected to only a few other buses
through the transmission lines. The Y Matrix is designated by Ybus and called the bus
admittance matrix. Y matrix is a symmetric and square matrix that completely describes the
configuration of power transmission lines. In realistic systems which are quite large containing
thousands of buses, the Y matrix is quite sparse. Each bus in a real power system is usually
connected to only a few other buses through the transmission lines. Ybus can be alternatively
assembled by use of singular transformation given by a graph theoretical approach. This
alternative approach is of great theoretical and practical significance.
clc
clear all
n=input (‘enter no. of buses:’);
nl= input (‘enter no. of lines:’);
sb= input (‘enter starting bus of each line:’);
eb= input (‘enter ending bus of each line’); % ending bus of a line
zser= input (‘enter resistance and reactance of each line:’); % line resistance and reactance (R, X)
yshty= input (‘enter shunt admittance of the bus:’); % shunt admittance
i=0;
k=1;
while i<nl
zser1(i+1)=zser(k)+j*zser(k+1); % impedance of a line (R+jX)
i=i+1;
k=k+2;
end
zser2=reshape(zser1,nl,1);
yser=ones(nl,1)./zser2;
ypri=zeros(nl+n,nl+n);
ybus=zeros(n,n);
a=zeros(nl+n,n);
for i=1:n
a(i,i)=1;
end
for i=1:nl
a(n+i,sb(i))=1;
a(n+i,eb(i))=-1;
ypri(n+i,n+i)=yser(i);
ypri(sb(i),sb(i))=ypri(sb(i),sb(i))+yshty(i);
ypri(eb(i),eb(i))=ypri(sb(i),eb(i))+yshty(i);
end
at=transpose(a);
ybus=at*ypri*a;
zbus=inv(ybus);
Input:
n= 3
nl=3
sb=[1 1 2 ]
eb=[2 3 3 ]
zser=[0.01 0.03 0.08 0.24 0.06 0.18 ]
yshty=[0.01 0.025 0.02]
Output:
bus admittance matrix
ybus:
11.29 + (-33.75)j -10.00 + ( 30.00)j -1.25 + ( 3.75)j
-10.00 + ( 30.00)j 11.70 + (-35.00)j -1.67 + ( 5.00)j
-1.25 + ( 3.75)j -1.67 + ( 5.00)j 2.94 + ( -8.75)j
Conclusion:
In conclusion, the formation of the Ybus matrix using singular transformation without
considering mutual coupling is a useful technique in power system analysis. By neglecting
mutual coupling, the Ybus matrix can be simplified, leading to computational efficiency while
still providing reasonably accurate results under certain conditions.
The Ybus matrix represents the admittance matrix of a power system network and is crucial for
various power system analysis applications, such as load flow studies, fault analysis, and stability
analysis. It characterizes the relationship between voltage and current at different nodes or buses
in the network.
The singular transformation technique involves transforming the Ybus matrix into a diagonal
matrix by diagonalizing the self-admittance elements. This simplifies the matrix and reduces the
number of operations required for subsequent analysis. However, mutual coupling, which
accounts for the interaction between different elements of the network, is ignored in this
approach.
While neglecting mutual coupling may lead to some loss of accuracy in certain scenarios, such as
when dealing with heavily interconnected networks or when analyzing systems with significant
coupling effects, the singular transformation technique can still provide reasonably accurate
results for many practical cases. It is particularly useful when studying systems with weak or
negligible mutual coupling.
It is important to note that the decision to neglect mutual coupling should be made after careful
consideration of the network characteristics and the specific analysis objectives. In cases where
mutual coupling cannot be ignored, alternative methods, such as the full Ybus formation with
mutual coupling considered, should be employed to ensure accurate results.
In summary, the Ybus formation using singular transformation without mutual coupling is a
valuable technique for simplifying the analysis of power system networks. While it may not
capture the full complexity of mutual coupling, it offers computational advantages and can
provide accurate results for systems with weak or negligible coupling effects.
Expt.No-1(iii)
Aim: Bus admittance matrix (Ybus) formation for power systems with mutual coupling using
singular transformation method.
Theory: The current flows through the coil and produces the flux in the same coil and this flux
also links with the neighboring coil. The amount of the flux linking with the second coil is called
mutual coupling. Then we say the two coils are mutually coupled. The amount of the energy
spent in mutual coupling is measured by an impedance called mutual impedance. Hence the
matlab program (Ybus formation using singular transformation without mutual coupling) should
be modified by considering mutual coupling.
MATLAB program:
clc
clear all
n=input(‘enter no. of buses:’); % no. of buses excluding reference
nl= input(‘enter no. of lines:’); % no. of transmission lines
sb= input(‘enter starting bus of each line:’); % starting bus of a line
eb= input(‘enter ending bus of each line’); % ending bus of a line
zser= input(‘enter resistance and reactance of each line:’); % line resistance and reactance (R, X)
nmc= input(‘enter no. of mutual couplings:’);
f1= input(‘enter first line no. of each mutual coupling:’);
s1= input(‘enter second line no. of each mutual coupling:’);
mz= input(‘enter mutual impedance between lines:’);
i=0
k=1
while i<nl
zser1(i+1)=zser(k)+j*zser(k+1) % impedance of a line (R+jX)
i=i+1
k=k+2
end
i=0
k=1
while i<nmc
mz1(i+1)=mz(k)+j*mz(k+1)
i=i+1
k=k+2
end
zser2=reshape(zser1,nl,1);
mz2=reshape(mz1,2,1);
zpri=zeros(nl,nl);
a=zeros(nl,n)
for i=1:nl
zpri(i,i)=zser2(i);
if sb(i)~=0
a(i,sb(i))=1;
end
if eb(i)~=0
a(i,eb(i))=-1
end
end
for i=1:nmc
zpri(f1(i),s1(i))=mz2(i);
zpri(s1(i),f1(i))=mz2(i);
end
ypri=inv(zpri)
at=transpose(a)
ybus=at*ypri*a
zbus=inv(ybus)
Input:
n=3
nl=5
sb=[0 0 2 0 1]
eb=[1 2 3 1 3]
zser=[0 0.6 0 0.5 0 0.5 0 0.4 0 0.2]
nmc=2
f1=[1 1]
s1=[2 4]
mz=[0 0.1 0 0.2 ]
Output:
Bus admittance matrix is
8.021 -0.208 -5.000
-0.208 4.083 -2.000
-5.000 -2.000 7.000
Bus impedance matrix is
0.271 0.126 0.230
0.126 0.344 0.189
0.230 0.189 0.361
Conclusion:
In conclusion, the Ybus formation using singular transformation with mutual coupling is a
comprehensive and accurate technique for analyzing power system networks. By incorporating
the effects of mutual coupling between network elements, this approach provides a more realistic
representation of the system behavior and enables a more detailed analysis of its performance.
The Ybus matrix, which represents the admittance matrix of the network, plays a crucial role in
various power system analysis tasks. It captures the complex interconnections and interactions
between different nodes or buses in the system, accounting for both self-admittance and mutual
admittance.
The singular transformation technique, when applied with consideration for mutual coupling,
involves diagonalizing the Ybus matrix to transform it into a diagonal matrix. This process takes
into account the mutual coupling effects, allowing for a more accurate representation of the
network behavior.
By considering mutual coupling, the Ybus matrix reflects the coupling between different
elements of the network, such as transmission lines, transformers, and generators. This enables a
more comprehensive analysis of power flow, fault analysis, stability studies, and other system-
wide assessments. It provides a better understanding of the system's dynamic behavior, voltage
distribution, and the impact of disturbances or faults.
While the Ybus formation using singular transformation with mutual coupling entails increased
computational complexity compared to neglecting mutual coupling, it offers a higher level of
accuracy and fidelity to real-world conditions. This is particularly important when studying
interconnected systems or analyzing scenarios where mutual coupling effects cannot be ignored.
In conclusion, the Ybus formation using singular transformation with mutual coupling is a robust
and reliable technique for power system analysis. By considering the mutual coupling effects, it
provides a more comprehensive and accurate representation of the network behavior, enabling
engineers and researchers to make informed decisions and ensure the reliable operation of power
systems.
Expt.No-2
Title: Determination of bus currents, bus power and line flows for a specified Bus system
profile.
Aim: To determine the bus currents, bus power and line flows for any power system
Theory: Calculating the power flows, bus currents, and bus power on the various network lines
is the final stage in the load flow analysis process.Examine the bus route that links busses I and
K. As seen here, a circuit with two shunt admittances, Y iko and Y kio, and a series admittance,
Y ik, can be used to represent the line and transformer at each end.
The current fed by bus I into the line can be expressed as,
I ik = I ik1 +I iko = (V i – V k )Y ik + V i Y ik0
The power fed into the line from bus i is
S ik = P ik +jQ ik = V i I ik * = V i (V i * - V k *) Y ik *+ V i V i *Y iko *
The power fed into the line from bus k is
S ki = P ki +jQ ki = V K I ki * = V k (V k * - Vi*) Y ik * + V k V k *Y kio *
The power loss in the (i-k) the line is the sum of the power flows determined from the last two
equations. The transmission loss can be computed by summing all line flows (i.e S ik + S ki for
all i, k).
The slack bus power can also be found by summing the flows on the lines terminating at the
slack bus.
MATLAB program
Clc
clear all
vbus=[1.05+0j; .98-0.06j; 1-0.05j];
yline= [0 0 10 -20 10 -30 10 -20 0 0 16 -32 10 -30 16 -32 0 0 ];
n= 3
i=0;
k=1;
while i<9
yline1(i+1)=yline(k)+j*yline(k+1);
i=i+1
k=k+2
end
yline2=reshape(yline1,3,3) % line admittance
Sp=zeros(n,n)
for i=1:n
for k=1:n
if i==k
continue
else
I(i,k)=yline2(i,k)*(vbus(i)-vbus(k));
I(k,i)= -I(i,k)
S(i,k)=vbus(i)*conj(I(i,k));
S(k,i)=vbus(k)*conj(I(k,i));
Sl(i,k)=S(i,k)+S(k,i);
Sp(i,i)=Sp(i,i)+S(i,k);
end
end
end
Sbus=[Sp(1,1); Sp(2,2); Sp(3,3)]
%bus curent calculation ;-
Ibus=conj((Sbus./vbus))
Input:
Yline= [0 0 10 -20 10 -30 10 -20 0 0 16 -32 10 -30 16 -32 0 0]
n= 3
Output:
line current matrix is
0.00+( 0.00)j 1.90+( -0.80)j 2.00+( -1.00)j
-1.90+( 0.80)j 0.00+( 0.00)j -0.64+( 0.48)j
-2.00+( 1.00)j 0.64+( -0.48)j 0.00+( 0.00)j
line flow matrix is
0.00+( 0.00)j 2.00+( 0.84)j 2.10+( 1.05)j
-1.91+( -0.67)j 0.00+( 0.00)j -0.66+( -0.43)j
-2.05+( -0.90)j 0.66+( 0.45)j 0.00+( 0.00)j
Line losses matrix is
0.00+( 0.00)j 0.09+( 0.17)j 0.05+( 0.15)j
0.09+( 0.17)j 0.00+( 0.00)j 0.01+( 0.02)j
0.05+( 0.15)j 0.01+( 0.02)j 0.00+( 0.00)j
Bus power matrix is
4.10+( 1.89)j 0.00+( 0.00)j 0.00+( 0.00)j
0.00+( 0.00)j -2.57+( -1.10)j 0.00+( 0.00)j
0.00+( 0.00)j 0.00+( 0.00)j -1.39+( -0.45)j
Sbus matrix is
4.10+( 1.89)j -2.57+( -1.10)j -1.39+( -0.45)j
Ibus matrix is
3.90+( -1.80)j -2.54+( 1.28)j -1.36+( 0.52)j
Conclusion:
In conclusion, the determination of bus currents, bus power, and line flows for a specified bus
system profile is crucial for assessing the performance and efficiency of power systems. Through
the utilization of MATLAB, we have successfully computed these parameters for the given
network configuration. The analysis involved calculating line currents, line flows, line losses,
bus power, and slack bus power, providing valuable insights into the distribution and
consumption of power within the system.
By understanding the flow of electricity through various network elements, engineers and
operators can optimize the operation of power systems, identify potential areas for improvement,
and ensure the reliability and stability of electrical grids. This analysis serves as a fundamental
step in load flow studies, aiding in the design, planning, and management of power networks to
meet the demands of modern society effectively.
Moving forward, advancements in computational techniques and modeling methodologies will
continue to enhance our ability to analyze and optimize power systems, contributing to the
development of sustainable and resilient energy infrastructures.
Expt.No-3
Aim: To determine the power angle diagram, reluctance power, excitation emf and regulation of
Theory: The steady state stability is basically concerned with the determination of the maximum
power flow possible through the power system, without loss of synchronism (stability).
The formation of power angle equation plays a vital role in the study of steady state stability. The
power angle equation for non-salient pole machine is given by,
|𝑉||𝐸| ∗ 𝑠𝑖𝑛𝛿
𝑃=
(𝑋𝑠 + 𝑋𝑒)
The above equation shows that the power P transmitted from the generator to the motor varies
with the sine of the displacement angle δ between the two rotors.
For salient pole machine,
MATLAB PROGRAM:
clc
clear all
p=input('power in MW = ');
pf_a=acos(pf);
q=p*tan(pf_a);
i=(p-j*q)*1000000/(3*vt_ph);
delta=0:1:180;
delta_rad=delta*(pi/180);
if xd==xq
ef=vt_ph+(j*i*xd);
excitation_emf=abs(ef)
reg=(abs(ef)-abs(vt_ph))*100/abs(vt_ph)
power_non=abs (ef)*vt_ph*sin(delta_rad)/xd;
net_power=3*power_non/1000000;
plot(delta,net_power);
xlabel('\delta (deg)');
end
if xd~=xq
eq=vt_ph+(j*i*xq);
del=angle(eq);
theta=del+pf_a;
id_mag=abs(i)*sin(theta);
ef_mag=vt_ph*cos(del)+id_mag*xd
reg=(ef_mag-abs(vt_ph))*100/abs(vt_ph)
pp=ef_mag*vt_ph*sin(delta_rad)/xd;
reluct_power=vt_ph^2*(xd-xq)*sin(2*delta_rad)/(2*xd*xq);
net_reluct_power=3*reluct_power/1000000;
power_sal=pp+reluct_power;
net_power_sal=3*power_sal/1000000;
plot(delta,net_reluct_power);
hold on
plot(delta,net_power_sal);
xlabel('\delta (deg)');
end
grid;
power in MW = 48
power factor = .8
xd in ohm = 13.5
xq in ohm = 9.333
Result:
ef_mag = 3.0000e+004
reg = 50.0032
Input data for non salient pole synchronous machine:
power in MW = 48
power factor = .8
xd in ohm = 10
xq in ohm = 10
Result:
excitation_emf = 2.7203e+004
reg = 36.0171
graph-1
The above graph shows the power angle curve. The maximum power transfer occurs at delta=900
Forvalues of delta >900, the power output of the machine reduces successively and finally the
machine may stall.The system is stable only if the displacement angle δ is in the range from -900
to +900 in which the slope dP/d δ is positive.
graph-2
The graph-2 consists of a second harmonic component of power. From the equation
P=V*E/Xd*sin(delta)+ V 2* (Xd-Xq)/2*Xd*Xq*sin(2*delta)
The first term is same as for non-salient pole machine with Xd=Xq. This constitutes the major
part of the power transfer. The second term is quite small (10-20%) compared to the first term
and known as reluctance power.
Conclusion:
We implemented the power angle curves using MATLAB. For non-salient pole machines, we
observed a linear relationship between power and rotor angle. For salient pole machines, the
curve exhibited nonlinearity due to the additional term involving (X_q).
In summary, understanding power angle diagrams is crucial for assessing the steady-state
stability of synchronous machines. These diagrams help us determine the maximum power
transfer capability without losing synchronism. Whether dealing with salient or non-salient pole
machines, these insights are essential for power system analysis and operation.
Expt.No-5
In summary, the experiment on the swing curve for a single synchronous machine connected to
an infinite bus revealed insights into transient stability. We explored the torque angle, swing
equation, and stability criteria. The swing curve, which plots rotor angle against time, helps
engineers design robust power systems that withstand disturbances. While no analytical solution
exists, numerical methods are used to approximate the swing equation’s solution. Understanding
transient stability is crucial for maintaining reliable power grids.
Expt.No-6
𝑃𝑖 = ∑|𝑉𝑖||𝑉𝑗||𝑌𝑖𝑗|cos(𝜃𝑖𝑗 − 𝛿𝑖 + 𝛿𝑗)
𝐽
𝑛
𝑄𝑖 = ∑|𝑉𝑖||𝑉𝑗||𝑌𝑖𝑗|sin(𝜃𝑖𝑗 − 𝛿𝑖 + 𝛿𝑗)
𝐽
Title: MATLAB Program to Solve Load Flow Equations using Gauss-Seidel Method
AIM: To find load flow solution of the given power system using Gauss-Seidel method
theoretically for one iteration and obtain full solution using MATLAB. (only PQ buses)
Theory: Load flow analysis is the study conducted to determine the steady state operating
condition of the given system under given conditions. A large number of numerical algorithms
have been developed and Gauss Seidel method is one of such algorithms.
Problem Formulation
Selecting one of the buses as the reference bus, we get (n-1) simultaneous equations. The bus
𝑃𝑖 = 𝑅𝑒[∑ 𝑉𝑖 ∗ 𝑌𝑖𝑘𝑉𝑘]
𝑘=1
𝑛
𝑄𝑖 = −𝐼𝑚[∑ 𝑉𝑖 ∗ 𝑌𝑖𝑘𝑉𝑘]
𝑘=1
The bus voltage can be written in form of
𝑛
1
𝑉𝑖 = ( ) [𝐼𝑖 − ∑ 𝑌𝑖𝑗𝑉𝑗]
𝑌𝑖𝑖
𝑗=1
j≠i (i=1,2,…………n)& i≠slack bus
Substituting Ii in the expression for Vi, we get
𝑛
1 𝑃𝑖 − 𝐽𝑄𝑖
𝑉𝑖 = ( ) [ − ∑ 𝑌𝑖𝑗𝑉𝑖𝑜
𝑌𝑖𝑖 𝑉𝑖0 ∗
𝑗=1
The latest available voltages are used in the above expression, we get
𝑛 𝑛
1 𝑃𝑖 − 𝐽𝑄𝑖
𝑉𝑖 = ( ) [ − ∑ 𝑌𝑖𝑗𝑉𝑗 𝑛 − ∑ 𝑌𝑖𝑗𝑉𝑖𝑜
𝑌𝑖𝑖 𝑉𝑖0 ∗
𝑗=1 𝑗=𝑖+1
The above equation is the required formula. this equation can be solved for voltages in
interactive manner. During each iteration, we compute all the bus voltage and check for
convergence is carried out by comparison with the voltages obtained at the end of previous
iteration. After the solutions is obtained. The stack bus real and reactive powers, the reactive
power generation at other generator buses and line flows can be calculated.
Algorithm:
Step1: Read the data such as line data, specified power, specified voltages, Q limits at the
generator buses and tolerance for convergences
Step2: Compute Y-bus matrix.
Step3: Initialize all the bus voltages.
Step4: Iter=1
Step5: Consider i=2, where i’ is the bus number.
Step6: check whether this is PV bus or PQ bus. If it is PQ bus goto step 8 otherwise go to next
step.
Step7: Compute Qi check for q limit violation. QGi=Qi+QLi.
If QGi>Qi max ,equate QGi = Qimax. Then convert it into PQ bus.
If QGi<Qi min, equate QGi = Qi min. Then convert it into PQ bus.
Step8: Calculate the new value of the bus voltage using gauss seidal formula. i=1
𝑛 𝑛
1 𝑃𝑖 − 𝐽𝑄𝑖
𝑉𝑖 = ( ) [ − ∑ 𝑌𝑖𝑗𝑉𝑗 𝑛 − ∑ 𝑌𝑖𝑗𝑉𝑖𝑜
𝑌𝑖𝑖 𝑉𝑖0 ∗
𝑗=1 𝑗=𝑖+1
Adjust voltage magnitude of the bus to specify magnitude if Q limits are not violated.
Step9: If all buses are considered go to step 10 otherwise increments the bus no. i=i+1 and Go to
step6.
Step10: Check for convergence. If there is no convergence goes to step 11 otherwise go to
step12.
Step11: Update the bus voltage using the formula.
Vi new =Vi old + α(vi new -Vi old ) (i=1,2,…..n) i≠ slackbus ,α is the acceleration factor=1.4
Step12: Calculate the slack bus power, Q at P-V buses real and reactive give flows real and
reactance line losses and print all the results including all the bus voltages and all the bus angles.
Step13: Stop.
Procedure:
PROBLEM:
a) The following figure shows the one-line diagram of a simple three-bus power system with
generation at bus 1. The magnitude of voltage at bus 1 is adjusted to 1.05 per unit. The
scheduled loads at buses 2 and 3 are as marked on the diagram. Line impedances are marked in
per unit on a 100-MVA base and the line charging susceptances are neglected. Using the Gauss-
Seidel method, determine the phasor values of the voltage at the load buses 2 and 3 (P-Q buses)
accurate to four decimal places and obtain full solution using MATLAB.
MATLAB PROGRAM:
clc
clear all
n=3;
V=[1.05 1 1];
Y=[20-j*50 -10+j*20 -10+j*30
-10+j*20 26-j*52 -16+j*32
-10+j*30 -16+j*32 26-j*62]
P=[inf -2.566 -1.386];
Q=[inf -1.102 -0.452];
diff=10;
iter=1;
Vprev=V;
while (diff>0.00001 | iter==1),
abs(V);
abs(Vprev);
Vprev=V;
for i=2:n
sumyv=0;
for k=1:n,
if(i~=k)
sumyv=sumyv+Y(i,k)*V(k);
end
end
V(i)=(1/Y(i,i))*((P(i)-j*Q(i))/conj(V(i))-sumyv);
end
diff=max(abs(abs(V(2:n))-abs(Vprev(2:n))));
V
iter=iter+1
end
SOLUTION:
At the P-Q buses, the complex loads expressed in per units are
(256.6+𝑗110.2)
✓ 𝑆2𝑠𝑐ℎ = − = −2.566 − 𝑗1.102𝑝𝑢
100
(138.6+𝑗45.2)
✓ 𝑆3𝑠𝑐ℎ = − = −1.386 − 𝑗0.452𝑝𝑢
100
MATLAB OUTPUT:
In this experiment, we focused on solving load flow equations using the Gauss-Seidel method.
The aim was to find the load flow solution for a given power system using the Gauss-Seidel
method. We specifically considered PQ buses (load buses) in the system.
Load flow analysis determines the steady-state operating conditions of a power system.
The performance equation relates bus currents to bus voltages using the admittance matrix.
The Gauss-Seidel method linearizes the nonlinear algebraic power equations. It iteratively
updates bus voltages based on the latest available values. Convergence is checked by comparing
with previous iteration results. The algorithm includes;
✓ Read system data (line data, specified power, voltages, Q limits).
✓ Compute the Y-bus matrix.
✓ Initialize bus voltages.
✓ Iteratively update bus voltages using the Gauss-Seidel method.
✓ Calculate stack bus real and reactive powers, reactive power generation at generator buses,
and line flows.
In summary, the Gauss-Seidel method provides an efficient way to solve load flow equations,
ensuring stable power system operation.