KETERKENDALIAN SISTEM BANDUL GANDA

J. Sains MIPA, Agustus 2008, Vol. 14, No. 2, Hal.: 90 - 100 ISSN 1978-1873

KETERKENDALIAN SISTEM BANDUL GANDA Amanto* dan La Zakaria Jurusan Matematika FMIPA Universitas Lampung korespondensi e-mail: [email protected]

*Alamat

Diterima 31 Juli 2008, disetujui untuk diterbitkan 19 September 2008

ABSTRACT Simple double pendulum system is a system consisting of two objects and two pieces of strong strand of metal but light. The system sway vertical with influenced by earth gravitation. There are three main consepts in mathematical system theory, which are : stability, controllability, and observability. In this paper we will focus on controllability. Given time-invariant differential systems x& = Ax + Bu , y = Cx + Du; with

x ∈ R n , u ∈ R m and y ∈ R p and will be indicated the system by (A,B,C,D). Controllability is completely determined by the matrices A and B and therefore we will discuss the controllability of the pair (B,A). The system (A,B,C,D) is controllable if and only if the controllability matrix, R = [B AB A 2 B ... A n −1 B ] has rank = n. In the case of a single output system, R is square matrix, then (B,A) is controllable if R is non singular. Moreover, we will investigate controllability of simple double pendulum system. The result, the system is controllable when the system is given by single and double input functions where m1 ≠ m2 . Keywords: controllability, simple double pendulum, controllability matrix, non singular matrix

1. PENDAHULUAN Suatu fenomena alam dapat dinyatakan dalam suatu model matematis, yang berupa sistem persamaan diferensial. Teori sistem merupakan cabang ilmu matematika yang berfokus pada kajian kontrol input dan output. Teori ini telah banyak diterapkan pada bidang teknik, fisika, biologi, dan bidang ilmu lain dan telah dikembangkan secara intensif oleh para ahli teknik teoritis, terutama teori sistem linear atas bilangan real 1). Kestabilan merupakan konsep dasar yang penting dalam teori sistem matematis, disamping konsep keterkendalian dan keterobservasian2) . Model-model dalam aplikasi matematika pada ilmu rekayasa, fisika, biologi dan lain sebagainya, banyak yang berbentuk sistem persamaan diferensial biasa. Perilaku solusi dari sistem seperti itu adalah contoh permasalahan yang dibahas dalam sistem dinamik. Sistem kerja pada Tim SAR (Search and Rescue) dan suplai makanan atau amunisi ke barak menggunakan pesawat helikopter merupakan dua contoh sistem kerja yang identik dengan sistem kerja pendulum. Dari contoh fenomena terakhir di atas, cara kerja permainan sirkus dan helikopter merupakan contoh sistem kerja yang dapat dikaitkan dengan sistem kerja bandul. Sistem kerja bandul tersebutlah yang lalu akan diubah bentuk ke dalam suatu persamaan diferensial. Sistem persamaan diferensial dari sistem bandul ganda dikategorikan ke dalam sistem Hamiltonian, karena sistem kerja bandul tersebut terkait dengan konservasi energi. Ia bersifat konservatif karena sistem tidak mengalami kehilangan energi seiring berjalannya waktu. Dalam pengilustrasiannya, sistem bandul ganda dapat kita gambarkan sebagai sebuah ayunan yang dipasangkan pada ayunan yang lainnya. Dengan demikian, sistem ini melibatkan 4 variabel, yaitu dua variabel yang melibatkan panjang dan dua variabel lainnya yang melibatkan massa. Secara teori fisika, dalam sistem ini berlaku hukum kekekalan energi. Oleh karena itu, bandul ganda memiliki sifat kualitatif. Sifat kualitatif pada Sistem Persamaan Diferensial Biasa perlu dipertahankan pada hasil integrasi numerik, karena sifat tersebut sangat berperan penting dalam menganalisa prilaku solusi sistem tersebut. Pembahasan terhadap sifat kualitatif sistem meliputi beberapa topik kajian, yaitu kestabilan (stability), keterkendalian (controllability), keterobservasian (observability), dan optimisasi (optimization). Analisis sifat ini perlu dilakukan karena menyangkut keberadaan sistem untuk suatu waktu tertentu. Misalnya pada analisis

