برنامه‌ریزی خطی

برنامه‌ریزی خطی، یا همان بهینه‌سازی خطی، روشی در ریاضیات است که به پیدا کردن مقدار کمینه یا بیشینه از یک تابع خطی روی یک چندضلعی (یا چندوجهی) محدب می‌پردازد. این چندضلعی محدب در حقیقت نمایش نموداری تعدادی محدودیت از نوع نامعادله روی متغیرهای تابع است. به بیان ساده‌تر، به وسیله برنامه‌سازی خطی می‌توان بهترین نتیجه (مثلاً بیشترین سود یا کمترین هزینه) را در شرایط خاص و با محدودیت‌های خاص به دست آورد. محل اصلی استفاده برنامه‌ریزی خطی در مدیریت و اقتصاد است، اما در مهندسی نیز کاربردهای فراوانی دارد. در واقع برنامه‌ریزی خطی بخشی از برنامه‌ریزی ریاضی و تحقیق در عملیات موسوم به علم مدیریت است که اولین بار توسط نیروی هوایی ارتش آمریکا بکار گرفته شد. می‌توان گفت حدود یک‌چهارم کل محاسبات علمی که بر روی رایانه انجام گرفته است، به برنامه‌ریزی خطی و مشتقات آن مربوط می‌شود.

تاریخچه

مسئلهٔ حل مجموعه‌ای از نامعادلات خطی از زمان فوریه مطرح بوده است. برنامه‌ریزی خطی به عنوان یک مدل ریاضی در زمان جنگ جهانی دوم شکل گرفت تا خرج‌ها و بازگشت‌های مالی را طوری سامان بخشد که به کاهش هزینه‌های ارتش و افزایش خسارات دشمن بینجامد. این طرح تا سال ۱۹۴۷ سری باقی ماند. پس از جنگ، بسیاری از صنایع به استفاده از آن پرداختند. پایه‌گذاران این حوزه جورج دانتزیگ منتشرکنندهٔ روش سیمپلکس در سال ۱۹۴۷، جان فون نویمان مطرح‌کننده نظریه دوگانگی در همان سال، و لئونید کانتروویچ ریاضیدان روس که از تکنیک‌های مشابهی پیش از دانتزینگ استفاده کرد و نوبل سال ۱۹۵۷ را برد هستند. نخستین بار در سال ۱۹۷۹ لئونید خاچیان نشان داد که مسئله برنامه‌ریزی خطی در مرتبه زمانی چندجمله‌ای قابل حل است. اما پیشرفت اساسی‌تر زمانی حاصل شد که نراندرا کارمارکار یک روش نقطه داخلی جدید برای حل این مسائل معرفی کرد. مثال دانتزینگ برای منتصب کردن هفتاد نفر به هفتاد شغل متمایز کارآمدی برنامه‌ریزی خطی را به نمایش می‌گذارد. توان محاسباتی لازم برای آزمودن همهٔ جایگشتهای ممکن این مسئله بسیار بالاست. این تعداد از تعداد ذرات موجود در عالم بیشتر است. با این حال، پیدا کردن پاسخ بهینه با تبدیل مسئله به یک مسئله برنامه‌ریزی خطی و حل آن با روش سیمپلکس تنها لحظه‌ای طول می‌کشد.

الگوریتم‌ها

مجموعه‌ای از محدودیت‌ها (خطوط مرزی) به صورت نامعادلات خطی روی دو متغیر منجر به ایجاد منطقه‌ای از مقادیر ممکن برای آن دو متغیر روی صفحه می‌شود. این منطقه برای مسائل حل‌شدنی به شکل یک چندضلعی محدب است.

