Pemrograman tujuan

Artikel ini perlu diwikifikasi agar memenuhi standar kualitas Wikipedia. Anda dapat memberikan bantuan berupa penambahan pranala dalam, atau dengan merapikan tata letak dari artikel ini. Untuk keterangan lebih lanjut, klik [tampil] di bagian kanan. Mengganti markah HTML dengan markah wiki bila dimungkinkan. Tambahkan pranala wiki. Bila dirasa perlu, buatlah pautan ke artikel wiki lainnya dengan cara menambahkan "[[" dan "]]" pada kata yang bersangkutan (lihat WP:LINK untuk keterangan lebih lanjut). Mohon jangan memasang pranala pada kata yang sudah diketahui secara umum oleh para pembaca, seperti profesi, istilah geografi umum, dan perkakas sehari-hari. Sunting bagian pembuka. Buat atau kembangkan bagian pembuka dari artikel ini. Susun header artikel ini sesuai dengan pedoman tata letak. Tambahkan kotak info bila jenis artikel memungkinkan. Hapus tag/templat ini.

Pemrograman tujuan (bahasa Inggris: Goal programming; disingkat dengan GP), adalah sebuah metode optimisasi multiobjektif dalam bidang analisis keputusan multi-kriteria (multi-criteria decision analysis, MCDA). Metode ini dapat dianggap sebagai perumuman dari pemrograman linear, agar dapat menangani beberapa objektif secara serentak yang umumnya saling berkonflik. Contoh permasalahan ini dapat berupa mencari solusi yang dapat meningkatkan keuntungan bersih sekaligus mengurangi besarnya ongkos usaha. Menggunakan informasi target nilai yang ingin dicapai setiap objektif, GP melakukan optimisasi dengan meminimumkan besar penyimpangan-penyimpangan yang tidak diinginkan solusi dengan nilai target-target tersebut. Fungsi objektif dalam GP disebut dengan fungsi pencapaian, dan digunakan untuk menghampiri preferensi pengambil keputusan terkait caranya meminimumkan penyimpangan-penyimpangan tersebut. Fungsi ini dapat bernilai skalar maupun vektor, tergantung varian GP yang digunakan.

Model GP menggunakan filosofi satisficing ketimbang optimising dalam pengambilan solusi.^[1] GP dapat digunakan untuk melakukan analisis terkait: daftar sumber daya yang diperlukan untuk memenuhi suatu kumpulan objektif, besarnya ketercapaian (atau kompromi) setiap objektif yang dihasilkan suatu kombinasi sumber daya, dan solusi dominan dari berbagai kendala dan tingkat kepentingan setiap objektif. Saat ini, GP telah diterapkan dalam pemilihan portofolio,^[2] perencanaan agrikultur,^[3] penjadwalan,^[4] dan banyak lagi.

Formulasi model GP paling awal diperkenalkan oleh Charnes et al. dalam permasalahan menentukan kompensasi untuk manajer.^[5] Walaupun pada saat itu model dianggap sebagai adaptasi dari pemrograman linear (LP), formulasi tersebut menampilkan konsep "regresi terbatas" yang menjadi dasar minimisasi dalam GP.^[6] Istilah goal programming baru muncul pada buku teks oleh Charnes dan Cooper pada tahun 1961 -- walau model masih disajikan sebagai perumuman pemrograman linear (LP) untuk menyelesaikan masalah LP yang infeasible (tidak memiliki solusi).^[7] Pengembangan model GP lebih lanjut yang dilakukan oleh oleh Ijiri,^[8] Lee,^[9] dan Ignizio,^[10] membuat model ini menjadi umum digunakan sebagai alat riset operasi.

Pada perkembangannya, salah satu aspek yang membuat model GP menjadi pilihan metode MCDA/MCDM yang umum digunakan, adalah faktor formulasi yang sederhana dan kemiripan dengan model LP.^[11] Dalam masa ini pula, berbagai varian GP juga dikembangkan untuk menyelesaikan masalah dalam banyak bidang. Akan tetapi banyak kritik disampaikan terkait GP pada tahun 1980-an, yang oleh sebagian pihak dianggap terjadi karena kesalahan mendasar yang disebabkan oleh kurangnya pemahaman praktek MCDA/MCDM.^[12][13] Perdebatan terkait fundamental model GP selanjutnya menghasilkan publikasi buku teks oleh Romero,^[14] yang menunjukkan masalah yang terjadi lebih diakibatkan oleh praktik pemodelan yang buruk, dan memberikan praktik-praktik GP yang baik.^[15]

Schniederjans merangkum perkembangan GP sampai tahun 1995, dan mencatat bibliografi banyak artikel yang berkaitan dengan GP.^[16] Selain itu, Jones dan Tamiz memberikan bibliografi serta deskripsi untuk artikel periode 1990-2000.^[17] Sebuah buku teks tahun 2010 oleh Jones dan Tamiz berisi penjelasan komprehensif terkait metode GP yang terbaru.^[18]

Karena model GP dibuat dari adaptasi model pemrograman linear (LP), perumusan model GP dapat dihasilkan dari mengubah notasi dan asumsi-asumsi model LP. Pada model LP, rumusan model dapat dinyatakan dalam bentuk kanonik sebagai

$${\displaystyle {\begin{aligned}\min &\,Z=\sum _{j=1}^{n}c_{j}x_{j}\\{\text{dengan kendala}}&\,f_{i}(\mathbf {x} )\geq b_{i};\,{\text{untuk }}i=1,\,\dots ,\,m\\&\,\quad \,x_{j}\geq 0;\,{\text{untuk }}j=1,\,\dots ,\,n.\end{aligned}}}$$

