Circuit Analysis I

with MATLAB® Computing and

Simulink®/SimPowerSystems® Modeling
Steven T. Karris

Orchard Publications
www.orchardpublications.com
 Students and working professionals will find Circuit
Circuit Analysis I Analysis I with MATLAB® Computing and
with MATLAB® Computing and Simulink®/SimPowerSystems Modeling to be a con-
cise and easy-to-learn text. It provides complete,
Simulink®/SimPowerSystems Modeling clear, and detailed explanations of the principal elec-
trical engineering concepts, and these are illustrated
with numerous practical examples.

This text includes the following chapters and appendices:

• Basic Concepts and Definitions • Analysis of Simple Circuits • Nodal and Mesh Equations -
Circuit Theorems • Introduction to Operational Amplifiers • Inductance and Capacitance
• Sinusoidal Circuit Analysis • Phasor Circuit Analysis • Average and RMS Values, Complex Power,
and Instruments • Natural Response • Forced and Total Response in RL and RC Circuits •
Introduction to MATLAB® • Introduction to Simulink® • Introduction to SimPowerSystems®
• Review of Complex Numbers • Matrices and Determinants
Each chapter and appendix contains numerous practical applications supplemented with detailed
instructions for using MATLAB, Simulink, and SimPowerSystems to obtain quick and accurate
results.

Steven T. Karris is the founder and president of Orchard Publications, has undergraduate and
graduate degrees in electrical engineering, and is a registered professional engineer in California
and Florida. He has more than 35 years of professional engineering experience and more than 30
years of teaching experience as an adjunct professor, most recently at UC Berkeley, California. His
area of interest is in The MathWorks, Inc. products and the publication of MATLAB® and
Simulink® based texts.
™

Orchard Publications
Visit us on the Internet
www.orchardpublications.com
or email us: [email protected]

ISBN-13: 978-1-934404-17-1

ISBN-10: 1-934404-17-9

$79.95 U.S.A.
Circuit Analysis I with MATLAB® Computing and Simulink® / SimPowerSystems® Modeling

Copyright  2009 Orchard Publications. All rights reserved. Printed in USA. No part of this publication may be
reproduced or distributed in any form or by any means, or stored in a data base or retrieval system, without the prior
written permission of the publisher.

Direct all inquiries to Orchard Publications, 39510 Paseo Padre Parkway, Fremont, California 94538, U.S.A.
URL: http://www.orchardpublications.com

Product and corporate names are trademarks or registered trademarks of the MathWorks, Inc., and Microsoft
Corporation. They are used only for identification and explanation, without intent to infringe.

Library of Congress Cataloging-in-Publication Data

Library of Congress Control Number: 2009923770

ISBN10: 1934404179
ISBN13: 9781934404171

TX 5737590

Disclaimer
The author has made every effort to make this text as complete and accurate as possible, but no warranty is implied.
The author and publisher shall have neither liability nor responsibility to any person or entity with respect to any loss
or damages arising from the information contained in this text.

This book was created electronically using Adobe Framemaker.

 Preface

This text is an introduction to the basic principles of electrical engineering. It is the outgrowth of
lecture notes prepared by this author while employed by the electrical engineering and computer
engineering departments as adjunct instructor at various colleges and universities. Many of the
examples and problems are based on the author’s industrial experience. The text is an expansion
of our previous publication, Circuit Analysis I with MATLAB® Applications, ISBN 9780
970951120, and this text, in addition to MATLAB scripts for problem solution, includes
several Simulink® and SimPowerSystems® models. The pages where these models appear are
indicated n the Table of Contents.
The book is intended for students of college grade, both community colleges and universities. It
presumes knowledge of first year differential and integral calculus and physics. While some
knowledge of differential equations would be helpful, it is not absolutely necessary. Chapters 9 and
10 include stepbystep procedures for the solutions of simple differential equations used in the
derivation of the natural and forces responses. Appendices D and E provide a thorough review of
complex numbers and matrices respectively.
In addition to several problems provided at the end of each chapter, this text includes multiple-
choice questions to test and enhance the reader’s knowledge of this subject. Moreover, the
answers to these questions and detailed solutions of all problems are provided at the end of each
chapter. The rationale is to encourage the reader to solve all problems and check his effort for
correct solutions and appropriate steps in obtaining the correct solution. And since this text was
written to serve as a selfstudy, primary, or supplementary textbook, it provides the reader with a
resource to test the reader’s knowledge.
A previous knowledge of MATLAB® would be very helpful. However he material of this text can
be learned without MATLAB, Simulink and SimPowerSystems. This author highly recommends
that the reader studies this material in conjunction with the inexpensive Student Versions of The
MathWorks™ Inc., the developers of these outstanding products, available from:
The MathWorks, Inc.
3 Apple Hill Drive
Natick, MA, 01760
Phone: 508-647-7000,
www.mathworks.com
[email protected].
Appendix A of this text provides a practical introduction to MATLAB, Appendix B is an
introduction to Simulink, and Appendix C is an introduction to SimPowerSystems. These
packages will be invaluable in later studies such as the design of analog and digital filters.
Preface

Like any other new book, this text may contain some grammar and typographical errors;
accordingly, all feedback for errors, advice and comments will be most welcomed and greatly
appreciated.
Orchard Publications
39510 Paseo Padre Parkway
Suite 315
Fremont, California 94538
www.orchardpublications.com
[email protected]
Table of Contents
1 Basic Concepts and Definitions 11
1.1 The Coulomb ........................................................................................................11
1.2 Electric Current and Ampere ...............................................................................11
1.3 Two Terminal Devices .........................................................................................14
1.4 Voltage (Potential Difference) .............................................................................15
1.5 Power and Energy .................................................................................................18
1.6 Active and Passive Devices ................................................................................111
1.7 Circuits and Networks ........................................................................................113
1.8 Active and Passive Networks .............................................................................113
1.9 Necessary Conditions for Current Flow .............................................................113
1.10 International System of Units ............................................................................114
1.11 Sources of Energy ...............................................................................................117
1.12 Summary .............................................................................................................119
1.13 Exercises .............................................................................................................121
1.14 Answers / Solutions to EndofChapter Exercises ............................................124
MATLAB Computing: Pages 16 through 18

2 Analysis of Simple Circuits 21

2.1 Conventions .........................................................................................................21
2.2 Ohm’s Law ............................................................................................................21
2.3 Power Absorbed by a Resistor ..............................................................................23
2.4 Energy Dissipated in a Resistor ............................................................................24
2.5 Nodes, Branches, Loops and Meshes ...................................................................25
2.6 Kirchhoff’s Current Law (KCL) ...........................................................................26
2.7 Kirchhoff’s Voltage Law (KVL) ............................................................................27
2.8 Single Mesh Circuit Analysis .............................................................................210
2.9 Single NodePair Circuit Analysis ....................................................................214
2.10 Voltage and Current Source Combinations .......................................................217
2.11 Resistance and Conductance Combinations .....................................................218
2.12 Voltage Division Expressions .............................................................................222
2.13 Current Division Expressions .............................................................................224
2.14 Standards for Electrical and Electronic Devices ................................................226
2.15 Resistor Color Code ...........................................................................................226
2.16 Power Rating of Resistors ...................................................................................228
2.17 Temperature Coefficient of Resistance ..............................................................228
2.18 Ampere Capacity of Wires .................................................................................229

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling i

Copyright © Orchard Publications
 2.19 Current Ratings for Electronic Equipment......................................................... 229
2.20 Copper Conductor Sizes for Interior Wiring ......................................................231
2.21 Summary .............................................................................................................236
2.22 Exercises ..............................................................................................................239
2.23 Answers / Solutions to EndofChapter Exercises ............................................247
Simulink / SimPowerSystems models: Pages 224, 226

3 Nodal and Mesh Equations - Circuit Theorems 31

3.1 Nodal, Mesh, and Loop Equations ....................................................................... 31
3.2 Analysis with Nodal Equations ............................................................................ 31
3.3 Analysis with Mesh or Loop Equations ................................................................ 38
3.4 Transformation between Voltage and Current Sources .................................... 320
3.5 Thevenin’s Theorem .......................................................................................... 323
3.6 Norton’s Theorem .............................................................................................. 333
3.7 Maximum Power Transfer Theorem .................................................................. 335
3.8 Linearity.............................................................................................................. 337
3.9 Superposition Principle....................................................................................... 338
3.10 Circuits with Non-Linear Devices...................................................................... 342
3.11 Efficiency ............................................................................................................ 344
3.12 Regulation........................................................................................................... 345
3.13 Summary ............................................................................................................. 347
3.14 Exercises ............................................................................................................. 349
3.15 Answers / Solutions to EndofChapter Exercises ............................................ 360
MATLAB Computing: Pages 3-4, 37, 311, 314, 316, 318, 332, 366,
370, 372, 374, 376, 380, 390
Simulink / SimPowerSystems models: Pages 38, 319

