Differential Equations & Integral Transforms
Differential Equations & Integral Transforms
B
INTEGRAL TRANSFORMS
C hapters
4. Fourier Transforms
5. Finite Fourier Transforms
7. Fourier Series
T-3
1
T he L aplace T ransform
is called the integral transform of the function F (t) and is denoted by T {F (t)}. The
function K ( p, t) is called the kernel of the transformation.
Remark: Some authors use the letter ‘s’ in place of ‘p’.
T-4
The function f ( p) defined by the integral (1) is called the Laplace transform of the
function F (t) and is also denoted by
L {F (t)} or F ( p).
Thus Laplace transform is a function of a new variable (or parameter) p given by (1).
Note: The Laplace transform of F (t) is said to exist if the integral (1) converges for some
values of p, otherwise it does not exist.
∞ − pt
∴ L { a1 F1 (t) + a2 F2 (t)} = ∫0 e { a1 F1(t) + a2 F2 (t)} dt
∞ ∞
= a1 ∫ e − pt F1 (t) dt + a2 ∫0 e − pt F2 (t) dt
0
t0 ∞
= ∫0 e − pt F (t) dt + ∫ t0 e − pt F (t) dt …(1)
t0
The integral ∫ e − pt F (t) dt exists since F (t) is piecewise continuous on every finite
0
interval 0 ≤ t ≤ t0 .
∞ ∞ ∞
Now ∫ t0 e − pt F (t) dt ≤ ∫ | e − pt F (t)| dt ≤ ∫ t0 e − pt Me at dt ,
t0
since | F (t)| ≤ Me at
∞ Me − ( p − a) t0
= ∫ t0 e −(p − a) t M dt = , p > a.
p− a
∞ Me −(p − a)t0
∴ ∫ t0 e − pt F (t) dt ≤ , p > a.
p− a
Me − ( p − a) t0
But can be made as small as we please by taking t0 sufficiently large.
p− a
Thus from (1), we conclude that L {F (t)} exists for all p > a.
Note 1: The above theorem of existence of Laplace transform can also be stated
as :
T-6
“If F (t) is a function which is piece-wise continuous on every finite interval in the range
t ≥ 0 and is of exponential order a as t → ∞, the Laplace transform of F (t) exists for all
p > a”.
Or
“If F (t) is a function of class A, the Laplace transform of F (t) exists for p > a”.
Note 2: Conditions in the theorem are sufficient but not necessary for the existence of
Laplace transform. If these conditions are satisfied, the Laplace transform must exist. If
these conditions are not satisfied, the Laplace transform may or may not exist.
We can show this by the following example.
Example: Consider the function F (t) = 1 / √ t.
∞ 1
= ∫0 e − pt ⋅ dt, which converges for p > 0
√t
2 ∞ 2 dt 2
= ∫0 e − x dx, putting √ ( pt) = x so that = dx
√p √t √ p
2 √π ∞ 2 √π
= ⋅ , since ∫0 e − x dx =
√p 2 2
= √ (π / p), p > 0.
Thus L {1/ √ t} exists for p > 0 even if 1/ √ t is not piecewise continuous in the
range t ≥ 0.
∞
∞ e − pt 1
∴ L {1} = ∫0 e − pt ⋅ 1 dt = − = , p> 0 .
p 0 p
Here the condition p > 0 is necessary, since the integral is convergent for p > 0 and
divergent for p ≤ 0.
T-7
(ii) Laplace transform of the function F (t) = t n, n is any real number greater than –1.
Solution: We have,
∞
L {F (t)} = ∫0 e − pt F (t) dt.
∞ ∞
∴ L { t n} = ∫0 e − pt t ndt = ∫0 e − pt t(n + 1) − 1dt …(1)
Now from Eulerian integral of the second kind known as ‘Gamma function’, we have
∞ − ax m −1 Γ(m)
∫0 e x dx = am , if a > 0 and m > 0.
Γ (n + 1)
∴ from (1), L { t n} = , if p > 0 and n + 1 > 0 i.e., n > − 1.
p n+1
Thus if n is any real number greater than –1, we have
Γ (n + 1)
L { t n} = , p > 0.
p n+1
Here the condition p > 0 is necessary for the convergence of the integral (1).
(iii) Laplace transform of the function F (t) = t n, n is a positive integer.
(Gorakhpur 2008)
Solution: We have
∞
L {F (t)} = ∫0 e − pt F (t) dt
∞ ∞
∴ L { t n} = ∫0 e − pt t n dt = ∫0 e − pt t(n + 1) −1 dt
Γ (n + 1)
= , if p > 0
p n+1 ∞ − ax m −1 Γ(m)
e x dx = m ,
∫ 0
∵
a
if a > 0 and m > 0. Here n + 1 > 0, n being a positive integer.
n!
= n +1 , p > 0.
p [ ∵ Γ (n + 1) = n !, n being a positive integer]
n!
Thus if n is any positive integer, we have L { t n} = n + 1 , p > 0 .
p
Here the condition p > 0 is necessary for the convergence of the integral defining the
Laplace transform of t n.
(iv) Laplace transform of the function F (t) = e at . (Rohilkhand 2009)
Solution: Here
∞ ∞
L { e at} = ∫0 e − pt ⋅ e at dt = ∫0 e − (p − a) t dt
T-8
∞
e −(p − a) t
= − , p≠ a
p− a
0
1
= , p > a.
p− a
Here the condition p > a is necessary, since the integral is convergent for p > a and
divergent for p ≤ a.
(v) Laplace transform of the function F (t) = sin at.
(Avadh 2010, 11)
∞
Solution: L {sin at} = ∫0 e − pt sin at dt
∞
e − pt
= 2 2
(− p sin at − a cos at)
p + a 0
ax e ax
∵ ∫ e sin bx dx = 2 2
(a sin bx − b cos bx)
a +b
a
= , p > 0.
p + a2
2
Here the condition p > 0 is necessary for the convergence of the integral (1).
(vi) Laplace transform of the function F (t) = cos at. (Rohilkhand 2006, 07, 10)
∞
Solution: L {cos at } = ∫0 e− pt cos at dt
∞
e − pt
= 2 2
(− p cos at + a sin at)
p + a 0
ax e ax
∵ ∫e cos bx dx = 2
a +b
(a cos bx + b sin bx)
2
p
= , p > 0.
p + a2
2
F (t) L {F (t)}
1
1. 1 , p> 0
p
tn n!
2. , p> 0
(n is a positive integer) p n +1
ta Γ (a + 1)
3. , p> 0
(a > − 1) p a +1
1
4. e at , p> a
p− a
a
5. sin at , p> 0
p + a2
2
p
6. cos at , p> 0
p + a2
2
a
7. sinh at , p > | a|
p2 − a2
p
8. cosh at , p > | a|
p − a2
2
If we know the Laplace transforms in the above table then nearly all the transforms can
be obtained by using the general theorems which we shall consider later on.
T-10
Solution: We have
L { (t2 + 1)2} = L { t 4 + 2 t2 + 1} = L { t 4} + 2 L {t2} + L {1}
2 4
4! 2! 1 24 + 4 p + p
= +2⋅ + = , p > 0.
p5 p3 p p5
Example 2: Find the L.T. of the function
F (t) = (sin t − cos t)2 . (Gorakhpur 2006, 10)
2 2 2
Solution: We have L { (sin t − cos t) } = L {sin t + cos t − 2 sin t cos t}
1 2
= L {1} − L {sin 2 t} = − , p> 0
p p2 + 22
p2 − 2 p + 4
= , p > 0.
p ( p2 + 4)
Solution: We have
1 1
L {sin2 at} = L { (1 − cos 2 at)} = [ L {1} − L {cos 2 at}]
2 2
1 1 p
= − 2 2 , p> 0
2 p p + (2 a)
2 a2
= , p > 0.
p ( p2 + 4 a2 )
Solution: We have
L {3 t 4 − 2 t 3 + 4 e −3 t − 2 sin 5 t + 3 cos 2 t}
= 3 L { t 4} − 2 L { t 3} + 4 L { e −3 t} − 2 L{sin 5 t} + 3 L {cos 2 t}
T-11
4! 3! 1 5 p
=3⋅ 5
−2⋅ 4
+4⋅ −2⋅ 2 2
+3⋅ 2 ,
p p p+3 p +5 p + 22
p > 0 and p > − 3 i. e. p > 0
72 12 4 10 3p
= − + − + , p > 0.
p5 p4 p + 3 p2 + 25 p2 + 4
Solution: We have
∞ 1 ∞
L { F (t)} = ∫0 F (t) e − pt dt = ∫ 0 ⋅ e − pt dt + ∫1 (t − 1)2 e − pt dt
0
∞ ∞ 2
= ∫1 (t − 1)2 e − pt dt = ∫0 x e − p ( x +1) dx, putting t − 1 = x ,
Solution: We have
(√ t)3 (√ t)5 (√ t)7
L {sin √ t} = L √ t − + − + ...
3! 5! 7!
1 /2 t3 /2 t5 /2 t7 /2
= L t − + − + ....
3! 5! 7!
1 /2 1 3 /2 1 5 /2 1
= L {t } − L {t } + L {t } − L { t7 /2} + ...
3! 5! 7!
3 5 7 9
Γ( ) Γ( ) Γ( ) Γ( )
2 1 2 1 2 1
= 3 /2 − ⋅ + ⋅ − ⋅ 2 + ...
p 3 ! p5 /2 5 ! p7 /2 7 ! p9 /2
T-12
1 3 1 5 3 1
√π ⋅ √π ⋅ ⋅ √π
1 1
= 2 3 /2 − ⋅2 2 + ⋅ 2 2 2 − ...
p 1⋅ 2 ⋅ 3 p5 /2 1⋅ 2 ⋅ 3 ⋅ 4 ⋅ 5 p7 /2
2 3
√π 1 1 1 1 1 1
= 3 /2
1 − + − + ...
2p 1! 4 p 2 ! 4 p 3 ! 4 p
√π
= 3 / 2 ⋅ e −1 / 4 p .
2p
1 1
= lim = lim , where 0 < m < 1
t→ ∞ t − n e at t→ ∞ t m e at
= 0 , if a > 0 .
n
∴ F (t) = t is of exponential order.
Since F (t) = t n is not sectionally continuous over every finite interval in the range t ≥ 0,
hence it is not a function of class A. But t n is integrable from 0 to any positive number
t0 .
∞ − pt ∞ − pt n ∞ − pt (n+1)−1
Now L { F (t)} = ∫0 e F (t) dt = ∫0 e t dt = ∫ e t dt
0
Γ(n + 1)
= , if p > 0 and n + 1 > 0 i. e., n > − 1.
p(n + 1)
Hence the Laplace transform of t n, − 1 < n < 0, exists, although it is not a function of
class A.
Comprehensive Exercise 1
et , 0 < t ≤ 1
3. (i) Find L {F (t)}, if F (t) =
0, t >1
0 , 0 < t < 2
(ii) Find L {F (t)}, where F (t) =
4, t>2
4. (i) Find the Laplace transform of the following function :
x / a, 0 < x < a
f ( x) =
1, x > a.
0 , 0 < t < 1
(ii) Find L {F (t)}, where F (t) = t, 1 < t < 2 (Lucknow 2010)
0 , t > 2.
5. (i) Find Laplace Transform of the function F (t), where
sin t, 0 < t < π
F (t) = .
0, t> π (Lucknow 2007)
e t, 0 < t < 5
(ii) Find L {F (t)}, if F (t) =
3, t > 5.
1, 0 < t < 2
6. (i) Find L {F (t)}, if F (t) =
t, t > 2.
t, 0 < t < 4
(ii) Find L { F (t)} if F (t) =
5, t > 4.
1 1
(ii) Prove that L = ⋅
√ (πt) √ p (Avadh 2014)
cos √ t π −1 / 4 p
8. Show that L = e .
√t p (Purvanchal 2007; Kashi 14)
2
9. Show that t is of exponential order 3.
A nswers 1
1
1. (i) 12 / p4 − 6 / p2 + 8 / p, p > 0 (ii) , p> 0
( p2 + 4)
p2 − 8 12 − 5 p
(iii) 2
, p> 4 (iv) , p> 0
p ( p − 16) p2 + 4
3 p − 20
2. (i) , p> 5
p2 − 25
T-14
16 p − 4 5p 42 + 2 p3 15
(ii) + + + , p> 2
p2 − 4 p2 + 1 p4 p2 + 9
p+9
(iii) , p> 3
p2 − 9
1
(iv) , p > 0 if a < 0 and p > a if a > 0
p ( p − a)
1 4 −2 p
3. (i) [1 − e − ( p − 1)], p ≠ 1 (ii) e , p> 0
( p − 1) p
1
4. (i) 2
(1 − e − ap ), p > 0
ap
2 1 1 1
(ii) − + 2 e −2 p + + 2 e − p , p ≠ 0 .
p p p p
e− p π + 1 1 − e −5( p −1) 3 −5 p
5. (i) 2
(ii) + e , p> 0
p +1 p−1 p
1 1 1 + ( p − 1) e −4 p
6. (i) [1 + e −2 p ] + 2 e −2 p , p > 0 (ii) , p> 0
p p p2
Obviously if the integral (1) converges when p > α, the integral (2) converges when
p − a > α i. e., p > α + a.
putting t − a = x so that dt = dx
∞ ∞
= e − pa ∫0 e − px F ( x) dx = e − pa ∫0 e − pt F (t) dt,
1 ∞
= ∫0 e − p ( x / a) F ( x) dx, putting at = x, so that dt = (1 / a) dx
a
1 ∞ ∵ b b
= ∫0 e −( p / a)t F (t) dt ∫a f ( x) dx = ∫a f (t) dt
a
1 p ∞ − pt
= f , since f ( p) = ∫0 e F (t) dt.
a a
1
Solution: We have L {sin2 t} = L { (1 − cos 2 t)}
2
1 1 p 2
= − 2 2 = 2
= f ( p) , say.
2 p p + 2 p ( p + 4)
p2 − p + 1
Example 12: Given L {F (t)} = ,
(2 p + 1)2 ( p − 1)
p
We have L { φ (t)} = L {cos t} = 2
= f ( p), say.
p +1
2 π /3 ∞ 2π
= ∫0 e − pt ⋅ 0 dt + ∫ 2 π /3 e − pt ⋅ cos t − dt
3
∞ 2π
= ∫ 2 π /3 e − pt ⋅ cos t − dt
3
∞
= ∫0 e − p [ x + (2 π /3)] cos x dx, putting t − 2 π / 3 = x and dt = dx
∞
= e − p (2 π /3) ∫0 e − px cos x dx
∞
= e − p (2 π /3) ∫0 e − pt cos t dt
Comprehensive Exercise 2
e −3 /( p + 1)
10. If L {F (t)} = (1 / p) e −1 / p , prove that L { e − t F (3 t)} = .
p+1
T-18
A nswers 2
1. (i) 6 /( p + 3)4 . (ii) 1 / ( p + a)n
(iii) n !/ ( p − a)n +1.
2. (i) 4 /( p2 − 6 p + 25). (ii) ( p − 3) /( p2 − 6 p + 34).
( p2 − 2 p + 3) p+4
3. (i) 2
⋅ (ii) 2
⋅
( p − 1) ( p − 2 p + 5) p + 8 p + 12
6 5 ( p + 1) 3 p − 24
4. (i) − ⋅ (ii) ⋅
p2 + 2 p + 5 p2 + 2 p − 3 p2 + 4 p + 40
9 p2 − 12 p + 5 1
5. 3
⋅ 6. [ f ( p − iω) + f ( p + iω)].
( p − 1) 2
3 p
7. (i) 2
. (ii) 2
, p > 0.
p −9 p + 25
e − ap e − π p /3
8. , p > 1. 9. , p > 0.
p−1 p2 + 1
∞ ∞
= [e − pt F (t)] 0 + p ∫ e − pt F (t) dt [Integrating by parts]
0
lim − pt
= e F (t) − F (0 ) + pL {F (t)}. …(2)
t→ ∞
Now | F (t)| ≤ Me at for all t ≥ 0 and for some constants a and M.
Case II: In case F ′ (t) is merely piece-wise continuous, the integtral (1) may be
broken as the sum of integrals in different ranges from 0 to ∞ such that F ′ (t) is
continuous in each of such parts.
Then proceeding as in case I, we shall have
L {F ′ (t)} = pL {F (t)} − F (0 ).
Note 1: If F (t) fails to be continuous at t = 0 but
lim
F (t) = F (0 + 0 ) exists, [F (0 + 0 ) is not equal to F (0 ),
t→0
which may or may not exist]
then L {F ′ (t)} = pL {F (t)} − F (0 + 0 ).
Note 2: If F (t) fails to be continuous at t = a, then
L { F ′ (t)} = pL {F (t)} − F (0 ) − e − ap [ F (a + 0 ) − F (a − 0 )]
where F (a + 0 ) and F (a − 0 ) are the limits of F at t = a, as t approaches a from the right
and from the left respectively.
The quantity F (a + 0 ) − F (a − 0 ) is called the jump at the discontinuity t = a.
∞ a ∞
Proof: L {F ′ (t)} = ∫0 e − pt F ′ (t) dt = ∫0 e − pt F ′ (t) dt + ∫a e − pt F ′ (t) dt
a a ∞
= [e − pt F (t)] 0 + p ∫ e − pt F (t) dt + [e − pt F (t)] a
0
∞
+ p∫ e − pt F (t) dt
a
lim − pt
= e − ap F (a − 0 ) − F (0 ) + e F (t)
t→ ∞
∞
− e − ap F (a + 0 ) + p ∫ e − pt F (t) dt
0
= pL {F (t)} − F (0 ) − e − ap [ F (a + 0 ) − F (a − 0 )].
lim − pt
∵ t → ∞ e F (t) = 0 , as shows in the theorem
Note 3: For more than one discontinuity of the function F (t), appropriate
modification can be made.
= p3 L {F (t)} − p2 F (0 ) − pF ′ (0 ) − F ′ ′ (0 ).
= pL { F (t)} − F (0 ). …(1)
Taking limit as p → 0 in (1), we have
lim ∞ − pt lim
e F ′ (t) dt = pL {F (t)} − F (0 )
p→ 0 ∫0 p→ 0
∞
or ∫0 F ′ (t) dt = lim pL { F (t) − F (0 )}
p→0
lim
or [ F(t)] 0∞ = p → 0 p L {F (t)} − F (0 )
lim lim
or F (t) − F (0 ) = pL {F (t)} − F (0 )
t→ ∞ p→ 0
lim lim
or F (t) = pL {F (t)}.
t→ ∞ p→ 0
lim
Note: If F (t) fails to be continuous at t = 0, but F (t) exists, the result still holds
t→0
by using Note (1), 1.13 theorem.
M at
∴ | G (t)| ≤ (e − 1), a > 0 . …(2)
a
Further G ′ (t) = F (t), except for points at which F (t) is discontinuous.
Therefore G ′ (t) is piece-wise continuous on each finite interval.
Hence from the theorem of 1.13, we have
L { G ′ (t)} = pL { G (t)} − G (0 ) = pL {G (t)} [Since G (0 ) = 0 from (2)]
1
∴ L { G (t)} = L { G ′ (t)}
p
t 1
L ∫ F ( x) dx = L { F (t)}
0 p
1.18 Multiplication By t
Theorem: If F (t) is a function of class A and if L {F (t)} = f ( p), then
L {t F (t)} = − f ′ ( p).
(Rohilkhand 2009, Gorakhpur 05, 09)
∞ − pt
Proof: We have f ( p) = L {F (t)} = ∫0 e F (t) dt.
d ∞ − pt
f ′ ( p) = e F (t) dt
dp ∫ 0
∴
∞ ∂ − pt
= ∫0 {e F (t)} dt,
∂p
by Leibnitz’s rule for differentiating under the sign of integral
∞
=− ∫0 te − pt F (t) dt
∞
=− ∫0 e − pt { t F (t)} dt = − L { t F (t)}.
1.19 Multiplication by t n
Theorem: If F (t) is a function of class A and if L {F (t)} = f ( p) then
n dn
L {t F (t)} = (− 1) n f ( p), where n = 1, 2, 3,......
dp n
(Rohilkhand 2000; Lucknow 06, 09, 11; Kanpur 08;
Avadh 10; Agra 01; Gorakhpur 11)
Proof: We shall prove this theorem by mathematical induction. By 1.18, we have
d
L { t F (t)} = (− 1)1 f ( p)
dp
T-23
∞ ∂ − pt m d m +1
or ∫0 {e t F (t)} dt = ( − 1)m f ( p),
∂p dp m + 1
∞ d m +1 f ( p )
or ∫0 e − pt ⋅ { t m +1
F (t)} dt = ( − 1) m + 1
dp m +1
m +1
d m + 1 f ( p)
or L {t F (t)} = ( − 1)m + 1
dp m +1
which shows that if the theorem is true for any particular value on n, it is true for the
next value of n. But we have already seen that the theorem is true for n = 1. Hence it is
true for n = 1 + 1 = 2 and n = 2 + 1 = 3, etc.
Therefore the theorem is true for every positive integral value of n.
1.20 Division by t
1 ∞
Theorem: If L {F (t)} = f ( p), then L F (t) = ∫p f ( x) dx
t
1
provided lim F (t) exists.
t→0 t
1
Proof: Let G (t) = F (t) i. e., F (t) = t G (t).
t
∴ L {F (t)} = L { t G (t)}
d
=− L { G (t)}, by the theorem of 1. 18
dp
T-24
d
or f ( p) = − L { G (t)}.
dp
Now integrating both sides with respect to p from p to ∞, we have
∞ ∞
− [ L { G(t)}]p = ∫p f ( p) dp
lim ∞
or − L { G (t)} + L { G (t)} = ∫p f ( p) dp
p→ ∞
[Note that L {G (t)} is a function of p]
∞
or 0 + L { G(t)} = ∫p f ( p) dp
∞ − pt
∵ plim
→∞
L { G (t)} = lim
p→ ∞
∫0 e G(t) dt = 0
1 ∞
or L F (t) = ∫p f ( x) dx
t
(iii) L {− a sin at }
Solution: We have
L {F ′ (t)} = pL {F (t)} − F (0 ) …(1)
(i) Here let F (t) = t, then F ′ (t) = 1 and F (0 ) = 0 .
∴ from (1), we have L {1} = pL { t } − 0
1 1 1 1
or L { t } = L {1} = ⋅ = 2 , p > 0 .
p p p p
a2
∴ L {− a sin at} = − ⋅
p2 + a2
2
2 p ( p2 − 3 a2 )
Example 15: Show that L { t cos at} = , p > 0.
( p2 + a2 )3
d2 p d − p2 + a2
∴ L { t2 cos at} = (− 1)2 =
dp2 ( p2 + a2 ) dp ( p2 + a2 )2
2 p ( p2 − 3 a2 )
= ⋅
( p2 + a2 )3
a2 d a 2 a3
=− ⋅ 2 = 2
p dp p + a ( p + a2 )2
2
a d p
= 2 2
+a 2
p +a dp p + a2
a a (a2 − p2 ) 2 a3
= + = .
p2 + a2 ( p2 + a2 )2 ( p2 + a2 )2
d2 3 d 3 3
= 2
2 +3 +2⋅ 2
dp p + 9 dp p2 + 9 p +9
18 p2 − 54 − 6p 6
= 2 3
+3⋅ 2 2
+ 2
( p + 9) ( p + 9) p + 9
6 p4 − 18 p3 + 126 p2 − 162 p + 432
= .
( p2 + 9)3
d 1 2p
∴ L { t sin t} = − 2 = 2
dp p + 1 ( p + 1)2
2 ( p + 1) 2p + 2
∴ L { t e − t sin t} = 2 2
= ⋅
[( p + 1) + 1] ( p + 2 p + 2)2
2
sin t −1 1 sin at
Example 19: Prove that L
= tan and hence find L ⋅ Does the Laplace
t p t
cos at
transform of exist?
t (Lucknow 2009; Gorakhpur 05; Rohilkhand 10; Avadh 13)
∞
sin at ∞ x 1
L = ∫p 2 2
dx = log ( x2 + a2 )
t x +a
2 p
1 lim 1
= log ( x 2 + a2 ) − log ( p2 + a2 ),
2 x→ ∞ 2
lim
which does not exist since log ( x2 + a2 ) is infinite.
x→ ∞
T-27
cos at
Hence L does not exist.
t
Example 20: If L {F (t), t → p} = f ( p),
t F (u) 1 ∞
show that L ∫ du, t → p = ∫ f ( y) dy .
0 u p p
t sin u cot −1 p
Hence show that L ∫ du, t → p = ⋅
0 u p
Solution: From 1.17, we have
t 1
L ∫ F (u) du = f ( p) …(1)
0 p
where f ( p) = L {F (t)}.
F (t)
Let G (t) = .
t
F (t) ∞
Then L { G(t)} = L = ∫p f ( y) dy [from article 1.20]
t
= g ( p), say.
∴ From (1), we have
t 1
L ∫ G (u) du = g ( p)
0 p
t 1 ∞
or L ∫ F (u) du = ∫ f ( y) dy …(2)
0 p p
Deduction. Let F (t) = sin t
1
so that f ( p) = L {sin t} = 2
.
p +1
∴ From (2), we have
t sin u 1 ∞ dy 1
L ∫ du = ∫ 2
= [tan−1 y]∞
p
0 u p p y +1 p
1 π 1
= − tan−1 p = cot −1 p.
p 2 p
Comprehensive Exercise 3
n!
(ii) Show that L { t n . e at} = , p > a.
( p − a)n +1
4. (i) Find L { t (3 sin 2 t − 2 cos 2 t)}.
(ii) Find L.T. of f (t) = sin α t + t cos α t.
2 a ( p − a)
5. Prove that L { t e at sin at} = 2 ⋅
( p − 2 ap + 2 a2 )2 (Kanpur 2010)
√π
6. Given L {sin √ t } = e −1 /4 p , show that
2 p3 /2
cos √ t π −1 / 4 p
L = e .
√t p
sinh t t sin x
7. Find L ⋅ 8. Find L ∫ dx .
t 0 x
9. Prove that if L {F (t)} = f ( p), then
∞ F (t) ∞
∫ 0 t dt = ∫ 0 f ( x) dx provided that the integral converges.
10. Find the Laplace transform of F (t) defined as
t + 1, 0 ≤ t ≤ 2
F (t) =
3, t>2.
Also determine L { F ′ (t)}.
A nswers 3
p2 − a2 p2 + 9
1. (i) , p > 0. (ii) , p > 0.
( p2 + a2 )2 ( p2 − 9)2
2 a (3 p2 − a2 ) 2
2. (i) 2 2 3
, p > 0. (ii) .
(p +a ) ( p − 2)3
8 + 12 p − 2 p2 (α + 1) p2 + (α − 1) α2
4. (i) . (ii) .
( p2 + 4)2 ( p2 + α2 )2
1 p+1
7. log ⋅ 8. (1 / p) cot −1 p.
2 p−1
1 1 − e− 2 p
10. [1 + p − e − 2 p ], ⋅
p2 p
∞ (e − at − e − bt)
Example 21: Evaluate ∫0 dt.
t
Solution: Let F (t) = e − at − e − bt .
1 1
Then L {F (t)} = L { e − at} − L { e − bt} = − = f ( p), say.
p+ a p+ b
∞
F (t) ∞ ∞ 1 − 1 dx = log ( x + a)
∴ L = ∫p f ( x) dx = ∫p
t x + a x + b ( x + b) p
lim x+a p+ a
= log − log
x→ ∞ x+b p+ b
lim 1 + (a / x) p+ a
= log − log
x→ ∞ 1 + (b / x) p+ b
p+ a p+ b
= 0 − log = log ⋅
p+ b p+ a
F (t) ∞ − pt e − at − e − bt p+ b
Thus L = ∫0 e ⋅ dt = log
t t p+ a
∞ e − at − e − bt b
∫0 dt = log ⋅
t a
∞
Example 22: Prove that ∫0 t3 e − t sin t dt = 0 .
(Gorakphur 2006)
3
Solution: Let F (t) = t sin t.
d3
∴ L {F (t)} = L { t3 sin t} = (−1)3 L {sin t}
dp 3
d3 1 d2 2p d 2 − 6 p2
=− = − − =
dp3 p2 + 1 dp2 ( p2 + 1)2 dp ( p2 + 1)3
∞ 24 ( p2 − 1) p
or ∫0 e − pt . t3 sin t dt = ⋅
( p2 + 1)4
Comprehensive Exercise 4
∞ sin t π
1. Show that ∫0 dt = ⋅
t 2 (Purvanchal 2007; Meerut 13, 13B; Rohilkhand 10, 14)
∞ e− t − e− 3 t
2. Evaluate ∫0 dt.
t (Rohilkhand 2010)
∞ − 3t 3
3. Show that ∫0 te sin t dt = ⋅
50
∞ 3
4. Show that ∫0 t e − 2 t cos t dt = ⋅
25
∞ e − x sin x π
5. Prove that ∫ dx = ⋅
0 x 4 (Rohilkhand 2009, 11)
∞ −3 t
6. Evaluate ∫0 te cos 4 t dt.
(Gorakhpur 2007, 09)
A nswers 4
5
3. log 3. 6. ⋅
169
T 2T 3T
= ∫0 e − pt F (t) dt + ∫T e − pt F (t) dt + ∫ 2T e − pt F (t) dt + ...
T T
= ∫0 e − pt F (t) dt + ∫0 e − p (u + T ) F (u + T ) du
T
+∫ e − p (u + 2T ) F (u + 2T ) du + ... ,
0
T
+ e −2 pT ∫ e − pu F (u) du + ...
0
T-31
T
= (1 + e − pT + e −2 pT + ...) ∫ e − pu F (u) du
0
T
∫0 e − pt F (t) dt
1 T − pu
= ∫0 e F (u) du = .
1 − e − pT 1 − e − pT
2 a2 p
and (ii) L { sinh at sin at } = ⋅
p4 + 4 a4
T-32
a
Solution: We have L {sinh at} = = f ( p), say.
p2 − a2
∴ L { e iat sinh at} = f ( p − ia)
a a
= 2 2
=
( p − ia) − a ( p − 2 a2 ) − 2 iap
2
a {( p2 − 2 a2 ) + 2 iap}
=
( p2 − 2 a2 )2 − (2 ipa)2
a ( p2 − 2 a2 ) + 2 ia2 p
or L {sinh at (cos at + i sin at)} =
p4 + 4 a4
a ( p2 − 2 a2 ) 2 a2 p
or L{sinh at cos at} + i L{sinh at sin at} = +i ⋅
p4 + 4 a4 p4 + 4 a4
a ( p2 − 2 a2 )
Hence L {sinh at cos at} =
p4 + 4 a4
2 a2 p
and L {sinh at sin at} = .
p4 + 4 a4
1 3 6 6
Example 24: Show that L {(1 + t e − t )3} = + + + ⋅
p ( p + 1)2 ( p + 2)3 ( p + 3)4
Solution: We have
L {(1 + t e − t )3} = L {1 + 3 t e − t + 3 t2 e −2 t + t3 e −3 t}
d d2 d3
= L {1} + 3 (−1) L { e − t} + 3 (−1)2 2 L { e −2 t} + (−1)3 3 L {e −3 t}
dp dp dp
1 1! 2! 3!
= + 3. 2
+ 3. 3
+
p ( p + 1) ( p + 2) ( p + 3)4
1 3 6 6
= + 2
+ 3
+ .
p ( p + 1) ( p + 2) ( p + 3)4
∞
Aliter: L {(1 + t e − t )3} = ∫0 (1 + t e − t )3 . e − pt dt
∞
= ∫0 [e − pt + 3 t e − ( p + 1)t + 3 t2 e − ( p + 2)t + t3 e − ( p + 3)t ] dt
∞ ∞ ∞
= ∫0 e − pt dt + 3 ∫ te − ( p + 1)t dt + 3 ∫ t2 e − ( p + 2)t dt
0 0
∞
+ ∫0 t 3 e − ( p +3)t dt
∞ ∞
= ∫0 e − pt t1−1 dt + 3 ∫ e − ( p + 1)t t2 − 1 dt
0
∞ ∞
+3∫ e − ( p + 2)t t3 − 1 dt + ∫0 e − ( p + 3)t t 4 − 1 dt
0
T-33
1 3 6 6
= + 2
+ 3
+ , p > 0.
p ( p + 1) ( p + 2) ( p + 3)4
∞ − at n−1 Γ(n)
e t dt = n , if a > 0 and n > 0
∫ 0
∵
a
1
Example 25: Prove that L { J0 (t)} =
√ (1 + p2 )
and hence deduce that (Avadh 2007)
1
(i) L { J0 (at)} =
√ ( p + a2 )
2
p
(ii) L { t J0 (at)} =
( p + a2 )3 /2
2 (Rohilkhand 2002)
1
(iii) L { e − at J0 (at)} =
√ ( p + 2 ap + 2 a2 )
2
∞
(iv) ∫0 J0 (t) dt = 1.
d d 1 p
(ii) L { t J0 (at)} = − L{ J0 (at)} = − = 2
dp dp √ ( p + a ) ( p + a2 )3 /2
2 2
1 1
∴ L { e − at J0 (at)} = 2 2 ∵ L { J0 (at)} = 2 2
√ [( p + a) + a ] √ ( p + a )
1
= .
√ ( p + 2 ap + 2 a2 )
2
T-34
∞ 1
(iv) We have L { J0 (t)} = ∫0 e − pt J0 (t) dt = .
√ (1 + p2 )
∞
∴ Putting p = 0, we have ∫0 J0 (t) dt = 1.
p
Example 26: Prove that L { J1 (t)} = 1 − where J1 (t) is the Bessel function of order
√ ( p 2 + 1)
one (Kanpur 2009)
1
and hence deduce that L {t J1 (t)} = .
( p + 1)3 /2
2
d p 1
=− 1 − = 2 .
dp √ ( p + 1) ( p + 1)3 /2
2
T-35
2 √ π 1 √π 1 1 √ π 1⋅ 3 1 √ π 1⋅ 3 ⋅ 5 1
= ⋅ 3 /2 − 5 / 2
+ ⋅ 7 / 2
− ⋅ 9 / 2
+ ...
√π 2 p 2 2 p 2 2 ⋅4 p 2 2 ⋅4 ⋅6 p
−1 / 2
1 1 1 1⋅ 3 1 1⋅ 3 ⋅ 5 1 1 1
= 3 /2
1− ⋅ + 2
− 3
+ ... = 3 /2 1 +
p 2 p 2 ⋅ 4 p 2 ⋅ 4 ⋅ 6 p p p
1
= ⋅ (Rohilkhand 2004)
p √ ( p + 1)
1 p
Since L { F (at)} = f , where f ( p) = L { F (t)},
a a
∴ L {erf (2 √ t)} = L {erf √ (4 t)}
1 1 2
= ⋅ = ⋅
4 p p p √ ( p + 4)
+ 1
4 4
d d 2 3p + 8
Hence L { t. erf (2 √ t)} = − L {erf (2 √ t)} = − = ⋅
dp dp p √ ( p + 4) p2 ( p + 4)3 /2
t3 t5 t7
=t− + − + ...
3 (3 !) 5 (5 !) 7 (7 !)
1 1 1
∴ L { Si (t)} = L { t } − L { t 3} + L { t 5} − ⋅ L { t7} + ...
3 (3 !) 5 (5 !) 7 (7 !)
1! 1 3! 1 5! 1 7!
= − ⋅ + ⋅ − ⋅ + ...
p2 3 (3 !) p4 5 (5 !) p6 7 (7 !) p8
1 1 1 1 1 1 1 1
= − ⋅ 3 + ⋅ 5 − ⋅ 7 + ...
p p 3 p 5 p 7 p
1 1
= tan−1 , by Gregory’s series.
p p
cos u
∞
Solution: L { Ci (t)} = L ∫ du ⋅
t u
∞ cos u t cos u
Let F (t) = ∫ du = − ∫ du
t u ∞ u
cos t
so that F ′ (t) = − or t F ′ (t) = − cos t.
t
∴ L { t F ′ (t)} = L { − cos t}
d p
or − L{F ′ (t)} = − 2
dp p +1
d p
or [ p f ( p) − F (0 )] = 2 where f ( p) = L {F (t)}
dp p +1
d p
or [ p f ( p)] = 2 , since F (0 ) is constant.
dp p +1
1
Integrating, p f ( p) = log ( p2 + 1) + C (constant) …(1)
2
But from the final-value theorem 1.16,
lim lim
p f ( p) = F (t) = 0 .
p→ 0 t→ ∞
∴ from (1) as p → 0 we have 0 = 0 + C or C = 0 .
1
∴ from (1), p f ( p) = log ( p2 + 1)
2
log ( p2 + 1)
or f ( p) = L {F (t)} = L { Ci (t)} = .
2p
3 t, 0 < t < 2
Example 30: If F (t) =
6, 2 < t < 4,
find L {F (t)} where F (t) has period 4.
Solution: Here F (t) is a periodic function with period T = 4.
T-37
1 6 −2 p 3 e −2 p 3 6 6
= − e − + 2 − e −4 p + e −2 p
1 − e −4 p p p2
p p p
3 − 3 e− 2 p − 6 p e− 4 p
= .
p2 (1 − e − 4 p )
Comprehensive Exercise 5
t 1 1 1
1. Given L 2 = 3 /2 , show that 1 /2 = L ⋅
π
p p √ (πt)
2. Find (i) L {F (t)} and (ii) L {F ′ (t)}, for the function given by
2 t , 0 ≤ t ≤ 1
F (t) =
t, t > 1.
25 30 9
3. Show that L {(5 e2 t − 3)2} = − + , p > 4.
p−4 p−2 p
sin2 t 1 p2 + 4
4. Show that L = log 2 ⋅ (Rohilkhand 2010; Kashi 14)
t 4 p
∞ sin2 t
π
5. Show that ∫ dt = . (Rohilkhand 2003)
0t2 2
∞ cos 6 t − cos 4 t 2
6. Show that ∫ dt = log ⋅
0 t 3
2
7. Prove that L { J0 (a √ t)} = (1 / p) e − (a /4 p)
.
8. If F (t) = t2 , 0 < t < 2 and F (t + 2) = F (t), find L {F (t)}.
[Hint. Here F (t) is a periodic function with period T = 2]
sin t, 0 < t < π
9. Compute L {F (t)}, if F (t) = where F (t) has period 2π.
0 , π < t < 2 π,
10. Find the Laplace transform of the Heaviside’s unit step function H (t − a).
t 1 − e −2 x
11. Find Laplace transform of ∫ dx.
0 x
T-38
∞ e− u log ( p + 1)
12. If E (t) = ∫t du, show that L { E (t)} = ⋅
u p
∞
13. Evaluate ∫0 e − t erf √ t dt.
A nswers 5
2 e −p
2. (i) 2 / p2 − (1/ p + 1/ p2 ) e − p (ii) −
p p
− (4 p2 + 4 p + 2) e −2 p + 2 1
8. 3 −2 p
9.
p (1 − e ) ( p + 1) . (1 − e − pπ )
2
e − ap 1 2
10. 11. log 1 +
p p p
12. 1/ √ 2
No. Operation F (t ) L { F (t )} = f ( p)
4. Change of scale 1 p
F (at) f
property a a
F ′ (t) pf ( p) − F (0 )
5. Differentiation
theorems n−1
F n (t) pn f ( p) − ∑ p n −1− r F r (0 )
r =0
T-39
t F (t) − f ′ ( p)
6. Multiplication
theorems n
t F (t) dn
(− 1)n n
f ( p)
dp
1 ∞
7. Division theorem F (t) ∫p f ( x) dx
t
t 1
8. Integral theorem ∫0 F ( x) dx f ( p)
p
lim lim
9. Initial value theorem F (t) = p L {F (t)}
t→0 p→ ∞
lim lim
10. Final-value theorem F (t) = p L {F (t)}
t→ ∞ p→ 0
T
∫0 e − pt F (t) dt
11. Fundamental L {F (t)} = ,
1 − e − pT
theorem for
F (t) is periodic function of period T
Periodic Functions
1 1
(c) (d) ⋅ (Rohilkhand 2003)
p − a2
2
p2 + a2
6. If L { F (t)} = f ( p) then L { t F (t)} is
(a) f ′ ( p) (b) − f ′ ( p)
∞ 1
(c) ∫p f ( x) dx (d) f ( p). (Rohilkhand 2002)
p
sin t
7. The value of Laplace transform of is
t
1 2
(a) tan−1 (b) tan−1
p p
(c) tan−1 p (d) None of these.
2 ap
(c) (d) None of these.
( p + a2 )2
2
T-41
∞
12. The value of ∫ J0 (t) dt is
0
(a) 1 (b) 0
(c) −1 (d) None of these.
∞ sin t
13. The value of ∫ dt is
0 t
π
(a) 0 (b)
2
(c) 1 (d) None of these.
1
14. If L {erf √ t} = then L { e3 t erf √ t} is
p √ ( p + 1)
1 1
(a) (b)
( p − 3) √ ( p − 2) ( p − 3) √ ( p + 2)
3
(c) (d) None of these.
( p − 3) √ ( p − 2)
the …… of the function F (t) and the function K ( p, t) is called the …… of the
transformation.
∞ − pt
2. The integral ∫0 e F (t) dt is called the …… transform of the function F (t).
True or False
Write ‘T’ for true and ‘F’ for false statement.
1. The Laplace transformation is a linear transformation.
2
2. The function e t is not of exponential order as t → ∞.
T-42
6
8. The value of L { t3 e −3 t } = ⋅
( p + 3)4
A nswers
Multiple Choice Questions
1. (c) 2. (c) 3. (b)
4. (a) 5. (a) 6. (b)
7. (a) 8. (c) 9. (d)
10. (b) 11. (c) 12. (a)
13. (b) 14. (a)
dn 1
5. (− 1)n n
f ( p) 6.
dp √ (1 + p2 )
π
7. 8. F (u)
2
True or False
1. T 2. T
3. F 4. F
5. F 6. F
7. T 8. T
¨
T-43
2
T he I nverse L aplace T ransform
function.
If we restrict ourselves to functions F (t) which are sectionally continuous in every finite
interval 0 ≤ t ≤ t0 and of exponential order for t > t0 , then the inverse Laplace transform of
f ( p) i.e., L−1 { f ( p)} = F (t) is unique.
Proof: We have
L { c1 F1 (t) + c2 F2 (t)} = c1 L {F1 (t)} + c2 L {F2 (t)}
= c1 f1( p) + c2 f2 ( p).
∴ L−1 { c1 f1 ( p) + c2 f2 ( p)} = c1 F1 (t) + c2 F2 (t)
= c1 L−1 { f1 ( p)} + c2 L−1 { f2 ( p)}.
p p
∵ L {cos at} = 2 2
, ∴ L−1 2 2
= cos at.
p +a p + a
(vi) Inverse L.T. of 1 /( p2 − a2 ), p > | a | .
a 1 1
∵ L {sinh at} = 2 2
, ∴ L−1 2 2
= sinh at.
p −a p − a a
T-45
1
1. 1
p
1 tn
2. n +1
, n is a positive integer
p n!
1
3. e at
p− a
1 1
4. 2 2
sin at
p +a a
p
5. cos at
p + a2
2
1 1
6. sinh at
p2 − a2 a
p
7. cosh at
p − a2
2
1 tn
8. n +1
,n> −1
p Γ (n + 1)
If we know the inverse L.T. given in the above table then nearly all the inverse Laplace
transforms can be obtained by using the general theorems, which we shall give later on.
1 1
Example 1: Find (i) L−1 4 (ii) L−1 2 .
p p + 4 (Meerut 2013B)
1 t 4 − 1 t3
Solution: (i) L−1 4 = =
p 3! 6
1 −1 1 1
(ii) L−1 2 =L 2 2
= sin 2 t
p + 4 p +2 2
T-46
5 p− 2 √ p+1 7 1
= L−1 2 + 2
− .
p p 3 p + (2 / 3)
1 1 1 7 1
= 6 L−1 2 + L−1 − 2 L−1 3 /2 − L−1
p p p 3 p + (2 / 3)
t 2 −1 t 3 /2 −1 7
=6 + 1− 2 . − e −2 t /3
1! Γ (3 / 2) 3
= 6 t + 1 − 4 √ (t / π) − (7 / 3) e −2 t /3 .
3 5 7
1 1 −1 1 1 (1/ p) (1/ p) (1/ p)
Solution: L−1 sin =L − + − + ...
p p p p 3! 5! 7!
1 1 −1 1 1 −1 1 1 −1 1
= L−1 2 − L 4 + L 6 − L 8 + ...
p 3! p 5! p 7! p
t3 t5 t7
=t− + − + ... .
(3 !)2 (5 !)2 (7 !)2
Comprehensive Exercise 1
3 3p + 2 3 p − 27 6 − 30 √ p
7. 2
+ 3
− 2
+ ⋅
p −3 p p +9 p4
3 ( p2 − 1)2 4 p − 18 ( p + 1) (2 − √ p)
8. 5
+ 2
+ ⋅
2p 9− p p5 /2
1 1 t2 t4 t6
9. Show that L−1 cos = 1− 2
+ 2
− + ...
p p (2 !) (4 !) (6 !)2
A nswers 1
1. 4 e 2 t. 2. cos √ 2 t + 6 cosh 4 t + 3 e 3 t .
3. e − t + 1. 4. 2 cosh 3 t − (5 / 3) sinh 3 t.
t 8 t 7
5. 6 − t − e −2 t /3 .
π 3 π 3
1 4t 4 4t 2 3t 3 3t
6. 3 e 3 t /2 − sinh − cosh + sin − cos ⋅
4 3 9 3 3 4 8 4
7. √ 3 sinh (√ 3 t) + 3 t + t 2 − 3 cos 3 t + 9 sin 3 t + t 3 − 16 t 2 √ (t / π).
1 3 2 1 4 8
8. − t + t − 4 cosh 3 t + 6 sinh 3 t + 4 √ (t / π) + t √ (t / π) − t.
2 2 16 3
∞ ∞ − pt
∴ f ( p − a) = ∫0 e −(p − a)t F (t) dt = ∫0 e .{ e at ⋅ F (t)} dt = L { e at F (t)}.
∞ − p (t + a)
∴ e − ap f ( p) = ∫0 e F (t) dt
T-48
∞ − px
= ∫0 e F ( x − a) dx, putting t + a = x, so that dt = dx
a − px ∞ − px
= ∫0 e ⋅ 0 dx + ∫a e F ( x − a) dx
a − pt ∞ − pt
= ∫0 e ⋅ 0 dt + ∫a e F (t − a) dt
∞ − pt
= ∫0 e G (t) dt = L { G (t)} ,
F (t − a), t > a
where G (t) =
0, t < a.
Note: In terms of Heaviside’s unit step function H (t − a) this theorem can be stated as :
If L−1 { f ( p)} = F (t),
then L−1 { e − ap f ( p)} = F (t − a). H (t − a).
∞ − apt 1 ∞ − px x
f (ap) = ∫0 e F (t) dt = e F dx .
a ∫0
∴
a
putting at = x, so that dt = (1 / a) dx
1 ∞ − pt t 1 t 1 t
= ∫ e F dt = L F = L F .
a 0 a a a a a
Hence L−1 { f (ap)} = (1 / a) F (t / a).
3p + 7
Example 4: Evaluate L−1 2 .
p − 2 p − 3 (Bundelkhand 2013)
3p + 7 −1
3 ( p − 1) + 10
Solution: L−1 2 =L 2
p − 2 p − 3 ( p − 1) − 4
T-49
3 ( p − 1) 10
= L−1 2
+ 2
( p − 1) − 4 ( p − 1) − 4
p−1 1
= 3 L−1 2 + 10 L
−1
2
( p − 1) − 4 ( p − 1) − 4
p 1
= 3 e t L−1 2 2
+ 10 e t L−1 2 2
p −2 p −2
= 3 e t cosh 2 t + 5 e t sinh 2 t = 4 e 3 t − e − t .
e −1 / p cos 2 √ t e−a / p
Example 5: If L−1 1 /2 = , find L−1 1 /2 , where a > 0.
p √ (πt) p
e −1 / p cos 2 √ t
Solution: Since L−1 1 /2 = ,
p √ (πt)
e −1 / pk 1 cos 2 √ (t / k ) e −1 / pk cos 2 √ (t / k )
∴ L−1 1 /2
= or L−1 1 /2 = .
( pk ) k √ (πt / k ) p √ (πt)
e − a / p cos 2 √ (at)
Taking k = 1 / a, we have L−1 1 /2 = ⋅
p √ (πt)
1 −1 1 1 −1 ( p + 2) − 6
= L − L 2
8 p 8 ( p + 2) + 4
1 1 −2 t −1 p − 6
= − e L 2
8 8 p + 4
1 1 −2 t −1 p 1
= − e L 2 2
−6 L−1 2 2
8 8 p +2 p + 2
1
= [1 − e −2 t (cos 2 t − 3 sin 2 t)].
8
p 2 t 1
Example 8: Prove that L−1 4 2 = sinh . sin . √ 3 t.
p + p + 1 √ 3 2 2
(Kanpur 2010; Meerut 13B; Kashi 14; Rohilkhand 14)
p p
Solution: We have L−1 4 2 =L
−1
2 2 2
p + p + 1 ( p + 1) − p
p
= L−1 2 2
( p − p + 1) ( p + p + 1)
1 ( p2 + p + 1) − ( p2 − p + 1)
= L−1 2 2
2 ( p − p + 1) ( p + p + 1)
1 1
= L−1 2
− 2
2 ( p − p + 1) 2 ( p + p + 1)
1 −1 1 1 −1 1
= L − L
2 1 2 3 2 1 2 3
( p − ) + ( p + ) +
2 4 2 4
1 t / 2 −1 1 1 − t / 2 −1 1
= e L − e L
2 2 1 2 2 2 1 2
p + ( √ 3) p + ( √ 3)
2 2
1 t /2 2 t 1 2 t
= e sin √ 3 − e − t /2 . sin √ 3
2 √3 2 2 √3 2
1 t 2 t t
= (e t /2 − e − t /2 ) sin √ 3 = sinh sin √ 3 ⋅
√3 2 √3 2 2
1
Example 9: Evaluate (i) L−1 2 2
,
( p + 4) ( p + 1) (Meerut 2013)
3 p3 − 3 p2 − 40 p + 36
(ii) L−1 ⋅
( p2 − 4)2
T-51
1
Solution: (i) We have L−1 2 2
( p + 4) ( p + 1)
2 1 2p + 3
= L−1 + 2
−
25 ( p + 1) 5 ( p + 1) 25 ( p2 + 4)
1 −1 1 −1 1
= 2 L +5L 2
25 p + 1 ( p + 1)
p −1 1
− 2 L−1 2 −3 L 2
p + 4 p + 4
1
= [2 e − t + 5 e − t L−1 {1/ p2 } − 2 cos 2 t − (3 / 2) sin 2 t]
25
1 −t
= [e (2 + 5 t) − 2 cos 2 t − (3 / 2) sin 2 t ].
25
3 p3 − 3 p2 − 40 p + 36
(ii) L−1
( p2 − 4)2
3 p3 − 3 p2 − 40 p + 36
= L−1 2 2
( p − 2) ( p + 2)
2 3 5
= L−1 − 2
+ + 2
( p − 2) p + 2 ( p + 2)
1 1 1
= − 2 L−1 2
+ 3 L−1 + 5L
−1
2
( p − 2) p + 2 ( p + 2)
= − 2 e2 t L−1 {1/ p2 } + 3 e −2 t + 5 e −2 t L−1 {1/ p2 }
= − 2 e2 t . t + 3 e −2 t + 5 e −2 t . t
= (5 t + 3) e −2 t − 2 t e2 t .
p2 1
Example 10: Show that L−1 4 4
= (cosh at sin at + sinh at cos at).
p + 4 a 2a
Solution: We have
p2 −1 p2
L−1 4 4 = L 2 2 2 2 2
p + 4a ( p + 2a ) − 4a p
p2
= L−1 2 2 2 2
( p − 2 ap + 2 a ) ( p + 2 ap + 2 a )
1 p ( p2 + 2 ap + 2 a2 ) − p ( p2 − 2 ap + 2 a2 )
= L−1 ⋅
4a ( p2 + 2 ap + 2 a2 ) ( p2 − 2 ap + 2 a2 )
T-52
p p
= L−1 2 2
− 2 2
4 a ( p − 2 ap + 2 a ) 4 a ( p + 2 ap + 2 a )
1 −1 ( p − a) + a 1 −1 ( p + a) − a
= L 2 2 − L 2 2
4a ( p − a) + a 4 a ( p + a) + a
1 at −1 p + a 1 − at −1 p − a
= e L 2 2
− e L 2 2
4a p + a 4 a p + a
1 at −1 p −1 1
= e L 2 2 +aL 2 2
4a p +a p + a
1 − at −1 p −1 1
e − L 2 2 −aL 2 2
4a p +a p + a
1 at 1 − at
= e (cos at + sin at) − e (cos at − sin at)
4a 4a
1 e at + e − at e at − e − at
= sin at + cos at
2 a 2 2
1
= (cosh at sin at + sinh at cos at).
2a
( p + 1) e − πp
Example 11: Find L−1 2 .
p + p+1
( p + 1) + 1
p+1
−1 −1 2 2
Solution: We have L 2 =L 1 3
p + p + 1 ( p + ) + 2
2 4
p+ 1
− t / /2 −1 2
=e L
3
p2 +
4
p 1 − t / 2 −1 1
= e − t /2 L−1 2 2
+ e L 2 2
p + (√ 3 / 2) 2 p + (√ 3 / 2)
1
= e − t /2 cos (√ 3 t / 2) + e − t /2 . (2 / √ 3) sin (√ 3 t / 2)
2
e − t /2
= [√ 3 cos (√ 3 t / 2) + sin (√ 3 t / 2)].
√3
− (t − π ) / 2
( p + 1 ) e − πp e √3 √3
∴ L−1 √ 3 cos 2 (t − π) + sin 2 (t − π) , t > π
2 = √3
p + p + 1
0 ,t< π
e − (t − π ) / 2 √3 √3
= √ 3 cos 2 (t − π) + sin 2 (t − π) H (t − π).
√3
T-53
Comprehensive Exercise 2
p 1 32 p
5. If L−1 2 2
= t sin t, find L−1 2 2
⋅
( p + 1) 2 (16 p + 1)
e −5 p −1 e
− 4p
6. Find (i) L−1 4 (ii) L 4
⋅
( p − 2) ( p − 3)
e4 − 3 p
7. Find L−1 5 /2
⋅
( p + 4)
8. Find the inverse L.T. of e − 3 p / p3 . (Lucknow 2009)
9. Prove the following :
e − pπ
(i) L−1 2 = − sin t. H (t − π)
p + 1 (Kanpur 2008)
p e −2 pπ /3
(ii) L−1 2 = cos 3 (t − 2 π / 3) . H (t − 2 π / 3).
p +9
10. Prove the following :
pe − ap
(i) L−1 2 2
= cos h ω (t − a). H (t − a), a > 0
p −ω
3 (1 + e − pπ )
(ii) L−1 2 = − sin 3 t. H (t − π) + sin 3 t.
p +9
11. Evaluate
p+2
(i) L−1 2 (Kanpur 2008; Lucknow 10; Avadh 13)
p − 2 p + 5
T-54
p+8
(ii) L−1 2 ⋅
p + 8 p + 5
4p + 5 1 t 1 −2 t
12. Prove that L−1 2
t
= 3t e + e − e .
( p − 1) ( p + 2) 3 3 (Purvanchal 2007)
4p + 5
13. Evaluate (i) L−1 2
( p − 4) ( p + 3)
5 p2 − 15 p − 11
(ii) L−1 3
.
( p + 1) ( p − 2) (Gorakhpur 2009)
2p + 1 1
14. Prove that L−1 2 2
= t (e t − e −2 t ).
( p + 2) ( p − 1) 3 (Gorakhpur 2007)
p 1
15. Prove that L−1 2 2 = sin t sinh t.
( p − 2 p + 2) ( p + 2 p + 2) 2
A nswers 2
1. (i) e 3 t sin t
4 − 3t 1 − t
(ii) − e + e (4 cos t − 3 sin t)
5 5
t n−1 2 −t t (3 − 2 t)
2. (i) e − at ⋅ (ii) e
(n − 1) ! 3 π
1 −t 1 −3 t /2
3. (i) e (4 t 3 − t 4 ) (ii) e (6 − 5 t)
24 8
1 3 1
4. (i) (t − ) e t + e − t (ii) [(3 t − 1) e t + e − 2 t ]
2 2 9
1 1
5. t sin ( t)
4 4
1 1
6. (i) (t − 5)3 e 2 (t − 5) . H (t − 5) (ii) (t − 4)3 e 3 (t − 4) . H (t − 4)
9 6
4 1
7. (1/ √ π) . (t − 3)3 /2 e − 4 (t − 4) . H (t − 3) 8. (t − 3)3 . H (t − 3)
3 2
3
9. (i) e t [cos 2 t + sin 2 t]
2
(ii) e − 4 t [cosh (√ 11t) + (4 / √ 11) sinh (√ 11t)]
1 − 3t 1 4t 1 7
13. (i) − e + e + 3 te 4 t (ii) (− e − t + e2 t ) + 4 te2 t − t2 e2 t
7 7 3 2
T-55
Proof: We have
L {F ′ (t)} = p L {F (t)} − F (0 ) = p f ( p). [ ∵ F (0 ) = 0 ]
−1
∴ L { p f ( p)} = F ′ (t).
Note 1: If F (0 ) ≠ 0 , then
L−1 { p f ( p) − F (0 )} = F ′ (t)
t
Then G ′ (t) = ∫ 0 F ( x) dx, G ′ ′ (t) = F (t). Also G (0 ) = G ′ (0 ) = 0 .
p
Example 12: Find L−1 2 2 2
⋅
( p + a ) (Purvanchal 2010)
p 1 d 1 1 −1 d 1
Solution: L−1 2
= L−1 − ⋅
2 2
2 2
= − L 2
2
( p + a ) 2 dp p +a 2 dp p + a
T-57
1 1 t
=− t . (− 1)1 L−1 2 2 = sin at.
2 p + a 2 a
2
p
Example 13: Find L−1 2 2
⋅
( p + 4)
p
Solution: Let f ( p) = .
( p + 4)2
2
p t
∴ L−1 { f ( p)} = L−1 2 2
= sin 2 t = F (t)
( p + 4) 4
[See Ex. 12, here a = 2]
and F (0 ) = 0 .
∴ L−1 { p f ( p)} = F ′ (t)
p2 d 1 1
or L−1 2 2
= ( t sin 2 t) = (sin 2 t + 2 t cos 2 t).
( p + 4) dt 4 4
p+2
Example 14: Find L−1 2 ⋅
p ( p + 3)
p+2 −1
( p + 3) − 1
Solution: We have L−1 2 =L 2
p ( p + 3) p ( p + 3)
1 1 −1 1 −1 1
= L−1 2 − 2 = L 2 − L 2 …(1)
p p ( p + 3) p p ( p + 3)
Now L−1 {1/ p2 } = t /1 ! = t.
1
Also L−1 =e
−3 t
= F (t), say.
p + 3
1 t
∴ we have L−1 = ∫ 0 F ( x) dx.
p ( p + 3) (by theorem I of article 2.14)
t −3 x 1 −3 t
=∫ e dx = (1 − e ) = F1 (t) say,
0 3
1 t 1 t
∴ L−1 2 = ∫ 0 F1 ( x) dx =
−3 x
∫ 0 (1 − e ) dx
p ( p + 3) 3
1 1
= t + (e −3 t − 1).
3 9
∴ From (1), we have
p+2 1 1 −3 t 2 1 1
L−1 2 = t − t − (e − 1) = t − e −3 t + ⋅
p ( p + 3) 3 9 3 9 9
1
Example 15: Find (i) L−1 log 1 + 2 (Agra 2001)
p
T-58
1 1
(ii) L−1 log 1 + 2 ⋅
p p
t 2
=∫ (1 − cos x) dx .
0 x
Comprehensive Exercise 3
p+3
8. (i) L−1 log (Kanpur 2007)
p+2
1 p+3
(ii) L−1 log ⋅
p p+2
p 1 1
9. If L−1 2 2
= t sin t , find L−1 2 2
⋅
( p + 1) 2 ( p + 1) (Kanpur 2011)
1
10. L−1 tan−1 ⋅
p
A nswers 3
1 1 2 − at
1. (i) t2 e at (ii) t e .
2 2
1 t sinh at 1
2. (i) (ii) te − t sin t.
2 a 2
1 2
3. 2 (1 − cosh t) / t. 4. 1 − t + t − e− t .
2
1 2 1
5. t + cos t − 1. 6. 1 − e − t (1 + t + t2 ).
2 2
t 1 −x
7. (i) (e − t − e − 2 t ) / t (ii) ∫ (e − e −2 x ) dx .
0 x
t 1 −2 x
8. (i) (e − 2 t − e − 3 t ) / t (ii) ∫ (e − e −3 x ) dx .
0 x
1 sin t
9. (sin t − t cos t). 10. ⋅
2 t
2.15 Convolution
Let F (t) and G (t) be two functions of class A, then the convolution of the two functions F (t) and
G (t) denoted by F * G is defined by the relation
t
F*G= ∫0 F ( x) G (t − x) dx.
t
= −∫ F (t − y) G ( y) dy, putting t − x = y so that dx = − dy
0
t
=∫ G ( y) F (t − y) dy = G * F .
0
t t
=∫ F ( x) G (t − x) dx + ∫0 F ( x) H (t − x) dx
0
= F*G+ F* H.
(Rohilkhand 2003, 07; Avadh 06, 09, 14; Lucknow 06, 10, 11; Kanpur 09)
Proof: Here we shall prove that the Laplace transform of the convolution of two
functions is equal to the product of their Laplace transforms
t
i. e., L ∫ F ( x) G (t − x) dx = f ( p) g ( p)
0
If we change the order of integration, then the strip is taken parallel to Ot so that the
limits of t are from x to ∞ and those of x are from 0 to ∞.
∴ changing the order of integration in the double integral (1), we have
t ∞ ∞ e − pt F ( x) G (t − x) dt dx
L ∫ F ( x) G (t − x) dx = ∫
0 x = 0 ∫ t = x
∞ ∞
= ∫ x =0 e − px F ( x) ∫ e px e − pt G (t − x) dt dx
t = x
∞ ∞
= ∫ x =0 e − px F ( x) ∫ e − p (t − x)G (t − x) dt dx
t = x
∞ ∞
= ∫ x =0 e − px F ( x) ∫ e − pz G (z ) dz dx,
z = 0
∵ b b b
∫a f ( x) dx = ∫a f (t) dt = ∫a f (z ) dz
F ( p) ⋅ ( p − α r )
Ar = lim
p→ α r G ( p)
(p − αr) Form 0
= F (α r ) ⋅ lim
p→ α r G ( p) 0
1
= F (α r ) ⋅ lim [by L’ Hospital’s rule]
p → αr G ′ ( p)
1
= F (α r ) ⋅ ⋅
G ′ (α r )
F ( p) F (α1) 1 F (α2 ) 1
∴ = ⋅ + . + ...
G ( p) G ′ (α1) ( p − α1) G ′ (α2 ) ( p − α2 )
F (α r ) 1 F (α n) 1
+ ⋅ + ... + ⋅ + ... .
G ′ (α r ) ( p − α r ) G ′ (α n) p − α n
F ( p) F (α1) −1 1 F (α2 ) −1 1
Hence L−1 = ⋅L + L + ...
G ( p ) G ′ (α1) p − α1 G ′ (α2 ) ( p − α2 )
F (α r ) −1 1 F (α n) −1 1
+ L + ... + .L
G ′ (α r ) p − α r G ′ (α n) p − αn
F (α1) α 1 t F (α r ) α r t F (α n) α n t
= ⋅e + ... + ⋅e + ... + ⋅e
G ′ (α1) G ′ (α r ) G ′ (α n)
n F (α r ) α r t
= Σ e .
r =1 G ′ (α r )
p2
Example 16: Use the convolution theorem to find L−1 2 2 2
⋅
( p + a )
p
Solution: We have L−1 2 2
= cos at.
p +a
∴ By the convolution theorem, we have
p2 p p
L−1 2 2 2
= L−1 2 2
⋅ 2 2
( p + a ) p + a p + a
T-63
t
= ∫ 0 cos ax cos a (t − x) dx
t
= ∫0 cos ax (cos at cos ax + sin at sin ax) dx
t 2 t
= cos at ∫ 0 cos ax dx + sin at ∫ cos ax sin ax dx
0
1 t 1 t
= cos at ∫ (1 + cos 2 ax) dx + sin at ∫ sin 2 ax dx
2 0 2 0
t t
1 1 1 1
= cos at x + sin 2 ax + sin at. − cos 2 ax
2 2a 0 2 2 a 0
1 1 1
= cos at t + sin 2 at + sin at (1 − cos 2 at)
2 2a 4 a
1 1 1
= t cos at + sin at + (sin 2 at cos at − sin at cos 2 at)
2 4a 4a
1 1 1
= t cos at + [sin at + sin (2 at − at)] = [at cos at + sin at].
2 4a 2a
1
Example 17: Find L−1 , by the convolution integral and deduce the value of
√ p . ( p − a)
1
L−1 ⋅
p √ ( p + a)
1 −1 1 t1 /2 − 1 1
Solution: We have L−1 = L 1 /2 = 1
= = F1 (t), say
√ p p Γ (2 ) √π √t
1
and L−1 at
= e = F2 (t), say.
p− a
∴ By the convolution theorem , we have
1 t 1
L−1 = F1 (t) * F2 (t) = ∫0 ⋅ e a(t − x)dx
√ p ⋅ ( p − a) √π √ x
e at √(at) √ a − u2 2 u 2 u du
= ∫0 e ⋅ du, putting ax = u2 , so that dx =
√π u a a
e at 2 √(at) − u2 e at
=⋅ ∫ e du = erf (√ (at)).
√a √π 0 √a
1 −1 1
Deduction. L−1 =L
p √ ( p + a) ( p + a − a) √ ( p + a)
at
1 − at e 1
= e − at L−1 =e ⋅ erf (√ (at)) = erf (√ (at)).
( p − a) √ p √a √a
Example 18: Apply convolution theorem to prove that
1 Γ (m) ⋅ Γ(n)
B (m, n) = ∫ x m −1 (1 − x)n−1 dx = , m > 0, n > 0 .
0 Γ(m + n)
(Gorakhpur 2008; Agra 02)
Hence deduce that
T-64
π /2 1 Γ (m) ⋅ Γ (n)
∫0 sin2 m −1 θ ⋅ cos2 n−1 θ dθ = B (m, n) = ⋅
2 2 Γ (m + n)
Solution: Consider the function
t
F (t) = ∫0 x m −1 ⋅ (t − x)n−1 dx.
t
We have, F (t) = ∫0 F1 ( x) ⋅ F2 (t − x) dx, where F1 (t) = t m −1 and F2 (t) = t n−1
= F1 * F2 .
∴ L {F (t)} = L {F1 * F2 }
= L {F1 (t)}. L {F2 (t)}, by convolution theorem
Γ (m) Γ (n) Γ (m) Γ (n)
= L { t m − 1} . L {t n − 1} = ⋅ n = ⋅
pm p pm+n
t m −1 Γ (m) Γ (n)
∴ F (t) = ∫0 x ⋅ (t − x)n − 1 dx = L−1 m+n
p
1 Γ (m) Γ (n) m + n − 1
= Γ (m) ⋅ Γ (n) ⋅ L−1 m + n = t .
p Γ (m + n)
Taking t = 1, we have
1
B (m, n) = ∫0 x m − 1 (1 − x)n − 1 dx
Γ (m) Γ (n)
= ⋅ …(1)
Γ (m + n)
Deduction: Taking x = sin2 θ, so that dx = 2 sin θ cos θ dθ. From (1), we have
π /2 Γ (m) Γ (n)
2 ∫0 sin2 m − 1 θ cos2 n − 1 θ dθ = ⋅
Γ (m + n)
π /2 Γ (m) Γ (n) 1
Hence ∫0 sin2 m − 1 θ ⋅ cos2 n − 1 θ dθ = = B (m, n)
2Γ (m + n) 2
3p + 1
Example 19: Using Heaviside’s expansion formula find L−1 2
( p − 1) ( p + 1)
(Gorakhpur 2008)
Solution: Here F ( p) = 3 p + 1
and G ( p) = ( p − 1) ( p2 + 1) = ( p − 1) ( p − i) ( p + i).
4et (3 i + 1) it (− 3 i + 1) − it
= + e + e
2 − (2 + 2 i) (− 2 + 2 i)
(3 i + 1) (1 − i) it (3 i − 1) (1 + i) − it
= 2e t − e + e
2 (1 + i)(1 − i) 2 (1 − i) (1 + i)
1 1
= 2et − (i + 2) e it + (i − 2) e − it
2 2
1
= 2 e t − i (e i t − e − i t ) − (e i t + e − i t )
2
1
= 2 e t − i . 2 i sin t − 2 cos t = 2 e t + sin t − 2 cos t.
2
Example 20: Using Heaviside’s expansion formula find
2 p2 + 5 p − 4
L−1 3 2 ⋅
p + p − 2 p (Avadh 2010)
Solution: Here F ( p) = 2 p2 + 5 p − 4
and G ( p) = p3 + p2 − 2 p = p ( p − 1) ( p + 2)
G ′ ( p) = 3 p2 + 2 p − 2.
G ( p) has 3 distinct zeros α1 = 0 , α2 = 1 and α3 = − 2.
∴ By the Heaviside’s expansion formula, we have
2 p2 + 5 p − 4 F (0 ) 0 t F (1) t F (− 2) −2 t
L−1 3 2 = e + e + e
p + p − 2 p G ′ (0 ) G ′ (1) G ′ (− 2)
= 2 + e t − e −2 t .
Comprehensive Exercise 4
1 1
3. (i) L−1 2 (ii) L−1 2 ⋅
( p − 2) ( p + 1) ( p + 4) ( p + 2)
1
4. (i) L−1 2 2 (Gorakhpur 2006)
p ( p + 1)
1
(ii) L−1 2 ⋅ (Kanpur 2011)
( p + 2) ( p − 2)
p2
5. Find L−1 2 2
⋅
( p + 4) (Kanpur 2010)
19 p + 37 p+5
9. (i) L−1 (ii) L−1 2 ⋅
( p + 1) ( p − 2) ( p + 3) ( p + 1) ( p + 1)
A nswers 4
1 t 1 2t
1. (i) (e − e − 2 t ) (ii) (e − e − t)
3 3
1 t
(iii) (e − e − 3 t )
4
1
2. (i) (1/ 2 a) t sin at (ii) (1 − t sin 2 t − cos 2 t)
16
1 2t 1 − 2t
3. (i) (e − 2 sin t − cos t) (ii) (e + sin 2 t − cos 2 t)
5 8
1 2t
4. (i) (t + 2) e − t + t − 2 (ii) [e − (4 t + 1) e 2 t ]
16
1 1 t 5
5. (2 t cos 2 t + sin 2 t). 7. e − e2t + e3t
4 2 2
5 − t 1 − 3t
8. −2+ e + e
2 2
9. (i) − 3 e − t + 5 e2 t − 2 e −3 t (ii) 2 e − t − 2 cos t + 3 sin t
T-67
e −√ p − x √ p
−1 e x
Example 21: Find L−1 and hence deduce that L = erfc ⋅
p p 2 √ t
Solution: Let f ( p) = e − √ p.
∴ F (t) = L−1 {e − √ p}
p p3 /2 p2 p5 /2
= L−1 1 − √ ( p) + − + − + …
2 ! 3 ! 4 ! 5 !
1 −1 1 −1 3 / 2
= L−1 {1} − L−1 { p1 /2} + L { p} − L {p }
2! 3!
1 −1 2 1 −1 5 / 2
+ L {p }− L { p } + ... …(1)
4! 5!
1
Now L−1{ pn + (1 /2)} = L−1 − n−(1 /2)
p
− n − (3 /2)
t
= , for n = 0 , 1, 2, ...
1
Γ (− n − )
2
(− 1)n + 1 1 3 5 2 n + 1 − n − (3 /2)
= ... t .
√ π 2 2 2 2
1 n+1 2 2 2 2
∵ Γ (− n − ) = (− 1) ... √ π
2 1 3 5 2 n + 1
e −√ p t 1
L−1 = ∫0 e −1 /(4 x)dx
p 2 √ (π) x3 /2
2 1 /(2 √ t) 2 1 dy
=− ∫∞ e − y dy, putting x = 2
so that dx = −
√π 4y 2 y3
2 ∞ − y2 1
= ∫1 / (2 √t) e dy = erfc ⋅
π 2 √ t
∵ Complementary Error Function erfc (t) = 1 − erf (t) = 1 − 2 t − x2
√π ∫0 e dx
2 ∞ − x2 0 2 2 ∞ − x2
e dx + e − x dx = e dx
√ π ∫0 ∫t ∫t
=
√π
e −√ p 1
Deduction: We have L−1 = erfc ⋅
p 2 √ t
e −√( x2 p) 1 1
∴ L−1 2 = 2 erfc 2
, by change of scale property
x p x 2 √ (t / x )
− x √p
e x ⋅
or L−1 = erfc
p 2 √ t
1 t2 t5 t8 t11
Example 22: (i) Prove that L−1 3 = − + − + ...
p + 1 2 ! 5 ! 8 ! 11!
(ii) Applying Heaviside’s Expansion formula , prove that
1 1 −t 1 1
L−1 3 t /2
= [e − e { cos ( √ 3 t) − √ 3 sin ( √ 3 t)}].
p + 1 3 2 2
1 1
Solution: (i) We have, 3 = 3 (1 + 1 / p3 )−1
p +1 p
1 1 1 1 1
= 3
1 − 3 + 6 − 9 + 12 − ...
p p p p p
1 1 1 1
= − + − + ...
p3 p6 p9 p12
1 t2 t5 t8 t11
∴ L−1 3 = − + − +…
p + 1 2 ! 5 ! 8 ! 11 !
(ii) Here F ( p) = 1 and G ( p) = p 3 + 1 = ( p + 1) ( p2 − p + 1), G ′ ( p) = 3 p2 .
1 1
G ( p) has 3 distinct zeros − 1, (1 + √ 3 i) and (1 − √ 3 i).
2 2
∴ By the Heaviside’s expansion formula , we have
T-69
1
1 F (− 1) − t F { 12 (1 + √ 3 i )} { (1 + √3 i)} t
L−1 3 = e + e 2
p + 1 G ′ (− 1) G ′ { 12 (1 + √ 3 i )}
1
F { 12 (1 − √ 3 i )} { (1 − √3 i)} t
+ e 2
G ′ { 12 (1 − √ 3 i )}
1 1
1 −t 1 { (1 + √3 i )}t 1 { (1 − √3 i)}t
= e + e2 + e2
3 3
4
(1 + √ 3 i )2 3
4
(1 − √ 3 i )2
1
1 2 (√ 3 i + 1) { (1 + √3 i )}t
= e− t + e 2
3 3 (√ 3 i − 1) (√ 3 i + 1)
1
2 (√ 3 i − 1) { (1 − √3 i )}t
+ e2
3 (√ 3 i + 1) (√ 3 i − 1)
1 − t 1 t /2
e − e
= [(√ 3 i + 1) e √3 it /2 − (√ 3 i − 1) e −√3 it /2 ]
3 6
1 1
= e − t − e t /2 [√ 3 i (e √3 it /2 − e − √3 it /2 ) + (e √3 it /2 + e −√3 it /2 )]
3 6
1 1 1
= [e − t − e t /2 {− √ 3 sin ( √ 3 t) + cos ( √ 3 t)}].
3 2 2
∞ 1
Example 23: Show that ∫ cos x2 dx = √ (π / 2) (Avadh 2008)
0 2
∞
Solution: Let F (t) = ∫0 cos tx2 dx.
∞ − pt ∞ − pt ∞ cos tx2 dx dt
L { F (t)} = ∫0 e F (t) dt = ∫0 e
∫ 0
∴
∞ ∞ − pt ∞ ∞ p
=∫ e cos tx2 dt dx = ∫ L {cos tx2 } dx = dx
0 ∫ 0 0 ∫0 p + x4
2
1 π /2 dθ
= ∫ , putting x = √ ( p tan θ) so that
2 √ p 0 √ (tan θ)
p sec2θ dθ
dx =
2 √ ( p tan θ)
1 π /2
= ∫ sin−1 /2 θ cos 1 /2 θ dθ
2√ p 0
1 3 1 1
1 Γ (4 ) Γ(4 ) 1 Γ (4 ) Γ (1 − 4 ) 1 π π
= = = ⋅ 1
= ⋅
2 √ p 2 Γ (1) 2√ p 2 4 √ p sin 4 π 2 √ (2 p)
π
∵ Γ ( p) Γ (1 − p) = , 0 < p < 1
sin pπ
π 1 π t(1 /2) −1 1 π
∴ F (t) = L−1 1 /2 = 1
=
2√2 p 2 2 Γ (2 ) 2 2t
∞ 1 π
or ∫0 cos tx2 dx = .
2 2t
Now taking t = 1, we have
T-70
∞ 1
∫0 cos x2 dx = √ (π / 2)
2
∞ 2 1
Example 24: Prove that ∫0 e − x dx = √ π.
2 (Meerut 2013B)
∞ 2
Solution: Let F (t) = ∫0 e − tx dx.
π −1 1 π 1 1 π
or F (t) = L = ⋅ =
2 √ p 2 √ (π t) 2 t
∞ 2 1 π
or ∫0 e − t x dx = .
2 t
∞ 2 1
Taking t = 1, we have ∫0 e − x dx = √ π.
2
8
Example 25: Prove that L−1 2 3
= (3 − t2 ) sint − 3 t cos t.
( p + 1)
1
Solution: We have L−1 2 = sin t .
+ 1
p
∴ by the convolution theorem, we have
1 1 t
L−1 2 ⋅ 2 = ∫ 0 sin x sin (t − x) dx
( p + 1) (( p + 1)
t
=∫ sin x (sin t cos x − cos t sin x) dx
0
t t
= sin t ∫ sin x cos x dx − cos t ∫ sin2 x dx
0 0
1 t 1 t
= sin t ∫ sin 2 x dx − cos t ⋅ ∫ (1 − cos 2 x) dx
2 0 2 0
1 1 1 1
= sin t . (1 − cos 2 t) − cos t . (t − sin 2 t)
2 2 2 2
1 1 1
= sin t.sin2 t − t cos t + cos t sin t cos t
2 2 2
1 1 t
or L−1 2 2
= sin t − cos t.
( p + 1) 2 2
1 1 1
∴ L−1 2 3
= 8 L−1 2 2
⋅ 2
( p + 1) ( p + 1) ( p + 1)
t 1 x
= 8 ∫ sin x − cos x sin (t − x) dx, by the convolution theorem
0 2 2
t
= 4 ∫ (sin x − x cos x) (sin t cos x − cos t sin x) dx
0
t t
= 4 sin t ∫ (sin x cos x − x cos2 x) dx − 4 cos t ∫0 (sin2 x − x sin x cos x) dx
0
T-71
t t
= 2 sin t ∫ {sin 2 x − x (1 + cos 2 x) } dx − 2 cos t ∫ {(1 − cos 2 x) − x sin 2 x} dx
0 0
t2 1 − cos 2 t t 1 − cos 2 t
= 2 sin t − + − sin 2 t +
2 2 2 4
1 t cos 2 t sin 2 t
−2 cos t t − sin 2 t + −
2 2 4
3 3
= − t 2 sin t + sin t − sin t cos 2 t − t sin t sin 2 t
2 2
3
− 2 t cos t + sin 2 t cos t − t cos t cos 2 t
2
3 3
= − t 2 sin t + sin t + (− sin t cos 2 t + sin 2 t cos t)
2 2
− t (cos 2 t cos t + sin t sin 2 t) − 2 t cos t
2
= (3 − t ) sin t − 3 t cos t.
Comprehensive Exercise 5
1
1. Find L−1 5 .
( p − 1) ( p + 2) (Kanpur 2009, 11)
1 t
2. Prove that L−1 2 2
= ∫0 J0 (ax) dx .
p √ ( p + a )
3. Apply Heaviside’s expansion formula, to prove that
1 1 t 1 1
L−1 3 = [e − e
−1 / 2
{cos ( √ 3 t) + √ 3 sin ( √ 3 t)}]
p − 1 3 2 2
∞ 1
4. Show that ∫ sin x2 dx = √ (π /2).
0 2
5. Apply the convolution theorem to show that
t t
∫ 0 sin u cos (t − u) du = 2 sin t . (Gorakhpur 2007, 11)
∞ sin2 x π
6. Use Laplace transform to prove that ∫0 dx = ⋅
x2 2
2 2
7. Prove that L−1 tan−1 2 = sin t sinh t. (Agra 2002)
p t
1
9. Find L−1 2 2 3 /2
⋅
( p + a )
t t t t (t − u)n−1
10. Prove that ∫0 ∫0 ...... ∫0 F (t) dt n = ∫0 F (u) du.
(n − 1) !
A nswers 5
1 t 4 4 3 4 2 8 8 1 −2 t
1. e t − t + t − t + − e
72 3 3 9 27 243
t
9. J1 (at)
a
5. Differentiation f n
( p) (− 1)n t n
F (t)
theorem
Integral theorem
∞ 1
6. ∫p f ( x) dx F (t)
t
7. Multiplication p f ( p) F ′ (t)
theorems
pn f ( p ) F n (t)
T-73
f ( p) t
∫0 F ( x) dx
p
8. Division theorems t t t
f ( p) F (t) (dt)n
∫0 ∫0 …∫0
pn
9. Convolution f ( p). g ( p ) t
theorem F*G= ∫0 F ( x) ⋅ G (t − x) dx
F ( p) n F (α r ) α r t
10. Heaviside’s , Σ e
G ( p) r =1 G ′ (α r )
expansion
theorem degree F ( p ) < degree where α r , r = 1, 2, ... n, are
G ( p) roots of
G ( p ) = 0 and are all
distinct.
(a) π (b) t
t π
1
(c) (d) √ (πt).
√ (πt)
1
3. The value of L−1 5 /2 is
p
4 3 /2 4 5 /2
(a) t (b) t [
3 3
4 t 4 π
(c) t (d) t ⋅
3 π 3 t
T-74
L−1 { f ( p ) g ( p )} = …… .
∞
5. The value of ∫ cos x2 dx = ......
0
True or False
Write ‘T’ for true and ‘F’ for false statement.
4
1. The value of L−1 − 4t
is 3 e .
p + 4
1
2. The value of L−1 2
is t2 e − t .
( p + 1)
1
3. The value of L−1 2 − 4t
is te .
p + 8 p + 16
4. The convolution of two functions F and G obeys the commutative law.
1 1 −t
5. The value of L−1 2t
is (e − e ).
( p − 1) ( p + 2) 3
A nswers
True or False
1. T 2. F 3. T
4 T 5. F
¨
T-77
3
A pplications of L aplace
T ransform
(To Solutions of Differential Equations and Integral Equations)
be the given initial or boundary conditions where A0 , A1, A2 ,..., An−1 are constants.
On taking the Laplace transform of both sides of equation (1) and using conditions (2)
we obtain an algebraic equation known as “subsidiary equation” from which
y ( p ) = L { y(t)} is determined. The required solution is then obtained by finding the
inverse Laplace transform of y ( p ).
T-78
Solution: (a) Taking Laplace transform of both sides of the given equation, we have
L { x ′ ′ } + m2 L { x } = a L {cos nt}
ap
or p2 L { x} − px (0 ) − x ′ (0 ) + m2 L { x} =
p + n2
2
ap
or ( p2 + m2 ) L {x} = px0 + x1 +
p2 + n2
p 1 a. p
or L { x} = x0 ⋅ + x1 ⋅ +
2 2 2 2
p +m p +m ( p + m2 ) ( p2 + n2 )
2
p 1 ap ( p2 + m2 ) − ( p2 + n2 )
= x0 . + x1. +
p2 + m2 p2 + m2 m2 − n2 ( p2 + m2 ) ( p2 + n2 )
p 1 a p p
= x0 . + x1 . + ⋅ − ⋅
p2 + m2 p2 + m2 (m2 − n2 ) p2 + n2 p2 + m2
Taking the inverse transform, we have
p 1
x = x0 . L−1 + x1 . L−1
2 2 2 2
p +m p +m
a p p
+ ⋅ L−1 − L−1
2
(m2 − n2 ) 2
p +n
2 2
p + m
x a
or x = x0 cos mt + 1 ⋅ sin mt + (cos nt − cos mt),
m (m − n2 )
2
(Rohilkhand 2000, 02, 04; Agra 01; Gorakhpur 07, 10, 11; Kashi 14)
Solution: Taking the Laplace transform of both sides of the given equation, we have
L { y ′ ′ } + 9 L { y } = L {cos 2 t}
or p2 L { y } − p y(0 ) − y ′ (0 ) + 9 L { y } = p /( p2 + 4)
or ( p2 + 9) L { y } − p − A = p / ( p2 + 4), where y ′ (0 ) = A
p+ A p
or L{ y} = +
p2 + 9 ( p2 + 9) ( p2 + 4)
p A p p
= + + − ⋅
p2 + 9 p2 + 9 5 ( p2 + 4) 5 ( p2 + 9)
p 1 1 −1 p 1 −1 p
∴ y = L−1 −1
+ AL 2 + L 2 − L 2
2
p + 9 p + 9 5 p + 4 5 p + 9
1 1 1
= cos 3 t + A sin 3 t + cos 2 t − cos 3 t
3 5 5
4 1 1
= cos 3 t + A sin 3 t + cos 2 t.
5 3 5
But y (π / 2) = − 1.
4 3 1 3 1
∴ − 1 = cos π + A sin π + cos π
5 2 3 2 5
1 1
or − 1= − A − or A = 12 / 5.
3 5
4 4 1
Hence the required solution is y = cos 3 t + sin 3 t + cos 2 t.
5 5 5
Example 4: Solve ( D + 1)2 y = t given that y = − 3, when t = 0 and y = − 1, when t = 1.
(Rohilkhand 2006)
2
Solution: The given equation can be written as ( D + 2 D + 1) y = t.
∴ L { y ′ ′ } + 2 L { y ′ } + L { y} = L { t}
or p2 L { y} − py (0 ) − y ′ (0 ) + 2 [ pL { y} − y (0 )] + L { y} = 1/ p2
or ( p2 + 2 p + 1) L { y } − p (− 3) − A − 2 (− 3) = 1/ p2 where y ′ (0 ) = A
or ( p + 1)2 L { y} = (1 / p2 ) − 3 p − 6 + A
1 3p + 6 A
or L{ y} = − +
2 2 2
p ( p + 1) ( p + 1) ( p + 1)2
1 3 ( p + 1) + 3 A
= − +
2 2 2
p ( p + 1) ( p + 1) ( p + 1)2
2 1 2 1 3 3 A
=− + + + − − +
p p2 p + 1 ( p + 1)2 ( p + 1) ( p + 1)2 ( p + 1)2
2 1 1 A−2
=− + − + ⋅
p p 2 ( p + 1) ( p + 1)2
T-80
1 1 1 1
∴ y = − 2 L−1 + L−1 − L−1 + ( A − 2) L−1
2 2
p p p + 1 ( p + 1)
. + t − e − t + ( A − 2) e − t L−1 {1/ p2} = − 2 + t − e − t + ( A − 2) te − t .
= − 21
Now, since y = − 1, when t = 1,
∴ − 1 = − 2 + 1 − e − 1 + ( A − 2) e − 1 or A = 3.
Hence the complete solution is y = − 2 + t − e − t + te − t .
dy
Example 5: Solve ( D2 + 1) y = t cos 2 t, y = 0 , = 0 when t = 0.
dt (Avadh 2006, 07; Meerut 13B)
Solution: Taking the Laplace transform of both the sides of the given equation, we
have L { y ′ ′ } + L { y } = L { t cos 2 t}
d
or p2 L { y } − p y (0 ) − y ′ (0 ) + L { y } = − ( L {cos 2 t})
dp
d p 1 2 p2
or ( p2 + 1) L { y} = − =− +
dp p2 + 4 p2 + 4 ( p2 + 4)2
p2 − 4 5 5 8
or L{ y} = =− + +
( p2 + 1) ( p2 + 4)2 9 ( p2 + 1) 9 ( p2 + 4) 3 ( p2 + 4)2
5 −1 1 5 −1 1 8 −1 1
∴ y=− L + L 2 + L 2
9 2 2
p + 1 9 p + 4 3 ( p + 4)
5 5 8 t1 1
= − sin t + sin 2 t + ∫ sin 2 x ⋅ sin 2 (t − x) dx ,
9 18 3 0 2 2
1 1
−1
By the convolution theorem since L 2 = sin 2t
p + 4 2
5 5 1 t
= − sin t + sin 2 t + {cos (2 t − 4 x) − cos 2 t} dx
9 18 3 ∫0
t
5 5 1 − 1 sin (2 t − 4 x) − x cos 2 t
=− sin t + sin 2 t +
9 18 3 4 0
5 5 1 1 1
=− sin t + sin 2 t + sin 2 t − t cos 2 t + sin 2 t
9 18 12 3 12
5 4 1
or y = − sin t + sin 2 t − t cos 2 t , which is the required solution.
9 9 3
Example 6: Solve ( D3 − D2 − D + 1) y = 8 t e − t if y = D2 y = 0 , Dy = 1 when t = 0.
Solution: Taking the Laplace transform of both sides of the given equation, we have
L { y ′ ′ ′ } − L { y ′ ′ } − L { y ′ } + L { y } = 8 L { te − t}
or p3 L { y } − p2 y (0 ) − py ′ (0 ) − y ′ ′ (0 ) − [ p2 L { y } − p y (0 ) − y ′ (0 )]
d
− [ pL { y } − y (0 )] + L { y } = − 8 [ L {e − t}]
dp
T-81
d 1
or ( p3 − p2 − p + 1) L { y } − p + 1 = − 8
dp p + 1
8
or ( p − 1)2 ( p + 1) L { y } = p − 1 +
( p + 1)2
1 8
or L{ y} = +
( p − 1) ( p + 1) ( p − 1)2 ( p + 1)3
1 1 1 3 1 3
= − − + +
2 p − 1 p + 1 2 ( p − 1) ( p − 1)2 2 ( p + 1)
2 2
+ +
( p + 1)2 ( p + 1)3
1 1 1 2 2
=− + + + + ⋅
2 2
p − 1 p + 1 ( p − 1) ( p + 1) ( p + 1)3
1 −1 1
1
∴ y = − L−1 +L
−1
+L 2
p − 1 p + 1 ( p − 1)
1 1
+ 2 L−1 + 2 L−1
2 3
( p + 1) ( p + 1)
1 1 1
= − e t + e − t + e t L−1 + 2 e − t L−1 + 2 e − t L−1
2 2 3
p p p
= − e t + e − t + e t t + 2 e − t t + 2 e − t (t 2 / 2 !)
= (1 + 2 t + t 2 ) e − t − (1 − t) e t , which is the required solution.
y ′ ′ ′ (0 ) = − 3.
Solution: Taking the Laplace transform of both sides of the given equation, we have
L { y iv} + 2 L { y ′ ′ } + L { y } = 0
or p4 L { y } − p3 y (0 ) − p2 y ′ (0 ) − p y ′ ′ (0 ) − y ′ ′ ′ (0 )
+ 2 [ p2 L { y } − p y (0 ) − y ′ (0 )] + L { y } = 0
or ( p4 + 2 p2 + 1) L { y } − p2 − 2 p + 3 + 2 (− 1) = 0
or ( p2 + 1)2 L { y } = p2 + 2 p − 1
p2 + 2 p − 1 1 2p − 2 1 2p 2
or L{ y} = = + = + − ⋅
2 2 2 2 2 2 2 2
( p + 1) p +1 ( p + 1) p +1 ( p + 1) ( p + 1)2
2
1 2p 1
∴ y = L−1 −1
+L 2 − 2 L−1 ⋅ …(1)
2 2 2 2
p + 1 ( p + 1) ( p + 1)
T-82
2p d 1 1
Now L−1 = − L−1 = − (− t) L−1 = t sin t
2 2 2 2
( p + 1) dp p + 1 p + 1
1 1 1
and L−1 = L−1 ⋅ = F (t) * F (t),
2 2 2 2
( p + 1) p + 1 p + 1
1
where F (t) = L−1 = sin t
2
p + 1
t 1 t
= ∫ sin x ⋅ sin (t − x) dx = ∫ [cos (2 x − t) − cos t] dx
0 2 0
t
1 1 1
= sin (2 x − t) − x cos t = (sin t − t cos t).
2 2 0 2
1
∴ x = L−1 f ( p) = (sin t) * F (t)
2
p +1
t
or x= ∫0 sin (t − x) ⋅ F ( x) dx , (by the convolution theorem)
Comprehensive Exercise 1
12. ( D3 + D) y = e2 t , y (0 ) = y ′ (0 ) = y ′ ′ (0 ) = 0 .
13. ( D3 − 2 D2 + 5 D) y = 0 , if y (0 ) = 0 , y ′ (0 ) = 1, y (π / 8) = 1.
14. ( D2 + D) y = t2 + 2 t, where y (0 ) = 4, y′ (0 ) = − 2.
T-84
A nswers 1
1. (i) y = e − t + 1 (ii) y = cos t
2. (i) y = 5 cos t + sin t − 2 cos 2 t (ii) y = 2 + 2 t + e − t
x a sin nt − n sin mt
3. (i) x = x0 cos mt + 1 sin mt +
m 2 2 m
m −n
(ii) y = e −2 t (2 t2 + 2 t − 1)
1 4 −2 t 8 3 t
4. (i) y = − + e + e
3 5 15
1
(ii) y = [(53 + 155 x) e −3 x − (3 cos x − 4 sin x)]
50
5. y = π sin 3 t + 2 t
3
6. y = 4 e3 t − e −3 t (4 cos 4 t + 7 sin 4 t)
4
7. y = 25 t2 + 40 t + 22 + 2 e2 t (2 sin t − 11 cos t)
1 −t 3 1 1 2t
8. y= e . t + 4 e − t + 6 te − t 9. y= + e − te2 t
2 2 2
15 1 1
10. y= cos t + t sin t + cos 3 t
16 4 16
11. y = 3 sinh t − sin t + cosh t + 2
1 1 2t 2 1
12. y=− + e + cos t − sin t
2 10 5 5
1 3
13. y = 1 + e t (sin 2 t − cos 2 t) 14. y = t + 2 (1 + e − t )
3
Solution: Taking the Laplace transform of both sides of the given equation, we have
L { t y ′ ′ } + L { y ′ } + 4 L { t y} = 0
T-85
d d
or − L{ y ′ ′ } + L{ y ′ } + 4. (− 1) L{ y } = 0
dp dp
d 2 d
or − [ p y − py (0 ) − y ′ (0 )] + [ p y − y(0 )] − 4 y = 0,
dp dp
where y = L{ y }
d d y
or − ( p2 y − 3 p) + ( p y − 3) − 4 =0
dp dp
d y d y p
or − ( p2 + 4) − p y = 0 or + dp = 0.
dp y 2
p +4
1 C1
Integrating, log y + log ( p2 + 4) = log C1 or y= ⋅
2 √ ( p2 + 4)
1
∴ y = L−1 { y } = C1 L−1
2
√ ( p + 4)
or y = C1 J0 (2 t) (See Ex. 25, after article 1.23 of Chapter 1)
Since y(0 ) = 3, therefore 3 = C1 J0 (0 ) = C1. [ ∵ J0 (0 ) = 1]
Hence y = 3 J0 (2 t), which is the required solution.
Example 11: Solve [t D2 + (1 − 2 t) D − 2] y = 0 if y(0 ) = 1, y ′ (0 ) = 2.
d 1
or p2 y − py (0 ) − y ′ (0 ) + [ L { y ′ }] + y =
dp p
d 1
or p2 y − p − 2 + [ p y − y(0 )] + y =
dp p
d 1
or p2 y − p − 2 + ( p y − 1) + y =
dp p
d y 1
or p + ( p2 + 2) y = p + 2 +
dp p
d y 2 2 1
or + p+ y = 1+ + 2
dp p p p …(1)
c p2 p4 1 2
= 1 1 − + −... + +
p2 2 4 ⋅2 ! p p2
(2 + c1) c c 1
= − 1 + 1 p2 − ... + ⋅
2 2 8 p
p
1 1
∴ y = (2 + c1) L−1 {1/ p2} − c1 L−1 {1} + c1 L−1 { p2} + ... + L−1 {1/ p}
2 8
= (2 + c1) t + 1, since L−1 { pn} = 0 , n = 0 , 1, 2, ...
But given y ′ (0 ) = 2. Therefore,
2 = 2 + c1 or c1 = 0 .
Hence y = 2 t + 1, which is the required solution.
T-87
Comprehensive Exercise 2
A nswers 2
1. y=t 2. y = 5e− t
1 sin t
3. y = t2 4. y =
2 t
at
5. y = A + e − t
3
(Rohilkhand 2010)
Solution: Taking the Laplace transform of both sides of the two equations, we have
L { x ′ } + L { y ′ } = L { t} and L { x ′ ′ } − L { y} = L { e − t}
or p x − x (0 ) + p y − y (0 ) = 1 / p2
and p2 x − p x (0 ) − x ′ (0 ) − y = 1 / ( p + 1)
or p x + p y = 3 + 1/ p2 and p2 x − y = 3 p − 2 + 1/( p + 1).
Solving for x and y, we have
3 p2 + 1 3p 2 1
x= + − −
3 2 2 2
p (1 + p ) 1+ p 1+ p ( p + 1) ( p2 + 1)
T-88
1 2 p4 3p 2 1 p 1
= 1 + 2 p2 − + − + − +
3 2 2 2
p 1+ p 1+ p 1+ p 2 ( p + 1) 2 ( p + 1) 2 ( p2 + 1)
2
2 1 1 p 3
= + + + −
p p 3 2 ( p + 1) 2 ( 1 + p ) 2 ( p2 + 1)
2
1 2
and y= +
2 2
p ( p + 1) ( p + 1) p +1
1 1 p 1 2
= − − − +
2 2 2
p 2 ( p + 1) 2 ( p + 1) 2 ( p + 1) p + 1
1 1 p 3
= − − + ⋅
p 2 ( p + 1) 2 ( p2 + 1) 2 ( p2 + 1)
1 1 1 −1 1
∴ x = 2 L−1 + L−1 + L
3
p p 2 p + 1
1 −1 p 3 −1 1
+ L − L 2
2 2
p + 1 2 p + 1
1 2 1 −t 1 3
=2+ t + e + cos t − sin t
2 2 2 2
1 1 1 1 −1 p 3 −1 1
and y = L−1 − L−1 − L 2 + L 2
p 2 p + 1 2 p + 1 2 p + 1
1 1 3
= 1 − e − t − cos t + sin t.
2 2 2
Example 14: Solve Dx + 2 D 2 y = e − t , ( D + 2) x − y = 1 if x (0 ) = y (0 ) = y ′ (0 ) = 0 .
Solution: Taking the Laplace transform of both sides of the given equations, we have
L { x ′ } + 2 L { y ′ ′ } = L { e− t }
and L { x ′ } + 2 L { x} − L { y} = L {1}
1
or p x − x (0 ) + 2 [ p2 y − py (0 ) − y ′ (0 )] = ,
p+1
where x = L { x } and y = L { y}
1
and p x − x (0 ) + 2 x − y =
p
1 1
or p x + 2 p2 y = and ( p + 2) x − y = ⋅
p+1 p
Solving for x and y, we have
1 2
x= +
p ( p + 1) (2 p2 + 4 p + 1) (2 p2 + 4 p + 1)
1 1
= +
2 − √ 2 2 + √ 2 2 − √ 2 2 + √ 2
2 p ( p + 1) p + p+ p+ p+
2 2 2 2
T-89
1 1 1 1
= + − −
p p+1 2 − √2 2 + √ 2
(2 − √ 2) p + (2 + √ 2) p +
2 2
1 1
+ −
2 − √ 2 2 + √ 2
√2 p+ √2 p+
2 2
1 1 1 1
= + − −
p p+1 2 − √ 2 2 + √ 2
p+ p+
2 2
1
and y=
p ( p + 1) (2 p2 + 4 p + 1)
1 1 1 1
= + − −
p p+1 2 − √ 2 2 + √ 2
(2 − √ 2) p + (2 + √ 2) p +
2 2
1 1 2 + √ 2 1 2 − √ 2 1
= + − −
p p+1 2 2 − √2 2 2 + √ 2
p+ p+
2 2
1 1 −1 1 −1 1
∴ x = L−1 + L−1 −L −L
p p + 1 p + a p + b
= 1 + e − t − e − at − e − b t
1 1 −1 1 −1 1
and y = L−1 + L−1 −bL − aL
p p + 1 p + a p + b
2 − √2 2 + √2
= 1 + e − t − be − a t − ae − b t , where a= and b = ⋅
2 2
Example 15: Solve ( D2 − 3) x − 4 y = 0 , x + ( D2 + 1) y = 0, t > 0
if x = y = D y = 0 , D x = 2 when t = 0.
Solution: Taking the Laplace transform of both sides of the two equations, we have
L { x ′ ′ } − 3 L { x } − 4 L { y } = 0 and L { x } + L { y ′ ′ } + L { y } = 0
or p2 x − px (0 ) − x ′ (0 ) − 3 x − 4 y = 0 , where x = L { x }
and y = L{ y}
2
and x + p y − py (0 ) − y ′ (0 ) + y = 0
or ( p2 − 3) x − 4 y = 2 and x + ( p2 + 1) y = 0 .
−2 1 1 1 1 1
and y= = − + − − ⋅
( p + 1)2 ( p − 1)2 2 p + 1 p − 1 ( p + 1)
2
( p − 1)2
1 1
∴ x = L−1 + L−1
2 2
( p − 1) ( p + 1)
= e t L−1 {1/ p2} + e − t L−1 {1/ p2} = (e t + e − t ) t
1 −1 1 1 1 1
and y= − L + L−1 −1
−L
−1
−L
2 p + 1 2 2
p − 1 ( p + 1) ( p − 1)
1 1 1
= [− e − t + e t − t e − t
− te t ] = (1 − t) e t − (1 + t) e − t .
2 2 2
Comprehensive Exercise 3
3. ( D2 − 1) x + 5 Dy = t, − 2 Dx + ( D2 − 4) y = − 2 ,
if x = 0 = Dx = y = D y when t = 0.
4. Applying Laplace transform solve the equations :
dx dy d2 x
+ = t and − y = e− t,
dt dt dt2
dx
given that x (0 ) = 0 = y (0 ) and = 0 when t = 0. (Rohilkhand 2007)
dt
5. Solve ( D − 2) x − ( D + 1) y = 6 e3 t , (2 D − 3) x + ( D − 3) y = 6 e3 t , if x = 3, y = 0
when t = 0.
6. Solve ( D − 2) x − ( D − 2) y = sin t, ( D2 + 1) x + 2 D y = 0 ,
if x = 0 = x ′ (0 ) = y (0 ).
A nswers 3
1. x = 5 e − t + 3 e4 t , y = 5 e − t − 2 e4 t .
1 1 1
2. x= − (2 cos t + cos 2 t), y = (sin 2 t − 2 sin t).
2 6 6
3. x = − t + 5 sin t − 2 sin 2 t, y = 1 − 2 cos t + cos 2 t.
1 1 1 1 1 1 1
4. x = t 2 − 1 + e − t + cos t + sin t, y = 1 − e − t − cos t − sin t.
2 2 2 2 2 2 2
T-91
5. x = e t + 2 t e t + 2 e 3 t, y = e t − t e t − 2 e 3 t.
1 4 2t 1 1
6. x= (1 + 3 t) + e − t + e − (cos t + 2 sin t), y = [(1 + 3 t) e − t − e2 t ].
9 45 5 9
∂y d y ∂2 y d 2 y
(c) L = and (d) L = where L { y ( x, t)} = y ( x, p).
2 2
∂x dx ∂ x dx
∂y ∞ ∂y s ∂y
Proof: (a) L = ∫0 e − pt dt = lim ∫0 e − pt dt
∂t ∂t s→ ∞ ∂t
s s
= lim e − pt y ( x, t)
{ + p ∫ e − pt y ( x, t) dt
}
s→ ∞ 0 0
∞
=p ∫0 e − pt y ( x, t) dt − y ( x, 0 ) = py ( x, p) − y ( x, 0 )
∂y
(b) Let V = ⋅
∂t
∂2 y
∂V
∴ L
2 = L = p . L {V } − V ( x, 0 )
∂t ∂t
= p [ pL { y} − y ( x, 0 )] − yt ( x, 0 ) ∂y
∵ V = = yt
∂ t
= p2 y ( x, p) − py ( x, 0 ) − yt ( x, 0 ).
∂y ∞ ∂y d ∞ − pt d y
(c) L = ∫0 e − pt dt = ∫ e y dt = ⋅
∂x ∂x dx 0 dx
∂2 y ∂y
∂U
(d) L = L , where U =
2
∂x ∂x ∂x
T-92
2
d d ∂y d d y d y
= L {U } = L = = ⋅
dx dx ∂x dx dx dx2
∂y ∂2 y
Example 16: Solve =2 where y (0 , t) = 0 = y (5, t) and y ( x, 0 ) = 10 sin 4 πx.
∂t ∂x 2
(Meerut 2008)
Solution: Taking the Laplace transform of both sides of given equation, we have
∂y ∂2 y d2 y
L = 2L or p y − y ( x, 0 ) = 2
2
∂t ∂x dx2
d2 y p
or − y = − 5 sin 4 πx whose general solution is
2 2
dx
5 sin 4 πx
y = C1e √( p /2 ) x + C2 e − √( p /2 ) x −
− (4 π)2 − p / 2
10
or y = C1e √( p /2 ) x + C2 e − √( p /2 ) x + sin 4 πx …(1)
32 π2 + p
= C1e5 √( p /2 ) + C2 e −5 √( p /2 ) + 0 .
Solving C1 = 0 = C2 .
10
∴ from (1), we have y= sin 4 πx.
32 π2 + p
10 2
∴ y = L−1 sin 4 πx or y = 10 e −32 π t sin 4 πx,
2
32 π + p
which is the required solution.
∂y ∂2 y
Example 17: Find the bounded solution of = , x > 0, t > 0,
∂t ∂x2
where y (0 , t) = 1, y ( x, 0 ) = 0 . (Meerut 2007)
Solution: Taking the Laplace transform of both the sides of the given equation, we
have
∂y ∂2 y d2 y
L = L or p y ( x, p) − y ( x, 0 ) =
2
∂t ∂x dx2
d2 y
or − p y =0
dx2
T-93
Comprehensive Exercise 4
(Meerut 2006)
T-94
∂y ∂2 y
3. = 2 , y ( x, 0 ) = 3 sin 2 π x , y (0 , t) = 0 = y (1, t), 0 < x < 1, t > 0 .
∂t ∂x
(Meerut 2006)
2
∂y ∂ y
4. = 2 2 , y (0 , t) = 0 , y (5, t) = 0 , y ( x, 0 ) = 10 sin 4 πx − 5 sin 6 πx.
∂t ∂x
(Meerut 2008)
∂u ∂u
5. − = 1 − e − t , 0 < x < 1, t > 0 , u ( x, 0 ) = x.
∂x ∂t (Meerut 2007)
A nswers 4
1. y ( x, t) = 6 e − 2 t − 3 x . 2. y = 30 e − 75 t cos 5 x.
2
3. y ( x, t) = 3 e − 4 π t
sin 2 πx.
2 2
4. y ( x, t) = − 5 e − 72 π t
sin 6 πx + 10 e − 32 π t
sin 4 πx.
5. u ( x, t) = ( x + 1) − e − t .
is called an “Integral equation” where y (t) and K (u, t) are known, a and b are either
constants or functions of t.
Here the function F (t) which appears under the sign of integration is to be determined.
2. Abel’s integral equation: An equation of the form
t F (u) du
G (t) = ∫
0 (t − u)n
is called Abel’s integral equation where F (u) is unknown and G (t) is known and n is a
constant between 0 and 1, i. e., 0 < n < 1.
3. Integral equation of Convolution type: An integral equation of the form
t
F (t) = y (t) + ∫ 0 K (t − u) F (u) du
which may also be expressed as
F (t) = y (t) + K (t) * F (t)
is called an integral equation of convolution type.
4. Integro-differential equation: An integral equation in which various derivatives of
the unknown function F (t) can also be present is called an Integro-differential equation.
t
For example, F ′ (t) = F (t) + y (t) + ∫ 0 sin (t − u) F (u) du.
T-95
1 1 p2 +1
or 1 − 2 L { F (t)} = or L { F (t)} = ⋅
p + 1 p p3
p 2 + 1 1 1
∴ F (t) = L−1 3 = L−1 + L−1 3
p p p
t2 t2
= 1+ = 1+ ⋅
Γ (3) 2
t2
Verification: We have F (t) = 1 + ⋅
2
Putting in the R.H.S. of the given equation, we have
t u2
R.H.S. = 1 + ∫ 1 + ⋅ sin (t − u) du
0 2
t
u2 t
= 1 + 1 + cos (t − u) − ∫0 u ⋅ cos (t − u) du
2 0
t2 t
= 1+ 1+ − cos t − [− u ⋅ sin (t − u)]0t − ∫0 sin (t − u) du
2
t2 t t2
=2+ − cos t − [cos (t − u)]0 = 2 + − cos t − (1 − cos t)
2 2
t2
= 1+ = F (t) = L.H. S.
2
Example 21: Solve the following equation for F (t) with the condition that F (0 ) = 0 ,
t
F ′ (t) = sin t + ∫0 F (t − u) cos u du .
1
∴ F (t) = ± 8 L−1 2 2
= ± 8 J0 (4 t).
√ ( p + 4 )
(See Q. 2 of Ex ercise 5 Chapter 2)
t F (u) du
Example 23: Solve the integral equation ∫0 = t (1 + t).
(t − u)1 /3
t
Example 24: Solve F ′ (t) = t + ∫0 F (t − u) cos u du, F (0 ) = 4.
Putting the values of F ′ ′ (t), F ′ (t) and F (t) from (3), (5) and (6) in (1), we get
t t
2 G (t) − 3 ∫0 G (u) du + 3 − 2 ∫ (t − u) G (u) du + 2 t − 8
0
= 4 e − t + 2 cos t
t
or 2 G (t) + ∫0 (− 2 t + 2 u − 3) G (u) du = 4 e − t + 2 cos t − 2 t + 5 ,
= 4 (1 − e − t ) + 2 sin t.
t
or 2 F ′ (t) − 3 F (t) − 2 ∫0 F (u) du = − 4 e − t + 2 sin t − 10 . [Using (2)]
= 4 (e − t − 1) − 2 (cos t − 1) − 10 t.
t
or 2 F (t) + ∫0 (− 2 t + 2 u − 3) F (u) du = 4 e − t − 2 cos t − 10 t + 6
1
or C1 = − ∫0 (1 − u) G (u) du.
t 1
= ∫0 (t − u) G (u) du + ∫0 (ut − t) G (u) du
t t 1
= ∫0 (t − u) G (u) du + ∫ 0 (ut − t) G (u) du + ∫ t (ut − t) G (u) du
t 1
= ∫0 (t − 1) u G (u) du + ∫t t (u − 1) G (u) du
1 (t − 1) u , if u < t
or F (t) = ∫0 K (t, u) G (u) du, where K (t, u) =
t (u − 1) , if u > t.
Hence from (1), we get
t (t − 1) u , u < t
G (t) + a ∫0 K (t, u) G (u) du = 0 , where K (t, u) =
(u − 1) t , u > t.
This is the required integral equation of the given differential equation.
Method 2: Integrating (1) between the limits 0 to t, we get
t t
∫0 F ′ ′ (u) du + a ∫0 F (u) du = 0
t
or [ F ′ (u)] 0t a + ∫ 0 F (u) du = 0
t
or F ′ (t) − F ′ (0 ) + a ∫ 0 F (u) du = 0.
Again integrating between the same limits, we get
t
[ F (u)] 0t − F ′ (0) [u]0t + a ∫ 0 (t − u) F (u) du = 0
t
or F (t) − F (0 ) − tF ′ (0 ) + a ∫ 0 (t − u) F (u) du = 0 …(1)
t
or F (t) − t F ′ (0 ) + a ∫ 0 (t − u) F (u) du = 0
For t = 1, this gives
1
F (1) − F ′ (0 ) + a ∫ 0 (1 − u) F (u) du = 0 [Using (2)]
T-101
1
or 0 − F ′ (0 ) + a ∫0 (1 − u) F (u) du = 0
1
or F ′ (0 ) = a ∫ 0 (1 − u) F (u) du.
Putting this value in (I), we get
1 t
F (t) − ta ∫0 (1 − u) F (u) du + a ∫0 (t − u) F (u) du = 0
1 t
or F (t) + a ∫ 0 t (u − 1) F (u) du + a ∫ 0 (t − u) F (u) du = 0
t 1
or F (t) + a ∫ 0 t (u − 1) F (u) du + a ∫ t t (u − 1) F (u) du
t
+a ∫0 (t − u) F (u) du = 0
t 1
or F (t) + a ∫ 0 u (t − 1) F (u) du + a ∫ t t (u − 1) F (u) du = 0
1 u (t − 1) , u< t
or F (t) + a ∫0 K (u, t) F (u) du = 0 , where K (u, t) =
t (u − 1) , u > t.
Example 27: Convert the given integral equation
t
F (t) − ∫0 (t − u) sec t ⋅ F (u) du = t
We know that,
d b (t) b (t) ∂K db da
∫ a (t) K (u, t) du = ∫ a (t) du + K (b, t) − K (a, t) ⋅ …(2)
dt ∂t dt dt
Differentiating (1) using the fact (2), we get
t
F ′ (t) − ∫ 0 [sec t + (t − u) sec t tan t] F (u) du = 1
t t
or F ′ (t) − ∫ 0 sec t ⋅ F (u) du − tan t ∫ 0 (t − u) sec t F (u) du = 1
t
or F ′ (t) − ∫ 0 sec t ⋅ F (u) du + tan t ⋅ [t − F (t)] = 1, [Using (1)] …(3)
Comprehensive Exercise 5
1 t
5. Solve y (t) = t + ∫0 (t − u)3 ⋅ y (u) du.
6
1 2 t
6. Solve the integral equation F (t) = t − ∫0 (t − u) F (u) du. (Agra 2003)
2
t
7. Solve 2 F (t) = 2 − t + ∫0 F (t − u) F (u) du.
t F (u) du 2
8. Solve ∫ 0 √ (t − u) = 1 + t + t .
A nswers 5
1 2 1 3
1. F (t) = at e − t 2. y (t) = − 1 + t − t + t
2 6
1
4. F (t) = 1 + 2 t 5. y (t) = (sinh t + sin t)
2
T-103
t 5 t4 t3
or F (t) + ∫ 0 (− 8t + 8u + 2) F (u) du = 12 − − t−2
2
10. F ′ ′ ′ (t) + 6 F (t) = 0 with F (0 ) = 4, F ′ (0 ) = − 3, F ′ ′ (0 ) = 2
∂ y
3. If y ( x, t) is a function of x and t, then L = …… .
∂t
2
∂ y
4. If y ( x, t) is a function of x and t, then L 2 = ……
∂x
b
5. An equation of the form F (t) = y (t) + ∫a K (u, t) F (u) du is called …… .
equation where y (t) and K (u, t) are known, a and b are either constants or
functions of t.
True or False
Write ‘T’ for true and ‘F’ for false statement.
1. If y ( x, t) is a function of x and t then
∂2 y
L 2 = y ( x, p) − py ( x, 0 ) − yt ( x, 0 ) where L { y ( x, t)} = y ( x, p).
∂t
∂y ∂2 y
2. The differential equation = 2 2 with the condition y ( x, 0 ) = 10 sin 4 π x is
∂t ∂x
d2 y p
takes the form − y = − 5 sin 4 π x, where L { y ( x, t)} = y ( x, p).
dx2 2
t F (u)
3. An equation of the form G (t) = ∫0 du is called Abel’s integral equation
(t − u)n
where F (u) is unknown and G (t) is known and n is a constant between 0 and 1,
i. e., 0 < n < 1.
4. An integral equation in which various derivatives of the unknown function F (t)
can also be present is called an Integro transform equation.
A nswers
Multiple Choice Questions
1. (a) 2. (b) 3. (a)
True or False
1. F 2. T 3. T 4. F
¨
T-105
4
F ourier T ransforms
1 l nπx
where an = ∫ −l f ( x) cos dx …(2)
l l
1 l nπx
and bn = ∫ − l f ( x) sin l dx. …(3)
l
The series (1) with coefficients an and bn given by (2) and (3) respectively is called the
Fourier series of f (x), and the coefficients an and bn are called the Fourier coefficients
corresponding to f ( x).
At a point of discontinuity
1
f ( x) = [ f ( x + 0 ) + f ( x − 0 )].
2
If the function f ( x) defined in the interval (− l, l ) be an even function of x i.e. if
f (− x) = f ( x), then
1 l nπx 2 l nπx
an = ∫ f ( x) cos dx = f ( x) cos dx
l −l l l ∫0 l
1 l nπx
and bn = ∫ f ( x) sin dx = 0 .
l −l l
Therefore in this case we get Fourier cosine series.
Again if f ( x) be an odd function of x i.e. if f (− x) = − f ( x), then
1 l nπx
an = ∫ f ( x) cos dx = 0
l − l l
1 l nπx 2 l nπx
and bn = ∫ f ( x) sin dx = ∫ f ( x) sin dx
l − l l l 0 l
and thus in this case we get Fourier sine series.
Note: If f ( x) is a function of period 2l but is defined only in (0 , l ), we can extend it to
(− l, 0 ) so as to be an even or an odd function of x in the interval (− l , l )
1 ∞ ∞
0= ∫ −∞ f (v) sin w ( x − v) dw dv …(2)
2π ∫ − ∞
putting w = − p so that dw = − dp
1 ∞ ~
− ipx
= ∫ −∞ e f ( p) dp.
√ (2 π)
Note: Some authors also define the Fourier transform in the following forms :
~ ∞
(1) f ( p) = ∫ e − ipx f ( x) dx
−∞
T-108
1 ∞ ~
and f ( x) = ∫ −∞ f ( p) ⋅ e ipx dp.
2π
~ ∞
(2) f ( p) = ∫ e ipx f ( x) dx
−∞
1 ∞ ~
and f ( x) = ∫ −∞ e − ipx f ( p) dp.
2π
~ 1 ∞
(3) f ( p) = ∫ −∞ e − ipx f ( x) dx
√ (2π)
1 ∞ ~
and f ( x) = ∫ −∞ f ( p) ⋅ e ipx dp.
√ (2π)
2 ∞ ~
f ( x) = ∫0 f s ( p) sin px dp
π
1 ∞ ∞
= ∫0 dp ∫ −∞ { f (v) cos px cos pv + f (v) sin px sin pv} dv, where w = p
π
1 ∞ ∞ 1 ∞ ∞
= ∫0 cos px dp ∫ −∞ f (v) cos pv dv + ∫0 sin px dp ∫ −∞ f (v) sin pv dv
π π
1 ∞ ∞
or f ( x) = ∫0 cos px dp ∫ −∞ f ( x) cos px dx
π
1 ∞ ∞
+ ∫0 sin px dp ∫ −∞ f ( x) sin px dx …(1)
π
Now define f ( x) in (− ∞, 0 ) such that f ( x) is an odd function of x in (− ∞, ∞). Then
obviously f ( x) cos px is an odd function of x and f ( x) sin px is an even function of x in
(− ∞, ∞).
∴ From (1), we have
2 ∞ 2 ∞
f ( x) = ∫0 sin px dp × ∫0 f ( x) sin px dx
π π
2 ∞ ~
or f ( x) = ∫0 f s ( p) sin px dp .
π
Note: According to the authors who define
~ ∞
f s ( p) = ∫0 f ( x) sin px dx,
2 ∞ ~
we have f ( x) = ∫0 f s ( p) sin px dp.
π
2 ∞ ~
f ( x) = ∫0 f c ( p) cos px dp
π
at every point of continuity of f ( x).
This is an inversion formula for infinite Fourier cosine transform.
Proof: Proceeding as in article 4.7, we have
1 ∞ ∞
f ( x) = cos px dp ∫ f ( x) cos px dx
π ∫0 −∞
1 ∞ ∞
+ ∫ sin px dp ∫ f ( x) sin px dx. …(1)
π 0 −∞
2 ∞ ~
we have f ( x) = ∫0 f c ( p) cos px dx.
π
2 ∞
Now Fc { f (ax)} = ∫0 f (ax) cos px dx
π
1 2 ∞ p
= ⋅ ∫0 f (t) cos t dt ,
a π a
putting ax = t so that dx = (1 / a) dt
1 ~ p
= f c , from (1).
a a
putting x − a = t so that dx = dt
ipa 1 ∞ ipt
=e e f (t) dt
√ (2 π) ∫ − ∞
⋅
~
= e ipa f ( p) , from (1).
1 1 ∞ 1 ∞
= ∫ −∞ e i( p + a ) x f ( x) dx + ∫ e i ( p − a ) x f ( x) dx
2 √ (2 π) √ (2 π) − ∞
1 ~ ~
= [ f ( p + a) + f ( p − a)].
2
T-113
1 2 ∞
f ( x) ⋅ sin ( p + a) x dx
2 π ∫ 0
=
2 ∞
+ ∫0 f ( x) ⋅ sin ( p − a) x dx
π
1 ~ ~
= [ f s ( p + a) + f s ( p − a)].
2
(ii) We have
2 ∞
Fc { f ( x) sin ax} = ∫0 f ( x) sin ax cos px dx
π
2 1 ∞
= ⋅ ∫ f ( x) [sin ( p + a) x − sin ( p − a) x] dx
π 2 0
1 2 ∞
= ∫0 f ( x) sin ( p + a) x dx
2 π
2 ∞
− ∫0 f ( x) ⋅ sin ( p − a) x dx
π
1 ~ ~
= [ f s ( p + a) − f s ( p − a)].
2
(iii) We have
2 ∞
Fs f { ( x) sin ax} = ∫ f ( x) sin ax sin px dx
π 0
2 1 ∞
= ⋅ ∫ f ( x) [cos ( p − a) x − cos ( p + a) x] dx
π 2 0
T-114
1 2 ∞
f ( x) cos ( p − a) x dx
2 π ∫0
=
2 ∞ f ( x) cos ( p + a) x dx
− ∫
π 0
1 ~ ~
= [ f c ( p − a) − f c ( p + a)].
2
4.15 Theorem
If φ ( p) is the Fourier sine transform of f ( x) for p > 0, then
Fs { f ( x)} = − φ (− p) for p < 0 .
Proof: We have
2 ∞
Fs { f ( x)} = ∫ F ( x) sin px dx
π 0
2 ∞
= − ∫0 f ( x) sin sx dx = − φ (s)
π
= − φ (− p ), for p < 0.
Hence in general
φ (| p |), p > 0
Fs { f ( x)} =
− φ (| p |), p < 0
or Fs { f ( x)} = φ (| p |). Sgn p,
+ 1, p > 0
where the symbol Sgn p =
− 1, p < 0 .
~ 1 ∞ ∞
or F ( p, q) = f ( x, y) e i ( px + qy )dx dy
2 π ∫ −∞ ∫ −∞
which is Fourier transform of f ( x, y).
Inversion formula: Using inversion formula for Fourier transforms, we have
1 ∞ ~
f ( x, y) = ∫ f ( p, y) e − ipx dp
√ (2π) − ∞
~ 1 ∞ ~
and f ( p, y) = ∫ −∞ F ( p, q) e − iqy dq.
√ (2π)
1 ∞ ~
Hence f ( x, y) = ∫ −∞ F ( p, q) e − i( px + qy )dp dq,
2π
which is the inversion formula for the Fourier transform of f ( x, y).
4.17 Convolution
The function
1 ∞
H ( x) = F * G = ∫ −∞ F (u). G ( x − u) du
√ (2π)
is called the convolution or Falting of two integrable functions F and G over the
interval (− ∞, ∞).
Note: Some authors also define
∞
F*G= ∫ −∞ F (u). G ( x − u) du.
1 ∞ ∞
f (u) ∫ g ( y) e ip ( u + y )dy du,
2π ∫ − ∞
=
− ∞
1 ∞ 1 ∞
= ∫ f (u) e ipu ∫ g ( x) e ipx dx du
√ (2 π) − ∞ √ (2 π) − ∞
1 ∞ ipu
f (u) [e F { g ( x)}] du
√ (2 π) ∫ − ∞
=
1 ∞
= ∫ f (u) e ipu du F { g ( x)}
− ∞
√ (2 π)
1 ∞
= ∫ f ( x) e ipx dx F { g( x)}
− ∞
√ (2 π)
= F { f ( x)} ⋅ F { g ( x)} .
0 ∞
=∫ 0 ⋅ e ipt dt + ∫0 e − xt g(t) ⋅ e ipt dt
−∞
∞ ∞
= ∫0 e( ip− x ) t g (t) dt = ∫0 e − st g (t) dt , putting x − ip = s
= L { g (t)}.
Hence the Fourier transform of the function f (t) defined by (1) is the Laplace transform of the
function g (t).
Proof: By definition
1 ∞
F { f ′ ( x)} = ∫ −∞ f ′ ( x) ⋅ e ipx dx
√ (2π)
1 ∞ f ( x + h) − f ( x) ipx
= ∫ lim ⋅ e dx
√ (2 π) − ∞ h→ 0 h
1 ∞ f ( x + h) ipx 1 ∞ f ( x) ipx
= lim ∫ ⋅ e dx − lim ∫ −∞ ⋅ e dx
h→ 0 √ (2 π) − ∞ h h→ 0 √ (2 π) h
~
1 ∞ f ( x + h) ip( x + h ) − iph f ( p)
= lim ∫ ⋅e e d ( x + h) − lim
h→ 0 √ (2π) − ∞ h h→ 0 h
~
e − iph ∞ 1 f ( p)
= lim ∫ ⋅ f ( y) e ipy dy − lim
h→ 0 √ (2 π) − ∞ h h→ 0 h
~ ~
e − iph f ( p) f ( p)
= lim − lim
h→ 0 h h→ 0 h
~ e − iph − 1 ~
= f ( p) ⋅ lim = (− ip) f ( p).
h→ 0 h
(b) The Fourier transform of f n( x), the nth derivative of f ( x) is (− ip)n times the Fourier
transform of f ( x) provided that the first (n − 1) derivatives of f ( x) vanish as x → ± ∞ .
Proof: By definition
1 ∞
F { f n ( x)} = ∫ −∞ f n ( x) . e ipx dx.
√ (2π)
~ n ~
2
f c2 n+1 ( p) = − ∑ (− 1)r α2 n−2 r p2 r + (−1)n p2 n+1 f s ( p);
π
r =1
~ n ~
2
f s2 n ( p) = − ∑ (− 1)r α2 n−2 r p2 r −1 + (−1)n+1 p2 n f s ( p);
π
r =1
~ n ~
2
and f c2 n+1 ( p) = − ∑ (− 1)r α2 n−2 r −1 p2 r −1 + (−1)n+1 p2 n+1 f c ( p);
π
r =1
Thus, we have
~ n ~
2
f s2 n( p) = − ∑ (− 1)r α2 n−2 r p2 r −1 + (−1)n+1 p2 n f s ( p)
π
r =1
~ 2 n
and f c2 n + 1( p) = − ∑ (− 1)r α2 n−2 r +1 p2 r −1
π
r =1
~
+ (− 1)n+1 p2 n+1 f c ( p).
Note. The infinite sine and cosine transforms can be applied when the range of the
variable selected for exclusion is 0 to ∞.
Solution: We have
1 ∞
F { f ( x)} = ∫ −∞ e ipx f ( x) dx
√ (2π)
1 a b ipx ∞
0 ⋅ e ipx dx + e ⋅ e iωx dx + ∫ 0 ⋅ e ipx dx
√ (2 π) ∫ − ∞ ∫a
=
b
b
1 b i ( p + ω)x 1 e i( p + ω ) x
e dx
√ (2 π) ∫ a
= = ⋅
√ (2 π) i ( p + ω)
a
1 ei ( p + ω ) a − ei ( p + ω ) b
= ⋅
√ (2 π) p+ω
T-120
1 e ipa e − ipa
= −
√ (2 π) ip ip
2 i sin pa 2 sin pa
= = , p ≠ 0.
ip √ (2 π) p √ (2 π)
~
For p = 0 , F ( p) = 2 a / √ (2 π).
Solution: Let f ( x) = e − x .
~
Then f s ( p) = 2 ∞ 2
f ( x) sin px dx =
∞ −x
e sin px dx
π ∫0 π ∫0
∞
2 e− x
= (− sin px − p cos px)
π 1 + p2 0
p 2
=
2 π
1+ p
~
and f c ( p) = 2 ∞ 2
f ( x) cos px dx =
∞
e − x cos px dx
π ∫0 π ∫0
∞
2 e− x
= (− cos px + p sin px)
π 1 + p2 0
1 2 ⋅
=
2 π
1+ p
Applying inversion to the sine transform, we have
2 ∞ ~
f ( x) = ∫0 f s ( p) ⋅ sin px dp
π
2 ∞ p sin px
= ∫0 dp …(1)
π 1 + p2
and applying inversion to the cosine transform, we have
2 ∞ ~
f ( x) = ∫0 f c ( p) cos px dp
π
2 ∞ cos px
= ∫0 dp. …(2)
π 1 + p2
Now from Fourier integral theorem, we have
1 ∞ ∞
f ( x) = ∫ dp ∫ f (v) cos p ( x − v) dv
π 0 − ∞
1 ∞ ∞
or f ( x) = cos px dp ∫ f (v) cos pv dv
π ∫0 −∞
1 ∞ ∞
+ ∫ sin px dp ∫ f (v) sin pv dv. …(3)
π 0 − ∞
T-122
∞ cos px π −x
∴ ∫0 dp = e .
2 2
1+ p
2 π −x
∴ from (2) we have f ( x) = ⋅ e = e− x .
π 2
Case II: Again defining f ( x) in (− ∞, 0 ) such that f ( x) is an odd function of x, from
(2), we have
2 ∞ ∞
f ( x) = ∫0 sin px dp ∫0 f (v) sin pv dv.
π
Taking f ( x) = e − x and simplifying, we have
∞ p sin px π −x
∫0 dp = e .
1 + p2 2
2 π −x
∴ from (1), f ( x) = ⋅ e = e− x .
π 2
1
Example 4: Find Fourier cosine transform of f ( x) = and hence find Fourier sine
1 + x2
x
transform of F ( x) = ⋅
1 + x2 (Kanpur 2007, 10; Meerut 13; Rohilkhand 14)
Solution: We have
~
f c ( p) = 2 ∞
f ( x) cos px dx
π ∫0
2 ∞ cos px
= ∫0 dx .
π 2
1+ x
Differentiating both sides w.r.t. p, we have
d ~ 2 ∞ x sin px
f c ( p) = − ∫0 dx
dp π 2
1+ x
T-123
2 ∞ ( x2 + 1 − 1) sin px
= − ∫0 dx
π 2
x (1 + x )
2 ∞ sin px 2 ∞ sin px
= − ∫0 dx + ∫0 dx
π x π x (1 + x2 )
2 π 2 ∞ sin px ∞ sin px π
= − ⋅ + ∫0 dx . ∵ ∫0 dx =
π 2 π 2
x (1 + x ) x 2
d ~ π
and f c ( p) = − .
dp 2
π
∴ from (1), we have = A + B
2
π
and − = A − B .
2
π
Solving, A = 0 , B = .
2
~ π
∴ from (1), we have f c ( p) = e − p.
2
1 dn −1 dn
= n (a − ip) − n
(a + ip)−1
2 i da da
1
= (− 1)n (n !) [(a − ip)−( n+1) − (a + ip)−( n+1)]
2i
1
= (− 1)n (n !) [2 i r −( n+1) sin (n + 1) θ], putting a = r cos θ, p = r sin θ
2i
= (− 1)n n !(1 / r)n+1 sin (n + 1) θ.
∞
∴ ∫0 x n e − ax sin px dx = (n !) . [1 / (a2 + p2 )( n+1)/2 ] sin {(n + 1) tan−1 ( p / a)}.
1
= (− 1)n (n !) [(a − ip)−( n+1) + (a + ip)−( n+1)]
2
= (− 1)n (n !) (1 / r)n+1 cos (n + 1) θ,
∞ ∞
We have ∫0 e − ax x n cos px dx + i ∫0 e − ax x n sin px dx
∞
=∫ e − ax x n (cos px + i sin px) dx
0
∞ ∞
= ∫0 e − ax x n e ipx dx = ∫0 e −( a − ip ) x x ( n + 1) −1 dx
Γ (n + 1) ∞ Γ (n)
= ∵ ∫0 e − az z n − 1 dz =
(a − ip ) n+1 a n
Γ (n + 1)
= , putting a = r cos θ and p = r sin θ
n+1
r (cos θ − i sin θ) n + 1
n! 1
= (cos θ + i sin θ) n + 1 ∵ = cos θ + i sin θ
n +1
r cos θ − i sin θ
n!
= [cos (n + 1) θ + i sin (n + 1) θ]. ...(1)
r n +1
Equating real and imaginary parts on both sides of (1), we have
∞ n!
∫0 e − ax x n cos px dx = cos (n + 1) θ
r n +1
∞ n!
and ∫0 e − ax x n sin px dx = sin (n + 1) θ ,
r n +1
where r2 = a2 + p2 i. e., r = (a2 + p2 )1 /2 and θ = tan−1 ( p / a).
T-126
∞
Hence ∫0 e − ax x n cos px dx
n!
= cos { (n + 1) tan−1 ( p / a)}
2 2 ( n+1)/2
(a + p )
∞
and ∫0 e − ax x n sin px dx
n!
= sin {(n + 1) tan−1 ( p / a)}.
2 2 ( n+1)/2
(a + p )
e ax + e − ax
Solution: If f ( x) = πx , then we have
e − e − πx
~
f s ( p) = 2 ∞
f ( x) sin px dx
π ∫0
2 ∞ e ax + e − ax
= ∫ sin px dx
π 0 πxe − πx
−e
2 ∞ e ax + e − ax e ipx − e − ipx
= ∫0 ⋅ dx
π e πx − e − πx 2i
2 1 ∞ e( a + ip) x − e −( a + ip) x
= ∫ dx
π 2 i 0
e π x − e − πx
1 ∞ e( a − ip) x − e −( a − ip) x
dx
2i ∫ 0
−
e π x − e− π x
2 1 1 a + ip 1 1 a − ip
= ⋅ tan − ⋅ tan
π 2 i 2 2 2i 2 2
∞ e az − e − az 1 a
∵ From definite integrals, ∫ 0 e π z − e− π z dz = 2 tan 2
a + ip a − ip
sin sin
2 1 2 1 2
= −
π 4i a + ip 4 i a − ip
cos cos
2 2
a + ip a − ip a − ip a + ip
sin cos − sin cos
2
= 2 2 2 2
π a + ip a − ip
4 i cos cos
2 2
T-127
Solution: We have
~
f s ( p) = 2 ∞ 1
sin px dx . …(1)
π ∫0 x (a + x2 )
2
∞ 1
Let I = ∫0 sin px dx. …(2)
x (a + x2 )
2
dI d ∞ sin px
Then = ∫0 dx
dp dp x (a2 + x2 )
∞∂ sin px
=∫
0 ∂p x (a2 + x2 )
dx
∞ cos px
= ∫0 dx. …(3)
a2 + x2
d2 I ∞ x sin px
∴ =− ∫0 dx
dp2 a2 + x2
∞ x2 sin px ∞ ( x2 + a2 ) − a2
=− ∫0 dx = − ∫0 sin px dx
x (a2 + x2 ) x (a2 + x2 )
∞ sin px ∞ sin px
=− ∫0 dx + a2 ∫0 dx
x x (a2 + x2 )
π ∞ sin px π
=− + a2 I. ∵ ∫0 dx = ⋅
2 x 2
2
d I π
∴ − a2 I = −
2 2
dp
π d
or ( D2 − a2 ) I = − , where D ≡ ⋅
2 dp
dI
∴ = − Aae − ap + Bae ap. …(5)
dp
Now from (2), when p = 0, we have I = 0 and from (3), when p = 0, we have
∞
dI ∞ 1 1 tan−1 x = π ⋅
= ∫0 dx =
dp a2 + x2 a a 0 2 a
Solution: We have
~ 2 2 ∞ 2
Fc { e − x } = ∫0 e − x cos px dx = I …(1)
π
1 2 ∞ 2
= ∫0 (− 2 xe − x ) ⋅ sin px dx
2 π
1 2 − x2 ∞ − x2
= (e sin px)0∞ − p ∫ e cos px dx
2
π 0
2 ∞ − x2 1
But when p = 0, from (1) I = ∫0 e dx = ⋅
π √2
∴ from (2), A = 1 / √ 2.
2 2
Hence I = Fc { e − x } = (1 / √ 2) e − p /4 .
~ p
Example 9: Use the sine inversion formula to obtain f ( x) if f s ( p) = ⋅
1 + p2
Solution: Using Fourier sine inversion formula, we have
2 ∞ p
f ( x) = ∫0 ⋅ sin px dp
π 1 + p2
2 ∞ p2 2 ∞ ( p2 + 1) − 1
= ∫ 0 p (1 + p2 ) ⋅ sin px dp = π ∫ ⋅ sin px dp
π 0 2
p (1 + p )
2 ∞ sin px 2 ∞ sin px
= ∫0 dp − ∫ dp
π p π 0 p (1 + p2 )
π 2 ∞ sin px
or f ( x) = − ∫0 dp …(1)
2 π p (1 + p2 )
∞ sin px π
∵ ∫0 dp =
p 2
df 2 ∞ cos px
∴ = − ∫0 dp …(2)
dx π 2
1+ p
d2 f 2 ∞ p sin px
and = ∫0 dp
2 π 2
dx 1+ p
2
d f
or − f =0
dx2
whose solution is f = A e x + B e − x . …(3)
df
∴ = A e x − B e− x . …(4)
dx
π
Now when x = 0 , f = , from (1)
2
df 2 ∞ dp π
and = − ∫0 = − , from (2).
dx π 1 + p2 2
π π
∴ from (3) and (4), = A + B and − = A − B .
2 2
π π
Solving, A = 0 , B = . Hence f ( x) = e − x .
2 2
T-130
Comprehensive Exercise 1
∞ Γ (m) Γ (m) mπ mπ
= ∫0 e ipx x m − 1 dx = m
= cos
m
− i sin
(ip) p 2 2
T-131
e ax + e − ax
(ii) Find the cosine transform of ⋅
e πx + e − πx (Purvanchal 2014; Kanpur 14)
1
6. Find the sine transform of and deduce that
e πx − e − πx
1
Fs (cosech πx) = tanh ( p / 2).
√ (2 π)
7. Find the Fourier sine transform of f ( x) , if
0 , 0 < x < a
f ( x) = x, a ≤ x ≤ b
0 , x > b.
8. Find f ( x) if its cosine transform is 1/(1 + p2 ).
~ 1
9. Find f ( x) if f c ( p) = πe − p.
2
10. Find f ( x) if its sine transform is π / 2.
11. Find f ( x) if (i) its sine transform is e − ap, (Kanpur 2012)
− ap
(ii) its cosine transform is e .
~
12. Find the inverse Fourier transform of f ( p) = e −| p | y .
~
13. Find f ( x) if f s ( p) = pn e − ap.
14. Find f ( x) if its cosine transform is
~ 1 p
a − , if p < 2 a
f c ( p) = √ (2 π) 2
0, if p ≥ 2 a
~ e − ap
15. Find f ( x) if f s ( p) = . Hence deduce Fs −1 {1/ p}.
p
A nswers 1
sin pε i 2
1. (ii) (ii) − ⋅ (ap cos ap − sin ap)
pε 2 π
p
1 2
2. (ii) ⋅ [1 − cos pa]
2 π
ap
2 sin p 2 cos p
4. (i) 2 ⋅ (1 − cos p); 2 ⋅ (1 − cos p).
π p2 π p2
Γ(m) 2 mπ Γ(m) 2 mπ
(ii) m
sin ; m cos
p π 2 p π 2
p /2
2 cos (a / 2) ⋅ (e + e − p /2 )
5. (i) √ (π / 2) . e − p (ii)
π 2 cos a + e + e − p
p
1 ep − 1
6.
2 √ (2 π) ep + 1
8. √ (π / 2) e − x
9. π ⋅ 1 10. √ (π / 2).(1 / x)
2 1 + x2
2 ⋅ x 2 ⋅ a y √2
11. 2 ; 2 12.
π a + x2 π a + x2 √ π( y2 + x2 )
−1
13. 2 ⋅ n !sin {(n + 1) tan ( x / a)} 14. π −1 x −2 sin2 ax
π (a2 + x2 )( n + 1)/2
15. √ (2 / π) tan−1 ( x / a) ; √ (π / 2)
True or False
Write ‘T’ for true and ‘F’ for false statement.
1. The Fourier transform is a linear transformation.
~
2. If f c ( p) is the infinite Fourier cosine transform of f ( x), then the infinite
~ p
Fourier cosine transform of f (ax) is f c ⋅
a
3. 2 ⋅ x
is the inverse Fourier sine transform of e − ap.
2
π a + x2
~ ~
4. If f s ( p) and f c ( p) are infinite Fourier sine and cosine transforms of f ( x)
1 ~ ~
respectively, then Fs { f ( x) sin ax} = [ f c ( p − a) + f c ( p + a).
2
T-134
~ ~
5. The Fourier transform of f ′ ( x), the derivative of f ( x) is (− ip) f ( p), where f ( p)
A nswers
True or False
1. T 2. F 3. T 4. F 5. T
6. T
¨
T-135
5
F inite F ourier T ransforms
1~ 2 ∞ ~ pπx
f ( x) =
l
f c (0 ) + ∑
l p=1
f c ( p) cos
l
,
~ l
where f c (0 ) = ∫0 f ( x) dx .
If π is taken as the upper limit for the finite Fourier cosine transform then the inversion formula is
given by
1 ~ 2 ∞ ~
f ( x) =
π
f c (0 ) + ∑ f c ( p) cos px ,
π p=1
~ π
where f c (0 ) = ∫0 f ( x) dx .
2 l pπx 2 ~
where ap = ∫0 f ( x) cos
dx = f c ( p).
l l l
2 l 2 ~
∴ a0 = ∫ 0 f ( x) dx = l f c (0).
l
1~ 2 ∞ ~ pπx
∴ f ( x) =
l
f c (0 ) + ∑ f c ( p) cos l .
l p=1
Note: The upper limit for Fourier sine or cosine transforms will be taken as x = π, if not
given in the problem.
~ π ~
Fs ( p, q) = ∫ f ( p, y) sin qy dy.
0 s
~ π π
∴ Fs ( p, q) = ∫0 ∫0 f ( x, y) sin px sin qy dx dy,
2 ∞ ~
and f ( x, y) = ∑ f s ( p, y) sin px.
π p=1
∞ ∞ ~
4
Hence f ( x, y) =
π 2 ∑ ∑ Fs ( p, q) sin px sin qy.
p=1q =1
Similarly we can find the inversion formula for double finite cosine transform of
f ( x, y).
π π
= [ f ′ ( x) sin px ]0 − p ∫ f ′ ( x) cos px dx
0
~
= − p f (π) cos pπ + p f (0 ) − p2 f s ( p).
~
Hence Fs { f ′ ′ ( x)} = − p2 f s ( p) + p { f (0 ) − (− 1) p f (π)}.
T-139
π π
= [ f ′ ( x) cos px ]0 + p ∫ f ′ ( x) sin px dx
0
(Integrating by parts taking f ′ ′ ( x) as second function)
π π
= f ′ (π) cos pπ − f ′ (0 ) + p [ f ( x) sin px]0 − p2 ∫ f ( x) cos px dx
0
~
= (− 1) p f ′ (π) − f ′ (0 ) − p2 f c ( p).
~
Hence Fc { f ′ ′ ( x)} = − p 2 f c ( p) − f ′ (0 ) + (− 1) p f ′ (π). ...(1)
Note: If f ′ ′ ( x) and f ′ ′ ′ ( x) are continuous and if f (4 ) ( x) is sectionally continuous,
then replacing f ( x) by f ′ ′ ( x) in (1), we have
Fc { f 4 ( x)} = − p2 Fc { f ′ ′ ( x)} − f ′ ′ ′ (0 ) + (−1) p f ′ ′ ′ (π)
~
= − p2 [− p2 f c ( p) − f ′ (0 ) + (− 1) p f ′ (π)] − f ′ ′ ′ (0 ) + (− 1) p f ′ ′ ′ (π)
T-140
~
= p 4 f c ( p) + p2 [ f ′ (0 ) − (− 1) p f ′ (π)] − f ′ ′ ′ (0 ) + (−1) p f ′ ′ ′ (π).
Similarly finite Fourier cosine transforms of other derivatives of even order can be
found.
~
Theorem 2: If f c ( p) is the cosine transform of a sectionally continuous function f ( x), 0 ≤ x ≤ π,
then
~ ~
−1 f c ( p) x π f (0 )
Fc = ∫0 ∫ t f (r) dr dt + c ( x − π)2 + A,
2 2π
p
π π
= [ f ( x) sin px]0 − p ∫ f ( x) cos px dx
0
(Integrating by parts, taking f ′ ( x) as second function)
= − p Fc { f ( x)}, ( p = 1, 2, 3,...).
(b) We have
π
Fc { f ′ ( x)} = ∫0 f ′ ( x) cos px dx
π π
= [ f ( x) cos px]0 + p ∫ f ( x) sin px dx
0
[Integrating by parts]
= f (π) cos pπ − f (0 ) + p Fs { f ( x)}, ( p = 0 , 1, 2,...)
= pFs { f ( x)} − f (0 ) + (−1) p f (π).
Note: If H ( x) is any sectionally continuous function, the above properties may be
written as
x
Fs { H ( x)} = − p Fc ∫ H (r) dr , ( p = 1, 2, ......)
0
1 ~ π x ~
and Fc H ( x) − Hc (0 ) = p Fs ∫ H (r) dr − Hc (0 ) ,
π 0 π
( p = 0 , 1, 2, ......)
T-141
5.9 Convolution
Let F ( x) and G ( x) be two functions defined on the interval −2 π < x < 2 π, then the function
π
F ( x) * G ( x) = ∫−π F ( x − y) G ( y) dy
Example 1: Find the finite Fourier sine and cosine transforms of the function
f ( x) = 2 x, 0 < x < 4. (Agra 2003)
Solution: We have
~ l
f s ( p) = ∫0 f ( x) sin ( pπx / l ) dx]
4
= ∫ 2 x sin ( pπx / 4) dx, as l = 4
0 (Given)
4
− 2 x cos ( pπx / 4) 4 cos ( pπx / 4)
= +2∫ dx
pπ / 4 0 pπ / 4
0
4
32 8 sin ( pπx / 4) 32
=− cos pπ + =− cos pπ.
pπ pπ pπ / 4 0 pπ
~ l 4
Also f c ( p) = ∫0 f ( x) cos ( pπx / l) dx = ∫ 2 x cos ( pπx / 4) dx, as l = 4
0
4
2 x sin ( pπx / 4) 4 sin ( pπx / 4)
=
p π / 4 − 2 ∫0 pπ / 4
dx
0
4
8 cos ( pπx / 4) 32
= = (cos pπ − 1) , if p > 0
pπ pπ / 4 0 p2 π2
~ 4
and if p = 0, then f c ( p) = ∫0 2 x ⋅ 1 dx = 16.
π
1 x 1 1 π1
= − − −1 + cos px − ∫ cos px dx
p
π p 0 p2 0 π
1 1 π 1
3 [
= − sin px]0 = , if p > 0
2
p p π p2
and when p = 0,
~ π π x2
f c ( p) = ∫0 3
− x + dx = 0
2π
~ π
(ii) f c ( p) = ∫0 sin nx cos px dx
1 π
= ∫0 [sin (n + p) x + sin (n − p) x] dx
2
π
1 cos (n + p) x cos (n − p) x
=
2
−
n+ p
−
n− p , if p ≠ n
0
~ 1 cos (n + p) π cos (n − p) π 1 1
∴ if p ≠ n, f c ( p) = − − + + ⋅
2 n + p n − p n + p n − p
If p = n,
~ π 1 π
then f c ( p) = ∫sin nx cos nx dx = sin 2 nx dx
∫
0 2 0
π
1 cos 2 nx
= − =0.
2 2 n 0
~ 2n
∴ f c ( p) = 0 or according as n − p is even or odd.
n − p2
2
cos k (π − x)
Example 3: Find the finite cosine transform of f ( x) if f ( x) = − ⋅
k sin kπ
Solution: We have
~ π cos { k (π − x)}
f c ( p) = − ∫ cos px dx
0 k sin k π
1 π
=− ∫0 [cos { k (π − x) + px} + cos { k (π − x) − px}] dx
2 k sin kπ
T-143
π
1 sin (kπ − kx + px) sin (kπ − kx − px)
=− −
2 k sin kπ p− k p+ k
0
1 sin pπ sin (− pπ) sin kπ sin kπ
=− − − +
2 k sin kπ p − k p+ k p− k p + k
1 1 1 1
= − = 2 , k ≠ 0 , 1, 2, 3, ......
2k p − k p + k p − k2
π
π sin kx x cos k (π − x)
=∫
0
2
− sin px dx
2 k sin kπ 2 k sin kπ
π π
= ∫0 sin kx sin px dx
2 k sin2 kπ
1 π
− ∫0 x cos k (π − x) sin px dx
2 k sin kπ
π π
= ∫0 [cos ( p − k ) x − cos ( p + k ) x] dx
2
4 k sin kπ
1 π
− ∫0 x [sin (kπ − kx + px) + sin ( px − kπ + kx)] dx
4 k sin kπ
π
π sin ( p − k ) x sin ( p + k ) x
= −
2 p− k p+ k
4 k sin kπ 0
π
1 cos (kπ − kx + px) cos ( px − kπ + kx)
− ⋅ x− −
4k sin kπ ( p − k) ( p + k) 0
1 π cos (kπ − kx + px) cos ( px − kπ + kx)
+
4 k sin kπ ∫ 0 1⋅ − ( p − k)
−
( p + k)
dx
π sin ( p − k ) π sin ( p + k ) π
= −
2 p− k p+ k
4 k sin kπ
π 1 1
+ cos pπ ⋅ +
4 k sin kπ p − k p + k
π
1 sin (kπ − kx + px) sin ( px − kπ + kx)
− +
4 k sin kπ ( p − k )2 ( p + k )2 0
T-144
π 1 1
= ⋅ sin pπ cos kπ ⋅ −
2 p − k p + k
4 k sin kπ
1 1 π 2p
− cos pπ sin kπ ⋅ + + ⋅ cos pπ
p − k p + k 4 k sin k π p − k2
2
1 1 1
+ ⋅ sin kπ ⋅ −
2
4 k sin kπ
( p − k ) ( p + k )2
2 2
1 ( p + k) − ( p − k)
= ⋅
4k ( p2 − k 2 )2
p
= , (| k | ≠ 0 , 1, 2,...).
( p − k 2 )2
2
Example 5: Find f ( x) if
pπ
~ 6 sin − cos pπ
2 2
f c ( p) = for p = 1, 2, 3,...... and for p = 0,
(2 p + 1) π π
where 0 < x < 4.
1~ 2 ∞ ~ pπx
Solution: We have f ( x) =
l
f c (0 ) + ∑ f c ( p) cos l
l p=1
pπ
6 sin ∞ − cos pπ
1 2 2 2 pπx
= ⋅ + ⋅ ∑ cos
4 π 4 p=1 (2 p + 1) π 4
Solution: We have
2 ∞ ~
f ( x) = ∑ f s ( p) sin px
π p=1
2 ∞ 1 − cos pπ
= ∑ sin px
π p = 1 p2 π2
∞
2 1 − cos pπ
= ∑
π3 p = 1 p2
sin px.
T-145
Comprehensive Exercise 1
A nswers 1
1
1. (i) [1 − (− 1) p]; 0
p
π (−1) p + 1 (−1) p − 1 π2
(ii) ; , if p = 1, 2, 3,... and , if p = 0 ⋅
p 2 2
p
1 (−1) p + 1
2. (i) ;
p 4p
1 1
(ii) [1 − (−1) p] ; [(−1) p − 1]
2
πp 4 πp2
3. (i) (2 / p2 ) sin ( pπ / 2) , (ii) (π / p) cos pc .
4. (i) (2 / p) sin ( pπ / 2), p > 0 and 0, if p = 0.
2 π
(ii) , if p > 0 and , if p = 0
2 3
πp
2 6π
5. (i) [1 − (− 1) p], (ii) (−1) p + 1
3
p p3
6. (i) 0, if p ≠ n , and π / 2, if p = n
π 1
(ii) {(− 1) p + 2} + {(− 1) p − 1}
6p πp3
p 6 π2
7. (i) [1 − (−1) p cos kπ], (ii) π (− 1) p − ,
2
p −k 2
p
3 p
p
(iii) [1 − (−1) p e cπ ]
c + p2
2
c p
8. 9. , k ≠ 0 , 1, 2,... .
2 2
c +p p − k2 2
64 128 128 64
10. − cos pπ + (cos pπ − 1) ; cos p π, if p > 0; , if p = 0
pπ 3 3 2 2 3
p π p π
∞
cos (2 pπ / 3)
11. 1 + 2 ∑ (2 p + 1)2
cos pπx
p=1
True or False
Write ‘T’ for true and ‘F’ for false statement.
~
1. If f s ( p) is the finite Fourier sine transform of f ( x) on the interval (0 , π) then
the
2 ∞ ~
inversion formula for sine transform is f ( x) = ∑ f s ( p)sin px
π p=1
~ 1 − cos pπ
4. If the finite sine transform is given by f s ( p) = , where 0 < x < π , then
p2 π2
∞ 1 − cos pπ
2
the value of f ( x) is
π 3 ∑
p2
sin px.
p =1
A nswers
Multiple Choice Questions
1. (c) 2. (a) 3. (c) 4. (d)
Fill in the Blank(s)
π
1. ∫0 f ( x ) sin px dx 2. inverse finite 3. Convolution
− 32
4. cos pπ
pπ
True or False
1. T 2. F 3. T 4. T
¨
T-149
6
A pplications
of F ourier
T ransforms in I nitial and
B oundary V alue P roblems
∞ ∂V ∂V
=−p
∫0 cos px ⋅ ∂x dx, if
∂x
→ 0 as x → ∞
T-150
∞
∞
= − p (V ⋅ cos px )0 + p
∫0 V sin px dx
~ ~
= p (V ) x = 0 − p2 V s where V s is the sine transform of V and
assuming that V → 0 as x → ∞ .
Again applying the cosine transform, we have
∞
∞ ∂2V ∂V ∞ ∂V
∫0 cos px ⋅ 2
∂x
dx =
∂x
cos px + p
0 ∫0 sin px ∂x
dx
∂V ∞
∞
+ p [V sin px]0 − p2
=−
∂x x = 0 ∫0 V cos px ⋅ dx
∂V
if → 0 as x → ∞
∂x
∂V ~
=− − p2 V c if V → 0 as x → ∞,
∂x x = 0
~
where V c is the cosine transform of V.
∂2V
Thus we see that for the exclusion of from a differential equation we require
∂x2
(V ) x = 0 in sine transform
∂V
and in cosine transform.
∂x x = 0
Note 1: By the Fourier sine or cosine transform we cannot exclude a derivative of odd
order from the given differential equation.
Note 2: When one of the variables in a differential equation ranges from − ∞ to ∞ then
that variable can be excluded with the help of complex Fourier transform.
2µ ∞ cos px 2
1 − e − kp t dp .
U ( x, t) =
π ∫0 p2
∂U
Solution: Since is given, so taking the Fourier cosine transform of both the
∂x x = 0
sides, we have
∞ ∞ ∂2U
2 ∂U 2
π ∫0 ∂t
cos px dx = k
π ∫0 ∂x2
cos px dx
∞
∞
2 d 2 ∂U
or
π dt ∫0 U cos px dx = k
π ∂x
cos px
0
∞
2 ∂U
+ kp
π ∫0 ∂x
sin px dx
d ~ 2 ∂U 2 ∞
or (U c ) = − k ⋅ + kp (U sin px )0
dt π ∂x x = 0
π
∞
2
− kp2
π ∫0 U cos px dx
∂U
if → 0 as x → ∞
∂x
~
dUc 2 ~
2
or = ⋅ kµ − kp U c if U → 0 as x → ∞
dt π
d ~ ~ 2
or U c + kp2 U c = kµ ,
dt π
which is a linear differential equation of first order.
kp2 dt 2
I.F. = e ∫ = e kp t .
∴ its solution is
2 ~ 2 2
e kp t
⋅ U c = A + k µ e kp t dt
π ∫
2 µ 2
= A + 2 e kp t .
π p
~
But when t = 0 ,Uc = 0 . [∵ U = 0 when t = 0 ]
2 µ
∴ 0 = A+ 2
π p
µ 2
or A= − 2
⋅√ .
p π
2 ~ 2 µ 2
∴ e kp t U c = ⋅ 2 − 1 + e kp t
π p
~ 2 µ 2
or U c = ⋅ 2 (1 − e − kp t ).
π p
∞ ~
2
Uc =
π ∫0 U c cos px dp
2µ ∞ cos px 2
⋅ 1 − e − kp t dp .
or Uc =
π ∫0 p 2
∂U ∂2U
Example 2: Solve =2 2
∂t ∂x
if U (0 , t) = 0 , U ( x, 0 ) = e − x , x > 0 , U ( x, t) is bounded where x > 0 , t > 0 . (Meerut 2013B)
Solution: Since U (0 , t) is given, so taking the Fourier sine transform of both the sides,
we have
∞ ∞ ∂2U
2 ∂U 2
π ∫0 ∂t
sin px dx = 2
π ∫0 ∂x2
sin px dx
∞ ∞
d 2 2 ∂U
or
dt
π ∫0 U sin px dx = 2 ⋅
π ∂x
sin px
0
∞
2 ∂U
−2 p
π ∫0 ∂x
cos px dx
∂U
if → 0 as x → ∞ .
∂x
d ~ 2 ∞ 2 ∞
U s = − 2 p (U cos px )0 − 2 p2
or
dt π π ∫0 U sin px dx
if U → 0 as x → ∞
2 ~
= 2 p U (0 , t) − 2 p2 U s
π
~
dUs ~ 2
or + 2 p2 U s = 2 p U (0 , t) = 0
dt π
~ 2
whose solution is U s = Ae −2 p t
…(1)
−x
But U ( x, 0 ) = e .
~ 2 ∞
e − x sin px dx
∴ when t = 0 , U s =
π ∫0
∞
−x
2 e
= 2
(− sin px − p cos px)
π 1 + p
0
p 2
= .
1 + p2 π
~ 2 p 2
Hence Us = e −2 p t .
π 1 + p2
Example 3: Using the Fourier sine transform, solve the partial differential equation
∂V ∂2V
= k 2 for x > 0 , t > 0 , under the boundary conditions V = V0 when x = 0 , t > 0 , and the
∂t ∂x
initial condition V = 0, when t = 0 , x > 0 .
Solution: Taking the Fourier sine transform of both the sides, we have
2 ∞ ∂V 2 ∞ ∂2V
π ∫0 ∂t
sin px dx = k
π ∫0 ∂x2
sin px dx
d 2 ∞ ~ 2
V sin px dx = − kp2 V s + kp
or
dt π ∫0 V (0 , t)
π
~ 2 2 V 2
or V s e kp t
= c + 0 e kp t …(1)
π p
~
But when t = 0 , V s = 0 . [∵ V = 0 when t = 0]
2
2V0 ∞ sin px ∞ e − kp t
or V =
π
∫0 p
dp −
∫0 p
⋅ sin px dp
2
2V0 π ∞ e − kp t
= −
π 2 ∫0 p
sin px dp
2
2 ∞ e − kp t
or V ( x, t) = V0 1 −
π ∫0 p
sin px dp ⋅
Note: In terms of Error function or Complementary error function the solution may be
written as
x x
V ( x, t) = V0 1 − erf = V0 erf c as shown below.
2 √ (kt) 2 √ (kt)
We know that
∞ 2 2 √ π −β2 / α 2
x
e−α
∫0 cos 2βx dx =
2α
e .
∞ 2 2 1 √π β/α 2
x
e−α e − u du,
or
∫0 2x
sin 2βx dx =
2α
α
∫0
β
putting = u so that dβ = α du
α
∞ 2 2 sin 2βx π 2 β/ α 2
x
e−α e − u du
or
∫0 ⋅
x
dx = ⋅
2 √π ∫0
π
= erf (β / α) .
2
2 ∞ 2 sin px
e − kp t
∴
π ∫0 ⋅
p
dp
2 ∞ 2 sin 2 ux
e −4 ktu ⋅
=
π ∫0 u
du , putting p = 2 u so that dp = 2 du
2 π x x
= ⋅ erf = erf .
π 2 2 √ (kt) 2 √ (kt)
d2 1 ∞
ye ipx dx = c 2 (− ip)2 ~
or 2
dt √ (2 π) ∫− ∞ y ( p, t)
d2 ~
y ( p, t)
or = − c 2 p2 ~
y ( p, t)
dt2
d2 ~
y ( p, t)
or + c 2 p2 ~
y ( p, t) = 0
dt2
whose solution is
~y ( p, t) = A cos cpt + B sin cpt. …(2)
Initially the string is at rest.
∂y
∴ = 0, at t = 0 .
∂t
1 ∞ ∂y d 1 ∞
⋅ e ipx dx = ye ipx dx
∴
√ (2 π) − ∞ ∂t ∫ dt √ (2 π) − ∞ ∫
d ~y ( p, t)
= = 0 , at t = 0 .
dt
∴ from (2), we have 0 = Bcp or B = 0 .
Also at t = 0, y = f ( x) .
1 ∞ ~
at t = 0 , ~ f (u) e ipu du = f ( p) .
∴ y ( p, 0 ) =
√ (2 π) ∫− ∞
~
∴ from (2), we have f ( p) = A .
~ ~
Hence y ( p, t) = f ( p) cos cpt.
Taking the inverse Fourier transform, we have
1 ∞ ~
f ( p) cos cpt e − ipx dp
y ( x, t) =
√ (2π) − ∞ ∫
1 ∞ ∞
f (u) e ipu du cos cpt e − ipx dp
=
2π ∫ − ∞ ∫ − ∞
T-156
1 ∞ ∞
f (u) e ipu du (e icpt + e − icpt ) e − ipx dp
=
4π ∫− ∞ ∫− ∞
1 1 ∞ ∞
f (u) e − iαudu (e − icαt + e icαt ) e iαx dα ,
=
2 2π ∫− ∞ ∫− ∞
putting p = − α so that dp = − dα
1 1 ∞ ∞
f (u) e − iαu e iα ( x + ct)dα du
=
2 2π ∫− ∞
∫− ∞
1 ∞ ∞
f (u) e − iαu e iα ( x − ct) dα du
+
2π ∫− ∞
∫− ∞
1
= [ f ( x + ct) + f ( x − ct)] .
2
(From Fourier integral formula, see article 4.3)
Example 5: If the flow of heat is linear so that the variation of θ (temperature) with z and y may
be neglected and if it is assumed that no heat is generated in the medium, then solve the differential
∂θ ∂2θ
equation = k 2 (one dimensional heat equation) where − ∞ < x < ∞ and θ = f ( x) where
∂t ∂x
t = 0, f ( x) being a given function of x.
Solution: Taking the Fourier transform of both the sides of the given equation, we have
1 ∞ ∂θ ipx 1 ∞ ∂2θ ipx
√ (2 π) ∫− ∞ ∂t
e dx = k
√ (2 π) ∫ − ∞ ∂x2 e dx
d 1 ∞ ~
θ ⋅ e ipx dx = k (− ip)2 θ
or
dt √ (2 π) ∫− ∞ [From article 4.20 (b)]
~ ~
where θ = θ ( p, t) is the Fourier transform of θ ( x, t)
d ~ ~
or θ = − kp2 θ
dt
whose solution is
~ 2
θ = A e − kp t
…(1)
Now at t = 0 , θ = f ( x) .
~ 1 ∞ ~
f ( x) e ipx dx = f ( p) .
∴ at t = 0, θ =
√ (2 π) ∫− ∞
~ ~
∴ from (1), at t = 0, θ = f ( p) = A ;
~ ~ 2
hence from (1), θ = f ( p) e − kp t .
Taking the inverse Fourier transform, we have
1 ∞ ~ 2
f ( p) e − kp t e − ipx dp
θ ( x, t) =
√ (2 π) − ∞ ∫
1 ∞ ~ 2
f ( p) e − kp t − ipx
or θ ( x, t) =
√ (2 π) ∫− ∞ dp.
T-157
∂4V ∂2V
Example 6: Solve 4
+ = 0 , − ∞ < x < ∞, y ≥ 0
∂x ∂ y2
2 ∞ ~
ip
∂ V e ipx
∞ ∂2V ipx d2
V
dx +
or −
√ (2 π) ∂x2
−∞
− ip
∫ −∞ ∂x2 e
dy2
=0 ,
∂3V
since → 0 as x → ± ∞
∂x3
∞ ~
(ip)2 ∂V ipx ∞ ∂V ipx d2 V
or
√ (2 π) ∂x
e
−∞
− ip
∫ −∞ ∂x
e dx +
dy2
= 0,
∂2V
since → 0 as x → ± ∞
∂x2
~
(ip)3 ∞ ∞ d2 V
or − Ve ipx
( )− ∞ − ip∫ −∞ V e ipx
dx + =0
√ (2 π) dy2
~
(ip)4 ∞
ipx d2 V
or
√ (2 π) ∫ −∞ V e dx +
dy2
=0
~
4
~ d2 V
or p V + =0
dy2
whose solution is
~
V = A cos p2 y + B sin p2 y. …(1)
∂V
Since on y = 0, V = f ( x) and = 0,
∂y
T-158
∞ ~
2 ∂U 2 ∞ d2 U c
or
π ∂x
⋅ cos px + p
0 π ∫0 sin px dx +
dy2
=0
∂U ∂U
→ 0 as x → ∞ and → 0 as x → 0 . By symmetry
∂x ∂x
T-159
~
2 ∞ 2 ∞ d2 U c
p [U sin px)] 0 − p2
or
π π ∫0 U cos px dx +
dy2
=0
~
d2 U c ~
or 2
− p2 U c = 0 (if U → 0 as x → ∞)
dy
~
whose solution is U c = Ae py + Be − py. …(2)
~
Now as y → ∞, U c → 0 . ∴ A = 0 .
~
∴ U c = Be − py. …(3)
Also when y = 0 ,
~ a ∞
2 2
Uc =
π ∫0 1⋅ cos px dx +
π ∫0 0 ⋅ cos px dx
a
2 sin px
=
π p
0
2 sin pa
= ⋅
π p
2 sin pa
∴ from (3), = B.
π p
~ 2 sin pa − py
Hence U c = ⋅e .
π p
1 ∞ e − py
=
π ∫0 p
⋅ [sin (a + x) p + sin (a − x) p] dp
1 ∞ 1 ∞
e − py ⋅ p−1 sin (a + x) p dp + e − py ⋅ p−1 sin (a − x) p dp
=
π ∫0 π ∫0
1 −1 a + x a − x
= tan + tan−1
π y y
∞ r
e − sx x −1 sin rx dx = tan−1
assuming the result
∫0 s
,r>0 ,s>0 .
l l
pπx pπ pπx (Integrating by parts)
= U ( x, t) sin
l 0
−
l ∫0 U cos l
dx
∂U pπ
or Fs =− Fc {U }
∂x l
l pπx
∂U ∂u
and Fc =−
∂x ∫0 ∂x cos l
dx
l l
pπx pπ pπx
= U ( x, t) cos
l 0
+
l ∫0 U sin l
dx
∂U pπ
or Fc = Fs {U } − {U (0 , t) − U (l, t) cos pπ} .
∂x l
∂2U
(b) The finite Fourier sine and cosine transforms of where U is a function of x and t, for
∂x2
0 < x < l, t > 0 .
l
∂2U l ∂2U pπx ∂U pπx pπ l ∂U pπx
Fs 2 =
∂x
∫0 ∂x2 sin
l
dx =
∂x
sin
l 0
−
l ∫0 ∂x cos l
dx
l
pπ pπx pπ l pπx
=−
l
U cos
l
+
0 l ∫0 U ⋅ sin
l
dx
pπ pπ
=− Fs {U } − {U (0 , t) − U (l, t)cos pπ}
l l
∂2U p2 π2 pπ
or Fs 2 = − Fs (U ) + {U (0 , t) − U (l, t) cos pπ}
∂x l l
l
∂2U l ∂2U pπx ∂U pπx pπ l ∂U pπx
and Fc 2 =
∂x
∫0 ∂x2 cos
l
dx =
∂x
cos
l 0
+
l ∫0 ∂x sin l
dx
pπ ∂U
= Fs − {U x (0 , t) − U x (l, t) cos pπ}
l ∂x
T-161
∂2U p2 π2
or Fc 2 = − 2 Fc {U } − {U x (0 , t) − U x (l, t) cos pπ}
∂x l
Solution: Taking the finite Fourier sine transform (with l = 4) of both sides of the given
partial differential equation, we have
4 ∂U pπx 4 ∂2U pπx
∫0 ∂t
sin
4
dx =
∫0 2
∂x
sin
4
dx
~
dUs p2 π2 ~ pπ
or = − 2 Us+ [U (0 , t) − U (4, t) cos pπ] [See article 6.4 (b)
dt 4 4
~
where U s is finite Fourier sine transform of U]
~
dUs p2 π2 ~
or =− Us
dt 16
~ 2 2
whose solution is U s = Ae − p π t /16
. …(1)
Since U ( x, 0 ) = 2 x, where 0 < x < 4 , taking finite Fourier sine transforms, we have
~ 4 pπx
at t = 0 , U s =
∫0 2 x ⋅ sin 4
dx
4
4 pπx 4 4 pπx
= − 2 x ⋅ cos +2⋅ ⋅ sin
pπ 4 pπ pπ 4 0
T-162
Solution: Taking the finite Fourier cosine transform (with l = π) of both the sides of the
given equation, we have
π ∂V π ∂2V
∫0 ∂t
cos px dx = k
∫0 ∂x2
cos px dx
~
dV c ~
or = k [− p2 V c − {V x (0 , t) − V x (π, t) cos pπ}],
dt
~
where V c is finite Fourier cosine transform of V [See article 6.4 (b)]
~
dV c ~
or = − k p2 V c
dt
whose solution is
~ 2
V c = A e − kp t . …(1)
Now taking finite Fourier cosine transform of initial condition, we have
T-163
~ π
at t = 0, V c =
∫0 f ( y) cos py dy .
π
∴ from (1), A =
∫0 f ( x) cos px dx .
~ 2 π
Hence V c = e − kp t
∫0 f ( y) cos py dy .
π ∞ π
1 2 2
=
π ∫0 f ( y) dy +
π ∑= e−kp t cos px. ∫0 f ( y) cos py dy ⋅
p 1
∂U ∂2U
Example 10: Solve = 2 , 0 < x < 6, t > 0
∂t ∂x
subject to the conditions
1, 0 < x < 3
U (0 , t) = 0 , U (6, t) = 0 , U ( x, 0 ) =
0 , 3 < x < 6
and interpret physically. (Kanpur 2011)
Solution: Taking the finite Fourier sine transform (l = 6) of both the sides of the given
partial differential equation, we have
6 ∂U pπx 6 ∂2U pπx
∫0 ∂t
⋅ sin
6
dx =
∫0 ∂x 2
sin
6
dx
~
dUs p2 π2 ~ pπ
or = − 2 Us + [U (0 , t) − U (6, t) cos pπ]
dt 6 6
~
dUs p2 π2 ~
or =− U s,
dt 36
whose solution is
~ 2 2
U s = Ae − p π t /36
. …(1)
1, 0 < x < 3
Now at t = 0, U =
0 , 3 < x < 6.
~ 6 pπx
∴ at t = 0, U s =
∫0 U ⋅ sin 6
dx
3 pπx 6 pπx 6 pπ
=
∫0 1⋅ sin 6
dx +
∫3 0 ⋅ sin 6
dx =
pπ
1 − cos
2
.
T-164
6 pπ
∴ from (1), when t = 0 , 1 − cos = A.
pπ 2
~ 6 pπ − p2 π 2 t /36
Hence Us = 1 − cos e .
pπ 2
~
dUs ~
or = − p2 U s + p [U (0 , t) − U (π, t) cos pπ]
dt
~
dUs ~
or = − p2 U s + p [1 − 3 cos pπ]
dt
~
dUs ~
or + p2 U s = p (1 − 3 cos pπ) …(1)
dt
which is a linear diff. equation of first order.
p2 dt 2
I. F. = e ∫ = ep t .
∴ solution of (1) is
~ 2
p2 t
U s⋅ ep t
= A+
∫ p (1 − 3 cos pπ) e dt
(1 − 3 cos pπ) p2 t
= A+ e
p
T-165
~ 2 (1 − 3 cos pπ)
or U s = Ae − p t
+ ⋅ …(2)
p
Now when t = 0 , U = 1 .
Taking finite Fourier sine transform
~ π 1 ⋅ (1 − cos pπ)
Us =
∫0 1 ⋅ sin px dx =
p
⋅
π
x π cos px
= − cos px +
p 0 ∫0 p
dx
π π
=− cos pπ = (− 1)p + 1.
p p
∞
2 1 − (− 1)p
∴ 1=
π ∑ p
sin px
p =1
∞ ∞
2 π (− 1)p +1 (− 1)p +1
and x=
π ∑ p
sin px = 2 ∑ p
sin px.
p =1 p =1
∞ ∞
2 (1 − 3 cos pπ) 2 1 − 3 (− 1)p
∴
π ∑ p
sin px =
π ∑ p
sin px
p =1 p =1
T-166
∞ ∞
2 1 − (− 1)p 4 (− 1)p +1
=
π ∑ p
sin px +
π ∑ p
sin px
p =1 p =1
2
= 1+ x.
π
Hence from (3),
∞
4 cos pπ − p2 t 2x
U ( x, t) =
π ∑ p
e sin px + 1 +
π
⋅
p=
1
Example 12: Using finite Fourier transform, find the solution of the wave equation
∂2U ∂2U
=4 ,
∂t2 ∂x2
subject to the conditions
U (0 , t) = 0 , U (π, t) = 0 , U ( x, 0 ) = (0 ⋅1) sin x + (0 ⋅ 01) sin 4 x
and U t ( x, 0 ) = 0 for 0 < x < π, t > 0 . (Meerut 2007)
Solution: Taking the finite Fourier sine transform of both the sides of the given
equation, we have
π ∂2U π ∂2U
∫0 ∂t2
sin px dx = 4
∫0 ∂x2
sin px dx
~
d2 U s ~
or = − 4 p2 U s + 4 p [U (0 , t) − U (π, t) cos pπ]
dt2
~
d2 U s ~
or + 4 p2 U s = 0
dt2
whose solution is
~
U s = A cos 2 pt + B sin 2 pt. …(1)
2 π
=
π
(0 ⋅ 1)
∫0 sin x ⋅ sin x dx cos 2 t sin x
2 π
+
π
(0 ⋅ 01)
∫0 sin 4 x ⋅ sin 4 x dx cos 8 t sin 4 x
(Since the integrals in first summation are all zero if p ≠ 1
and the integrals in second summation are all zero if p ≠ 4)
or U ( x, t) = (0 ⋅ 1) cos 2 t sin x + (0 ⋅ 01) cos 8 t sin 4 x.
Comprehensive Exercise 1
∂U ∂2U
1. Solve = 2 , x > 0 , t > 0 subject to the conditions
∂t ∂x
1, 0 < x < 1
(i) U = 0 , when x = 0 , t > 0 (ii) U = when t = 0
0 , x ≥ 1
(iii) U ( x, t) is bounded.
2. If the function U ( x, y) is determined by the differential equation
∂U ∂2U
= for x ≥ 0, − ∞ < y < ∞ and U = f ( y) when x = 0,
∂x ∂ y2
1 ∞ ~ 2
f ( p) e − p x − ipydp
show that U ( x, y) =
2π − ∞ ∫
~
where f ( p) is the Fourier transform of f ( y).
[Hint: Replace k → 1, θ → U , t → x and x → y in Ex. 5]
T-168
A nswers 1
2 ∞1− cos p − p2 t
1. U ( x, t) =
π ∫0 p
e sin px dp
∞ cos pπ − 1 − p2 π 2 t /36
24 pπx
π2 ∑
3. V ( x, t) = 6 + e cos
p2 6
p =1
∞
(1 ⋅ 6) 1 3 (2 n − 1) πt (2 n − 1) πx
4. U ( x, t) =
π 3 ∑
n =1
(2 n − 1)3
cos
2
sin
2
∞
4V0 sinh (2 n + 1) y sin (2 n + 1) x
5. V ( x, y) =
π ∑
n=0
(2 n + 1) sinh (2 n + 1) π
∞
2g 1
6. U ( x, t) = 2 ∑ 3
[1 − (− 1)p ] (1 − cos pat) sin px .
πc p =1
p
¨
T-169
7
F ourier S eries
= bn π.
T-171
1 2π
∴ bn = ∫0 f ( x) sin nx dx. …(4)
π
By taking n = 1, 2 , … we get the values of b1, b2 , …… .
1
Note. To get similar formula of a0 , has been written with a0 in Fourier series.
2
∞
If we write the Fourier series as f ( x) = a0 + Σ [an cos nx + bn sin nx],
n =1
At the point of discontinuity, x = α, the Fourier series gives the value of f ( x) as the
arithmetic mean of left and right limits.
1
∴ At x = α, f ( x) = [ f (α − 0 ) + f (α + 0 ) ].
2
Comprehensive Exercise 1
A nswers
π 1 1 1 1 1
1. e x = − cos x − cos 2 x + cos 3 x − cos 4 x + ...
2 sinh π 2 2 2 10 17
1 2 3 4
+ sin x − sin 2 x + sin 3 x − sin 4 x +...
2 5 10 17
π 4 cos x cos 3 x sin x sin 3 x
2. f ( x) = − 2 + +... + 4 + +. ..
2 π 1 32 1 3
−2 π
1− e 1 1 1 1
3. e− x = + cos x + cos 2 x + cos 3 x +...
π 2 2 5 10
1 2 3
+ sin x + sin 2 x + sin 3 x + ...
2 5 10
1 2
4. f ( x) = 2 sin x − 2 sin 2 x + sin 3 x + … .
π 3
π x, 0 ≤ x ≤1
Example 4: Obtain Fourier Series for the function f ( x) =
π(2 − x), 1 ≤ x ≤ 2.
2
π 1 1 1
Deduce that = + + + ….
8 12 32 52
Solution: Here the function is defined in the interval (0 , 2). Let
a ∞
f ( x) = 0 + Σ (an cos n π x + bn sin nπx) . …(1)
2 n =1
1 2c
Then, we have a0 = f ( x) dx , where c = 1
c ∫0
1 2 1 2
f ( x) dx = πx dx + π (2 − x) dx
1 ∫0 ∫0 ∫1
=
1 2
x2 x2 π π
= π + π 2 x − = + = π;
2 0 2
1
2 2
1 2c nπ x 1 2
an = ∫ f ( x) cos dx = ∫ f ( x) cos nπ x dx
c 0 c 1 0
1 2
= ∫0 π x cos nπ x dx + ∫1 π (2 − x) cos nπ x dx
1 2
sin nπx cos nπx sin nπx cos nπx
= π x + 2
+ π − (2 − x) −
nπ (nπ) 0 nπ (nπ)2 1
Comprehensive Exercise 2
l l ∞1 2 nπx
1. Prove that − x = Σ sin , 0 < x < l.
2 π 1 n l
1 + 2 x , − 1 < x < 0
l
2. Obtain Fourier series of the function F ( x) =
2x
1− , 0 < x < l.
l
x + 1 for − 1 < x < 0 ,
3. Given f ( x) =
x − 1 for 0 < x < 1.
Expand f into a Fourier series on (− 1, 1).
4. Obtain the Fourier series expansion of the periodic function of period 1
1 + x, − 1 < x ≤ 0
f ( x) = 2 2
1 1
− x, 0 < x < .
2 2
0 , 0 < x < c
5. Given f ( x) = ; expand f ( x) in a Fourier series of period 2c.
1, c < x < 2 c
A nswers 2
8 1 πx 1 3 πx 1 5 πx
2. F ( x) = 2 2
cos + 2 cos + 2 cos + ...
π l l 3 l 5 l
2 1 1 1
3. f ( x) = − sin πx + sin 2 πx + sin 3 πx + …
π 1 2 3
1 2 cos 2 πx cos 6 πx cos 10 πx
4. f ( x) = + 2 + + +…
4 π 12 32 52
1 2 1 sin πx + 1 sin 3 πx + …
5. f ( x) = −
2 π 1 c 3 c
Example 6: Obtain Fourier’s series for the expansion of f ( x) = x sin x in the interval
− π < x < π. Hence deduce that
π −2 1 1 1
= − + − …⋅
4 1. 3 3 . 5 5 . 7
Solution: Here f ( x) = x sin x being an even function of x. The Fourier series is
∞
ƒ( x) = a0 + Σ an cos nx
n =1
T-180
1 π 2 ∞ π
ƒ( x) = ∫0 ƒ(u) du + Σ cos n x ∫0 ƒ(u) cos nu du,
π π n =1
2 (− 1)n 1 2 0 , n is even
=− 2 − 2 = 2
[ 1 − (− 1)n ] = 2
π n n πn 4 / πn , n is odd ⋅
π 4 cos x cos 3 x cos 5 x
∴ The required series is f ( x) = − + + + +...
2 π 12 32 52
Comprehensive Exercise 3
A nswers 3
8 cos x cos 3 x cos 5 x
1. f ( x) = + + +... .
π2 12 32 52
π2 ∞ cos nx π2
2. x2 = + 4 Σ (− 1)n ; ⋅
3 n =1 n2 12
π2 6 π2 6 π2 6
3. x3 = 2 − − + 3 sin x + − + 3 sin 2 x − − + 3 sin 3 x …
1 1 2 2 3 3
1 4 sin 2 x 6 sin 3 x
4. x cos x = − sin x + 2 − 2 +…
2 2 −1 3 −1
Again for the interval (0 , T ) the Half Range Fourier sine series is
πx 2 πx nπx
f ( x) = b1 sin + b2 sin + … + bn sin + …,
T T T
2 T nπx
where bn = f ( x) sin dx.
T ∫0 T
Thus when f ( x) is an odd function with period T , we have
∞ nπx 2 T nπx
f ( x) = Σ bn sin , where bn = ∫0 f ( x) sin dx.
n =1 T T T
In this case a0 = 0 = an .
Remark: While expanding a function in (0 , T ) as a series of sines or cosines, we only
see if it is an odd or even function of period 2T . It makes no matter if the function is
odd or even or neither.
2
2 nπt 4 nπt
+ (4 − 2 t) sin − (−2) − 2 2 cos
nπ 2 n π 2 1
4 nπ 8 nπ 8
= sin + 2 2 cos − 2 2
nπ 2 n π 2 n π
8 4 nπ 8 nπ
+ 0 − 2 2 cos nπ − sin + 2 2 cos
n π nπ 2 n π 2
8 nπ 8 8 8 cos nπ
= 2 2
cos − 2 2 − 2 2 cos nπ + 2 2
n π 2 n π n π n π 2
T-184
∴ a1 = − 8 / π2 , a2 = 0 , a3 = − 8 / 32 π2 , a4 = 0 , a5 = − 8 / 52 π2 , etc.
Hence, the required half-range Fourier series for f ( x) in (0, 2) is
8 cos πx /2 cos 3 πx /2 cos 5 πx /2
f ( x) = 1 − 2 + + +…
π 1 2
3 2
5 2
Comprehensive Exercise 4
x , for 0 ≤ x≤ π /2,
4. Given f ( x) =
π − x , for π / 2 ≤ x ≤ π.
Express this function by a sine series and also by a cosine series.
5. Develop sin ( πt / l ) in half range cosine series in the range 0 < t < l.
6. Expand f ( x) = e x in a cosine series over (0 , 1).
1 − x, if 0 < x < 1 ,
7. Expand f ( x) = 4 2
3 1
x − , if < x < 1.
4 2
A nswers 4
1 4 1 πx 1 3 πx 1 5 πx
1. f ( x) = + cos + 2 cos + 2 cos … ⋅
2 π2 12 l 3 l 5 l
− e − 1 e −1 − e −1
6. ex = e − 1+ 2 2 cos πx + 2
cos 2 πx + 2
cos 3 πx + …
π +1 4 π +1 9 π +1
1 4 1 4 1 − 4 sin 5 π x + …
7. f ( x) = − 2 sin π x + + sin 3 π x +
π π 3 π 32 π2 5 π 52 π2
∞ 1 π4 ∞ 1 π
(c) Σ 4
= . (d) Σ 6
= ⋅
n =1 n 90 n =1 n 945
Solution: Let f ( x) = x (π − x), 0 < x < π.
(a) The half range cosine series for f ( x) is
a ∞
f ( x) = 0 + Σ an cos nx.
2 n =1
π
2 π 2 πx2 x3 2 π3 π3 π2
Here a0 = ∫0 x (π − x) dx = − = − =
π π 2 3 π2 3 3
0
2 π
and an = ∫0 x (π − x) cos nx dx
π
π
2 sin nx − cos nx − sin nx
(πx − x2 )
= − (π − 2 x) + (− 2)
π n n2 n3 0
2 π (−1)n π 2 π
= 0 − 2
+ 0 − 2 = 2 [− (−1)n − 1]
π n n π n
−4 / n2 , when n is even
=
0 , when n is odd.
π2 cos 2 x cos 4 x cos 6 x
∴ x (π − x) = −4 + + +… ⋅
6 22 42
6 2
2 π 1 π4 1 1 1
or ∫0 (π2 x2 − 2 πx3 + x4 ) dx = + 16 4 + 4 + 4 + …
π 2 9 2 4 6
π
2 π2 x3 2 πx4 x5 π4 1 1 1
or − + = + 4 + 4 + 4 + …
π 3 4 5
18 1 2 3
0
2 x5 2 π5 π5 π4 1 1 1
or − + = + + + + …
π 3 4 5 18 14 24 34
π4 π4 ∞ 1 ∞ 1 π4
or = + Σ 4
or Σ 4
= ⋅
15 18 n =1 n n =1 n 90
∞
(b) The half range sine series for f ( x) is f ( x) = Σ bn sin nx .
n =1
2 π
Here bn = ∫0 x (π − x) sin nx dx
π
π
2 − cos nx − sin nx cos nx
= (πx − x2 ) − (π − 2 x) + (− 2)
π n n2 n3 0
8 , when n is odd
2 (−1)n 2 4 n
= − 2 3 + 3 = 3 [− (−1) + 1] = n3 π
π n n πn 0 , when n is even.
T-188
π2 64 1 1 1 π4 1 1 1
or = 2 6 + 6 + 6 + … or = 6 + 6 + 6 ⋅
15 π 1 3 5
960 1 3 5
1 1 1 1 1 1 1 1 1 1
Let S = 6 + 6 + 6 + 6 + … = 6 + 6 + 6 + … + 6 + 6 + 6 + …
1 2 3 4 1 3 5 2 4 6
4 4 4
π 1 1 1 π 1 1 1 1 π S
= + 6 + 6 + 6 + … = + 6 6 + 6 + 6 + … = + ⋅
960 2 4 6 960 2 1 2 3
960 64
S π4 63 S π4
∴ S− = or = ⋅
64 960 64 960
π4 64 π4 ∞ 1 π4
or S= × = or Σ 6
= ⋅
960 63 945 n =1 n 945
π2 1 π2 16 1
or − π2 0 − = + 2 1+ + …
3 3 2
π 81
2
2π π2 16 1 1
or − = 2 1 + 4 + 4 + …
3 2 π 3 5
π2 16 1 + 1 + 1 + … π4 1 1
or = 2 or = 1+ 4 + 4 + ….
6 π 34 54 96 3 5
This is the complex form of Fourier series and its coefficients are given by (3).
If the function is defined in the interval (− T , T ), then the coefficients are given by
1 T
cn = f ( x) e − i n πx /T dx. (n = 0 , ± 1, ± 2,...)
2T ∫ −T
The complex form of a Fourier series is especially useful in problems on electrical
circuits having periodic voltage.
Example 14: Obtain the complex form of the Fourier series of the function
0, − π≤ x≤0
f ( x) =
1, 0 ≤ x ≤ π.
Solution: The complex form of Fourier series for f ( x) is
∞
f ( x) = Σ c n e inπ x /T
.
n= − ∞
1 T 1 π
Here cn = ∫ −T f ( x) e − inπ x /T
dx = ∫−π f ( x) e − inx dx
2T 2π
1 0 π
0 . e − inx dx + 1 . e − inx dx
2 π ∫ − π ∫0
=
π
1 π − inx 1 e − inx 1
= ∫0 e dx = =− [e − inπ − 1]
2π 2 π − in 2 nπ i
0
1 , n is odd
1 1
=− [cos nπ − 1] = − [(−1)n − 1] = inπ
2 nπ i 2 nπ i 0 , n is even
1 T 1 π 1 π 1 1
and c0 =
2T ∫ −T
f ( x) dx =
2π ∫ −π
f ( x) dx =
2π ∫0
dx =
2π
[ x ]0π = 2 ⋅
1 1 e ix e 3 ix e 5 ix 1 e − ix e − 3 ix e − 5 ix
∴ f ( x) = + + + + … + + + + …
2 iπ 1 3 5 i π − 1 − 3 − 5
1 1 ix 1 1
= − (e − e
− ix
) + (e3 ix − e − 3 ix ) + (e5 ix − e −5 ix ) + … ⋅
2 iπ 3 5
Comprehensive Exercise 5
A nswers 5
4 4 4 4
1. 1= sin x + sin 3 x + sin 5 x + sin 7 x + …
π 3π 5π 7π
n2 ∞ 4 ∞ (− 1)n (1 − in π)
2. f ( x) = + Σ 2 (− 1)n cos nx 3. e − x = Σ sinh 1 ⋅ e inπx .
3 n =1 n n= − ∞ 1 + n2 π2
True or False
Write ‘T’ for true and ‘F’ for false statement.
1. A function f ( x) is even if f (− x) = − f ( x).
2. A function f ( x) is odd if f (− x) = f ( x).
3. Most functions are neither even nor odd.
4. A function f ( x) can always be expressed as an arithmetic mean of an even and
1 1
odd function as f ( x) = [ f ( x) + f (− x)] + [ f ( x) − f (− x)].
2 2
5. For an even function f ( x), the graph of f ( x) is symmetric about x-axis.
6. If f ( x) = x2 in (− π, π) then the Fourier series of f ( x) contains only sine terms.
7. The period of f ( x) = sin 2 x is π.
8. For an odd function f ( x), the graph of f ( x) is symmetric about the origin.
A nswers
True or False
1. F 2. F 3. T 4. T 5. F
6. F 7. T 8. T