TEKNIK ITERASI VARIASIONAL UNTUK MENYELESAIKAN PERSAMAAN NONLINEAR ABSTRACT

TEKNIK ITERASI VARIASIONAL UNTUK MENYELESAIKAN PERSAMAAN NONLINEAR Koko Saputra1∗ , Supriadi Putra2 , Zulkarnain2 1

Mahasiswa Program Studi S1 Matematika Laboratorium Matematika Terapan, Jurusan Matematika Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Riau Kampus Binawidya Pekanbaru (28293), Indonesia 2

∗ [email protected]

ABSTRACT This article discusses the variational iteration technique to solve nonlinear equations. This iterative technique is obtained by estimating the second derivative that appears in the variational iteration method. Then by choosing certain functions in the form of exponents, three iterative methods of order three are found. Because the three methods require two evaluations of the function and one evaluation of the derivative function, then the efficiency index of the three methods that have been suggested is 1.442, greater than the efficiency index of Newton’s method that is 1, 414. Keywords: variational iteration, iterative method, order of convergence, Newton method, Taylor expansion. ABSTRAK Artikel ini membahas teknik iterasi variasional, untuk menyelesaikan persamaan nonlinear. Teknik iterasi ini diperoleh dengan menaksir turunan kedua yang muncul di metode iterasi variasional dan memilih beberapa fungsi tertentu yang berbentuk eksponen, sehingga ditemukan tiga metode iterasi yang memiliki orde konvergensi kubik. Karena setiap metode yang ditemukan ini memerlukan dua kali evaluasi fungsi dan satu kali evaluasi turunan fungsi, maka indek efisiensi ketiga metode yang dikemukankan ini adalah 1.442, lebih besar dari pada indek efisiensi metode Newton, yaitu 1.414. Kata kunci: iterasi variasional, metode iterasi, orde konvergensi, metode Newton, ekspansi Taylor.

JOM FMIPA Volume 1 No. 2 Oktober 2014

53

1. PENDAHULUAN Perkembangan dunia yang semakin maju membuat persoalan semakin kompleks, tidak terkecuali persoalan yang melibatkan persoalan matematika. Salah satu persoalan yang sering muncul adalah bagaimana menyelesaikan persamaan nonlinear f (x) = 0.

(1)

Dalam matematika, terdapat banyak teknik ataupun metode yang dapat digunakan untuk menyelesaikan persamaan (1), salah satunya dengan menggunakan metode Newton, dengan bentuk iterasi xn+1 = xn −

f (xn ) , f 0 (xn )

n = 0, 1, 2, · · · ,

(2)

dengan f (xn ) 6= 0, metode ini merupakan metode iterasi dengan tebakan awal x0 . Metode Newton memiliki orde konvergensi kuadratik apabila tebakan awal cukup dekat ke akar (α) [1, h. 55]. Dalam perkembangannya, metode Newton banyak mengalami modifikasi yang bertujuan untuk mempercepat iterasi dan meningkatkan orde konvergensi. Salah satu modifikasi dari metode Newton adalah modifikasi yang telah dikembangkan oleh He [2], yang merupakan modifikasi dari metode Newton dengan persamaan koreksi. Metode yang dihasilkan disebut metode iterasi He, dalam proses iterasinya metode ini menggunakan nilai turunan kedua dari fungsi f . Selanjutnya, Noor [3] memodifikasi metode iterasi yang telah dikembangkan oleh He dengan melakukan pendekatan menggunakan ekspansi Taylor yang bertujuan untuk menghindari bentuk turunan kedua yang terdapat pada metode iterasi He. Metode yang dihasilkan disebut metode iterasi Noor yang memuat fungsi sembarang g(x). Kemudian dengan melihat beberapa kasus dari fungsi sembarang g(x) yang terdapat pada metode iterasi Noor, diperoleh metode iterasi baru dengan orde konvergensi kubik. Pada artikel ini, di bagian dua dibahas proses terbentuknya metode iterasi baru yang merupakan review dari artikel [3] yang berjudul ”Variational Iteration Technique for Solving Nonlinear equations”. Pada bagian tiga dilakukan analisa kekonvergenan dari metode yang diusulkan, kemudian dilanjutkan dengan melakukan uji komputasi menggunakan program Maple 13. 2. TEKNIK ITERASI VARIASIONAL Pada bagian ini diberikan beberapa definisi dasar untuk pembahasan selanjutnya, kemudian dilanjutkan dengan proses terbentuknya beberapa metode iterasi baru. Definisi 1 (Orde Konvergensi) [5]. Misalkan barisan bilangan real {xn }∞ n=0 konvergen ke α dan nyatakan en = xn − α untuk n ≥ 0. Jika terdapat konstanta positif p > 0, maka xn+1 − α = C 6= 0, n→∞ (xn − α)p lim


