LINEAR QUADRATIC REGULATOR (LQR) UNTUK SISTEM DESKRIPTOR BERINDEKS SATU

Jurnal Konvergensi Vol. 4, No. 1, April, 2014

LINEAR QUADRATIC REGULATOR (LQR) UNTUK SISTEM DESKRIPTOR BERINDEKS SATU Muhammad Wakhid Musthofa Program Studi Matematika Universitas Islam Negeri Sunan Kalijaga Yogyakarta, [email protected]

Abstrak Dalam makalah ini dikaji masalah Linear Quadratic Regulator(LQR) pada interval waktu berhingga untuk sistem deskrirptor berindeks satu. Berdasarkan prinsip minimum Pontryagin akan dikonstruksikan persamaan differensial Riccati untuk masalah LQR tersebut yang berperan penting dalam penyelesaiannya.Selanjutnya, akan diturunkan teorema yang menyatakan keberadaan solusi masalah LQR pada sistem linear kuadratik. Kata kunci: linear quadratic regulator, sistem deskriptor berindeks satu, persamaan differensial riccati.

1. Pendahuluan Desain kendali optimal pada suatu sistem bekerja dengan tujuan mengendalikan sistem tersebut dengan biaya seminimal mungkin. Dalam hal sistem disajikan dengan beberapa persamaan diferensial dan fungsi yang menyatakan biaya yang dikeluarkan oleh kendali dinyakatan dalam bentuk fungsi kuadratik maka masalah kendali tersebut disebut kendali linear kuadratik. Pada umumnya kendali didesain untuk membawa state dari sistem yang dikendalikan menuju ke state nol ataupun ke sebuah persekitaran yang sangat dekat dengan nol. Masalah ini disebut sebagai masalah regulator. Sehingga masalah Linear Quadratic Regulator dapat dipandang sebagai masalah mencari kendali pada suatu sistem linear yang akan meminimumkan biaya pengendalian yang disajikan dalam bentuk fungsi kuadratik dengan tujuan pengendalian membawa state dari sistem tersebut menuju ke state nol ataupun ke sebuah persekitaran yang sangat dekat dengan nol. Masalah Linear Quadratic Regulator merupakan bagian dari masalah kendali optimal yang paling banyak dikaji dan dikembangkan baik dari sisi pengembangan teori maupun aplikasinya. Referensi dari masalah tersebut untuk sistem nonsingular dapat dilihat pada Bryson dan Ho (1975), Lewis (1995), Kirk (1998) dan Vinter (2000). Makalah ini akan mengkaji desain kendali optimal pada masalah Linear Quadratic Regulator untuk sistem deskriptor berindeks satu. Sistem deskriptor adalah generalisasi dari sistem biasa (sistem nonsingular). Sistem ini memuat persamaan diferensial dan sekaligus persamaan aljabar. Banyak permasalahan yang disajikan dalam sistem ini, diantaranya adalah proses-proses kimia (Kumar dan Dauotidis 1996), sistem sirkuit listrik (Newcomb 1981, Newcomb dan Dziurla 1989), sistem ekonomi (Luenberger 1977), interkoneksi antar sistem berskala besar (Luenberger dan Arbel 1977, Singh dan Liu 1973), sistem pada teknik mesin (Hemami dan Wyman 1979), sistem pembangkit daya (Scott 1979) dan sistem robot (Mills dan Goldenberg 1989).

Linear Quadratic Regulator (LQR) Untuk Sistem Deskriptor Berindeks Satu Muhammad Wakhid Musthofa

1


Makalah ini disajikan dengan runtutan alur sebagai berikut. Setelah pendahuluan, bagian kedua dari makalah ini menyajikan rumusan masalah yang akan diselesaikan, yaitu mencari kendali optimal pada masalah Linear Quadratic Regulator untuk sistem deskriptor berindeks satu. Pada bagian ini juga dipaparkan transformasi masalah Linear Quadratic Regulator dari sistem deskriptor ke sistem nonsingular. Selanjutnya pada bagian ketiga akan dikonstruksikan persamaan diferensial Riccati yang dibutuhkan dalam mencari kendali optimal dengan memanfaatkan hasil transformasi bada bagian sebelumnya. Bagian keempat adalah bagian inti dari makalah ini. Pada bagian ini disajikan teorema yang memberikan eksistensi solusi dari masalah Linear Quadratic Regulator. Terakhir, bagian kelimamenyajikan contoh numerik yang menggambarkan aplikasi dari sebagian teori yang telah dipaparkan pada bagian-bagian sebelumnya.

