Residues and Its Applications

Download as pdf or txt
Download as pdf or txt
You are on page 1of 60

EE202 - EE MATH II Jitkomut Songsiri

12. Residues and Its Applications

• isolated singular points

• residues

• Cauchy’s residue theorem

• applications of residues

12-1
Isolated singular points

z0 is called a singular point of f if

• f fails to be analytic at z0
• but f is analytic at some point in every neighborhood of z0

a singular point z0 is said to be isolated if f is analytic in some punctured disk

0 < |z − z0| < 

centered at z0 (also called a deleted neighborhood of z0)

example: f (z) = 1/(z 2(z 2 + 1)) has the three isolated singular points at

z = 0, z = ±j

Residues and Its Applications 12-2


Non-isolated singular points

1
example: the function has the singular points
sin(π/z)

1
z = 0, z= , (n = ±1, ±2, . . .)
n

• each singular point except z = 0 is isolated


• 0 is nonisolated since every punctured disk of 0 contains other singularities
• for any ε > 0, we can find a positive integer n such that n > 1/ε
• this means z = 1/n always lies in the punctured disk 0 < |z| < ε

Residues and Its Applications 12-3


Residues
assumption: z0 is an isolated singular point of f , e.g.,
there exists a punctured disk 0 < |z − z0| < r0 throughout which f is analytic
consequently, f has a Laurent series representation

X b1 n bn
f (z) = an(z − z0) + + ··· + n
+ ··· , (0 < |z − z0| < r0)
n=0
z − z 0 (z − z 0 )

let C be any positively oriented simple closed contour lying in the disk

0 < |z − z0| < r0

the coefficient bn of the Laurent series is given by


Z
1 f (z)
bn = dz, (n = 1, 2, . . .)
j2π C (z − z0)−n+1

Residues and Its Applications 12-4


the coefficient of 1/(z − z0) in the Laurent expansion is obtained by
Z
f (z)dz = j2πb1
C

b1 is called the residue of f at the isolated singular point z0, denoted by

b1 = Res f (z)
z=z0

this allows us to write

Z
f (z)dz = j2π Res f (z)
C z=z0

which provides a powerful method for evaluating integrals around a contour

Residues and Its Applications 12-5


2
e1/z dz when C is the positive oriented circle |z| = 1
R
example: find C

1/z 2 is analytic everywhere except z = 0; 0 is an isolated singular point

the Laurent series expansion of f is

1/z 2 1 1 1
f (z) = e =1+ 2 + + + ··· (0 < |z| < ∞)
z 2!z 4 3!z 6

the residue of f at z = 0 is zero (b1 = 0), so the integral is zero

remark:
R the analyticity of f within and on C is a sufficient condition for
C
f (z)dz to be zero; however, it is not a necessary condition

Residues and Its Applications 12-6


1 R
example: compute C 3
dz where C is circle |z + 2| = 1
z(z + 2)

f has the isolated singular points at 0 and −2


choose an annulus domain: 0 < |z + 2| < 2
on which f is analytic and contains C

f has a Laurent series on this domain and is given by

1 1 1 1
f (z) = =− · ·
(z + 2 − 2)(z + 2)3 2 1 − (z + 2)/2 (z + 2)3
∞ ∞
1 X (z + 2)n X (z + 2)n−3
=− = − , (0 < |z + 2| < 2)
2(z + 2)3 n=0 2n n=0
2n+1

the residue of f at z = −2 is −1/23 which is obtained when n = 2

therefore, the integral is j2π(−1/23) = −jπ/4 (check with the Cauchy formula)

Residues and Its Applications 12-7


Cauchy’s residue theorem

let C be a positively oriented simple closed contour

Theorem: if f is analytic inside and on C except for a finite number of singular


points z1, z2, . . . , zn inside C, then
Z n
X
f (z)dz = j2π Res f (z)
C z=zk
k=1

Proof.
• since zk ’s are isolated points, we can find small
circles Ck ’s that are mutually disjoint
• f is analytic on a multiply connected domain
• from the Cauchy-Goursat theorem:
R Pn R
C
f (z)dz = k=1 C f (z)dz
k

Residues and Its Applications 12-8


example: use the Cauchy residue theorem to evaluate the integral
Z
3(z + 1)
dz, C is the circle |z| = 2, in counterclockwise
C z(z − 1)(z − 3)