Pada bentuk di atas, ${\displaystyle m}$ menyatakan banyaknya kendala, ${\displaystyle x_{j}}$disebut variabel keputusan, dan ${\displaystyle Z}$ disebut fungsi objektif. Fungsi linear$${\displaystyle f_{i}(\mathbf {x} )=\sum _{j=1}^{n}a_{ij}x_{j}}$$menyatakan persamaan kendala ke-${\displaystyle i}$, yang dapat bernilai lebih besar (memiliki deviasi positif) daripada nilai target ${\displaystyle b_{i}}$. Dengan bantuan aritmetika dasar, model LP tersebut juga dapat menyertakan kendala dalam bentuk "kurang dari" dan "sama dengan." Akan tetapi, terlepas dari jenis kendala yang digunakan, model LP hanya mempertimbangkan solusi yang memenuhi semua kendala. Ketika kendala-kendala dalam model LP tidak semuanya dapat saling dipenuhi, model LP disebut infeasible dan tidak memiliki solusi. Agar tetap menghasilkan solusi yang wajar dalam kondisi tersebut, Charnes dan Cooper menganggap masing-masing kendala^{[note 1]} yang menyusun model LP sebagai sebuah fungsi:^[7]$${\displaystyle h_{i}(\mathbf {x} )={\Big |}f_{i}(\mathbf {x} )-b_{i}{\Big |};\;\;{\text{untuk }}i=1,\,\dots ,\,Q.}$$Fungsi-fungsi ini dianggap sebagai sebuah tujuan (goal) yang perlu dipenuhi agar menghasilkan solusi feasible. Dengan kata lain, solusi feasible dihasilkan ketika ${\displaystyle h_{i}(\mathbf {x} )=0}$. Pada keadaan hal tersebut tidak mungkin terjadi, solusi "terbaik" selanjutnya adalah solusi yang terletak sedekat mungkin dengan tujuan-tujuan yang ditetapkan, dan dapat dicari dengan melakukan optimisasi$${\displaystyle {\begin{aligned}\min &\,Z=\sum _{i=1}^{Q}(d_{i}^{+}+d_{i}^{-})\\{\text{dengan kendala}}&\,f_{i}(\mathbf {x} )-d_{i}^{+}+d_{i}^{-}=b_{i};\,{\text{untuk }}i=1,\,\dots ,\,Q\\&\,d_{i}^{+},\,d_{i}^{-}\geq 0;\,\,{\text{untuk }}i=1,\,\dots ,\,Q\\&\,x_{j}\geq 0;\,\,{\text{untuk }}j=1,\,\dots ,\,n\end{aligned}}}$$dengan ${\displaystyle d_{i}^{+}}$ disebut variabel deviasi positif dan ${\displaystyle d_{i}^{-}}$ disebut variabel deviasi negatif. Variabel ${\displaystyle Q}$ digunakan untuk menyatakan banyaknya tujuan yang dipertimbangkan model. Bentuk di atas selanjutnya menjadi dasar dari model GP.^[19] Kendala-kendala pada model GP yang melibatkan variabel-variabel deviasi disebut dengan kendala lunak (soft constraint, juga disebut goal constraint). Istilah 'kendala lunak' digunakan untuk membedakannya dengan kendala pada model LP, yang selanjutnya dirujuk sebagai kendala tegas (hard constraint). Kendala tegas dapat disertakan dalam model GP namun disarankan tidak dalam jumlah yang banyak. Penggunaan kendala tegas yang berlebihan dapat meniadakan solusi-solusi alternatif yang mungkin menarik bagi pengambil keputusan. Tanpa mengurangi keumuman, semua kendala pada model GP dalam artikel ini dianggap sebagai kendala lunak -- dengan kata lain, ${\displaystyle Q:=m}$.

Lebih lanjut, fungsi ${\displaystyle Z}$ pada bentuk di atas disebut dengan fungsi pencapaian, untuk membedakannya dengan model LP. Fungsi ini hanya menyertakan variabel-variabel deviasi yang ingin diminimumkan dan tidak melibatkan variabel-variabel keputusan.^{[note 2]} Variabel-variabel deviasi dapat dihasilkan dari ${\displaystyle f_{i}(\mathbf {x} ),}$ menggunakan hubungan

$${\displaystyle d_{i}^{+}={\frac {1}{2}}{\Big [}|h_{i}(\mathbf {x} )|+h_{i}(\mathbf {x} ){\Big ]}={\begin{cases}h_{i}(\mathbf {x} )-b_{i}&{\text{jika }}h_{i}(\mathbf {x} )\geq b_{i}\\0&{\text{lainnya.}}\\\end{cases}}}$$dan

$${\displaystyle d_{i}^{-}={\frac {1}{2}}{\Big [}|h_{i}(\mathbf {x} )|-h_{i}(\mathbf {x} ){\Big ]}={\begin{cases}b_{i}-h_{i}(\mathbf {x} )&{\text{jika }}h_{i}(\mathbf {x} )\leq b_{i}\\0&{\text{lainnya.}}\\\end{cases}}}$$

Pendefinisian variabel-variabel deviasi di atas mengartikan sebuah kendala tidak dapat mengalami pencapaian melebihi target (nilai deviasi positif ${\displaystyle d_{i}^{+}>0}$) sekaligus pencapaian dibawah target (nilai deviasi ${\displaystyle d_{i}^{-}>0}$) secara bersamaan.

Model GP dapat diklasifikasikan ke dalam beberapa varian, yang masing-masing dicirikan oleh metrik jarak yang mendasari atau fungsi utilitas yang berbeda. Tiga varian umum dari model GP akan dideskripsikan di sini.^[17]

Pada varian model GP yang pertama, fungsi pencapaian akan mengukur jumlah berbobot besar penyimpangan yang tidak diinginkan dari nilai targetnya. Bobot tersebut disesuai dengan ukuran kepentingan relatifnya bagi pengambil keputusan. Metrik jarak yang mendasari model ini adalah jarak Manhattan. Varian ini dikenal sebagai GP berbobot (weighted GP atau non-preemptive GP) dan selanjutnya disingkat dengan WGP. Formulasi WGP dengan ${\displaystyle Q}$ tujuan dapat dinyatakan sebagai

$${\displaystyle {\begin{aligned}\min &\,Z=\sum _{i=1}^{Q}{\frac {1}{k_{i}}}{\Big (}w_{i}^{+}d_{i}^{+}+w_{i}^{-}d_{i}^{-}{\Big )}\\{\text{dengan kendala}}&\,f_{i}(\mathbf {x} )-d_{i}^{+}+d_{i}^{-}=b_{i};\,{\text{untuk }}i=1,\,\dots ,\,Q\\&\,d_{i}^{+},\,d_{i}^{-}\geq 0;\,\,{\text{untuk }}i=1,\,\dots ,\,Q\\&\,x_{j}\geq 0;\,\,{\text{untuk }}j=1,\,\dots ,\,n\end{aligned}}}$$

Konstanta ${\displaystyle w_{i}^{+}}$ dan ${\displaystyle w_{i}^{-}}$ masing-masing menyatakan bobot (bernilai non-negatif) untuk variabel-variabel deviasi dalam fungsi pencapaian. Bobot akan diberi nilai nol jika pengambil keputusan menganggap variabel deviasi yang bersangkutan tidak penting untuk diminimumkan. Setiap variabel deviasi dalam fungsi pencapaian selanjutnya perlu dinormalisasi agar memiliki satuan ukuran yang sama. Beberapa pilihan normalisasi dapat dilihat di bagian Normalisasi.

