Syarat Fritz John ... (Caturiyati)
SYARAT FRITZ JOHN PADA MASALAH OPTIMASI BERKENDALA KETAKSAMAAN Caturiyati1 Himmawati Puji Lestari2 1,2 Jurusan Pendidikan Matematika FMIPA UNY 1
[email protected] 2
[email protected]
Abstrak Syarat-syarat optimasi masalah optimasi berkendala merupakan perluasan syaratsyarat optimasi masalah optimasi tak berkendala. Pada paper ini diuraikan syarat Fritz John yang merupakan perluasan syarat optimasi masalah optimasi tak berkendala untuk menyelesaikan masalah optimasi berkendala. Kata kunci : syarat Fritz John, masalah optimasi berkendala, masalah optimasi tak berkendala
PENDAHULUAN Suatu masalah tak berkendala adalah suatu masalah berbentuk: meminimalkan π(π) tanpa vektor π sebagai suatu kendalanya. Masalah tak berkendala jarang muncul dalam penggunaan. Masalah yang lebih sering muncul adalah masalah optimasi berkendala. Masalah optimasi berkendala terbagi berdasarkan kendalanya, yaitu masalah optimasi berkendala ketaksamaan dan masalah optimasi berkendala campuran (kendala persamaan dan ketaksamaan). Syarat optimasi masalah optimasi berkendala merupakan perluasan syarat optimasi masalah optimasi tak berkendala. Fritz John mengemukakan syarat optimasi masalah optimasi berkendala ketaksamaan. Pembahasan pada paper ini terbagi menjadi beberapa bagian, yaitu mengenai syarat perlu dan cukup optimasi tak berkendala dan syarat Fritz John untuk masalah optimasi berkendala ketaksamaan.
Syarat Perlu dan Cukup Optimasi Tak Berkendala Masalah optimasi tak berkendala adalah suatu masalah dengan model matematika berbentuk meminimumkan π π , tanpa suatu kendala pada vektor π. Untuk masalah optimasi dengan model matematika tersebut, nilai optimum masalah ada pada peminimal
111
Vol. 5, No. 2, Desember 2009: 111-120
atau pemaksimal. Berikut diberikan definisi peminimal, untuk pemaksimal didefinisikan secara analog dengan tanda pertidaksamaannya adalah
.
Definisi 1. (Mital, 1979: 42) Pandang masalah meminimalkan π(π) atas πΈπ , dan diberikan π β πΈπ . Jika π(π) β€ π(π) untuk setiap π β πΈπ , maka π disebut peminimal global. Jika terdapat suatu persekitaran π, ππ (π), disekitar π sehingga π(π) β€ π(π) untuk setiap π β ππ (π), maka π disebut peminimal lokal. Jelaslah suatu peminimal global juga suatu peminimal lokal.
Syarat Perlu Keoptimalan Suatu titik π di πΈπ , dapat menjadi peminimal lokal atau peminimal global suatu fungsi π jika memenuhi karakteristik peminimal sebagai berikut.
Teorema 1. (Bazaraa and Shetty, 1990: 124) Misalkan π: πΈπ β πΈ1 terdiferensial pada π. Jika terdapat suatu vektor π
sehingga ππ(π)π π
< 0, maka terdapat suatu πΏ > 0 sehingga π π + ππ
< π(π) untuk setiap π β (0, πΏ), dan π
adalah arah descent π pada π. Bukti: Karena π terdiferensial pada π, diperoleh π π + ππ
= π π + πππ(π)π π
+ π π
πΌ π + ππ
untuk πΌ π + ππ
β 0 pada π β 0. Persamaan di atas dapat dinyatakan menjadi π π + ππ
β π π = π ππ π π π
+ π
πΌ π + ππ
. Membagi persamaan dengan π, diperoleh π π + ππ
β π π = ππ π π π
+ π
πΌ π + ππ
. π Karena ππ(π)π π
< 0 dan πΌ π + ππ
β 0 pada π β 0, terdapat suatu πΏ > 0 sehingga ππ π π π
+ π
πΌ π + ππ
< 0 untuk semua π β (0, πΏ). Dengan kata lain terdapat suatu πΏ > 0 sehingga π π + ππ
< π π untuk setiap π β (0, πΏ). qed. Sebagai akibat dari Teorema 1 setiap π peminimal lokal mengakibatkan derivatif pertama fungsi terhadap π adalah vektor nol, sebagai berikut:
112
Syarat Fritz John ... (Caturiyati)
Corollary 1. (Bazaraa and Shetty, 1990: 125) Misalkan π: πΈπ β πΈ1 terdiferensial pada π. Jika π peminimal lokal maka ππ(π) = π. Bukti: ππ(π) β π
Misalkan
ππ(π)π π
= ππ π
π
dan
diberikan
βππ(π) = β ππ(π)
2
π
= βππ(π).
Berlaku
< 0, dan dengan Teorema 1, terdapat
suatu πΏ > 0 sehingga π π + ππ
< π(π) untuk π β (0, πΏ), yaitu π pemaksimal, kontradiksi dengan asumsi bahwa π peminimal lokal. Sehingga ππ π = π. qed.
Syarat optimasi pada Teorema 1 dan corollary Teorema 1 menggunakan vektor gradien yang komponennya adalah derivatif parsial pertama dari π, disebut dengan syarat order pertama. Syarat perlu dapat juga dinyatakan dalam bentuk matriks Hessian π― dengan elemennya adalah derivatif parsial kedua dari π, disebut dengan syarat order kedua sebagai berikut:
Teorema 2. (Bazaraa and Shetty, 1990: 125) Misalkan π: πΈπ β πΈ1 terdiferensial dua kali pada π. Jika π peminimal lokal maka ππ(π) = π, dan π― π matriks semidefinit positif. Bukti: Pandang suatu arah sebarang π
. Karena π terdiferensial dua kali pada π, diperoleh π π + ππ
= π π + πππ π π π
+ 12π2 π
π π― π π
+ π2 π
2 πΌ π; ππ
1
dengan πΌ π; ππ
β 0 untuk π β 0. Karena π peminimal lokal, dari corollary Teorema 1, diperoleh ππ π = π, sehingga persamaan (1) menjadi π π + ππ
β π π = π2 (12π
π π― π π
+ π
2 πΌ π; ππ
)
1π
Persamaan (1a) dibagi dengan π2 , diperoleh π π+ππ
βπ π π2
1
= 2 π
π π― π π
+ π
2 πΌ π; ππ
.
Karena π peminimal lokal berlaku π π + ππ
β₯ π π π π+ππ
βπ π π2
2 untuk π cukup kecil, atau
β₯ 0. Sehingga berlaku 12π
π π― π π
+ π
2 πΌ π; ππ
β₯ 0 untuk π cukup kecil.
Dengan mengambil limitnya untuk π β 0, limπβ0 (12π
π π― π π
+ π
2 πΌ π; ππ
) = 1 2
π
π π― π π
maka π
π π― π π
β₯ 0, dan akibatnya π― π matriks semidefinit positif. qed.
113
Vol. 5, No. 2, Desember 2009: 111-120
Syarat Cukup Keoptimalan Berikut adalah syarat cukup optimasi suatu masalah optimasi tak berkendala, yang memberikan syarat cukup untuk peminimal lokal.
Teorema 3. (Bazaraa and Shetty, 1990: 126) Misalkan π: πΈπ β πΈ1 terdiferensial dua kali pada π. Jika ππ π = π dan π― π matriks definit positif, maka π peminimal lokal. Bukti: Karena π terdiferensial dua kali pada π, dapat diperoleh untuk setiap π β πΈπ : π π = π π + ππ π
π
1
π β π + 2 π β π ππ― π π β π
+ π β π 2 πΌ π; π β π
(3)
dengan πΌ π; π β π β 0 pada π β π. Misalkan π bukan peminimal lokal, yaitu misalkan terdapat barisan {ππ } konvergen ke π sehingga π ππ < π π untuk setiap π. Sehingga berlaku ππ π = 0 π π βπ π πβπ
2
1 πβπ π
=2
πβπ
2
notasikan π
π = ππ β π π π βπ π πβπ 2
π― π
πβπ πβπ 2
+ πΌ π; π β π
3π
ππ β π , persamaan (3a) menjadi
1
= 2 π
ππ π― π π
π + πΌ π; π β π
(3π)
dan juga berlaku π ππ < π π , yaitu π π β π π < 0 akibatnya
π π βπ π πβπ 2
< 0,
sehingga (3b) menjadi 1 2
π
ππ π― π π
π + πΌ π; ππ β π < 0
untuk setiap π.
4
Namun π
π =
ππ β π ππ β π
=
ππ β π =1 ππ β π
untuk setiap π, dan akibatnya terdapat suatu himpunan indeks πΎ sehingga {π
π }πΎ konvergen ke π
, dengan π
= 1. Dari subbarisan ini dan fakta bahwa πΌ π; π β π β 0 pada π β β β πΎ, maka (4) berakibat π
π π― π π
β€ 0. Kontradiksi dengan asumsi bahwa π― π definit positif dan fakta bahwa π
= 1. Sehingga, terbukti π₯ peminimal lokal. Pada teorema berikut ditunjukkan jika fungsi π adalah fungsi pseudokonveks.
114
Syarat Fritz John ... (Caturiyati)
Teorema 4. (Bazaraa and Shetty, 1990: 126) Diberikan π: πΈπ β πΈ1 adalah pseudokonveks pada π. π peminimal global jika dan hanya jika ππ π = π. Bukti: (β) Dari corollary Teorema 1, jika π peminimal global, maka ππ π = π. (β) Misalkan ππ π = 0, maka ππ π
π
π β π = 0 untuk setiap π β πΈπ . Karena π
pseudokonveks pada π, berakibat π π β₯ π π untuk setiap π β πΈπ . qed. Syarat Fritz John Pada Masalah Optimasi dengan Kendala Ketaksamaan Masalah optimasi berkendala mempunyai model matematika sebagai berikut: Meminimumkan π(π) dengan kendala π β π. Lebih lanjut π merupakan daerah layak masalah pemrograman linear berbentuk: meminimumkan π(π) dengan kendala π(π) β€ 0 dan π β π.
Syarat Geometrik Keoptimalan Pada bagian ini akan dibangun syarat perlu keoptimalan untuk masalah meminimumkan π(π) dengan kendala π β π, menggunakan kerucut arah layak yang didefinisikan berikut ini.
Definisi 2. (Bazaraa and Shetty, 1990: 127) Diberikan π himpunan tak kosong pada πΈπ dan π β ππ π. Kerucut arah layak π pada π, adalah π· = π
π
β π, dan π + ππ
β π untuk setiap π β 0, πΏ untuk suatu πΏ > 0 . Setiap vektor taknol π
β π· disebut arah layak. Jelas bahwa sedikit pergerakan pada π sepanjang vektor π
β π· akan berakibat titik layak cepat dicapai. Lebih lanjut, dari Teorema 1, jika ππ π π π
< 0, maka π
adalah suatu arah improving, yaitu mulai dari π, terjadi pergerakan kecil sepanjang π
yang akan mereduksi nilai π. Seperti ditunjukkan dalam Teorema 5, jika π peminimal lokal dan jika ππ π π π
< 0, maka π
β π·, yaitu syarat perlu keoptimalan lokal adalah setiap arah
115
Vol. 5, No. 2, Desember 2009: 111-120
improving bukan arah layak. Ilustrasinya pada Gambar 1, dengan sudut kerucut πΉ0 dan π· ditranslasi dari titik origin ke π.
Gambar 1
Teorema 5. (Bazaraa and Shetty, 1990: 128) Pandang masalah meminimumkan π(π) dengan kendala π β π, dengan π: πΈπ β πΈ1 , dan π himpunan tak kosong pada πΈπ . Misalkan π terdiferensial pada titik π β π. Jika π adalah solusi optimal lokal, maka πΉ0 β© π· = β
, dengan πΉ0 = π
ππ π π π
< 0 dan π· adalah kerucut arah layak π pada π. Bukti: Misalkan terdapat vektor π
β πΉ0 β© π·. Dengan Teorema 1, terdapat suatu πΏ1 > 0 sehingga π + ππ
< π π untuk setiap π β 0, πΏ1 .
5
Lebih lanjut, dari Definisi 2, terdapat suatu πΏ2 > 0 sehingga, π + ππ
β π untuk setiap π β 0, πΏ2
(6)
(5) dan (6) kontradiksi dengan hipotesis bahwa π solusi optimal lokal (minimum). Akibatnya πΉ0 β© π· = β
. qed. Dikhususkan daerah layak π adalah: π = π β π ππ π β€ 0, untuk π = 1, β¦ , π dengan ππ : πΈπ β πΈ1 untuk π = 1, β¦ , π, dan π adalah himpunan terbuka tak kosong dalam πΈπ . Diperoleh masalah pemrograman nonlinear dengan kendala ketaksamaan berikut: 116
Syarat Fritz John ... (Caturiyati)
Masalah P: Meminimumkan π(π), dengan kendala ππ π β€ 0, untuk π = 1, β¦ , π, π β π. Syarat perlu keoptimalan lokal pada π adalah πΉ0 β© π· = β
, dengan πΉ0 adalah setengah ruang terbuka didefinisikan dalam bentuk vektor gradien ππ π , dan π· adalah kerucut arah layak, yang tidak perlu didefinisikan dalam bentuk gradien fungsi terkait.
Teorema 6. (Bazaraa and Shetty, 1990: 129) Diberikan ππ : πΈπ β πΈ1 untuk π = 1, β¦ , π, dan π adalah himpunan terbuka tak kosong dalam πΈπ . Pandang masalah P: meminimumkan π(π) dengan kendala ππ π β€ 0, untuk π = 1, β¦ , π, dan π β π. Diberikan π adalah titik layak, dan diberikan πΌ = π ππ π = 0 . Lebih lanjut, misalkan π dan ππ untuk π β πΌ terdiferensial pada π dan ππ untuk π β πΌ kontinu pada π. Jika π solusi optimal lokal, maka πΉ0 β© πΊ0 = β
, dengan πΉ0 = π
ππ π π π
< 0 dan πΊ0 = π
πππ π π π
< 0 untuk setiap π β πΌ Bukti: Diberikan π
β πΊ0 . Karena π β π, dan π terbuka, maka terdapat suatu πΏ1 > 0 sehingga π + ππ
β π untuk setiap π β 0, πΏ1 .
(7)
Karena ππ π < 0 dan karena ππ kontinu pada π untuk π β πΌ, terdapat suatu πΏ2 > 0 sehingga, ππ π + ππ
< 0 untuk setiap π β 0, πΏ2 dan untuk π β πΌ.
(8)
Akhirnya karena π
β πΊ0 maka πππ π π π
< 0 untuk setiap π β πΌ, dan dengan Teorema 1, terdapat suatu πΏ3 > 0 sehingga, ππ π + ππ
< ππ π = 0 untuk setiap π β 0, πΏ3 dan untuk π β πΌ
(9)
Dari (7), (8), dan (9), jelas bahwa titik-titik berbentuk π + ππ
layak untuk masalah P untuk setiap π β 0, πΏ , dengan πΏ = minimum (πΏ1 , πΏ2 , πΏ3 ). Sehingga π
β π·, dengan π· adalah kerucut arah layak dari daerah layak pada π. Telah ditunjukkan bahwa π
β πΊ0 berakibat π
β π·, sehingga πΊ0 β π·. Dengan Teorema 5, karena π solusi lokal masalah P, πΉ0 β© π· = β
. Karena πΊ0 β π·, akibatnya πΉ0 β© πΊ0 = β
. qed. Terdapat beberapa kasus khusus dengan syarat perlu Teorema 6 dipenuhi secara trivial oleh titik non optimal yang mungkin.
117
Vol. 5, No. 2, Desember 2009: 111-120
Kasus-kasus yang mungkin: 1.
Jika
π
adalah
titik
layak
sehingga
ππ π = π.
Jelaslah
πΉ0 = π
ππ π π π
< 0 = β
, dan akibatnya πΉ0 β© πΊ0 = β
. Akibatnya suatu titik π dengan ππ π = π memenuhi syarat perlu keoptimalan. 2.
Jika suatu titik π dengan πππ π = π untuk setiap π β πΌ. Perhatikan masalah meminimumkan dengan kendala persamaan berikut: Meminimumkan π π , dengan kendala π π = 0 Kendala persamaan π π = 0 dapat digantikan oleh kendala ketaksamaan a. π1 π = π π β€ 0, b. π2 π = βπ π β€ 0. Misalkan π adalah suatu titik layak. Berlaku π1 π = π2 π = 0. Untuk vektor gradien π π adalah ππ π , maka ππ1 π = βππ2 π , dan akibatnya tidak terdapat vektor π
sehingga ππ1 π π π
< 0 dan ππ2 π π π
< 0. Sehingga, πΊ0 = β
, dan berakibat πΉ0 β© πΊ0 = β
.
Dengan kata lain, syarat perlu dari Teorema 6 dipenuhi oleh semua solusi layak. Teorema 7. (Syarat Fritz John) (Bazaraa and Shetty, 1990: 133) Diberikan π suatu himpunan terbuka tak kosong dalam πΈπ , dan diberikan π: πΈπ β πΈ1 , dan ππ : πΈπ β πΈ1 untuk π = 1, β¦ , π. Pandang masalah P untuk meminimumkan π(π) dengan kendala π β π dan ππ (π) β€ 0 untuk π = 1, β¦ , π. Diberikan π adalah solusi layak, dan diberikan πΌ = π ππ π = 0 . Lebih lanjut, misalkan π dan ππ untuk π β πΌ terdiferensial pada π dan bahwa ππ untuk π β πΌ kontinu pada π. Jika π solusi lokal masalah P, maka terdapat π’0 dan π’π untuk π β πΌ, sehingga π’0 ππ π +
πβπΌ π’π πππ
π = π, π’0 , π’π β₯ 0 untuk π β πΌ, π’0 , ππΌ β (0, π)
dengan ππΌ adalah vektor yang komponennya adalah π’π untuk π β πΌ. Lebih lanjut, jika ππ untuk π β πΌ juga terdiferensial pada π, maka syarat Fritz John dapat ditulis dalam bentuk yang ekuivalen berikut ini: π’0 ππ π +
π π=1 π’π πππ
π = π, π’π ππ π = 0 untuk π = 1, β¦ , π,
π’0 , π’π β₯ 0 untuk π = 1, β¦ , π, π’0 , π β (0, π) dengan π adalah vektor yang komponennya adalah π’π β₯ 0 untuk π = 1, β¦ , π.
118
Syarat Fritz John ... (Caturiyati)
Bukti: Karena π solusi lokal masalah P, maka dengan Teorema 6, tidak terdapat vektor π
sehingga ππ π π π
< 0 dan πππ π π π
< 0 untuk setiap π β πΌ. Sekarang, misalkan π¨ adalah matriks yang baris-barisnya adalah ππ π
π
dan πππ π
π
untuk setiap π β πΌ.
Syarat optimalitas dari Teorema 6 ekuivalen dengan pernyataan bahwa sistem π¨π
< π tak konsisten. Akibatnya terdapat vektor tak nol π β₯ π sehingga π¨π π = π. Notasikan komponen dari π dengan π’0 dan π’π untuk π β πΌ, sehingga bagian pertama teorema dipenuhi. Bentuk ekuivalen dari syarat perlu teorema diperoleh dengan memisalkan π’π = 0 untuk π β πΌ. qed. Dalam syarat Fritz John, skalar π’0 dan π’π untuk π = 1, β¦ , π sering disebut dengan pengali Lagrange. Syarat π’π ππ π = 0 untuk π = 1, β¦ , π disebut syarat komplemen slack. Ini menghasilkan π’π = 0 jika ππ π < 0. π’π > 0 hanya untuk kendala batas. Syarat Fritz John dapat juga ditulis dalam notasi vektor sebagai berikut: π’0 ππ π + ππ π π = π, ππ π π = 0, π’0 , π β₯ (0, π), π’0 , π β (0, π) dengan ππ π adalah matriks π Γ π yang kolom ke- π nya adalah πππ π , dan π adalah π-vektor pengali Lagrange. Seperti pada kasus Teorema 6, terdapat titik-titik yang memenuhi syarat Fritz John secara trivial. Kasus-kasus yang mungkin: 1.
Jika suatu titik π memenuhi ππ π = π. Dimisalkan pengali Lagrange yang terkait adalah suatu bilangan positif, kemudian himpun semua pengali yang sama dengan nol, dan memenuhi syarat Teorema 7.
2.
πππ π = π untuk suatu π β πΌ. Syarat Fritz John dari Teorema 7 juga dipenuhi benar secara trivial pada setiap titik layak untuk masalah dengan kendala persamaan yang digantikan oleh dua ketaksamaan. Secara khusus, jika π π = 0 digantikan oleh π π β€ 0 dan βπ π β€ 0, maka syarat Fritz John benar dengan π’1 = π’2 = πΌ dan menghimpun semua pengali yang sama dengan nol, dengan πΌ adalah skalar positif.
119
Vol. 5, No. 2, Desember 2009: 111-120
SIMPULAN 1.
Untuk menentukan optimasi masalah optimasi tak berkendala, digunakan vektor gradien dengan komponennya adalah derivatif pertama fungsi atau matriks Hessian yang komponennya adalah derivatif kedua fungsi.
2.
Syarat perlu keoptimalan masalah optimasi tak berkendala adalah bila dapat dicari vektor arah descent π
sehingga π π + ππ
< π(π) untuk setiap π β (0, πΏ) untuk πΏ > 0.
3.
Syarat cukup keoptimalan masalah optimasi tak berkendala adalah bila dipenuhi ππ(π) = π, dan π― π matriks semidefinit positif, maka π peminimal lokal.
4.
Untuk menentukan optimasi masalah optimasi berkendala, digunakan kerucut arah layak.
5.
Syarat perlu keoptimalan masalah optimasi berkendala diperoleh jika π solusi optimal lokal, maka πΉ0 β© πΊ0 = β
, dengan πΉ0 = π
ππ π π π
< 0
dan πΊ0 =
π
πππ π π π
< 0 untuk setiap π β πΌ . 6.
Syarat Fritz John untuk keoptimalan masalah optimasi berkendala yaitu jika π solusi lokal masalah P, maka terdapat π’0 dan π’π untuk π β πΌ, sehingga π’0 ππ π +
πβπΌ π’π πππ
π = π, π’0 , π’π β₯ 0 untuk π β πΌ, π’0 , ππΌ β (0, π)
dengan ππΌ adalah vektor yang komponennya adalah π’π untuk π β πΌ.
DAFTAR PUSTAKA Bazaraa, M. S. & Shetty, C.M. 1990. Nonlinear Programming: Theory and Algorithms . Toronto: John Wiley and Sons Inc. Mital, K.V. 1979. Optimization Methods: In Operations Research And Systems Analysis. New Delhi : Wiley Eastern Limited.
120