2. Rumusan Masalah Ex  t   Ax  t   Bu  t  , x  0   x0 (1) dengan  nr m . Vektor u    adalah rank  E   n , B  R  s

Diberikan sistem deskriptor

E, A  R

n  r  n  r 

,

kendali yang diberikan oleh desainer kendali pada sistem (1). Masalah yang ingin diselesaikan dalam makalah ini adalah masalah Linear Quadratic Regulator yaitu mendesain kendali u     s pada sistem (1) yang akan meminimalkan fungsi ongkos kuadratik tf

J  u     xT  t  Qx  t   u T  t  Ru  t   dt  x T  t f  Qt f x  t f  , (2) 0

dengan matriks-matriks matriks R

Q, R dan Qt f adalah matriks simetri dan diasumsikan

mempunyai invers. Selanjutnya dikarenakan kendali u     s yang

didesain meminimalkan fungsi ongkos (2) maka kendali u  s disebut kendali optimal dan dilambangkan dengan u     s .Sistem (1) dikatakan regular jika

det   E  A   0. Sistem (1) mempunyai solusi untuk setiap nilai awal yang konsisten jika dan hanya jika sistem (1) regular (Dai 1989). Jika sistem (1) regular maka berlaku teorema transformasi bentuk kanonik Weierstrass berikut. Teorema 1. (Gantmacher, 1959)Jika sistem deskriptor (1) regular maka terdapat dua matriks nonsingular X   X1 X 2  dan Y  Y1 Y2  sedemikian sehingga

0  A1 0  I T dan Y AX  Y T EX   n  0 I  (3)   0 N  r dengan A1 adalah matriks dalam bentuk Jordan yang elemen-elemennya nilai-nilai eigen dari A, I k adalah matriks identitas dan N adalah matriks nilpoten juga dalam bentuk Jordan. Jika sistem (1) regular, maka solusi dari sistem (1) berbentuk (Engwerda dan Salmah 2009)

2



x  t   X 1x1  t   X 2 x2  t  dengan X   X1

X 2  , Y  Y1T 

Y2T  , X 1 , Y1T  R  n  r n , X 2 , Y2T  R  n  r n 

, dan t J s t x1  t   e Jt x1  0    e  Y1 f  s  ds

x1  0    I n

0 X

0 1

x0

k 1

di

i 0

dt i

x2  t     N iY2

f t  ,

nilai awal bagi persamaan di atas diberikan oleh

0

k 1

di

i 0

dt i

I r  X 1x0    N iY2

f 0 .

Bilangan k adalah derajat kenilpotenan matriks N. Yaitu bilangan bulat k dengan k 1

k

sifat N  0 and N  0 . Indeks dari sistem deskriptor (1) dinyatakan dengan derajat kenilpotenan k dari matriks N. Jika matriks E nonsingular maka didefinisikan sistem tersebut berindeks nol. Selanjutnya berdasarkan Teorema 1. di atas, dengan

 x1  t   1 n  : X x  t  dengan x1  t   R dan  x2  t 

mendefinisikan variabel state yang baru 

x2  t   R r , maka sistem (1) dapat ditulis sebagai

In 0 

0  x1  t    A1   0  x2  t    0

0   x1  t   T    Y Bu  t  ,  I r   x2  t  

 x1  0   1    X x0  x2  0  

(4) atau

x1  t   A1x1  t   Y1Bu  t  , x1  0    I n

x2  t     0 I r Y T Bu  t   Y2 Bu  t  .

0  X 1 x0 (5)

(6)

Sedangkan fungsi ongkos (2) dapat ditulis sebagai tf

J u  

 T    x1  t  0 

 x t   x2T  t   X T QX  1  u T  t  Ru  t   dt  x1T t f Qt f x1 t f .    x2  t  

 

 

(7) Berdasarkan hubungan (6), fungsi ongkos (7) dapat disajikan sebagai


3


tf

  X1T QX1  X1T QX 2Y2 B   x1  t   T T   J  u     x1  t  u  t       BT Y T X T QX BT Y T X T QX Y B   u  t    0   2 2 1 2 2 2 2    uT  t  Ru  t   dt  x1T t f Qt f x1 t f 

 

