Gasal2324 Bab 7 PD Order Dua Koef Var

Unduh sebagai pdf atau txt
Unduh sebagai pdf atau txt
Anda di halaman 1dari 21

Bab 7 PD Order 2 Koef.

Variabel

Pada bab ini akan dipelajari mengenai PD order dua yang mempunyai bentuk umum:
d2y dy
2
+ R( x) + S ( x) y = Q( x)
dx dx

Penyelesaian:

A. PD ini diselesaikan dengan merubah variabel tak bebas y,


Misalkan y = u × v , maka:
dy d dv du
= (u × v) = u +v
dx dx dx dx
2
d y d æ dv du ö
2
= çu + v ÷
dx dx è dx dx ø
d 2v
du dv d 2 u dv du
= u 2 + +v 2 +
dx dx dx dx dx dx
2 2
d v du dv d u
= u 2 +2 +v 2
dx dx dx dx
dengan mensubtitusikan ke PD, didapat :
æ d 2v du dv d 2 u ö÷ æ dv du ö
çu +2 + v 2 + R( x)ç u + v ÷ + S ( x)u × v = Q( x)
ç dx 2 dx dx dx ø ÷ è dx dx ø
è
d 2 v æ du ö dv æç d u ö
2
du ÷ v = Q( x)
u 2 + ç2 + uR( x) ÷ + + R ( x ) + S ( x )u
dx è dx ø dx çè dx 2 dx ÷
ø
d 2 v æ 2 du ö dv æç 1 æç d u öö
2
du ÷ ÷ v = Q( x)
+ ç + R ( x ) ÷ + + R ( x ) + S ( x )u
2
èudx ç ç 2 ÷÷
dx  ø dx èu è dx

dx øø u

R1 ( x ) Q ( x)
1
S1 ( x )

d 2v dv
+ R1 ( x) + S1 ( x) v = Q1 ( x)
dx 2 dx

Selanjutnya, kemungkinan-kemungkinan yang terjadi antara lain:


1. Jika u solusi partikelir dari PD yang diberikan semula, maka PD-nya menjadi:
d 2v dv
2
+ R1 ( x) = Q1 ( x)
dx dx
dv d 2 v dp
dan dengan permisalan = p dan = , didapat PD order 1
dx dx 2 dx
dp
+ R1 ( x) p = Q1 ( x)
dx

7-1
Bab 7 PD Oredr 2 Koef. Variabel

2. Jika u dipilih sedemikian sehingga R1 ( x) = 0 , maka:


2 du
+ R( x) = 0
u dx
du R( x)dx
=-
u 2
1
ln u = - ò R( x) dx
2
1
- ò R ( x ) dx
u= e 2
1
du 1 1 - ò R ( x ) dx
= - uR( x) = - e 2 R( x)
dx 2 2
d 2u 1 dR( x) 1 du
2
=- u - R( x)
dx 2 dx 2 dx

sehingga dengan memsubtitusi ke S1 ( x) , didapat :


1 æç d 2 u du ö
÷
S1 ( x) = + R ( x ) + S ( x )u
u çè dx 2 dx ÷
ø
1 d 2u 1 du
= 2
+ R( x) + S ( x)
u dx u dx
1 du 1 d 2 u
= S ( x) + R( x) +
u dx u dx 2
1 æ 1 ö 1 æ 1 dR( x) 1 du ö
= S ( x) + R( x)ç - uR( x) ÷ + ç - u - R( x) ÷
u è 2 ø uè 2 dx 2 dx ø
1 1 dR( x) 1 du
= S ( x) - R 2 ( x) - - R( x)
2 2 dx 2u dx
1 1 dR( x) 1 æ 1 ö
= S ( x) - R 2 ( x) - - R( x)ç - uR( x) ÷
2 2 dx 2u è 2 ø
1 2 1 dR( x) 1 2
= S ( x) - R ( x) - + R ( x)
2 2 dx 4
1 1 dR( x)
= S ( x) - R 2 ( x) -
4 2 dx
1 1 dR
\ S1 = S - R 2 -
4 2 dx
> Jika S1 = A , maka PD menjadi :
d 2v
+ Av = Q1 ( x) (merupakan PDlinier order 2 koefisien konstanta)
dx 2
A
> Jika S1 = , maka PD menjadi :
x2
d 2v
x2 + Av = Q1 ( x) (merupakan PD Cauchy)
dx 2

7-2
Bab 7 PD Order 2 Koef. Variabel

Teorema 7.1:
Misalkan diketahui PD order dua sebagai berikut:
( D 2 + RD + S ) y = 0
Maka :
a). Jika R + Sx = 0 , maka y = x adalah solusi komplementer.
b). Jika 1 + R + S = 0 , maka y = e x adalah solusi komplementer.
c). Jika 1 - R + S = 0 , maka y = e - x adalah solusi komplementer.
d). Jika m 2 + mR + S = 0 , maka y = e mx adalah solusi komplementer.

Bukti :
a). ( D 2 + RD + S ) y = 0
( D 2 + RD + S ) x = 0
D 2 x + RDx + Sx = 0
0 + R + Sx = 0
R + Sx = 0

