METODE RAYLEIGH-RITZ

STUDI LITERATUR
Disusun untuk memenuhi syarat kelulusan pada mata kuliah
studi literatur di Jurusan Matematika
Oleh
Riad Tauk Lazwardi
207700255
JURUSAN MATEMATIKA
FAKULTAS SAINS DAN TEKNOLOGI
UNIVERSITAS ISLAM GUNUNG DJATI
BANDUNG 2012
HALAMAN PENGESAHAN
METODE RAYLEIGH-RITZ
STUDI LITERATUR
Disusun untuk memenuhi syarat kelulusan pada mata kuliah
studi literatur di Jurusan Matematika
Oleh :
Riad Tauk Lazwardi
207700255
Telah Diperiksa dan Disetujui oleh
Pembimbing Pada Tanggal 18 Februari 2012
Dosen Pembimbing
Esih Sukaesih ,[Link].
197903132009012007
Mengetahui,
Ketua Jurusan Matematika
Elis Ratna Wulan ,[Link].,M.T.
197301122000032001
i
Abstraksi
Masalah nilai batas pada persamaan diferensial merupakan salah satu
masalah yang menarik untuk dikaji terutama pada permasalahan yang kom-
pleks sehingga dibutuhkan solusi [Link] berbagai metode untuk
[Link] satunya adalah metode [Link] literatur
ini menjelaskan bagaimana metode Rayleigh-Ritz menyelesaikan masalah nilai
batas pada persamaan
d
dx
_
p(x)
dy
dx
_
+ q(x)y = f(x) , dimana 0 x 1
dengan fungsi p(x) = 1 , q(x) =
2
, f(x) = 2
2
sin(x), panjang subin-
terval h
i
= h = 0.1 sehingga x
i
= 0, 1ii = 1, 2, . . . , 9 dan fungsi basis berupa
PWL(Piecewise Linier) yang berbentuk:

i
(x) =
_

_
0, if 0 x x
i1
;
(xx
i1
)
h
i1
, if x
i1
< x x
i
;
(x
i+1
x)
h
i
, if x
i
< x x
i+1
;
0, if x
i+1
< x 1;
Metode Rayleigh Ritz dengan penggunaan fungsi basis di atas akan
menghasilkan sistem persamaan linier yang jika dibentuk matriks Ac = b a-
kan menghasilkan matriks A tridiagonal simetri yang denit positif sehingga
nonsingular (solusi spl unik) dan diperoleh komputasi (penyelesaian spl) yang
stabil, solusi persamaan balok di atas unik, stabil dan mempunyai orde error
O(h
2
).
Penggunaan PWL sebagai fungsi basis dalam pendekatan solusi ke persa-
maan di atas mempunyai kelemahan yaitu kontinyu tetapi tidak diferensiabel
di [0, 1].
ii
KATA PENGANTAR
Assalamualaikum warahmatullahiwabarakatuh.
Terima kasih kepada Allah SWT karena dengan rahman dan rahim-
Nya saya dapat menyelesaikan studi literatur [Link] lupa kepada Ibu Esih
Sukaesih, [Link]. selaku pembimbing yang telah dengan sabar membimbing saya
sehingga mata kuliah studi literatur ini dapat [Link] saya
ucapkan terima kasih kepada orang tua yang telah mendukung dan mendoakan
tanpa henti.
Saya akui studi literatur ini jauh dari sempurna, tetapi semoga dapat
bermanfaat.
Wassalamualaikumwarahmatullahiwabarakatuh.
iii
Daftar Isi
BAB I PENDAHULUAN 1
1.1 Latar Belakang Masalah . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Rumusan Masalah . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.3 Batasan Masalah . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.4 Tujuan Penelitian . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.5 Metode Penelitian . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.6 Sistematika Penulisan . . . . . . . . . . . . . . . . . . . . . . . . 2
1.7 Kerangka Berkir . . . . . . . . . . . . . . . . . . . . . . . . . . 3
BAB II LANDASAN TEORI 4
2.2 Persamaan Diferensial Biasa . . . . . . . . . . . . . . . . . . . . 4
2.3 Piecewise Linier . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.4 Masalah Nilai Batas . . . . . . . . . . . . . . . . . . . . . . . . 5
2.4.1 Persamaan Homogen . . . . . . . . . . . . . . . . . . . . 5
2.4.2 Persamaan nonhomogen dengan nilai batas homogen . . 6
2.5 Variational Principle . . . . . . . . . . . . . . . . . . . . . . . . 6
BAB III METODE RAYLEIGH-RITZ 9
3.2 Langkah-Langkah Metode Rayleigh-Ritz . . . . . . . . . . . . . 11
3.3 Penerapan Metode Rayleigh-Ritz pada Persamaan Balok . . . . 13
BAB IV KESIMPULAN dan SARAN 21
4.2 Kesimpulan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
4.3 Saran . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
iv
Daftar Gambar
1 Gambar Tekanan Pada Balok dengan Penampang Persegi Panjang 9
2 Gambar Fungsi Basis . . . . . . . . . . . . . . . . . . . . . . . . 11
3 Gambar y = sin(x) . . . . . . . . . . . . . . . . . . . . . . . . 14
4 Gambar y = (x) . . . . . . . . . . . . . . . . . . . . . . . . . . 16
5 Gambar y = sin(x) dan y = (x) . . . . . . . . . . . . . . . . 20
Daftar Tabel
1 Tabel y = sin(x) dan y = (x) . . . . . . . . . . . . . . . . . . 20
v
BAB I PENDAHULUAN
1.1 Latar Belakang Masalah
Topik persamaan diferensial dalam matematika merupakan objek kajian
yang sangat [Link] dengan perkembangan mesin hitung pada jaman
dulu, para matematikawan mulai mencari berbagai metode untuk memecahkan
masalah-masalah pada topik ini dengan cepat dan akurat.
Salah satu masalah dalam persamaan diferensial adalah masalah nilai
[Link] nilai batas adalah masalah tersendiri dalam memperoleh solusi
pada persamaan [Link] sekali metode untuk menyelesaikan masa-
lah tersebut sehingga diperlukan pengetahuan untuk mengetahui kelemahan
dan kelebihan suatu metode yang digunakan.
Tujuan SL ini adalah untuk memahami salah satu metode untuk menye-
lesaikan masalah nilai batas yaitu metode Rayleigh-Ritz dan penerapannya.
1.2 Rumusan Masalah
Pada studi literatur ini akan dibahas persamaan diferensial linier nonho-
mogen dengan nilai batas homogen dan akan ditelaah:
1. Bagaimana langkah-langkah metode Rayleigh-Ritz?
2. Apakah penggunaan metode Rayleigh-Ritz menghasilkan solusi unik dan
stabil ?
1.3 Batasan Masalah
Studi literatur ini hanya membahas metode Rayleigh-Ritz dan contoh
penerapannya dalam persamaan tekanan pada balok,yaitu :

