Created By Aristastory.Wordpress.com BAB I PENDAHULUAN. Teori sistem dinamik adalah bidang matematika terapan yang digunakan untuk

Created By Aristastory.Wordpress.com

BAB I PENDAHULUAN

1.1

Latar Belakang Teori sistem dinamik adalah bidang matematika terapan yang digunakan untuk

memeriksa

kelakuan

sistem

dinamik

kompleks,

biasanya

dengan

menggunakan persamaan diferensial ataupun persamaan beda. Bila digunakan persamaan diferensial, teori tersebut dinamakan sistem dinamik kontinu. Bila digunakan persamaan beda, teori tersebut dinamakan sistem dinamik diskret. Bila variabel waktu berjalan dalam himpunan yang diskret pada beberapa selang dan kontinu pada selang lain, atau himpunan-waktu lain seperti himpunan Cantor, maka kita mendapatkan persamaan dinamik pada skala waktu. Beberapa keadaan juga mungkin dimodelkan oleh operator campuran seperti persamaan diferensial-beda. Teori ini membahas kelakuan kualitatif jangka panjang sistem dinamik, dan studi pemecahan persamaan gerak dari sistem yang terutama bersifat mekanis di alam. Ini termasuk

orbit

planet,

kelakuan

rangkaian

elektronik,

dan

pemecahan

terhadap persamaan diferensial parsial yang muncul dalam biologi. Kebanyakan penelitian modern memusatkan perhatian kepada studi sistem khaos (Wikipedia). Sistem dinamik kontinu seperti disebutkan diatas, sangat membantu untuk menyelesaikan persamaan-persamaan dengan variabel parameter yang saling berhubungan dan sering digunakan sebagai solusi penyelesaian dari beberapa model

1


matematika, khusunya bidang fisika, kimia, dan juga biologi. Dalam sistem dinamik, akan diketahui perilaku sistem yang diberikan. Beberapa metode yang dapat digunakan untuk menganalisa kestabilan sistem yaitu melalui nilai eigen, metode Lyapunov, limit cycle dan bifurkasi. Melalui nilai eigen akan diperoleh informasi kestabilan sistem, baik menggunakan sistem persamaan diferensial linier maupun sistem almost linier. 1.2

Rumusan Masalah Rumusan masalah dalam penelitian ini adalah analisa kestabilan suatu persamaan

linier 2 dimensi melalui nilai eigennya. 1.3

Tujuan Penelitian Adapun tujuan yang ingin dicapai dalam penelitian ini adalah mengetahui

kestabilan suatu persamaan linier 2 dimensi melalui nilai eigennya. 1.4

Batasan Masalah Dalam penelitian ini, penulis hanya memfokuskan pada analisa kestabilan sistem

dengan persamaan linier 2 dimensi dengan sistem dinamik kontinu. Penelitian ini juga hanya menganalisa kestabilan sistem melalui nilai eigennya.

2


BAB II TINJAUAN PUSTAKA 2.1

Sistem Dinamik Secara umum sistem dinamik adalah pemodelan sebuah masalah nyata secara

matematis dengan menggunakan persamaan diferensial yang didalam persamaannya mengandung parameter-parameter yang saling berhubungan. Perubahan pada parameter-parameter tersebut akan mengakibatkan perubahan kestabilan pada titik equilibrium. Sistem dinamik pada 𝐸 adalah pemetaan ∅ = 𝐶 1 𝜙: ℝ × 𝐸 → 𝐸 dengan 𝐸 adalah himpunan bagian terbuka dari ℝ𝑛 dan jika 𝜙𝑡 (𝑥) = 𝜙(𝑡, 𝑥), maka 𝜙𝑡 memenuhi : (i)

𝜙0 (𝑥) = 𝑥, ∀𝑥 ∈ 𝐸 dan

(ii)

𝜙𝑡 ∘ 𝜙𝑠 (𝑥) = 𝜙𝑡+𝑠 (𝑥), ∀𝑥 ∈ 𝐸 dan 𝑡, 𝑠 ∈ ℝ (Perko, 2000:182) Secara geometri, sistem dinamik menggambarkan pergerakan titik-titik di dalam

ruang fase sepanjang kurva-kurva solusi dari sistem persamaan diferensialnya. Secara grafik, sistem dinamik akan memunculkan orbit. 1.1.1 Orbit Orbit melalui 𝑥0 , dinotasikan sebagai 𝑂𝑟(𝑥0 ), adalah himpunan titik-titik 𝑥 dalam ruang keadaan 𝑋 yang berada pada suatu flow sehingga 𝑥 = 𝜑 𝑡 𝑥0 , yakni 𝑂𝑟(𝑥0 ) = {𝑥 ∈ 𝑋 ∶ 𝑥 = 𝜑 𝑡 𝑥0 , 𝑡 ∈ 𝑇} (Wiggins, 1990:2)