tf





T

   x1  t  0



Q uT  t    T  V 

 

V   x1  t    T    dt  x1 t f Qt f x1 t f . u t R      

 

 

(8) Dengan demikian masalah LQR (1,2) ekuivalen dengan masalah LQR (5,8).

3. Persamaan Diferensial Riccati Pada bagian ini akan dikonstruksikan persamaan diferensial Riccati yang memegang peranan penting dalam menentukan eksistensi solusi dari masalah LQR (1,2) atau LQR (5,8). Transformasi masalah LQR sistem deskriptor (1,2) ke dalam masalah LQR sistem nonsingular (5,8) yang telah berhasil dilakukan mengakibatkan metode pengkonstruksian persamaan diferensial Riccati pada masalah LQR sistem deskriptor dapat menggunakan metode pengkonstruksian persamaan diferensial Riccati pada masalah LQR sistem nonsingular. Untuk kepentingan tersebut, dibentuk fungsi Hamiltonian (Lewis 1995)

 Q V   x1  t   T H  2  x1T  t  uT  t    T    1  t   A1x1  t   Y1Bu  t     V u t   R      2  x1T  t  Qx1  t   uT  t V T x1  t    x1T  t Vu  t   uT  t  Ru  t   

1T  A1 x1  t   Y1Bu  t   . (9) Berdasarkan prinsip maksimum Pontriagyn (Lewis 2006) syarat perlu penyelesaian optimal masalah di atas adalah

H  A1 x1  t   Y1Bu  t   H ii.   t      2Qx1  t   2Vu  t   A1T  t  ,  t f  2Qt f x1 t f x1 H iii. 0   2V T x1  t   2 Ru  t   BT Y1T  t  atau u 2 Ru  t   2V T x1  t   BT Y1T  t  . i.

x1  t  





 

 

Dengan mengasumsikan matriks R mempunyai invers diperoleh

u  t    R 1V T x1  t   1 BT Y1T  t  . 2

4

(10)



Selanjutnya didefinisikan  t   1   t  yang akan menghasilkan persamaanstate

2

and costate

x1  t   A1 x1  t   Y1Bu  t 

  t   Qx1  t   Vu  t  

(11)

A1T

t  ,

 

 

 t f  Qt f x1 t f

(12) Berikutnya dengan menggunakan metode sweep akan dicari persamaan diferensial Riccati untuk masalah LQR (5,8) sebagai berikut. Diasumsikan

  t   K  t  x1  t  dengan K  t f   Qt f .

(13)

Turunkan persamaan (13) terhadap variabel t dihasilkan   t   K  t  x1  t   K  t  x1  t  ,

(14)

Substitusikan persamaan (10) – (13 )ke dalam persamaan (14) dihasilkan





Qx1  t   V  R1V T x1  t   BT Y1T K  t  x1  t   A1T K  t  x1  t   K  t  x1  t   K  t  A1x1  t 



  

Y1B  R1V T x1  t   BT Y1T K  t  x1  t  , K t f  Qt f .

(15)

Persamaan (15) dapat pula ditulis sebagai

K  t  x1  t    A1T K  t  x1  t   K  t  A1x1  t    K  t  Y1Bx1  t   Vx1  t   R 1





 BT Y1T K  t  x1  t   V T x1  t   Qx1  t  ,

(16)

 

K t f  Qt f .

Dikarenakan persamaan (16) berlaku untuk semua nilai awal x1  0  maka diperoleh

K  t    A1T K  t   K  t  A1





  K  t  Y1B  V  R 1 BT Y1T K  t   V T  Q,

(17)

 

K t f  Qt f .

Persamaan (17) adalah persamaan diferensial Riccati yang diturunkan dari syarat perlu penyelesaian optimal masalahLQR (5,8).

3.1. Eksistensi Solusi Setelah berhasil dikonstruksikan persamaan diferensial Riccati pada masalah LQR, berikut disajikan teorema eksistensi solusi dari masalah LQR (1,2) yang ekuivalen dengan (5,8). Teorema 2. Diberikan masalah Linear Quadratic Regulator (1,2) yang ekuivalen dengan (5,8). Masalah tersebut mempunyai solusi untuk setiap nilai awal x1 jika dan hanya 0

jika persamaan diferensial Riccati (17) mempunyai solusi simetri K  