Model WGP digunakan ketika semua tujuan dapat dibandingkan secara langsung, dan pengambil keputusan bersedia dan mampu memberikan bobot yang mencerminkan kepentingan relatif dari tujuan-tujuan dalam situasi tersebut. Selain itu, WGP sebaiknya digunakan ketika pengambil keputusan tertarik pada solusi yang memberikan total deviasi terbobot terkecil dari tujuan-tujuan, dan bukan keseimbangan dari keseluruhan pencapaian tujuan-tujuan tersebut. Dalam situasi seperti ini, WGP merupakan alat yang ampuh yang tidak hanya memberikan solusi, namun juga memberikan informasi mengenai tarik-ulur (trade-off) antar tujuan-tujuan.

Varian umum lain dari GP dihasilkan ketika variabel deviasi ditetapkan ke dalam sejumlah tingkat prioritas, lalu diminimumkan secara berurutan (leksikografis). Lebih spesifik, minimisasi dilakukan bertahap dari prioritas tertinggi ke terendah, sambil mempertahankan nilai minimum yang dicapai oleh semua minimisasi pada tahap-tahap sebelumnya. Varian GP ini dikenal sebagai model GP leksikografik (lexicographic GP atau pre-emptive GP) dan selanjutnya akan dirujuk dengan LGP. Varian LGP diperkenalkan dan dikembangkan (terutama) oleh Ijiri,^[8] Lee,^[9] dan Ignizio.^[20]

Model LGP dapat dituliskan sebagai berikut:

$${\displaystyle {\begin{aligned}{\text{lex min}}&\,\mathbf {a} ={\Big [}g_{1}(\mathbf {d} ^{+},\,\mathbf {d} ^{-}),\,\dots ,g_{L}(\mathbf {d} ^{+},\,\mathbf {d} ^{-}){\Big ]}\\{\text{dengan kendala}}&\,f_{i}(\mathbf {x} )-d_{i}^{+}+d_{i}^{-}=b_{i};\,{\text{untuk }}i=1,\,\dots ,\,Q\\&\,d_{i}^{+},\,d_{i}^{-}\geq 0;\,\,{\text{untuk }}i=1,\,\dots ,\,Q\\&\,x_{j}\geq 0;\,\,{\text{untuk }}j=1,\,\dots ,\,n\end{aligned}}}$$

Model ini memiliki ${\displaystyle L}$ tingkat prioritas, dan ${\displaystyle Q}$ tujuan. Fungsi pencapaian ${\displaystyle \mathbf {a} }$ adalah fungsi bernilai vektor terurut dari ${\displaystyle L}$ tingkat prioritas tersebut. Setiap kendala tegas ditempatkan, berdasarkan konvensi, pada tingkat prioritas pertama. Fungsi ${\displaystyle g}$ umumnya (dalam setiap tingkat prioritas) memiliki bentuk

$${\displaystyle g_{l}(\mathbf {d} ^{+},\,\mathbf {d} ^{-})=\sum _{i=1}^{Q}{\frac {1}{k_{i}}}{\Big (}w_{il}^{+}d_{i}^{+}+w_{il}^{-}d_{i}^{-}{\Big )},}$$

dengan ${\displaystyle w_{il}^{+}}$ dan ${\displaystyle w_{il}^{-}}$ menyatakan bobot relatif intra-prioritas. Seperti WGP, bobot bernilai nol diberikan pada variabel-variabel deviasi yang tidak penting untuk diminimumkan. Ketika variabel-variabel deviasi pada suatu tingkat prioritas memiliki satuan pengukuran yang berbeda, teknik-teknik normalisasi WGP dapat diterapkan agar aritmatika wajar dilakukan.

Model LGP adalah varian yang paling banyak diperdebatkan dalam aspek fungsi utilitas yang mendasarinya. Romero^[14] memberikan diskusi tentang topik ini serta praktik-praktik pemodelan ketika menggunakan pemrograman tujuan leksikografis. Varian ini digunakan ketika pengambil keputusan memiliki tujuan-tujuan yang dapat diurutkan secara alami, dan bukan yang berbentuk perbandingan relatif. Hal ini juga digunakan ketika pengambil keputusan tidak dapat, atau tidak mau memberikan kepentingan relatif yang relevan, dari tujuan-tujuan dengan menggunakan bobot. Pemrograman tujuan leksikografis secara historis merupakan varian pemrograman tujuan yang paling banyak digunakan.^[21]

Varian lain yang kurang banyak digunakan tetapi secara teoritis signifikan adalah varian GP Chebyshev, juga dikenal sebagai GP Minmax. Dalam varian ini, model akan memimumkan nilai deviasi terbesar dari kumpulan deviasi-deviasi, bukan jumlah dari deviasi-deviasi itu sendiri. Formulasi aljabar dari model GP Chebyshev diberikan sebagai:

$${\displaystyle {\begin{aligned}\min &\,\lambda \\{\text{dengan kendala}}&\,\sum _{i=1}^{Q}{\frac {1}{k_{i}}}{\Big (}w_{i}^{+}d_{i}^{+}+w_{i}^{-}d_{i}^{-}{\Big )}\leq \lambda ;{\text{untuk }}i=1,\,\dots ,\,Q\\&\,f_{i}(\mathbf {x} )-d_{i}^{+}+d_{i}^{-}=b_{i};\,{\text{untuk }}i=1,\,\dots ,\,Q\\&\,d_{i}^{+},\,d_{i}^{-}\geq 0;\,\,{\text{untuk }}i=1,\,\dots ,\,Q\\&\,x_{j}\geq 0;\,\,{\text{untuk }}j=1,\,\dots ,\,n\end{aligned}}}$$

Model tersebut mempertimbangkan ${\displaystyle Q}$ tujuan. Konstanta bobot relatif ${\displaystyle w_{i}^{+}}$ dan ${\displaystyle w_{i}^{-}}$, dan normalisasi ${\displaystyle k_{i}}$, memiliki sifat yang sama dengan yang telah disampaikan pada varian WGP. Akan tetapi, berbeda dengan WGP, GP Chebyshev meminimumkan nilai deviasi terbesar, ketimbang total nilai deviasi. Akibatnya, GP Chebyshev sebisa mungkin akan menghasilkan keseimbangan antar besar penyimpangan pada setiap tujuan. Pengambil keputusan sebaiknya menggunakan GP Chebyshev jika kendala-kendala mereka didefinisikan dalam aspek keseimbangan dan keadilan.

Dalam proses penyelesaian model GP, peran variabel-variabel deviasi dalam fungsi pencapaian perlu dinormalisasi. Normalisasi menyebabkan variabel-variabel deviasi memiliki satuan dan magnitudo yang sama, sehingga perhitungan antar variabel-variabel deviasi masuk akal dilakukan dan hasilnya dapat dibandingkan. Akan tetapi, ada banyak pilihan metode normalisasi beserta untung-rugi dalam menggunakan mereka; pilihan yang sesuai tergantung pada situasi masalah yang terjadi dan preferensi pengambil keputusan.