C encloses the two singular points of the integrand, so


Z Z
3(z + 1) h i
I= f (z)dz = dz = j2π Res f (z) + Res f (z)
C C z(z − 1)(z − 3) z=0 z=1

• calculate Resz=0 f (z) via the Laurent series of f in 0 < |z| < 1
• calculate Resz=1 f (z) via the Laurent series of f in 0 < |z − 1| < 1

Residues and Its Applications 12-9


1 3 2
rewrite f (z) = − +
z z−1 z−3

• the Laurent series of f in 0 < |z| < 1

1 3 2 1 2
f (z) = + − = +3(1+z+z +. . .)− (1+(z/3)+(z/3)2+. . .)
2
z 1 − z 3(1 − z/3) z 3

the residue of f at 0 is the coefficient of 1/z, so Resz=0 f (z) = 1

• the Laurent series of f in 0 < |z − 1| < 1

1 3 1
f (z) = − −
1 + z − 1 z − 1 1 − (z − 1)/2
 2 !
3 z−1 z−1
= 1 − (z − 1) + (z − 1)2 + . . . − − 1+ + + ...
z−1 2 2

the residue of f at 1 is the coefficient of 1/(z − 1), so Resz=0 f (z) = −3

Residues and Its Applications 12-10


therefore, I = j2π(1 − 3) = −j4π
alternatively, we can compute the integral from the Cauchy integral formula

!
1 3 2
Z
I= − + dz
C z z−1 z−3
= j2π(1 − 3 + 0) = −j4π

Residues and Its Applications 12-11


Residue at infinity
f is said to have an isolated point at z0 = ∞ if

there exists R > 0 such that f is analytic for R < |z| < ∞

singular points of

C is a positive oriented simple closed contour

Theorem: if f is analytic everywhere except for a finite number of singular


points interior to C, then
Z   
1 1
f (z)dz = j2π Res 2 f
C z=0 z z

(see a proof on section 71, Churchill)

Residues and Its Applications 12-12


z−3 R
example: find I = C dz, C is the circle |z| = 2 (counterclockwise)
z(z − 1)

2
 
I = j2π Res (1/z )f (1/z)
z=0
 
1 − 3z
= j2π Res , j2π Res g(z)
z=0 z(1 − z) z=0

find the residue via the Laurent series of g in 0 < |z| < 1
 
1
write g(z) = − 3 (1 + z + z 2 + · · · ) =⇒ Res g(z) = 1
z z=0

compare the integral with other methods .

• Cauchy integral formula (write the partial fraction of f )


• Cauchy residue theorem (have to find two residues; hence two Laurent series)

Residues and Its Applications 12-13


Principal part

f has an isolated singular point at z0, so f has a Laurent seires


X b1 n b2 bn
f (z) = an(z − z0) + + 2
+ ··· + n
+ ···
n=0
(z − z 0 ) (z − z 0 ) (z − z 0 )

in a punctured disk 0 < |z − z0| < R

the portion of the series that involves negative powers of z − z0

b1 b2 bn
+ + ··· + + ···
(z − z0) (z − z0)2 (z − z0)n

is called the principal part of f

Residues and Its Applications 12-14


Types of isolated singular points

three possible types of the principal part of f

• no principal part

z2 z4
f (z) = cos z = 1 − + + · · · , (0 < |z| < ∞)
2! 4!

• finite number of terms in the principal part

1 1 1 2
f (z) = 2
= 2
− + 1 − z + z + ··· , (0 < |z| < 1)
z (1 + z) z z

• infinite number of terms in the principal part

1 1 1
f (z) = e1/z = 1 + + 2
+ 3
+ ··· , (0 < |z| < ∞)
z 2!z 3!z

Residues and Its Applications 12-15


classify the number of terms in the principal part in a general form

• none: z0 is called a removable singular point



X
f (z) = an(z − z0)n
n=0

• finite (m terms): z0 is called a pole of order m



X b1 b2 bm
f (z) = an(z − z0)n + + + · · · +
n=0
z − z0 (z − z0)2 (z − z0)m

• infinite: z0 is said to be an essential singular point of f



X b1 b2 bn
f (z) = an(z − z0)n + + 2
+ · · · + n
+ ···
n=0
z − z0 (z − z0) (z − z0)