4 Introduction to Operational Amplifiers 41

4.1 Signals................................................................................................................... 41
4.2 Amplifiers ............................................................................................................. 41
4.3 Decibels ................................................................................................................ 42
4.4 Bandwidth and Frequency Response.................................................................... 44
4.5 The Operational Amplifier................................................................................... 45
4.6 An Overview of the Op Amp............................................................................... 45
4.7 Active Filters ...................................................................................................... 413
4.8 Analysis of Op Amp Circuits.............................................................................. 416
4.9 Input and Output Resistance ............................................................................. 428
4.10 Summary ............................................................................................................. 432
4.11 Exercises ............................................................................................................. 434

ii Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling

Copyright © Orchard Publications
 4.12 Answers / Solutions to EndofChapter Exercises............................................ 443
MATLAB Computing: Page 4-47
Simulink / SimPowerSystems models: Pages 452

5 Inductance and Capacitance 51

5.1 Energy Storage Devices....................................................................................... 51
5.2 Inductance........................................................................................................... 51
5.3 Power and Energy in an Inductor...................................................................... 511
5.4 Combinations of Inductors in Series and in Parallel......................................... 514
5.5 Capacitance....................................................................................................... 516
5.6 Power and Energy in a Capacitor ...................................................................... 521
5.7 Combinations of Capacitors in Series and in Parallel ....................................... 524
5.8 Nodal and Mesh Equations in General Terms.................................................. 526
5.9 Summary............................................................................................................ 529
5.10 Exercises ............................................................................................................ 531
5.11 Answers / Solutions to EndofChapter Exercises ........................................... 536
MATLAB Computing: Pages 513, 523, 540

6 Sinusoidal Circuit Analysis 61

6.1 Excitation Functions ........................................................................................... 61
6.2 Circuit Response to Sinusoidal Inputs ................................................................ 61
6.3 The Complex Excitation Function ..................................................................... 63
6.4 Phasors in R, L, and C Circuits........................................................................... 68
6.5 Impedance ......................................................................................................... 614
6.6 Admittance ....................................................................................................... 617
6.7 Summary ........................................................................................................... 623
6.8 Exercises ............................................................................................................ 626
6.9 Solutions to EndofChapter Exercises............................................................ 632
MATLAB Computing: Pages 621, 632, 635
Simulink / SimPowerSystems models: Pages 622, 637, 638

7 Phasor Circuit Analysis 71

7.1 Nodal Analysis .................................................................................................... 71
7.2 Mesh Analysis...................................................................................................... 75
7.3 Application of Superposition Principle ............................................................... 76
7.4 Thevenin’s and Norton’s Theorems ................................................................. 710
7.5 Phasor Analysis in Amplifier Circuits ............................................................... 714

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling iii
Copyright © Orchard Publications
 7.6 Phasor Diagrams ................................................................................................ 717
7.7 Electric Filters .................................................................................................... 722
7.8 Basic Analog Filters ........................................................................................... 723
7.9 Active Filter Analysis ........................................................................................ 728
7.10 Summary ............................................................................................................ 731
7.11 Exercises............................................................................................................. 732
7.12 Answers to EndofChapter Exercises.............................................................. 739
MATLAB Computing: Pages 7-4, 76, 78, 712, 713, 715, 717, 721, 730,
744, 745, 746, 748, 750, 751, 755, 756, 758
Simulink models: Pages 79, 710

8 Average and RMS Values, Complex Power, and Instruments 81

8.1 Periodic Time Functions ...................................................................................... 81
8.2 Average Values .................................................................................................... 82
8.3 Effective Values.................................................................................................... 84
8.4 Effective (RMS) Value of Sinusoids..................................................................... 85
8.5 RMS Values of Sinusoids with Different Frequencies ......................................... 87
8.6 Average Power and Power Factor ........................................................................ 89
8.7 Average Power in a Resistive Load .................................................................... 811
8.8 Average Power in Inductive and Capacitive Loads ........................................... 811
8.9 Average Power in NonSinusoidal Waveforms................................................. 814
8.10 Lagging and Leading Power Factors................................................................... 815
8.11 Complex Power  Power Triangle...................................................................... 816
8.12 Power Factor Correction .................................................................................... 818
8.13 Instruments......................................................................................................... 820
8.14 Summary............................................................................................................. 831
8.15 Exercises ............................................................................................................. 834
8.16 Answers to EndofChapter Exercises .............................................................. 840
MATLAB Computing: Page 8-3 continued on Page 84

9 Natural Response 91

9.1 Natural Response of a Series RL circuit ............................................................... 91
9.2 Natural Response of a Series RC Circuit.............................................................. 99
9.3 Summary ............................................................................................................. 917
9.4 Exercises.............................................................................................................. 918
9.5 Answers to EndofChapter Exercises............................................................... 925
MATLAB Computing: Page 927
Simulink / SimPowerSystems models: Page 924

iv Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling

Copyright © Orchard Publications
10 Forced and Total Response in RL and RC Circuits 101
10.1 Unit Step Function .............................................................................................101
10.2 Unit Ramp Function ..........................................................................................10 6
10.3 Delta Function ...................................................................................................10 7
10.4 Forced and Total Response in an RL Circuit ...................................................1014
10.5 Forced and Total Response in an RC Circuit...................................................1021
10.6 Summary............................................................................................................1031
10.7 Exercises ............................................................................................................1033
10.8 Answers to EndofChapter Exercises .............................................................1040
MATLAB Computing: Pages 1018, 1030
Simulink / SimPowerSystems models: Page 1051

A Introduction to MATLAB A1

A.1 Command Window ..............................................................................................A1
A.2 Roots of Polynomials ............................................................................................A3
A.3 Polynomial Construction from Known Roots ......................................................A4
A.4 Evaluation of a Polynomial at Specified Values ..................................................A5
A.5 Rational Polynomials ...........................................................................................A8
A.6 Using MATLAB to Make Plots ..........................................................................A9
A.7 Subplots .............................................................................................................A18
A.8 Multiplication, Division and Exponentiation ...................................................A19
A.9 Script and Function Files ..................................................................................A26
A.10 Display Formats .................................................................................................A31
MATLAB Computations: Entire Appendix A

B Introduction to Simulink B1

B.1 Simulink and its Relation to MATLAB ............................................................... B1
B.2 Simulink Demos ................................................................................................. B20
Simulink Modeling: Entire Appendix B

C Introduction to SimPowerSystems C1

C.1 Simulation of Electric Circuits with SimPowerSystems ...................................... C1
SimPowerSystems Modeling: Entire Appendix C

D A Review of Complex Numbers D1

D.1 Definition of a Complex Number ........................................................................ D1

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling v

Copyright © Orchard Publications
 D.2 Addition and Subtraction of Complex Numbers.................................................D2
D.3 Multiplication of Complex Numbers ...................................................................D3
D.4 Division of Complex Numbers.............................................................................D4
D.5 Exponential and Polar Forms of Complex Numbers ...........................................D4
MATLAB Computing: Pages D6, D7, D8
Simulink Modeling: Page D7

E Matrices and Determinants E1

E.1 Matrix Definition ................................................................................................ E1
E.2 Matrix Operations............................................................................................... E2
E.3 Special Forms of Matrices ................................................................................... E6
E.4 Determinants..................................................................................................... E10
E.5 Minors and Cofactors ........................................................................................ E12
E.6 Cramer’s Rule .................................................................................................... E17
E.7 Gaussian Elimination Method........................................................................... E19
E.8 The Adjoint of a Matrix.................................................................................... E21
E.9 Singular and NonSingular Matrices................................................................ E21
E.10 The Inverse of a Matrix..................................................................................... E22
E.11 Solution of Simultaneous Equations with Matrices .......................................... E24
E.12 Exercises ............................................................................................................ E31
MATLAB Computing: Pages E3, E4, E5, E7, E8, E9, E10, E12,
E15, E16, E18, E22, E25, E6, E29
Simulink Modeling: Page E3

Excel Spreadsheet: Page E27

References R1
Index IN1

vi Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling

Copyright © Orchard Publications
 Chapter 1
Basic Concepts and Definitions

T
his chapter begins with the basic definitions in electric circuit analysis. It introduces the
concepts and conventions used in introductory circuit analysis, the unit and quantities used
in circuit analysis, and includes several practical examples to illustrate these concepts.
Throughout this text, a left justified horizontal bar will denote the beginning of an example, and
a right justified horizontal bar will denote the end of the example. These bars will not be shown
whenever an example begins at the top of a page or at the bottom of a page. Also, when one
example follows immediately after a previous example, the right justified bar will be omitted.