90

 2008 FMIPA Universitas Lampung

J. Sains MIPA, Agustus 2008, Vol. 14, No. 2

keterkendalian sistem, suatu sistem dikatakan terkendali apabila pada keadaan tertentu, katakan x1 untuk suatu waktu tertentu (t >0) yang dapat dimulai dari keadaan awal x0 dengan cara mengatur fungsi masukan, pembuat model dapat menginginkan keadaan tertentu pula3) . Berdasar latar belakang tersebut di atas permasalahan yang akan diselesaikan dalam penelitian ini adalah mengkaji sifat kualitatif dari sistem Hamiltonian pada bandul ganda. Solusi sistem bandul ganda secara teori memiliki sifat kualitatif, karena bandul melibatkan hukum kekekalan energi. Karenanya proses penyelesaian sistem bandul ganda harus mempertahankan sifat kualitatifnya. Dua sifat yang dapat dicermati langsung dalam bandul ganda adalah Kestabilan dan Keterkendalian. Masalah solusi eksak dan kestabilan telah dibahas4) . Penyelidikan keterkendalian sistem memiliki prosedur yang berbeda dengan kestabilan sistem. Keterkendalian memerlukan fungsi masukan agar sistem tersebut dapat diketahui keterkendaliannya. Pada penelitian ini penulis membahas syarat dan sifat keterkendalian sistem bandul ganda. Sistem dinamik yang cukup realistis adalah sistem bandul ganda sederhana. Sistem bandul ini dapat merepresentasikan sistem kerja tim SAR (Search and Rescue) yang sedang mengupayakan penyelamatan dengan membawa perlengkapan penyelamatan atau korban bencana alam atau kecelakaan dengan menggunakan helikopter. Sistem bandul ganda sederhana adalah sistem yang terdiri dari dua benda B1 dan B2 dengan massa masing-masing benda adalah m1 dan m2 . Selain itu benda tersebut masing-masing dihubungkan dengan dua helai kawat yang kuat tapi ringan L1 dan L2 dengan panjang masing-masing kawat adalah l1 dan l2 . Benda B1 terpasang pada ujung kawat L1 (ujung kawat L1 lainnya terpasang mantap pada sebuah bidang). Sementara itu benda B2 terpasang pada ujung kawat L2, di bawah pengaruh gravitasi 5) (ujung kawat L2 lainnya terpasang mantap pada benda pertama) (lihat Gambar 1).

θ1

L1

l = l1 + l2

l1

B1

m1

θ2

L2

B2

m1

l2

m2 m2

m2

Gambar 1. Sistem Bandul Ganda Sederhana Dari Gambar 1 sistem bandul ganda memiliki 4 (empat) parameter yakni l1, l2, m1, dan m2. Dengan Dipengaruhi oleh grafitasi bandul gandaberosilasi pada bidang vertikal dengan sudut perpindahan untuk suatu waktu adalah θ1 ( t ) dan θ 2 ( t ) . Model matematika untuk sistem bandul ganda yang diilustrasikan dalam Gambar 1. adalah sebagai berikut 6) :


91

Amanto dan La Zakaria… Keterkendalian Sistem Bandul Ganda

2 ″ ″ ″ 2 ″ ( m1 + m 2 ) l1 θ 1 + m 2 l1l 2θ 2 + ( m1 + m 2 )l1 gθ 1 = 0 , m 2 l1l 2θ 1 + m 2 l 2 θ 2 + m 2 l 2 gθ 2 = 0

(1)

dengan g adalah kecepatan gravitasi. Kajian tentang keterkendalian penting artinya, sebab kita akan tahu bahwa keadaan tertentu dapat dicapai dengan memulai dari keadaan awal dan mengatur fungsi masukan. 1.1. Keterkendalian Sistem Persamaan Diferensial Linier Keterkendalian (controllability) merupakan konsep dasar dalam teori sistem matematis. Konsep ini memegang peranan penting dalam desain dan kontrol suatu sistem. Misalkan diberikan sistem persamaan diferensial time-invariant:

x& = Ax + Bu

(2)

y = Cx + Du

dengan x,u dan y adalah fungsi vektor-fungsi vektor yang bernilai real, untuk fungsi y=fungsi keluaran (output) yang tergantung pada x ∈ R n , u ∈ R m dan y ∈ R p , u=fungsi masukan (input). Dalam sistem Persamaan (2) A,B,C,D masing-masing adalah matriks berukuran n × n, n × m, p × n, p × m atas lapangan bilangan R. Sistem Persamaan (2) disebut sistem linier kontinu berdimensi n dengan m-masukan dan p-keluaran yang invariant terhadap waktu (elemen-elemen dari matriks A,B,C dan D bebas waktu t). Sistem Persamaan (2) dengan syarat awal x (0) = x 0 yang kemudian ditulis dengan x (t , x 0 , u ) dan fungsi keluarannya y (t , x 0 , u ) . t

x(t , x0 , u ) = e At x0 + ∫ e A (t − s ) Bu ( s ) ds

(3)

0

t

y (t , x0 , u ) = Ce x 0 + ∫ Ce A (t − s ) Bu ( s ) ds + Du (t ) At

(4)

0

1.2. Definisi3) (Kerkendalian) Sistem (A,B,C,D) dikatakan terkendali (controllable) jika untuk setiap vektor awal x0 ∈ R n dan untuk vektor lain x1 ∈ R n , terdapat t1 > 0 dan u (t ) ∈ R n sehingga berlaku x(t , x 0 , u ) = x1 . Keadaan tertentu x1 akan tercapai dengan memasukkan syarat awal x0 . Dari definisiketerkendalian di atas dapat diinterpretasikan bahwa suatu sistem dikatakan terkendali jika untuk sebarang titik (keadaan) x1 dapat dicapai (reachable) dengan memulai dari titik (keadaan) x0 dan mengatur fungsi masukan u. Pada prinsipnya, keterkendalian sistem Persamaan (4) dapat dikaji dari apa yang dinamakan matriks keterkendalian (controllability matrix).

ℜ = [ B | AB | A 2 B | K | A n −1 B ]

(5)

yang merupakan matriks berukuran n × nm yang memuat n blok AjB,j=0,1,…,n-1. Image dari ℜ dinotasikan dengan Im( ℜ ) adalah himpunan bagian dari Rn yang direntang oleh vektor kolom dari ℜ dan dinamakan ruang bagian keterkendalian (controllability subspace). 1.3. Uji Keterkendalian Untuk memeriksa keterkendalian sistem linear dapat digunakan definisi keterkendalian, tetapi hal ini rumit, karena harus mencari penyelesaian dalam bentuk integral yang kadang-kadang dihadapkan pada bentuk integral yang tidak sederhana. Untuk mengatasai masalah ini digunakan teorema uji keterkendalian

92



yaitu dengan cukup memeriksa matriks keterkendaliannya. Cara ini lebih mudah dan merupakan metode modern dalam analisis keterkendalian sistem. 1.4. Teorema3) (Uji Keterkendalian) Sistem ( A, B, C , D ) terkendali jika dan hanya jika rank (ℜ) = n 3). Diketahui bahwa rank ( ℜ ) = n jika dan hanya jika det ( ℜ ) ≠ 0. Dengan demikian teorema di atas dapat digunakan untuk memeriksa keterkendalian suatu sistem.

2. METODE PENELITIAN 2.1. Bahan dan Alat Bahan penelitian yang digunakan adalah buku-buku dan jurnal matematika dan analisis komputasi. Sedangkan alat yang digunakan adalah seperangkat komputer PC Pentium IV yang berisikan software Mathematica 4.0.

