Jurnal Matematika Murni dan Terapan “εpsilon” Vol.10 No.1 (2016) Hal. 14 - 20

## POLINOMIAL CHEBYSHEV PADA SYARAT BATAS SERAP GELOMBANG AKUSTIK DUA DIMENSI

## Mohammad Mahfuzh Shiddiq
Program Studi Matematika Fakultas MIPA Universitas Lambung Mangkurat E-mail: mmahfuzhs@unlam.ac.id

## Abstrak
Masalah syarat batas pada persamaan gelombang mempunyai banyak jenis dan metode penyelesaian. Salah satu masalah syarat batas adalah syarat batas serap pada persamaan gelombang. Syarat batas serap ini muncul sebagai akibat domain natural pada permasalah perambatan gelombang tidak terbatas dan membutuhkan perhitungan yang besar. Metode penyelesaian secara numerik tidak bisa dihindari pada persamaan gelombang jenis ini. Solusi numerik yang akan dibahas pada makalah ini adalah dengan mendekati solusi masalah perambatan gelombang akustik dua dimensi dengan menggunakan polinomial chebyshev. Beberapa perbandingan hasil solusi dengan mengg8unakan pendekatan yang lain yang sudah dikerjakan juga diberikan sehingga menunjukkan efektifitas solusi mana yang lebih baik.
Kata kunci : Polinomial Chebyshev, Solusi Numerik, Gelombang Akustik dua dimensi

## 1. PENDAHULUAN
Masalah perambatan gelombang sangat menarik untuk diteliti. Hal ini menjadi dasar teori pada cabang ilmu aplikasi yang menggunakan gelombang sebagai alat. Salah satu bidang yang menggunakan adalah pemodelan gelombang air laut, tsunami, geofisika, seismik, pencitraan medik dan lain sebagainya [5,3]. Pemodelan perambatan gelombang biasa menggunakan domain natural yang tidak terbatas serta melibatkan perhitungan yang besar. Hal ini menyebabkan banyak peneliti menggunakan syarat batas buatan untuk menyederhanakan tapi masih bisa dipakai mensimulasikan keadaan sebenarnya [1,4]  Syarat batas buatan pada perambatan gelombang untuk mengatasi permasalahan tersebut diantaranya adalah syarat batas serap. Syarat batas serap ini memungkinkan gelombang merambat di domain yang terbatas akan meneruskan perambatannya di batas- batas domain, dengan kata lain meminimkan gelombang pantulan di batas domain atau bahkan menghilangkan gelombang pantulan di batas.
Pada penelitian sebelumnya yang dilakukan oleh Shiddiq [7]  telah menunjukkan bahwa pada masalah perambatan gelombang akustik dua dimensi bisa meminimalkan gelombang pantul di batas domain. Selanjutnya Shiddiq [8] menyelesaikan masalah tersebut dengan menggunakan pendekatan numerik dengan menggunakan persamaan beda hingga. Beberapa peneliti yang lain juga pernah melakukan dengan pendekatan yang berbeda [2,10]  Pada paper ini akan ditentukan solusi dari masalah pada Shiddiq [7] dengan pendekatan numerik dan memakai polynomial Chebysev untuk mendekati syarat batas serap yang didefinisikan pada masalah perambatan gelombang.