54

maka barisan tersebut dikatakan konvergen ke α dengan orde kekonvergenan p. Konstanta C disebut konstanta kesalahan asimtotik (asymptotic error constant). Jika p = 2 dan 3, maka orde kekonvergenan dengan barisan bilangan real berturutturut dikenal dengan istilah kuadratik dan kubik. Definisi 2 (Persamaan Error) [5]. Apabila notasi en = xn −α merupakan notasi untuk error pada iterasi ke−n, maka . (3) en+1 = Cepn + O ep+1 n Persamaan (3) disebut sebagai persamaan error, dengan p merupakan orde konvergensi dari metode iterasi. Definisi 3 [6, h. 12]. Misalkan p adalah orde suatu metode iterasi dan m adalah banyaknya fungsi yang dievaluasi pada setiap iterasinya, maka indek efisiensi dari 1 metode iterasi itu adalah p m . 2.1 Metode Iterasi He Diberikan persamaan koreksi dengan bentuk xn+1 = xn + λf (xn ),

(4)

dengan λ merupakan pengali Lagrange. Untuk memperoleh nilai optimal dari persamaan (4) digunakan dxn+1 = 0, dxn

(5)

sehingga diperoleh λ=−

1 f 0 (x

n)

.

Kemudian, dengan mensubstitusikan nlai λ ke persamaan (4), sehingga diperoleh metode Newton. Untuk meningkatkan orde konvergensi, He [2] menambahkan fungsi sembarang g(xn ) ke dalam persamaan (4), perhatikan persamaan berikut xn+1 = φ(xn ) + λ(f (xn )g(xn ))r ,

(6)

dengan φ(xn ) adalah fungsi iterasi metode Newton dan r adalah orde dari φ(xn ). Selanjutnya dari persamaan (6) dapat dicari nilai λ, yaitu λ=−

φ0 (xn ) , r[f (xn )g(xn )]r−1 (f 0 (xn )g(xn ) + f (xn )g 0 (xn ))

kemudian, dengan mensubstitusikan nilai λ ke persamaan (6), maka diperoleh xn+1 = φ(xn ) −

φ0 (xn )(f (xn )g(xn )) . r[f 0 (xn )g(xn ) + f (xn )g 0 (xn )]


(7)

55

Kemudian, dengan menurunkan φ(xn ) terhadap xn diperoleh φ0 (xn ) =

f 00 (xn )f (xn ) , [f 0 (xn )]2

(8)

selanjutnya, dengan mensubstitusikan persamaan (2) dan (8) ke persamaan (7), setelah disederhanakan diperoleh  f (xn )   y n = xn − 0    f (xn ) 2 00 [f (xn )] g(xn )f (xn ) (9) xn+1 = yn −   r[f 0 (xn )]2 [f 0 (xn )g(xn ) + f (xn )g 0 (xn )]    n = 0, 1, 2, · · · . Persamaan (9) merupakan metode Iterasi He, metode ini menggunakan turunan kedua dalam proses perhitungannya, oleh karena itu Noor mengembangkan metode iterasi bebas dari bentuk turunan kedua. 2.2 Metode Iterasi Noor Untuk menghindari bentuk turunan kedua yang terdapat pada metode iterasi He yang terdapat pada persamaan (9), dilakukan pendekatan dengan menggunakan ekspansi Taylor untuk f (yn ) disekitar yn = xn sampai pangkat dua dan dengan mengabaikan pangkat yang lebih tinggi, dan setelah disederhanakan diperoleh f (yn ) =

1 [f (xn )]2 f 00 (xn ) . 2 [f 0 (xn )]2

(10)

