## ALGORITMA FILTER KALMAN PADA PENGUKURAN BEBAN MOTOR DC DENGAN METODE LABVIEW

## Sayuti Syamsuar
1

## Abstrak
Pemodelan matematika dan sistem identifikasi merupakan suatu proses yang lengkap dalam industry atau system dinamika. Demikian juga halnya struktur-struktur model linier mempunyai banyak kegunaan sebagai langkah awal dalam pengetahuan pemodelan  dan aplikasi sistem kontrol Lab View menggunakan algoritma Filter Kalman pada eksperimen motor DC di laboratorium.    Analisa dalam domain waktu (t) , kita menggunakan fungsi ramp, step dan sinusoidal sebagai fungsi input (t) terhadap model matematika orde 2 dari pembebanan dinamik motor DC. Dan, kita memperoleh respons waktu (t) yang berosilasi dengan faktor redaman ( ζ ) yang berbeda-beda pada kasus ini. Metodologi ini dapat dikembangkan ke arah analisa stabilitas dan pengendalian pada kapal dan pesawat udara .
Kata kunci : sistem identifikasi, motor DC,Filter Kalman, rasio redaman ( ζ ), fungsi transfer

## Abstract
The mathematical modeling and system identification are the complete development process on the industrial or the dynamics system. Although linear model structures have several advantages as a starting point to know the modeling and Lab View control system application by using Kalman Filter algorithm on the motor DC experiment in the Laboratory. In the time domain (t) analysis, we use the input ramp, step and sinusoidal on the second order of motor DC loading dynamic system. And, we found the oscillation time responses (t) with several damping ratio ( ζ ) on this cases. This methodology could be developed to stability and control analysis of ship and aircraft.
Keywords : identification system, motor DC, Kalman Filter, damping ratio( ζ ), transfer function

## PENDAHULUAN
Proses estimasi yang sering digunakan orang untuk mengetahui parameter dinamik dari suatu sistem yang dianalisa seperti metoda Maximum Likelihood, Extended Kalman Filter, Least Square dan lain-lain. Pada kasus ini digunakan metoda Filter Kalman untuk mengidentifikasi proses pembebanan pada motor listrik DC (Direct Current).
Perlunya perancangan lebih awal terhadap keperluan besarnya pembebanan optimum dari suatu
motor DC sehingga diperoleh ukuran kerja mesin yang dijual di pasaran dari suatu kapal laut atau pengendalian putaran propeller dari mesin pesawat udara. Sistem identifikasi ini dapat dikembangkan untuk analisa kestabilan dan kontrol dari wahana transportasi laut atau udara. Akibat banyaknya jenis dan ukuran motor DC yang dijual di pasaran, maka diperlukan juga perhitungan analisis dari kinerja mesin yang digunakan apakah sudah sesuai dengan keperluan. Ketidak sesuaian besarnya tenaga mesin yang menggerakkan propeller sering menyebabkan
in-efisiensi. Beberapa teknik pengukuran yang dilakukan pada eksperimen ini menggunakan optimalisasi algoritma Filter Kalman dari persamaan fungsi transfer sistem motor DC skala laboratorium. Pemahaman tentang motor DC, seperti pengendalian over speed atau lower speed pada propeller pesawat udara. Pembahasan tentang komponen serta jenis motor DC telah disampaikan pada beberapa paper sebelumnya (Syamsuar, dkk, 2011).
Tujuan dari penelitian ini adalah memperkenalkan teknik atau metoda perhitungan dan pengukuran skala laboratorium untuk mengetahui respons pembebanan motor DC terhadap input yang diberikan pada suatu sistem dinamik atau model matematik motor tersebut. Penggunaan Filter Kalman sebagai proses identifikasi model matematika pada eksperimen motor DC di laboratorium  Seperti diketahui bahwa motor DC digunakan sebagai penggerak propeller pada kapal laut dan mesin pesawat udara.
Adapun kegunaan dari penelitian ini adalah memperkenalkan beberapa contoh simulasi dalam bentuk model matematika persamaan orde 2 disampaikan pada tulisan ini. Beberapa parameter sistem dinamik diperkenalkan, demikian juga pengaruh input ramp, step dan sinusoidal yang diberikan terhadap sistem menghasilkan respons output akibat beberapa variasi parameter damping ratio ( ζ ). Proses identifikasi ini juga dapat dikembangkan ke sistem dinamik pesawat udara yang mempunyai persamaan diferensial orde tinggi.