1.1 The Coulomb

Two identically charged (both positive or both negative) particles possess a charge of one coulomb
–7
when being separated by one meter in a vacuum, repel each other with a force of 10 c 2 newton
8
where c = velocity of light  3  10 m  s . The definition of coulomb is illustrated in Figure 1.1.

q 1m q
Vacuum –7
F = 10 c 2 N
q=1 coulomb

Figure 1.1. Definition of the coulomb

The coulomb, abbreviated as C , is the fundamental unit of charge. In terms of this unit, the
– 19 18
charge of an electron is 1.6  10 C and one negative coulomb is equal to 6.24  10 electrons.
Charge, positive or negative, is denoted by the letter q or Q .

1.2 Electric Current and Ampere

Electric current i at a specified point and flowing in a specified direction is defined as the instan-
taneous rate at which net positive charge is moving past this point in that specified direction, that
is,
dq q
i = ------ = lim ------ (1.1)
dt t  0 t

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling 11
Copyright © Orchard Publications
Chapter 1 Basic Concepts and Definitions

Example 1.4
It is assumed a 12volt automotive battery is completely discharged and at some reference time
t = 0 , is connected to a battery charger to trickle charge it for the next 8 hours. It is also assumed
that the charging rate is
 8e – t  3600 A
it = 
0  t  8 hr
0 otherwise

For this 8hour interval compute:

a. the total charge delivered to the battery
b. the maximum power (in watts) absorbed by the battery
c. the total energy (in joules) supplied
d. the average power (in watts) absorbed by the battery
Solution:
The current entering the positive terminal of the battery is the decaying exponential shown in
Figure 1.11 where the time has been converted to seconds.

i(t)
(A)
8 – t § 3600
i = 8e

t (s)
28800
Figure 1.11. Decaying exponential for Example 1.4
Then,
a.

0 0
15000 28800 28800
15000 – t  3600 8 – t  3600
q = i dt = 8e dt = ---------------------- e
t=0 – 1  3600 0

–8
= – 8  3600  e – 1   28800 C or 28.8 kC
b.
i max = 8 A (occurs at t=0)
Therefore,
p max = vi max = 12  8 = 96 w

 0 0
c.
28800 28800 28800
– t  3600 96 – t  3600
W = p dt = vi dt = 12  8e dt = ---------------------- e
– 1  3600 0

5 –8
= 3.456  10  1 – e   345.6 KJ.

110 Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling
Copyright © Orchard Publications
 Chapter 2
Analysis of Simple Circuits

T
his chapter defines constant and instantaneous values, Ohm’s law, and Kirchhoff’s Current
and Voltage laws. Series and parallel circuits are also defined and nodal, mesh, and loop
analyses are introduced. Combinations of voltage and current sources and resistance com-
binations are discussed, and the voltage and current division formulas are derived.

2.1 Conventions
We will use lower case letters such as v , i , and p to denote instantaneous values of voltage, cur-
rent, and power respectively, and we will use subscripts to denote specific voltages, currents,
resistances, etc. For example, v S and i S will be used to denote voltage and current sources
respectively. Notations like v R1 and i R2 will be used to denote the voltage across resistance R 1
and the current through resistance R 2 respectively. Other notations like v A or v 1 will represent
the voltage (potential difference) between point A or point 1 with respect to some arbitrarily cho-
sen reference point taken as “zero” volts or “ground”.
The designations v AB or v 12 will be used to denote the voltage between point A or point 1 with
respect to point B or 2 respectively. We will denote voltages as v  t  and i  t  whenever we wish
to emphasize that these quantities are time dependent. Thus, sinusoidal (AC) voltages and cur-
rents will be denoted as v  t  and i  t  respectively. Phasor quantities, to be introduced in Chapter
6, will be represented with bold capital letters, V for phasor voltage and I for phasor current.

2.2 Ohm’s Law

We recall from Chapter 1 that resistance R is a constant that relates the voltage and the current
as:
v R = Ri R (2.1)

This relation is known as Ohm’s law.

The unit of resistance is the Ohm and its symbol is the Greek capital letter  . One ohm is the
resistance of a conductor such that a constant current of one ampere through it produces a volt-
age of one volt between its ends. Thus,

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems  Modeling 21
Copyright © Orchard Publications
 Resistance and Conductance Combinations

 RK
n
R eq = R 1 + R 2 + R 3 +  + R n =
k=1
(2.36)
For Resistors in Series

Example 2.4
For the circuit of Figure 2.30, find the value of the current i after combining the voltage sources
to a single voltage source and the resistances to a single resistor.

4 64 V 6 
v + R 


+
v S1 + R v S3
1 S2 2
+ +
i  
200 V  R 4 + 80 V
 R3 +
A 10  8
Figure 2.30. Circuit for Example 2.4
Solution:
We add the values of the voltage sources as indicated in Definitions 8 and 9, we add the resis-
tances in accordance with (2.36), and we apply Ohm’s law. Then,

v
i = ------- 200 –  64 + 80  56
- = -------------------------------------- = ------ = 2 A (2.37)
R 28 28

Next, we consider the case where n resistors are connected in parallel as shown in Figure 2.31.
R1

i1
iT iT
A i2 R2 B

in Rn
v AB
Figure 2.31. Addition of resistances in parallel
By KCL,
iT = i1 + i2 +  + in (2.38)

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems  Modeling 219
Copyright © Orchard Publications
Chapter 2 Analysis of Simple Circuits

R2

4 +
R1 iS
R3 5 v 12 
i2 i1 3A
20 

R4
Figure 2.41. Application of current division expressions for the circuit of Example 2.7

and observing the passive sign convention, the voltage v is

87 1044
v = – i 1 R 1 = – ------  12 = – ------------ V
41 41
or
v = – 25.46 V

Resistances in Ohms -0.88 -2.12

Display 2 Display 3
CCS = Current Controlled Source
+v -25.46
-
VM Display 4
VM = Voltage Measurement
+

R1= 29 R2 = 12
CCS
s
-

-3 + -i + -i Continuous
+ -i CM 2 CM 3
Constant powergui
CM 1

CM = Current Measurement 3.00

Display 1

Figure 2.42. Simulink / SimPower Systems model for Example 2.7

2.14 Standards for Electrical and Electronic Devices

Standardization of electronic components such as resistors, capacitors and diodes is carried out by
various technical committees. In the United States, the Electronics Industries Association (EIA)
and the American National Standards Institute (ANSI) have established and published several
standards for electrical and electronic devices to provide interchangeability among similar prod-
ucts made by different manufacturers. Also, the U.S. Department of Defense or its agencies issue

226 Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling
Copyright © Orchard Publications
 Chapter 3
Nodal and Mesh Equations  Circuit Theorems

T
his chapter begins with nodal, loop and mesh equations and how they are applied to the
solution of circuits containing two or more nodepairs and two or more loops or meshes.
Other topics included in this chapter are the voltagetocurrent source transformations
and vice versa, Thevenin’s and Norton’s theorems, the maximum power transfer theorem, linear-
ity, superposition, efficiency, and regulation.

3.1 Nodal, Mesh, and Loop Equations

Network Topology is a branch of network theory concerned with the equations required to com-
pletely describe an electric circuit. In this text, we will only be concerned with the following two
theorems.

Theorem 3.1

Let N = number of nodes in a circuit ; then N – 1 independent nodal equations are required to
completely describe that circuit. These equations are obtained by setting the algebraic sum of the
currents leaving each of the N – 1 nodes equal to zero.
Theorem 3.2
Let L = M = number of loops or meshes , B = number of branches , N = number of nodes
in a circuit; then L = M = B – N + 1 independent loop or mesh equations are required to com-
pletely describe that circuit. These equations are obtained by setting the algebraic sum of the
voltage drops around each of the L = M = B – N + 1 loops or meshes equal to zero.

3.2 Analysis with Nodal Equations

In writing nodal equations, we perform the following steps:
1. For a circuit containing N nodes, we choose one of these as a reference node assumed to be zero
volts or ground.
2. At each nonreference node we assign node voltages v 1 v 2  v n – 1 where each of these volt-
ages is measured with respect to the chosen reference node, i.e., ground.
3. If the circuit does not contain any voltage sources between nodes, we apply KCL and write a
nodal equation for each of the node voltages v 1 v 2  v n – 1 .

Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling 31
Copyright © Orchard Publications
 Analysis with Nodal Equations

Check with Excel®:

The spreadsheet of Figure 3.3 shows the solution of the equations of (3.7). The procedure is dis-
cussed in Appendix A.

