## Analisis Faktorisasi Matriks Tak Negatif

## Abstrak
Hasil perkalian matriks tak negatif W yang berukuran n × r dengan matriks tak negatif H yang berukuran r × m akan menghasilkan sebuah matriks tak negatif V berukuran n × m , sedemikian sehingga dapat ditulis dalam bentuk persamaan matriks W × H = V . Pada tulisan ini, kondisi seperti itu akan kita balik, yaitu dari sebuah matriks tak negatif V berukuran n × m , akan dicari dua faktor matriks tak negatif yaitu matriks W yang berukuran n × r dengan matriks H yang berukuran r × m , sedemikian sehingga berlaku V = WH . Pencarian tersebut mungkin tidak tunggal ataupun mungkin tidak ada, sehingga yang terjadi adalah hanya pendekatan terhadap matris V , yang didasarkan pada pertimbangan dimensi matriks-matriks faktornya dan norm matriks V − WH adalah seminimum mungkin.
Kata kunci : Faktorisasi , matriks tak negatif , dimensi , norm .

## 1. Pendahuluan
Penelitian ini mengkaji aspek teoritis dari faktorisasi matriks tak negatif (fmtn) yang telah banyak digunakan oleh para peneliti di bidang terapan untuk menganalisis dan mereduksi dimensi gambar ataupun signal suara seperti yang telah dilakukan oleh Guillamet, Schiele dan Vitria, 2001 yang menggunakannya untuk klasifikasi gambar, Smaragdis dan Brown, 2003 [1] menggunakannya untuk transkripsi musik polyphonic, hal yang sama juga dikemukakan oleh Abdallah  dan Plumbley pada tahun 2004 [2]. Hoyer  di tahun 2002 [3] membahas pengkodean sparse tak negatif dan menghubungkannya dengan dekomposisi yang dilakukan pada sebuah matriks tak negatif.
Prinsip sederhana pada penggunaan fmtn setiap data asli suara atau gambar disusun dalam bentuk matriks V yang tak negatif berukuran n × m . Kita akan mencari dua matriks W dan H yang ranknya adalah r ≤ min{ n , m } dan hasil kali kedua matriks tersebut dapat mendekati keorisinalan matriks V . Setiap kolom dari matriks W merupakan vektor basis sedangkan matriks H menyatakan pembobotan yang diinginkan dapat mengaproksimasi kolom-kolom yang bersesuaian pada matriks V . Pencarian kedua matriks ini akan dilakukan berulang-ulang sehingga akan didapatkan hasil kali W dan H yang dapat mengaproksimasi V . Sebuah Matriks berukuran n × m disebut tak negatif jika entri-entrinya adalah bilangan riil tak negatif dan disebut positif jika entri-entrinya adalah bilangan riil positif. Adapun karakteristik yang melekat dari bentuk matriks tak negatif ini jika ukurannya n × n , adalah tracenya atau jumlah diagonal matriksnya tak negatif dan terdapat nilai eigen yang positif sedemikian hingga absolut nilai eigen yang lain tidak lebih besar dari nilai eigen tersebut, dan disebut nilai eigen maksimal [5]. Namun untuk matriks tak negatif n × m dapat dilihat pada nilai singular terbesarnya. Berikut ini untuk menentukan nilai skalar V − WH akan digunakan bentuk nilai singularnya sebagai salah satu bentuk normnya. Norm ini berguna
Staf Pengajar pada Jurusan Matematika FMIPA Universitas Hasanuddin Makassar
Andi Kresna Jaya †
†

## Jurnal Matematika, Statistika & Komputasi Edisi Khusus Juli 2007
157
Vol 13, No 2, 142-146, Januari 2017 Vol. 13, No. 2, 147-151, January 2017
sebagai fungsi tujuan agar proses pencarian berulang-ulang atau iterasinya menjadi konvergen ke faktor-faktor yang dapat mendekati keorisinalitas matriks V .