## TINJAUAN PUSTAKA
Berikut ini diberikan beberapa contoh persamaan model matematika orde 2 yang juga merupakan model matematika yang sering ditemukan pada motor DC.
Beberapa contoh respons model sistem dinamik
Memperhatikan karakteristik respons output terhadap input yang diberikan pada suatu persamaan model matematika dengan fungsi transfer, H (s) orde 2. Persamaan diferensial orde dua:
dimana, u(t) :    input Sistem dinamik,
dimana,
n ω :  frekuensi natural
ς :  damping ratio Jika, ς 〈 1   :  under damped ς 〉 1   :  over damped ς =1    :  critical damping
Misalkan, bentuk input step, ramp dan fungsi sinusoidal yang diberikan kepada kedua fungsi transfer, H 1 (s) dan    H 2 (s) sebagai orde 2, seperti terlihat pada Gambar 1.
Contoh fungsi transfer, H 1 (s) orde 2 (pertama) adalah
;  (kurva warna kuning)
Contoh fungsi transfer, H 2 (s) orde 2 (kedua) adalah
;  (kurva warna merah)

## Penggunaan MATLAB dalam penyelesaian persamaan model matematika
Algoritma model matematika yang digunakan untuk menyelesaikan persoalan simulasi adalah penggunaan SIMULINK pada MATLAB. Pada Gambar 1 di bawah ini diperlihatkan penggunaan simulasi data tersebut terhadap dua model matematika orde 2 yang dipilih, sebagai berikut:
18 2s  +16s+5 2 Transfer Fcn6 18 2s  +20s+5 2 Transfer Fcn5 18 2s  +20s+5 2 Transfer Fcn4 18 2s  +16s+5 2 Transfer Fcn3 18 2s  +20s+5 2 Transfer Fcn2 18 2s  +16s+5 2 Transfer Fcn1 Step2 Step1 Sine Wave2 Sine Wave1 Scope3 Scope2 Scope Ramp2 Ramp1
Gambar 1    Diagram blok fungsi transfer, H (s) orde 2.
Untuk mengetahui bagaimana bentuk respons output (t) yang diberikan dan membandingkan kedua grafik tersebut di atas dan kenapa mempunyai respons waktu (t) yang berbeda. Perhatikan respons output (t) yang diperoleh dari input ramp untuk kedua fungsi transfer, H (s) orde dua tersebut. Fungsi transfer yang diwakili oleh kurva berwarna merah mempunyai damping ratio ( ζ ) lebih kecil dari kurva berwarna kuning pada Gambar 2, Gambar 3 dan Gambar 4.
Gambar 2    Respons ouput dengan input ramp terhadap kedua fungsi transfer H (s).
Respons output yang diperoleh dari input step (t) untuk kedua fungsi transfer H (s) orde 2 tersebut:
Gambar 3    Respons ouput dengan input step terhadap kedua fungsi transfer, H (s).
Respons output yang diperoleh dari input fungsi sinusoidal (t) untuk kedua fungsi transfer, H (s) orde 2 tersebut:
Gambar 4    Respons ouput dengan input fungsi sinusoidal terhadap kedua fungsi transfer, H (s).

## TEORI
Beberapa teori umum yang diuraikan merupakan proses stokastik dan probabilistik dalam persamaan matematika untuk mendekati keadaan yang sebenarnya yang diperoleh dari data pengukuran pembebanan motor DC di laboratorium.

## Beberapa aturan dalam algoritma Filter Kalman
Filter Kalman hampir sama dengan dengan filter lain, seperti band-pass filter, low-pass filter ataupun high-pass filter yang bekerja dalam respons frekuensi
( ω ). Bedanya hanya ada proses estimasi. Beberapa pengertian tentang PRBS (Pseudo Random Binary Sequence) dan Finite Impulse Response (FIR) dalam domain waktu (t), memberikan pengertian yang hampir sama, dimana sinyal input dalam bentuk fungsi ramp. Sebelum membahas teknik pemograman Filter Kalman, berikut ini diberikan beberapa uraian atau tahapan pekerjaan pemograman.