Figure 3.3. Spreadsheet for the solution of (3.7)

Example 3.2
For the circuit of Figure 3.4, write nodal equations in matrix form and solve for the unknowns
using matrix theory, Cramer’s rule, or Gauss’s elimination method. Verify your answers with
Excel or MATLAB. Please refer to Appendix A for procedures and examples. Then construct a
table showing the voltages across, the currents through and the power absorbed or delivered by
each device. Verify your answer with a Simulink / SimPowerSystems model.

+
8 10 V
+ 4 6

12 V 18 A 24 A

Figure 3.4. Circuit for Example 3.2

Solution:
We observe that there are 4 nodes and we denote these as as , , , and G (for ground) as
shown in Figure 3.5. We assign voltages v 1 v 2 , and v 3 at nodes , , and  respectively; these
are to be measured with respect to the ground node G . We observe that v 1 is a known voltage,
that is, v 1 = 12 V and thus our first equation is

v 1 = 12 (3.8)

Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling 35
Copyright © Orchard Publications
 Analysis with Mesh or Loop Equations

Now, we can find the voltage drop across the 5 A current source by application of KVL around
Mesh 1 using the following relation:
2  3.3975 + 4   3.3975 + 1.6025  + v 5A + 10   3.3975 – 1.2315  – 12 = 0
This yields
v 5A = – 36.455

We can verify this value by application of KVL around Mesh 2 where beginning with the lower
end of the 6 w resistor and going counterclockwise we obtain
 6 + 8   1.6025 – 36 + 4   3.3975 + 1.6025  – 36.455 + 16   1.6025 + 0.2737  = 0
With these values, we can also compute the power delivered or absorbed by each of the voltage
sources and the current source.

3.40 C VS 2 36
36 V
I1 K3
+ i s
- +
-
CM 1
2 8
+

C VS 1
12 V
4 6
s
-

-1.60
12 I2 i +
+

-
CCS 1
K1 -5 5A CM 2
s
-

K2

10 16
+

C VS 3
12 V1
R1=3 20
s

0.27
i +
I4 -
12 1.23 CM 4
C VS 4 24
K4 I3 24 V Continuous
i +
s K5
+
- - powergui
CM 3 18

Figure 3.19. Simulink / SimPower Systems model for Example 3.6

Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling 319
Copyright © Orchard Publications
 Chapter 4
Introduction to Operational Amplifiers

T
his chapter is an introduction to amplifiers. It discusses amplifier gain in terms of decibels
(dB) and provides an overview of operational amplifiers, their characteristics and applica-
tions. Numerous formulas for the computation of the gain are derived and several practical
examples are provided.

4.1 Signals
A signal is any waveform that serves as a means of communication. It represents a fluctuating
electric quantity, such as voltage, current, electric or magnetic field strength, sound, image, or
any message transmitted or received in telegraphy, telephony, radio, television, or radar. A typical
signal which varies with time is shown in figure 4.1 where f  t  can be any physical quantity such
as voltage, current, temperature, pressure, and so on.

f t

t
Figure 4.1. A signal that changes with time

4.2 Amplifiers
An amplifier is an electronic circuit which increases the magnitude of the input signal. The sym-
bol of a typical amplifier is a triangle as shown in Figure 4.2.

vin
vout

Electronic Amplifier
Figure 4.2. Symbol for electronic amplifier
An electronic (or electric) circuit which produces an output that is smaller than the input is
called an attenuator. A resistive voltage divider is a typical attenuator.

Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling 41
Copyright © Orchard Publications
 An Overview of the Op Amp

120 K

20 K Rf

v in+ R in v out
 + +


Figure 4.12. Circuit for Example 4.4

Solution:
This is the same circuit as that of the previous example except that the input is a sine wave with
unity amplitude and the voltage gain G v is the same as before, that is,

R 120 K
G v = – -------f- = – -------------------- = – 6
R in 20 K
and the output voltage is
v out = G v v in = – 6  sin t = – 6 sin t mV

The voltages v in and v out are plotted as shown in Figure 4.13.

6
v (mv)
v out = – 6 sin t
4
v in = sin t
2

-2

-4

-6
0 2 4 6 8 10 12

Figure 4.13. Input and output waveforms for the circuit of Example 4.4

An op amp is said to be connected in the noninverting mode when an input signal is connected
to the noninverting () input through an external resistor R which serves as a current limiter,
and the inverting () input is grounded through an external resistor R in as shown in Figure 4.14.
In our subsequent discussion, the resistor R will represent the internal resistance of the applied
voltage v in .

Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling 49
Copyright © Orchard Publications
Chapter 4 Introduction to Operational Amplifiers

- = – ------ R 4 + R 2  -----4- + 1
 
v out 1 R
G v = ---------
v in R1 R3
9.

 Rf +
+ R in ------- v in v out
 + R in 
v in

v in = 1 v peak f = 0.2 Hz Phase = 0 deg R in = 1 K R f = 10 K

VM = Voltage Measurement Vin

Vin = 1 volt peak, frequency 0.2 Hz
Vout
+v
+

- Scope
CVS +v Bus
VM 1 - Continuous
Creator
s
-

Vin VM 2
Rin Product powergui

-Rf/Rin CVS=Controlled Voltage Source

Constant Rin=1 K, Rf= 10 K, entered at the MATLAB command prompt

452 Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling
Copyright © Orchard Publications
 Chapter 5
Inductance and Capacitance

his chapter is an introduction to inductance and capacitance, their voltagecurrent rela-

T tionships, power absorbed, and energy stored in inductors and capacitors. Procedures for
analyzing circuits with inductors and capacitors are presented along with several examples.

5.1 Energy Storage Devices

In the first four chapters we considered resistive circuits only, that is, circuits with resistors and
constant voltage and current sources. However, resistance is not the only property that an elec-
tric circuit possesses; in every circuit there are two other properties present and these are the
inductance and the capacitance. We will see through some examples that will be presented later
in this chapter, that inductance and capacitance have an effect on an electric circuit as long as
there are changes in the voltages and currents in the circuit.

The effects of the inductance and capacitance properties can best be stated in simple differential
equations since they involve the changes in voltage or current with time. We will study induc-
tance first.

5.2 Inductance
Inductance is associated with the magnetic field which is always present when there is an electric
current. Thus, when current flows in an electric circuit the conductors (wires) connecting the
devices in the circuit are surrounded by a magnetic field. Figure 5.1 shows a simple loop of wire
and its magnetic field represented by the small loops.

Figure 5.1. Magnetic field around a loop of wire

The direction of the magnetic field (not shown) can be determined by the lefthand rule if con-
ventional current flow is assumed, or by the righthand rule if electron current flow is assumed.
The magnetic field loops are circular in form and are referred to as lines of magnetic flux. The unit
of magnetic flux is the weber (Wb).

Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling 51
Copyright © Orchard Publications
Chapter 5 Inductance and Capacitance

vC  t 
(V)
1.00
0.75
0.50
0.25
0 t (ms)
0.25 1 2 3 4
0.50
0.75
1.00

Figure 5.23. Voltage waveform for 0  t  2 ms of Example 5.5

There is no need to calculate the values of the capacitor voltage v c at t = 3 ms and at t = 4 ms

because the waveform of the current i c starts repeating itself at t = 2 ms , and the initial condi-
tions and the areas are the same as before. Accordingly, the capacitor voltage v c waveform of fig-
ure (b) starts repeating itself also as shown in Figure 5.24.
vC  t 
(V)
1.00
0.75
0.50
0.25
t (ms)
0
0.25 1 2 3 4
0.50
0.75
1.00
Figure 5.24. Voltage waveform for 0  t  4 ms of Example 5.5

Example 5.5 has illustrated the well known fact that the voltage across a capacitor cannot change
instantaneously. Referring to the current and voltage waveforms for this example, we observe
that the current through the capacitor can change instantaneously as shown by the discontinui-
ties at t = 1 2 3 and 4 ms in Figure 5.21. However, the voltage across the capacitor never
changes instantaneously, that is, it displays no discontinuities since its value is explicitly defined
at all instances of time as shown in Figure 5.24.

520 Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling
Copyright © Orchard Publications
Chapter 6
Sinusoidal Circuit Analysis

T
his chapter is an introduction to circuits in which the applied voltage or current are sinu-
soidal. The time and frequency domains are defined and phasor relationships are developed
for resistive, inductive and capacitive circuits. Reactance, susceptance, impedance and
admittance are also defined. It is assumed that the reader is familiar with sinusoids and complex
numbers. If not, it is strongly recommended that Appendix B is reviewed thoroughly before read-
ing this chapter.

6.1 Excitation Functions