b). ( D 2 + RD + S ) y = 0
( D 2 + RD + S )e x = 0
D 2 e x + RDe x + Se x = 0
e x + R e x + Se x = 0
e x (1 + R + S ) = 0
karena e x ¹ 0 , maka 1 + R + S = 0

c). ( D 2 + RD + S ) y = 0
( D 2 + RD + S )e -x = 0
D 2 e - x + RDe - x + Se - x = 0
e - x - R e - x + Se - x = 0
e - x (1 - R + S ) = 0
karena e x ¹ 0 , maka 1 - R + S = 0

d). ( D 2 + RD + S ) y = 0
( D 2 + RD + S )e mx = 0
D 2 e mx + RDe mx + Se mx = 0
m 2 e mx + Rme mx + Se mx = 0
e mx (m 2 + mR + S ) = 0
karena e x ¹ 0 , maka m 2 + mR + S = 0

7-3
Bab 7 PD Oredr 2 Koef. Variabel

Contoh 7.2:
1. Selesaiakan PD berikut:
xy ¢¢ - ( x + 2) y ¢ + 2 y = 0

Penyelesaian:

Dengan membagi x pada setiap sukunya, didapat:


( x + 2) 2
y ¢¢ - y¢ + y = 0
x x
( x + 2) 2 2 2
Þ 1+ R + S = 1- + = 1-1- + = 0
x x x x
x
\ y = e adalah solusi komplementer.

Transformasi y = u × v = e x v sehingga dengan subtitusi ke PD, didapat :


d 2v dv
2
+ R1 ( x)
= Q1 ( x)
dx dx
Q( x) 0
dengan : Q1 ( x) = = =0
u u

2 du ( x + 2) 2 x 2 2
R1 ( x) = R( x) + =- + x e = -1 - + 2 = 1 -
u dx x e x x
d 2v æ 2 ö dv
\ + ç1 - ÷ = 0
dx 2 è x ø dx

dv d 2 v dp
Misalkan : = p , maka = , sehingga:
dx dx 2 dx
dp æ 2 ö
+ ç1 - ÷ p = 0
dx è x ø
dp æ 2 ö
= ç - 1÷dx
p èx ø
ln p = 2 ln x - x + C
x 2C
p=
ex

Selanjutnya, langkah dilanjutkan dengan mengembalikan pemisalan p:


dv
= Cx 2 e - x
dx
v = ò Cx 2 e - x dx
y
x
= ò Cx 2 e - x dx
e
y = e x ò Cx 2 e - x dx

7-4
Bab 7 PD Order 2 Koef. Variabel

2. Selesaikan PD berikut:
x 2 y ¢¢ - x(2 x + 3) y ¢ + ( x 2 + 3 x + 3) y = (6 - x 2 )e x

Penyelesaian:
Dengan setiap sukunya dengan x2, didapat:
2x + 3 x 2 + 3x + 3 6 - x2 x
¢¢
y - ¢
y + y= e
x x2 x2
Jika u dipilih sedemikian sehingga R1 ( x) = 0 , maka:
1 1 dR
S1 ( x) = S - R 2 -
4 2 dx
2x + 3 (2 x + 3) 2 dR 3
dengan R = - , R2 = 2
dan = 2
x x dx x
x 2 + 3x + 3 (2 x + 3) 2 3 3
\ S1 ( x) = 2
- 2
- 2
=-
x 4x 2x 4x 2
æ A 3ö
ç S1 = 2 ; A=- ÷
è x 4ø

Sehingga PD menjadi :
d 2v dv
+ R1 ( x ) + S1 ( x)v = Q1 ( x) ; R1 ( x) = 0
dx 2 dx
d 2v 3 1 6 - x2
- v=
; ex
dx 2 4 x 2 x2 u
dengan u sedemikian sehingga R1 ( x) = 0 ;
2 du
R( x) + =0
u dx
du 2
= - R( x)
u dx
du 1
= - R( x)dx
u 2
1
ln u = - ò R( x)dx
2
1 1 3 1
- ò R ( x ) dx - ò ( -2- ) dx - ( -2 x -3 ln x ) 32
u= e 2 = e 2 x = e 2 = e x +ln x = e x × x3 2

d 2v 3 6 - x2 1
2
- 2
v= 2
ex
dx 4x x e × x3 2
x

d 2v 3
2
- 2
v = (6 - x 2 ) x -7 2
dx 4x
2
d v 3
x 2 2 - v = ( 6 - x 2 ) x -3 2 (PD Linear Cauchy order 2)
dx 4

7-5
Bab 7 PD Oredr 2 Koef. Variabel

Selanjutnya, misalkan x = e z , maka:


dv d 2v
x = Dv dan x 2 2 = D( D - 1)v
dx dx
Sehingga dengan mensubstitusikan pemisalan ini ke PD Cauchy, didapat:
3
D( D - 1)v - v = (6 - e 2 z )e -3 z 2
4
(4 D - 4 D - 3)v = 4(6e -3 z 2 - e z 2 )
2


dan seterusnya dapat dilanjutkan sendiri sebagai latihan.

3. Selesaikan PD berikut:
d 2 y 2 dy æ 2 ö
2
- + ç1 + 2 ÷ y = xe x
dx x dx è x ø