Metode normalisasi pada bagian ini difokuskan pada varian GP berbobot. Walaupun demikian, konsep setiap normalisasi dapat diterapkan atau disesuaikan untuk banyak varian model GP lainnya. Menggunakan masalah WGP di bawah ini sebagai contoh

$${\displaystyle {\begin{aligned}\min &\,d_{1}^{+}+d_{2}^{-}+d_{3}^{-}+d_{4}^{-}\\{\text{dengan kendala}}&\,3x_{1}+4x_{2}-d_{1}^{+}+d_{1}^{-}=560\\&\,100x_{1}+150x_{2}-d_{2}^{+}+d_{2}^{-}=7800\\&\,x_{1}-d_{3}^{+}+d_{3}^{-}=80\\&\,x_{2}-d_{4}^{+}+d_{4}^{-}=90\\&\,x_{1}+x_{2}\leq 160\\&\,d_{i}^{+},\,d_{i}^{-}\geq 0;{\text{untuk }}i=1,\,\dots ,\,m\\&\,x_{j}\geq 0;{\text{untuk }}j=1,\,\dots ,\,n\end{aligned}}}$$

Berikut adalah beberapa teknik normalisasi yang umum dilakukan.^[22][23]

Normalisasi ini akan menskala setiap variabel deviasi sehingga kontribusinya di fungsi pencapaian menyatakan persen penyimpangan dari nilai target yang ditetapkan. Secara matematis, fungsi pencapaian pada model di atas akan berubah menjadi

$${\displaystyle \min {\Bigg (}{\frac {d_{1}^{+}}{560}}+{\frac {d_{2}^{-}}{7800}}+{\frac {d_{3}^{-}}{80}}+{\frac {d_{4}^{-}}{90}}{\Bigg )}}$$

Normalisasi jenis ini mudah diterapkan dan konsep total persentase penyimpangan sederhana untuk diartikan. Akan tetapi, metode ini tidak cocok jika pengambil keputusan ingin membandingkan secara langsung beberapa tujuan yang diukur dalam unit satuan yang sama. Sebagai contoh, misalkan masalah GP dengan dua tujuan berikut yang dinyatakan dalam satuan rupiah:

$${\displaystyle {\begin{aligned}x_{1}+3x_{2}-d_{1}^{+}+d_{1}^{-}&=1000\\2x_{1}+x_{2}-d_{1}^{+}+d_{2}^{-}&=500\end{aligned}}}$$

Menerapkan normalisasi persentase akan menghasilkan fungsi pencapaian$${\displaystyle \min {\Bigg (}{\frac {d_{1}^{-}}{1000}}+{\frac {d_{2}^{-}}{500}}{\Bigg )}}$$yang mengartikan penyimpangan seribu rupiah dari tujuan pertama, setara dengan setengah penyimpangan seribu rupiah dari tujuan kedua. Sudut pandang ini benar jika persentase penyimpangan yang ingin dibandingkan, namun tidak sesuai jika penyimpangan kedua tujuan dibandingkan dalam satuan rupiah. Beberapa masalah lain terkait normalisasi adalah formulasinya yang tidak kokoh (robust), karena perubahan nilai target akan menghasilkan persentase penyimpangan yang berbeda (walaupun besar penyimpangannya sama); tidak dapat digunakan ketika target bernilai nol; dan perlu modifikasi pada fungsi pencapaian untuk target-target bernilai negatif.^[24]

Normalisasi nol-satu akan menskala setiap variabel deviasi agar bernilai diantara 0 sampai 1 dalam fungsi pencapaian. Dengan kata lain, normalisasi ini mengukur penyimpangan terburuk yang dapat terjadi untuk setiap variabel deviasi. Suatu nilai batas yang realistis dapat ditetapkan bagi variabel deviasi yang unbounded (tidak terbatas). Fungsi pencapaian pada model akan berubah menjadi

$${\displaystyle \min {\Bigg (}{\frac {d_{1}^{+}}{\max d_{1}^{+}}}+{\frac {d_{2}^{-}}{\max d_{2}^{-}}}+{\frac {d_{3}^{-}}{\max d_{3}^{-}}}+{\frac {d_{4}^{-}}{\max d_{4}^{-}}}{\Bigg )},}$$

dengan masing-masing nilai ${\displaystyle \max }$ didapatkan dari optimisasi satu-objektif (model LP standar) dari kendala-kendala model GP di atas. Variasi lain dari normalisasi ini adalah membagi setiap variabel deviasi dengan rentang nilainya. Hal ini mengakibatkan fungsi pencapaian memiliki bentuk

$${\displaystyle \min {\Bigg (}{\frac {d_{1}^{+}}{\max d_{1}^{+}-\min d_{1}^{+}}}+{\frac {d_{2}^{-}}{\max d_{2}^{-}-\min d_{2}^{-}}}+{\frac {d_{3}^{-}}{\max d_{3}^{-}-\min d_{3}^{-}}}+{\frac {d_{4}^{-}}{\max d_{4}^{-}-\min d_{4}^{-}}}{\Bigg )},}$$

Metode ini cocok pada kasus permasalahan dengan setiap objektif memiliki interval nilai yang jelas, dan setiap solusi pada ruang solusi ingin diperhatikan pengambil keputusan. Tapi, tujuan yang tidak terbatas (tidak memiliki batas penyimpangan terburuk), atau masalah dengan banyak ruang solusi yang unbounded, dapat menghasilkan solusi-solusi yang tidak relevan ketika menggunakan normalisasi ini. Disamping itu, normalisasi nol-satu memerlukan optimisasi satu-objektif sebanyak variabel deviasi yang disertakan dalam fungsi pencapaian. Hal ini mungkin tidak praktis pada masalah yang kompleks dengan waktu komputasi yang lama.

Pada beberapa kasus, model GP perlu mengukur jarak geometris antara nilai suatu tujuan dengan targetnya, ketimbang mengukur besar deviasi antara keduanya. Sebagai contoh, titik ${\displaystyle (20,\,20)}$ pada tujuan pertama, akan memiliki jarak geometris sebesar ${\displaystyle 12}$ dari target ${\displaystyle 560}$. Perubahan koefisien normalisasi dalam tujuan ini tidak akan mengubah jarak geometris tersebut. Oleh karena itu, normalisasi perlu dilakukan agar setiap variabel deviasi dapat mengukur jarak Euklides nilai tujuan dari targetnya.^[23] Untuk contoh masalah GP di atas, fungsi pencapaian pada model akan berubah menjadi$${\displaystyle \min {\Bigg (}{\frac {d_{1}^{+}}{\sqrt {4^{2}+3^{2}}}}+{\frac {d_{2}^{-}}{\sqrt {100^{2}+150^{2}}}}+{\frac {d_{3}^{-}}{\sqrt {1^{2}}}}+{\frac {d_{4}^{-}}{\sqrt {1^{2}}}}{\Bigg )}.}$$Normalisasi Euklides kokoh secara komputasi; karena dapat dilakukan untuk setiap tujuan dan nilai target, dan tidak memerlukan optimisasi atau perhitungan yang kompleks untuk mendapatkan konstanta normalisasi. Akan tetapi, dua jenis permasalahan terkait metode ini terlihat jelas dari fungsi pencapaian yang dihasilkan. Pertama, normalisasi ini tidak mempertimbangkan nilai target sehingga menghasilkan konstanta yang rendah untuk tujuan ke-3 dan ke-4. Sedangkan permasalahan kedua, berbeda dengan kedua normalisasi sebelumnya, nilai optimal yang dihasilkan oleh fungsi pencapaian tidak memiliki makna yang mudah diintepretasikan. Karena dua hal tersebut, normalisasi ini sebaiknya digunakan untuk kasus permasalahan yang tidak praktis untuk menggunakan normalisasi persentase maupun nol-satu.^[25]

