EAS PDP 2024
1. Diberikan persamaan diferensial parsial sebagai berikut:
∂u ∂ 2u
=k 2
∂t ∂x
dengan k adalah konstanta positif. Temukan solusi umum dari persamaan diferensial parsial
tersebut dengan menggunakan pemisalan eksponensial.
Asumsikan u(x, t) = eαx+βt adalah solusi untuk PD tersebut. Maka didapat
ut = βeαx+βt dan kuxx = kα2 eαx+βt
Tulis ulang persamaan sehingga menjadi
ut = kuxx
βeαx+βt = kα2 eαx+βt
n
β = kα2
re
Substitusi β = kα2 sehingga untuk semua α didapatkan solusi khusus
uα (x, t) = eαx+kα
2t
la
Karena persamaan linear, maka solusi umum adalah superposisi semua solusi yang mungkin
yaitu Z ∞
sC
2
u(x, t) = A(α)eαx+kα t dα
−∞
Ubah batasan domain menjadi [0, L] sehingga u(0, t) = T0 , u(L, t) = T1 , u(x, 0) = f (x)
didapat deret Fourier:
Ja
∞
T1 − T0 nπx −k( nπ
2
L )
t
X
u(x, t) = T0 + x+ Bn sin e
L n=1 L
dengan
2ZL T1 − T0 nπx
Bn = f (x) − T0 + x sin dx
L 0 L L
sehingga jika u(0, t) = 0, u(L, t) = 0, u(x, 0) = f (x) didapat deret Fourier:
∞
nπx −k( nπ
2
L )
t
X
u(x, t) = Bn sin e
n=1 L
dengan
2ZL nπx
Bn = f (x) sin dx
L 0 L
2. Diketahui PD: x2 uxx − 2xyuxy + y 2 uyy = 0
a) Ditinjau dari persamaan karateristik termasuk tipe apa?
1
b) Tentukan fungsi s(x, y) dan t(x, y) yang mengubah persamaan menjadi bentuk kanonik
dengan transformasi u(x, y) = w(s, t) dan dapatkan bentuk kanonik PD tersebut.
c) Dengan metode pemisahan variabel, selesaikan PD tersebut.
a) Bentuk umum PDP linier orde 2 dengan dua variabel adalah
Auxx + Buxy + Cuyy + Dux + Eut + F u = G
Berdasarkan bentuk umum didapat bahwa A = x2 , B = −2xy, C = y 2 , D = E = F =
G = 0, kemudian klasifikan tipe dengan cara
B 2 − 4AC
Substitusi A = x2 , B = −2xy, C = y 2 , D = E = F = G = 0 kedalam persamaan
(−2xy)2 − 4x2 y 2 = 0
Karena B 2 − 4AC = 0 maka tergolong PDP parabolik
n
b) Langkah-langkah penyelesaian:
i. Cari persamaan karakteristik
re
dy
dx
=
B
2A
la
dy −2xy
=
dx 2x2
sC
dy y
=−
dx x
Z
1 Z
1
dy = − dx
y x
Ja
ln y = − ln x + c1
ln y + ln x = c1
ln xy = ln C
xy = C
|{z}
ϕ1 (x, y)
sehingga diperoleh s(x, y) = xy, pilih t(x, y) = x maka
s s y x
J = t x t y = 1 0 = −x ̸= 0, karena x > 0
x y
ii. Cari turunan s, t terhadap x, y
sx =y tx = 1
sy =x ty = 0
sxx =0 txx = 0
sxy =1 txy = 0
syy =0 tyy = 0
2
iii. Substitusi (ii) ke persamaan
A∗ wss + B ∗ wst + C ∗ wtt + D∗ ws + E ∗ wt + F ∗ w = G∗
dengan
A∗ = As2x + Bsx sy + Cs2y
B∗ = 2Asx tx + B (sx ty + sy tx ) + 2Csy ty
C∗ = At2x + Btx ty + Ct2y
D∗ = Asxx + Bsxy + Csyy + Dsx + Esy
E∗ = Atxx + Btxy + Ctyy + Dtx + Ety
F∗ =F
G∗ =G
Didapat
A∗ = x2 y 2 − 2x2 y 2 + y 2 x2 = 0
B∗ = 2x2 y − 2x2 y + 0 = 0
n
C∗ = x2 − 0 + 0 = x2
D∗ = 0 − 2xy + 0 + 0 + 0 = −2xy
E∗
F∗
G∗
=0
=0
re
=0−0+0+0+0=0
la
Substitusi A∗ , B ∗ , C ∗ , D∗ , E ∗ , F ∗ , G∗ ke persamaan
A∗ wss + B ∗ wst + C ∗ wtt + D∗ ws + E ∗ wt + F ∗ w = G∗
sC
x2 wtt − 2xyws = 0
y
wtt = 2 ws
x
s
iv. Selesaikan hasil (iii) dengan substitusi t = x dan y = sehingga diperoleh
Ja
t
s
wtt = 2 2 ws
t
2s
Persamaan wtt = w
t2 s
memiliki solusi umum
Z
w(s, t) = 1
Aλ sλ/2 tr1 (λ) + Bλ sλ/2 tr2 (λ) dλ
λ>− 4
√
+ s−1/8 t (C + D ln t)
Z
λ/2
√
+ 1
s t (Eλ cos (ω ln t) + Fλ sin (ω ln t)) dλ
λ<− 4
+ G + H ln t
dengan
√ √
1+
1 + 4λ 1 + 1 + 4λ
r1 (λ) = , r1 (λ) =
q 2 2
|1 + 4λ|
ω=
2
3
Kembalikan ke variabel awal
Xh i
u(x, y) = Cλ (xy)λ/2 xr1 (λ) + Dλ (xy)λ/2 xr2 (λ)
λ
Atau dapat ditulis menjadi
u(x, y) = C1 + C2 x + C2 y + C3 x3 y + C4 xy 3 + . . .
dengan
u(x, y) = xa y b , (a − b)2 = a + b
c) Asumsikan u(x, y) = X(x)Y (t) adalah solusi untuk PD tersebut. Maka didapat
uxx = X ′′ (x)Y (y), uxy = X ′ (x)Y ′ (y), uyy = X(x)Y ′′ (y)
Tulis ulang persamaan sehingga menjadi
x2 uxx − 2xyuxy + y 2 uyy = 0
x2 X ′′ (x)Y (y) − 2xyX ′ (x)Y ′ (y) + y 2 X(x)Y ′′ (y) = 0
n
x2 X ′′ (x)Y (y) + y 2 X(x)Y ′′ (y) = 2xyX ′ (x)Y ′ (y)
Bagi kedua ruas dengan X(x)Y (Y )
x2
X ′′ (x)
X(x)
+ y2 re
Y ′′ (y)
Y (y)
= 2xy
X ′ (x) Y ′ (y)
X(x) Y (y)
la
Karena terdapat xy, maka Misal
sC
X ′ (x) k Y ′ (y) m
= , =
X(x) x Y (y) y
sehingga
X ′ (x)Y ′ (Y ) km
Ja
=
X(x)Y (y) xy
Substitusi ke persamaan sehingga
X ′′ (x) Y ′′ (y) km
x2 + y2 = 2xy
X(x) Y (y) xy
′′
X (x) km Y ′′ (y)
x2 = 2xy − y2 =λ
X(x) xy Y (y)
Solusi umum PD adalah
X
u(x, y) = Cλ Xλ (x)Yλ (y)
λ
Atau dapat ditulis menjadi
u(x, y) = C1 + C2 x + C2 y + C3 x3 y + C4 xy 3 + . . .
dengan
u(x, y) = xa y b , (a − b)2 = a + b
4
3. Selesaikan PDP konduksi panas
ut = uxx ; 0 ≤ x ≤ π, t > 0
dengan kondisi awal: u(x, 0) = x(π − x), 0 ≤ x ≤ π dan kondisi batas u(0, t) = 0 dan
u(π, t) = 0, t > 0
Andaikan u(x, t) = X(x)T (t) adalah solusi untuk PD tersebut. Maka didapat
uxx = X ′′ (x)T (t), ut = X(x)T ′ (t)
Tulis ulang persamaan sehingga menjadi
X(x)T ′ (t) = X ′′ (x)T (t)
Bagi kedua ruas dengan X(x)T (t) sehingga menjadi
X ′′ (x) T ′ (t)
= =λ
n
X(x) T (T )
X ′′ (x) T ′ (t)
= λ dan =λ
X(x)
re
Ambil λ negatif supaya didapat solusi non-trivial
T (t)
la
λ = −β 2 , β > 0
X ′′ (x)
= −β 2
X(x)
sC
X(x) = A cos (βx) + B sin (βx)
Masukan batasan soal
Ja
u(0, t) = 0 →
− AT (t) = 0
A=0
u(π, t) = 0 →
− (A cos (πβ) + B sin (πβ)) T (t) = 0
B sin (πβ) = 0
nπ
β= =n
π
Diperoleh
Xn (x) = Bn sin nx
Substitusi β = n ke T (t)
T ′ (t)
= −β 2
T (t)
2t
T (t) = Ce−β
2t
Tn (t) = Cn e−n
5
Solusi PD adalah
un (x, t) = Xn (x)Tn (t)
2t
= Bn sin (nx)Cn e−n
2
= an sin (nx)e−n t , n = 1, 2, 3, . . .
∞
2t
an sin (nx)e−n
X
u(x, t) =
n=1
Masukan kondisi
u(x, 0) = x (π − x)
∞
X
an sin (nx) = x (π − x)
n=1
Sehingga diperoleh
n
2Zπ
an = x (π − x) sin (nx) dx
π 0
Solusi umum PD adalah
u(x, t) =
∞
X
n=1
re 2t
an sin (nx)e−n
la
dengan
2Zπ
an = x (π − x) sin (nx) dx
sC
π 0
Ja