d
dx
_
p(x)
dy
dx
_
+ q(x)y = f(x) , dimana 0 x 1
dengan fungsi p(x) = 1 , q(x) =
2
, f(x) = 2
2
sin(x), panjang su-
binterval h
i
= h = 0.1 sehingga x
i
= 0, 1ii = 1, 2, . . . , 9 dan fungsi basisnya
:
1

i
(x) =
_

_
0, if 0 x x
i1
;
(xx
i1
)
h
i1
, if x
i1
< x x
i
;
(x
i+1
x)
h
i
, if x
i
< x x
i+1
;
0, if x
i+1
< x 1;
1.4 Tujuan Penelitian
Tujuan penelitian dalam studi literatur ini adalah :
Mengetahui langkah-langkah metode Rayleigh-Ritz dan penerapannya
pada persamaan balok (21).
Mengetahui apakah penggunaan metode Rayleigh-Ritz dalam persamaan
balok (21) dengan PWL (16), partisi menjadi 9 subinterval dan panjang
subintervalnya 0, 1 menghasilkan solusi unik dan stabil?
1.5 Metode Penelitian
Metode Penelitian studi literatur ini adalah tinjauan pustaka.
1.6 Sistematika Penulisan
1. Bab I : Latar Belakang Masalah, Rumusan Masalah, Batasan Ma-
salah, Tujuan Penelitian, Metode Penelitian, Sistematika Penulisan, Ke-
rangka Berkir
2. Bab II : Landasan Teori
3. Bab III : Metode Rayleigh-Ritz
4. Bab IV : Kesimpulan dan Saran
2
1.7 Kerangka Berkir
Dasar pemilihan topik masalah nilai batas adalah begitu banyaknya me-
tode yang digunakan dalam penyelesaiannya sehingga harus ditelaah kelema-
han dan kelebihannya .Langkah-langkah penyelesaian masalah nilai batas de-
ngan metode Rayleigh-Ritz adalah :
1. Partisi interval menjadi beberapa subinterval.
2. Tentukan fungsi basis.
3. Bentuk sistem persamaan linier.
4. Cari konstanta c
i
untuk
n
i=1
c
i

i
.
3
BAB II LANDASAN TEORI
2.2 Persamaan Diferensial Biasa
Persamaan diferensial biasa adalah persamaan yang memuat x, y(x) be-
serta turunan-turunan dari y(x).Bentuk umumnya adalah :
F(x, y, y

, . . . , y
n
) = 0
Persamaan diferensial biasa dikatakan linier jika fungsi F linier terha-
dap y, y