Terdapat beberapa normalisasi lainnya yang dipertimbangkan dan digunakan dalam literatur. Normalisasi penjumlahan akan membagi masing-masing variabel deviasi dengan jumlah mutlak koefisien kendala lunak.^[24] Pada contoh masalah, variabel ${\displaystyle d_{2}^{-}}$ akan dinormalisasi menjadi

${\displaystyle {\frac {d_{2}^{-}}{|100|+|150|}}}$

Metode ini memiliki pembagi yang lebih besar daripada metode Euklides dan terbukti lebih baik ketika menskalakan masalah yang sangat tidak dapat dibandingkan. Metode ini sama kokohnya dengan metode Euklides, namun tidak memberikan makna berarti pada nilai fungsi pencapaian. Lebih lanjut, Tamiz dan Jones mengembangkan metode hibrida untuk menentukan jenis normalisasi -- Euklides atau penjumlahan -- yang cocok untuk setiap tujuan dalam model.^[26] Di lain pihak, Romero membahas metode yang menggabungan normalisasi persentase dan Euklides, dan konsekuensi dari penggunaannya.^[23]

Filosofi yang mendasari model GP adalah teori Simon tentang satisficing.^[27][28] Filosofi ini diterapkan dengan mencoba sebisa mungkin memenuhi sekumpulan nilai target yang ditetapkan oleh pengambil keputusan. Akan tetapi, model GP yang juga mencakup elemen-elemen dari filosofi optimising, karena dikembangkan dari sudut pandang pemrograman linear. Secara matematis, model GP merupakan sebuah optimasi tunggal atau serangkaian optimasi yang saling terkait. Menggabungkan kedua filosofi pengambilan keputusan mengakibatkan model GP dapat menghasilkan solusi tak-efisien Pareto -- yakni, solusi dengan nilai satu atau lebih objektif yang dapat dibuat lebih baik, tanpa membuat nilai objektif lain pada solusi tersebut menjadi lebih buruk.^[29]

Beberapa metode telah dikembangkan untuk mendeteksi dan memperbaiki solusi tak-efisien Pareto yang dihasilkan model GP. Sebelum metode restorasi dapat dilakukan, setiap tujuan harus dijaga agar nilainya tidak mengalami degradasi. Hal ini dapat dicapai dengan menempatkan batas atas atau bawah pada variabel deviasi. Sebagai contoh, misalkan masalah LGP

$${\displaystyle {\begin{aligned}\min &\,{\Big [}d_{1}^{-},\,5d_{2}^{-}+2d_{3}^{+}{\Big ]}\\{\text{dengan kendala}}&\,f_{1}(\mathbf {x} )-d_{1}^{+}+d_{1}^{-}=b_{1}\\&\,f_{2}(\mathbf {x} )-d_{2}^{+}+d_{2}^{-}=b_{2}\\&\,f_{3}(\mathbf {x} )-d_{3}^{+}+d_{3}^{-}=b_{3}\\&\,d_{i}^{+},\,d_{i}^{-}\geq 0;{\text{untuk }}i=1,\,\dots ,\,m\\&\,x_{j}\geq 0;{\text{untuk }}j=1,\,\dots ,\,n\end{aligned}}}$$

menghasilkan solusi (yang mungkin tak-efisien Pareto) dengan variabel-variabel deviasi bernilai

$${\displaystyle {\Big (}(d_{1}^{+})^{*},\,(d_{1}^{-})^{*},\,(d_{2}^{+})^{*},\,(d_{2}^{-})^{*},\,(d_{3}^{+})^{*},\,(d_{3}^{-})^{*}{\Big )}}$$

Nilai tujuan-tujuan yang dihasilkan model dapat dijaga dengan menambahkan kendala-kendala berikut pada kendala model tersebut

$${\displaystyle {\begin{aligned}d_{1}^{-}\leq (d_{1}^{-})^{*}\\d_{2}^{-}\leq (d_{2}^{-})^{*}\\d_{3}^{+}\leq (d_{3}^{+})^{*}\\d_{1}^{+}\geq (d_{1}^{+})^{*}\\d_{2}^{+}\geq (d_{2}^{+})^{*}\\d_{3}^{-}\geq (d_{3}^{-})^{*}\\\end{aligned}}}$$

Dengan adanya batasan-batasan ini, setiap perbaikan dalam tujuan apapun akan menjadi langkah Pareto yang dapat diterima menuju batas efisien. Mungkin terdapat banyak solusi pada batas efisien yang merupakan hasil dari proses restorasi. Metode restorasi yang baik akan memilih satu titik di antara sekumpulan solusi tersebut yang optimal sehubungan dengan preferensi DM. Ada beberapa metode restorasi yang mungkin dilakukan, dan diuraikan di bawah ini.^[30]

Pada metode ini, variabel-variabel deviasi yang tidak diperhatikan pada setiap tujuan dalam penyelesaian model GP akan dimaksimumkan.^[14][31][32] Hal ini dilakukan dengan menambahkan satu prioritas tambahan pada fungsi pencapaian LGP, atau dengan mengonversi varian WGP menjadi LGP dengan dua tingkat prioritas. Sebagai contoh, fungsi pencapaian pada contoh di atas berubah menjadi (dengan asumsi ketiga tujuan tak-efisien),

$${\displaystyle \min \,{\Big [}d_{1}^{-},\,5d_{2}^{-}+2d_{3}^{+},-(d_{1}^{+}+d_{2}^{+}+d_{3}^{-}){\Big ]}}$$

Metode ini akan memastikan solusi yang dihasilkan bersifat efisien Pareto, asalkan model masalah awal terbatas Pareto -- tidak ada tujuan yang selalu dapat dibuat menjadi lebih baik. Akan tetapi, restorasi ini tidak mempertimbangkan bobot relatif pada setiap variabel deviasi.