## Aturan umum pemberhentian program Filter Kalman
Untuk mengidentifikasi sistem dinamik suatu objek dengan beberapa persamaan matematik yang tidak diketahui, dan juga dengan koefisien yang tidak diketahui dan harus diprediksi; seperti seperti persoalan parameter. Hal itu juga selalu diikut sertakan dengan persoalan untuk menentukan kapan proses komputasional harus berhenti atau jumlah yang dihitung perlu estimasi dan ketelitian berapa yang merupakan bagian dari estimasi komputasional. Hasil dari teori matematika secara keseluruhan dari pemberhentian proses komputasional akhirnya masuk ke dalam proses optimalisasi aturan pemberhentian program. Hal ini akhirnya dikenal dengan istilah chi square distribution atau error ellipsoid.
Berikut ini disampaikan prosedur identifikasi untuk pemberhentian program yang diaplikasikan pada multi-dimensional parameter pada persoalan identifikasi.
dimana, i D θ Η :    matriks covariance dari estimasi Ө i dan Ө i-1 ε :  nilai kecil yang didefinisikan sebelumnya.
Untuk menguji konsistensi dari estimasi, adopsi
tingkat ketepatan α 1 dimana hubungan tersebut diperlihatkan oleh koefisien ketepatan
1 1 1 α β − = .
Definisikan threshold dari 2 β χ dalam  bentuk probabilitas menggunakan distribusi investigasi statistik 2 r .
[ ] 1 2 1 2 β χ χ β = 〈 P ; 1 0 1 〈 〈 β