, . . . , y
(n)
, namun tidak perlu linier terhadap variabel [Link] umum
persamaan diferensial biasa linier berbentuk:
a
0
(x)y
(n)
+ a
1
(x)y
(n1)
+ . . . + a
n
(x)y = g(x) (1)
[10].
Jika y = f(x) memenuhi persamaan diferensial, maka f(x) dikatakan so-
lusi dari persamaan diferensial [Link] umum suatu persamaan diferen-
sial adalah bentuk umum solusi persamaan diferensial [Link] umum
bisa menjadi solusi khusus dengan adanya informasi/syarat tambahan , disebut
syarat awal/syarat batas.
Dua buah fungsi dikatakan bebas linier jika tidak saling berkelipatan.
Teorema untuk mengetahui fungsi yang bebas linier yaitu:
Teorema Kebebaslinieran dua buah fungsi 2.2.1 Wronskian dari dua
fungsi yang bergantung linier /konjugat pada selang I bernilai nol.
W(f, g)

f g
f

= 0 , x I
2.3 Piecewise Linier
Denisi [Link] I = [a, b] adalah [Link] fungsi g : I R
dikatakan pieciwise linier pada I jika I adalah union pada sejumlah interval
I
1
, . . . , I
m
dimana fungsi g pada setiap I
k
adalah fungsi linier [6].
4
Teorema Piecewiselinier 2.3.1 Misal I interval terbatas yang tertutup dan
misal f : I R kontinyu di [Link] > 0 maka terdapat fungsi piecewise linier
yang kontinyu g

: I R dimana |f(x) g

(x)| < untuk semua x I.

2.4 Masalah Nilai Batas
Masalah nilai batas pada persamaan diferensial biasa dibagi 3 menurut
nilai batasnya [4] yaitu :
1. persamaan diferensial homogen
2. persamaan diferensial nonhomogen dengan nilai batas homogen
3. persamaan diferensial nonhomogen dengan nilai batas nonhomogen
2.4.1 Persamaan Homogen
Contoh masalah nilai batas sederhana pada persamaan homogen yaitu :
y

+ y = 0
dengan kondisi batas y(0) = 1 dan y(

2
) = 2 mempunyai solusi umum:
y(x) = c
1
cos(x) + c
2
sin(x)
untuk memenuhi kondisi y(0) = 1 dan y(

2
) = 2 ,kita tinggal memasukan
syarat ke solusi umum sehingga diperoleh solusi khusus y(x) = cosx + 2sinx.
Teorema Persamaan Diferensial Homogen 2.4.1 Persamaan diferensial
linier homogen y

+p(x)y

+q(x)y = 0 mempunyai koesien kontinyu p(x) dan

q(x) pada interval a x b dengan nilai batas y(a) = dan y(b) = agar
mempunyai solusi yang unik adalah:
1. solusi yang hanya memenuhi y(a) = y(b) = 0 adalah solusi nol
2. untuk sepasang solusi bebas linier misal y
1
, y
2
determeniannya yaitu
y
1
(a)y
2
(b) y
2
(a)y
1
(b) = 0
Jika tidak bebas linier maka disebut konjugat dan tidak mempunyai so-
lusi yang unik [4].
5
2.4.2 Persamaan nonhomogen dengan nilai batas homogen
Persamaan nonhomogen dengan nilai batas homogen dapat diselesaikan
dengan formula variasi parameter yaitu y = y
h
(x) + y
p
(x) dimana y
h
adalah
solusi homogennya dan solusi partikularnya:
y
p
(x) = y
1
(x)
_
x
b
y
2
(t)f(t)
w(t)
dt + y
2
(x)
_
x
a
y
1
(t)f(t)
w(t)
dt
y
p
(x) =
_
x
a
y
2
(x)y
1
(t)f(t)
w(t)
dt +
_
b
x
y
1
(x)y
2
(t)f(t)
w(t)
dt
dimana w(t) = y
1
(t)y

2
(t) y
2
(t)y

1
(t) adalah wronksian dari y
1
, y
2
.Selan-
jutnya adalah mensubstitusikan syarat nilai batas sehingga diperoleh solusi
khusus [4].
Contoh sederhana dari penyelesaian masalah nilai batas nonhomogen de-
ngan menggunakan metode variasi parameter yaitu menyelesaikan persamaan:
y

(x) = 1, y(0) = y(1) = 0