Residues and Its Applications 12-16


examples:

z2 z4
f1(z) = cos z = 1 − + + · · ·
2! 4!
 
3 1
f2(z) = =− + 1 + (z − 2) + (z − 2)3 + · · ·
(z − 1)(z − 2) z−2
1 1 1 2
f3(z) = = − + 1 − z + z + ···
z 2(1 + z) z 2 z
1 1 1
f4(z) = e1/z = 1 + + 2
+ 3
+ ···
z 2!z 3!z
• 0 is a removeable singular point of f1
• 2 is a pole of order 1 (or simple pole) of f2
• 0 is a pole of order 2 (or double pole) of f3
• 0 is an essential singular point of f4

note: for f2, f3 we can determine the pole/order from the denominator of f

Residues and Its Applications 12-17


Characterization of poles

an isolated singular point z0 of a function f is a pole of order m if and only if

φ(z)
f (z) =
(z − z0)m

where φ(z) is analytic and nonzero at z0

Proof. since φ is analytic at z0, it has Taylor series about z = z0


φ(m−1)(z0)(z − z0)m−1 X φ(k)(z0)(z − z0)k
φ(z) = φ(z0) + · · · + +
(m − 1)! k!
k=m

φ(z0) φ(m−1)(z0) X (k)
φ (z0)(z − z0)k−m
f (z) = m
+ ··· + +
(z − z0) (m − 1)!(z − z0) k!
k=m

f has a pole at z0 of order m when φ is nonzero at z0

Residues and Its Applications 12-18


Residue formula

if f has a pole of order m at z0 then

1 dm−1 m
Res f (z) = lim m−1
(z − z 0 ) f (z)
z=z0 (m − 1)! z→z0 dz

Proof. if f has a pole of order m, its Laurent series can be expressed as


X b1 n b2 bm
f (z) = an(z − z0) + + 2
+ ··· + m
n=0
(z − z 0 ) (z − z 0 ) (z − z 0 )

X
(z − z0)mf (z) = an(z − z0)m+n + b1(z − z0)m−1 + b2(z − z0)m−2 + · · · + bm
n=0

to obtain b1, we take the (m − 1)th derivative and take the limit z → z0

Residues and Its Applications 12-19


(z + 1)
example 1: find Resz=0 f (z) and Resz=2 f (z) where f (z) = 2
z (z − 2)
 
d z+1
Res f (z) = lim = −3/4 (0 is a double pole of f )
z=0 z→0 dz z−2
z+1
Res f (z) = lim 2 = 3/4
z=2 z→2 z

z+1
example 2: find Resz=0 g(z) where g(z) =
1 − 2z
g is analytic at 0 (0 is a removable singular point of g), so Resz=0 g(z) = 0

check . apply the results from the above two examples to compute

(z + 1)
Z
2 (z − 2)
dz, C is the circle |z| = 3 (counterclockwise)
C z

by using the Cauchy residue theorem and the formula on page 12-12

Residues and Its Applications 12-20


sometimes the pole order cannot be readily determined

sinh z
example 3: find Resz=0 f (z) where f (z) =
z4

use the Maclaurin series of sinh z

z3 z5
   
1 1 1 z
f (z) = 4 · z + + + · · · = + + + ···
z 3! 5! z 3 3!z 5!

0 is the third-order pole with residue 1/3!

here we determine the residue at z = 0 from its definition (the coeff. of 1/z )

no need to use the residue formula on page 12-19

Residues and Its Applications 12-21


L’Hôpital’s rule for complex functions

let f (z) and g(z) be analytic in a region containing z0 and

• f (z0) = g(z0) = 0

• g 0(z0) 6= 0

then the L’Hôpital’s rule states that

f (z) f 0(z0)
lim = 0
z→z0 g(z) g (z0)

in case f 0(z0) = g 0(z0) = 0, the rule may be extended

example: limz→0 sinz z = 1

Residues and Its Applications 12-22


when the pole order (m) is unknown, we can

• assume m = 1, 2, 3, . . .
• find the corresponding residues until we find the first finite value

1+z
example 4: find Resz=0 f (z) where f (z) =
1 − cos z
• assume m = 1
z(1 + z) 1 + 2z
Res f (z) = lim = 0/0 = lim = 1/0 = ∞ =⇒ (not 1st order)
z=0 z→0 1 − cos z z→0 sin z