## 2. TINJAUAN PUSTAKA
Persamaan gelombang akustik adalah tipe gelombang longitudinal yang merambat dengan kompresi dan dekompresi adibiatik. Persamaan gelombang merambat dengan kecepatan suara yang bergantung pada media dimana gelombang tersebut merambat. Hal ini terjadi misalkan pada defraksi, refleksi dan interfensi.
Misalkan diberikan persamaan gelombang dua dimensi sebagai berikut
2 2 2 1 0 2 2 2 2 u u u c t x y ∂ ∂ ∂ − + = ∂ ∂ ∂         (1)
Didefinisikan di dalam media terbatas  ( ) { } , , , 0 , 0 D x y t L x L y M t T = − ≤ ≤ ≤ ≤ ≤ ≤ dan c k ω = menotasikan kecepatan di media, ω adalah frekuensi sudut dan k adalah bilangan gelombang. Persamaan gelombang (1) dengan nilai awal
( ) ( ) ( ) , , 0 , , , 0 0 u u x y F x x y t ∂ = = ∂ (2)
Masalah nilai awal persamaan gelombang (1)  dan  (2)  dilengkapi dengan syarat batas Dirichlet, yaitu
( ) ( ) , , 0, , , 0
u L y t u x M t ± = =
atau syarat batas Neumann, yaitu
( ) ( ) , , , , 0, 0 u L y t u x M t x y ∂ ± ∂ = =
∂ ∂
menghasilkan gelombang pantul pada batas dengan amplitudo sama dengan gelombang datang, yaitu koefisien pantul sama dengan 1. Koefisien pantul ini didapatkan dari mensubtitusikan solusi dari persamaan gelombang (1), yaitu
( ) ( ) cos sin cos sin , 0 2 i t kx ky i t kx ky u e Re ω θ θ ω θ θ π θ − ± + ±   = + ≤ ≤     (3)
dengan θ adalah sudut diantara bidang gelombang dengan sumbu x dan R adalah koefisien pantul.
Polinomial Chebyshev memiliki banyak macam dan jenis. Salah satunya adalah polinomial jenis pertama dan kedua, ( ) n T x dan ( ) n U x . Polinomial Chebyshev jenis pertama biasa digunakan untuk fungsi yang didefinisikan di dalam interval [ ] 1,1 − . Sedangkan untuk polinomial Chebyshev jenis kedua digunakan untuk interval yang lebih umum, yaitu [ ] , a b .
Definisi 1 Polinomial Chebyshev ( ) n T x jenis pertama adalah polinomial dalam x dengan derajat n dan didefinisikan dengan
Polinomial Chebyshev pada  Syarat Batas Serap Gelombang Akuatik Dua Dimensi
( ) cos , cos n
T x n ketika x θ θ = = (4)
Jika nilai x berada di dalam interval [ ] 1,1 − maka nilai variabel θ berada di dalam interval [ ] 0, π . Berdasarkan persamaan (4)  maka dapat dilihat empat polinomial Chebyshev pertama sebagai berikut
( ) ( ) ( ) ( ) ( ) 2 0 1 2 3 4 2 3 4 1, , 2 1 4 3 , 8 8 1,
T x T x x T x x
T x x x T x x x = = = − = − = − + 
Polinomial Chebyshev juga bisa didefinisikan secara rekursif jika dilihat dari (4) dan sifat identitas trigonometri, yaitu sebagai berikut
( ) ( ) ( ) 1 2 2 , 2,3, , n n n T x xT x T x n − − = − =  dengan nilai awal ( ) ( ) 0 1 1, T x T x x = =
[6].
Pendekatan polinomial dengan derajat n untuk fungsi kontinu ( ) f x yang diberikan pada [ ] 1,1 − adalah dengan menginterpolasi diantara nilai dari ( ) f x di 1 n + titik yang dipilih secara benar di dalam interval [ ] 1,1 − . Misalkan untuk interpolasi di titik
1 2 1 , , , n x x x + 
dengan polinomial ( ) 0 1 n n n p x c c x c x = + + + 
dibutuhkan
( ) ( ) 0 1 , 1, , n 1 n k n k k c c x c x f x k + + + = = +  
Sedangkan pilihan yang baik dalam interpolasi titik untuk menjamin kekonvergenan seragam adalah himpunan zeroes dari polinomial Chebyshev ( ) 1 n T x + , yaitu
( ) 1 2 1 cos , 1, , 1 k k n x x k n π   −     + = = = =  (5)
Jadi polinomial derajat n , ( ) n p x , menginterpolasi ( ) f x di dalam titik (5) dalam bentuk jumlahan polinomial Chebyshev adalah sebagai berikut
( ) ( ) 0 n n i i i p x c T x = = ∑
dengan koefisien i c diberikan oleh ( ) ( ) 1 1 2 1 n i k i k k c f x T x n + = = + ∑
Teorema 2.  Jika ( ) f x adalah sebarang fungsi di dalam interval [ ] 1,1 − dan jika koefisien , 0,1, , 1 j c j n = −  didefinisikan dengan ( ) ( ) 1 0 1 1 1 2 2 0 2 2 cos cos n j k j k k n k j k n n k c f x T x n f n π π − =         − + +                         =     =   =       ∑ ∑ Maka formula aproksimasi fungsi ( ) f x ( ) ( ) 1 0 0 1 2 n k k k f x c T x c − =   ≈ −     ∑
adalah eksak untuk x sama dengan ke semua n zeroes ( ) n T x .

## 3. PEMBAHASAN