The applied voltages and currents in electric circuits are generally referred to as excitations or driv-
ing functions, that is, we say that a circuit is “excited” or “driven” by a constant, or a sinusoidal, or
an exponential function of time. Another term used in circuit analysis is the word response; this
may be the voltage or current in the “load” part of the circuit or any other part of it. Thus the
response may be anything we define it as a response. Generally, the response is the voltage or cur-
rent at the output of a circuit, but we need to specify what the output of a circuit is.
In Chapters 1 through 4 we considered circuits that consisted of excitations (active sources) and
resistors only as the passive devices. We used various methods such as nodal and mesh analyses,
superposition, Thevenin’s and Norton’s theorems to find the desired response such as the voltage
and/or current in any particular branch. The circuit analysis procedure for these circuits is the
same for DC and AC circuits. Thus, if the excitation is a constant voltage or current, the
response will also be some constant value; if the excitation is a sinusoidal voltage or current, the
response will also be sinusoidal with the same frequency but different amplitude and phase.
In Chapter 5 we learned that when the excitation is a constant and steadystate conditions are
reached, an inductor behaves like a short circuit and a capacitor behaves like an open circuit.
However, when the excitation is a timevarying function such as a sinusoid, inductors and capac-
itors behave entirely different as we will see in our subsequent discussion.

6.2 Circuit Response to Sinusoidal Inputs

We can apply the circuit analysis methods which we have learned in previous chapters to circuits
where the voltage or current sources are sinusoidal. To find out how easy (or how difficult) the
procedure becomes, we will consider the simple series circuit of Example 6.1.

Example 6.1
For the circuit shown in Figure 6.1, derive an expression for v C  t  in terms of V p , R , C , and 

Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystem Modeling 61
Copyright © Orchard Publications
Chapter 6 Sinusoidal Circuit Analysis

I 1 = 10 45 and I 2 = 5 – 135 

and by addition,

I = I 1 + I 2 = 10 45 + 5 – 135  = 10  ------- + j ------- + 5  – ------- – j -------

 2 2  2 2
2 2 2 2

I = 5  ------- + j ------- = 5 45

or

 2 2
2 2

and finally transforming the phasor I into the t – domain , we obtain

i  t  = 5 cos  120  t + 45 

Also, for brevity, in our subsequent discussion we will designate resistive, inductive and capaci-
tive circuits as R , L , and C respectively.

6.4 Phasors in R, L, and C Circuits

The circuit analysis of circuits containing R , L , and C devices, and which are excited by sinusoi-
dal sources, is considerably simplified with the use of phasor voltages and phasor currents which
we will represent by the boldface capital letters V and I respectively. We will now derive V and
I phasor relationships in the j – domain . We must always remember that phasor quantities
exist only in the j – domain .

1. V and I phasor relationship in R branches

Consider circuit 6.4 (a) below where the load is purely resistive.

vS  t  + VS +
iR  t  vR  t  IR


VR
 

v R  t  = Ri R  t  = V p cos   t +   V R = RI R
 a  t – domain network  b  j – domain  phasor  network
Figure 6.4. Voltage across a resistive load in t – domain and j – domain

We know from Ohm’s law that v R  t  = Ri R  t  where the resistance R is a constant. We will
show that this relationship also holds for the phasors V R and I R shown in circuit 6.4 (b), that
is, we will prove that

68 Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling
Copyright © Orchard Publications
Chapter 6 Sinusoidal Circuit Analysis

v a = 5 sin t v b = – 15 sin t v ab = 20 sin t  = 2f = 0.4

Vs = AC Voltage Source, Is = AC Current Source, f= 0.2 Hz

VA VB
2 1

4 8
Vs +v +v Is
- -
1 VM 1 2 VM 2

Continuous

Scope1 Scope 2
Resistances in Ohms VB powergui
VA
Capacitors in microfarads
Inductor in millihenries
VM = Voltage Mezasurement

Scope 3
VA , VB
Subtract VA-VB
Bus
Creator

Figure 6.19. Simulink / SimPowerSystems model for the circuit in Example 6.13

Figure 6.20. Waveforms for v a v b and v a – v b

622 Circuit Analysis I with MATLAB Computing and Simulink / SimPowerSystems Modeling
Copyright © Orchard Publications
Chapter 7
Phasor Circuit Analysis

T
his chapter begins with the application of nodal analysis, mesh analysis, superposition, and
Thevenin’s and Norton’s theorems in phasor circuits. Then, phasor diagrams are intro-
duced, and the inputoutput relationships for an RC lowpass filter and an RC highpass
filter are developed.

7.1 Nodal Analysis

The procedure of analyzing a phasor* circuit is the same as in Chapter 3, except that in this chap-
ter we will be using phasor quantities. The following example illustrates the procedure.

Example 7.1
Use nodal analysis to compute the phasor voltage V AB = V A – V B for the circuit of Figure 7.1.

VA VB
2 j3 
 
4 8

5 0 A – j 6  – j 3  
10 0 A

Figure 7.1. Circuit for Example 7.1

Solution:
As in Chapter 3, we choose a reference node as shown in Figure 7.2, and we write nodal equa-
tions at the other two nodes A and B . Also, for convenience, we designate the devices in series
as Z 1 Z 2 and Z 3 as shown, and then we write the nodal equations in terms of these impedances.

Z 1 = 4 – j6 = 7.211 – 56.3

Z 2 = 2 + j3 = 3.606 56.3

Z 3 = 8 – j3 = 8.544 – 20.6

* A phasor is a rotating vector

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling 71
Copyright © Orchard Publications
Chapter 7 Phasor Circuit Analysis

– Z 1 I 1 +  Z 1 + Z 2 + Z 3 I 2 – Z 3 I 3 = 0

–  4 – j6  I 1 +  14 – j6 I 2 –  8 – j3 I 3 = 0 (7.10)

Also, by inspection for Mesh 3,

I 3 = 10 0 (7.11)

and in matrix form, (7.9), (7.10), and (7.11) are written as

1 0 0 I1 5
–  4 – j6   14 – j6  ––  8 – j3  I 2 = 0 (7.12)
0 0 1 I3 10

We use MATLAB for the solution of 7.12.*

Z=[1 0 0; (4j*6) 14j*6 (8j*3); 0 0 1];
V=[5 0 10]';
I=Z\V; i1=I(1); i2=I(2); i3=I(3); fprintf(' \n');
disp('i1 = '); disp(i1); disp('i2 = '); disp(i2); disp('i3 = '); disp(i3); fprintf(' \n');
i1 = 5 i2 = 7.5862 - 1.0345i i3 = 10
Therefore, the voltage across the 10 0 A current source is
V 10A = Z 3  I 2 – I 3  =  8 – j3   7.586 – j1.035 – 10  = – 22.417 – j1.038

We observe that this is the same value as that of the voltage V B in the previous example.

7.3 Application of Superposition Principle

As we know from Chapter 3, the superposition principle is most useful when a circuit contains
two or more independent voltage or current sources. The following example illustrates the appli-
cation of the superposition principle in phasor circuits.

Example 7.3
Use the superposition principle to find the phasor voltage across capacitor C 2 in the circuit of
Figure 7.5.

* As we experienced with Example 7.1, the computation of phasor voltages and currents becomes quite tedious. Accordingly,
in our subsequent discussion we will use MATLAB for the solution of simultaneous equations with complex coefficients.

76 Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems  Modeling
Copyright © Orchard Publications
 Chapter 8
Average and RMS Values, Complex Power, and Instruments

T
his chapter defines average and effective values of voltages and currents, instantaneous and
average power, power factor, the power triangle, and complex power. It also discusses elec-
trical instruments that are used to measure current, voltage, resistance, power, and energy.

8.1 Periodic Time Functions

A periodic time function satisfies the expression
f  t  = f  t + nT  (8.1)
where n is a positive integer and T is the period of the periodic time function. The sinusoidal and
sawtooth waveforms of Figure 8.1 are examples of periodic functions of time.

cos t cos  t +  

T T T T

T T T T

Figure 8.1. Examples of periodic functions of time

Other periodic functions of interest are the square and the triangular waveforms.

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems  Modeling 81
Copyright © Orchard Publications
Chapter 8 Average and RMS Values, Complex Power, and Instruments

Thus, the total average power supplied by the current sources is equal to the total average power
absorbed by the resistors. The small difference is due to rounding of fractional numbers.

8.9 Average Power in NonSinusoidal Waveforms

If the excitation in a circuit is nonsinusoidal, we can compute the average power absorbed by a
resistor from the relations

0 0 0 i R dt
T T T
1 1 v 2- 1 2
P ave = --- p dt = --- ----- dt = --- (8.35)
T T R T