Pada metode ini, variabel-variabel deviasi yang tidak diperhatikan akan dimaksimumkan, dengan cara yang berbeda dengan restorasi langsung. Restorasi ini menggunakan kembali bobot-bobot relatif dan/atau tingkat prioritas yang digunakan model awal.^[33] Sebagai contoh, restorasi berbasis preferensi untuk pencapaian pada contoh di atas adalah

$${\displaystyle \min \,{\Big [}d_{1}^{-},\,5d_{2}^{-}+2d_{3}^{+},-(d_{1}^{+}),\,-(5d_{2}^{+}+2d_{3}^{-}){\Big ]}}$$

Metode ini memastikan kontinuitas preferensi, baik di-bawah maupun di-atas nilai target, sehingga memastikan solusi efisien Pareto yang ditemukan sesuai dengan preferensi awal pengambil keputusan.

Dalam metode ini, serangkaian tujuan yang tidak efisien disajikan kepada pengambil keputusan. Selanjutnya mereka diminta untuk memilih tujuan yang paling mendesak untuk diperbaiki. Proses ini berlanjut secara berulang hingga solusi yang efisien dihasilkan.^[33] Restorasi interaktif mempertimbangkan fakta bahwa preferensi pengambil keputusan dapat tidak kontinu. Hal ini sering terjadi dalam penerapan model GP, sebagai contoh ketika beberapa tujuan mewakili standar yang harus dipenuhi (persyaratan hukum, standar lingkungan, dll.) sementara yang lain mewakili keinginan pengambil keputusan (biaya rendah, keuntungan tinggi, dll.). Metode restorasi interaktif memiliki keuntungan dalam melibatkan pengambil keputusan dalam proses restorasi. Sebagai akibatnya, restorasi ini juga termasuk dalam metode multi-objektif interaktif.^[34] Untungnya, jumlah iterasi interaktif yang diperlukan untuk memperbaiki efisiensi solusi umumnya kecil.^[33]

Dalam model GP, sifat matematis dari tujuan dan/atau variabel keputusan, maupun hubungan antar keduanya, dapat diperumum menyesuaikan masalah nyata yang diteliti. Model yang dihasilkan ini berbeda dengan model-model dalam bagian Varian, yang mencatat variasi metrik dan fungsi utilitas model yang umum. Pada kenyataannya, banyak implementasi model GP disusun dari gabungan dari dua atau lebih variasi, dan dapat memiliki beberapa nama berbeda. Berikut adalah beberapa perumuman model GP secara garis besar.^[17][35][36]

GP interval adalah salah satu cara tertua untuk memperbesar struktur preferensi dalam model GP. Charnes dan Collomb (1972) mengusulkan sebuah model untuk memperlonggar syarat setiap nilai target dapat didefinisikan dengan pasti.^[37] Dalam model ini, mereka memungkinkan pengambil keputusan untuk menentukan interval kepuasan dan menetapkan penalti untuk nilai deviasi yang tidak termasuk dalam interval tersebut.

Metode ini berdasar pada LP interval, dengan suatu nilai batas bawah ${\displaystyle _{L}b}$ dan suatu nilai batas atas ${\displaystyle _{U}b}$ menjadi batas interval untuk nilai ekspresi di ruas-kanan; secara matematis:

$${\displaystyle {_{L}b_{i}}\leq f_{i}(\mathbf {x} )\leq {_{U}b_{i}}.}$$

Bentuk GP interval yang bersesuaian dapat dihasilkan dengan menggunakan dua pasang variabel deviasi:

$${\displaystyle {\begin{aligned}f_{i}(\mathbf {x} )-{_{U}d_{i}}^{+}+{_{U}d_{i}}^{-}&={_{U}b_{i}}\\f_{i}(\mathbf {x} )-\,{_{L}d_{i}}^{+}+\,{_{L}d_{i}}^{-}&=\,{_{L}b_{i}}.\\\end{aligned}}\,}$$

Variabel deviasi ${\displaystyle {_{U}d_{i}}^{+}}$ dan ${\displaystyle {_{L}d_{i}}^{-}}$ perlu diminimumkan dalam fungsi pencapaian, sedangkan variabel deviasi yang lain diperbolehkan mengambil nilai yang dihasilkan dari hasil optimasi.^[38] GP Interval dapat digunakan untuk menyelesaikan banyak masalah formulasi yang umumnya digunakan untuk mengkritik model-model GP, seperti penetapan besar nilai (yang spefisik) target untuk setiap tujuan.^[39]

Dalam aspek teoritis, pembaca dapat membaca bukti-bukti dasar yang diberikan oleh Charnes dan Cooper.^[40] Beberapa penelitian juga telah dilakukan terkait penggunaan interval pada nilai pembobotan dan/atau target,^[41] penggunaan relasi preferensi,^[42] dan penggunaan analisis sensitivitas.^[43] Sedangkan beberapa penelitian lain telah dilakukan untuk merangkum dan menyatukan banyak macam penggunaan GP interval.^[44][45]

GP Samar (fuzzy GP) salah satu cara selain GP Interval untuk menangani ketidakpastian nilai dalam model GP. Ketidakpastian ini umumnya berhubungan dengan besar nilai target ${\displaystyle b_{i}}$, namun juga dapat dikembangkan untuk aspek-aspek GP yang lain seperti struktur prioritas. Konsep GP samar didasarkan pada fungsi keanggotaan samar (fuzzy membership function) dalam teori himpunan samar,^[46] menyatakan besar bobot maupun rentang realisasi nilai tujuan. Sebagai contoh, jika kita ingin memaksimumkan fungsi keuntungan pada rentang nilai lima sampai dengan enam juta rupiah, kita dapat memberikan bobot ${\displaystyle 0}$ pada keuntungan dibawah lima juta, bobot ${\displaystyle 0.5}$ untuk keuntungan lima setengah juta, dan ${\displaystyle 1.0}$ untuk keuntungan di atas enam juta rupiah. Hubungan antara pembobotan dan fungsi keuntungan dapat bersifat non-linear. Hal ini memungkinkan pengambil keputusan yang tidak dapat menyatakan tujuan-tujuan dengan pasti, untuk setidaknya menyatakannya dalam bentuk bobot dan rentang realisasi nilai. Hal ini membuat GP samar sebagai pendekatan yang baik ketika tujuan-tujuan dalam model GP memiliki sifat fungsi utilitas.^[47]

Fungsi keanggotaan samar yang digunakan model GP memiliki nilai ${\displaystyle 1}$ untuk keadaan "puas seluruhnya" dan nilai ${\displaystyle 0}$ untuk keadaan "tidak-puas seluruhnya". Bentuk aljabar dari salah satu bentuk fungsi keanggotaan samar linear yang umum digunakan adalah fungsi trapesium:^[48]