## 2. Norm-norm untuk Matriks
Misalkan A adalah matriks berukuran n × m maka yang dimaksud dengan norm matriks adalah sebuah skalar yang menyatakan ukuran dari matriks yang bergantung pada entri- entrinya. Notasinya adalah || A ||, sedemikian sehingga memenuhi sifat - || A || ≥ 0 dan akan sama dengan nol jika A = 0 .
- || cA ||=| c | || A || untuk suatu c skalar - || A + B || ≤ || A || + || B || untuk matriks B berukuran n × m
Bentuk norm matriks ini dapat diklasifikasikan atas
a. Jumlah nilai mutlak dari entri-entri pada baris yang terbesar, disebut norm- ∞ dan dinotasikan oleh
⎭ ⎬ ⎫ ⎩ ⎨ ⎧ = = ∑ = ∞ m j ij i n i a A 1 , , 3 , 2 ,1 , max L b. Jumlah nilai mutlak entri-entri pada kolom yang terbesar, disebut norm-1 dan dinotasikan oleh
⎭ ⎬ ⎫ ⎩ ⎨ ⎧ = = ∑ = n i ij j m j a A 1 1 , , 3 , 2 ,1 , max L c. Nilai singular terbesar dari A, dinotasikan oleh ( ) ( ) t t svd AA A A A max max λ λ = = dengan () . max λ adalah nilai eigen terbesar.
d. Akar dari jumlah kuadrat tiap entri matriks disebut norm frobenius dan dinotasikan oleh bentuk
2 1 2 2 ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ = ∑ ij ij a A Contoh 1: Misalkan A matriks berukuran 2×3 berikut ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = 0 1 1 1 2 1 A , Maka norm-normnya adalah { } 4 2 , 4 max = = ∞ A { } 3 1, 3 , 2 max 1 = = A ( ) ( ) 7578 . 2 6056 . 7 max max = = = = t t svd AA A A A λ λ ( ) 2 2 0 1 1 1 2 1 2 1 2 2 2 2 2 2 2 = + + + + + = A
Masih banyak bentuk norm matriks yang lain, namun yang akan kita gunakan sebagai bentuk fungsi tujuan pada fmtn hanya keempat bentuk norm di atas.
Jurnal Matematika, Statistika & Komputasi Edisi Khusus Juli 2007
Jurnal Matematika, Statistika, & Komputasi
Vol. 13  No.2 Januari  2017

## 3. Faktorisasi Matriks Tak Negatif
Untuk mendapatkan hampiran faktorisasi V ≈ WH adalah pemilihan dua matriks tak negatif sebagai faktor kemudian harus dilihat bentuk matriks V – WH apakah normnya adalah yang paling kecil.
Contoh 2: Misalkan V = A pada contoh 1, kemudian kita pilih dua buah faktor matriks W dan H yang ranknya adalah r =1 ≤ min{2,3} sebagai berikut
⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = 0 1 W dan , [ ] 1 1 1 = H
maka bentuk norm matriks V – WH adalah { } 2 2 ,1 max = = − ∞ WH V { } 2 0 , 2 ,1 max 1 = = − WH V ( ) ( ) ( ) WH V WH V WH V t svd − − = − max λ ( )( ) ( ) 6180 . 1 6180 . 2 max = = − − = t WH V WH V λ ( ) ( ) 3 0 1 1 0 1 0 2 1 2 2 2 2 2 2 2 = + + + + + = − WH V
Proses iterasi yang dilakukan untuk memperoleh norm minimum untuk matriks V – WH dilakukan dengan cara menentukan terlebih dahulu dua buah faktor awal W 0 dan H 0 yang ranknya r ≤ min{ n , m }, kemudian iterasi pertama akan menghasilkan dua buah faktor berikutnya dan seterusnya. Adapun iterasi yang dilakukan pada tiap entri-entri kedua faktor yang kita nyatakan dengan teorema berikut:
Teorema 1: Norm matriks V – WH merupakan barisan nilai-nilai yang monoton tak naik untuk bentuk iterasi entri-entri matriks faktor
( ) ( ) ik ik T ik T ik w WHH VH w → , untuk n i , , 2 , 1 L = dan r k , , 2 , 1 L = dan ( ) ( ) kj kj T kj T kj h WH W V W h → , untuk r k , , 2 , 1 L = dan m j , , 2 , 1 L = ,
asalkan ( ) 0 ≠ ik T WHH , untuk n i , , 2 , 1 L = dan r k , , 2 , 1 L = dan ( ) 0 ≠ kj T WH W , untuk dan . r k , , 2 , 1 L = m j , , 2 ,1 L = Contoh 3:
Untuk V pada contoh 2, misalkan matriks faktor awal diberikan oleh