• assume m = 2
2
 
d z (1 + z)
Res f (z) = lim = 2 (finite) =⇒ 0 is a double pole
z=0 z→0 dz 1 − cos z

note: use L’Hôpital’s rule to compute the limit

Residues and Its Applications 12-23


Summary
R
many ways to compute a contour integral ( C f (z)dz)

• parametrize the path (feasible when C is easily described)


• use the principle of deformation of paths (if f is analytic in the region
between the two contours)
• use the Cauchy integral formula (typically requires the partial fraction of f )
• use the Cauchy’s residue theorem on page 12-8 (requires the residues at
singular points enclosed by C)
• use the theorem of the residue at infinity on page 12-12 (find one residue at 0)

to find the residue of f at z0

• read from the coeff of 1/(z − z0) in the Laurent series of f


• apply the residue formula on page 12-19

Residues and Its Applications 12-24


Application of the residue theorem

• improper integrals

• improper integrals from Fourier series

• inversion of Laplace transforms

• integrals involving sines and cosines

ingredients: residue theorem, upper bound of contour integral, Jordan inequality

Residues and Its Applications 12-25


Improper integrals

let’s first consider a well-known improper integral


Z ∞
dx
I= 2

−∞ 1 + x

of course, this can be evaluated using the inverse tangent function

we will derive this kind of integral by means of contour integration


some poles of the integrand lie in the upper half plane
let CR be a semicircular contour with radius R → ∞
Z R Z X
f (x)dx + f (z)dz = j2π Res f (z)
−R CR z=zk
k

R
and show that CR
f (z)dz → 0 as R → ∞

Residues and Its Applications 12-26


Theorem: if all of the following assumptions hold

1. f (z) is analytic in the upper half plane except at a finite number of poles

2. none of the poles of f (z) lies on the real axis

M
3. |f (z)| ≤ k when z = Rejθ ; M is a constant and k > 1
R

then the real improper integral can be evaluated by a contour integration, and
Z ∞  
sum of the residues of f (z) at the poles
f (x)dx = j2π
−∞ which lie in the upper half plane

• assumption 2: f is analytic on C1
R
• assumption 3: C f (z)dz → 0 as R → ∞
R

Residues and Its Applications 12-27


Proof. consider a semicircular contour with radius R large enough to include all
the poles of f (z) that lie in the upper half plane

• from the Cauchy’s residue theorem


Z hX i
f (z)dz = j2π Res f (z) at all poles within C1 ∪ CR
C1 ∪CR

(to apply this, f (z) cannot have singular points on C1, i.e., the real axis)
• the integral along the real axis is our desired integral
Z R Z Z
lim f (x)dx + lim f (z)dz = lim f (z)dz
R→∞ −R R→∞ CR R→∞ C1 ∪CR

• hence, it suffices to show that


Z
lim f (z)dz = 0 by using |f (z)| ≤ M/Rk , where k > 1
R→∞ CR

Residues and Its Applications 12-28


Upper bounds for contour integrals

setting: C denotes a contour of length L and f is piecewise continuous on C

. Theorem: if there exists a constant M > 0 such that

|f (z)| ≤ M

for all z on C at which f (z) is defined, then


Z
b
f (z)dz ≤ M L


a