pada

5


 0, t f  . Jika masalah LQR tersebut mempunyai solusi maka solusi tersebut adalah   tunggal dan kendali optimal feedback tersebut berbentuk

u *  t    R 1V T x1  t   1 BT Y1T K  t  x1  t  . 2

(18)

Dalam struktur lingkar terbuka, bentuk kendali (18) adalah

u *  t    R 1V T x1  t   1 BT Y1T K  t    t , 0  x10 2

(19)

dengan matriks  memenuhi persamaan transisi





  t , 0   A  Y BR 1 V T  BT Y T K  t     t , 0  ,   0, 0   I .  1  1 1  Lebih lanjut, nilai optimal yang diberikan oleh kendali optimal u *  t  adalah

 

J u*  x1T0 K  0  x10 .

(20)

Pembuktikan dari teorema di atas dilakukan dengan beberapa langkah berikut. Pertama, akan ditunjukkan bahwa jika persamaan diferensial Riccati (17) mempunyai solusi simetri maka masalah LQR (1,2) atau (5,8) mempunyai solusi pada  0, t f  .

Bukti Teorema 2, bagian “  ” Menggunakan hubungan tf

 0

d T x1  t  K  t  x1  t  dt  x1T t f K t f x1 t f  x1T  0  K  0  x1  0  , dt





     

maka fungsi ongkos (8) dapat ditulis sebagai tf

V   x1  t       dt  R   u  t   

 Q J  u      x1T  t  uT  t    T   V   0  

x1T

 0  K  0  x1  0  

tf

J u  

x1T

 T Q  x1  t  uT  t        V T  0

t f  Qt

f

tf

d

 dt x1  t  K t  x1 t  dt T

0

   x1 t f 

 K tf

V   x1  t   d T x1  t  K  t  x1  t    R   u  t   dt



  Qt

 x1T  0  K  0  x1  0   x1T t f

f



  dt 

   x1 t f .

K tf

Kemudian, dengan menggunakan persamaan (5) dan (17), integrand dari persamaan di atas dapat ditulis sebagai (ketergantungan terhadap waktu dihilangkan)  xT  1

Q uT   T  V 

V   x1  d T x1 Kx1  u  R 1V T x1  1 BT Y1T Kx   2 R   u  dt



 

T

 R u  R

1 T

V x1  1 BT Y1T Kx 2

. Akibatnya, fungsi ongkos (8) dapat disajikan dengan persamaan

6





tf

u  R J 0

1 T

V x1  1 BT Y1T Kx 2

T

 R u  R



1 T

V x1  1 BT Y1T Kx 2 .

 x1T  0  K  0  x1  0 

Berdasarkan persamaan di atas didapat J  x1T  0  K  0  x1  0  untuk semua u dan kesamaan akan terjadi jika u    dipilih pada persamaan (18). Untuk membuktikan arah sebaliknya dari Teorema 2 di atas, diasumsikan bahwa masalah LQR (1,2) atau (5,8) mempunyai solusi u *    dengan x*  t  adalah



trayektori dari state yang dihasilkan oleh kendali u *    dan J * 0, x1

0

 adalah nilai

dari ongkos minimum yang dihasilkan oleh kendali u *    . Selanjutnya akan



ditunjukkan bahwa jika J * 0, x1



nilai minimum J * t , x1

0

0

 ada untuk sembarang nilai awal maka terdapat

 dari masalah optimisasi

tf

J u  

 T Q  x1  t  uT  t        V T  t 

V   x1  t    T    dt  x1 t f Qt f x1 t f R   u  t   

 

 

(21)

dengan kendala sistem dinamik

x1  s   A1x1  s   Y1Bu  s  , x1  t   x10 .

(22)

Kemudian dengan mengasumsikan bahwa terdapat nilai infimum yang dinyatakan dengan J inf  t , x0  maka dipenuhi





J inf t , x10  t   x1T0  t  P  t  x10

(23)

dengan tanpa mengurangi keumuman diasumsikan P  t  adalah matriks simetri. Hal di atas akan ditunjukkan dengan menggunakan bantuan beberapa lemma berikut.

Lemma 1. (Molinari, 1975) Jika

memenuhi identitas parallelogram, V  x  y   V  x  y   2 V  x   V  y  , untuk semua x, y (24) dan untuk setiap y, fungsi

suatu

fungsi

V