Solusi homogennya:
y
h
(x) = c
1
+ c
2
(x)
Substitusikan syarat nilai [Link] memenuhi y
1
(0) = 0, harus mempu-
nyai c
1
= 0 maka dapat memilih dengan mudah c
2
= [Link] x = 1 dan harus
mempunyai c
1
+ c
2
= 0 , sehingga dapat memilih c
1
= 1, c
2
= [Link] :
y
1
(x) = x, y
2
(x) = 1 x
untuk solusi homogen bebas [Link] adalah
W(f, g)

x 1 x
1 1

= 1
y
p
(x) =
x
2
2

x
2
. y = x + 1 x +
x
2
2

x
2
.
2.5 Variational Principle
Pada variational principle (Kalkulus variasi) fokus utamanya yaitu men-
cari fungsi yang memaksimasi atau meminimasi nilai suatu [Link] pa-
6
ling sederhana adalah mencari y(x) pada interval a x b yang memaksimasi
atau meminimasi nilai integral tentu :
I[y] =
_
b
a
F(x, y(x), y

(x))dx , dimana ()

=
d()

dx
(2)
min
yS
I[y] | y(a) = A, y(b) = B
dimana S adalah himpunan semua fungsi yang memenuhi syarat batas.
Pada umumnya digunakan pendekatan numerik terhadap solusi karena
solusi eksak untuk masalah variational principle hanya ada untuk masalah
yang simpel [2].
Contoh :
I[y] =
_
1
0
1
2
[(y

)
2
y
2
2sin(x)y]dx , dimana y(0) = y(1) = 0 (3)
Langkah-langkah dalam masalah variational principle:
1. Tentukan partisi
2. Tentukan fungsi basis
3. Cari turunannya terhadap koesien sehingga diperoleh fungsi yang me-
minimasi nilai I[y].
Misal interval I dipartisi menjadi n subinterval dengan fungsi basisnya
i
(x) =
x
i
(1 x).Pada subinterval ke-1.
y
1
= a
1

1
(x)
dimana a
1
tidak [Link] diketahui bahwa
1
(0) = (1) = 0 sehing-
ga y(x) memenuhi syarat nilai [Link] menemukan y
1
, dinotasikan
y
1
,yang menyebabkan I[y] minimum ini sama saja dengan menemukan kon-
stanta optimal a
1
.
I[y
1
] =
_
1
0
1
2
[(a
1

1
)
2
(a
1

1
)
2
2sin(x)a
1

1
]dx =
1
2
a
2
1

11
a
1

10
(4)
dimana
11
=
_
1
0
[(

1
)
2
(
1
)
2
]dx =
3
10
,
10
=
_
1
0
sin(x)
1
dx =
4

3
I[y
1
] =
1
2
a
2
1

11
a
1

10
=
1
2
a
2
1
3
10
a
1
4

3
(5)
7
Selanjutnya mencari konstanta yang tepat agar nilai integral minimal yaitu
dengan menganggap I sebagai fungsi dari a
1
.Jelas, supaya a
1
minimal turunan
pertamaya harus nol.
I

a
1
= a
1

11

10
= 0. (6)
I

a
1
=
3
10
(7)
I(y
1
) =
1
2
a
2
1

11
a
1

10
=
1
2
a
2
1
3
10
a
1
4

3
dengan a
1
=
40
3
3
adalah minimum karena turunan keduanya positif.
Tapi I(y
1
) bisa jadi bukan minimum yang sebenarnya (minimum lokal)
dimana ini diperoleh dari subinterval ke-1 .Maka pilih
2
= x
2
(1 x) dan me-
lakukan cara yang [Link] dibandingkan sehingga diperoleh fungsi
y yang minimum.
8
BAB III METODE RAYLEIGH-RITZ
Gambar 1: Gambar Tekanan Pada Balok dengan Penampang Persegi Panjang
Salah satu masalah persamaan diferensial linier dengan 2 nilai batas pada
analisis tekanan-balok adalah pada persamaan :

d
dx
_
p(x)
dy
dx
_
+ q(x)y = f(x) , dimana 0 x 1 (8)
dengan kondisi batas
y(0) = y(1) = 0 (9)
Persamaan dierensial ini menggambarkan deeksi y(x) pada balok de-
ngan panjang 1, luas penampang-lintang q(x) dan tekanan f(x) dan p(x).
Teorema Penyelesaian Persamaan Balok 3.1 Misal p C
1
[0, 1] dan
p(x) > 0, q(x) 0 , dimana 0 x 1
fungsi y C
2
0
[0, 1] adalah solusi unik pada persamaan diferensial

d
dx
_
p(x)
dy
dx
_
+ q(x)y = f(x) , dimana 0 x 1
jika dan hanya jika y adalah fungsi unik di C
2
0
[0, 1] yang meminimasi integral
I[u] =
_
1
0
p(x)[u

(x)]
2
+ q(x)[u(x)]
2
2f(x)u(x)dx (10)
Bukti teorema ini dapat dilihat di [7].
Variational principle (Kalkulus Variasi) untuk persamaan (8) adalah ide
dasar dari metode Rayleigh-Ritz sehingga solusinya adalah fungsi yang memi-
nimasi nilai integral tentu I[u].
9
Metode Rayleigh-Ritz mendekati solusi fungsi y dengan fungsi
(x) =
n

