Suatu matriks disebut matriks Monge jika entri-entrinya memenuhi: ... Kata Kunci:
Aljabar max- plus, Matriks invers Monge, Nilai eigen, Vektor eigen. 1.
Seminar Sains dan Penididkan Sains VI UKSW Salatiga, 11 Juni 2011 _____________________________________________________________________________________
MASALAH VEKTOR EIGEN MATRIKS INVERS MONGE DI ALJABAR MAX-PLUS Farida Suwaibah1, Subiono2, Mahmud Yunus3 Jurusan Matematika FMIPA Institut Teknologi Sepuluh Nopember Surabaya1,2,3 e-mail:1
[email protected] ,
ABSTRAK Pada penelitian ini akan dibahas masalah nilai eigen dan vektor eigen dari matriks invers Monge di Aljabar Max-Plus. Suatu matriks disebut matriks Monge jika entri-entrinya memenuhi: , untuk semua dan . Sedangkan matriks invers Monge atau sering juga disebut sebagai matriks anti Monge adalah matriks yang entri-entrinya memenuhi: , untuk semua dan . Dalam paper ini dibahas tentang path monoton. Dengan menggunakan sifat path yang monoton selanjutnya ditentukan matriks yang akan berguna dalam penentuan vektor eigen matriks invers Monge. Pada akhir pembahasan akan diterapkan satu Algoritma untuk mendapatkan nilai eigen dan vektor eigen dari matriks invers Monge Kata Kunci: Aljabar max- plus, Matriks invers Monge, Nilai eigen, Vektor eigen.
1. PENDAHULUAN Pada beberapa permasalahan, matriks digunakan untuk memodelkan suatu sistem dan sistem tersebut di selesaikan sehingga didapatkan solusinya. Dari sebuah matriks dapat ditentukan nilai eigen dan vektor eigen yang mana keduanya memainkan peranan yang sangat penting dalam suatu sistem. Peranan nilai eigen dan vektor eigen di Aljabar Max-Plus diantaranya menggambarkan kestabilan sistem dan menganalisis kedinamikan suatu sistem. Sebagai contoh adalah sistem produksi sederhana yang bisa dilihat di Subiono (2003) dan analisis kedinamikan sistem pada masalah penjadwalan flow shop di Nur Shofianah dan Subiono (2008). Seperti pada matriks biasa, matriks invers Monge juga mempunyai nilai eigen dan vektor eigen. Banyak permasalahan yang diakui dapat diselesaikan dengan mudah jika melibatkan matriks Monge atau matriks invers Monge. (Aleksey A. Imaev, Robert P.Judd, 2010). Sebagai contoh pada klas travelling salesmen problem (TSP) dengan cost matrix merupakan matriks invers Monge maka akan didapatkan rute terpendek dengan pola tertentu. Sebuah problem spectral pada matriks Monge dan matriks invers Monge di Aljabar Max-Plus telah dipelajari di Gavalec dan Plavka. Di Gavalec dan J.Plavka (2003) ditunjukkan bahwa nilai eigen dari suatu matriks invers Monge yang diberikan dapat ditentukan dengan mengambil nilai maximum dari elemen-elemen pada diagonal utama. Sedangkan di Aleksey, A. Imaev (2009) ditunjukkan peran path monoton dalam menentukan yang digunakan untuk mendapatkan vektor eigen di Aljabar Max-Plus. Selanjutnya diterapkan Algoritma 1, yang dipakai untuk menentukan vektor eigen matriks invers Monge di Aljabar Max-Plus. Pada penelitian ini akan dibahas mengenai permasalahan yang diberikan oleh Aleksey A. Imaev dan Robert P. Judd yaitu : bagaimana menentukan vektor eigen matriks invers Monge di Aljabar Max-Plus dengan menerapkan satu algoritma.
2. ALJABAR MAX-PLUS. Aljabar max-plus merupakan contoh struktur aljabar yang disebut dioid (semiring idempotent). Aljabar max-plus adalah Himpunan dengan operasi dan yang , dimana dengan dinyatakan dengan dan . Adapun struktur Aljabar dari dijelaskan dalam definisi berikut: Definisi 1. Struktur Aljabar
(Bacelli, dkk, 1992)
Simbol menyatakan himpunan yang mempunyai dua operasi biner yaitu maksi mum yang dinotasikan dengan dan penjumlahan yang dinotasikan dengan . Sedemikian hingga untuk setiap
berlaku :
dan Pada operasi mempunyai elemen netral , karena untuk setiap mempunyai elemen satuan , karena . operasi
berlaku:
dan
Himpunan semua matriks ukuran di Aljabar Max-Plus dinyatakan dengan . Dengan cara yang sama pada aljabar biasa, elemen pada matriks baris ke kolom ke ditulis “ ”, dimana dan . Sebagaimana biasa Matriks dapat ditulis sebagai:
Operasi pada matriks dijelaskan berikut ini:
adalah penjumlahan dan perkalian, seperti yang
Penjumlahan matriks oleh Matriks
dinotasikan dengan dan didefinisikan sebagai: , dengan dan . Adapun perkalian Skalar dengan matriks , didefinisikan = . Sedangkan Perkalian Matriks dengan , didefinisikan sebagai:
,
dan
.
Suatu matriks dapat direpresentasikan ke dalam bentuk graf dan sebaliknya dari suatu graf dapat dinyatakan ke dalam bentuk suatu matriks. dimana elemen-elemen matriks merupakan bobot busur dari graph, seperti disebutkan dalam defnisi berikut: Definisi 2. Communication Graph (Aleksey Imaev, 2009) Untuk sebarang matriks dapat dikaitkan dengan sebuah graph ( dinotasikan dengan dan disebut dengan “ communication graph” pada mempunyai node dengan . Untuk sebarang dua node sebuah busur dengan bobot jika hanya jika
Page 2
.
yang ada
Communication graph dinotasikan dengan Sebuah path
di
dapat dinyatakan sebagai urutan dari busur-busurnya yaitu : atau sebagai urutan nodenya yaitu:
Suatu path dikatakan elementer jika tidak ada node terjadi dua kali. Sedangkan suatu sirkuit adalah suatu path dengan atau disebut path elementer tertutup. Dan suatu sirkuit yang mempunyai bobot rata-rata maksimum disebut sirkuit kritis. Pengertian panjang, bobot dan bobot rata-rata untuk suatu sirkuit sama seperti path. (Subiono, 2010) Panjang suatu path adalah banyaknya busur dalam suatu path yang dinyatakan dengan , , yaitu jumlah bobot dari seluruh sedangkan bobot dari suatu path dinyatakan dengan busur pada suatu path. Sedangkan bobot rata-rata dari path adalah . Suatu graph dikatakan strongly connected jika untuk sebarang dua node ada suatu path dari ke . Dan suatu matrik dikatakan irreducible jika communication graphnya adalah strongly connected. Selanjutnya nilai eigen dan vektor eigen dalam Aljabar Max-Plus diberikan melalui definisi sebagai berikut: Definisi 3. ( Aleksey, A.Imaev, 2009) Penentuan suatu nilai dan suatu vektor satu elemen berhingga sedemikian hingga memenuhi dari . Maka disebut nilai eigen dan dengan nilai eigen .
yang memuat paling sedikit disebut suatu eigenproblem adalah vektor eigen yang bersesuaian
Interpretasi grafis nilai eigen dari matriks irreducible adalah sirkuit elemnter mempunyai bobot rata-rata maksimum, atau dengan kata lain:
di
,
dimana
(1)
merupakan himpunan sirkuit elementer di
Selanjutnya, untuk sebarang
A
...
A
.
didefinisikan beberapa hal sebagai berikut:
A
k
A
yang
( k
1)
,
k 1
Page 3
...
k
1
A
k
3. MATRIKS INVERS MONGE Definisi 4. (Aleksey, A.Imaev, 2009) Suatu matriks
berukuran
disebut matriks Monge jika dan hanya jika memenuhi: , untuk
(2)
dan disebut matriks invers Monge jika memenuhi: untuk Di Aljabar Max-Plus, suatu matriks memenuhi:
(3)
disebut matriks
invers Monge jika entri-entrinya
Untuk mengetahui suatu matriks ukuran memenuhi matriks invers Monge atau yang berukuran matriks Monge, yaitu dengan memeriksa disetiap submatrik apakah memenuhi pertidaksamaan (2) atau pertidaksamaan (3) (Gerhard J. Woeginger, 2006). Pertidaksamaan (2) atau (3) berlaku untuk baris dan kolom yang berdekatan. Contoh 1: Matriks
dengan ukuran
Ada 9 submatriks ukuran ( tersebut memenuhi (3), yaitu:
, yaitu:
yang dapat diperiksa apakah masing-masing submatriks , ,
3. PATH MONOTON Pada bagian ini akan dijelaskan pengertian suatu path monoton melalui satu definisi, dan beberapa contoh. Dari pengertian path yang monoton dapat di gunakan untuk menentukan “ ” dengan merupakan matriks invers Monge . Definisi 4.(Aleksey. A.Imaev, 2009) Diberikan matriks misalkan suatu path
dalam
dengan graf dari matriks yaitu: ,
Titik
disebut sebagai puncak, jika dan
Sedangkan titik
. disebut sebagai lembah, jika
dan Page 4
adalah
, selanjutnya di
Definisi 5. (Aleksey. A.Imaev, 2009, hal 97) Suatu path elementer atau sirkuit elementer dalam graf path tersebut tidak terdapat puncak atau lembah.
dikatakan monoton, jika dalam
Dari definisi di atas, berikutnya diturunkan satu Teorema untuk mendapatkan matriks Teorema 1. (Aleksey, A.Imaev, 2009) Misalkan
merupakan matriks invers Monge dengan
Dimana
, maka:
menyatakan himpunan semua path mono di
dari titik ke
Contoh 2: Diberikan matriks
Matriks
yang irreducible. Akan ditentukan matriks
dengan memperhatikan path yang monoton di
Pertama, merepresentasi matriks
ke graf
merupakan matriks invers Monge
sebagai berikut: 3
1
8
2 2
5 4
2
1 3
7
Gambar 1. Graf Seperti yang telah disebutkan pada bagian sebelumnya, bahwa eigenvalue untuk suatu matriks sebarang di Aljabar Max-Plus diperoleh dengan menentukan bobot rata-rata maksimum dari sebarang sirkuit di graf Adapun bobot rata-rata dari sebarang sirkuit di graf
diberikan pada tabel berikut ini:
Tabel (1). Tabel sirkuit elementer dari Gambar 1. Cirkuit
Panjang 1 1 1 2 2 3 3
Bobot 8 2 7 8 5 9 9
Bobot rata-rata 8 2 7 4 3
Berdasarkan Tabel (1) di atas, nilai maksimum bobot rata-rata adalah 8, maka eigenvalue dari adalah 8. dan didapatkan yaitu: Page 5
Selanjutnya untuk memperoleh dapat dilakukan dengan memperhatikan path mono dari Terlebih dahulu dibuat graf representasi dari yang ditunjukkan pada gambar berikut: -5 0
1
-6 2
-3 -6
-4 -7 3
-1
Gambar 2. Graf Adapun path mono dari graf matriks
diberikan pada tabel berikut:
Tabel (2). Tabel path monoton dari Gambar 2. Path Mono dari node ke node
Panjang path mono 1 1 1 2
Bobot path mono 0
Bobot maksimum node ke node 0
1 1 1 1 2 1 1
Berdasarkan Tabel (2) diatas, diperoleh matriks
Page 6
dari
4. Eigenvalue dan Eigenvektor Matriks Invers Monge Matriks invers Monge seperti halnya matriks biasa juga mempunyai nilai eigen dan vektor eigen. Misalkan dan memenuhi matriks invers Monge maka nilai eigen dari yang dinotasikan dengan merupakan nilai maksimum dari elemen-elemen pada baris ke kolom ke , sebagaimna dijelaskan dalam teorema berikut: Teorema(2): (Martin Gavalec, dkk, 2003). Jika
max
memenuhi matriks invers Monge maka:
i N
berikut:
Sedangkan eigenvector dari matriks invers Monge diberikan melalui definisi
Definisi 6. (Aleksey, A Imaev, 2009, hal 101) n
Misalkan
dimana
adalah eigenvalue dari
dan
i 1
misalkan merupakan bilangan bulat sedemikian hingga dari adalah eigenvector matriks dan , dengan kata lain:
maka vektor kolom ke –
Dari definisi diatas, eigenvector dari suatu matriks dapat ditentukan dengan memperhatikan kolom ke – dari matriks yang memenuhi . Berikut ini diberikan satu teorema yang menjelaskan tentang suatu matriks invers Monge yang irreducible dengan . Teorema 3. (Aleksey, A.Imaev, 2009) Misalkan Matriks irreducible Maka:
dengan
, untuk
dan
memenuhi matriks invers Monge,
,
i 1
j 1
k
untuk
j 1
k
j 1
untuk
Contoh 3: Dari Contoh 2 telah didapatkan matriks Definisi 6 maka vektor eigen dari matriks
adalah:
. Berdasarkan .
Selanjutnya, dalam menentukan vektor eigen matriks invers Monge di Aljabar Max-Plus diberikan Algoritma sebagai berikut:
Page 7
Algoritma 1 Input :
suatu matriks invers Monge irreducible
Output: Eigenvalue dan eigenvector dari n
1: i 1
2: 3: Tentukan , sedemikian hingga 4: 5: 6: 7: for
to
do
l 1
8: k c 1
9: end for 10: for
to 1 do c 1
11: k l 1
12: end for Contoh 4. Misalkan matriks invers Monge pada Contoh 3 di atas, akan ditentukan eigenvalue dan eigenvector dari matriks . Karena matriks memenuhi matriks invers Monge maka dengan menggunakan Teorema (2) diperoleh eigenvalue dari adalah : 3
i 1
Dari Contoh 3 diperoleh matriks
, sehingga
adalah 1. Selanjutnya ditentukan eigenvector dari Algoritma 1 di atas, melalui tahapan-tahapan sebagai berikut: Untuk 1: 2:
, , maka diperoleh
,
, maka diperoleh
Untuk
Page 8
yang memenuhi dengan menerapkan
l 1
3:
, diperoleh: k c 1 2
k 2
Sehingga eigenvector dari
adalah:
. Hasil perhitungan dari menerapkan Algoritma (1) di atas
apabila di perhatikan menyatakan kolom ke-1 dari matriks
5. Kesimpulan 1. Vektor eigen di Aljabar Max-Plus memainkan peran penting dalam mempelajari perilaku steady state dari Sistem Event Diskret Deterministic Cyclic. 2. Path monoton memainkan peran dalam menentukan vektor eigen matriks invers Monge di Aljabar Max-Plus sehingga: , dimana yang memenuhi vektor eigen digunakan Algoritma 1.
. Selanjutnya untuk mendapatkan
6. DAFTAR PUSTAKA [1] Baccelli, F., Cohen, G., Olsder, G.J. dan Quadrat, J.P. (1992), Synchronization and Linearity An Algebra for Discrete Event Systems, John Wiley & Sons, New York. [2] Gavalec, M.,Plavka J, An algorithm for maximum cycle mean of Monge matrices in Max-Algebra, Discrete Applied Mathematics 127(2003) 651-656 [3] Imaev, A.A. (2009). Hierarchical Modeling of Manufacturing Systems Using Max-Plus Algebra, Disertasi P.Hd,(114pp), Ohio University, Ohio. [4] Imaev, A.A., Judd, P.R, (2010), Computing an eigenvector of an inverse Monge matrix in Max-Plus Algebra, Discrete Applied Mathematics 158 (2010) Ohio university, Athens, OH 45701, United States. (1701-1707). [5] Subiono, (2003), Terapan Aljabar Max-Plus pada Sistem Produksi Sederhana serta Simulasinya dengan Menggunakan Matlab, Jurusan Matematika FMIPA ITS Surabaya [6] Subiono, (2010), Aljabar Max-Plus Dan Terapannya, versi 1.00, Jurusan Matematika FMIPA ITS, Surabaya [7] Woeginger, G.J, (2006), Some Problem Around Travelling Salesman Problem, The Algorithmics Column by Gerhard J Woeginger, Departement of Mathematematics and Computer Science Eindhoven University of Technology
Page 9