Penyelesaian:
Jika u dipilih sedemikian sehingga R1 ( x) = 0 , maka:
1 1 dR
S1 ( x) = S - R 2 -
4 2 dx
2 4 dR 2
dengan R = - , R 2 = 2 dan = 2 , maka:
x x dx x
2 1 1
S1 ( x) = 1 + 2 - 2 - 2 = 1
x x x

Sehingga PD menjadi:
d 2v dv
2
+ R1 ( x) + S1 ( x)v = Q1 ( x) ; R1 ( x) = 0
dx dx
1 1 2 dx
d 2v x 1 - ò R ( x ) dx - ò - dx ò
2
+ v = xe ; u = e 2 = e 2 x = e x = e ln x = x
dx u
d 2v
2
+ v = e x (P Linear order 2 koefisien konstanta)
dx
( D 2 + 1)v = e x
dan seterusnya dapat dilanjutkan sendiri sebagai latihan.

B. Dengan merubah variabel bebas


Pandang PD berikut:
d2y dy
2
+ R( x) + S ( x) y = Q( x)
dx dx
dy dy dz
Misalkan z = θ (x) dan juga = , sehingga:
dx dz dx

7-6
Bab 7 PD Order 2 Koef. Variabel

d2y d æ dy ö d æ dy dz ö
= ç ÷ = ç ÷
dx 2 dz è dx ø dz è dz dx ø
d2y 1 dy dz dy d 2 z
= +
dx 2 dz dz dx dz dx 2
d dz dy dz dy d 2 z
= +
dz dx dz dx dz dx 2
2
d 2 y æ dz ö dy d 2 z
= ç ÷ +
dz 2 è dx ø dz dx 2
Subtitusikan pemisalan di atas ke PD, sehingga didapat:
æ d 2 y æ dz ö 2 dy d 2 z ö
ç ç ÷ + ÷ + R( x) dy dz + S ( x) y = Q( x)
ç dz 2 è dx ø dz dx 2 ÷ø dz dx
è
2 2
æ dz ö d y æç d z dz ö÷ dy
2
ç ÷ + + R ( x ) + S ( x) y = Q( x)
è dx ø dz 2 çè dx 2 dx ÷ø dz
æ d 2z dz ö÷
ç + R ( x )
2
d y ç dx 2 dx ÷ dy S ( x) Q( x)
2
+ç 2 ÷ dz + 2
y= 2
dz ç æ dz ö ÷ æ dz ö æ dz ö
ç ç ÷ ÷ ç ÷ ç ÷
è è dx ø ø è dx ø è dx ø

2
æ dz ö ±S
Pilih z = θ (x) sedemikian sehingga ç ÷ = 2
è dx ø a
2
d z dz
2
+ R( x)
• Jika dx dx
2
= A ; dengan A konstanta
æ dz ö
ç ÷
è dx ø
d2y dy 2 a 2Q
Maka PD menjadi + A +a y =
dz 2 dz S
2
d z dz
2
+ R ( x )
dx
• Jika dx 2
¹ A , maka cari cara lain
æ dz ö
ç ÷
è dx ø
Contoh 7.3:
Selesaikan PD berikut:
3
xy¢¢ - 3 y¢ + y = x + 2
x
Penyelesaian:
Dengan membagi setiap sukunya dengan x, maka didapat:
3 3 x+2
y¢¢ - y¢ + 2 y =
x x x

7-7
Bab 7 PD Oredr 2 Koef. Variabel

3 3 x+2
R( x) = - , S ( x) = 2 , Q( x) =
x x x
Pilih z = θ (x) sedemikian sehingga:
2
æ dz ö ±S
ç ÷ = 2
è dx ø a
3
2
æ dz ö 2 1 dz 1 d 2z 1
dengan ç ÷ = x = 2 , = , dan 2
=- 2
è dx ø 3 x dx x dx x
d 2z dz 1 3 1
2
+ R( x) - 2- ×
dx dx = x x x = -4
æ dz ö
2 1
ç ÷ x2
è dx ø
maka didapat:
d2y dy Q
2
+ A + a2 y = 2
dx dx æ dz ö
ç ÷
è dx ø
2
d y dy x+2 2
2
- 4 + 3y = ×x
dx dx x
( D 2 - 4 D + 3) y = x 2 + 2 x (PDL order 2 koef. Konstanta),
dan seterusnya kerjakan sendiri sebagai latihan.

C. Dengan Faktorisasi.

Contoh 7.4:
Selesaikan PD berikut:
d2y dy
( x + 1) 2 - (3 x + 4) + 3 y = (3 x + 2)e3 x
dx dx

Penyelesaian:
PD di atas dapat ditulis sebagai berikut:
( )
( x + 1) D 2 - (3 x + 4) D + 3 y = (3 x + 2)e3 x
(( x + 1) D - (3x + 3 + 1) D + 3)y = (3x + 2)e
2 3x

(( x + 1) D - 3( x + 1) - D + 3)y = (3x + 2)e


2 3x

(( x + 1)(D - 3D) - ( D - 3))y = (3x + 2)e


2 3x

(( x + 1)( D - 3) D - ( D - 3) ) y = (3x + 2)e3 x