i=1
c
i

i
(x)
.
Karena termasuk dalam variational principle, langkah-langkah dalam
me-tode Rayleigh-Ritz nilai integral I diminimasi tidak pada semua fungsi
di C
2
0
[0, 1] tetapi pada himpunan fungsi yang lebih kecil yang terdiri dari kom-
binasi linier fungsi basis
1
,
2
, . . . ,
n
.Fungsi basis {
i
}
n
i=1
harus bebas linier
dan memenuhi:

i
(0) =
i
(1) = 0 , dimana i = 1, 2, . . . , n
Pendekatan solusi y(x) dengan menggunakan fungsi kombinasi linier (x) =

n
i=1
c
i

i
(x) diperoleh dengan cara mencari konstanta c
1
, c
2
, . . . , c
n
untuk me-
minimasi I[
n
i=1
c
i

i
(x)].
Dari persamaan (10) misalkan u = .
I[] = I
_
n

i=1
c
i

i
_
=
_
1
0
_
p(x)
_
n

i=1
c
i

i
(x)
_
2
+ q(x)
_
n

i=1
c
i

i
(x)
_
2
2f(x)
n

i=1
c
i

i
(x)
_
dx (11)
Untuk meminimumkan nilai integral I[u], misalkan I adalah fungsi yang
variabel bebasnya c
1
, c
2
, . . . , c
n
kemudian turunkan secara parsial yaitu:
I
c
j
= 0 untuk setiap j = 1, 2, . . . , n (12)
Mendiferensialkan (11) diperoleh :
I
c
j
=
_
1
0
_
2p(x)
n

i=1
c
i

i
(x)

j
(x) + 2q(x)
n

i=1
c
i

i
(x)
j
(x)
2f(x)
j
(x)
_
dx (13)
10
dan mensubstitusikan ke persamaan (12) diperoleh :
0 =
n

i=1
__
1
0
p(x)

i
(x)

j
(x) + q(x)
i
(x)
j
(x)dx
_
c
i

_
1
0
f(x)
j
(x)dx untuk setiap j = 1, 2, . . . , n (14)
sehingga akan diperoleh sistem linier Ac = b dengan A merupakan ma-
triks yang entrinya :
a
ij
=
_
1
0
[p(x)

i
(x)

j
(x) + q(x)
i
(x)
j
(x)]dx (15)
dan matriks b yang entrinya : b
i
=
_
1
0
f(x)
i
(x)dx.
3.2 Langkah-Langkah Metode Rayleigh-Ritz
Langkah pertama dalam metode Rayleigh-Ritz adalah membentuk par-
tisi pada [0, 1] dengan memilih titik x
0
, x
1
, . . . , x
n+1
dengan 0 = x
0
<
x
1
< . . . , x
n
, x
n+1
= 1 .Misal h
i
= x
i+1
x
i
untuk setiap i = 0, 1, . . . , n
adalah lebar setiap subinterval.
Langkah kedua adalah menentukan fungsi basis, dimana pada SL ini
fungsi basisnya:
Gambar 2: Gambar Fungsi Basis
11

i
(x) =
_

_
0, if 0 x x
i1
;
(xx
i1
)
h
i1
, if x
i1
< x x
i
;
(x
i+1
x)
h
i
, if x
i
< x x
i+1
;
0, if x
i+1
< x 1;
(16)
Maka turunannya diperoleh :

i
(x) =
_

_
0, if 0 < x < x
i1
;
1
h
i1
, if x
i1
< x < x
i
;
1
h
i
, if x
i
< x < x
i+1
;
0, if x
i+1
< x < 1;
untuk setiap i = 1, 2, . . . , n. Karena
i
dan

i
tidak nol hanya pada
(x
i1
, x
i+1
),maka :

i
(x)
j
(x) = 0
dan

i
(x)

j
(x) = 0
kecuali ketika j = i 1 atau i + 1 sitem persamaan linier (14) apabila
dirubah ke dalam matriks n x n akan diperoleh matriks [Link]
bukan nol di A diperoleh dari :
a
ii
=
_
1
0
[p(x)[

i
(x)]
2
+ q(x)[
i
(x)]
2
]dx
=
_
x
i
x
i1
(
1
h
i1
)
2
p(x)dx +
_
x
i+1
x
i
(
1
h
i
)
2
p(x)dx
+
_
x
i
x
i1
(
1
h
i1
)
2
(x x
i1
)
2
q(x)dx
+
_
x
i+1
x
i
(
1
h
i
)
2
(x
i+1
x)
2
q(x)dx (17)
untuk setiap i = 1, 2, . . . , n.
a
i,i+1
=
_
1
0
p(x)