2.2. Prosedur Penelitian Metode yang digunakan dalam penelitian ini adalah studi pustaka dan praktik komputasi. Langkah-langkah yang di ambil dalam penelitian ini adalah: 1. Mengkaji model matematika sistem persamaan diferensial dari fenomena pendulum ganda (persamaan gerakan). 2. Mengubah model matematika pada gerak pendulum ganda ke dalam sistem persamaan diferensial biasa yang terdiri dari 4 variabel, yaitu m1, m2, l1, dan l2 dimana m adalah massa benda dan l adalah panjang tali pendulum. 3. Mencari matriks keterkendalian dari sistem persamaan diferensial biasa linear yang telah disusun pada point (b), yaitu : ℜ = [ B | AB | K | A n−1 B ] . 4. Memeriksa keterkendalian pada matriks pada (c) dengan menggunakan rank matriks atau determinan matriks. 5. Melakukan analisis atas keterkendalian sistem pada (d).

3. HASIL DAN PEMBAHASAN Pembahasan mengenai persamaan gerak dari sistem Hamiltonian bandul ganda akan dilakukan. Tinjau persamaan gerak sistem bandul ganda pada sistem ( 1) dan dimisalkan: 2 a = (m1 + m2 )l1 g , b = (m1 + m2 )l1 (6)

p = m 2 l1l 2 , q = m2 l 2 g , r = m 2 l 2

2

dengan a, b, p, q dan r > 0 .

Setelah dimisalkan, sistem (1) menjadi:

bθ&&1 + pθ&&2 = − aθ1 pθ&&1 + rθ&&2 = − qθ 2

(7)

Untuk lebih menyederhanakan sistem, kembali dimisalkan:

θ1 = x1 ;θ&1 = x 2 ;θ&& = x& 2 ; θ 2 = x3 ;θ&2 = x 4 ;θ&&2 = x& 4

(8)

Sehingga Persamaan (7) menjadi:

bx& 2 + px& 4 = − ax1 px& 2 + rx& 4 = − qx3

(9)

Dari sistem yang sudah disederhanakan di atas, diperoleh sistem persamaan diferensial sebagai berikut:


93


x&1 = θ1 ' = x2 − ar pq − ar pq θ1 + 2 θ2 = x1 + x3 2 2 2 − p + br − p + br − p + br − p + br x&3 = θ 2 ' = x4 x&2 = θ1 " =

x&4 = θ 2 " =

(10)

− ap bq − ap bq x1 + 2 x3 θ1 + 2 θ2 = 2 p − br p − br p − br p − br 2

Dalam bentuk matriks, sistem Persamaan (9) menjadi:  − ar  T T T 0  x1  0  x1   − p 2 + br   x1  1  x  0   x    x  0   2  , x&1 =    2  , x&1 =    2  , x& 2 =  pq   x  0   x 3  0   x 3   3 2           − p + br 0   x 4  1  x 4  x4    0   atau

 − ap  T  p 2 − br   x1     0  x2  x& 2 =  bq    x3   p 2 − br   x   0   4  

x& = Ax + Bu

 x&1   x&   2 =  x& 3     x& 4 

 x1  x  A 2  + Bu  x3     x4 

(11)

Pada sistem gerak bandul ganda akan diberikan masukan pada tiap variabel x.

0   − ar  2  − p + br x& =  0   − ap   p 2 − br

1

0

pq 2 − p + br 0 bq 2 p − br

0 0 0

0  0   1  0 

 x1  x   2  + Bu  x3     x4 

(12)

Dengan

0    x&1   −2ar  x&   − p + br x& =  2  , A =   x&3  0     −ap  x&4   2  p − br

1 0 0 0

0 pq 2 − p + br 0 bq 2 p − br

0  0  , x= 1   0 

 x1  x   2  , dan B adalah matriks dari fungsi  x3     x4 

masukan (u) yang akan diberikan terhadap variabel x. Berikut ini akan ditinjau 2 kasus fungsi masukan, yaitu diberikan satu dan dua fungsi masukan (u) pada sistem gerak bandul ganda Kasus 1: a. Diberikan satu masukan pada variable x1.

94



1  0 B = ( u1 ) =   0   0