(( x + 1) D - 1)( D - 3) y = (3x + 2)e3 x
misalkan ( D - 3) y = v , maka didapat:
( x + 1) Dv - v = (3 x + 2)e3 x

7-8
Bab 7 PD Order 2 Koef. Variabel

1
Dv - v = (3 x + 2)e3 x (PD linier order satu),
x +1
dan seterusnya kerjakan sendiri sebagai latihan.

PD CAUCHY-EULER ORDER 2

Yang dimaksud dengan PD Cauchy – Euler adalah PD yang mempunyai bentuk umum:
2
2d y dy
P0 x 2
+ P1x + P2 y = α ( x)
dx dx
dengan P0 , P1, P2 konstanta.

Penyelesaian:
dz 1
Misalkan x = e z , maka z = ln x dan =
dx x
d
Definisikan D = , sehingga didapat:
dx
dy dy dz 1 dy dy
= = Þ x = Dy
dx dz dx x dz dx
d2y d æ dy ö d æ 1 dy ö
2
= ç ÷= ç ÷
dx dx è dx ø dx è x dz ø
1 dy 1 d 2 y 1 dy 1 d 2 y
=- + = - +
x 2 dz x dxdz x 2 dz x 2 dz 2
1 æ d 2 y dy ö
= 2 çç 2 - ÷÷
x è dz dz ø
d 2 y d 2 y dy æç d 2 d ö
x2 2 = 2 - = - ÷ y = ( D 2 - D) y = D( D - 1) y
dx dz dz çè dz 2 dz ÷ø

subtitusikan ke PD, maka didapat PD linier non homogen koefisien konstanta:

P0 D( D - 1) y + P1Dy + P2 = α ( x)
(P D
0
2
)
+ ( P1 - P0 ) D + P2 y = α ( x)

selanjutnya analog dengan metode yang lalu .....!

Contoh 7.5:

1. Selesaikan PD berikut:
( x 2 D 2 - 3 xD + 4) y = x + x 2 ln x
Penyelesaian:
d
Dengan pendefinisian D = , sehingga didapat:
dx

7-9
Bab 7 PD Oredr 2 Koef. Variabel

x 2 D 2 = D( D - 1)
3xD = 3D
misalkan x = e z , maka z = ln x
(D( D - 1) - 3D + 4) y = e z + ze2 z
Solusi homogen:
D 2 - D - 3D + 4 = 0
D2 - 4D + 4 = 0
( D - 2) 2 = 0
D12 = 2
\ yh = C1e 2 z + C2 ze 2 z

Solusi non homogen :


( D - 2) 2 y = e z + ze 2 z
1 1
y = 2
ez + ze 2 z
( D - 2) ( D - 2) 2
1
= e z + e2 z z
( D + 2 - 2) 2
1
= e z + e2 z 2 z
D
1
y = 2 z Þ D2 y = z
D
1
= e z + e2 z z
6
1 1
Dy = z 2 Þ y = z 3
2 6
Solusi umum :
1 3 1
y = yh + yk = C1e 2 z + C2 ze 2 z + e z + e 2 z z = C1x 2 + C2 x 2 ln x + x + x 2 ln 3 x
6 6

2. Selesaiakan PD berikut:
( x 2 D 2 - 2 xD + 2) y = (ln x) 2 - ln x 2

Penyelesaian:
Misalkan x = e z , maka ln x = z , ln x 2 = 2 z , dan:
x 2 D 2 = D( D - 1)
2 xD = 2 D
sehingga setelah disubstitusikan:
( D( D - 1) D - 2 D + 2) y = z 2 - 2 z
( D 2 - D - 2 D + 2) y = z 2 - 2 z
( D 2 - 3D + 2) y = z 2 - 2 z

7 - 10
Bab 7 PD Order 2 Koef. Variabel

( D - 2)( D - 1) y = z 2 - 2 z

Solusi homogen :
\ yh = C1e 2 z + C2e z

Solusi non homogen :


( D - 2)( D - 1) y = z 2 - 2 z
z2 - 2z
y=
( D - 2)( D - 1)

Deret McLaurint:
1 1
f ( x) = 2
Þ f (0) =
x - 3x + 2 2
3
f ¢( x) = -( x 2 - 3 x + 2) - 2 (2 x - 3) Þ f ¢(0) =
4
7
f ¢¢( x) = 2( x 2 - 3 x + 2) -3 (2 x - 3) 2 - 2( x 2 - 3 x + 2) - 2 Þ f ¢¢(0) =
8
f ¢(0) f ¢¢(0) 2 1 3 7
\ f ( x) = f (0) + x+ x = + x + x2
1! 2! 2 4 8
æ1 3
è2 4
7 ö
yk = ç + D + D 2 ÷ z 2 - 2 z
8 ø
( )
(
1 2 3
) ( 7
)
= z - 2z + D z 2 - 2z + D2 z 2 - 2z
2 4 8
( )
1
( 3
)
= z 2 - 2 z + (2 z - 2 ) +
2 4
14
8
1 1 1
= z2 + z +
2 2 4

Solusi Umum:
1 2 1 1
y = yh + yk = C1e 2 z + C2e z + z + z+
2 2 4