Rb Rb
Proof sketch: need lemma: | a
w(t)dt| ≤ a
|w(t)|dt for complex
R b R
f (z)dt = f (z(t))z (t)dt ≤ b |f (z(t)z 0(t)|dt ≤ b M |z 0(t)|dt ≤ M · L
R 0
R
C a a a

Residues and Its Applications 12-29


(continue the proof of applying residue theorem)

• apply the modulus of the integral and use |f (z)| ≤ M/Rk


Z
M M πR
f (z)dz ≤ k · length of CR =


CR R Rk
R
hence, limR→∞ CR
f (z)dz = 0 if k > 1

remark: an example of f (z) that satisfies all the conditions in page 12-27

p(x)
f (x) = , p and q are polynomials
q(x)

q(x) has no real roots and deg q(x) ≥ deg p(x) + 2

(relative degree of f is greater than or equal to 2)

Residues and Its Applications 12-30


example: show that Z
f (z)dz = 0
CR
as R → ∞ where CR is the arc z = Rejθ , 0 ≤ θ ≤ π

• f (z) = (z + 2)/(z 3 + 1) (relative degree of f is 2)

|z + 2| ≤ |z| + 2 = R + 2, |z 3 + 1| ≥ ||z 3| − 1| = |R3 − 1|

R+2
hence, |f (z)| ≤ R3 −1
and apply the modulus of the integral

1 + R22
Z Z

f (z)dz ≤ R + 2
|f (z)|dz ≤ 3 · πR = π ·

C

C R − 1 R − R12

the upper bound tends to zero as R → ∞

Residues and Its Applications 12-31


• f (z) = 1/(z 2 + 2z + 2)

z 2 + 2z + 2 = (z − (1 + j))(z − (1 − j)) , (z − z0)(z − z¯0)



hence, |z − z0| ≥ ||z| − |z0|| = R − |1 + j| = R − 2 and similarly,


|z − z¯0| ≥ ||z| − |z¯0|| = R − 2

then it follows that

2
√ 2 1
|z + 2z + 2| ≥ (R − 2) ⇒ |f (z)| ≤ √
(R − 2)2
Z Z
1 π
f (z)dz ≤ |f (z)|dz ≤ √ · πR = √
(R − 2)2
(1 − R2 )2

C C

the upper bound tends to zero as R → ∞

Residues and Its Applications 12-32



dx
Z
example: compute I = 2
−∞ 1 + x

1
• define f (z) = 2
and create a contour C = C1 ∪ CR as on page 12-26
1+z
R
• relative degree of f is 2, so C f (z)dz = 0 as R → ∞
R

• f (z) has poles at z = j and z = −j (no poles on the real axis)


• only the pole z = j lies in the upper half plane
• by the residue’s theorem

X I Z R Z
j2π · Res f (z) = f (z)dz = f (x)dx + f (z)dz
z=zk C | −R {z } | CR {z }
=I as R→∞ =0 as R→∞

I = j2π Res f (z) = j2π lim (z − j)f (z) = π


z=j z→j

Residues and Its Applications 12-33


example: compute

x2
Z
I= 2 + a2 )(x2 + b2 )
dx
−∞ (x

z2
• define f (z) = 2 2 2 2
and create C = C1 ∪ CR as on page 12-26
(z + a )(z + b )
R
• relative degree of f is 2, so C f (z)dz = 0 as R → ∞
R

• f (z) has poles at z = ±ja and z = ±jb (no poles on the real axis)
• only the poles z = ja and z = jb lie in the upper half plane
• by the residue’s theorem
X I Z R Z
j2π · Res f (z) = f (z)dz = f (x)dx + f (z)dz
z=zk C | −R {z } | CR {z }
=I as R→∞ =0 as R→∞
   
a b π
I = j2π Res f (z) + Res f (z) = j2π + =
z=ja z=jb j2(a2 − b2) j2(b2 − a2) a+b

Residues and Its Applications 12-34


Applications of residue theorem

• improper integrals

• improper integrals from Fourier series

• inversion of Laplace transforms

• integrals involving sines and cosines

Residues and Its Applications 12-35


Improper integrals from Fourier analysis

we can use residue theory to evaluate improper integrals of the form


Z ∞ Z ∞ Z ∞
f (x) sin mx dx, f (x) cos mx dx, or ejmxf (x) dx
−∞ −∞ −∞

we note that ejmz is analytic everywhere, moreover

|ejmz | = ejm(x+jy) = e−my < 1 for all y in the upper half plane

therefore, if |f (z)| ≤ M/Rk with k > 1, then so is |ejmz f (z)|

hence, if f (z) satisfies the conditions in page 12-27 then


Z ∞  jmz

sum of the residues of e f (z) at the poles
ejmxf (x)dx = j2π
−∞ which lie in the upper half plane

Residues and Its Applications 12-36


denote
jmz
 
sum of the residues of e f (z) at the poles
S=
which lie in the upper half plane

and note that S can be complex

by comparing the real and imaginary part of the integral


Z ∞ Z ∞
jmx
e f (x)dx = (cos mx + j sin mx)f (x)dx = j2πS
−∞ −∞

we have
Z ∞
cos mxf (x) dx = Re(j2πS) = −2π · Im S
−∞
Z ∞
sin mxf (x) dx = Im(j2πS) = 2π · Re S
−∞

Residues and Its Applications 12-37



cos mx dx
Z
example: compute I =
−∞ 1 + x2

ejmz
• define f (z) = 1+z 2
and create C = C1 ∪ CR as on page 12-26
R
• relative degree of f is 2, so CR
f (z)dz = 0 as R → ∞

• f has poles at z = j and z = −j (no poles on the real axis)

• the pole z = j lies in the upper half plane

• by residue’s theorem

X I Z R Z
j2π · Res f (z) = f (z)dz = f (x)dx + f (z)dz
z=zk C | −R {z } | CR {z }
=I as R→∞ =0 as R→∞

Residues and Its Applications 12-38


• therefore,


ejmx ejmz
Z
2
dx = j2π Res
−∞ 1 + x z=j 1 + z 2

(z − j)ejmz −m
= j2π lim = πe
z→j 1 + z2

• our desired integral can be obtained by


Z ∞
cos mx dx
= Re(πe−m) = πe−m,
−∞ 1 + x2
Z ∞
sin mx dx
2
= Im(πe−m) = 0
−∞ 1 + x

Residues and Its Applications 12-39


Summary of improper integrals

the examples of f we have seen so far are in the form of

p(x)
f (x) =
q(x)

where p, q are polynomials and deg p(x) ≥ deg q(x) + 2

Residues and Its Applications 12-40


the assumption on the degrees of p, q is sufficient to guarantee that
Z
f (z)ejaz dz = 0 (a > 0)
CR

as R → ∞ where CR is the arc z = Rejθ , 0 ≤ θ ≤ π

we can relax this assumption to consider function f such as

z 1
, (relative degree is 1)
z 2 + 2z + 2 z+1

and obtain the same result by making use of Jordan’s inequality

Residues and Its Applications 12-41


Jordan inequality

for R > 0, Z π
−R sin θ π
e dθ <
0 R
Proof.

π
sin θ ≥ 2θ/π, 0≤θ≤
2
π
e−R sin θ ≤ e−2Rθ/π , R > 0, 0 ≤ θ ≤
2
Z π/2
π
e−R sin θ dθ ≤
0 2R

the last line is another form of the Jordan inequality

because the graph of y = sin θ is symmetric about the line θ = π/2

Residues and Its Applications 12-42


z
show that C f (z)ejaz dz = 0 for a > 0 as
R
example: let f (z) = 2
z + 2z + 2 R
R→∞

• first note that |ejaz | = |eja(x+jy)| = |ejax · e−ay | = e−ay < 1 (since a > 0)


• similar to page 12-32, we see that |f (z)| ≤ R/(R − 2)2 , MR and
Z Z
jaz
R π
f (z)e dz ≤ √ · πR = √
2
CR (R − 2) (1 − R2 )2

CR

which does not tend to zero as R → ∞

• however, for z that lies on CR, i.e., z = Rejθ

jaRejθ
f (z)e jaz
= f (z)e = f (z)ejaR(cos θ+j sin θ) = f (z)e−aR sin θ · ejaR cos θ

Residues and Its Applications 12-43


• if we find an upper bound of the integral, and use Jordan’s inequality:
Z Z π
jaz −aR sin θ jaR cos θ jθ


f (z)e dz =
f (z)e ·e jRe dθ
CR 0
Z π
−aR sin θ jaR cos θ jθ

≤ f (z)e ·e jRe dθ
0
Z π
= RMR e−aR sin θ dθ
0
πMR
<
a
the final term approach 0 as R → ∞ because MR → 0

conclusion: then we can apply the residue’s theorem to integrals like


Z ∞
x cos(ax)
2 + 2x + 2
dx
−∞ x

Residues and Its Applications 12-44


Applications of residue theorem

• improper integrals

• improper integrals from Fourier series

• inversion of Laplace transforms

• integrals involving sines and cosines

Residues and Its Applications 12-45


Inversion of Laplace transforms
recall the definitions:
Z ∞
F (s) , L[f (t)] , f (t)e−stdt
0
a+j∞
1
Z
−1
f (t) = L [F (s)] = F (s)estds
j2π a−j∞

Theorem: suppose F (s) is analytic everywhere except at the poles

p1, p2, . . . , pn,

all of which lie to the left of the vertical line Re(s) = a (a convergence factor)
if |F (s)| ≤ MR and MR → 0 as s → ∞ through the half plane Re(s) ≤ a then
n
X
L−1[F (s)] = Res F (s)est
s=pi
i=1

Residues and Its Applications 12-46


Proof sketch.

parametrize C1 and C2 by

C1 = {z | z = a + jy, −R ≤ y ≤ R }
 
π 3π
C2 = z | z = a + Rejθ , ≤θ≤
2 2

1. create a huge semicircle that is large enough to contain all the poles of F (s)
2. apply the Cauchy’s residue theorem to conclude that
Z n
X Z
estF (s)ds = j2π Res [estF (s)] − estF (s)ds
C1 s=pk C2
k=1

3. prove that the integral along C2 is zero when the circle radius goes to ∞

Residues and Its Applications 12-47


choose a and R: choose the center and radius of the circle

• a > 0 is so large that all the poles of F (s) lie to the left of C1

a> max Re(pk )


k=1,2,...,n

• R > 0 is large enough so that all poles of F (s) are enclosed by the semicircle
if the maximum modulus of p1, p2, . . . , pn is R0 then

∀k, |pk − a| ≤ |pk | + a ≤ R0 + a =⇒ pick R > R0 + a

C1 = {z | z = a + jy, −R ≤ y ≤ R }
 
π 3π
C2 = z | z = a + Rejθ , ≤θ≤
2 2

Residues and Its Applications 12-48


integral along C2 is zero

C1 = {z | z = a + jy, −R ≤ y ≤ R }
 
π 3π
C2 = z | z = a + Rejθ , ≤θ≤
2 2

• for s = a + Rejθ and ds = jRejθ dθ, the integral becomes

Z Z 3π/2

estF (s)ds = eat · eRt cos θ+jRt sin θ F (a + Rejθ )Rjejθ dθ



C2 π/2

Residues and Its Applications 12-49


• apply the modolus of the integral
Z Z 3π/2
st at Rt cos θ jRt sin θ jθ jθ


e F (s)ds ≤ e e ·e F (a + Re )Rje dθ
C2 π/2

• since |F (s)| ≤ MR for s that lies on C2


Z Z 3π/2
estF (s)ds ≤ MRReat eRt cos θ dθ



C2 π/2

• make change of variable φ = θ − π/2 and apply the Jordan inequality


Z Z π at
st
at −Rt sin φ πMR e
e F (s)ds ≤ MRRe e dφ <

C2 |0 {z } t
<π/Rt

the last term approaches zero as R → ∞ because MR → 0 (by assumption)

Residues and Its Applications 12-50


−1 s
example: find L [F (s)] where F (s) = 2 and c > 0
(s + c2)2

π 3π
 jθ

C2 = z | z = a + Re , 2 ≤θ≤ 2
poles of F (s) are s = ±jc so we choose a > 0
the semicircle must enclose all the pole
so we have R > a + c

first we verifty that |F (s)| ≤ MR and MR → 0 as s → ∞ for s on C2


we note that |s| = |a + Rejθ | ≤ a + R and |s| ≥ |a − R| = R − a
since |s2 + c2| ≥ ||s|2 − c2| ≥ (R − a)2 − c2 > 0, then

|s| (R + a)
|F (s)| = 2 2 2
≤ 2 2 2
, MR → 0 as R → ∞
|s + c | [(R − a) − c ]

Residues and Its Applications 12-51


therefore, we can apply the theorem on page 12-46

−1
X
st sest sest
L [F (s)] = Res [e F (s)] = Res 2 + Res 2
s=sk s=jc (s + c2 )2 s=−jc (s + c2 )2

poles of F (s) are s = ±jc (double poles)


st
   st st

d se e (1 + ts) 2se
Res estF (s) = lim = −
s=jc s→jc ds (s + jc)2 (s + jc)2 (s + jc)3 s=jc
tejct
=
j4c
st st st
   
st d se e (1 + ts) 2se
Res e F (s) = lim = −
s=−jc s→−jc ds (s − jc)2 (s − jc)2 (s − jc)3 s=−jc

te−jct
=−
j4c

−1 t jct −jct
t sin ct
hence L [F (s)] = 4jc (e −e )=
2c
Residues and Its Applications 12-52
−1 1
example: find L [F (s)] where F (s) =
(s + a)2 + b2

F (s) has poles at s = −a ± jb (simple poles)

L−1[F (s)] = Res estF (s) + Res estF (s)


s=−a+jb s=−a−jb

(provided that |F (s)| ≤ MR and MR → 0 as s → ∞ on C2 ... please check .)

est e(−a+jb)t
Res = lim =
s=−a+jb s=−a+jb s + a + jb j2b
est e(−a−jb)t
Res = lim =
s=−a−jb s=−a−jb s + a − jb −j2b

−at jbt −jbt −at


e (e − e ) e sin(bt)
hence, L−1[F (s)] = =
2jb b

Residues and Its Applications 12-53


Applications of residue theorem

• improper integrals

• improper integrals from Fourier series

• inversion of Laplace transforms

• integrals involving sines and cosines

Residues and Its Applications 12-54


Definite integrals involving sines and cosines

we consider a problem of evaluating definite integrals of the form


Z 2π
F (sin θ, cos θ)dθ
0

since θ varies from 0 to 2π, we can let θ be an argument of a point z

z = ejθ (0 ≤ θ ≤ 2π)

this describe a positively oriented circle C centered at the origin

make the substitutions

z − z −1 z + z −1 dz
sin θ = , cos θ = , dθ =
j2 2 jz

Residues and Its Applications 12-55


this will transform the integral into the contour integral
−1 −1
 
z−z
Z
z+z dz
F ,
C j2 2 jz

• the integrand becomes a function of z

• if the integrand reduces to a rational function of z, we can apply the Cauchy’s


residue theorem

example:
Z 2π Z Z Z
dθ 1dz dz
= (z−z −1 ) jz
= 2 + j5z − 2
, g(z)dz
0 5 + 4 sin θ C 5+4 C 2z C
2j
Z  
dz
= = j2π Res g(z) = 2π/3
C 2(z + 2j)(z + j/2) z=−j/2

where C is the positively oriented circle |z| = 1

Residues and Its Applications 12-56


the above idea can be summarized in the following theorem

Theorem: if F (cos θ, sin θ) is a rational function of cos θ and sin θ which is


finite on the closed interval 0 ≤ θ ≤ 2π, and if f is the function obtained from
F (·, ·) by the substitutions

z + z −1 z − z −1
cos θ = , sin θ =
2 j2

then !
Z 2π X f (z)
F (cos θ, sin θ) dθ = j2π Res
C z=zk jz
k

where the summation takes over all zk ’s that lie within the circle |z| = 1

Residues and Its Applications 12-57



cos 2θ
Z
example: compute I = 2
dθ, −1 < a < 1
0 1 − 2a cos θ + a

make change of variables

ej2θ + e−j2θ z 2 + z −2 z 4 + 1
• cos 2θ = = =
2 2 2z 2
2 2
2 −1 2
az − (a + 1)z + a
• 1 − 2a cos θ + a = 1 − 2a(z + z )/2 + a = −
z
R 2π R f (z) R
we have 0
F (θ)dθ = C jz
dz , C
g(z)dz where

(z 4 + 1)z (z 4 + 1)
g(z) = − =
jz · 2z (az − (a + 1)z + a) j2z 2(1 − az)(z − a)
2 2 2

we see that only the poles z = 0 and z = a lie inside the unit circle C

Residues and Its Applications 12-58


therefore, the integral becomes
Z  
I= g(z)dz = j2π Res g(z) + Res g(z)
C z=0 z=a

• note that z = 0 is a double pole of g(z), so

d 2 1 a2 + 1
Res g(z) = lim (z g(z)) = − ·
z=0 z=0 dz j2 a2

1 a4 + 1
• Res g(z) = lim (z − a)g(z) = · 2
z=a z=a j2 a (1 − a2)

2πa2
hence, I =
1 − a2

Residues and Its Applications 12-59


References

Chapter 6-7 in
J. W. Brown and R. V. Churchill, Complex Variables and Applications, 8th
edition, McGraw-Hill, 2009

Chapter 7 in
T. W. Gamelin, Complex Analysis, Springer, 2001

Chapter 22 in
M. Dejnakarin, Mathematics for Electrical Engineering, CU Press, 2006

Residues and Its Applications 12-60

You might also like