W  x, y  : V  x  y   V  x  y  mempunyai sifat W  x, y  adalah

fungsi kontinu dalam  pada saat   0 (25) Maka V  x  adalah berbentuk kuadratik.



Lemma berikut berguna untuk menunjukkan bahwa J inf t , x1

0



memenuhi

identitas parallelogram (24). Lemma 2. (Engwerda, 2005)


7




 ada maka J inf  t,  x1  juga ada dan J inf  t ,  x1    2 J inf  t , x1  . Jika J inf t , x1 0

0

0

0

Berdasarkan Lemma2 di atas dihasilkan akibat berikut.

Akibat 1. Diasumsikan terdapat nilai infimum pada persamaan (21) untuk setiap nilai awal x1

0

. Maka, terdapat matriks simetri P  t  sedemikian sehingga





J inf t , x10  x1T0  t  P  t  x10 . Berdasarkan lemma-lemma dan akibat di atas dapat dikonstruksikan bukti arah kanan dari Teorema 2 sebagai berikut. Bukti Teorema 2, bagian “  ” Diberikan persamaan diferensial Riccati (17). Berdasarkan teorema fundamental eksistensi dan ketunggalan solusi persamaan diferensial terdapat interval waktu



maksimal t1 , t f  dengan persamaan diferensial Riccati (17) mempunyai solusi tunggal. Diasumsikan t1  0 . Dari Teorema 2 bagian “  ” didapat bahwa pada interval waktu

 t1, t f 

masalah optimisasi (21,22) mempunyai solusi yang T

merupakan solusi minimum. Solusi minimum tersebut adalah x10  t  P  t  x10 . Dikarenakan hal ini berlaku untuk sembarang nilai awal x1

0

maka dapat



disimpulkan bahwa P  t   K  t  pada interval waktu t1 , t f  . Namun dikarenakan

P  t  terbatas pada  0, t f  maka demikian juga K  t  juga harus terbatas pada t1 , t f  . Hal ini berakibat bahwa persamaan diferensial Riccati (17) juga mempunyai solusi pada suatu interval waktu t2 , t f  untuk t2  t1 . Sehingga keberadaan interval t1 , t f  bukanlah interval yang maksimal. Hal ini menunjukkan pengasumsian t1  0 tidaklah tepat dan harusnya t1  0 .







4.

Contoh Numerik

Diberikan

masalah

Ex  t   Ax  t   Bu  t  ,

Linear

x  0   x0

Quadratic

Regulator

(26)

1

J  u     xT  t  Qx  t   u T  t  Ru  t   dt  x T  t f  Qt f x  t f  ,

(27)

0

8



0 0  0 1  1   2 0 , A , B   ,Q   E    , R  1 dan 1 1 1 2 0  0 1 1 0  1  , x0    . Menggunakan transformasi bentuk kanonik Weierstrass, Qt f    0 1   0 1 1  1 1 dapat dicari matriks nonsingular Y T   dan X     sedemikian sehingga 1 0   0 1 dengan

masalah LQR (26,27) ekuivalen dengan

x1  t   A1x1  t   Y1Bu  t  , x1  0    I n 1

Q J  u      x1T  t  uT  t    T   V   0

0  X 1 x0

(28)

V   x1  t    T    dt  x1 t f Qt f x1 t f R   u  t   

 

1 

1

Q

 





dengan X1  Y2T    , X 2  Y1T    , A1  1 ,  0 1 VT

 

(29)

V   2 2   . R   2 3 

Substitusikan matriks-matriks di atas ke dalam persamaan (17) menghasilkan persamaan diferensial Riccati 2 K  t   7 K  t   2K  t   3 , K 1  1

(30)

yang mempunyai solusi K  t   t . Selanjutnya, dengan menggunakan persamaan (18) didapat bentuk kendali optimal feedback bagi masalah di atas adalah *

u  t   5.5e

7 t

dan trayektori

optimal x1*

Gambar 1. kendali optimal u *  t 

t 

diberikan

* oleh x1

t 

 7 t t 2   4   e .

Gambar 2. trayektori optimal

x1*  t  Gambar 1 mengilustrasikan kendali optimal yang dilakukan untuk meminimalkan fungsi ongkos (29) terhadap sistem (28), sedangkan Gambar 2 menyatakan


9