PD LAGRANGE- LEGENDRE ORDER 2

Yang dimaksud dengan PD Lagrange Euler adalah PD yang mempunyai bentuk umum:
d2y dy
P0 (ax + b) 2 2 + P1 (ax + b) + P2 y = Q( x)
dx dx
dengan P0 , P1, P2 konstanta.

Penyelesaian:
dz a
Misalkan (ax + b) = e z , maka z = ln(ax + b) dan =
dx ax + b

7 - 11
Bab 7 PD Oredr 2 Koef. Variabel

d
Definisikan D = , sehingga didapat:
dz
dy dy dz a dy
• = =
dx dz dx ax + b dz
dy dy
(ax + b) =a
dx dz
dy
(ax + b) = aDy
dx
2
d y d æ dy ö
• 2
= ç ÷
dx dx è dx ø
d æ a dy ö
= ç ÷
dx è ax + b dz ø
a2 dy a d2y
=- +
(ax + b) 2 dz ax + b dxdz
a2 dy a d2y
=- +
(ax + b) 2 dz ax + b ax + b (dz ) 2
a
2 2
a dy a d2y
=- +
(ax + b) 2 dz (ax + b) 2 dz 2
a2 æ d 2 y dy ö
= ç - ÷
(ax + b) 2 çè dz 2 dz ÷ø

d2y
(ax + b) 2 2
= a 2 ( D 2 - D) y = a 2 D( D - 1) y
dx

Subtitusikan pemisalan tersebut ke PD, didapat PD linier non homogen koefisien


konstanta :
P0 a 2 D( D - 1) y + P1aDy + P2 y = Q( x)
( P0 a 2 D 2 + ( P1a - P0 a 2 ) D + P2 ) y = Q( x)

Selanjutnya gunakan metode penyelesaian PD linier non homogen koefisien konstanta.

Contoh 7.6:
1. Selesaikan PD berikut:
( )
( x + 1) 2 D 2 + ( x + 1) D - 1) y = ln( x + 1) 2 + x - 1

Penyelesaian:
Misalkan x + 1 = e z , maka x = e z - 1 , ( x + 1) 2 = e 2 z , ln( x + 1) 2 = 2 z
sehingga : ( D( D - 1) + D - 1) y = 2 z + e z - 2

7 - 12
Bab 7 PD Order 2 Koef. Variabel

• Solusi homogen :
D( D - 1) + D - 1 = 0
D2 - D + D - 1 = 0
D2 - 1 = 0
( D + 1)( D - 1) = 0
\ yh = C1e z + C2e - z

• Solusi non homogen :


( D 2 - 1) y = 2 z + e z - 2
1
yk = 2 (2 z + e z - 2)
D -1
1 1 2
= 2 2z + ez - 2
D -1 ( D + 1)( D - 1) D -1
1 1
= (-1 - D 2 )2 z + ez + 2
2 ( D - 1)
1
= -2 z + ze z + 2
2

• Solusi umum :
1 z
y = yh + yk = C1e z + C2e - z - 2 z + ze + 2
2

2. Tunjukkan solusi dari:


(1 + x) 2 y¢¢ - 3(1 + x) y¢ + 4 y = (1 - x)3
adalah
y = C1 ( x + 1) 2 + C2 ( x + 1) 2 ln( x + 1) + ( x + 1)3
Kerjakan sendiri sebagai latihan

Reduksi Order PD Linear dengan Metode Integral

• Misalkan diketahui bahwa y = Φ1 adalah solusi dari PD linier homogen order 2,


maka solusi ini dapat digunakan untuk mereduksi order menjadi order 1, sehingga
dengan menyederhanakan seperlunya, didapat satu solusi lain yang bebas linear.

• Jika solusi umumnya telah ditentukan, maka dapat digunakan untuk mencari solusi
partikelir dengan metode yang sudah dipelajari.

• Untuk jelasnya, perhatikan contoh berikut:

Contoh 7.7:
Tentukan solusi umum dari:

7 - 13
Bab 7 PD Oredr 2 Koef. Variabel

d2y dy
a2 ( x) + a1 ( x ) + a0 ( x) y = 0
dx 2 dx
dengan y1 ( x) adalah solusi yang diketahui.

Jawab:
Misalkan y = uy1 , maka:
dy du dy
= y1 + u 1
dx dx dx
d 2 y d 2u du dy1 du dy1 d 2 y1
= y1 + + + u
dx 2 dx 2 dx dx dx dx dx 2
Subtitusikan pemisalan di atas ke PD, maka didapat:
d 2 y a1 ( x) dy a0 ( x)
+ + y=0
dx 2 a2 ( x) dx a2 ( x)

æ d 2u du dy1 du dy1 d 2 y1 ö÷ a1 ( x) æ du dy ö a ( x)
ç y + + +u 2 + ç y1 + u 1 ÷ + 0 uy1 = 0
ç dx 2 1 dx dx dx dx ÷
dx ø a2 ( x) è dx dx ø a2 ( x)
è
d 2u æ dy a ( x) ö du æç d 2 y1 a1 ( x) dy1 a0 ( x) ö÷
y1 2 + çç 2 1 + 1 y1 ÷÷ + + + y u=0
ç 2 a ( x) dx a ( x) 1 ÷
dx è dx a2 ( x) ø dx èdx ø
2 2
0
karena y1 solusi