## Tinjau kembali persamaan gelombang pada (1), yaitu
2 2 2 1 0 2 2 2 2 u u u c t x y ∂ ∂ ∂ − + = ∂ ∂ ∂        
dengan salah satu kelas dari solusi (1)  yang direpresentasikan oleh gelombang yang merambat ke kanan diberikan oleh
( ) ( ) cos sin cos sin i t kx ky i t kx ky u e Re ω θ θ ω θ θ − ± + ± = +
Jika variabel ( ) , y t dibuat tetap, maka salah satu syarat batas orde satu
cos 0 x L ik u x θ = ∂   +
=   ∂   (6)
akan menyebabkan gelombang diteruskan pada batas x L = dan semua gelombang (3) tidak menghasilkan gelombang pantul, yaitu 0 R = .
Pada kasus yang umum, tidak mungkin mendesain syarat batas sedemikian sehingga pantulan tidak terjadi, sehingga syarat batas yang mungkin dikontruksi adalah syarat batas yang mereduksi sebanyak mungkin sejumlah pantulan yang terjadi. Oleh karena itu, akan dicari dan ditaksir syarat batas pada (6) sedemikian sehingga koefisien pantul dapat dibuat seminim mungkin atau dengan kata lain amplitudo gelombang pantul mendekati nol.

## Polinomial Chebyshev pada  Syarat Batas Serap Gelombang Akuatik Dua Dimensi
Ambil simbol syarat batas pada (6) yang diberikan dengan
2 cos 1 sin ik ik x x θ θ ∂ ∂ + = + − ∂ ∂ (7)
Jika menggunakan aproksimasi polinomial Chebyshev derajat nol untuk fungsi akar pangkat dua yaitu 1 1 x − = maka diperoleh 2 1 sin ik ik x x θ ∂ ∂ + − = + ∂ ∂
dengan mengingat bahwa 1 ik i c ω = berkaitan dengan 1 c t ∂ ∂ sehingga didapatkan hampiran pertama syarat batas serap
1 0 u x c t ∂ ∂   + =
  ∂ ∂   (8)
Selanjutnya akan ditentukan koefisien pantul dari syarat batas serap (8)  dengan mensubstitusikan ke (3) sehingga diperoleh
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) cos sin cos sin cos sin cos sin cos sin cos sin cos sin cos sin 1 1 0 cos i t kx ky i t kx ky i t kx ky i t kx ky i t kx ky i t kx ky i t kx ky i t kx ky u e Re x c t x e Re c t ik e Re i e Re c ω θ θ ω θ θ ω θ θ ω θ θ ω θ θ ω θ θ ω θ θ ω θ θ θ ω − ± + ± − ± + ± − ± + ± − ± + ± ∂ ∂ ∂   + = +   ∂ ∂ ∂   ∂ + + ∂ = − + + +
Kedua ruas dibagi dengan ik dan didapatkan
( ) 1 1 cos 1 cos R R θ θ θ − = = + (9)
Syarat batas serap (8)  dapat diketahui bahwa bekerja hanya pada kasus satu dimensi. Jadi akan ditentukan lagi syarat batas hampiran dengan menggunakan derajat yang lebih tinggi.
Jika ditinjau kembali aproksimasi polinomial Chebyshev derajat satu
untuk fungsi yaitu 1 0.9239 0.5412 x x − = − maka dapat dilihat bahwa ( ) 2 2 2 2 2 2 1 sin 0.9239 0.5412sin 0
0.9239 0.5412 sin 0 0.9239 0.5412 sin ik ik x x ik ik x ik k k x θ θ θ θ ∂ ∂ + − = + − ∂ ∂ ∂ = + − ∂ ∂ = − + ∂
Vol.10 No.1 (2016) Hal. 14 - 20
Diketahui bahwa ik , 2 k − dan 2 2 sin k θ − berturut-turut berkaitan dengan 1 c t ∂ ∂ ,
2 2 2 1 c t ∂ ∂ dan 2 2 y ∂ ∂ sehingga menghasilkan hampiran kedua syarat batas serap 2 2 2 2 2 2 1 1 0.9239 0.5412 0 u c x t c t y   ∂ ∂ ∂ + − =   ∂ ∂ ∂ ∂   (10)
Dengan mensubtitusi persamaan (10) ke dalam persamaan (3) diperoleh rumusan koefisien pantul sebagai berikut
( ) 2 2 2 cos 0.9239 0.5412sin cos 0.9239 0.5412sin R R θ θ θ θ θ − + = = + − (11)
Jika ditinjau kembali aproksimasi polinomial Chebyshev derajat dua untuk fungsi akar kuadrat, yaitu
( ) 2 2 1 0.9107 0.5774 0.0893 2 1 1 0.5774 0.1786 x x x x x − = − − − = − −
dan diperoleh
( ) ( ) ( ) 2 2 4 2 4 4 5 3 2 2 4 4 1 sin 1 0.5774sin 0.1786sin
0
0.5774 sin 0.1786 sin 0 0.5774 sin 0.1786 sin
ik ik x x ik ik ik x k ik ik k ik k x θ θ θ θ θ
θ θ ∂ ∂ + − = + − − ∂ ∂ ∂ = + − − ∂ ∂ = − + + ∂
Dengan cara yang serupa diperoleh hampiran ketiga syarat batas serap
5 5 5 5 4 4 5 5 3 3 2 4 1 1 1 1
0.5774 0.1786 0 u c t x c t c t y c t y 
 ∂ ∂