trayektori optimal pada sistem (28) jika menggunakan kendali u *  t  untuk mencapai nilai optimal.

5. Kesimpulan Dalam makalah ini telah dikaji konsep Linear Quadratic Regulator (LQR) untuk sistem deskriptor berindeks satu. Dengan menggunakan bentuk kanonik Weierstrass, masalah LQR untuk sistem deskriptor dapat dibawa ke dalam bentuk LQR sistem nonsingular. Hal tersebut mengakibatkan persamaan diferensial Riccati yang memegang peranan penting dalam mencari eksistensi solusi masalah LQR tersebut dapat dikonstruksikan menggunakan metode pada sistem nonsingular. Selanjutnya telah dibuktikan pula teorema yang menjamin eksistensi solusi dari permasalahan tersebut. Contoh numerik sederhana yang menggambarkan penggunaan teorema yang telah dikonstruksikan telah pula disajikan di akhir makalah ini. Namun demikian kajian dalam makalah ini masih terbatas pada sistem deskriptor yang berindeks satu. Sehingga pengembangan masalah LQR untuk sistem deskriptor berorde tinggi merupakan obyek penelitian yang masih harus dikaji lebih lanjut.

10



6. Daftar Pustaka [1] Bryson, A.E. dan Ho, Y-C., (1975). Applied Optimal Control, Taylor and Francis, New York. [2] Dai, L. (1989). Singular Control Systems, Springer-Verlag, Berlin. [3] Engwerda, J.C.(2005) Linear Quadratic Dynamic Optimization and Differential Games, John Wiley & Sons, West Sussex. [4] Engwerda, J.C. dan Salmah, (2009). The Open-Loop Linear Quadratic Differential Game for Index One Descriptor Systems, Automatica, vol. 45, pp 585-592. [5] Gantmacher, F., (1959). Theory of Matrices, Chelsea Publishing Company, New York. [6] Hemami, H. dan Wyman, B. F. (1979). Modeling and Control of Constrained Dynamic Systems with Application to Biped Locomotion in the Frontal Plane, IEEE Transactions on Automatic Control, vol. 24, 526 – 535.

[7] Kirk, D.E., (1998), Optimal Control Theory: An Introduction, Dover Publications, New York. [8] Kumar, A. danDaoutidis, P. (1996). State-Space Realizations of Linear Differential Algebraic-Equation Systems with Control-Dependent State Space, IEEE Transactions on Automatic Control, vol. 41, 269 – 274. [9] Lewis, A. D. The Maximum Principle of Pontryaginin Control and in Optimal Control, catatan kuliah, tersedia dihttp://www.mast.queensu.ca/~andrew/ [10] Lewis, F.L. dan Syrmos, V.L., (1995), Optimal Control, John Wiley and Sons, New York. [11] Luenberger, D. G., (1977). Dynamic Equation in Descriptor Form, IEEE Transaction on Automatic Control, vol. 22, 312 – 321. [12] Luenberger, D. G. danArbel. (1977). Singular Dynamic Leontief systems, Econometrica, vol. 45, 991 – 995. [13] Mills, J. K. dan Goldenberg, A. A., (1989). Force and Position Control of Manipulators During Constrained Motion Tasks, IEEE Transactions on Robotics and Automation, vol. 5, 30 – 46. [14] Molinari, B.P. (1977). Solving Polynomial Sistem Using Continuation for Scientific and Engineering Problems, Prentice-Hall, New Jersey. [15] Newcomb, R. W. (1981). The Semistate Description of Nonlinear TimeVariable Circuits, IEEE Transactions on Circuits Systems, vol. vol. 28, 62 – 71. [16] Newcomb, R. W. dan Dziurla, B. (1989). Some Circuits and Systems Applications of Semistate Theory, Circuits Systems Signal Processes, vol. 8, 235 – 260.

[17] Singh, S. dan Liu, R. W. (1973). Existence of State Equation Representation of Linear Large-Scale Dynamical Systems, IEEE Transaction Circuits Systems, vol. 20, 239 – 246. [18] Scott, B., (1979). Power system Dynamic Response Calculations, IEEE Proceeding, vol. 67, 219 – 247. [19] Vinter, R.B., (2000). Optimal Control, Birkhauser, Boston.


11

LINEAR QUADRATIC REGULATOR (LQR) UNTUK SISTEM DESKRIPTOR BERINDEKS SATU

Recommend Documents