0   −ar  2  − p + br AB =  0   −ap  2  p − br  − ar  − p 2 + br   0  2  A B=  −ap  p 2 − br   0  

1 0 0 0

0 pq 2 − p + br 0 bq 2 p − br

0

pq − p + br

−ar − p 2 + br

0

0

bq p 2 − br

−ap p 2 − br

0

2

0  0  1   0 

0    1  −ar   0  2   =  − p + br   0   0       − ap 0    2   p − br 

   pq  2 − p + br    1   bq  p 2 − br  0

 −ar  2 1   − p + br    0  0    =  0  − ap     2 0   p − br    0  

−ar pq   0 0  −p2 + br −p2 + br     a2r2  −ap2q −arpq pbq2 0 0 + +  2  2 2 2 2 2 2 2 ( − p + br ) ( p − br )( − p + br ) ( − p + br ) ( − p + br )( p − br )  A3B =    −ap 0 0 0   2 p − br     a2 pr −apbq −ap2q b2q2 + 2 0 + 2 0   2 2 2 2 2 2 ( p −br)(−p + br) ( p − br)  (−p + br)( p − br) ( p − br) 

1 0  0  0

0     2 2 2 −ap q  ar  + (−p2 + br)2 ( p2 − br)(−p2 + br)  =  0    a2 pr −apbq  + 2   2 2 2  (−p + br)( p − br) ( p − br) 

Matriks keterkendaliannya:

ℜ1 = [ B | AB | A 2 B | A 3 B] Berdasarkan teorema uji keterkendalian, untuk mengetahui suatu sistem linier terkendali atau tidak adalah apabila rank matriks keterkendaliannya sama dengan banyaknya variabel bebas dalam sistem persamaan (rank( ℜ) = n) , atau dengan kata lain determinan matriks keterkendaliannya tidak sama dengan nol.


95


Matriks keterkendalian yang pertama adalah ℜ1 yang berukuran 4 × 4 . − ar   0 0 2 1  − p + br     − ar a2r 2 − ap 2 q 0 + 2 0  2 2 2 2 − p + br ( − p + br ) ( p − br )( − p + br )   ℜ1 =   − ap 0 0 0  2 p − br    − ap a 2 pr − apbq  0 + 0  p 2 − br ( − p 2 + br )( p 2 − br ) ( p 2 − br ) 2   Determinan matriks ℜ1 yakni:

a3 p2q (− p 2 + br ) 3

; p 2 ≠ br

Karena det ≠ 0, berarti bahwa saat diberikan satu fungsi masukan terhadap variabel x1 , sistem linier untuk matriks ℜ1 terkendali. b. Diberikan satu masukan pada variabel x2 , yaitu

0 1  B = (u2 ) =   0   0 Maka sistem akan terkendali. c. Diberikan satu masukan pada variabel x3 , yaitu:

0  0  B = (u3 ) =   1    0  Maka sistem terkendali. d. Diberikan masukan pada variabel x4 , yaitu: 0  0  B = (u4 ) =   0    1  Maka sistem akan terkendali. Kasus 2 : Selanjutnya akan diberikan dua fungsi masukan (u) terhadap dua variabel dari x1 , x 2 , x3 dan x 4 secara bersamaan, matriks keterkendalian yang dapat dibentuk dari 4C2 adalah enam buah. a. Diberikan masukan pada fungsi x1 dan x2 .

96



1  1  B = (u5 ) =   0   0 Maka sistem akan terkendali jika p 2 ≠ (a + b)(q + r ) atau pada saat nilai

l1 l 2 m2 ≠ ( g + l1 )( g + l 2 )(m1 + m2 ) . b. Diberikan masukan pada fungsi x1 dan x3 , yaitu: 1  0  B = (u6 ) =   1    0  Maka sistem akan terkendali untuk (b + p ) q ≠ a ( p + r ) atau pada saat nilai g l1 l 22 m1 m2 ≠ 0 . Sistem pasti terkendali karena g l1 l 22 m1 m 2 pasti memiliki nilai. c. Diberikan masukan pada fungsi x1 dan x4 , yaitu: 1  0 B = (u 7 ) =   0   1  Maka sistem akan terkendali jika

