Scribd
Upload
23 views3 pages

Lab 1

This document provides instructions for completing a MATLAB worksheet to analyze a three node circuit using matrix methods. It defines the circuit components and their admittances. It provides 6 steps to: 1) build an admittance matrix Y; 2) calculate currents from voltages; 3) calculate power from voltages and currents; 4) calculate the impedance matrix Z; 5) calculate voltages from currents and Z; 6) calculate power flows and verify they sum to the source power. Lowercase y's denote component admittances while uppercase Y denotes the admittance matrix.

Uploaded by

karlTronxo
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
23 views3 pages

Lab 1

This document provides instructions for completing a MATLAB worksheet to analyze a three node circuit using matrix methods. It defines the circuit components and their admittances. It provides 6 steps to: 1) build an admittance matrix Y; 2) calculate currents from voltages; 3) calculate power from voltages and currents; 4) calculate the impedance matrix Z; 5) calculate voltages from currents and Z; 6) calculate power flows and verify they sum to the source power. Lowercase y's denote component admittances while uppercase Y denotes the admittance matrix.

Uploaded by

karlTronxo
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

Worksheet for Lab1 Power Calculations

In the figure below you are presented with a three node circuit that shall be
analyzed by matrix methods using MATLAB.
I2
I1

I3
2

1
y12

3
y23

y12
y23
y1g
y2g
y3g

=
=
=
=
=

0.1
0.5
0.03
0.05
0.02

+
+
+

0.2i
0.15i
0.05i
0.05i
0.04i

The MATLAB program Lab1.m starts by defining a matrix, Y which has


dimension 3x3 and initialized to a matrix of all zeros. Then statements in
MATLAB are used to define the terms in Y. Note that in the figure above,
the admittances of each individual component are given in lower case y. We
shall use this nomenclature throughout the course: lower case y's and z's will
designate individual component admittances and impedances, capital Y's
and Z's will be used to designate an admittance matrix or an impedance
matrix. When a particular term in a matrix is to be designated we will use a
capital with numbers to indicate row and column. Therefore y12 will
designate the admittance between nodes 1 and 2 as in the figure, whereas
Y12 will designate the row 1 column 2 term in the Y matrix. We shall keep
to this convention inside all MATLAB programs.

y3g

y2g

y1g

All quantities in this circuit are complex (voltages, currents, and


admittances). The admittances for the circuit are:

Step 1) Add the lines of code in Lab1.m to complete building the Y matrix.
Step 2) Given the voltage vector (note that all vectors will be capital and
underlined):
V a 100 + 0 j
V = 100 + 10 j
b

Vc 100 10 j

Step 3 continued)
Find the current vector by multiplying Y times V, that is:

I = YV
Step4) Given the I and V vectors, we can calculate the complex power
delivered by each current source. That is, find:

P1 + jQ 1 = V 1I *1
I = YV
This can be done in MATLAB by using the statement:
S=V.*conj(I)
The .* tells MATLAB to perform the operations pairwise for each element
and not sum over the products. Conj is MATLAB's notation for complex
conjugate. The result is a vector S each element a complex number
representing the P +jQ for the current sources.
Step 5) Calculate the inverse of the Y matrix which is called the impedance
matrix Z by using the MATLAB statement: Z = inv(Y) then multiply Z
times I (obtained previously) to get a new voltage vector Vc. Vc should be
the same as the original V vector. That is:

Vc = ZI

Step6) Use MATLAB to calculate the P+jQ flowing on the individual


components y1g and y12 using the formulas below:

flow 1g = V 1 ( V1 *y1g ) *

flow12 = V 1 (( V 1 V2 ) *y12 )*
Confirm that the sum of these flows equals the P+jQ entering node 1 from
the I1 current source.

You might also like