## Jurnal Matematika, Statistika & Komputasi
205 Edisi Khusus Juli 2007 Jurnal Matematika, Statistika, & Komputasi Vol. 13  No.2 Januari  2017 149 Andi Kresna Jaya
⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = 1 1 1 W dan [ ] 1 1 1 1 = H
Maka dengan menggunakan program dalam Matlab 7 diperoleh 3 iterasi pertama
⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = 6667 . 0 3333 . 1 2 W dan ; [ ] 5 . 0 5 . 1 1 2 = H ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = 7143 . 0 2857 . 1 3 W dan ; [ ] 6 . 0 5 . 1 9 . 0 3 = H ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = 7018 . 0 3158 . 1 4 W dan [ ] 5943 . 0 5189 . 1 9245 . 0 4 = H
Nilai keempat norm untuk masing-masing faktor iterasi ditunjukkan pada tabel .
Tabel 1 .  Norm fmtn untuk 4 I terasi P ertama
W dan H 1 WH V − ∞ − WH V svd WH V − 2 WH V − W 1 , H 1 1 1 1 1.4142 W 2 , H 2 0.6667 0.6667 0.6667 0.6667 W 3 , H 3 0.6571 0.8571 0.6299 0.6312 W 4 , H 4

## 0.6351 0.8342 0.6282 0.6293
Tabel di atas menunjukkan bahwa norm svd WH V − mempunyai pendekatan yang lebih baik, walaupun untuk kasus lain hal itu belum tentu benar. Kasus di atas juga memperlihatkan kepada kita bahwa iterasinya berjalan lambat untuk faktor awal yang diberikan. Namun dengan proses iterasi yang terus-menerus maka kita akan mendapatkan hampiran terbaik untuk fmtn.

## 4.  Kesimpulan
Bentuk fmtn yang dapat dilakukan adalah meminimumkan norm WH V − dan mengasumsikan norm sebagai bentuk fungsi tujuan, maka fungsi kendalanya diberikan adalah W dan H adalah positif.
Algoritma untuk faktorisasi matriks tak negatif ini adalah
( ) ( ) kj kj T kj T kj h WH W V W h → dan
( ) ( ) ik ik T ik T ik w WHH VH w → untuk meminimumkan norm.
Jurnal Matematika, Statistika & Komputasi Edisi Khusus Juli 2007
206 Jurnal Matematika, Statistika, & Komputasi Vol. 13  No.2 Januari  2017
[2]   Mark D Plumbley. Samer A Abdallah, 2004. Polyphonic Music Transcription by Non-Negative Sparse Coding of Power Spectra, Centre for Digital Music , University of London.
[3]  Patrik O Hoyer, 2002, Non-Negative Sparse Coding, Proc. IEEE Workshop on Neural Networks for Signal Processing .
[4] Hiroki Asari. 2004. Non-negative Matrix Factorization: A possible way to learn sound dictionaries. Paper Tony Zador Lab , Watson school of Biological Sciences cold Spring Harbor Laboratory.
[5] Andi Kresna Jaya. 2001. Masalah Nilai Karakteristik Invers Real Tak Negatif, Tesis Magister . PPs Magister Matematika, Institut Teknologi Bandung.
[1] J udith C Brown. Paris Smaragdis. 2003. Non-Negative Factorization for Polyphonic Music Transcription, IEEE Workshop on Application of Signal Proccessing to Audio and Acoustics , pp 177 – 178.
Edisi Khusus Juli 2007

## Daftar Pustaka
Jurnal Matematika, Statistika, & Komputasi Vol. 13  No.2 Januari  2017 151 Andi Kresna Jaya 146 151