d 2u æ dy1 a1 ( x) ö du
y1 + ç2 + y1 ÷ =0
dx 2 çè dx a2 ( x) ÷ø dx
d 2u æ 2 y1¢ a1 ( x) ö du
+ç + ÷ =0
dx 2 çè y1 a2 ( x) ÷ø dx
Selanjutnya misalkan:
du d 2u dv
=v, =
dx dx 2 dx
maka didapat:
dv æ 2 y1¢ a1 ( x) ö
+ç + ÷v = 0
dx çè y1 a2 ( x) ÷ø
æ 2 y1¢ a1 ( x ) ö
ò çç + ÷ dx
÷
faktor integrasi v = e è y1 a2 ( x ) ø
y1¢ a1 ( x )
2ò dx òa dx
y1 2 ( x)
=e ×e
a1 ( x )
2 ln y12 òa dx
=e ×e 2 ( x)

a1 ( x )
ò a ( x ) dx
= y12 × e 2

setelah setiap suku dikalikan dengan v, didapat :

7 - 14
Bab 7 PD Order 2 Koef. Variabel

a1 ( x ) a1 ( x )
òa dx
dv æ 2 y1¢ a1 ( x) ö 2 ò a2 ( x ) dx
y12 × e 2 ( x) +ç + ÷v y e =0
dx çè y1 a2 ( x) ÷ø 1
æ a1 ( x ) ö
ç 2 ò a2 ( x ) dx ÷
d ç v y1 e ÷=0
ç ÷
è ø
a1 ( x )
òa dx
v y12 e 2 ( x) = C1
a1 ( x )
-ò dx
a2 ( x )
e
v = C1
y12
a1 ( x )
-ò dx
a2 ( x )
e
u = C1 ò dx + C2
y12
Kemudian subtitusikan ke solusi umum:
y = uy1

Contoh 7.8:
Tentukan dua solusi bebas linear dari : x 2 y¢¢ + xy¢ - y = 0 jika diberikan Φ1 = x adalah
salah satu solusinya. Selanjutnya, gunakan solusi ini untuk menentukan solusi dari
x 2 y¢¢ + xy¢ - y = 2 x 2 + 1 .

Jawab :
Misal : y = xv , maka y¢ = xv¢ + v dan y¢¢ = xv¢¢ + v¢ + v¢
subtitusikan pemisaklan tersebut ke PD homogen, didapat :
x 2 ( xv¢¢ + v¢ + v¢) + x( xv¢ + v) - xv = 0
x3v¢¢ + 3 x 2v¢ = 0
3
v¢¢ + v¢ = 0 (misal : w = v¢ , w¢ = v¢¢ )
x
3
3 - ò dx -3
w¢ + w = 0 (faktor integrasi : e x = e -3 ln x = eln x = x -3 )
x
w = C1x -3
C1
v= + C2
x2
• Solusi komplementer :
C
yh = xv = 1 + xC2
x

• Solusi partikelir (dengan metode koefisien tak tentu) :


yk = ax 2 + b
yk¢ = 2ax

7 - 15
Bab 7 PD Oredr 2 Koef. Variabel

yk¢¢ = 2a

subtitusikan ke PD non homogen :


2ax 2 + x × 2ax - ax 2 - b = 2 x 2 + 1
3ax 2 - b = 2 x 2 + 1
2
a=
3
b = -1
2
yk = x 2 - 1
3
• Solusi umum :
C 2
y = yh + yk = 1 + xC2 + x 2 - 1
x 3

Latihan 7.10:
Kerjakan seperti contoh 7.9 di atas, untuk soal berikut:
1. ( x + 1) 2 y¢¢ - 2( x + 1) y¢ + 2 y = 0 ; y1 = x + 1
1 cos x
2. x 2 y¢¢ + xy¢ + ( x 2 - ) y = 0 ; y1 = ; x>0
4 x
3. ( x 2 + 1) y¢¢ + 2 xy¢ = 0 ; y1 = C1
4. (2 x + 1) y¢¢ - 4( x + 1) y¢ + 4 y = 0 ; y1 = e 2 x
5. y¢¢ - y = e 2 x ; y1 = e x
1
6. xy¢¢ + 2 y¢ = ; y1 = 1
x
7. y¢¢ - 4 y¢ + 5 y = 0 ; Φ1 = e 2 x sin x
8. xy¢¢ - 2 y¢ + 9 x 2 y = 0 ; Φ1 = sin( x3 )

Reduksi ke Bentuk Normal : u¢¢ - f ( x)u = 0

Teorema 7.11

Misalkan y adalah solusi dari:


y¢¢ + p ( x) y¢ + q ( x) y = 0
1

- p dx
dengan mendefinisikan y = uv ; v = e ,
maka didapat persamaan differensial
dalam bentuk normal:
1 1
u¢¢ - f ( x)u = 0 ; f ( x) = p¢ + p 2 - q
2 4

Bukti :

Diketahui y adalah solusi dari:

7 - 16
Bab 7 PD Order 2 Koef. Variabel

y¢¢ + p ( x) y¢ + q ( x) y = 0
misalkan : y = uv , maka y¢ = uv¢ + vu¢ dan y¢¢ = uv¢¢ + u¢v¢ + v¢u¢ + vu¢¢
subtitusikan ke PD, didapat :
(uv¢¢ + u¢v¢ + v¢u¢ + vu¢¢) + p (uv¢ + vu¢) + uv = 0
vu¢¢ + (2v¢ + pv)u¢ + (v¢¢ + pv¢ + qv)u = 0 (*)
dari persamaan (*) ini, u¢ akan hilang jika 2v¢ + pv = 0
1

- p dx
2v¢ + pv = 0 Û v = e
1
dari 2v¢ + pv = 0 , maka v¢ = - pv
2
1 1 1 1 1 1 1
v¢¢ = - p¢v - pv¢ = - p¢v - p (- pv) = - p¢v + p 2v
2 2 2 2 2 2 4
1 1 1 1
vu¢¢ + (2 × - pv + pv)u¢ + (- p¢v + p 2v + p × - pv + qv)u = 0
2 2 4 2
subtitusi ke (*), didapat :
1 1
u¢¢ - ( p¢ + p 2 - q )u = 0
2 4
u¢¢ - f ( x)u = 0
terbukti bahwa bentuk awal menjadi :
1 1
f ( x ) = p¢ + p 2 - q
2 4

Contoh 7.12:
1
1.Diketahui : x 2 y¢¢ + xy¢ + ( x 2 - ) y = 0
4
a. Reduksi PD di atas ke dalam bentuk normal
b. Tentukan 2 solusi yang bebas linear dari bentuk normal tersebut
c. Tentukan solusi umum y

Jawab :

1
a. x 2 y¢¢ + xy¢ + ( x 2 - ) y = 0
4
1 1
y¢¢ + y¢ + (1 - 2 ) y = 0
x 4x
1 1 1
p = dan q = 1 - 2 , maka p¢ = - 2
x 4x x
Menurut teorema 7.11, didapat :
1 1 1 1 1 1
2ò 2ò x
- p dx - dx - ln x - -
ln x 2
v=e =e =e 2 =e =x 2

dan :
1 1 1 1 1
f ( x) = p¢ + p 2 - q = - 2 + 2 - 1 + 2 = -1
2 4 2x 4x 4x

7 - 17
Bab 7 PD Oredr 2 Koef. Variabel

sehingga bentuk normalnya adalah :


u¢¢ + u = 0
( D 2 + 1)u = 0
u = C1 cos x + C2 sin x

b. Jadi dua solusi u yang bebas linear adalah u1 ( x) = cos x dan u2 ( x) = sin x
c. Karena solusi PD-nya adalah y = uv , maka dua solusi y yang bebas linear:
1 1
- -
y1 ( x) = x 2 cos x dan y2 ( x) = x 2 sin x

Solusi umumnya adalah :


1
-
y = uv = x 2 (C cos x + C sin x )
1 2

2. Diketahui : xy¢¢ - 2 xy¢ - 4 y = 0 . Tentukan bentuk normalnya

Jawab :
xy¢¢ - 2 xy¢ - 4 y = 0
4
y¢¢ - 2 y¢ - y = 0
x
4
p = -2 dan q = - , maka p¢ = 0
x

menurut teorema (7.11), didapat :


1

- p dx
= eò
dx
v=e = ex
dan :
1 1 1 1 4 4
f ( x) = p¢ + p 2 - q = × 0 + × 4 + = 1 +
2 4 2 4 x x

sehingga bentuk normalnya adalah :

u¢¢ - f ( x)u = 0
æ 4ö
u¢¢ - ç1 + ÷u = 0
è xø

Catatan :
Persamaan ini tidak dapat diintegralkan dari fungsi elementer, akan tetapi solusi
pendekatannya sangat sederhana untuk x ® ¥ . Solusi untuk kasus ini disebut solusi
asimtotis.

3. Untuk latihan, tentukan bentuk normal dari persamaan berikut:


a. y¢¢ - 4 y¢ + 3 y = 0
b. y¢¢ + 2 sin x y¢ + cos x (1 - cos x) y = 0

7 - 18
Bab 7 PD Order 2 Koef. Variabel

c. x 2 y¢¢ - 2 y¢ + 4 xy = 0

Fungsi Green

