Filter Kalman adalah suatu algoritma yang digunakan untuk mengestimasi variabel keadaan pada sistem linier. Sedangkan untuk sistem taklinier Filter Kalman tidak dapat digunakan secara langsung. Untuk beberapa tahun kemudian muncul metode baru yang dikenal dengan nama Unscented Kalman Filter (UKF) yang menggunakan teknik transformasi unscented. Akar Kuadrat merupakan suatu skema yang dapat diterapkan dalam UKF untuk mengestimasi suatu model dinamik taklinear. Dalam makalah ini dilakukan suatu kajian mengenai skema Akar Kuadrat yang diterapkan pada Unscented Kalman Filter (UKF) sehingga terbentuk suatu algoritma baru yang dinamakan dengan Akar Kuadrat-Unscented Kalman Filter (AK-UKF). Dan kemudian algoritma ini diimplementasikan pada model sistem deteksi kerak pada alat penukar panas

Erna, Apriliani

Tholib, M.

Indonesian

Lumbung Pustaka UNY  (UNY Repository)

Prosiding Seminar Nasional Penelitian, Pendidikan dan Penerapan MIPA,  Fakultas MIPA, Universitas Negeri Yogyakarta, 14 Mei 2011 M - 9 SKEMA AKAR KUADRAT DALAM UNSCENTED KALMAN FILTER   UNTUK MENDETEKSI KERAK PADA ALAT PENUKAR PANAS  M. Tholib1, Erna Apriliani2  Jurusan Matematika FMIPA Institut Teknologi Sepuluh Nopember Surabaya Email: 101tholib@gmail.com  Abstrak  Filter Kalman adalah suatu algoritma yang digunakan untuk mengestimasi variabel keadaan pada sistem linier. Sedangkan untuk sistem taklinier Filter Kalman tidak dapat digunakan secara langsung. Untuk beberapa tahun kemudian muncul metode baru yang dikenal dengan nama Unscented Kalman Filter (UKF) yang menggunakan teknik  transformasi unscented. Akar Kuadrat merupakan suatu skema yang dapat diterapkan dalam UKF untuk mengestimasi suatu model dinamik taklinear. Dalam makalah ini dilakukan suatu kajian mengenai skema Akar Kuadrat yang diterapkan pada Unscented Kalman Filter (UKF) sehingga terbentuk suatu algoritma baru yang dinamakan dengan Akar Kuadrat-Unscented Kalman Filter (AK-UKF). Dan kemudian algoritma ini diimplementasikan pada model sistem deteksi kerak pada alat penukar panas. Kata kunci: Filter Kalman, Unscented Kalman Filter (UKF), Akar Kuadrat-Unscented Kalman Filter (AK-UKF), kerak,  alat penukar panas.  PENDAHULUAN Filter Kalman merupakan salah satu metode untuk menaksir variabel keadaan dari sebuah sistem linier dengan meminimumkan kovariansi kesalahan estimasi. Salah satu pendekatan Filter Kalman yang bisa digunakan untuk sistem taklinier, yaitu Extended Kalman Filter (EKF) dan Unscented Kalman Filter (UKF)[4]. Selama kurang lebih 20-30 tahun EKF diakui secara umum sebagai metode untuk menaksir variabel keadaan sistem taklinier, hingga ditemukan suatu pendekatan baru yang dikenal dengan Unscented Kalman Filter [4]. Perbedaan kedua metode ini yaitu pada metode pelinieran yang digunakan, pada EKF sistem dilinierkan dengan menggunakan pendekatan deret Taylor, sementara pada UKF menggunakan Transformasi Unscented. Metode dasar Transformasi Unscented pertama kali diperkenalkan Uhlman dan Julier [5], mereka membentuk sebuah kerangka yang mewakili variabel random yang dinamakan Titik-Titik Sigma. Skema akar kuadrat merupakan salah satu skema yang dapat diimplementasikan pada UKF. Skema ini dapat mempengaruhi pada hasil estimasi menjadi lebih baik, baik dalam hal tingkat akurasi maupun waktu komputasi yang digunakan[6]. Dalam penelitian ini akan dilakukan suatu kajian mengenai implementasi skema akar kuadrat pada UKF, yang selanjutnya diterapkan untuk mengestimasi variabel keadaan pada sistem dengan model pengukuran taklinear yaitu sistem/model pada alat penukar panas. Hasil estimasi dengan metode Akar Kuadrat Unscented Kalman Filter (AK-UKF) selanjutnya akan dibandingkan metode UKF standar, sehingga diperoleh metode yang memiliki tingkat ketelitian yang lebih baik.                                                   1 Mahasiswa Pasca Sarjana Matematika FMIPA ITS Surabaya 2 Staf Pengajar Jurusan Matematika ITS Surabaya M. Tholib / Skema Akar Kuadrat M - 10 UNSCENTED KALMAN FILTER (UKF) Proses estimasi pada metode Filter Kalman menyajikan bentuk umum dari sistem yang digunakan dan beberapa tahapan dari proses estimasinya. Lewis  memberikan suatu sistem dinamik linear, secara umum berbentuk sebagai berikut[3]: kkkkkkk wGuBxAx ++=+1  kkkk vxHz +=     (1) ( )00 0~ , xx x P ; ( )kk Qw ,0~ ; ( )kk Rv ,0~  dengan kx  variabel keadaan sistem pada waktu k yang nilai estimasi awalnya 0x  dan kovarian awal 0xP ,  nkx ∈ℜ , ku  variabel input deterministik pada waktu k, mku ∈ℜ . kw  gangguan (noise) pada sistem yang mempunyai mean 0kw =  dan kovarian kQ , kz  variabel pengukuran, pkz ∈ℜ , kv  gangguan (noise) pada pengukuran dengan mean 0kv =  dan kovarian kR . kkkk HGBA ,,, adalah matriks-matrik dengan ukuran yang bersesuaian. Proses estimasinya ada dua tahap, yaitu tahap prediksi (time update) dipengaruhi oleh dinamika sistem, dan tahap koreksi (measurement update) dipengaruhi oleh informasi dari pengukuran. Kedua tahap ini akan berulang terus-menerus sampai pada waktu k yang ditentukan. Unscented Kalman Filter adalah pengembangan dari Filter Kalman untuk sistem yang non-linear dengan menggunakan teknik Transformasi Unscented. Misalkan diberikan suatu fungsi kepadatan peluang diskrit yk=f(xk,k) mempunyai variabel random x dari sebuah model taklinear dengan dimensi L mempunyai mean  dan kovarian . Fungsi yk=f(xk,k)  didekati dengan transformasi unscented. Mean dan kovarian tersebut digunakan untuk menentukan penyebaran 2L+1 titik-titik sigma disekitar . Titik-titik sigma dalam bentuk vektor sigma  diperoleh dengan menggunakan persamaan berikut [4]:    ,                      i = 1,........., L    i = L + 1,......., 2L     (2) dengan:  adalah parameter penskalaan, α adalah sebuah konstanta yang digunakan untuk menentukan sebaran dari titik sigma di sekeliling , dimana α selalu bernilai positif kecil dan  adalah skala penskalaan tambahan, dimana nilai ≥ 0. Nilai yang paling sering digunakan yaitu  = 0. Misalkan diberikan variabel keadaan: [ ]1 2 TLx x x x= L         (3) Jika dinyatakan dalam bentuk matriks sigma points bisa dituliskan menjadi: [ ]0 1 1 2 2 Ti L L L Lχ χ χ χ χ χ+ +=χ L L      (4) Karena ( )k ky f x= , maka penyebaran vektor sigma yk adalah:  ,        i = 0,........, 2L,  Pembobot mean dan kovarian dihitung berdasarkan persamaan [6]:                   (5)  ,   i = 1, ...., 2L Prosiding Seminar Nasional Penelitian, Pendidikan dan Penerapan MIPA,  Fakultas MIPA, Universitas Negeri Yogyakarta, 14 Mei 2011 M - 11 Dengan menggunakan titik-titik sigma persamaan (2) dan persamaan pembobot mean-kovarians pada persamaan (5), maka diperoleh mean:  ( )2 1 ( )0ˆLmi iiy W f χ+= =  ∑              (6) Sedangkan untuk menghitung kovarians dari y menggunakan persamaan: ( )( ) ( )( )2 ( )0ˆ ˆL Tcy i i iP W f y f yχ χ = − − ∑                      (7) Secara ringkas algoritma Unscented Kalman Filter dapat dituliskan sebagai berikut[5]: Inisialisasi Pada saat k=0      Untuk k=1,2,3,…,N :  Hitung titik sigma   dengan:     dan   Tahap prediksi (time update)          Tahap koreksi (measurement update)      Dengan  kQ = kov. eror proses dan kR = kov. eror pengukuran  AKAR KUADRAT UNSCENTED KALMAN FILTER (AK-UKF) Algoritma UKF standar dapat dimodifikasi dengan cara menyebarkan matriks akar kuadrat/faktor Cholesky  secara langsung, untuk menghindari kebutuhan pemfaktoran ulang pada setiap time step sehingga bisa mengurangi beban komputasi. Untuk memperoleh nilai estimasi, kovarian error dan penghitungaan faktor Cholesky tahap prediksi dan koreksi adalah dengan menggunakan teknik Transformasi Unscented. Secara urut algoritma akar kuadrat pada UKF dapat disusun sebagai berikut:   M. Tholib / Skema Akar Kuadrat M - 12 Inisialisasi Pada saat k = 0    Faktor Cholesky S0 :     Faktor Cholesky tambahan   :   Misal diberikan variabel keadaan:            (8) Dengan kovarian awal: 11 12 121 22 21 2LLkL L LLP P PP P PPP P P   =    LLM M ML                                                        (9) Berdasarkan model yang diberikan variabel random xk dengan dimensi L mempunyai mean ˆkx  , kovarian kP  dan juga faktor Cholesky Sk  yang diperoleh dari persamaan  , atau bisa ditulis: . Kemudian didefinisikan variabel keadaan tambahan: [ ]Tak k kx x w=                                                           (10) Sedangkan kovarian awal tambahan: xwa k kk xwk kP PPP Q =                                           (11) Dan faktor Cholesky awal tambahan adalah:             (12) Tk adalah factor Cholesky dari Qk . Selanjutnya mean, kovarian dan faktor Cholesky tersebut digunakan untuk menentukan penyebaran 2L+1 titik-titik sigma disekitar . Titik-titik sigma dalam bentuk vektor sigma   diperoleh dengan menggunakan persamaan berikut[7]:   ,  i = 0  ,  i  = 1,2,…L  ,  i = L+1,…,2L                           (13) Dengan  adalah skala parameter   dan  ,  adalah elemen baris ke-i dari , sedangkan L adalah dimensi variabel tambahan. Berdasarkan variabel keadaan diatas persamaan matriks sigma poin bisa dituliskan menjadi:                        (14) Sedangkan titik-titik sigma untuk xk+1 diperoleh dari: ( ), 1 , ,,i kx x wi k i kFχ χ χ+ =                         (15) Tahap Prediksi: Dengan menggunakan titik-titik sigma persamaan (15) dan pembobot mean dan kovarian persamaan (5) diperoleh hasil estimasi:              (16) Prosiding Seminar Nasional Penelitian, Pendidikan dan Penerapan MIPA,  Fakultas MIPA, Universitas Negeri Yogyakarta, 14 Mei 2011 M - 13 Langkah berikutnya adalah menentukan faktor Cholesky prediksi ( )  dengan cara melakukan faktorisasi QR terlebih dulu dari persamaan:                 (17) Lalu menghitung Update Faktor Cholesky  dari dengan persamaan:   cholupdate                    (18) Sigma poin dari model pengukuran, yaitu:                     (19) Dengan menggunakan persamaan bobot (5) dan persamaan (19) diperoleh persamaan estimasi dari model pengukuran:                  (20)  Tahap Koreksi  Menghitung faktor Cholesky  dengan melakukan faktorisasi QR  dari persamaan:               (21) Lalu menghitung Update Faktor Cholesky dari:   cholupdate ( )              (22) Kovarian error model proses dan pengukuran (cross covarian) diperoleh dari persamaan:         (23) Kalman Gain diperoleh dari persamaan:                              (24) Persamaan estimasi tahap koreksi adalah:                    (25) Menghitung matrik U dengan persamaan:                       (21) Menghitung faktor Cholesky tahap koreksi   dari persamaan:  cholupdate                               (22) MODEL /  SISTEM PADA ALAT PENUKAR PANAS Alat penukar panas atau dalam industri kimia populer dengan istilah Heat Exchanger (HE), merupakan suatu alat yang berfungsi untuk memindahkan panas antara dua fluida yang berbeda temperatur dan dipisahkan oleh suatu sekat pemisah. Perpindahan panas antara dua fluida dapat dipastikan mengakibatkan kerak[2]. Cukup banyak kerugian yang dapat ditimbulkan oleh kerak tersebut. Sehingga apabila terjadi perpindahan panas yang besar maka penting untuk dapat mendeteksi kerak yang terbentuk pada alat penukar panas.   Gambar 1. Alat Penukar Panas Tipe counter-flow  Dari gambar diatas persamaan / model yang mewakili alat penukar panas adalah[2]: M. Tholib / Skema Akar Kuadrat M - 14 ,1 ,1,2 ,2,1 ,1,2 ,211 . 02 2 21 11 . 1 . 02 2 211 . 02 2 21 10 1 . 1 .2 2 2h h hh hh hh h hc cc cc c cc c cT TT TdT TdtT Tα α ατ τ τα α ατ τ τβ β βτ τ τβ β βτ τ τ  − +           − − +            =    − +               − − +        ,,11 . 020210 1 .202hh inhc inccTTατατβτβτ  −        +      −                      (23) dengan 1,hT  adalah temperatur pada bagian panas 1( C° ), 2,hT  adalah temperatur pada bagian panas 2 = temperatur keluar (outlet) dari fluida panas( C° ), 1,cT  adalah temperatur pada bagian dingin 1( C° ), 2,cT  adalah temperatur pada bagian dingin 2 = temperature keluar (outlet) dari fluida dingin( C° ),  inhT ,  adalah temperatur inlet (masuk) pada bagian panas( C° ), incT ,  adalah temperatur inlet (masuk) pada bagian dingin( C° ).  Model tersebut memiliki 4 parameter yang dapat dinyatakan dalam bentuk berikut:  ( ) ( ) hhhctmUAt&=α , ( ) ( )tmMthhh&=τ , ( ) ( ) cccctmUAt&=β , ( ) ( )tmMtccc&=τ ,  dengan:  dan  adalah  jumlah unit perpindahan panas pada fluida panas dan dingin,  dan  adalah waktu yang dibutuhkan perpindahan panas fluida panas dan dingin,  dan  adalah  area perpindahan panas pada fluida panas dan dingin (m2), U adalah koefisian perpindahan panas menyeluruh, diasumsikan U konstan, M adalah massa fluida,  dan  adalah konstanta pada bagian panas dan bagian dingin. Misalkan model state: ( ) ( )0, , indx t f m x Tdt =   & ( ),1 ,2 ,1 ,2h h c cx t T T T Tα β Τ =   [ ] [ ]ΤΤ == incinhinch TTTmmm ,.,&&&  adalah laju alir massa pada bagian panas (kg/s),  adalah laju alir massa pada bagian dingin (kg/s). Didapatkan persamaan model state sebagai berikut:    Prosiding Seminar Nasional Penelitian, Pendidikan dan Penerapan MIPA,  Fakultas MIPA, Universitas Negeri Yogyakarta, 14 Mei 2011 M - 15 ,1 ,1 ,2 ,,1,1 ,2 ,1 ,,2,1,1 ,2 ,1 ,,2,1 ,1 ,2001 2 1 22 21 2 1 22 21 2 1 22 21 2 1 22 2h c c h inh h h hhh h c c inh h h h hch h c c inc c c cch c cc c c cT T T TTd T T T TTdtTT T T TTT T Tα α α α αβ τ τ τ τα α α ατ τ τ τβ β β βτ τ τ τβ β β βτ τ τ τ  + −− + + +   − +− + +=    + −+ − +   − ++ − +,h inT                 Selanjutnya dari 6 parameter pada model di atas, variabel yang akan diestimasi pada maklah ini adalah  dan  , dengan tujuan untuk mengetahui jumlah unit perpindahan panas pada fluida panas dan dingin sebagai dasar deteksi adanya kerak. Dan dengan cara yang sama, dapat diestimasi pula nilai dari 4 parameter yang lain yaitu ,  ,   dan  untuk mengetahui pola distribusi temperatur pada masing-masing bagian. Model diatas disimulasikan dengan nilai awal yang definisikan sebagai berikut[2]: ( )1170ˆ 0 ;603570x    =             ( )0.1 0 0 0 0 00 0.1 0 0 0 00 0 0.1 0 0 000 0 0 0.1 0 00 0 0 0 0.1 00 0 0 0 0 0.1P    =        Kemudian Q (kovarian noise pada sistem) dan R  (kovarian noise pada pengukuran) dapat didefinisikan sebagai matriks diagonal berikut: 0.01 0 0 0 0 00 0.01 0 0 0 00 0 0.01 0 0 0 0.01 0;0 0 0 0.01 0 0 0 0.010 0 0 0 0.01 00 0 0 0 0 0.01Q R      = =          HASIL SIMULASI Pada gambar 2 terlihat bahwa dengan 60 iterasi grafik estimasi  dengan metode AK-UKF relatif lebih mendekati grafik realnya dibandingkan dengan grafik dengan metode UKF standar. Kemudian grafik estimasi  dengan metode AK-UKF juga jauh lebih mendekati grafik realnya dibandingkan dengan grafik dengan metode UKF. Hal menunjukkan bahwa  dari 60 kali iterasi estimasi  dan  dengan metode AK-UKF lebih akurat daripada metode UKF standar. M. Tholib / Skema Akar Kuadrat M - 16 0 10 20 30 40 50 600.20.40.60.811.21.41.61.8Estimasi Nilai Alphawaktu ke k Nilai Alpha  realAK-UKFUKF0 10 20 30 40 50 6000.20.40.60.811.21.41.61.8Estimasi Nilai Betawaktu ke kNilai Beta  realAK-UKFUKF Gambar 2. Hasil Estimasi  dan  menggunakan UKF dan AK-UKF   Perbandingan nilai error kedua metode pada estimasi  dan  terlihat pada gambar 3, dapat dikatakan bahwa untuk waktu sampai ke-60 satuan waktu,  nilai error estimasi  dan  dengan menngunakan metode AK-UKF cenderung stabil pada kisaran 0.01–0.75 relatif lebih kecil daripada dengan menggunakan metode UKF yang berada pada kisaran 0.25-1.00. Dapat dikatakan pula bahwa nilai estimasi  dan  dengan AK-UKF memiliki error maksimal kurang dari 0.75, yang relatif kecil dibanding dengan metode UKF standar yang mencapai 1.00. Artinya pada estimasi  dan  nilai error dengan metode UKF standar, jauh lebih besar dibandingkan dengan metode AK-UKF. Jadi, dapat disimpulkan bahwa metode AK-UKF lebih akurat daripada metode UKF standar. 0 10 20 30 40 50 6000.51Error pada α  0 10 20 30 40 50 6000.511.5Error pada β  AK-UKFUKFAK-UKFUKF Gambar 3. Error Estimasi  dan   dengan Metode UKF dan AK-UKF  Selanjutnya akan dilakukan estimasi sampai dengan 400 satuan waktu. Dengan waktu yang lebih lama, diprediksikan kerak telah timbul atau terakumulasi. Gambar 4 menunjukkan bahwa ketika memasuki waktu ke-200 telah terjadi penurunan drastis nilai estimasi dari  dan  yaitu dari kisaran (-0.5)oC–(1.5) oC sampai pada kisaran (-1.3) oC – (-2.6) oC. Maka secara teoritis ataupun fisis dapat dikatakan telah terdeteksi adanya kerak pada alat penukar panas.   Prosiding Seminar Nasional Penelitian, Pendidikan dan Penerapan MIPA,  Fakultas MIPA, Universitas Negeri Yogyakarta, 14 Mei 2011 M - 17 0 50 100 150 200 250 300 350 400-3-2.5-2-1.5-1-0.500.511.5Estimasi Nilai Alpha dan Betawaktu ke k Nilai Alpha dan Beta  Estimasi AlphaEstimasi Beta Gambar  4.  Estimasi dari  dan  sampai waktu ke -400  Untuk  pada waktu ke-75 satuan waktu, nilai dari estimasinya sudah terlihat mulai mengalami penurunan dibawah 0oC, dan secara terus menerus menurun. Mulai terjadi penurunan drastis setelah waktu ke-200 sampai pada suhu (-1.4)oC. Sedangkan untuk  sampai pada waktu ke-75 satuan waktu nilai estimasinya terus mengalami penurunan sampai mendekati suhu (-0.5)oC, meskipun kemudian sempat mengalami kenaikan suhu pada waktu ke-180, namun berikutnya setelah waktu ke-200 mengalami penurunan drastis sampai pada titik  (-2.6) oC. Hal ini menunjukkan adanya kerak pada alat tersebut ketika memasuki waktu ke-200. Perbandingan waktu komputasi algoritma UKF dan AK-UKF pada masing-masing iterasi dapat dilihat pada Tabel 1 berikut:  Tabel.1 Perbandingan waktu komputasi algoritma UKF dan AK-UKF Jumlah iterasi Waktu Komputasi (detik) Metode UKF Metode AK-UKF k=60 2.3879 1.5740 k=100 4.7372 2.9715 k=600 10.5155 5.8328 k=1000 35.4859 14.6148 Nilai  Root Mean Square Error (RMSE) dari tiap-tiap variabel pada iterasi k=60, k=100 dan k=400  ditunjukkan pada Tabel 2. Tabel 2 Perbandingan RMSE Algoritma UKF dan AK-UKF tiap variabel.   Variabel keadaan Root Mean Square Error (RMSE) k=60 k=100 k=400 UKF AK- UKF UKF AK- UKF UKF AK- UKF α 2.5646 1.8775 3.5354 2.1375 10.4661 2.4596 β  3.7862 1.1688 4.9112 1.4220 37.4690 1.9113  KESIMPULAN DAN SARAN Dari tabel perbandingan rata-rata error antara UKF standar dan AK-UKF dapat disimpulkan bahwa metode AK-UKF jauh lebih baik dan akurat daripada UKF standar. Sedangkan dari tabel perbandingan waktu komputasi yang digunakan dari kedua algoritma tersebut dapat disimpulkan pula bahwa metode AK-UKF membutuhkan waktu yang lebih sedikit  daripada UKF standar. Jadi M. Tholib / Skema Akar Kuadrat M - 18 secara keseluruhan dapat disimpulkan berarti metode AK-UKF lebih akurat dan efisien dibanding dengan metode UKF standar. Penulis menyarankan untuk mengkaji lebih lanjut penerapan metode AK-UKF pada model sistem yang taklinier dan berorde yang lebih tinggi.  DAFTAR PUSTAKA [1] Golub, H. G. dan Loan, V. F. Charles. (1996), “Matrix Computations (Third Edition)”, The John Hopkins University Press, Baltimore and London. [2] Jonsson, G.R., Lalot, S., Palsson, O.P., dan Desmet, B. (2007). “Use of Extended Kalman Filtering in Detecting Fouling in Heat Exchangers”. University of Iceland, France. [3] Lewis, L Frank. (1986), “Optimal Estimation, With An Introduction To Stochastic Control Theory”, John Wiley and Sons, New York. [4] Rudi. (2006), ”Penerapan Extended Kalman Filter dan Unscented  Kalman Filter pada Estimasi Variabel Keadaan Sistem dengan Model Pengukuran Taklinier”, Tesis Magister Jurusan Matematika FMIPA Institut Teknologi Sepuluh Nopember, Surabaya. [5] Terejanu, Gabriel.A.(2003), “Unscented Kalman Filter Tutorial”, Departement of  Computer and Engineering University at Buffalo. [6] Wan, Eric dan Merwe, Rudolph V.D. (2001), “ The Square Root Unscented Kalman Filter For State and Parameter Estimation”. Oregon Graduate Institute of Science and Technology USA.   

SKEMA AKAR KUADRAT DALAM UNSCENTED KALMAN FILTER UNTUK MENDETEKSI KERAK PADA ALAT PENUKAR PANAS

Abstract

Similar works

Full text

Available Versions

Lumbung Pustaka UNY (UNY Repository)