i
(x)

i+1
(x) + q(x)
i
(x)
i+1
(x)dx
=
_
x
i+1
x
i
(
1
h
i
)
2
p(x)dx
+
_
x
i+1
x
i
(
1
h
i
)
2
(x
i+1
x)(x x
i
)q(x)dx (18)
12
untuk setiap i = 1, 2, 3, . . . , n 1.
a
i,i1
=
_
1
0
p(x)

i
(x)

i1
(x) + q(x)
i
(x)
i1
(x)dx
=
_
x
i
x
i1
(
1
h
i1
)
2
p(x)dx
+
_
x
i
x
i1
(
1
h
i1
)
2
(x
i
x)(x x
i1
)q(x)dx (19)
untuk setiap i = 2, 3, . . . , n. Nilai entri b diperoleh dari:
b
i
=
_
1
0
f(x)
i
(x)dx
=
_
x
i
x
i1
1
h
i1
(x x
i1
)f(x)dx +
_
x
i+1
x
i
1
h
i
(x
i+1
x)f(x)dx (20)
untuk setiap i = 1, 2, . . . , n.
Langkah ketiga adalah memasukan entri-entri di atas ke dalam persama-
an Ac = b.
Langkah keempat adalah menyelesaikan Ac = b sehingga diperoleh kon-
stanta c
1
, c
2
, . . . , c
n
.Selanjutnya masukan ke dalam persamaan:
y(x) =
n

i=1
c
i

i
(x)
3.3 Penerapan Metode Rayleigh-Ritz pada Persamaan
Balok
y

+
2
y = 2
2
sin(x), 0 x 1, y(0) = y(1) = 0 (21)
Penyelesaian dengan Metode Variasi Parameter.
y

2
y = 2
2
sin(x)
y
h
(x) = c
1
e
x
+ c
2
e
x
13
y
p
(x) = uy
1
+ vy
2
y
p
(x) =
_
x
a
y
2
(x)y
1
(t)f(t)
w(t)
dt +
_
b
x
y
1
(x)y
2
(t)f(t)
w(t)
dt
y
p
(x) = sin(x)
y(x) = sin(x) + c
1
e
x
+ c
2
e
x
masukan kondisi batas, diperoleh:
y(x) = sin(x)
Gambar 3: Gambar y = sin(x)
Penyelesaian dengan Metode Rayleigh-Ritz.
Langkah [Link]. Misal h
i
= h = 0.1 sehingga x
i
= 0.1i untuk
setiap i = 0, 1, 2, . . . , 9
Langkah [Link] fungsi [Link] basisnya adalah (16).
Langkah [Link] sistem persamaan [Link] h
i
= h = 0.1,q(x)
=
2
dan p(x) = 1 maka :
a
ii
= (
1
h
i1
)
2
_
x
i
x
i1
dx + (
1
h
i
)
2
_
x
i+1
x
i
dx
+ (
1
h
i1
)
2

2
_
x
i
x
i1
(x x
i1
)
2
dx
+ (
1
h
i
)
2

2
_
x
i+1
x
i
(x
i+1
x)
2
dx
= 20 +

2
15
(22)
14
a
i,i+1
= (
1
h
i
)
2
_
x
i+1
x
i
dx
+ (
1
h
i
)
2

2
_
x
i+1
x
i
(x
i+1
x)(x x
i
)dx
= 10x|
x
i+1
x
i
+ 100
2
_
x
i+1
x
i
(x
i+1
+ x
i
) x
i+1
x
i
x
2
(23)
untuk i = 1, . . . , 8.
a
i,i1
= (
1
h
i1
)
2
_
x
i
x
i1
dx
+(
1
h
i1
)
2

2
_
x
i
x
i1
(x
i
x)(x x
i1
)dx
= 10x|
x
i
x
i1
+ 100
2
_
x
i
x
i
1
(x
i
+ x
i1
) x
i1
x
i
x
2
(24)
untuk i = 2, . . . , [Link] diperoleh matriks tridiagonal simetri
dengan entri :
a
i,i+1
untuk i = 1, . . . , 8 sama dengan a
i,i1
untuk i = 2, . . . , 9 yaitu
10 +

