Qualifying Exams Study Sheet: Jesse Adams August 10, 2014
Qualifying Exams Study Sheet: Jesse Adams August 10, 2014
Jesse Adams
August 10, 2014
Contents
1) Previous Material
1.a) Calculus . . . . . . . . . . . .
1.a.i) Fundamental Theorem
1.a.ii) Sums . . . . . . . . .
1.a.iii) Trig Substitutions . .
1.a.iv) Integrals . . . . . . . .
1.a.v) Sequences . . . . . . .
1.a.vi) Series . . . . . . . . .
1.a.vii) Vectors . . . . . . . .
1.a.viii)Coordinate Systems .
1.a.ix) Integration Elements .
1.a.x) Multiple Integrals . .
1.a.xi) Vector Fields . . . . .
1.a.xii) Trig Identities . . . . .
1.b) Differential Equations . . . .
. . . . . . .
of Calculus
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
3
3
3
3
3
3
3
4
4
4
4
5
5
5
6
2) Numerical Analysis
2.a) Linear Algebra . . . . . . . . . . . . .
2.a.i) Basics . . . . . . . . . . . . . .
2.a.ii) Norms . . . . . . . . . . . . . .
2.a.iii) Singular Value Decomposition .
2.a.iv) Projectors . . . . . . . . . . . .
2.a.v) QR Factorization . . . . . . . .
2.a.vi) Least Squares . . . . . . . . . .
2.a.vii) Conditioning . . . . . . . . . .
2.a.viii)Stability . . . . . . . . . . . . .
2.a.ix) Gaussian Elimination . . . . .
2.a.x) Cholesky Factorization . . . . .
2.a.xi) Eigenvalues . . . . . . . . . . .
2.a.xii) Iterative Methods . . . . . . .
2.b) Numerical Methods . . . . . . . . . . .
2.b.i) Functional Iteration . . . . . .
2.b.ii) Polynomial Interpolation . . .
2.b.iii) ODEs . . . . . . . . . . . . . .
2.b.iv) Basic Concepts . . . . . . . . .
2.b.v) Methods . . . . . . . . . . . . .
2.b.vi) BVPs . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
6
6
6
6
7
7
7
8
8
9
9
9
10
10
11
11
11
12
12
12
13
3) Analysis
3.a) Metric Spaces . . . . . . . . .
3.a.i) `p Spaces . . . . . . .
3.a.ii) Lebesgue (Lp ) Spaces
3.a.iii) Normed Linear Spaces
3.b) Topology . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
14
14
14
15
15
15
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
3.c)
3.d)
3.e)
3.f)
3.g)
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
15
16
16
16
17
17
18
18
18
18
19
19
19
19
20
20
20
20
21
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
21
21
21
21
22
22
22
22
23
23
23
23
23
24
24
24
25
25
25
25
25
26
26
26
1)
Previous Material
1.a)
Calculus
1.a.i)
n
X
k=1
2.
Sums
n
X
n(n + 1)
k=
;
2
k=1
n
X
n(n + 1)(2n + 1)
;
k =
6
2
k =
k=1
n(n + 1)
2
2
X
2
1
=
2
n
6
n=1
a
1. b2 x2 a2 x = sec ()
b
a
2. a2 b2 x2 x = sin ()
b
a
3. a2 + b2 x2 x = tan ()
b
4. x = tan (/2) sin() =
1.a.iv)
2x
1 x2
2 dx
,
cos()
=
, d =
2
2
1+x
1+x
1 + x2
Integrals
R
x=b
dy
dx
2
y=d
dx
dy
2
t=f
s
dx
dt
2
dy
dt
2
1. Arc Length: L = ds =
1+
dx =
1+
dy =
+
x=a
y=c
t=e
s
2
Z =h
dr
2
r +
d
d
=g
R
R
2. Surface Area: A = 2y ds about x-axis, A = 2x ds about y-axis, with ds as defined in arc length.
1.a.v)
dt =
Sequences
R
1. Integral test: Given continuous, positive,
and decreasing f (x) on [k, ) with f (n) = an , then if k f (x) dx
P
is convergent/divergent, then so is n=k an .
P
P
P
2. P
Comparison test: Given two series with 0 an bn n, then
bn <
an < ;
an =
bn = .
3. Limit comparison test: c = lim
an
. If 0 < c < , then either both converge or both diverge.
bn
4. Alternating series test: Given an = (1)n bn or an = (1)n+1 bn with bn 0, if lim bn = 0 and {bn }
n
P
decreasing, then
an is convergent.
an+1
. L < 1 convergence, L > 1 divergence.
5. Ratio test: L = lim
n
an
p
6. Root test: L = lim n |an |. L < 1 absolute convergence, L > 1 divergence.
n
7. Absolute convergence:
1.a.vi)
Series
1. Power series:
n
1
1 X f (x)
provided |f (x)| < |a|, and a 6= 0.
=
a f (x)
a n=0
a
X
f (n) (x0 )
(x x0 )n
n!
n=0
n
3. Binomial series: (a + b) =
n
X
n
k=0
ank bk
X
k n
k
k(k 1) (k n + 1)
k
4. (1 + z) =
.
z for |z| < 1, and
=
n!
n
n
n=0
5. ez =
1.a.vii)
X
zn
;
n!
n=0
sin(z) =
z 2n+1
(1)n ;
(2n
+
1)!
n=0
cos(z) =
X
z 2n
(1)n ;
(2n)!
n=0
ln(z) =
X
(z 1)n
(1)n+1
n
n=1
Vectors
r0 (t)
.
||r0 (t)||
T0 (t)
.
||T0 (t)||
Coordinate Systems
x = r cos(), y = r sin(), z = z.
2 = r 2 + z 2 ;
Integration Elements
x1
u1
(x)
du = x2
1. In general, take the determinant: dx =
u1
u
x3
u1
2. Cylindrical:
(a)
(b)
(c)
(d)
Constant radius: dA = r d dz
Constant angle: dA = dr dz
Constant height: dA = r dr d.
Volume: dV = r dr d dz
3. Spherical:
(a)
(b)
(c)
(d)
x1
u2
x2
u2
x3
u2
x1
u3
x2
u3
x3
u3
du
1.a.x)
Multiple Integrals
ZZ
ZZ
1. Change of variables:
f (x, y) dA =
D
(x, y)
dudv (similarly for triple integrals).
f (g(u, v), h(u, v))
(u, v)
S
Vector Fields
dA.
P dx + Q dy =
x
y
C
D
Then
j
k
z
R
P
Q R
+
+
x
y
z
7. Stokes Theorem: Given smooth surface S bounded by simple, closed, smooth curve C and vector field F,
Z
ZZ
F dr =
curl(F) dS
C
8. Divergence Theorem: Given simple solid region E, boundary surface S, and vector field F with continuous
1st order partials,
ZZ
ZZZ
F dS =
div(F)dV
S
1.a.xii)
Trig Identities
1.b)
2)
2.a)
2.a.i)
Differential Equations
Numerical Analysis
Linear Algebra
Basics
Norms
1. Must satisfy
(a) ||x|| 0, with ||x|| = 0 x = 0
(b) ||x|| = || ||x||
(c) ||x + y|| ||x|| + ||y||
2. All vector norms on Cn are equivalent, i.e. 0 < c1 < c2 < : c1 |||| |||| c2 ||||
3. Induced Matrix Norms: ||A||(m,n) = sup
xCn
x6=0
||Ax||(m)
||x||(n)
sup
||x||(n) =1
||Ax||(m)
1/2
X
2
5. Frobenius Norm: ||A||F =
|aij |
=
i,j
!1/2
X
2
||ai ||2
p
= trace(A A) =
sX
i2
P
6. Nuclear Norm: ||A|| = trace( A A) = i i
7. Unitary matrices Q satisfy ||QA||2 = ||A||2 , ||QA||F = ||A||F
2.a.iii)
V , with U
Cmn ,
Rnn , and V Cnn ; U
, V unitary, and
diagonal with
1. Reduced SVD: A = U
0
decreasing elements.
2. Full SVD: A = U V with
U=
V =
V
(V )
,
0(mn)n
1 .
=
1 AV , or V = A U
(c) Find the other matrix: U
4. Properties:
(a) r = rank(A) = #{j > 0}.
(b) range(A) = spanhu1 , . . . , ur i, null(A ) = spanhur+1 , . . . , um i.
(c) null(A) = spanhvr+1 , . . . , vn i, range(A ) = spanhv1 , . . . , vr i.
Pr
(d) A = j=1 j uj vj (for rank k < r approx, use the first k; gives min error in 2 and Frobenius norms).
2.a.iv)
Projectors
1. Projector: P 2 = P
2. Complimentary Projector: I P ; range(I P ) = null(P ) and null(I P ) = range(P )
3. Orthogonal Projector: P 2 = P and P = P .
aa
; Pa = I Pa .
a a
(b) Onto range of A (arbitrary basis): P = A(A A)1 A = AA+
(a) Rank 1 orthogonal projectors (a 6= 0): Pa =
2.a.v)
QR Factorization
,
R=
0(mn)n
4. Modified Gram-Schmidt: A R1 R2 Rn = Q
|
{z
}
1
R
vv
I 0
where F = I 2 .
5. Householder Triangularization: Qn Qn1 Q1 A = R, with Qk =
0 F
{z
}
|
v v
6. Algorithms: Pseudocode.
Classical Gram-Schmidt:
Modified Gram-Schmidt:
for j = 1:n
Q[ : , j ] = A[ : , j ]
f o r i = 1 : j 1
R[ i , j ] = Q [ : ,
Q[ : , j ] = Q[ : ,
R[ j , j ] = norm (Q [ :
Q[ : , j ] = Q[ : , j ]
for i = 1:n
R[ i , i ] = norm (A [ :
Q[ : , i ] = A[ : , i ]
f o r j = i +1:n
R[ i , j ] = Q [ : ,
A[ : , j ] = A[ : ,
i ] A[ : , j ]
j ] R[ i , j ] Q [ : ,
, j ] , 2)
/ R[ j , j ]
i]
, i ] , 2)
/ R[ i , i ]
i ] A[ : , j ]
j ] R[ i , j ] Q [ : ,
i]
Householder QR:
% Note t h a t t h e o ut pu t i s R = A, and W m a t r i c e s
for i = 1:n
x = A[ i : , i ]
x [ 1 ] = s i g n ( x [ 1 ] ) norm ( x ) + x [ 1 ]
W[ i : , i ] = x / norm ( x , 2 )
A[ i : , i : ] = A[ i : , i : ] 2 W[ i : , i ] (W[ i : , i ] A[ i : , i : ] )
% To s o l v e , n o t e t h a t
% Rx = Qb ,
for i = 1:n
b [ i : ] = b [ i : ] 2 W[ i : , i ] (W[ i : , i ] b [ i : ] )
% Use back s u b s t i t u t i o n t o s o l v e f o r x
x = R\b
2.a.vi)
Least Squares
1. Normal Equations: A Ax = A b
= V
1 U
1 Q
2. Pseudoinverse: A+ = (A A)1 A = R
3. Solution methods:
(a) Cholesky: A A = R R, where R is upper triangular (requires full rank A). Best for speed, bad for error.
(b) QR: Reduces to Rx = Q b, use back substitution. Good method unless A is close to rank deficient.
V x = U
b, then solve. Stable even if A close to rank deficient, but more time/memory consuming.
(c) SVD:
2.a.vii)
Conditioning
||f ||
= ||J(x)||
||x||
||f || ||x||
||J(x)||
2. Relative condition number: = sup
/
=
||f
(x)||
||x||
||f
(x)||
/ ||x||
x
3. Matrix-vector condition: =
||A||||x||
||Ax||
4. Matrix condition number: (A) = ||A|| A1 or ||A|| ||A+ || in rank deficient case.
2.a.viii)
Stability
f
(
x
)
||
x x||
= O(m ) and
= O(m )
2. Stability:
x:
||x||
||f (x)||
3. Backward stability: f(x) = f (
x) for some x
as above.
4. Both Householder triangularization is backward stable, and so is MGS provided Q b is formed implicitly.
2.a.ix)
Gaussian Elimination
1
1
1. No pivoting: Use for hand calcs, not stable. A = LU where L = L1
1 L2 Lm1 is lower triangular, U is
xjk
upper triangular. We have `jk = xkk for k < j m, and
Lk =
..
.
1
`k+1,k
..
.
1
..
`mk
`21
L=
.
..
`m1
0
1
..
.
..
.
..
.
`m,m1
0
..
.
0
1
3. Full pivoting: (L0m1 L02 L01 )(Pm1 P2 P1 )A(Q1 Q2 Qm1 ) = L1 P AQ = U . Here, L and P are as
before, and Q = Q1 Q2 Qm is another permutation matrix. Solve for P AQ = LU .
4. Algorithms:
GE No Pivot:
GE Partial Pivot:
U = A; L = I
f o r k = 1 :m1
f o r j = k +1:m
L [ j , k ] = U[ j , k ] / U[ k , k ]
U[ j , k : ] = U[ j , k : ] L [ j , k ] U[ k , k : ]
U = A; L = I ; P = I
f o r k = 1 :m1
f o r j = k +1:m
S e l e c t i >= k : abs (U[ i , k ] ) i s maximized
U[ k , k : ] , U[ i , k : ] = U[ i , k : ] , U[ k , k : ]
P[ k , : ] , P[ i , : ] = P[ i , : ] , P[ k , : ]
f o r j = k +1:m
L [ j , k ] = U[ j , k ] / U[ k , k ]
U[ j , k : ] = U[ j , k : ] L [ j , k ] U[ k , k : ]
2.a.x)
Cholesky Factorization
0 1
= R1 Rn I Rn R1
A = 11
=
w K
w/ I
0
I
0 K ww
| {z } | {z }
a11
{z
}|
{z
}
|
R
{z
}|
R
R1
where =
R1
A1
a11 .
2. Algorithm:
Cholesky:
R = A
for k = 1:n
f o r j = k +1:n
R[ j , j : ] = R[ j , j : ] R[ k , j : ] R[ k , j : ] / R[ k , k ]
R[ k , k : ] = R[ k , k : ] / s q r t (R[ k , k ] )
2.a.xi)
Eigenvalues
Iterative Methods
TJ
4. Algorithms:
Jacobi:
D = d i a g ( d i a g (A) )
for i = 1: cIter
x = D \ ( (D A) x + b )
Gauss-Seidel:
D = d i a g ( d i a g (A ) ) ; L = t r i l (A ) ; U = t r i u (A)
for i = 1: cIter
x = (D L) \ ( b + Ux)
10
cGS
2.b)
Numerical Methods
2.b.i)
Functional Iteration
Casting nonlinear system as a fixed point problem: x = g(x) for x Rn , and g(x) = (g1 (x), g2 (x), . . . , gn (x))> .
Iterating x(k+1) = g(x(k) ).
1. Contraction
Mapping:
If ||g(x) g(y)|| ||x y|| x, y : x x(0) , y x(0) with 0 <
1, and g(x(0) x(0) (1 ), then limk x(k) = a where g(a) = a, and a unique (in this region).
gi (x)
st
i, j = 1 : n and x B (a) = {x Rn : ||x a|| },
2. If gi (x) has continuous 1 order s:
xj
n
then x(0) B (a) x(k) a (unique).
3. Nonlinear Systems: g(x) = x A(x)f (x).
(a) Easy: A(x) = A.
(b) Newton: A(x) = J 1 (x) (inverse Jacobian of f ).
2.b.ii)
Polynomial Interpolation
1. Lagrange: Pn (x) =
n
X
j=0
(x xi )/
i6=j
Y
(xj xi ), where xj , j = 0 : n are known
i6=j
points.
def
(xxj ) (n+1)
()
j=0 (n+1)! f
Qn
n
X
n (x) (n+1)
f
()
(n + 1)!
j=0
|
f (k) (xi )
n
X
{z
approx
n (x) =
n
Y
(x xj )
j=0
n0 (xi ) (n+1)
f
((xi ))
(n + 1)!
{z
}
|
err
(k)
j=0
f (x) dx
a
Z
n
X
j=0
nj (x) dx f (xj )
Rb
Pn
5. Weighted Least Squares: f L2 ([a, b]). Then Qn (x) = j=0 cj Pj (x) with cj = a f (x)Pj (x)w(x) dx miniRb
Rb
2
2
mizes ||f (x) Qn (x)||2 = a |f (x) Qn (x)| w(x) dx, with inner product hf (x), g(x)iw = a f (x)g(x)w(x) dx.
Rb
w(x) 0 on [a, b], and a w(x) dx > 0.
6. G-S Orthonormalization: Given {gi (x)}ni=0
d0 = 1/ ||g0 ||L2
f0 (x) = d0 g0 (x)
f1 (x) = d1 [g1 (x) c01 f0 (x)]
n1
X
fn (x) = dn gn (x)
cjn fj (x)
j=0
11
n
X
1
cn einx + cn einx +
cj eijx
2
j=n
1
cj =
2n
2.b.iii)
n
X
1
f (xn )eijxn + f (xn eijxn +
f (xk )eijxk
2
!
j = n : n
k=n
ODEs
Basic Concepts
Methods
hn
(f (tn1 , yn1 ) + f (tn , yn ))
2
0
4. Taylor Series Method: y0 = f (t, y), yn = yn1 + hyn1
+
12
h2 00
2! yn1
5. Explicit Midpoint:
(
n1/2 = yn1 + h2 f (tn1 , yn1 )
y
n1/2 )
yn
= yn1 + hf (tn1/2 , y
6. RK Methods: For 1 i s
(
Yi
yn
Ps
= yn1 + h j=1 aij f (tn1 + cj , h, Yj )
Ps
= yn1 + h i=1 bi f (tn1 + ci h, Yi )
Ki
yn
Ps
= f tn1 + ci h, yn1 + h j=1 aij Kj
Ps
= yn1 + h i=1 bi Ki
(a) Tableau:
c1 a11
c2 a21
..
..
.
.
cs as1
b1
(b) Require: ci =
a12
a22
..
.
as2
b2
Ps
j=1
a1s
a2s
..
.
ass
bs
Ps
j=1 bj
=1
j ynj = h
j=0
k
X
j fnj
j=0
k
X
i yni = h0 f (tn , yn )
i=0
k
k
k
X
X
X
1
1
C0 =
j ,
Ci = (1)i
j i j +
j i1 j
i!
(i
1)!
j=0
j=0
j=0
(d) Char. Poly.s: () =
Pk
j=0
j kj , () =
Pk
j=0
j kj .
(e) 0-Stable: All roots of () satisfy |i | 1, and if |i | = 1 it is simple. Strongly stable if all |i | < 1
(except = 1), weakly stable if 0-stable but not strongly stable.
(f) Stability Region: z =
(ei )
(ei )
8. Predictor Corrector Methods: Use explicit multistep method to predict, then implicit method to correct.
2.b.vi)
BVPs
(
y0 = A(t)y + q(t)
0<t<b
B0 y(0) + Bb y(b) = b
13
Q = B0 + Bb Y (b)
3)
Analysis
1. Pointwise Convergence: > 0, x I, N (, x) : |fn (x) f (x)| < , n > N (x, ).
2. Uniform Convergence: > 0, N () : |fn (x) f (x)| < x and n > N (). ALT: limn ||fn f || = 0.
3.a)
Metric Spaces
`p Spaces
X
p
|xj | ,
||x|| = sup |xj |
||x||p =
0j<
j=0
|xj yj |
j=0
1/p
p
|xj |
j=0
1/q
q
|yj |
j=0
n
n
X
X
pj f (xj )
f
pj xj
j=1
j=1
3.a.ii)
1/p
,
|f (x)| dx
p
A norm is derived from an inner product iff the parallelogram law holds, i.e. ||x + y|| + ||x y|| =
2
2
2 ||x|| + 2 ||y|| .
3.b)
Topology
3.b.i)
Basic Definitions
(b) A B A B
(c) A = A
(d) A B = A B
(e) =
9. G : Countable intersection of open sets.
10. F : Countable union of closed sets.
3.b.ii)
3.c)
Distributions
Topological linear space of functions D. Distributions are complex-valued continuous (with respect to D) linear
functions D0 . I.e. if n in D, then hT, n i hT, i.
3.c.i)
Basics
Z
N
N
X
X
xN +1 1
1 (j)
1 (j)
(0)xj +
(1 u)N (N +1) (xu) du =
(0)xj + xN +1 (x)
j!
N
!
j!
0
j=0
j=0
lim
A B
Z
f (x) dx
or
lim
f (x) dx + lim
A
16
f (x) dx
C
lim
f (x) dx
R
Specific Distributions
"Z
#
Z
Z
(x)
(0)
(x) (0)
1
def
dx = lim
dx +
dx
() = lim
1. Principal Value: P V
0 |x|>
0
x
x
x
|x|> x
|x|>
1
2.
= PV
x i0
3.c.iii)
1
i(x)
x
Function Spaces
(b) Derivatives n
sup
n x[N,N ]
(m)
n (x) (m) (x) 0
k
ALT: T D0 iff N , B(N ), k(N ): |hT, i| B ||||N,k CN
.
7. Norm: ||||N,k =
k
X
sup (j) (t)
j=0 [N,N ]
def
8. Topology: d(1 , 2 ) =
||1 2 ||N,k
k=0
1 + ||1 2 ||N,k
1
2k
10. Derivatives: T (n) , = (1)n T, (n) .
17
3.d)
Measure Theory
3.d.i)
Basic Definitions
Bn
and
n=1
volume of Bn
n=1
Measurable Spaces
Probability Spaces
P (Em ) < ,
m=1
n=1
4. Borel Zero-One Law: En independently often (i.o.) can only have probability 0 or 1.
18
3.e)
3.e.i)
3.f )
3.f.i)
Uniformity
Basic Definitions
1. Uniformly Continuous: f : (M, d) (M 0 , d0 ). If > 0, > 0 : d(x, y) < d0 (f (x), f (y)) < .
If A M is compact and f : M M 0 is continuous, then f is uniformly continuous on A.
2. Uniform Convergence: fn : (M, d) (M 0 , d0 ).
(a) fn converges to f if , x, N (x, ) : n > N (x, ) d0 (f (x), fn (x)) < .
(b) fn converges uniformly to f if , N () : n > N () d0 (f (x), fn (x)) < x.
(c) {fn } is Cauchy if , x, N () : n, m > N (x, ) d0 (fn (x), fm (x)) < .
(d) {fn } is uniformly Cauchy if , N () : n, m > N () d0 (fn (x), fm (x)) < x.
3. fn : M M 0 , M 0 complete. fn converges Uniformly iff it is uniformly Cauchy.
P
4. A M , fn : A R. n=0 fn (x) converges uniformly on A if the sequence of partial sums converges uniformly.
P
5. Weirstrass
M-test: (for uniform convergence). If Mn 0 : |fn (x)| Mn x A, and n=0 Mn < ,
P
then n=0 fn converges uniformly.
19
3.f.ii)
More on Integrals
1. Riemann integrable: Define Mi = max{f (t) : ti t ti+1 }, mi = min{f (t) : ti t ti+1 }, and U (f ; ) =
n
n
X
X
Mi (ti+1 ti ), L(f ; ) =
mi (ti+1 ti ) for some partition . If inf all U (f ; ) = supall L(f ; ) then
i=0
i=0
f is Riemann integrable.
(a) f on [a, b] is R-integrable if it is continuous.
(b) f on [a, b] is R-integrable iff the set of points of discontinuity of f has Lebesgue measure zero.
(c) Riemann integral exists Lebesgue integral exists, and they are equal (not converse).
2. Lebesgue integrable: Suppose f bounded, M, f (t) M and partition range: M 1 = y0 < y1 < <
n
X
yn < yn+1 = M + 1. Let i be the length of {t : yi f (t) < yi+1 }. Then the Lebesgue sum is:
yi i .
i=0
Alternatively, define
(
f (x)
f+ (x) =
0
if f (x) 0
else
f (x) =
Z
Then f (x) = f+ (x) f (x), and
Z
f+ d
f d =
X
n
X
j=1
Z
non-negative simple function, and has integral
d =
X
n
X
j (Aj ).
j=1
3.g)
3.g.i)
Convergence Theorems
Basic Definitions
def
mn
def
inf xm , alternatively the leftmost limit point (or ).
sup xm , alternatively the rightmost limit point (or ).
mn
1
20
Z
|f | d.
X
1
2
|f | d.
X
3.g.ii)
Theorems
Z
fn d.
X
Z X
gn d =
X n=1
X n
X
Z
X
n=1
gn d.
4)
4.a)
4.a.i)
X n
d
x
=
1. Set x = L
x, t = T t, u = C u
, where is dimensionless. Then x = x dx
1
,
L x
t = t ddtt =
1
T
t.
2. Discrete Symmetries: Example: sign invariance, i.e. given solution u to differential equation, u is also a
solution.
3. Continuous Symmetries: Translation invariance, e.g. t t + , x x + .
4. Traveling Wave: Both time and spatially invariant. u(x, t) = u(x ct) = u(z).
5. Scaling Symmetry: u (x, t) = c u(a x, b t), with z(x, t) = z(a x, b t).
6. Buckingham Pi: Given f (x1 , . . . , xn ) = 0 of n physical variables in k physical units, can restate as F (1 , . . . , p ) =
n
Y
a
0, where i =
xj j , and p = n k.
j=0
7. Symmetries Reductions:
(a) Time translation (t t + ) steady solutions (u(x, t) = u(x))
(b) Space translation (x x + ) homogeneous solutions (u(x, t) = u(t))
(c) Rescaling symmetries (u u ) self similar solutions (u(x, t) = tc/a u(x/ta/b ))
4.a.ii)
Phase Plane
21
Dispersion Relations
4.b)
Contour Integration
4.b.i)
Complex Basics
Integration
I
f (z)
dz.
k+1
C (z z0 )
I
X
f (z)
1
n
2. Laurent Series: f (z) =
cn (z z0 ) where cn =
dz.
2i
(z
z0 )n+1
C
n=
k!
2i
I
5. Cauchys Residue Theorem: For a simple closed contour C,
f (z) dz = 2i
C
N
X
j=1
CR
22
4.c)
Fourier Series
n
n
n
X
a0 X
an cos
x + bn sin
x
or
=
x
cn exp
2 n=1
L
L
L
n=
Z
Z
Z L
1 L
1
1 L
f (x) cos(nx) dx
bn =
f (x) sin(nx) dx
cn =
f (x)einx dx
an =
L L
L L
2L L
f (x) =
a0 = 2c0 ;
4.c.i)
an = (cn + cn );
bn = i(cn cn );
Hilbert Space
A complete, normed linear space whose norm comes from an inner product.
1. Inner Product: Satisfies hf, f i 0 with = iff f = 0, hf, gi = hf, gi, hf + g, hi = hf, hi + hg, hi,
hf, gi = hg, f i.
Rb
2. For L2 [a, b], hf, gi = a f (x)g(x) dx.
2
Theorems
1. Riemann-Lebesgue: Given f L1 [a, b], limn |cn | = 0. More specifically, for f C n [a, b] (actually only
c
need f (n) L1 [a, b]), |ck |
n for |k| 1, and constant c.
|k|
2
|ck | .
k=
2
3. Carlsons Theorem: f L Sn (f ) f pointwise almost everywhere, where Sn is the partial sum of the
Fourier series.
4. For f (x) periodic, piecewise smooth on [a, b], and integrable, then in (a, b), limn Sn (f ) = f (x) where f is
continuous, and limn Sn (f ) = 12 (f (x+ ) + f (x )) where f is discontinuous (convergence may break down at
endpoints).
5. For f (x) periodic, continuous, and piecwise smooth on [a, b], then Sn (f ) f uniformly. The convergence is at
the rate ||Sn (f ) f || = O(n1p ), where f (p) L1 [a, b] (or O(np ) for f C p [a, b]).
6. Gibbs Phenomenon: Apporximately 10% error near discontinuities in function.
7. Can integrate term by term to get Fourier series of F (x) (but need to make it periodic), and differentiate term
by term to get Fourier series of f 0 (x) on the condition that f C, f 0 L1 .
4.d)
Distributions
Basics
n
d
1. Null sequence: In C0 (R), {m (x)} such that lim
sup n m (x) = 0 for all n Z+ . In S(R),
m x[K,K] dx
k dn
{m (x)} such that lim sup x
m (x) = 0 for all k, n Z+ .
m xR
dxn
dm
dm
2. Seminorms: Examples are pm,k () = sup xk m , and qm,k () = sup (1 + x2 )k m .
dx
dx
xR
xR
23
5. If f Lloc
1 , then Tf is a distribution in C0 and S.
4.e)
Fourier Transforms
Z
Z
h
i
1
fb(k)eikx dk
F 1 fb(k) = f (x) =
2
f (x)eikx dx
F [f (x)] = fb(k) =
4.e.i)
Common Functions/Distributions
1. Sinc and box: F [sinc(x)] = [1,1] , and F [1,1] = 2 sinc(k)
h
i
2
2
2. Gaussian: F ex /a = a exp ak
4
3. Sin and cos: F [sin(ax)] = i [(k + a) (k a)] and F [cos(ax)] = [(k + a) + (k a)]
4. Delta: F [(x)] = 1, and F [1] = 2(k)
5. Heaviside: F [H(x)] = (k) iP V
6. Principal Value and sign: F P V
4.e.ii)
1
k
1
x
1
k
Other stuff
1. Poisson Sum:
(x 2n) =
n=
1 X ikx
e
2
|a|
hR
x
i
f (t) dt =
(nT ) =
n=
k=
1 X
2k
b
T
T
k=
b
fb(x) = 2f (x)
1 b
ik f (k)
+ fb(0)(k)
D
E
2
fb, gb = 2 hf, gi, so ||f ||2 =
f L2 (R)).
h
i
9. Convolution: F [f ? g] = fb gb, and F 1 f ? g = 2f g.
Z
10. Multi dimensions: F [f (x)] = fb(k) =
f (x)eikx dx;
Rn
cf [] = T b[] = Tf []
11. FT of Distributions: T
b
f
24
h
i
F 1 fb(k) = f (x) =
1
(2)n
Z
Rn
1
2
2
b
f (for
2
fb(k)eikx dk
4.e.iii)
Sampling
Other Transforms
4.f )
1
f
Spectral Theory
4.f.i)
Basics
3. Adjoint: An operator L such that hLu, vi = hu, L vi. If L is bounded and linear, then L exists and is
bounded and linear. Also, L = L, and ||L|| = ||L ||.
4. Normal: L L = LL .
5. Self Adjoint: L = L. formally if they have different domains. All eigenvalues are real. No residual
spectrum.
6. Resolvent: Complete normed space X 6= {0} and linear operator L : D X, D X. The resolvent is
R (L) = (L I)1 .
7. Regular value: such that R (L) exists, is bounded, and is defined on a dense set in X.
8. Point Spectrum: p (L) = { : (L I)1 d.n.e.}. Eigenvalues of L, Lu = u. p (L) p (L ).
9. Continuous Spectrum: c (L) = { : (L I)1 exists and is unbounded}.
10. Residual Spectrum: c (L) = { : (LI)1 exists, possibly bounded, but not defined on a dense set of X}.
11. Rayleigh Quotient: =
4.f.ii)
hLu, ui
.
hu, ui
Sturm-Liouville
1. Lu =
1
((p(x)ux )x + q(x)u(x)) on some domain [a, b]
(x)
25
4.g)
Greens Functions
4.g.i)
Variation of Parameters
u u2
= u1 u02 u2 u01 , for solutions Lu1 = Lu2 = 0 of the S-L operator.
1. Wronskian: W (x) = det 10
u1 u02
2. v10 (x) =
u2 (x)f (x)
;
p(x)W (x)
u1 (x)f (x)
.
p(x)W (x)
v20 (x) =
Z
u2 (t)f (t)
dt + u2 (x)
p(t)W (t)
3. up (x) = u1 v1 + u2 v2 = u1 (x)
0
Z
0
u1 (t)f (t)
dt
p(t)W (t)
Greens Functions
1
p(t)
Z
G(x, t)f (t) dt =
26