الگوریتم سیمپلکس که توسط جورج دانتزینگ شکل گرفت، مسائل برنامه‌ریزی خطی را به این ترتیب حل می‌کند که یک جواب قابل قبول در یکی از رئوس چندضلعی فراهم می‌کند و سپس در راستای اضلاع چندضلعی به طرف رئوسی با مقدار بالاتری از تابع هدف حرکت می‌کند تا این که به نقطه بهینه برسد. اگرچه در عمل این الگوریتم بسیار کارآمد است و می‌تواند با در نظر گرفتن برخی پیش‌گیری‌های مربوط به جلوگیری از ایجاد دور، با اطمینان جواب بهینه مطلق را بیابد، اما در حالاتی که به اصطلاح بدترین حالت نامیده می‌شوند عملکرد بدی دارد. تا حدی که می‌توان مسائل برنامه‌ریزی خطی طراحی کرد که روش سیمپلکس برای حلشان در برخی مراحل زمانی از مرتبه زمانی نمایی نیاز داشته باشد. حتی در دورانی دانشمندان نمی‌دانستند که این مسائل راه حل چندجمله‌ای هم دارند.

سرانجام این مسئله را لئونید خاچیان در سال ۱۹۷۹ با ارائه روش بیضوی حل کرد. این روش در بدترین حالت هم دارای زمان اجرای چندجمله‌ای بود. این روش تأثیر چندانی در جنبهٔ عملی مسئله نداشت چرا که همچنان روش سیمپلکس در همه موارد به جز تعداد محدودی از مسائل بهتر عمل می‌کرد. اما اهمیت نظری روش خاچیان غیرقابل‌انکار بود. این روش الهام‌بخش به وجود آمدن نسل جدیدی از راه‌حل‌ها شد که به آن‌ها روش نقطه داخلی گفته می‌شود. در این روش‌ها نقاط داخلی محدوده قابل بررسی متغیرها پیموده می‌شود و به سمت نقطه بهینه حرکت انجام می‌گیرد.

برنامه‌ریزی خطی استوار

در بسیاری از مسایل بهینه‌سازی، تابع هدف یا قیود، تصادفی و در حال تغییر هستند؛ بنابراین پاسخ بهینهٔ مسئله می‌بایست به گونه‌ای بدست آورده شود که برای تمامی تغییرات توابع بهینه باشد.

یک مسئله برنامه‌ریزی خطی (LP) به صورت زیر فرض می‌شود: الگو:چپ‌چین $minimize{C}^{T}x$الگو:سخ $subjectto{a}_i^{T}x\le b$ الگو:پایان چپ‌چین اگر در مسئلهٔ فوق پارامترهای $a_i$ ،$b_i$ و C ثابت و یقینی نباشند، حل مسئله، به حل یک مسئلهٔ مقاوم در برابر تغییرات پارامترها تبدیل می‌شود.

برای حل مسائل بهینه‌سازی به صورت مقاوم دو رهیافت وجود دارد:

1. بدترین حالت
2. مدل تصادفی

بدترین حالت

در این روش، بازهٔ تغییرات پارامترها مشخص است و مسئله برای بدترین حالتی که ممکن است رخ دهد، حل می‌شود.

مثال

در حالت خاص فرض می‌شود تنها پارامتر $a_i$ متغیر تصادفی باشد. در صورتی که$a_i\in {\epsilon }_i$ فرض شود، ($a_i$ عضو یک بیضی باشد).

الگو:چپ‌چین $minimize{C}^{T}x$الگو:سخ $subjectto{a}_i^{T}x\le b_{i}for{a}_i\in {\epsilon }_i$ الگو:پایان چپ‌چین قید نامساوی را می‌توان به صورت زیر نوشت: الگو:چپ‌چین $Sup(a_i^{T}x|a_i\in {\epsilon }_i)\le b_i$, $a_i=\{{\overline{a}}_i^T+P_{i}u|{\Vert u\Vert }_2\le 1\}$ الگو:پایان چپ‌چینالگو:سخ با جایگذاری $a_i$ در $Sup$ خواهیم داشت: الگو:چپ‌چین $Sup(a_i^{T}x|a_i\in {\epsilon }_i)$ $={\overline{a}}_i^{T}x+Sup\{u^T{P}_i^{T}x|{\Vert u\Vert }_2\le 1\}$