3


Namun, tidak semua sistem persamaan diferensial dapat ditentukan solusi dari sistemnya. Sehingga, tujuan utama sistem dinamik adalah mempelajari perilaku sistem di sekitar titik equilibrium. 1.1.2 Titik Equilibrium Titik 𝑥0 ∈ 𝑋 dikatakan titik equilibrium jika memenuhi 𝜑 𝑡 𝑥0 = 𝑥0 untuk semua 𝑡 ∈ 𝑇 (Kuznetsov, 1990:9). Analisis kestabilan adalah suatu pendekatan yang dapat digunakan untuk mempelajari perilaku dari sistem. Analisis kestabilan dapat dilakukan dengan berbagai cara seperti penyelidikan terhadap perilaku titik setimbang persamaan diferensialnya. 1.1.3 Potret Fase Gabungan dari beberapa orbit yang digambarkan dalam satu bidang disebut potret fase atau biasa disebut bidang fase. Potret fase digambarkan dengan memanfaatkan nilai eigen dan vektor eigen dari persamaan yang dianalisis. 2.2

Nilai Eigen dan Vektor Eigen Nilai eigen dan vektor eigen sistem dinamik dapat digunakan untuk mempelajari

keadaan dinamik dari suatu sistem khususnya sistem linear. Misalkan 𝐴 adalah matriks 𝑛 × 𝑛, maka vektor 𝑥 yang tidak nol di ℝ𝑛 disebut vektor eigen (eigen vektor) dari 𝐴 jika 𝐴𝑥 adalah kelipatan skalar dari 𝑥, yaitu 𝐴𝑥 = 𝜆𝑥

4


untuk 𝜆 suatu skalar. Skalar 𝜆 dinamakan nilai eigen (eigen value) dari 𝐴 (Anton, 1988:277). Persamaan 𝐴𝑥 = 𝜆𝑥 dapat dituliskan sebagai berikut 𝐴𝑥 = 𝜆𝑥 𝐴𝑥 − 𝜆𝑥 = 0 (𝐴 − 𝜆𝐼)𝑥 = 0 Persamaan (𝐴 − 𝜆𝐼)𝑥 = 0 memiliki pemecahan tak nol jika dan hanya jika (𝐴 − 𝜆𝐼) = 0 tidak memiliki invers, akibatnya 𝑑𝑒𝑡(𝐴 − 𝜆𝐼) = 0. 2.3

Kestabilan Sistem Kestabilan suatu sistem linear dapat dilihat dari nilai eigen sistem tersebut. Pada persamaan diferensial orde satu 𝑥̇ = 𝑓(𝑥) dengan solusi awal 𝑥(𝑡, 𝑥0 ) pada

waktu 𝑡 dan dengan kondisi awal 𝑥(0) = 𝑥0 , pernyataan berikut bernilai benar a. Suatu nilai 𝑥 dimana memenuhi 𝑓(𝑥) = 0 maka nilai 𝑥̅ disebut sebagai titik ekuilibrium. b. Titik ekuilibrium 𝑥̅ dikatakan stabil jika untuk setiap 𝜀 > 0 dan 𝛿 > 0 , sedemikian hingga jika ‖𝑥0 − 𝑥̅ ‖ < 𝛿 maka ‖𝑥(𝑡, 𝑥0 ) − 𝑥̅ ‖ < 𝜀 untuk setiap 𝑡 ≥ 0 c. Titik ekuilibrium 𝑥̅ dikatakan stabil asimtotis jika titik ekuilibrium tersebut stabil dan selain itu untuk 𝛿1 > 0, sedemikian hingga lim ‖𝑥(𝑡, 𝑥0 ) − 𝑥̅ ‖ = 0 dengan 𝑡→∞

ketentuan bahwa ‖𝑥0 − 𝑥̅ ‖ < 𝛿1 .

5


d. Titik ekuilibrium 𝑥̅ tidak stabil jika untuk setiap 𝜀 > 0 ada 𝛿 > 0 sedemikian sehingga, jika ‖𝑥0 − 𝑥̅ ‖ < 𝛿, maka ‖𝑥(𝑡, 𝑥0 ) − 𝑥̅ ‖ > 𝜀 untuk semua 𝑡 ≥ 0 (Olsder, 2004:57)

6


