Engineering Mathematics Overview
Engineering Mathematics Overview
Prof N. Magaji
Syllabus
1 Linear algebra: linear vector spaces: transformations and matrices;
spectra theory; geometric concepts; structure of linear transformations;
quadratic forms, etc; network topology and graphs; matrix analysis
techniques; cut-set; incidence and nodal admittance matrix;
2 Complex variables. Phasor rotation by π/2 Complex numbers and
operations Polar form of a complex number Applications of complex
numbers to AC linear circuits Circular motionThe importance of being
exponential
3 Optimization: Introduction & formations of the problems, graphical
method, simplex method, duality concept in Linear Programing.
4 Statistical Quality Control: Introduction, control charts of all types, ISO
9000 series & their importance, OC curves, advantages & limitations of SQL
in industries.
5 Statistical distributions: Discrete & continuous frequency distributions,
Stochastic process, Markov chain, simulation techniques.
ELE8301(chapter 1) Prof. N. Magaji
Basic concepts
Vector in Rn is an ordered 1
set of n real numbers.
6
e.g. v = (1,6,3,4) is in R4 3
“(1,6,3,4)” is a column
4
vector:
as opposed to a row
vector: 1 6 3 4
m-by-n matrix is an object
with m rows and n columns,
1 2 8
each entry fill with a real
number: 4 78 6
9 3 2
3
ELE8301(chapter1) Prof. N. Magaji
Basic concepts
Transpose: reflect vector/matrix on line:
T T
a a b a c
a b
b c d b d
• Measure the magnitude of a vector
– Is the error in x small or large?
• General class of p-norms:
1
• 1-norm: x 1 i 1 xi
n
n p p
x x
• 2-norm: x2 n
i 1
xi
2 1/ 2
p i 1
• -norm: x
max i xi
4
ELE8301(chapter1) Prof. N. Magaji
Properties of matrices Norm 1
The length of a vector with n components,
[a1 a2 … an], is defined as ,
which is also called the 1 norm of [a1 a2 … an].
2 Norm 12 22 32 42 52
Example 1 Compute 2 2 1
a) 1 norm
b) infinity Norm
A 1 3 1
2 4 1
x 1 max 2 1 2 , 2 3 4 , 1 1 1 9
The maximum of the row sums
x
max 2 2 1 , 1 3 1 , 2 4 15 8
Properties of matrices (Linear Independence)
ELE8301(chapter1) Prof. N. Magaji
Non-square matrices
m=3
1 0 n=2 1 0 1
r=2 x1
0 1 System Ax=b may not 0 1 1
2 3 have a solution (x has x2
2 variables but 3 2 3 1
constraints).
m=2
1 0 2 n=3 x1
r=2 1 0 2 1
0 1 3 System Ax=b is x2
underdetermined (x 0 1 3 1
has 3 variables and 2 x3
constraints). 10
Orthogonal matrices
ELE8301(chapter1) Prof. N. Magaji
Orthogonal matrices
12
ELE8301(chapter1) Prof. N. Magaji
Determinants
If det(A) = 0, then A is singular.
If det(A) ≠ 0, then A is invertible.
To compute:
Simple example:
a b
Matlab: det(A) det ad bc
c d
m-by-n matrix A is rank-deficient if it has
rank r < m (≤ n)
Thm: rank(A) < r iff det(A) = 0 for all t-by-t
submatrices, r ≤ t ≤ m
13
Determinants
ELE8301(chapter1) Prof. N. Magaji
14
Computation of Eigenvalues and
ELE8301(chapter1) Prof. N. Magaji
Eigenvectors
Let A be an n n matrix with eigenvalue and corresponding
eigenvector x. Thus Ax = x. This equation may be written
Ax – x = 0
given
(A – In)x = 0
Solving the equation |A – In| = 0 for leads to all the eigenvalues
of A.
On expending the determinant |A – In|, we get a polynomial in .
This polynomial is called the characteristic polynomial of A.
The equation |A – In| = 0 is called the characteristic equation of
A.
15
ELE8301(chapter1) Prof. N. Magaji
1 4 5
A 0 3 / 4 6
0 0 1 / 2
1 4 5
det( A I ) 0 3/ 4 6 (1 )(3 / 4 )(1 / 2 )
0 1 / 2
0
1, 3 / 4, 1 / 2
16
Eigenvalues & eigenvectors
ELE8301(chapter1) Prof. N. Magaji
Example 6 (cont.)
1 = 10
We get ( A 10 I 3 )x 0
5 4 2 x1
4 5 2 x2 0
2 2 8 x3
2 = 1
Let = 1 in (A – I3)x = 0. We get
( A 1I 3 )x 0
4 4 2 x1
4 4 2 x2 0
2 2 1 x3
The solution to this system of equations can be shown to be
x1 = – s – t, x2 = s, and x3 = 2t, where s and t are scalars.
Thus the eigenspace of 2 = 1 is the space of vectors of the
form.
s t
s
2t
21
ELE8301(chapter1) Prof. N. Magaji
Definition
Let A and B be square matrices of the same size. B is said to be
similar to A if there exists an invertible matrix C such that
B = C–1AC. The transformation of the matrix A into the matrix B
in this manner is called a similarity transformation.
Theorem 2.1
Similar matrices have the same eigenvalues.
Definition
A square matrix A is said to be diagonalizable if there exists a
matrix C such that D = C–1AC is a diagonal matrix.
23
Example 7
ELE8301(chapter1) Prof. N. Magaji
Let A be an n n matrix.
(a) If A has n linearly independent eigenvectors, it is
diagonalizable. The matrix C whose columns consist of n
linearly independent eigenvectors can be used in a similarity
transformation C–1AC to give a diagonal matrix D. The
diagonal elements of D will be the eigenvalues of A.
(b) If A is diagonalizable, then it has n linearly independent
eigenvectors
Example 8
4 6
A
(a) Show that the matrix 3 5 is diagonalizable.
(b) Find a diagonal matrix D that is similar to A.
(c) Determine the similarity transformation that diagonalizes A.
25
Example 8(cont.)
ELE8301(chapter1) Prof. N. Magaji
Solution
(a) The eigenvalues and corresponding eigenvector of this
matrix were found to be :
1 2
1 2, v1 r 2 1, and v 2 s
1 1
Since A, a 2 2 matrix, has two linearly independent
eigenvectors, it is diagonalizable.
(b) A is similar to the diagonal matrix D, which has diagonal
elements 1 = 2 and 2 = –1. Thus
4 6 2 0
A is similar to D
3 5 0 1
26
(c) Select two convenient linearly independent eigenvectors, say
1 2
ELE8301(chapter1) Prof. N. Magaji
v1 and v 2
1 1
Let these vectors be the column vectors of the diagonalizing
matrix C. 1 2
C
1 1
We get
1
1 1 2 4 6 1 2
C AC
1 2 3 5 1 1
1 2 4 6 1 2 2 0
D
1 1 3 5 1 1 0 1
A9 CD9C 1
1
1 2 512 0 1 2 514 1026
1 1 0 1 1 1 513 1025
28
ELE8301(chapter1) Prof. N. Magaji
Special matrices
Matrix A is symmetric if A = AT
A is positive definite if xTAx>0 for all non-zero x (positive semi-
definite if inequality is not strict)
1 0 0 a
a b c 0 1 0 b a 2 b 2 c 2
0 0 1 c
1 0 0 a
a b c 0 1 0 b a 2 b 2 c 2
0 0 1 c
Useful fact: Any matrix of form ATA is positive semi-
[Link] see this, xT(ATA)x = (xTAT)(Ax) = (Ax)T(Ax) ≥ 0
29
Quadratic Forms.
ELE8301(chapter1) Prof. N. Magaji
j 1 k 1
Solution:
(b) x – y = 0 x = y
Suppose y = t, z = s
x t 1 0 1 0
Then y t 1 t 0 s Hence B = 1 , 0
z 0 1
0 1
s
x 2t 2 2
(c) y t 1 t Hence B = 1
4
z
4t
4
a a 1 0 1 0
b a c 1 a 1 c Hence B 1 , 1
(d)
0 1
c
c
0
1
TUTORIAL
ELE8301(chapter1) Prof. N. Magaji
2 1 3
Q2 If A
4 8 6 , find A
1
0 7 5
Solution
2 16 18
From the example 2.23, adj A 20 10 0 A 60
and
28 14 12
2 16 18
1
Then, A1
20 10 0
60
28 14 12
z T = zH
Complex conjugate of a matrix A is a matrix derived from A by taking the complex conjugate
of all the elements in A.
3 i 2 3 i 2
A= A = 2i 1 i
2i 1 i
5 i t i 0 5 i t i 0
A= A =
12 5i t 2 12
5i t 2
5 i 12
AH = A* = A = t i 5i
T
0 t 2
Q5 : Compute the 1,-2-, -, and Frobenious norms of the matrix
TUTORIAL
ELE8301(chapter1) Prof. N. Magaji
0 3 1
(1) Prove that A 3 0 5 is a skew-symmetric matrix.
1 5 0
(2) (a) Show that if A is a 3x3 matrix such that A t A , then detA=0.
1 2 74
(b) Given that B 2 1 67 ,
74 67 1
use (a) , or otherwise , to show det(I B) 0 .
Hence deduce that det(I B 4 ) 0 .
(3) (a) If , and are the roots of x 3 px q 0 , find a cubic equation whose
roots are 2 , 2 and 2 .
x 2 3
(b) Solve the equation 2 x 3 0 .
2 3 x
Hence, or otherwise, solve the equation
x 3 38x 2 361x 900 0 .
ELE8301(chapter 1) Prof. N. Magaji
ELE8301
Chapter 1
Lecture 2
Outline
1 Network topology and graphs;
2 matrix analysis techniques;
3 cut-set matrix analysis
4 Loop incidence matrix and nodal admittance matrix;
Network topology
ELE8301(chapter1) Prof. N. Magaji
37
Network topology
ELE8301(chapter1) Prof. N. Magaji
38
ELE8301(chapter1) Prof. N. Magaji 1 Incidence matrix
Incidence matrix is that matrix which
represents the graph such that with the
help of that matrix we can draw a graph.
This matrix can be denoted as [AC] As in
every matrix, there are also rows and
columns in incidence matrix [AC].
The rows of [AC] represent the number of
nodes
Column of the matrix [AC] represent the
number of branches in the given graph.
The entries of incidence matrix is always -
1, 0, +1.
Steps to draw the incidence matrix :-
If a given kth node has outgoing branch, Fig. 3
then we will write +1.
If a given kth node has incoming branch,
then we will write -1. Figure 3 has 4 nodes and
Rest other branches will be considered 0. 6 branches. Thus the
incidence matrix for this
figure will be 4 rows and
6 columns 39
ELE8301(chapter1) Prof. N. Magaji 1 Incidence matrix(cont.)
For Figure 3 graph
the incidence matrix
is written as:
Fig. 3
40
ELE8301(chapter1) Prof. N. Magaji 2 Loop Analysis
Consider a connected graph with b branches and nt nodes.
Pick a tree T There are n = nt-1 tree branches and l = b-nt
links. Number the links first to be 1,2….l and number the tree
from l+1 to b . Every link and a unique path of tree
branches
The graphdefines
of Fig. a4 fundamental
illustrates loop.
Assign the direction of loop
fundamental loop for the chosen current to the same as the
Tree direction of the link the KVL for
each fundamental loop are.
42
ELE8301(chapter1) Prof. N. Magaji 2 Loop Analysis(cont.)
In matrix form
LV=0
Properties of the loop impedance matrix
For a RLC networks in sinusoid steady state the loop impedance
matrix
Z i(jw) = BZb(jw)BT and has the following properties
If there is no coupling element the matrix Zb(jw) is diagonal and
the loop impedance matrix is symmetric.
If there is no coupling element the matrix Zb(jw) can be written by
inspection
Zii ( jw) is the sum of impedance in the loop i and
Zjk(jw) is the sum or negative sum of impedance of branch i
impedance common to loop j the plus sign applied
if the branch i direction agree with the loop j direction
43
ELE8301(chapter1) Prof. N. Magaji 3 Cut set Analysis
Cut set analysis is a dual of loop analysis
Every tree branch defines a unique cut set
The fundamental cut set of the circuit of
Fig.5 is shown in Fig.6
KCL can be written for each
cut set as shown
45
Tutorial
ELE8301(chapter1) Prof. N. Magaji
1 2
d 3 c
a
b
Solution
(a) No. of branches, b = n – 1 = 4 – 1 = 3
No. of links, _ e b e (n 1) e n
1 5 3 2
46
Tutorial
ELE8301(chapter1) Prof. N. Magaji
47
Tutorial
ELE8301(chapter1) Prof. N. Magaji
Example 2 Fig. 2
Write the fundamental loop equation
for the circuit shown in Fig.2
Fig. 2
48
Tutorial
ELE8301(chapter1) Prof. N. Magaji
Example 3
Write the cut set equation of Fig. 3 by inspection.
Fig. 3
49