$${\displaystyle \mu {\Big [}f_{i}(\mathbf {x} ){\Big ]}{\begin{cases}0&f_{i}(\mathbf {x} )\leq {_{L}b_{i}}-n_{\max }\\1-{\frac {_{L}b_{i}-f_{i}(\mathbf {x} )}{n_{\max }}}&{_{L}b_{i}}-n_{\max }\leq f_{i}(\mathbf {x} )\leq {_{L}b_{i}}\\1&{_{L}b_{i}}\leq f_{i}(\mathbf {x} )\leq {_{U}b_{i}}\\1-{\frac {f_{i}(\mathbf {x} )-\,_{U}b_{i}}{p_{\max }}}&{_{U}b_{i}}\leq f_{i}(\mathbf {x} )\leq {_{U}b_{i}}+p_{\max }\\0&{_{U}b_{i}}+p_{\max }\leq f_{i}(\mathbf {x} ).\\\end{cases}}}$$Fungsi ${\displaystyle \mu {\big [}f_{i}(\mathbf {x} ){\big ]}}$ tersebut menyatakan fungsi keanggotaan samar untuk tujuan ke-${\displaystyle i}$, dengan rentang nilai ${\displaystyle {\big [}{_{L}b_{i}},\,{_{U}b_{i}}{\big ]}}$ menyatakan daerah nilai fungsi ${\displaystyle f_{i}(\mathbf {x} )}$ "memuaskan seluruhnya". Sedangkan, nilai fungsi di luar rentang ${\displaystyle {\big [}{_{L}b_{i}}-n_{\max },\,{_{U}b_{i}}+p_{\max }{\big ]}}$ dianggap "*tidak memuaskan sama sekali*". Yaghoobi berhasil menyusun model yang menyatukan empat fungsi keanggotaan samar yang umum digunakan: segitiga, sisi-kiri, sisi-kanan, dan trapesium.^[49] Sebagai contoh, jika semua tujuan dalam permasalahan dapat dianggap sebagai fungsi trapesium tersebut, model GP yang perlu diselesaikan dapat ditulis secara aljabar sebagai berikut:

$${\displaystyle {\begin{aligned}\min &\,\sum _{i=1}^{Q}w_{i}{\Bigg (}{\frac {d_{i}^{+}}{p_{\max }}}+{\frac {d_{i}^{-}}{n_{\max }}}{\Bigg )}\\{\text{dengan kendala}}&\,f_{i}(\mathbf {x} )-d_{i}^{+}\leq {_{U}b_{i}}\\&\,f_{i}(\mathbf {x} )+d_{i}^{-}\geq {_{L}b_{i}}\\&\,\mu _{i}+{\frac {d_{i}^{+}}{p_{\max }}}+{\frac {d_{i}^{-}}{n_{\max }}}=1\\&\,\mu _{i},\,d_{i}^{+},\,d_{i}^{-}\geq 0;{\text{untuk }}i=1,\,\dots ,\,Q\\&\,x_{j}\geq 0;{\text{untuk }}j=1,\,\dots ,\,n.\end{aligned}}}$$

Pada model tersebut, ${\displaystyle \mu _{i}}$ menyatakan besar nilai keanggotaan samar (kepuasan) untuk fungsi tujuan ke-${\displaystyle i}$. Formulasi model oleh Yaghoobi memungkinkan GP berbobot samar diselesaikan menggunakan sembarang software pemrograman linear standar.

• Romero, Carlos (1991). Handbook of critical issues in goal programming. Oxford New York: Pergamon Press. ISBN 978-0-08-040661-9. 
• Jones, Dylan; Tamiz, Mehrdad (2010). "Practical Goal Programming". International Series in Operations Research & Management Science (dalam bahasa Inggris). doi:10.1007/978-1-4419-5771-9. ISSN 0884-8289. 
• Jones, Dylan; Tamiz, Mehrdad; Ries, Jana, ed. (2010). "New Developments in Multiple Objective and Goal Programming". Lecture Notes in Economics and Mathematical Systems (dalam bahasa Inggris). doi:10.1007/978-3-642-10354-4. ISSN 0075-8442. 
• Simon, Herbert A. (1955-02). "A Behavioral Model of Rational Choice". The Quarterly Journal of Economics. 69 (1): 99. doi:10.2307/1884852.  Periksa nilai tanggal di: |date= (bantuan)
• Ignizio, Jas P.; Ignizio, James P. (1982). Linear programming in single- and multiple- objective systes. Prentice-Hall fnternational series in industrial and systems engineering. Englewood Cliffs, N.J: Prentice-Hall. ISBN 978-0-13-537027-8. 
• Ignizio, James P.; Cavalier, Tom M. (1994). Linear programming. Prentice Hall international series in industrial and systems engineering. Englewood Cliffs, NJ: Prentice Hall. ISBN 978-0-13-183757-7. 
• Ignizio, James P. (1979). Goal programming and extensions: Tab. Lexington books (edisi ke-2. print). Lexington: Lexington Books. ISBN 978-0-669-00021-4. 
• Schniederjans, Marc J. (1995). "Goal Programming: Methodology and Applications". SpringerLink (dalam bahasa Inggris). doi:10.1007/978-1-4615-2229-4. 
• Ehrgott, Matthias; Gandibleux, Xavier, ed. (2002). "Multiple Criteria Optimization: State of the Art Annotated Bibliographic Surveys". International Series in Operations Research & Management Science (dalam bahasa Inggris). doi:10.1007/b101915. ISSN 0884-8289. 
• Lee, Sang M. (1972). Goal programming for decision analysis. Management and communications series (edisi ke-1. ed). Philadelphia, Pa.: Auerbach. ISBN 978-0-87769-144-0. Pemeliharaan CS1: Teks tambahan (link)