(

) (

)

g l2 (m1 + m2 ) ⋅ (l1 − l2 ) m1 + l2 (l1 + l2 ) m2 ≠ − l1 l23 m1 m22 + g 2 l12 (m1 + m2 ) . d. Diberikan masukan pada fungsi x2 dan x3 , yaitu: 2

2

3

0  1  B = (u8 ) =   1    0  Maka sistem akan terkendali jika 2 g l1 (m1 + m2 ) (l1 − l2 ) m1 + l1 (l1 + l2 ) m2 ≠ g 2 l13 l23 m23 (m1 + m2 ) g 2 l22 m22 + l13 l2 m1 (m1 + m2 )

(

)

(

)

e. Diberikan masukan pada fungsi x 2 dan x 4 , yaitu:

0  1  B = (u 9 ) =   0    1  Sistem akan terkendali jika (b + p )q ≠ a ( p + r ) atau pada saat − g l 1 l 22 m1 m 2 ≠ 0 . Sistem pasti terkendali karena

− g l 1 l 22 m1 m2 pasti memiliki nilai. f. Diberikan masukan pada fungsi x3 dan x 4  2008 FMIPA Universitas Lampung

97


0  0  B = (u10 ) =   1    1  Untuk mengetahui keterkendalian pada sistem bandul ganda, diberikan satu contoh berikut ini: Contoh untuk Kasus 1.a. : Diketahui m1 = 937.5kg , m2 = 312.5kg , l1 = 16m, l 2 = 16m, g = 9.8 m

s2

dan

θ1 (0) = −0.5rad , θ 2 (0) = 1.0rad , θ1′ (0) = −1.0 rad sec , θ 2 ′ (0) = 2.0 rad sec . Dengan mensubstitusikan nilai-nilai tersebut di atas ke dalam Pers. (4.2), didapatkan a = 64, b = 32, p = 8, q = 16, dan r = 8. Dari nilai-nilai tersebut dapat ditentukan nilai determinan untuk

ℜ1 . Dalam bentuk matriks: T

0  x1  1   x  &x1 =    2  , 0   x 3     0   x 4 

 − 512  x T  192   1    x  &x 2 =  0   2  , 128  x    3  192   x 4   0 

T

0  x1  0   x  &x1 =    2  , 0   x 3     1  x 4 

 − 512  x T  192   1    x  &x 2 =  0   2  − 152  x    3  192   x 4   0 

Sehingga

 0  − 512  A =  192 0  512   192  1   0 ℜ1 =   0   0  dengan

1 0 0 0

0 0 128  0 192 0 1  − 512 0 192 

0

− ar − p 2 + br

− ar − p 2 + br

0

0

−ap p 2 − br

−ap p 2 − br

0

    a2r 2 − ap 2 q + 2  2 2 2 (− p + br ) ( p − br )(− p + br )   0   a 2 pr −apbq  +  (− p 2 + br )( p 2 − br ) ( p 2 − br ) 2  0

a = (m1 + m 2 )l1 g = 64 b = (m1 + m2 )l12 = 32 p = m 2 l1l 2 = 8 q = m 2 l 2 g = 16 r = m2 l 22 = 8 Substitusi nilai di atas ke dalam determinan ℜ1

a3 p 2q 643 ⋅ 82 ⋅ 16 = (− p 2 + br )3 (−82 + 32 ⋅ 8)3 = 0,333611 98



Karena determinan yang dihasilkan dari sistem ini tidak sama dengan nol, maka sistem dapat dikendalikan pada saat diberikan masukan seperti di atas. Kemudian akan diperlihatkan keterkendalian sistem bandul ganda melalui gambar. Gambar-gambar berikut merupakan visualisasi dari keterkendalian sistem bandul ganda dengan satu fungsi masukan (u1); dengan nilai τ = 0,01 ; dan banyaknya iterasi (N) adalah 100.000. Gambar berikut adalah Poincare keterkendalian sistem antara θ1 dan p1 . Gambar 2 menunjukkan