2
60
dan matriks b mempunyai entri :
b
1
= 3.297.10
6
b
5
= 1.6483.10
5
b
2
= 6.594.10
6
b
6
= 1.9778.10
5
b
3
= 9.89.10
6
b
7
= 2.3072.10
5
b
4
= 1.3187.10
5
b
8
= 2.6367.10
5
b
9
= 1.096
Langkah [Link] konstanta c
i
.
Teorema Eliminasi Gauss 3.3.1 Jika A adalah matriks n x n
yang denit positif, maka A [Link] Eliminasi Ga-
uss dapat digunakan pada sembarang sistem linier Ax = B untuk
memperoleh solusi unik tanpa perubahan baris atau kolom, dan kom-
putasi(perhitungan spl) stabil .
Denisi 3.2.1 .Denit positif adalah jika x
t
Ax > 0n dimensi.
Karena matriks A denit positif maka dengan metode Eliminasi
Gauss diperoleh [7] :
15
c
1
= 0.3102866742 c
6
= 0.9549641893
c
2
= 0, 5902003271 c
7
= 0.8123410598
c
3
= 0.8123410598 c
8
= 0.5902003271
c
4
= 0.9549641893 c
9
= 0.3102866742
c
5
= 1.004108771
sehingga
(x) =
9
i=1
c
i

i
(x)
Gambar 4: Gambar y = (x)
Penyelesaian dengan Algoritma Rayleigh-Ritz
input integer n sehingga x
0
= 0 < x
1
< . . . < x
n
< x
n+1
= 1.
output koesien c
1
, . . . , c
n
.
For i = 0, . . . , n , h
i
= x
i+1
x
i
.
For i = 0, . . . , n, denisikan fungsi basis (16).
Buat variabel untuk nilai inisial.
IOLD1 = (
1
h
0
)
2
_
x
1
x
0
p(x)dx
= 100
_
0.1
0
dx
= 10 (25)
16
IOLD2 = (
1
h
0
)
2
_
x
1
x
0
(x x
0
)
2
q(x)dx
= 100
_
0.1
0
x
2

2
dx
=

2
30
(26)
For i = 1, . . . , n [Link] variabel :
INEW1 = (
1
h
i
)
2
_
x
i+1
x
i
p(x)dx
= 100
_
0.1i+0.1
0.1i
dx
= 10 (27)
INEW2 = (
1
h
i
)
2
_
x
i+1
x
i
(x x
i
)
2
q(x)dx
= 100
_
0.1i+0.1
0.1i
(x 0.1i)
2

2
dx
=

2
30
(28)
I3 = (
1
h
i
)
2
_
x
i+1
x
i
(x
i+1
x)
2
q(x)dx
= 100
_
0.1i+0.1
0.1i
(0.1i + 0.1 x)
2

2
dx
=

2
30
(29)
I4 = (
1
h
i
)
2
_
x
i+1
x
i
(x
i+1
x)(x x
i
)q(x)dx
= 100
_
0.1i+0.1
0.1i
(0.1i + 0.1 x)(x 0.1i)
2
dx
=

2
60
(30)
17
I5 =
1
h
i1
_
x
i
x
i1
(x x
i1
)f(x)dx
= 10
_
0.1i
0.1i0.1
(x 0.1i + 0.1)2
2
sin(x)dx
= 2 cos(0.1i) + 20[sin(0.1i) sin((0.1i 0.1))](31)
I6 =
1
h
i
_
x
i+1
x
i
(x
i+1
x)f(x)dx
= 10
_
0.1i+0.1
0.1i
(0.1i + 0.1)
2

2
sin(x)dx
= 2 cos(0.1i) 20[sin((0.1i + 0.1)) sin(0.1i)] (32)
a
ii
= IOLD1 + INEW1 + IOLD2 + I3 = 20 +

2
15
a
i,i+1
= INEW1 + I4 = 10 +

2
60
b
i
= I5 + I6 = 40 sin(0.1i)(1 cos(0.1))
IOLD1 + INEW1 = 10
IOLD2 + INEW2 =

2
30
Buat variabel lagi untuk menentukan a
99
dan b
9
.
INEW1 = (
1
h
n
)
2
_
x
n+1
x
n
p(x)dx
= 100
_
1
0.9
dx
= 10 (33)
I3 = (
1
h
n
)
2
_
x
n+1
x
n
(x
n+1
x)
2
q(x)dx
= 100
_
1
0.9
(1 x)
2

2
dx
=

2
30
(34)
18
I5 =
1
h
n1
_
x
n
x
n1
(x x
n1
)f(x)dx
= 2 cos(0.9) + 20[sin(0.9) sin(0.8)] (35)
I6 =
1
h
n
_
x
n+1
x
n
(x
n+1
x)f(x)dx
= 2 cos(0.9) + 20 sin(0.9) (36)
a
n,n
= IOLD1 + INEW1 + IOLD2 + I3 =
9
= 20 +