$={\overline{a}}_i^{T}x+{\Vert P_i^{T}x\Vert }_2$ الگو:پایان چپ‌چینالگو:سخ در نتیجه مسئله برنامه‌ریزی خطی مقاوم به یک مسئلهٔ SOCP تبدیل می‌شود: الگو:چپ‌چین $minimize{C}^{T}x$الگو:سخالگو:سخ $subjectto{\overline{a}}_i^{T}x+{\Vert P_i^{T}x\Vert }_2$ الگو:پایان چپ‌چین

مدل تصادفی

در این روش، متغیر در بازه مشخص نبوده و دارای توزیع آماری است بنابراین بهینه‌سازی برای بدترین حالت امکان‌پذیر نبوده و مسئله با توجه به توزیع آماری متغیر حل می‌شود.

مثال

الگو:چپ‌چین $minimize{C}^{T}x$الگو:سخ $subjecttoprob(a_i^{T}x\le b_i)\ge \eta$ الگو:پایان چپ‌چین $a_i$ یک بردار گوسی با میانگین ${\overline{a}}_i^T$ و تابع کوواریانس ${\mathrm{\Sigma }}_i$ فرض می‌شود. الگو:چپ‌چین $a_i\sim N({\overline{a}}_i^T,{\mathrm{\Sigma }}_i)$الگو:سخ $a_{i}x\sim N({\overline{a}}_i^{T}x,x^T{\mathrm{\Sigma }}_{i}x)$الگو:سخ الگو:پایان چپ‌چین قید مسئله به صورت زیر خواهد شد: الگو:چپ‌چین $prob(a_i^{T}x\le b_i)=1-Q(\left(\frac{b_i-a_i^{T}x}{{\Vert {\mathrm{\Sigma }}_i^{1/2}x\Vert }_2}\right)$الگو:سخ الگو:پایان چپ‌چین $\varphi (x)=1-Q(x)$ تابع CDF متغیر گوسی است. الگو:چپ‌چین

$prob(a_i^{T}x\le b_i)=\varphi (\left(\frac{b_i-a_i^{T}x}{{\Vert {\mathrm{\Sigma }}_i^{1/2}x\Vert }_2}\right)$$\ge \eta$الگو:سخ $b_i-a_i^{T}x\ge {\varphi }^{-1}(\eta ){\Vert {\mathrm{\Sigma }}_i^{1/2}x\Vert }_2$

$a_i^{T}x+{\varphi }^{-1}(\eta ){\Vert {\mathrm{\Sigma }}_i^{1/2}x\Vert }_2\le b_i$ الگو:پایان چپ‌چین با جایگذاری قید بدست آمده در مسئله اصلی، مسئله به صورت زیر خواهد شد: الگو:چپ‌چین $minimize{C}^{T}x$الگو:سخ $subjectto{a}_i^{T}x+{\varphi }^{-1}(\eta ){\Vert {\mathrm{\Sigma }}_i^{1/2}x\Vert }_2\le b_i$ الگو:پایان چپ‌چین مسئلهٔ برنامه‌ریزی خطی مقاوم حاصل به شکل یک مسئلهٔ SOCP است.

منابع

• فردریک س. هیلیر- جرالد ج. لیبرمن، ترجمه محمد مدرس و اردوان آصف‌وزیری، تحقیق در عملیات، چاپ دهم تهران: نشر جوان، ۱۳۸۲

• سلیمانی دامنه- مجید، بهینه‌سازی خطی ، انتشارات دانشگاه تهران، ۱۴۰۰

پیوند به بیرون

• Guidance on Formulating LP problems
• 0-1 Integer Programming Benchmarks with Hidden Optimum Solutions الگو:Webarchive
• Mathematical Programming Glossary
• A Tutorial on Integer Programming
• The linear programming FAQ
• Linear Programming Survey OR/MS Today
• Linear Programming: Guide to Formulation, Simplex Algorithm, Goal Programming and Excel Solver examples
• George Dantzig

جستارهای وابسته

• برنامه‌ریزی پویا
• بهینه‌سازی
• بهینه‌سازی استوار
• سیمپلکس