Example 8.7
Compute the average power absorbed by a 5  resistor when the voltage across it is the half
wave rectification waveform shown in Figure 8.14.
Voltage (v)

10 sin t

1 2
t (ms)
Figure 8.14. Waveform for Example 8.7
Solution:
We first need to find the numerical value of  . It is found as follows:

= 2------ = 10 
–3 3
T = 2 ms = 2  10 s T = 2 
T
and thus 3
10 sin t = 10 sin 10 

0 0  10
Then,
–3 2 3 –3
1 T
v 2- 1 - 10 2
10 sin 10 -t + 2  10
P ave = --- ----- dt = -------------------
–3
--------------------------------- 0 dt
T R 2  10 5 –3

0 0 0
or
–3 –3 –3
100 10 1 10 10 3
P ave = ----------------------
- ---  1 – cos 2  10 3 t dt = 5  10 3 dt – cos 2  10  t dt
–3 2
10  10

 –3 
–3

=  5  10 t – ------------------------------   10 
10
3 3

 
3 sin 2  10 -t = 5  10 10
3 –3 sin 2  10
– ----------------------------------------------
 
3 3
2  10  2  10 
0

814 Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems  Modeling
Copyright © Orchard Publications
 Chapter 9
Natural Response

T
his chapter discusses the natural response of electric circuits.The term natural implies that
there is no excitation in the circuit, that is, the circuit is sourcefree, and we seek the cir-
cuit’s natural response. The natural response is also referred to as the transient response.

9.1 Natural Response of a Series RL circuit

Let us find the natural response of the circuit of Figure 9.1 where the desired response is the cur-
rent i, and it is given that at t = 0 , i = I0 , that is, the initial condition is i  0  = I 0 .

 +
R L
+ i 

Figure 9.1. Circuit for determining the natural response of a series RL circuit
Application of KVL yields
vL + vR = 0
or
di
L ----- + Ri = 0 (9.1)
dt
Here, we seek a value of i which satisfies the differential equation of (9.1), that is, we need to find
the natural response which in differential equations terminology is the complementary function. As
we know, two common methods are the separation of variables method and the assumed solution
method. We will consider both.
1. Separation of Variables Method

Rearranging (9.1), so that the variables i and t are separated, we obtain

di R
----- = – ---- dt
i L

Next, integrating both sides and using the initial condition, we obtain

0 
i1 t
--- di = – R
---- d
I i L 0

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling 91
Copyright © Orchard Publications
 Natural Response of a Series RC Circuit

and we sketch it as shown in Figure 9.11.

– 1   RC  t
e
v C  t   V 0 = – 1   RC  t + 1
Percent v C  t   V 0 
36.8%

13.5%
5%

Time constants
Figure 9.11. Circuit for determining the natural response of a series RC circuit
From Figure 9.11 we observe that at t = 0 , v C  V 0 = 1 , and i  0 as t  

The initial rate (slope) of decay is found from the derivative of v C  t   V 0 evaluated at t = 0 ,
that is,
d-  -----
vC 
dt  V 0 
---- - 1 - e – 1  RC  t
= – ------- 1-
= – -------
RC t=0
RC
t=0

and thus the slope of the initial rate of decay is – 1   RC 

Next, we define the time constant  as the time required for v C  t   V 0 to drop from unity to zero
assuming that the initial rate of decay remains constant. This constant rate of decay is repre-
sented by the straight line equation
vC  t  1
------------- = – -------- t + 1 (9.12)
V0 RC

and at t =  , v C  t   V 0 = 0 . Then,
1
0 = – --------  + 1
RC
or
 = RC
(9.13)
Time Cons tan t for R C Circuit

Evaluating (9.11) at t =  = RC , we obtain

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling 911
Copyright © Orchard Publications
 Chapter 10
Forced and Total Response in RL and RC Circuits

T
his chapter discusses the forced response of electric circuits.The term “forced” here implies
that the circuit is excited by a voltage or current source, and its response to that excitation
is analyzed. Then, the forced response is added to the natural response to form the total
response.

10.1 Unit Step Function u 0  t 

A function is said to be discontinuous if it exhibits points of discontinuity, that is, if the function
jumps from one value to another without taking on any intermediate values.

A wellknown discontinuous function is the unit step function u 0  t  * which is defined as

0
u0  t  = 
t0
1
(10.1)
t0

It is also represented by the waveform in Figure 10.1.

u0  t 
1
0
Figure 10.1. Waveform for u 0  t 

In the waveform of Figure 10.1, the unit step function u 0  t  changes abruptly from 0 to 1 at
t = 0 . But if it changes at t = t 0 instead, its waveform and definition are as shown in Figure
10.2.

0 t  t0
u0  t – t0  = 
1
u0  t – t0 
t 1 t  t0
0 t0
Figure 10.2. Waveform and definition of u 0  t – t 0 

* In some books, the unit step function is denoted as u  t  , that is, without the subscript 0. In this text we will reserve this des-
ignation for any input.

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling 101
Copyright © Orchard Publications
 Delta Function

The proofs of (10.22) through (10.24) and additional properties of the delta function are beyond
the scope of this book. They are provided in Signals and Systems with MATLAB Computing and
Simulink Modeling, ISBN 9781934404119.
MATLAB has two builtin functions for the unit step and the delta functions. These are desig-
nated by the names of the mathematicians who used them in their work. The unit step u 0  t  is
called Heavyside(t) and the delta function   t  is called Dirac(t). Shown below are examples of
how they are being used.
syms k a t
u=k*sym('Heaviside(ta)') % Create unit step function at t=a
u =
k*Heaviside(t-a)
d=diff(u) % Compute the derivative of the unit step function
d =
k*Dirac(t-a)
int(d) % Integrate the delta function
ans =
Heaviside(t-a)*k

Example 10.5
For the circuit shown in Figure 10.14, the inputs are applied at different times as indicated.

50 K
Rf

+ +
3 K 6 K
5 K v out
+ +

v in 1 
i in
v in 2


v in 1 = 0.8u 0  t – 3  V
v in 2 = 0.5u 0  t – 1  V
i in = 0.14  u 0  t + 1  + u 0  t – 2   mA
Figure 10.14. Circuit for Example 10.5

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling 109
Copyright © Orchard Publications
 Answers / Solutions to EndofChapter Exercises

R 1 + R 2
 R1  R2 
------------------
-
dv C
dt  1 – ----- 10-   C
--------- + ------------------------------ v C = 0
 R1 

and from this differential equation we see that the coefficient of v C is

1 1 30  216 30  216 15  9 135 5

 1 – 10 ------   1
------------- = --- = ----------------------------- = ------------------- = ------------------ = --------- = ------ = 0.3125
  
 18 
R eq C 8 18 4 108 432 16

and thus
– 0.3125t
v C  t  = Ae


and with the given initial condition v C  0  = V 0 = A = 5 V we obtain
– 0.3125t
v C  t  = 5e
Then, using the relation
dv C
i C = C ---------
dt
we find that for t  0
– 0.3125t – 0.3125t
i C  t  =  1   – 0.3125  5e  = – 1.5625e

and the minus () sign indicates that the i C  t  direction is opposite to that shown.

9.

CVS 10^(6) Continuous

s
+ +
Step - v powergui
Step time: 0 -
10^(-6) Scope
Initial value:0 VM
Final value: 12
VM=Voltage Measurement
CVS=Controlled Voltage Source

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems Modeling 1051
Copyright © Orchard Publications
Appendix D
A Review of Complex Numbers

T
his appendix is a review of the algebra of complex numbers. The basic operations are
defined and illustrated by several examples. Applications using Euler’s identities are pre-
sented, and the exponential and polar forms are discussed and illustrated with examples.

D.1 Definition of a Complex Number

In the language of mathematics, the square root of minus one is denoted as i , that is, i = – 1 .
In the electrical engineering field, we denote i as j to avoid confusion with current i . Essentially,
j is an operator that produces a 90degree counterclockwise rotation to any vector to which it is
applied as a multiplying factor. Thus, if it is given that a vector A has the direction along the
right side of the xaxis as shown in Figure D.1, multiplication of this vector by the operator j will
result in a new vector jA whose magnitude remains the same, but it has been rotated counter-
clockwise by 90  .

y
jA

j  j A  = j2 A = –A A
x
j  –j A  = –j A = A
2

j  –A  = j 3 A = –j A
Figure D.1. The j operator

Also, another multiplication of the new vector jA by j will produce another 90  counterclock-
wise direction. In this case, the vector A has rotated 180  and its new value now is – A . When
this vector is rotated by another 90  for a total of 270  , its value becomes j  – A  = – j A . A
fourth 90  rotation returns the vector to its original position, and thus its value is again A .
2 3 4
Therefore, we conclude that j = – 1 , j = – j , and j = 1 .