Pada bagian ini akan dibahas cara mancari solusi umum dari PD non homogen order 2
koefisien variabel :
a ( x) y¢¢ + b( x) y¢ + c( x) y = f ( x) .... (1)
Solusi umum ini memuat fungsi Green dari persamaan homogennya :
a ( x) y¢¢ + b( x) y¢ + c( x) y = 0 .... (2)
dengan a ( x), b( x), c( x) fungsi kontinu
Jika a ( x) ¹ 0 , maka fungsi Green dapat dibentuk dengan menggunakan 2 solusi dari (2)
yang bebas linear Φ1 dan Φ2 .
Berdasarkan teori integral tak tentu :
x
ò f ( x)dx = ò f (t )dt + C ;untuk x ³ a Î I
a
dan dari metode variasi parameter :
Φ2 f Φ f
y = C1Φ1 ( x) + C2Φ2 ( x) - Φ1 ( x) ò dx + Φ2 ( x) ò 1 dx
Wa Wa
maka didapat :
x
Φ2 f Φ (t ) f (t )
ò Wa dx = ò W2 (t )a(t ) + k1
x 0
x
Φ f Φ (t ) f (t )
ò Wa
1 dx =
ò W1 (t )a(t ) + k2
x0
Sehingga :
æ x Φ (t ) f (t ) ö æ x Φ (t ) f (t ) ö
ç ÷ ç ÷
y = C1Φ1 ( x) + C2Φ2 ( x) - Φ1 ( x)ç ò 2 dt + k1 ÷ + Φ2 ( x)ç ò 1 dt + k2 ÷
çx W (t ) a (t ) ÷ çx W (t ) a (t ) ÷
è 0 ø è 0 ø
x x
- Φ1 ( x)Φ2 (t ) f (t ) Φ ( x)Φ1 (t ) f (t )
= C1Φ1 ( x) + C2Φ2 ( x) - k1Φ1 ( x) + k2Φ2 ( x) + ò dt + ò 2 dt
x
W (t ) a (t ) x
W (t ) a (t )
0 0
x
(Φ1 ( x)Φ2 (t ) - Φ2 ( x)Φ1 (t ))
= (C1 - k1 )Φ1 ( x) + (C2 - k2 )Φ2 ( x) - ò W (t )a (t )
f (t )dt
x0
x
= C1Φ1 ( x) + C2Φ2 ( x) + ò G ( x, t ) f (t )dt
x0
dengan :

7 - 19
Bab 7 PD Oredr 2 Koef. Variabel

Φ1 ( x) Φ2 ( x)
1 Φ1 (t ) Φ2 (t )
G ( x, t ) = - ; x0 £ t £ x disebut fungsi Green.
a (t ) Φ1 (t ) Φ2 (t )
Φ1¢ (t ) Φ2¢ (t )
yc = C1Φ1 ( x) + C2Φ2 ( x) solusi komplementer
x
y p = ò G ( x, t ) f (t )dt solusi partikelir
x0
Catatan :
x x
¶G ( x, t )
Jelas G ( x, x) = 0 , maka y p = ò G ( x, t ) f (t )dt , y¢p = G ( x, x) f ( x) + ò f (t )dt
x0 x0
¶x
sehingga y¢p ( x0 ) = G ( x0 , x0 ) f ( x0 ) = 0
\Solusi partikelirnya y p memenuhi y p ( x0 ) = y¢p ( x0 ) = 0
Untuk memperoleh 2 solusi yang bebas linear Φ1 dan Φ2 bisa dari faktorisasi (koefisien
konstan) atau dengan bentuk normal ( koefisien Variabel).

Contoh 7.13:
Tentukan fungsi Green dari y¢¢ + y¢ = 0 , kemudian gunakan fungsi tersebut untuk
mencari solusi umum dari y¢¢ + y¢ = 3 sin x .

Jawab :
Dari bentuk a ( x) y¢¢ + b( x) y¢ + c( x) y = f ( x) , maka dari soal tersebut didapat :
a ( x) = 1 , b( x) = 0 , c( x) = 1 dan f ( x) = 3 sin x
Dan persamaan karakteristiknya adalah λ 2 + 1 = 0 dan akar karakteristiknya λ = ± i
Jadi dua solusi bebasnya adalah Φ1 = cos x dan Φ2 = sin x
Sehingga fungsi Green-nya adalah :
Φ1 ( x) Φ2 ( x) cos x sin x
1 Φ1 (t ) Φ2 (t ) cos t sin t sin x cos t - cos x sin t
G ( x, t ) = - = - =
a (t ) Φ1 (t ) Φ2 (t ) cos t sin t cos 2 t + sin 2 t
Φ1¢ ( x) Φ2¢ ( x) - sin t cos t
= sin( x - t ) ; x0 £ t £ x
Solusi umum :
x
y = C1Φ1 ( x) + C2Φ2 ( x) + ò G ( x, t ) f (t )dt
x0
x
= C1 cos x + C2 sin x + 3 ò sin( x - t ) sin t dt
x0
3 3 3
= C1 cos x + C2 sin x + sin x - sin(2 x0 - x) + ( x0 - x) cos x
4 4 2
3 3 3
= C1 cos x + C2 sin x + sin x + sin x - x cos x
4 4 2

7 - 20
Bab 7 PD Order 2 Koef. Variabel

3 3
= C1 cos x + C2 sin x + sin x - x cos x
2 2
Contoh soal 7.14:
Kerjakan soal berikut sebagaimana contoh 7.13 di atas:
1. y¢¢ + 2 y¢ - 8 y = 2e3 x 4. y¢¢ + y = cos x
2. y¢¢ - 9 y = e3 x 5. xy¢¢ - 3 y¢ = 4 x - 6
3. y¢¢ + 2 y¢ + y = e 2 x 6. x 2 y¢¢ - 2 xy¢ + 2 y = 2 x - 8

7 - 21

Anda mungkin juga menyukai