BAB III HASIL DAN PEMBAHASAN 3.1

Hasil Tulisan ini merupakan kajian terhadap bentuk-bentuk kestabilan sistem melalui

persamaan diferensial. Adapun sistem persamaan diferensial yang digunakan dalam penelitian ini adalah sistem persamaan diferensial linier 2 dimensi. Berikut beberapa contoh persamaan diferensial linier 2 dimensi dan pengecekan nilai eigen serta kestabilannya. 1.

𝑑𝑥 𝑑𝑡

𝑑𝑦

= 2𝑥 + 𝑦 dan

𝑑𝑡

= 2𝑥 + 3𝑦

Penyelesaian : Titik kritis : (𝑥, 𝑦) = (0,0) 𝐴=(

2 1 ) 2 3 1 ) −1

𝜆1 = 1 → ⃗⃗⃗⃗ 𝑣1 = (

1 𝜆2 = 4 → ⃗⃗⃗⃗ 𝑣2 = ( ) 2 [

𝑥(𝑡) 1 1 ] = 𝐶1 𝑒 𝑡 ( ) + 𝐶2 𝑒 4𝑡 ( ) 𝑦(𝑡) −1 2

∴ Tidak stabil dengan tipe titik kestabilan adalah node

7


Gambar 3.1 Potret fase (tidak stabil dengan tipe titik kestabilan node)

2.

𝑑𝑥 𝑑𝑡

= 𝑥 + 2𝑦 dan

𝑑𝑦 𝑑𝑡

= 2𝑥 − 2𝑦


1 2 ) 2 −2 1 ) −2

𝜆1 = −3 → ⃗⃗⃗⃗ 𝑣1 = (

2 𝜆2 = 2 → ⃗⃗⃗⃗ 𝑣2 = ( ) 1 [

𝑥(𝑡) 1 2 ] = 𝐶1 𝑒 −3𝑡 ( ) + 𝐶2 𝑒 2𝑡 ( ) 𝑦(𝑡) −2 1

∴ Tidak stabil dengan tipe titik kestabilan adalah titik pelana (saddle point)

8


Gambar 3.2 Potret fase (tidak stabil dengan tipe titik kestabilan pelana)

3.

𝑑𝑥 𝑑𝑡

= 3𝑥 dan

𝑑𝑦 𝑑𝑡

= 3𝑦


3 0 ) 0 3

1 𝜆1 = 𝜆2 = 3 → ⃗⃗⃗⃗ 𝑣1 = ⃗⃗⃗⃗ 𝑣2 = ( ) 0 [

𝑥(𝑡) 1 1 ] = 𝐶1 𝑒 3𝑡 ( ) + 𝐶2 𝑒 3𝑡 𝑥 ( ) 𝑦(𝑡) 0 0

∴ Tidak stabil dengan tipe titik kestabilan adalah node

9


Gambar 3.3 Potret fase (tidak stabil dengan tipe titik kestabilan node)

4.

𝑑𝑥 𝑑𝑡

= −5𝑥 − 𝑦 dan

𝑑𝑦 𝑑𝑡

= 𝑥 − 3𝑦


−5 −1 ) 1 −3 1 ) −9

𝜆1 = 𝜆2 = −4 → ⃗⃗⃗⃗ 𝑣1 = ⃗⃗⃗⃗ 𝑣2 = (

[

𝑥(𝑡) 1 1 ] = 𝐶1 𝑒 −4𝑡 ( ) + 𝐶2 𝑒 −4𝑡 𝑥 ( ) 𝑦(𝑡) −9 −9

∴ Stabil asimtotik dengan tipe titik kestabilan adalah node

10


Gambar 3.4 Potret fase (stabil asimtotik dengan tipe titik kestabilan node)

5.

𝑑𝑥 𝑑𝑡

= 𝑦 dan

𝑑𝑦 𝑑𝑡

= −50𝑥 − 15𝑦


0 1 ) −50 −15 1 ) −5

𝜆1 = −5 → ⃗⃗⃗⃗ 𝑣1 = (

1 ) −10

𝜆2 = −10 → ⃗⃗⃗⃗ 𝑣2 = (

[

𝑥(𝑡) 1 1 ] = 𝐶1 𝑒 −5𝑡 ( ) + 𝐶2 𝑒 −10𝑡 ( ) 𝑦(𝑡) −5 −10

∴ Stabil asimtotik dengan tipe titik kestabilan adalah node

11


Gambar 3.5 Potret fase (stabil asimtotik dengan tipe titik kestabilan node)

6.

𝑑𝑥 𝑑𝑡

= −4𝑥 + 5𝑦 dan

𝑑𝑦 𝑑𝑡

= −5𝑥 + 2𝑦


−4 5 ) −5 2 3 − 4𝑖 ) 5