Circuit Analysis I with MATLAB  Computing and Simulink / SimPower Systems Modeling D1
Copyright © Orchard Publications
A Review of Complex Numbers

Example D.4
It is given that A = 3 + j 5 . Find A  A
Solution:
A  A =  3 + j 5   3 – j 5  = 3 + 5 = 9 + 25 = 34
2 2

D.4 Division of Complex Numbers

When performing division of complex numbers, it is desirable to obtain the quotient separated
into a real part and an imaginary part. This procedure is called rationalization of the quotient, and it
is done by multiplying the denominator by its conjugate. Thus, if A = a + jb and B = c + jd ,
then,
A a + jb  a + jb   c – jd  A B  ac + bd  + j  bc – ad 
---- = -------------- = ------------------------------------- = ----  ------- = ------------------------------------------------------
B c + jd  c + jd   c – jd  B B 2
c +d
2
(D.2)
 ac + bd - j ----------------------
bc – ad 
= ----------------------
2 2
+ 2 2
c +d c +d
In (D.2), we multiplied both the numerator and denominator by the conjugate of the denomina-
tor to eliminate the j operator from the denominator of the quotient. Using this procedure, we see
that the quotient is easily separated into a real and an imaginary part.

Example D.5
It is given that A = 3 + j 4 , and B = 4 + j 3 . Find A  B
Solution:
Using the procedure of (D.2), we obtain
A 3 + j 4  3 + j 4   4 – j 3  12 – j 9 + j 16 + 12 24 + j 7 24 7
---- = -------------- = -------------------------------------- = -------------------------------------------- = ----------------- = ------ + j ------ = 0.96 + j 0.28
B 4 + j3 4 + j34 – j3 4 +3
2 2 25 25 25

D.5 Exponential and Polar Forms of Complex Numbers

The relations

j
e = cos  + j sin  (D.3)

D4 Circuit Analysis I with MATLAB  Computing and Simulink / SimPower Systems Modeling
Copyright © Orchard Publications
Appendix E
Matrices and Determinants

T
his appendix is an introduction to matrices and matrix operations. Determinants, Cramer’s
rule, and Gauss’s elimination method are reviewed. Some definitions and examples are not
applicable to the material presented in this text, but are included for subject continuity,
and academic interest. They are discussed in detail in matrix theory textbooks. These are
denoted with a dagger (†) and may be skipped.

E.1 Matrix Definition

A matrix is a rectangular array of numbers such as those shown below.

1 3 1
2 3 7 or –2 1 –5
1 –1 5
4 –7 6

In general form, a matrix A is denoted as

a 11 a 12 a 13  a 1 n
a 21 a 22 a 23  a 2 n
A = a 31 a 32 a 33  a 3 n (E.1)
    
a m 1 a m 2 a m 3  a mn

The numbers a ij are the elements of the matrix where the index i indicates the row, and j indi-
cates the column in which each element is positioned. For instance, a 43 indicates the element
positioned in the fourth row and third column.
A matrix of m rows and n columns is said to be of m  n order matrix.
If m = n , the matrix is said to be a square matrix of order m (or n ). Thus, if a matrix has five
rows and five columns, it is said to be a square matrix of order 5.

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems  Modeling E1
Copyright © Orchard Publications
 Matrix Operations

Check with MATLAB:

A=[1 2 3; 0 1 4]; B=[2 3 0; 1 2 5]; % Define matrices A and B
A+B, AB % Add A and B, then Subtract B from A
ans =
3 5 3
-1 3 9
ans =
-1 -1 3
1 -1 -1
Check with Simulink:

A
3 5 3 Note: The elements of matrices
Constant 1
-1 3 9 A and B are specified in
MATLAB's Command prompt
B
Sum 1 Display 1 (A+B)
Constant 2
-1 -1 3
1 -1 -1

Sum 2 Display 2 (A-B)

If k is any scalar (a positive or negative number), and not  k  which is a 1  1 matrix, then mul-
tiplication of a matrix A by the scalar k is the multiplication of every element of A by k .

Example E.2
Multiply the matrix

A = 1 –2
2 3
by
a. k 1 = 5

b. k 2 = – 3 + j2

Circuit Analysis I with MATLAB  Computing and Simulink / SimPowerSystems  Modeling E3
Copyright © Orchard Publications
Index
Symbols and Numerics complex excitation function 6-3, 6-23 delta function
complex number(s) defined 10-7
% (percent) symbol in MATLAB A-2 addition B-2 sampling property 10-8
3-dB down 4-4 conjugate A-3, B-3 sifting property 10-9
defined A-3, B-2 demos in MATLAB A-2
A division B-4 dependent source(s)
exponential form B-5 current 1-11, 3-38
abs(z) MATLAB function A-23 multiplication B-3 voltage 1-11, 3-38
admittance 6-17 polar form B-5 determinant C-9
ampere 1-2, 1-19 rectangular form B-5 device(s)
ampere capacity of wires 2-30 subtraction B-2 active 1-11, 1-20
amplifier 4-1, 4-32 complex power 8-16 passive 1-11, 1-20
buffer 4-20 conductance 2-2 dielectric 5-16, 5-29
unity gain 4-13, 4-20 conj(A) MATLAB function C-8 differential input amplifier 4-5
analog-to-digital converter 8-28, 8-33 conjugate of a complex number B-3 digital filter 7-21
angle(z) MATLAB function A-23 conv(a,b) MATLAB function A-6 diode(s) 1-12
attenuation 4-13, 4-33 conversion factors 1-16 Dirac function 10-9
attenuator 4-1 conductor sizes for interior wiring 2-33 direct current 1-4
average value 8-2, 8-31 coulomb 1-1, 1-19 discontinuous function 10-1
axis MATLAB command A-16 Cramer’s rule 3-2, C-16, C-17 disp(A) MATLAB function 7-19, A-32
critical frequency 4-13, 4-33 display formats in MATLAB A-31
B current 1-1 division in MATLAB A-18
current division expressions 2-25 dot multiplication operator in MATLAB A-20
bandwidth 4-4 current flow driving functions 6-1
box MATLAB command A-12 conventional 1-2 duality 6-18, 6-25
branch 2-5 electron 1-2
current gain 4-2 E
C current limiting devices 2-2
current ratings for editor window in MATLAB A-1
capacitance 5-1, 5-17 electronic equipment 2-30 editor/debugger in MATLAB A-1
capacitance combinations 5-24 current source effective (RMS) value of sinusoids 8-5
capacitor(s) 1-11, 1-20, 5-16 combinations 2-14 effective values 8-4
in parallel 5-25 ideal 1-11 efficiency 3-44
in series 5-24 independent 1-11 eight-to-three line encoder 8-29
chemical processes 1-17, 1-20 practical 3-21 electric field 5-16, 5-17, 5-29
circuit(s) cutoff frequency electric filters - see filters
defined 1-13, 1-20 band-elimination filter 4-15 energy dissipated in a resistor 2-4
analysis with loop equations 3-8 band-pass filter 4-15 energy stored in a capacitor 5-21
analysis with mesh equations 3-8 high-pass filter 4-14 energy stored in an inductor 5-12
analysis with nodal equations 3-1 low-pass filter 4-13 eps in MATLAB A-22
with non-linear devices 3-42 lower 4-4 Euler’s identities B-4
clc MATLAB command A-2 upper 4-4 excitations 6-1
clear MATLAB command A-2 exit MATLAB command A-2
combined mesh 3-17 D exponential form of complex numbers B-5
combined node 3-6 exponentiation in MATLAB A-18
command screen in MATLAB A-1 data points in MATLAB A-14 eye(n) in MATLAB C-7
command window in MATLAB A-1 DC (Direct Current) 1-4 eye(size(A)) in MATLAB C-7
commas in MATLAB A-8 decibel 4-2, A-13
comment line in MATLAB A-2 deconv(c,d) MATLAB function A-6, A-7 F
comparators 8-29 default color in MATLAB A-15
complementary function 9-1 default in MATLAB A-12 Farad 5-17, 5-29
complete response 10-16 default line in MATLAB A-15 Faraday’s law of
complex conjugate A-4, B-3 default marker in MATLAB A-15 electromagnetic induction 5-2