Perhatikan kembali persamaan (9), dengan menggunakan bentuk persamaan (10) dan r = 2, maka persamaan (9) dapat dituliskan dalam bentuk  f (xn )   yn = xn − 0    f (xn ) f (yn )g(xn ) (11) xn+1 = yn − 0   (g (xn )f (xn ) + g(xn )f 0 (xn ))    n = 0, 1, 2, · · · . Persamaan (11) dikenal dengan metode Iterasi Noor. Metode ini memiliki orde konvergensi kubik dan setiap iterasinya memerlukan tiga kali evaluasi fungsi dan dua kali evaluasi turunan fungsi, sehingga dengan menggunakan Definisi 3 maka indek efisiensinya adalah 1.245. 2.3 Teknik Iterasi Variasional Metode Noor Metode Iterasi Noor pada persamaan (11) adalah ide pokok untuk memperoleh metode iterasi untuk menyelesaikan persamaan nonlinear dengan menggunakan beberapa kasus dari fungsi sembarang g(xn ) diantaranya, 1. g(xn ) = ±e−βxn JOM FMIPA Volume 1 No. 2 Oktober 2014

56

2. g(xn ) = ±e

β f 0 (xn )

βf (xn )   0 3. g(xn ) = ±e f (xn ) 



Berikut ini akan ditunjukkan proses terbentuknya metode iterasi baru, dengan menggunakan kasus ke−1, yaitu g(xn ) = ±e−βxn .

(12)

Persamaan (12) dapat ditulis dalam 2 bentuk, yaitu g(xn ) = e−βxn ,

(13)

g(xn ) = −e−βxn .

(14)

dan

Apabila persamaan (13) diturunkan terhadap xn diperoleh g 0 (xn ) = −β e−βxn .

(15)

Selanjutnya, dengan mensubstitusikan persamaan (13) dan (15) ke persamaan (11), setelah disederhanakan diperoleh xn+1 = yn −

f 0 (x

f (yn ) . n ) − βf (xn )

(16)

Kemudian, apabila persamaan (14) diturunkan terhadap xn diperoleh g 0 (xn ) = β e−βxn .

(17)

Selanjutnya, dengan mensubstitusikan persamaan (14) dan (17) ke persamaan (11), maka akan diperoleh xn+1 = yn − atau yn xn+1

f 0 (x

f (yn ) , n ) − βf (xn )

f (xn ) = xn − 0 , f (xn ) f (yn ) = yn − 0 , n = 0, 1, 2, · · · . f (xn ) − βf (xn )

(18)    

(19)

  

Persamaan (19) merupakan metode iterasi yang dikenal dengan Teknik Iterasi Variasional Tipe 1. Apabila persamaan (15) ataupun persamaan (17) disubstitusikan ke persamaan (11) akan memberikan hasil yang sama seperti yang terdapat pada persamaan (16) dan (18). Selanjutnya untuk kasus ke−2 dan ke−3 digunakan persamaan dalam JOM FMIPA Volume 1 No. 2 Oktober 2014

57

bentuk positif, karena akan menghasilkan bentuk yang sama seperti yang terdapat pada kasus ke−1. Selanjutnya, dengan menggunakan kasus ke−2, yaitu β

g(xn ) = e f 0 (xn ) .

(20)

Apabila persamaan (20) diturunkan terhadap xn diperoleh β  00 βf (xn ) f 0 (xn ) 0 g (xn ) = − 0 . e [f (xn )]2 

 

(21)

Selanjutnya dengan mensubstitusikan persamaan (20) dan (21) ke persamaan (11), setelah disederhanakan diperoleh  f (xn )   , y n = xn − 0    f (xn ) f (yn )f (xn ) (22) ,  xn+1 = yn − 0  f (xn )f (xn ) − 2βf (yn )    n = 0, 1, 2, · · · . Persamaan (22) merupakan metode iterasi yang dikenal dengan Teknik Iterasi Variasional Tipe 2. Kemudian, dengan menggunakan kasus ke−3 yaitu βf (xn )   0 (x ) f n g(xn ) = e . 



(23)