Abraham Charnes and William Wager Cooper, Management Models and Industrial Applications of Linear Programming (Wiley, 1961), [OpenLibrary:OL173690W](https://openlibrary.org/works/OL173690W/Management_models_and_industrial_applications_of_linear_programming).

• Charnes, A.; Cooper, W. W.; Ferguson, R. O. (1955-01). "Optimal Estimation of Executive Compensation by Linear Programming". Management Science (dalam bahasa Inggris). 1 (2): 138–151. doi:10.1287/mnsc.1.2.138. ISSN 0025-1909.  Periksa nilai tanggal di: |date= (bantuan)
• Masud, Abu S.; Hwang, C. L. (1981). "Interactive Sequential Goal Programming". The Journal of the Operational Research Society. 32 (5): 391–400. doi:10.2307/2581556. ISSN 0160-5682. 
• Wheeler, B. M.; Russell, J. R. M. (1977-04). "Goal Programming and Agricultural Planning". Journal of the Operational Research Society (dalam bahasa Inggris). 28 (1): 21–32. doi:10.1057/jors.1977.2. ISSN 0160-5682.  Periksa nilai tanggal di: |date= (bantuan)
• Chanas, Stefan; Kuchta, Dorota (1996-11-08). "Multiobjective programming in optimization of interval objective functions — A generalized approach". European Journal of Operational Research. 94 (3): 594–598. doi:10.1016/0377-2217(95)00055-0. ISSN 0377-2217. 
• Charnes, A.; Collomb, B. (1972-08-01). "Optimal economic stabilization policy: Linear goal-interval programming models". Socio-Economic Planning Sciences. 6 (4): 431–435. doi:10.1016/0038-0121(72)90016-X. ISSN 0038-0121. 
• Charnes, A.; Cooper, W. W. (1977-01-01). "Goal programming and multiple objective optimizations: Part 1". European Journal of Operational Research. 1 (1): 39–54. doi:10.1016/S0377-2217(77)81007-2. ISSN 0377-2217. 

• Hannan, Edward (1980-11). "Nondominance In Goal Programming". INFOR: Information Systems and Operational Research (dalam bahasa Inggris). 18 (4): 300–309. doi:10.1080/03155986.1980.11731801. ISSN 0315-5986.  Periksa nilai tanggal di: |date= (bantuan)
• Hannan, Edward L. (1981-07). "ON FUZZY GOAL PROGRAMMING*". Decision Sciences (dalam bahasa Inggris). 12 (3): 522–531. doi:10.1111/j.1540-5915.1981.tb00102.x. ISSN 0011-7315.  Periksa nilai tanggal di: |date= (bantuan)
• Inuiguchi, Masahiro; Kume, Yasufumi (1991-06-17). "Goal programming problems with interval coefficients and target intervals". European Journal of Operational Research. 52 (3): 345–360. doi:10.1016/0377-2217(91)90169-V. ISSN 0377-2217. 
• Lu, Hao-Chun; Chen, Tzu-Li (2013-02-01). "Efficient model for interval goal programming with arbitrary penalty function". Optimization Letters (dalam bahasa Inggris). 7 (2): 325–341. doi:10.1007/s11590-011-0422-z. ISSN 1862-4480. 
• Tamiz, M.; Jones, D. F. (1997). Caballero, Rafael; Ruiz, Francisco; Steuer, Ralph, ed. "An Example of Good Modelling Practice in Goal Programming: Means for Overcoming Incommensurability". Advances in Multiple Objective and Goal Programming. Lecture Notes in Economics and Mathematical Systems (dalam bahasa Inggris). Berlin, Heidelberg: Springer: 29–37. doi:10.1007/978-3-642-46854-4_3. ISBN 978-3-642-46854-4. 
• Tamiz, M.; Jones, D. F. (1996-05). "Goal programming and Pareto efficiency". Journal of Information and Optimization Sciences (dalam bahasa Inggris). 17 (2): 291–307. doi:10.1080/02522667.1996.10699283. ISSN 0252-2667.  Periksa nilai tanggal di: |date= (bantuan)
• Tamiz, M.; Jones, D. F. (1997-01). "Interactive Frameworks for Investigation of Goal Programming Models: Theory and Practice". Journal of Multi‐Criteria Decision Analysis. 6 (1): 52–60. doi:10.1002/(SICI)1099-1360(199701)6:1<52::AID-MCDA124>3.0.CO;2-3.  Periksa nilai tanggal di: |date= (bantuan)
• Tamiz, M.; Jones, D. F.; El-Darzi, E. (1995-01-01). "A review of Goal Programming and its applications". Annals of Operations Research (dalam bahasa Inggris). 58 (1): 39–53. doi:10.1007/BF02032309. ISSN 1572-9338. 
• Azaiez, M. N.; Al Sharif, S. S. (2005-03-01). "A 0-1 goal programming model for nurse scheduling". Computers & Operations Research. 32 (3): 491–507. doi:10.1016/S0305-0548(03)00249-1. ISSN 0305-0548. 
• Tamiz, Mehrdad; Jones, Dylan; Romero, Carlos (1998-12-16). "Goal programming for decision making: An overview of the current state-of-the-art". European Journal of Operational Research. 111 (3): 569–581. doi:10.1016/S0377-2217(97)00317-2. ISSN 0377-2217. 
• Min, Hokey; Storbeck, James (1991-04). "On the Origin and Persistence of Misconceptions in Goal Programming". Journal of the Operational Research Society (dalam bahasa Inggris). 42 (4): 301–312. doi:10.1057/jors.1991.68. ISSN 0160-5682.  Periksa nilai tanggal di: |date= (bantuan)
• Narasimhan, Ram (1980-04). "GOAL PROGRAMMING IN A FUZZY ENVIRONMENT". Decision Sciences (dalam bahasa Inggris). 11 (2): 325–336. doi:10.1111/j.1540-5915.1980.tb01142.x. ISSN 0011-7315.  Periksa nilai tanggal di: |date= (bantuan)
• Steuer, Ralph E. (1979-12-01). "Goal programming sensitivity analysis using interval penalty weights". Mathematical Programming (dalam bahasa Inggris). 17 (1): 16–31. doi:10.1007/BF01588223. ISSN 1436-4646. 
• Tiwari, R. N.; Dharmar, S.; Rao, J. R. (1987-10-01). "Fuzzy goal programming — An additive model". Fuzzy Sets and Systems. 24 (1): 27–34. doi:10.1016/0165-0114(87)90111-4. ISSN 0165-0114. 
• Vitoriano, B.; Romero, C. (1999). "Extended Interval Goal Programming". The Journal of the Operational Research Society. 50 (12): 1280–1283. doi:10.2307/3010638. ISSN 0160-5682. 
• Yaghoobi, M. A.; Jones, D. F.; Tamiz, M. (2008-10). "WEIGHTED ADDITIVE MODELS FOR SOLVING FUZZY GOAL PROGRAMMING PROBLEMS". Asia-Pacific Journal of Operational Research (dalam bahasa Inggris). 25 (05): 715–733. doi:10.1142/S0217595908001973. ISSN 0217-5959.  Periksa nilai tanggal di: |date= (bantuan)
• Zadeh, L. A. (1965-06-01). "Fuzzy sets". Information and Control. 8 (3): 338–353. doi:10.1016/S0019-9958(65)90241-X. ISSN 0019-9958. 

• Yuji Ijiri, Management Goals and Accounting for Control (Amsterdam : North-Holland Pub. Co., 1965), https://openlibrary.org/works/OL4814197W/Management_goals_and_accounting_for_control.
• D. F. Jones, “The Design and Development of an Intelligent Goal Programming System.” (Ph.D., University of Portsmouth, 1995), https://ethos.bl.uk/OrderDetails.do?uin=uk.bl.ethos.282556