Poincare keterkendalian sistem pada saat diberikan masukan pada θ1 , nilai θ1 berada pada

− 0,3 ≤ θ1 ≤ 0,2 sedangkan p1 berada pada − 1,75 ≤ θ1 ≤ 1 . YS 4;t

0 . 0 1;

N10 0 0 0 0

0 -0.25 -0.5

p1

-0.75 -1 -1.25 -1.5 -1.75 -0.3

-0.2

-0.1

0

0.1

0.2

q1

Gambar 2. Poincare keterkendalian sistem antara θ1 dan p1 Gambar 3. berikut adalah visusalisasi keterkendalian sistem antara θ1 dan θ 2 YS 4;t

0 .0 1;

N1 0 00 00

0.4

q2

0.2

0

-0.2

-0.4

-0.3 -0.2

-0.1

0

0.1

0.2

0.3

0.4

q1

Gambar 3. Hubungan keterkendalian sistem antara θ1 dan θ 2 Gambar 3 menunjukkan keterkendalian sistem pada saat diberikan masukan pada θ1 , ternyata

sistem terkendali untuk θ1 berada pada − 0,3 ≤ θ 1 ≤ 0,4 sedangkan θ 2 bernilai − 0,4 ≤ θ 2 ≤ 0,4.


99


Gambar 4 berikut adalah visusalisasi keterkendalian sistem antara θ1 dan p1 . Gambar 4 menunjukkan keterkendalian sistem pada saat diberikan masukan pada p1 , ternyata sistem terkendali untuk p1 bernilai

− 1,75 ≤ θ1 ≤ 1 sedangkan θ1 berada pada − 0,3 ≤ θ1 ≤ 0,2 . YS 4;t

0 . 0 1;

N1 0 0 0 0 0

0.4 0.3 0.2

q1

0.1 0 -0.1 -0.2 -0.3 -1.5

-1

-0.5

0 p1

0.5

1

1.5

Gambar 4. Hubungan keterkendalian sistem antara θ1 dan p1

4. KESIMPULAN Sifat keterkendalian sistem bandul ganda dapat diketahui melalui i matriks keterkendalian dari sistem PDB linear yang terbentuk dengan menggunakan teorema uji keterkendalian sistem. Sistem bandul ganda yang dapat merepresentasikan sistem kerja Tim SAR dan suplai makanan atau amunisi ke barak menggunakan Helikopter merupakan sistem yang terkendali baik dengan satu fungsi masukan maupun dua fungsi masukan.

UCAPAN TERIMAKASIH Penelitian ini dapat dilaksananakan dengan baik berdasarkan bantuan dana dari Hibah Penelitian (Research Grant) yang diprogramkan Program Hibah Kompetisi (PHK) A2 Jurusan Matematika FMIPA Universitas Lampung. Oleh karena itu, penulis mengucapkan terima kasih yang sebesar-besarnya kepada pengelola PHK A2 Jurusan Matematika atas segala bantuan yang telah diberikan.

DAFTAR PUSTAKA 1. Amanto, 2005. Analisis Keterobservasian Sistem Gerak Kereta Sederhana Dengan Dua Pendulum Terbalik. J. Sains Tek.., 11(2) : 101 – 104. 2. Cobb, J.D. 1984. Controllability, observability and duality in singular syastem. IEEE Trans Aut. Control, AC-29 (12) : 1076 – 1082. 3. Olsder, G.J. 1994. Mathematical Systems Theory First Edition. Delftse Uitgevers Maatschappij, Delft, Netherlands. 4. Amanto dan La Zakaria, 2007. Solusi Eksak dan Kestabilan Sistem Bandul Ganda. Jurnal Sains MIPA, Edisi Khusus 14 (1) : 23 – 32. 5. Shinbrot. 1992. Chaos in a double pendulum. Am. J. Phys. 6 (60) : 491 – 499. 6. Giordano dan Weir. 1994. Differential Equations A Modeling Approach. Addison-Wesley, Publishing Company. 100


KETERKENDALIAN SISTEM BANDUL GANDA

Recommend Documents