Apabila persamaan (23) diturunkan terhadap xn diperoleh βf (xn )   0 2 00 β([f (x )] − f (x )βf (x ) n n n f 0 (xn ) , g 0 (xn ) = e [f 0 (xn )]2 



(24)

dengan mensubstitusikan persamaan (23) dan (24) ke persamaan (11), setelah disederhanakan diperoleh  f (xn )   yn = xn − 0 ,    f (xn ) f (yn ) (25) xn+1 = yn − 0 ,   f (xn ) + β[f (xn ) − 2f (yn )]    n = 0, 1, 2, · · · . Persamaan (25) merupakan metode iterasi yang dikenal dengan Teknik Iterasi Variasional Tipe 3. Selanjutnya, akan dilakukan analisa kekonvergenan Teknik Iterasi Variasional Tipe 1, Teknik Iterasi Variasional Tipe 2 dan Teknik Iterasi Variasional Tipe 3.


58

3. ANALISA KEKONVERGENAN Teorema 4 Asumsikan bahwa fungsi f : D ⊂ R → R untuk interval buka D yang memiliki akar sederhana α ∈ D. Misalkan f adalah fungsi mulus yang berada dilingkungan akar, maka Teknik Iterasi Variasional Tipe 1 yang diberikan oleh persamaan (19) memiliki orde konvergensi kubik. Bukti: Misalkan α adalah akar dari persamaan nonlinear f (x) = 0, maka f 0 (α) 6= 0 dan en = xn − α. Dengan menggunakan ekspansi Taylor untuk f (xn ) disekitar xn = α, diperoleh f 00 (α) f 00 (α) (α − xn )2 + (α − xn )3 + O(e4n ) 2! 3! 1 1 = f (α) + f 0 (α)en + f 00 (α)en 2 + f 000 (α)en 3 + O(e4n ). (26) 2! 3!

f (xn ) = f (α) + f 0 (α)(α − xn ) +

Karena f (α) = 0, maka setelah disederhanakan persamaan (26) menjadi 1 f 000 (α)e3n 1 f 00 (α)e2n 4 0 + + O(en ) . f (xn ) = f (α) en + 2! f 0 (α) 3! f 0 (α)

(27)

Misalkan (j) 1 f (α) Cj = , j = 2, 3, · · · , j! f 0 (α)

(28)

maka persamaan (27) dapat ditulis dalam bentuk f (xn ) = f 0 (α)(en + C2 e2n + C3 e3n + O(e4n )).

(29)

Kemudian dengan menggunakan ekspansi Taylor terhadap f 0 (xn ) disekitar xn = α, maka setelah disederhanakan diperoleh f 0 (xn ) = f 0 (α) 1 + 2C2 en + 3C3 e2n + O(e3n ) . (30) Apabila persamaan (29) dibagi dengan persamaan (30), setelah disederhanakan diperoleh f (xn ) en + C2 e2n + C3 e3n + O(e4n ) = . f 0 (xn ) 1 + 2C2 en + 3C3 e2n + O(e3n )

(31)

Selanjutnya dengan menggunakan identitas geometri 1 = 1 − s + s2 − s3 + O(s4 ), 1+s

(32)

dengan s = 2C2 en + 3C3 e2n + O(e3n ), maka setelah disederhanakan persamaan (31) dapat ditulis menjadi f (xn ) = en − C2 e2n + (2C22 − 2C3 )e3n + O(e4n ). f 0 (xn ) JOM FMIPA Volume 1 No. 2 Oktober 2014

(33)

59

Selanjutnya, dengan mensubstitusikan persamaan (33) ke persamaan (19), setelah disederhanakan diperoleh yn = α + C2 e2n + 2(C3 − C22 )e3n + O(e4n ).

(34)

Kemudian dengan melakukan ekspansi Taylor untuk f (yn ) disekitar yn = α diperoleh f (yn ) = f 0 (α)(C2 e2n − (2C22 − 2C3 )e3n ) + O(e4n ).

(35)

Perhatikan kembali persamaan (19), dengan menggunakan persamaan (29), (30) dan (35) diperoleh f 0 (x

f (yn ) = C2 e2n + (βC2 − 4C22 + 2C3 )e3n + O(e4n ), n ) − βf (xn )

(36)

dengan mensubstitusikan persamaan (34) dan (36) ke persamaan (19) diperoleh xn+1 = α + (2C22 − βC2 )e3n + O(e4n ), karena en+1 = xn+1 − α, maka persamaan (36) menjadi en+1 = (2C22 − βC2 )e3n + O(e4n ).

(37)

Persamaan (37) merupakan persamaan error untuk Teknik Iterasi Variasional Tipe 1. Berdasarkan Definisi 2 maka metode ini memiliki orde konvergensi kubik. Teorema 5 Dengan menggunakan asumsikan seperti pada Teorema 4, maka Teknik Iterasi Variasional Tipe 2 seperti yang diberikan oleh persamaan (22) memiliki orde konvergensi kubik. Bukti: Perhatikan kembali persamaan (22). Dengan menggunakan persamaan (29), (30) dan (35) diperoleh 2βC22 f (yn )f (xn ) 2 2 = C2 en + − 3C2 e3n + O(e4n ), (38) f 0 (xn )f (xn ) − 2βf (yn ) f 0 (α) kemudian, dengan mensubstitusikan persamaan (34) dan (38) ke persamaan (22), setelah disederhanakan diperoleh 2βC22 3 2 xn+1 = α + 2C3 + C2 − 0 e + O(e4n ), f (α) n karena en+1 = xn+1 − α, maka 2βC22 3 2 en+1 = 2C3 + C2 − 0 e + O(e4n ). f (α) n

(39)

Berdasarkan Definisi 2, maka Teknik Iterasi variasional Tipe 2 yang diberikan oleh persamaan (22) memiliki orde konvergensi kubik. JOM FMIPA Volume 1 No. 2 Oktober 2014

60

Teorema 6 Dengan menggunakan asumsi seperti pada Teorema 4, maka Teknik Iterasi Variasional Tipe 3 seperti yang diberikan oleh persamaan (25) memiliki orde konvergensi kubik. Bukti: Perhatikan kembali persamaan (25). Dengan menggunakan persamaan (29), (30) dan (35) diperoleh f 0 (x

f (yn ) = C2 e2n + (2C3 − 4C22 − βC2 )e3n + O(e4n ). n ) + β(f (xn ) − 2f (yn ))

(40)

kemudian, dengan mensubstitusikan persamaan (34) dan (40) ke persamaan (25), setelah disederhanakan diperoleh xn+1 = α + (2C22 + βC2 )e3n + O(e4n ), karena en+1 = xn+1 − α, maka en+1 = (2C22 + βC2 )e3n + O(e4n ).

(41)

Berdasarkan Definisi 2, maka metode iterasi pada persamaan (25) memiliki orde konvergensi kubik. 4. UJI KOMPUTASI Berikut ini akan dilakukan uji komputasi untuk membandingkan kecepatan dalam menemukan akar antara metode Newton (MN), Teknik Iterasi Variasional Tipe 1 (TIV1), Teknik Iterasi Variasional Tipe 2 (TIV2) dan Teknik Iterasi Variasional Tipe 3 (TIV3). Adapun fungsi-fungsi yang digunakan adalah 1. f1 (x) = x3 + 4x2 − 10 2. f2 (x) = sin2 (x) + x2 − 1 3. f3 (x) = x2 − ex − 3x + 2 4. f4 (x) = cos(x) − x 2

5. f5 (x) = xex − sin2 x + 3 cos(x) + 5 6. f6 (x) = ex

2 +7x−30

−1

Dalam menentukan solusi numerik dari contoh-contoh fungsi nonlinear di atas, digunakan program Maple 13 dengan memakai 128 digit dan toleransi ε = 1.0 × 10−32 , dengan tebakan awal yang sama untuk masing-masing fungsi yang telah diberikan. Dalam menentukan solusi numerik juga ditentukan kriteria pemberhentian jalannya program komputasi yang sama untuk semua metode, yaitu : 1. Jika selisih nilai mutlak antara dua iterasi yang berdekatan bernilai lebih kecil dari toleransi yang diberikan. 2. Jika nilai mutlak fungsi lebih kecil dari toleransi yang diberikan. 3. Jika jumlah iterasi mencapai maksimum iterasi yang ditetapkan.


61

Hasil komputasi dapat dilihat pada tabel berikut ini. Tabel 1: Perbandingan Komputasi untuk MN, TIV1, TIV2 dan TIV3 β fi x0 xn MN TIV1 TIV2 TIV3 f1 0.7 1.3652300134140968 7 4 6 6 f2 2.3 0.7390851332151606 5 5 4 5 f3 0.0 0.2575302854398608 5 4 3 4 1.0 f4 1.7 0.7390851332151606 6 5 4 4 f5 −1.0 −1.2076478271309189 7 5 5 4 f6 3.5 3.0000000000000000 13 9 9 9 f1 0.7 1.3652300134140968 7 5 7 6 f2 2.3 0.7390851332151606 5 5 4 5 f3 0.0 0.2575302854398608 5 4 3 3 0.5 f4 1.7 0.7390851332151606 6 4 4 4 f5 −1.0 −1.2076478271309189 7 5 5 5 f6 3.5 3.0000000000000000 13 9 9 9

Tabel 1, β merupakan parameter dengan dua nilai yang berbeda, sedangkan xn merupakan akar hampiran dari fungsi. Secara umum ke empat metode yang dibahas berhasil menemukan akar yang diharapkan, akan tetapi metode baru yang diusulkan diantaranya Teknik Iterasi Variasional Tipe 1, Teknik Iterasi Variasional Tipe 2 dan Teknik Iterasi Variasional Tipe 3 lebih cepat menemukan akar yang diharapkan dibandingkan dengan metode Newton. Hal ini ditunjukkan dari jumlah iterasi metode baru lebih sedikit dibandingkan metode Newton. Dari Tabel 1 terlihat bahwa, parameter berpengaruh terhadap jumlah iterasi yang dihasilkan, baik itu untuk Teknik Iterasi Variasional Tipe 1, Teknik Iterasi Variasional Tipe 2 maupun untuk Teknik Iterasi Variasional Tipe 3. Hal ini dapat dilihat pada f1 untuk β = 1.0 maupun untuk β = 0.5. Akan tetapi, pengaruh parameter akan sangat kecil sekali, bahkan mungkin tidak ada apabila tebakan awal diambil cukup dekat ke akar, contohnya dapat dilihat pada fungsi f6 . Kemudian, dari analisa kekonvergenan yang telah dilakukan pada bagian sebelumnya dapat dilihat bahwa, ketiga metode yang diusulkan memiliki orde konvergensi kubik. Ketiga metode ini memerlukan dua kali evaluasi fungsi dan satu kali evaluasi turunan fungsi pada setiap iterasinya, sehingga indeks efisiensi dari ketiga metode yang diusulkan adalah 1.442. Sedangkan metode Newton memiliki orde konvergensi kuadratik dengan satu kali evaluasi fungsi dan satu kali evaluasi turunan fungsi pada setiap iterasinya, sehingga indeks efisiensi dari metode Newton adalah 1.414. Secara umum keempat metode yang dibahas berhasil menemukan akar yang diharapkan, akan tetapi metode baru yang diusulkan lebih unggul dari pada metode Newton, baik itu dilihat dari jumlah iterasi, indeks efisiensi maupun jika dilihat dari orde konvergensi metode itu sendiri.


62

DAFTAR PUSTAKA [1] Faires, J.D. & R.L. Burden. 2002. Numerical Methods, Third Edition. Brooks Cole, New York. [2] He, J.H. 2007. Variational Iteration Method-Some Recent Results and New Interpretations. Applied Mathematics Computation, 207: 3-17. [3] Noor, M.A. 2007. Variational Iteration Technique for Solving Nonlinear Equations. Applied Mathematics Computation, 321: 247-131. [4] Ostrowski, A.M. 1966. Solution of Equations and Systems of Equations, Third Edition. Academic Press, New York. [5] Sharma, J.R., R.K. Guha & R. Sharma. 2011. Some Modified Newton’s Methods with fourth-order Convergence. Applied Science Research, 1:240-247. [6] Traub, J.F. 1964. Iterative Methods for the Solution of Equations, Prentice Hall Inc. Englewood Cliffs. New Jersey.


63

TEKNIK ITERASI VARIASIONAL UNTUK MENYELESAIKAN PERSAMAAN NONLINEAR ABSTRACT

Recommend Documents