2
15
b
n
= I5 + I6 = b
9
= 40 sin(0.9) 20 sin(0.8)
Selesaikan sistem persamaan [Link]:
c
1
= 0.3102866742 c
6
= 0.9549641893
c
2
= 0, 5902003271 c
7
= 0.8123410598
c
3
= 0.8123410598 c
8
= 0.5902003271
c
4
= 0.9549641893 c
9
= 0.3102866742
c
5
= 1.004108771
sehingga
(x) =
9
i=1
c
i

i
(x)
sedangkan dalam solusi analitik diperoleh :
y(x) = sin x
Error dari setiap aproksimasi di x
i
untuk setiap i = 1, 2, . . . , 9 dapat
dilihat dari Tabel 1 sehingga diperoleh solusi stabil dan :
|(x) y(x)| = O(h
2
) 0 x 1
Karena O(h
2
) maka apabila subinterval dibagi 2, errornya adalah error sebe-
lumnya dibagi 4 [9].Pembuktian ini terdapat pada buku Schultz [7].
19
i x
i
(x
i
) y(x
i
) |(x
i
) y(x
i
)|
1 0.1 0.3102866742 0.3090169943 0.00127
2 0.2 0.5902003271 0.5877852522 0.00242
3 0.3 0.8123410598 0.8090169943 0.00332
4 0.4 0.9549641896 0.9510565162 0.00391
5 0.5 1.004108771 1.0000000000 0.00411
6 0.6 0.9549641893 0.9510565162 0.00391
7 0.7 0.8123410598 0.8090169943 0.00332
8 0.8 0.5902003271 0.5877852522 0.00242
9 0.9 0.3102866742 0.3090169943 0.00127
Tabel 1: Tabel y = sin(x) dan y = (x)
Gambar 5: Gambar y = sin(x) dan y = (x)
20
BAB IV KESIMPULAN dan SARAN
4.2 Kesimpulan
Langkah-langkah Metode Rayleigh-Ritz :
1. Partisi interval menjadi beberapa subinterval.
2. Tentukan fungsi basis.
3. Bentuk sistem persamaan linier.
4. Cari konstanta c
i
untuk
n
i=1
c
i

i
.
Penerapan langkah-langkah Metode Rayleigh-Ritz pada persamaan (21) :
1. Partisi interval menjadi 9 subinterval dengan panjang 0, 1.
2. Fungsi basis yang dipilih yaitu (16).
3. Bentuk sistem persamaan linier Ac = b dari persamaan (10).
4. Gunakan Eliminasi Gauss untuk mencari konstanta c
i
,setelah diketaui
masukkan pada persamaan
n
i=1
c
i

i
.
Penggunaan PWL (16) sebagai fungsi basis dalam pendekatan solusi per-
samaan (21) dengan partisi menjadi 9 subinterval dan panjang subintervalnya
0, 1 menghasilkan :
Sistem persamaan linier yang jika dibentuk ke dalam matriks Ac = b
maka akan menghasilkan matriks A yang berbentuk tridiagonal simetri
dan denit positif sehingga komputasi (penyelesaian spl) stabil.
Solusi persamaan (21) unik dan stabil.
4.3 Saran
Penggunaan PWL pada studi literatur ini mempunyai kelemahan yaitu
kontinyu tetapi tidak diferensiabel di [0, 1] sehingga diperlukan fungsi basis
yang lebih kompleks misalnya
i
C
2
0
[0, 1] yang diharapkan memperoleh solusi
aproksimasi yang lebih baik.
21
Daftar Pustaka
[1] Francis Scheild,Ph.D, 2000 Solved Problems In Numerical Analysis,
Boston University, McGraw-Hill Publishing Company.
[2] Frederick [Link] , Introduction to the calculus of variations and its
application, Chapman Hall Mathematics.
[3] O.C Zienkiewicz and [Link] , Finite Elements And Approxima-
tion, University of Wales, Swansea, United Kingdom, John Wiley and
Sons.
[4] Richard .[Link] , Introduction Dierental equation and Dyna-
mic System.
[5] Richard L. Burden and [Link] Faires , Numerical Analysis third
edition, Priadle Weber and Schmidt, Boston.
[6] Robert [Link] and Donald [Link] , Introduction to Real Ana-
lysis, Third Edition, John Wiley and Sons,Inc.
[7] Schultz,M.H , Spline Analysis, Prentice Hall, Englewood Clis,N.J.
[8] Sri Redjeki P,Dr , Metoda Matematika, Fakultas MIPA Institut Tek-
nologi Bandung,2009.
[9] Steven Chapra , Numerical Methods , McGraw-Hill.
[10] William E Boyce and Richard C Diprima , Elementary Dierential
Equations and Boundary Value Problem, John Wiley and [Link].
22