𝜆1 = −1 + 4𝑖 → ⃗⃗⃗⃗ 𝑣1 = (

3 + 4𝑖 ) 5

𝜆2 = −1 − 4𝑖 → ⃗⃗⃗⃗ 𝑣2 = (

𝑥(𝑡) = 3𝑒 −𝑡 cos 4𝑡 + 4𝑒 −𝑡 sin 4𝑡 𝑦(𝑡) = 5𝑒 −𝑡 cos 4𝑡 ∴ Stabil asimtotik dengan tipe titik kestabilan adalah spiral point

12


Gambar 3.6 Potret fase (stabil asimtotik dengan tipe titik kestabilan spiral point)

7.

𝑑𝑥 𝑑𝑡

= 4𝑥 + 5𝑦 dan

𝑑𝑦 𝑑𝑡

= −5𝑥 − 2𝑦

Penyelesaian : Titik kritis : (𝑥, 𝑦) = (0,0) −4 5 𝐴=( ) −5 −2 −3 − 4𝑖 ) 5

𝜆1 = 1 + 4𝑖 → ⃗⃗⃗⃗ 𝑣1 = (

−3 + 4𝑖 𝜆2 = 1 − 4𝑖 → ⃗⃗⃗⃗ 𝑣2 = ( ) 5 𝑥(𝑡) = −3𝑒 −𝑡 cos 4𝑡 − 4𝑒 −𝑡 sin 4𝑡 𝑦(𝑡) = 5𝑒 −𝑡 cos 4𝑡 ∴ Tidak stabil dengan tipe titik kestabilan adalah spiral point

13


Gambar 3.7 Potret fase (tidak stabil dengan tipe titik kestabilan spiral point)

8.

𝑑𝑥 𝑑𝑡

= 𝑦 dan

𝑑𝑦 𝑑𝑡

= −𝑥


0 1 ) −1 0

1 𝜆1 = 𝑖 → ⃗⃗⃗⃗ 𝑣1 = ( ) 𝑖 1 ) −𝑖

