

Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Matrix of a lin. transf T with respect to bases B and C: For every vector v in B, evaluate T(v), and express. T(v) as a linear combination of vectors in C.
Typology: Lecture notes
1 / 2
This page cannot be seen from the preview
Don't miss anything!


Subspace: If u and v are in W , then u + v are in W , and cu is in W
Nul(A): Solutions of Ax = 0. Row-reduce A.
Row(A): Space spanned by the rows of A: Row-reduce A and choose the rows that
contain the pivots.
Col(A): Space spanned by columns of A: Row-reduce A and choose the columns of A
that contain the pivots
Rank(A): = Dim(Col(A)) = number of pivots
Rank-Nullity theorem: Rank(A) + dim(N ul(A)) = n, where A is
m × n
Linear transformation: T (u + v) = T (u) + T (v), T (cu) = cT (u),
where c is a number.
T is one-to-one if T (u) = 0 ⇒ u = 0
T is onto if Col(T ) = R
m .
Linearly independence:
a 1 v 1
To show lin. ind, form the matrix of the vectors, and show that N ul(A) = { 0 }
Linear dependence: a 1 v 1
a 1 , a 2 , · · · , an, not all zero.
Span: Set of linear combinations of v 1 , · · · vn
Basis B for V : A linearly independent set such that Span (B) = V
To show sthg is a basis, show it is linearly independent and spans.
To find a basis from a collection of vectors, form the matrix A of the vectors, and find
Col(A).
To find a basis for a vector space, take any element of that v.s. and express it as a linear
combination of ’simpler’ vectors. Then show those vectors form a basis.
Dimension: Number of elements in a basis.
To find dim, find a basis and find num. elts.
Theorem: If V has a basis of vectors, then every basis of V must have n vectors.
Basis theorem: If V is an n−dim v.s., then any lin. ind. set with n elements is a basis,
and any set of n elts. which spans V is a basis.
Matrix of a lin. transf T with respect to bases B and C: For every vector v in B,
evaluate T (v), and express T (v) as a linear combination of vectors in C. Put the
coefficients in a column vector, and then form the matrix of the column vectors you
found!
Coordinates: To find [x] B
, express x in terms of the vectors in B.
x = P B
[x] B
, where P B
is the matrix whole columns are the vectors in B.
Invertible matrix theorem: If A is invertible, then: A is row-equivalent to I, A has n
pivots, T (x) = Ax is one-to-one and onto, Ax = b has a unique solution for
every b, A
is invertible, det(A) 6 = 0, the columns of A form a basis for R
n ,
N ul(A) = { 0 }, Rank(A) = n [
a b
c d
ad−bc
d −b
−c a
Change of basis: [x] C
[x] B (think of C as the new, cool basis)
is the matrix whose columns are [b] C
, where b is in B
Diagonalizability: A is diagonalizable if A = P DP
for some diagonal D and
invertible P.
A and B are similar if A = P BP
for P invertible
Theorem: A is diagonalizable ⇔ A has n linearly independent eigenvectors
Theorem: IF A has n distinct eigenvalues, THEN A is diagonalizable, but the opposite
is not always true!!!!
Notes: A can be diagonalizable even if it’s not invertible (Ex: A =
). Not
all matrices are diagonalizable (Ex:
Consequence: A = P DP
n = P D
n P
How to diagonalize: To find the eigenvalues, calculate det(A − λI), and find the
roots of that.
To find the eigenvectors, for each λ find a basis for N ul(A − λI), which you do by
row-reducing
Rational roots theorem: If p(λ) = 0 has a rational root r =
a
b
, then a divides
the constant term of p, and b divides the leading coefficient.
Use this to guess zeros of p. Once you have a zero that works, use long division! Then
, where D= diagonal matrix of eigenvalues, P = matrix of
eigenvectors
Complex eigenvalues If λ = a + bi, and v is an eigenvector, then
, where P =
Re(v) Im(v)
a b
−b a
C is a scaling of
det(A) followed by a rotation by θ, where:
det(A)
cos(θ) sin(θ)
− sin(θ) cos(θ)
u, v orthogonal if u · v = 0.
‖u‖ =
u · u
{u 1 · · · un} is orthogonal if u i
· u j
= 0 if i 6 = j, orthonormal if
u i
· u i
: Set of v which are orthogonal to every w in W.
If {u 1
· · · un} is an orthogonal basis, then:
y = c 1
u 1
y·u j
u j
·u j
Orthogonal matrix Q has orthonormal columns! Consequence:Q
= Orthogonal projection on Col(Q).
‖Qx‖ = ‖x‖
(Qx) · (Qy) = x · y
Orthogonal projection: If
u 1 · · · u k
is a basis for W , then orthogonal projection
of y on W is: ˆy =
y·u 1
u 1 u 1
u 1
y·u 1
u k u k
u k
y − ˆy is orthogonal to yˆ, shortest distance btw y and W is ‖y − ˆy‖
Gram-Schmidt: Start with B = {u 1
, · · · un}. Let:
v 1
= u 1
v 2
= u 2
u 2
·v 1
v 1 ·v 1
v 1
v 3
= u 3
u 3 ·v 1
v 1 ·v 1
v 1
u 3 ·v 2
v 2 ·v 2
v 2
Then {v 1
· · · vn} is an orthogonal basis for Span(B), and if w i
v i ∥ ∥ v i
then {w 1
· · · wn} is an orthonormal basis for Span(B).
QR-factorization: To find Q, apply G-S to columns of A. Then R = Q
Least-squares: To solve Ax = b in the least squares-way, solve A
Ax = A
b.
Least squares solution makes ‖Ax − b‖ smallest.
ˆx = R
b, where A = QR.
Inner product spaces f · g =
b
a
f (t)g(t)dt. G-S applies with this inner product
as well.
Cauchy-Schwarz: |u · v| ≤ ‖u‖ ‖v‖
Triangle inequality: ‖u + v‖ ≤ ‖u‖ + ‖v‖
Has n real eigenvalues, always diagonalizable, orthogonally diagonalizable
, P is an orthogonal matrix, equivalent to symmetry!).
Theorem: If A is symmetric, then any two eigenvectors from different eigenspaces are
orthogonal.
How to orthogonally diagonalize: First diagonalize, then apply G-S on each eigenspace
and normalize. Then P = matrix of (orthonormal) eigenvectors, D = matrix of
eigenvalues.
Quadratic forms: To find the matrix, put the x
i
-coefficients on the diagonal, and evenly
distribute the other terms. For example, if the x 1 x 2 −term is 6 , then the (1, 2)th and
(2, 1)th entry of A is 3.
Then orthogonally diagonalize A = P DP
Then let y = P
x, then the quadratic form becomes λ 1 y
n
where λ i are the eigenvalues.
Spectral decomposition: λ 1 u 1 u 1
Homogeneous solutions: Auxiliary equation: Replace equation by polynomial, so y
becomes r
etc. Then find the zeros (use the rational roots theorem and long division,
see the ‘Diagonalization-section). ’Simple zeros’ give you e
rt , Repeated zeros
(multiplicity m) give you Ae
rt
rt
m− 1 e
rt , Complex
zeros r = a + bi give you Ae
at cos(bt) + Be
at sin(bt).
Undetermined coefficients: y(t) = y 0 (t) + yp(t), where y 0 solves the hom.
eqn. (equation = 0), and yp is a particular solution. To find yp:
If the inhom. term is Ct
m e
rt , then:
yp = t
s (Amt
m · · · + A 1 t + 1)e
rt , where if r is a root of aux with
multiplicity m, then s = m, and if r is not a root, then s = 0.
If the inhom term is Ct
m e
at sin(βt), then: yp = t
s (Amt
m · · · +
t + 1)e
at cos(βt) + t
s (Bmt
m · · · + B 1 t + 1)e
rt sin(βt),
where s = m, if a + bi is also a root of aux with multiplicity m (s = 0 if not).
cos always goes with sin and vice-versa, also, you have to look at a + bi as one
entity.
Variation of parameters: First, make sure the leading coefficient (usually the coeff. of
y
) is = 1.. Then y = y 0
yp(t) = v 1 (t)y 1 (t) + v 2 (t)y 2 (t), where y 1 and y 2 are your hom.
solutions. Then
y 1 y 2
y
y
v
v
f (t)
. Invert the matrix and solve for v
and v
, and integrate to get v 1
and v 2
, and finally use:
yp(t) = v 1
(t)y 1
(t) + v 2
(t)y 2
(t).
Useful formulas:
a b
c d
ad−bc
d −b
−c a
sec(t) = ln |sec(t) + tan(t)|,
tan(t) = ln |sec(t)|,
∫
tan
(t) = tan(x) − x,
ln(t) = t ln(t) − t
Linear independence: f, g, h are linearly independent if
af (t) + bg(t) + ch(t) = 0 ⇒ a = b = c = 0. To show linear
dependence, do it directly. To show linear independence, form the Wronskian:
W (t) =
f (t) g(t)
f
(t) g
(t)
(for 2 functions),
W (t) =
f (t) g(t) h(t)
f
(t) g
(t) h
(t)
f
(t) g
(t) h
(t)
(for 3 functions). Then pick a point
t 0 where det(
W (t 0 )) is easy to evaluate. If det 6 = 0, then f, g, h are linearly
independent! Try to look for simplifications before you differentiate.
Fundamental solution set: If f, g, h are solutions and linearly independent.
Largest interval of existence: First make sure the leading coefficient equals to 1. Then
look at the domain of each term. For each domain, consider the part of the interval which
contains the initial condition. Finally, intersect the intervals and change any brackets to
parentheses. Harmonic oscillator: my
damping, k = stiffness)
To solve x
= Ax: x(t) = Ae
λ 1
t v 1
λ 2
t v 2
λ 3
t v 3 (λ i
are your eigenvalues, v i are your eigenvectors)
Fundamental matrix: Matrix whose columns are the solutions, without the constants
(the columns are solutions and linearly independent)
Complex eigenvalues If λ = α + iβ, and v = a + ib. Then:
x(t) = A
e
αt cos(βt)a − e
αt sin(βt)b
e
αt sin(βt)a + e
αt cos(βt)b
Notes: You only need to consider one complex eigenvalue. For real eigenvalues, use the
formula above. Also,
a+bi
a−bi
a
+b
Generalized eigenvectors If you only find one eigenvector v (even though there are
supposed to be 2 ), then solve the following equation for u: (A − λI)(u) = v
(one solution is enough).
Then: x(t) = Ae
λt v + B
te
λt v + e
λt u
Undetermined coefficients First find hom. solution. Then for xp, just like regular
undetermined coefficients, except that instead of guessing
xp(t) = ae
t
t
a 1
a 2
is
a vector. Then plug into x
= Ax + f and solve for a etc.
Variation of parameters First hom. solution x h
(t) = Ax 1
(t) + Bx 2
(t). Then
sps xp(t) = v 1
(t)x 1
(t) + v 2
(t)x 2
(t), then solve
W (t)
v
v
= f ,
where
W (t) =
x 1 (t) | x 2 (t)
. Multiply both sides by
W (t)
integrate and solve for v 1 (t), v 2 (t), and plug back into xp. Finally,
x = x h
Matrix exponential e
n=
n t
n
n!
. To calculate e
At , either
diagonalize: A = P DP
⇒ e
At = P e
Dt P
, where e
Dt is a
diagonal matrix with diag. entries e
λ i t
. Or if A only has one eigenvalue λ with
multiplicity m, use e
At = e
λt
m− 1
n=
(A−λI)
n t
n
n!
. Solution of
x
= Ax is then x(t) = e
At c, where c is a constant vector.
Case N = 2
Equation: x
= Ax, A =
Proper frequencies: Eigenvalues of A are: λ = − 1 , − 3 , then proper frequencies
±i, ±
3 i (± square roots of eigenvalues)
Proper modes: v 1
sin
π
3
sin
π
3
v 2
sin
π
3
sin
π
3
Case N = 3
Equation: x
= Ax, A =
Proper frequencies: Eigenvalues of A: λ = − 2 , − 2 −
2 , then
proper frequencies ±
2 i, ±
i, ±
i
Proper modes: v 1
sin
π
4
sin
π
4
sin
π
4
, v 2
sin
π
4
sin
π
4
sin
π
4
, v 3
sin
π
4
sin
π
4
sin
π
4
General case (just in case!)
Equation: x
= Ax,
Proper frequencies: ± 2 i sin
kπ
2(N +1)
, k = 1, 2 , · · · N
Proper modes: v k
sin
kπ
N +
sin
2 kπ
sin
N kπ
Full Fourier series: f defined on (−T , T ):
f (x) ˜
m=
am cos
πmx
T
πmx
T
, where:
a 0
f (x)dx
f (x) cos
πmx
b 0
f (x) sin
πmx
Cosine series: f defined on (0, T ): f (x) ˜
m=
am cos
πmx
where:
a 0
f (x)dx (not a typo)
f (x) cos
πmx
Sine series: f defined on (0, T ): f (x) ˜
m=
bm sin
πmx
, where:
b 0
f (x) sin
πmx
Tabular integration: (IBP:
f
g = f g −
f g
) To integrate
f (t)g(t)dt
where f is a polynomial, make a table whose first row is f (t) and g(t). Then
differentiate f as many times until you get 0 , and antidifferentiate as many times until it
aligns with the 0 for f. Then multiply the diagonal terms and do + first term − second
term etc.
Orthogonality formulas:
cos
πmx
sin
πnx
dx = 0
cos
πmx
T
cos
πnx
T
dx = 0 if m 6 = n
sin
πmx
T
sin
πnx
T
dx = 0 if m 6 = n
Convergence: Fourier series F goes to f (x) is f is continuous at x, and if f has a
jump at x, F goes to the average of the jumps. Finally, at the endpoints, F goes to
average of the left/right endpoints.
Heat/Wave equations:
Step 1: Suppose u(x, t) = X(x)T (t), plug this into PDE, and group X-terms
and T -terms. Then
(x)
X(x)
= λ, so X
= λX. Then find a differential
equation for T. Note: If you have an α-term, put it with T.
Step 2: Deal with X
= λX. Use boundary conditions to find X(0) etc. (if you
have
∂u
∂x
, you might have X
(0) instead of X(0)).
Step 3: Case 1: λ = ω
, then X(x) = Ae
ωx
−ωx , then find
ω = 0, contradiction. Case 2: λ = 0, then X(x) = Ax + B, then eihter find
X(x) = 0 (contradiction), or find X(x) = A. Case 3: λ = −ω
, then
X(x) = A cos(ωx) + B sin(ωx). Then solve for ω, usually ω =
πm
T
Also, if case 2 works, should find cos, if case 2 doesn’t work, should find sin.
Finally, λ = −ω
, and X(x) = whatever you found in 2) w/o the constant.
Step 4: Solve for T (t) with the λ you found. Remember that for the heat equation:
= λT ⇒ T (t) =
Ame
λt
. And for the wave equation:
= λT ⇒ T (t) =
Am cos(ωt) +
Bm sin(ωt).
Step 5: Then u(x, t) =
m=
T (t)X(x) (if case 2 works),
u(x, t) =
m=
T (t)X(x) (if case 2 doesn’t work!)
Step 6: Use u(x, 0), and plug in t = 0. Then use Fourier cosine or sine series or
just ‘compare’, i.e. if u(x, 0) = 4 sin(2πx) + 3 sin(3πx), then
= 3, and
Am = 0 if m 6 = 2, 3.
Step 7: (only for wave equation): Use
∂u
∂t
u(x, 0): Differentiate Step 5 with respect
to t and set t = 0. Then use Fourier cosine or series or ‘compare’
Nonhomogeneous heat equation:
∂u
∂t
= β
u
∂x
u(0, t) = U 1
, u(L, t) = U 2
u(x, 0) = f (x)
Then u(x, t) = v(x) + w(x, t), where:
v(x) = [
z
0
β
P (s)dsdz
x
L
x
0
z
0
β
P (s)dsdz
and w(x, t) solves the hom. eqn:
∂w
∂t
= β
w
∂x
w(0, t) = 0, w(L, t) = 0
u(x, 0) = f (x) − v(x)
D’Alembert’s formula: ONLY works for wave equation and −∞ < x < ∞:
u(x, t) =
(f (x + αt) + f (x − αt)) +
2 α
x+αt
x−αt
g(s)ds,
where u tt
= α
uxx, u(x, 0) = f (x),
∂u
∂t
u(x, 0) = g(x). The
integral just means ‘antidifferentiate and plug in’.
Laplace equation:
Same as for Heat/Wave, but T (t) becomes Y (y), and we get
(y) = −λY (y). Also, instead of writing
Y (y) =
Ame
ωy
Bme
−ωy , write
Y (y) =
Am cosh(ωy) +
Bm sinh(ωy). Remember cosh(0) = 1,
sinh(0) = 0