## Pendeteksian kesalahan algoritma Filter Kalman
Perhatikan ekspresi persamaan berikut:
[ ] 2 2 2 2 β χ χ β = 〈 P ; 1 0 2 〈 〈 β
Apabila,
2 2 2 β χ ≤ i r :Filter Kalman masih bekerja normal.
2 2 2 β χ ≥ i r : kesalahan yang terjadi pada program.
Pada Gambar 5 terlihat kurva daerah kemungkinan Filter Kalman saat bekerja, merupakan fungsi densitas dari distribusi 2 χ dengan n = 4 derajat kebebasan, limit kepercayaan komputasional β 1 = 1- α 1 = 0.2 dan β 2 = 1- α 2 = 0.95 dimana α 2 merupakan tingkat kepercayaan paling tinggi.
Gambar 5    Kurva faktor kemungkinan ( probability ) dari Filter Kalman saat bekerja.
dimana, Daerah 1  :  operasi berhenti, Daerah 2  :  observasi, Daerah 3  :  pendeteksian kesalahan Gambar 5 tersebut di atas merupakan manifestasi dari
persamaan berikut ini:
2 0 2
1 ) ( 2 1 χ χ β β d f x ∫ = 2 0 2 2 ) ( 2 2 χ χ β β d f x ∫ = 2 2 1 2 ) ( 2 2 2 1 χ χ β β β β d f x x ∫ = −
Komputasional pada matriks covariance merupakan penyebaran antara dua estimasi yang diperhitungkan. Perhatikan bentuk persamaan linier sebagai berikut,
n i x y T i i , 1 = = θ
dimana,
] , , , 1 [ 2 m i i i T i p p p x = :    vektor dari variabel input ] ,...... , , [ 2 1 0 m T a a a a = θ :  vector dari parameter unknown sebagai parameter yang akan diestimasi. Sinyal output dari i y diukur oleh peralatan instrumentasi Lab View.
i i i h y z + = Apabila, plant mempunyai persamaan diferensial orde
2, maka
) 1 ( ) 2 ( ) 1 ( ) ( 0 2 1 − + − − − − = ∧ t u b t y a t y a t y Perhatikan persamaan: Bu Ax x + = . Cx y =
Maka, proses identifikasi menggunakan Kalman Filter, adalah
∧ ∧ ∧ = ψ θ . ) ( t y
⎥ ⎦ ⎤ ⎢ ⎣ ⎡ − − − − ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = ) 2 ( ) 2 ( ) 1 ( ) 1 ( t u i y t u i y b a Sehingga, ∧ ∧ ∧ = y P T . ψ θ Atau, ∧ − ∧ = y T T . . ] . [ 1 ψ ψ ψ θ

## Proses estimasi menggunakan algoritma

## Filter Kalman
Kalman filter awalnya digunakan untuk estimasi state x pada time discrete controlled process yang diatur oleh persamaan diffrensial linear:
x k+1 =  A k x k + B u k + v k
dimana,
A dan B adalah matrix, u input yang berfungsi sebagai kendali, dan v derau proses. Derau ini secara statistis juga bersifat tak gayut satu sama lain, terdistribusi normal, dengan nilai rata-rata nol dan kovarians Q=cov(u k ). Dengan perkataan lain, v ~ N(0,Q k ).
Jika P k = cov(x k ), maka operasi kovarians persamaan di atas hasilnya adalah
P k = A k P k-1 A k T + Q k
Nilai z yang ditampilkan oleh sebuah alat ukur yang mempunyai dua komponen, yaitu x yang memang harus diukur dan nilai lain w yang datang ke dalam alat ukur itu sebagai gangguan/derau pengukuran. Secara matematis z merupakan kombinasi linear dari x dan w:
z k = H k x k + w k Derau atau gangguan w k pada umumnya disepakati memiliki sifat statistis tak gayut satu sama lain, terdistribusi normal dengan nilai rata-rata nol dan memiliki  matrix kovarians R k = cov(w k ). Secara singkat biasanya ditulis w k ~ N(0,R k ). Operasi kovarians atas persamaan ini menghasilkan:
S k = P k H k P k T + R k
dimana,
P k = cov(x k ) = kovarians vector keadaan, dan S k = cov(z k ) = kovarians (data) pengukuran
Penggunaan Filter Kalman dalam aspek praktikal menggunakan sinyal input PRBS seperti diperlihatkan oleh Gambar 6 [3].
Gambar 6    Sistem dinamik PRBS dan Filter Kalman.

## ANALISIS DATA

## Penggunaan Filter Kalman pada model sistem dinamik motor DC
Proses identifikasi dari seperangkat motor DC di laboratorium menggunakan Lab View. Hal ini sebagai cara menentukan persamaan dari orde 2 untuk n a , n b dan d (time delay) dari sistem dinamik motor DC. Pemilihan data dalam bentuk grafik respons waktu (t) sebagai output, maka perlu diperhatikan besarnya frekuensi sampling, f s . Hasil pengukuran 999 data praktikum dengan faktor noise seperti terlihat pada Gambar 7.

## Gambar 7    Data Output, respons waktu (t) orde 2 dari pengukuran LabView
Fungsi input dari PRBS merupakan sinyal sequential dengan bentuk pulsa rectangular sebanyak 49 data. Data digital atau pengkodean dengan 0 dan 1 sebagai bilangan binary digunakan pada bagian input
ini. Ambil 49 data saja sebagai input seperti pada Gambar 8.

## Gambar 8    Data input PRBS pada motor DC
Respons output sebagai fungsi (t) dengan input PRBS pada motor DC yang diberi beban dengan setting di angka nomor 9 sebagai output transient dan steady state dengan 49 data yang diseleksi. Output transient dan steady state seperti pada Gambar 9.
Gambar 9    Data pengukuran Output, respons waktu (t) persamaan orde 2 [3] .
Apabila, plant mempunyai persamaan diferensial Orde 2, maka:
) 1 ( ) 2 ( ) 1 ( ) ( 0 2 1 − + − − − − = ∧ t u b t y a t y a t y
Maka, pada Gambar 7 merupakan manifestasi dari persamaan:
Bu Ax x + = . Cx y =
Maka, proses identifikasi menggunakan Filter Kalman, adalah:
∧ ∧ ∧ = ψ θ . ) ( t y
⎥ ⎦ ⎤ ⎢ ⎣ ⎡ − − − − ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = ) 2 ( ) 2 ( ) 1 ( ) 1 ( t u i y t u i y b a Sehingga,
∧ ∧ ∧ = y
P T . ψ θ
Atau, ∧ − ∧ = y T T . . ] . [ 1 ψ ψ ψ θ
Berdasarkan Tabel 1, jika dikelompokkan data dari 0 sampai dengan 48; maka diperoleh persamaan matriks dari no. (0 – 6); (7 – 13); (14 – 20); (21 – 27); (28 – 34); (35 – 41) dan (42 – 48) dalam beberapa kelompok. Data yang diambil mempunyai 7 kelompok data. Contoh pengambilan data:
⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ = ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ = ∧ 299 . 1 889 . 0 474 . 0 ) 5 ( ) 4 ( ) 3 (
y y y y ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ = ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ = 5 . 1 474 . 0 889 . 0 5 . 1 103 . 0 474 . 0 75 . 0 005 . 0 103 . 0 ) 4 ( ) 3 ( ) 4 ( ) 3 ( ) 2 ( ) 3 ( ) 2 ( ) 1 ( ) 2 ( u y y u y y u y y ψ
Hasil dari perhitungan dengan menggunakan MATLAB. Harga rata-rata yang diperoleh dari 7 kelompok data tersebut adalah:
a 1 rata-rata = (1.3985-0.138+0.0+0.0+0.0+0.0+0.0)/7=0.18007
a 2 rata-rata =
(-0.6448+0.1360+0.0+0.0+0.0+0.0+0.0)/7=-0.07268
b 0 rata-rata = (0.2120+1.000+1.000+1.000+1.000+1.000+1.000)/7= 0.88740
Hasil perhitungan koefisien a1, a2 dan b0 dari persamaan matematik orde-2 menggunakan algoritma ARMAX menggunakan beban nomor 9 pada motor DC seperti terlihat pada Gambar 10.
Gambar 10    Hasil perhitungan koefisien a1, a2 dan b0 sebagai fungsi (t) menggunakan algoritma ARMAX.

## Hasil perhitungan menggunakan algoritma Filter Kalman di MATLAB
Terdapat perbedaan antara parameter yang diestimasi dengan data simulasi. Hal ini disebabkan adanya time lag, noise dan kelembaman dari sistem pembebanan motor DC. Time lag, noise dan kelembaman yang terjadi pada sistem motor DC diperlihatkan oleh kurva error. Kurva hasil perhitungan algoritma Filter Kalman yang merupakan estimasi data pengukuran dan simulasi diperlihatkan oleh Gambar 11 yang merupakan persamaan diferensial orde 2 dari motor DC dengan menggunakan program MATLAB.
0 10 20 30 40 50 60 70 80 90 100 -4 -2 0 2 4 No. of samples Ou tp u t Kalman filter response 0 10 20 30 40 50 60 70 80 90 100 -4 -2 0 2 4 No. of samples Er ro r

## Gambar 11 Hasil perhitungan estimasi dan data pengukuran.
Pada Gambar 12 terlihat penggunaan motor DC pada sistem propulsi Remote Control (RC) pesawat udara jenis amphibian.

## Gambar 12 RC model pesawat udara amphibian

## KESIMPULAN
Terdapat beberapa kesimpulan yang dapat ditarik dari eksperimen motor DC menggunakan metodologi Lab View sebagai system control sebagai berikut:
• Terlihat adanya perbedaan antara simulasi data dengan estimasi dari data pengukuran menggunakan algoritma Filter Kalman.
Terlihat respons output orde 2 yang dihasilkan oleh proses estimasi menggunakan algoritma Filter Kalman adalah tidak smooth. Hal ini disebabkan adanya time delay, lag dan noise dari rangkaian sistem motor DC dari bagian awal sampai akhir sistem pada saat eksperimen dilaksanakan.
Untuk mengetahui solusi identifikasi parameter aerodinamika pesawat udara yang mempunyai persamaan diferensial orde tinggi, maka perlu melakukan penguraian dalam respons waktu (t) terhadap parameter atau state variable yang diseleksi dan data pengukuran sistem sensor peralatan instrumentasi yang dipasang pada sumbu koordinat pesawat udara tersebut.

## DAFTAR PUSTAKA
Zoran Vukic, Ljubomir Kuljaca, Dali Donlagic, Sejid Tesnjak, Non Linear Control Systems , Printed in the United States of America, New York, 2003. Chingiz Hajiyev, Stopping of algorithms and faults detection in Kalman Filter application , Istanbul Technical University, Turkey, proceeding of the 7 th. Mediterranean Conference on Control and Automation (MED99), Haifa, June, 28-30, 1999. Ari Santoso, Catatan Kuliah Estimasi Parameter , Departemen Teknik Elektro, Fakultas Teknik, Institut Teknologi Sepuluh Nopember, ITS, Surabaya, 2011. S. B. Bozic, Digital and Kalman Filtering , 2 nd.
Edition, University of Birmingham, Great Britain,
1979, Second Edition, 1994.
Ioan Dore Landau, System Identification and Control Design , Prentice Hall, New Jersey, USA, 1990. S. Syamsuar, R. Wibawaningrum, H. Makarim, Cara kerja dan penggunaan Motor Direct Current (DC) pada Kapal Selam , hal. 509-523, Jurnal Warta Penelitian, Kementerian Perhubungan, Badan Penelitian dan Pengembangan, ISSN. 0852-1824, Volume 23, Nomor 5, Jakarta, Mei
2011.
S. Syamsuar, Progress Report pada program S3, Jurusan Teknik Sistem dan Pengaturan Kelautan, Fakultas Teknologi Kelautan, Institut Teknologi
Sepuluh Nopember, ITS, Surabaya, 2011.