𝜆2 = −𝑖 → ⃗⃗⃗⃗ 𝑣2 = ( [

𝑥(𝑡) 1 1 ] = 𝐴 cos 𝑥 ( ) + 𝐵 sin 𝑥 ( ) 𝑦(𝑡) 𝑖 −𝑖

∴ Stabil dengan tipe titik kestabilan adalah center

14


Gambar 3.8 Potret fase (stabil dengan tipe titik kestabilan center)

3.2

Pembahasan Beberapa contoh persamaan diferensial diatas memberikan nilai eigen dan

vektor eigen yang berbeda-beda, dimana untuk persamaan pertama, diperoleh nilai eigen positif yang berbeda sehingga mengakibatkan persamaan tersebut tidak stabil dengan tipe titik kritis node. Dari potret fase yang dihasilkan ditunjukkan bahwa titiktitik uji menjauh dari titik kritis. Untuk persamaan kedua, salah satu nilai eigen yang dihasilkan positif dan satunya lagi negatif, sehingga mengakibatkan persamaan tersebut tidak stabil dengan tipe titik kritisnya adalah pelana atau saddle point. Dari potret fase yang dihasilkan ditunjukkan bahwa titik-titik uji menjauh dari titik kritis dan membentuk pola seperti pelana kuda. Untuk persamaan ketiga, nilai eigen yang diperoleh adalah positif sama. Sehingga mengakibatkan persamaan tersebut tidak stabil dengan tipe titik kritisnya

15


adalah node. Dari potret fase yang dihasilkan ditunjukkan bahwa titik-titik uji menjauh dari titik kritis. Untuk persamaan keempat diperoleh nilai eigen yang sama, namun bernilai negatif. Hal ini mengakibatkan persamaan tersebut stabil asimtotik dengan tipe titik kritisnya adalah node. Dari potret fase yang dihasilkan ditunjukkan bahwa titik-titik uji masuk ke titik kritis. Dari persamaan kelima diperoleh nilai eigen yang berbeda dan juga kedua nilai eigen yang dihasilkan bernilai negatif, sehingga mengakibatkan persamaan tersebut stabil asimtotik dengan tipe titik kritisnya adalah node. Dari potret fase yang dihasilkan ditunjukkan bahwa titik-titik uji masuk ke titik kritis. Persamaan keenam diperoleh nilai eigen yang imaginer, dimana 𝛼 dari nilai eigen tersebut bernilai negatif. Hal ini mengakibatkan persamaan tersebut stabil saimtotik dengan tipe titik kritisnya adalah spiral point. Dari potret fase yang dihasilkan ditunjukkan bahwa titik-titik uji masuk ke titik kritis dan berbentuk seperti spiral. Persamaan ketujuh diperoleh nilai eigen yang imaginer, dimana 𝛼 dari nilai eigen tersebut bernilai positif. Hal ini mengakibatkan persamaan tersebut tidak stabil dengan tipe titik kritisnya adalah spiral point. Dari potret fase yang dihasilkan ditunjukkan bahwa titik-titik uji menjauh dari titik kritis dan berbentuk seperti spiral.

16


Persamaan kedelapan diperoleh nilai eigen yang imaginer, dimana 𝛼 dari nilai eigen tersebut bernilai 0. Hal ini mengakibatkan persamaan tersebut stabil dengan tipe titik kritisnya adalah center. Dari potret fase yang dihasilkan ditunjukkan bahwa titiktitik uji masuk ke titik kritis dan membentuk pola center.

17


BAB IV KESIMPULAN DAN SARAN 4.1

Kesimpulan Berikut beberapa jenis nilai eigen yang diperoleh dan kestabilan yang dihasilkan.

1. Nilai eigen positif berbeda menyebabkan persamaan tidak stabil dengan tipe titik kritisnya node. 2. Nilai eigen positif negatif menyebabkan persamaan tidak stabil dengan tipe titik kritisnya titik pelana (saddle point). 3. Nilai eigen positif sama menyebabkan persamaan tidak stabil dengan tipe titik kritisnya node. 4. Nilai eigen negatif sama menyebabkan persamaan stabil asimtotik dengan tipe titik kritisnya node. 5. Nilai eigen negatif berbeda menyebabkan persamaan stabil asimtotik dengan tipe titik kritisnya node. 6. Nilai eigen imaginer dengan 𝛼 negatif menyebabkan persamaan stabil asimtotik dengan tipe titik kritisnya spiral point. 7. Nilai eigen imaginer dengan 𝛼 positif menyebabkan persamaan tidak stabil dengan tipe titik kritisnya spiral point. 8. Nilai eigen imaginer dengan 𝛼 = 0 menyebabkan persamaan stabil dengan tipe titik kritisnya center.

18


4.2

Saran Untuk proses pencarian nilai eigen dan khususnya penggambaran potret fase

akan lebih efektif jika menggunakan program Maple 13 dari pada menggunakan cara manual. Sehingga sangatlah penting untuk setiap mahasiswa menguasai program tersebut. Dalam penyusunan tulisan ini, penulis menyadari masih terdapat banyak kekurangan, sehingga penulis mengharapkan kritik dan saran yang membangun untuk memperbaiki kinerja penulis di kemudian hari.

19


DAFTAR PUSTAKA Anton, H. 1988. Aljabar Linier Elementer. (Alih bahasa: Pantur Silaban, dan I. Nyoman Susila. Jakarta: Penerbit Erlangga. Tersedia http://eprints.ung.ac.id/5536/5/2013-1-84202-411409045-bab201082013012017.pdf. (diakses tanggal 29 November 2015)

pada

G.J. Olsder & J.W.Van der Woude. 2004. Mathematical Systems Theory. VSSD. Tersedia pada http://eprints.ung.ac.id/5536/5/2013-1-84202-411409045-bab201082013012017.pdf. (diakses tanggal 29 November 2015) Kustnetsov, Y.A. 1998. Element of Applied Bifurcation Theory. Springer-Verlag: New York. Tersedia pada http://eprints.ung.ac.id/5536/5/2013-1-84202411409045-bab2-01082013012017.pdf. (diakses tanggal 29 November 2015) Perko, L. 2000. Differential Equations and Dynamical System. Springer-Verlag: New York. Tersedia pada http://eprints.ung.ac.id/5536/5/2013-1-84202411409045-bab2-01082013012017.pdf. (diakses tanggal 29 November 2015) Wiggins, S. 1990. Introduction to Applied Nonlinier Dynamic System and Chaos. Springer-Verlag: New York. Tersedia pada http://eprints.ung.ac.id/5536/5/20131-84202-411409045-bab2-01082013012017.pdf. (diakses tanggal 29 November 2015) https://id.wikipedia.org/wiki/Teori_sistem_dinamik

20

Created By Aristastory.Wordpress.com BAB I PENDAHULUAN. Teori sistem dinamik adalah bidang matematika terapan yang digunakan untuk

Recommend Documents