IN-1
feedback 4-4 imaginary M
negative 4-5 axis B-2
positive 4-5 number B-2 magnetic field 5-1, 5-16, 5-29
figure window in MATLAB A-13 impedance 6-14 magnetic flux 5-2, 5-29
filter inductance 5-2 matrix, matrices
active 4-13 inductive adjoint C-20
all-pass 7-22 reactance 6-15, 6-23 cofactor of C-12
analog 7-23 susceptance 6-18, 6-23 conformable for addition C-2
band-elimination 4-15, 4-33, 7-22 inductor(s) conformable for multiplication C-4
band-pass 4-14, 4-33, 7-22 defined 1-11, 1-20, 5-2 congugate of C-8
band-rejection 4-15, 4-33, 7-22 in parallel 5-15 defined C-1
band-stop 4-15, 4-33, 7-22 in series 5-14 diagonal of C-1, C-6
high-pass 4-14, 4-33, 7-22 initial condition 5-3 Hermitian C-9
low-pass 4-13, 4-33, 7-22 initial rate of decay 9-3, 9-11 identity C-6
passive 4-13, 7-23 instantaneous values 2-1 inverse of C-21
phase shift 7-22 int(f,a,b) MATLAB function 1-7 left division in MATLAB C-24
RC high-pass 7-25 International System of Units 1-14 lower triangular C-6
RC low-pass 7-23 minor of C-12
stop-band 4-15, 4-33, 7-22 J multiplication using MATLAB A-20
flash converter 8-28 non-singular C-21
flux linkage 5-2, 5-29 j operator B-1 singular C-21
fmax(f,x1,x2) MATLAB function A-29 scalar C-6
fmin(f,x1,x2) MATLAB function A-29 K skew-Hermitian C-9
forced response 6-4, 10-16, 10-22 skew-symmetric C-9
format command in MATLAB A-31 KCL 2-6 square C-1
format in MATLAB A-31 Kirchhoff’s Current Law 2-6 symmetric C-8
fplot MATLAB command A-27 Kirchhoff’s Voltage Law 2-7 theory 3-2
fplot(fcn,lims) KVL 2-7 trace of C-2
MATLAB command A-27 transpose C-7
fprintf(format,array) L upper triangular C-5
MATLAB command 7-19, A-32 zero C-2
frequency response A-12 left-hand rule 5-1 maximum power
frequency-domain to time-domain lims = MATLAB function A-27 transfer theorem 3-35, 7-35
transformation 6-6, 6-23 linear mechanical forms of energy 1-17, 1-20
full-wave rectification circuit 3-38 mesh
function file in MATLAB A-26 devices 1-11 combined 3-18
fzero(f,x) MATLAB function A-26 factor A-9 defined 2-6
inductor 5-2 equations 2-10, 3-1, 5-25, 7-5
G passive element 3-37 generalized 3-17
linearity 3-37 mesh(x,y,z) MATLAB function A-18
Gaussian elimination method C-19 lines of magnetic flux 5-1, 5-29 meshgrid(x,y) MATLAB function A-18
grid MATLAB command A-12 linspace(values) MATLAB command A-14 metric system 1-14, 1-20
ground ln (natural log) A-13 m-file in MATLAB A-1, A-26
defined 2-1, 2-14 load mho 2-2
virtual 4-17 capacitive 8-15, 8-32 Military Standards 2-27
gtext(‘string’) MATLAB function A-13 inductive 8-15, 8-32 MINVERSE in Excel C-26
lighting 2-33 MMULT in Excel C-26, C-27
H resistive 8-11 multiplication of complex numbers B-3
log (common log) A-13 multiplication in MATLAB A-18
half-power points 4-4 log(x) MATLAB function A-13 multirange ammeter/milliammeter 8-24
half-wave rectification 8-3 log10(x) MATLAB function A-13
Heavyside function 10-9 log2(x) MATLAB function A-13 N
Henry 5-3, 5-29 loglog(x,y) MATLAB function A-13
loop NaN in MATLAB A-26
I defined 2-5 National Electric Code (NEC) 2-30
equations 3-1, 3-13 natural response
imag(z) MATLAB function A-23 circuits with single 2-10 9-1, 9-9, 10-16, 10-22

IN-2
NEC 2-30 complex 8-16, 8-17 series connection 2-8, 2-16, 2-17
negative charge 5-16 gain 4-2 short circuit 2-2
network in a capacitor 5-22 SI Derived Units 1-17
active 1-13, 1-20 in an inductor 5-11 siemens 2-2
passive 1-13, 1-20 in a resistor 2-3, 2-4, 2-28 signal 4-1, 4-32
topology 3-1 instantaneous 8-4 single ended output amplifier 4-5
newton 1-1, 1-19 power factor 8-10 single node-pair parallel circuit 2-14
nodal analysis 2-14, 3-1, 7-1 defined 8-10 slope converter 8-28
node lagging 8-15 solar energy 1-17, 1-20
combined 3-6 leading 8-15 sources of energy 1-17, 1-20
defined 2-5 power factor correction 8-18 standard prefixes 1-15
generalized 3-6 power triangle 8-16 Standards for Electrical and
equations 2-14, 3-2, 5-25, 7-1 prefixes 1-15, 1-16 Electronic Devices 2-26
non-reference 3-1 principle of superposition 3-41 steady-state conditions 5-12
reference 3-1 string in MATLAB A-18
non-linear devices 1-11 Q subplot(m,n,p) MATLAB command A-18
Norton’s theorem 3-33, 7-10 substitution method of solving a system
nuclear energy 1-17, 1-20 quad MATLAB function 1-8 of simultaneous equations 3-2
quad(‘f’,a,b,tol) MATLAB function 1-8 supermesh 3-17
O quad8 MATLAB function 1-8 supernode 3-6
quadratic factors A-9 superposition principle 3-38, 7-6
Ohm 2-1 quit MATLAB command A-2 susceptance
Ohm’s law 2-1 capacitive 6-18, 6-25
Ohm’s law for AC circuits 6-14 R inductive 6-18, 6-25
Ohmmeter 8-26
parallel type 8-26 rational polynomials A-8 T
series type 8-26 reactance
shunt type 8-26 capacitive 6-15, 6-24 temperature scales equivalents 1-16
op amp 4-5 inductive 6-15, 6-24 text(x,y,’string’) MATLAB function A-14
inverting mode 4-6 real text(x,y,z,’string’) MATLAB function A-16
non-inverting mode 4-9 axis B-2 Thevenin’s theorem 3-23, 7-10
open circuit 2-2 number B-2 time constant 9-3, 9-11, 10-18, 10-24
operational amplifier - see op amp real(z) MATLAB function A-23 time-domain to frequency-domain
regulation 3-45 transformation 6-5, 6-23
P resistance 2-1 time-window converter 8-28
input 4-28 title(‘string’) MATLAB command A-12
parallel connection 2-8, 2-17, 2-18 negative 2-3 total response 10-1, 10-14
particular solution 6-4 output 4-28 tracking converter 8-28
passive sign convention 1-9, 1-19 resistive network 8-29 transient response 9-1
periodic functions of time 8-1 resistors 1-11, 2-2 transistors 1-11
phasor analysis in amplifier circuits 7-14 color code 2-27 trivial solution 9-2
phasor diagram 7-17 failure rate 2-27 two-terminal device 1-4, 1-19
plot(x,y) MATLAB command A-10, A-12 shunt (parallel) 8-22
plot3(x,y,z) MATLAB command A-15 tolerance 2-27 U
polar plot in MATLAB A-24 response 6-1, 6-23
polar(theta,r) MATLAB function A-23 right-hand rule 5-1 unit impulse function 10-7
poly(r) MATLAB function A-4 RMS value of sinusoids 8-5 unit ramp function 10-6
polyder(p) MATLAB function A-6 RMS values of sinusoids with unit step function 10-1
polynomial construction from different frequencies 8-7
known roots in MATLAB A-4 roots(p) MATLAB function A-3, A-8 V
polyval(p,x) MATLAB function A-6 round(n) MATLAB function A-24
potential difference 1-4 virtual ground 4-17
power S volt 1-5, 1-19
absorbed 1-8, 1-19 voltage
average 8-9, 8-14 script file in MATLAB A-26 defined 1-4
in capacitive loads 8-11 semicolons in MATLAB A-8 dividers 2-2
in inductive loads 8-11 semilogx(x,y) MATLAB command A-12 division expressions 2-22
in a resistive loads 8-11 semilogy(x,y) MATLAB command A-12 drop 1-5

IN-3
 follower 4-20
gain 4-2
instantaneous 1-6
rise 1-5
voltage source
combinations 2-14
ideal 1-11
independent 1-11
practical 3-20
voltmeter 8-24

watt 1-8
watt-hour meter 8-28
wattage 2-4, 2-29
wattmeter 8-28
weber 5-1, 5-29
Wheatstone bridge 8-27, 8-32

xlabel(‘string’) MATLAB command A-12

ylabel(‘string’) MATLAB command A-12

zero potential 2-14

IN-4