∂ ∂ − − + =   ∂ ∂ ∂ ∂ ∂ ∂ ∂   (12)
dan koefisien pantul mempunyai rumusan sebagai berikut
( ) 2 3 2 cos 1 0.756sin cos 1 0.756sin R R θ θ θ θ θ − + = = + − (13)
Nilai dari koefisien pantul untuk aproksimasi polinomial Chebyshev derajat 0,1, 2 n = pada sebagian sudut dapat dilihat pada tabel 1 sebagai berikut
Tabel 1. Koefisien Pantul dengan polinomial Chebyshev
θ 0 6 π 4 π 3 π 2 π ( ) 1 R θ 0 0.0718  0.1716  0.3333 1 ( ) 2 R θ 0.0396  0.0468  0.0396  -0.0177  -1 ( ) 3 R θ 0 0.0328  0.0640  0.0718  -1

## Polinomial Chebyshev pada  Syarat Batas Serap Gelombang Akuatik Dua Dimensi
Sedangkan nilai koefisien dari hampiran syarat batas serap dengan penaksir Padé [9] dapat dilihat pada tabel 2 sebagai berikut
Tabel 2. Koefisien Pantul dengan penaksir Padé
θ 0 6 π 4 π 3 π 2 π  ( ) 1 R θ 0  0.0718  0.1716  0.3333  1
 ( ) 2 R θ 0  0.0004  0.0051  0.0370  1
 ( ) 3 R θ 0  0.1221  0.2886  0.4999  1

## 4. KESIMPULAN
Berdasarkan tabel 1 diperoleh rata-rata terkecil dari nilai koefisien pantul terjadi pada hampiran ketiga syarat batas serap yaitu dengan menggunakan polinomial Chebyshev derajat 2. Sedangkan jika dibandingkan dengan menggunakan penaksir Padé pada tabel 2, hampiran dengan menggunakan polinomial Chebyshev derajat dua tetap memiliki rata-rata terkecil.
5. DAFTAR RUJUKAN
[1] Brunner H dan Han H. 2014. Artificial Boundary Condition and Finite Difference Approximations for A Time-Fractional Wave Equationon on A two-dimensional unbounded spatial domain. Journal of Computational Physics, 276: 541-562
[2] Higdon R. 1987. Numerical absorbing boundary conditions for wave equation. Math. Comp.  49. 65–90
[3] Hu, J., dan Jia X. 2016. Numerical Modeling of Seismic Wave Using Frequency-Adaptive Meshes. Journal of Applied Geophysics. 131:41-53
[4] Lin X dan Yu X. 2015. A finite difference method for effective treatment of mild-slope wave equation subject to non-reflecting boundary conditions. Applied Ocean Research, 53: 179-189
[5] Mousavi, M. R. Karimi, M. dan Jamshidi, A. 2016. Probability Distribution of Acoustic Scattering from Slightly Rough Sea Surface. Ocean Engineering, 112: 134-144
[6] Mason J.C dan Handscomb D. 2002. Chebyshev Polynomials. A CRC Press Company. Florida
[7] Shiddiq, M. Mahfuzh. 2011. Syarat batas serap pada gelombang akustik dua dimensi. Jurnal Matematika Murni dan Terapan Epsilon. 05 (2)
[8] Shiddiq, M. Mahfuzh. 2013. Numerical solution of absorbing boundary conditions on two dimensional acoustic wave. Jurnal Matematika Murni dan Terapan Epsilon. 07 (1)
[9] Shiddiq, M. 2015. Numerical Solution of Absorbing Boundary Condition in Padé Approximation. The 2015 International Conference on Mathematics, Its Application, and Mathematics Education. Yogyakarta.
[10] Rowley C and Colonius T. 2000. Discretely nonreflecting boundary condition for linear hyperbolic system. J. Comput. Phys. 15. 500–538