Stochastic programming — Στοχαστικός προγραμματισμός

From Systems analysis wiki
Jump to navigation Jump to search

Στοχαστικός προγραμματισμός (αγγλ. stochastic programming) — κλάδος του μαθηματικού προγραμματισμού που αναπτύσσει μοντέλα και μεθόδους επίλυσης βελτιστοποιητικών προβλημάτων υπό συνθήκες αβεβαιότητας, όταν ορισμένες παράμετροι του μοντέλου δεν είναι ακριβώς γνωστές, αλλά αναπαρίστανται ως τυχαίες μεταβλητές με γνωστές ή εκτιμώμενες κατανομές πιθανότητας[1][2].

Σε αντίθεση με τα ντετερμινιστικά προβλήματα, όπου όλα τα δεδομένα θεωρούνται σταθερές ποσότητες, ο στοχαστικός προγραμματισμός αποσκοπεί στην εύρεση λύσης (ή πολιτικής λήψης αποφάσεων) που είναι βέλτιστη με κάποια στατιστική έννοια. Τις περισσότερες φορές αυτό σημαίνει ελαχιστοποίηση ή μεγιστοποίηση της μαθηματικής προσδοκίας της αντικειμενικής συνάρτησης[1]. Η βασική ιδέα είναι να βρεθεί μια πολιτική λήψης αποφάσεων που να είναι η καλύτερη «κατά μέσο όρο» για όλες τις πιθανές πραγματοποιήσεις των τυχαίων παραμέτρων, κάτι που είναι ιδιαίτερα κρίσιμο για προβλήματα όπου οι αποφάσεις λαμβάνονται επανειλημμένα υπό παρόμοιες συνθήκες (π.χ. στη διαχείριση αποθεμάτων ή ενεργειακών συστημάτων)[3].

Μαθηματική διατύπωση του προβλήματος

Σε γενική μορφή, το πρόβλημα στοχαστικού προγραμματισμού μπορεί να διατυπωθεί ως: minxX𝔼[f(x,ξ)] όπου:

  • x — διάνυσμα μεταβλητών ελέγχου (αποφάσεων) που πρέπει να προσδιοριστεί.
  • X — σύνολο αποδεκτών λύσεων για x, ορισμένο από ντετερμινιστικούς περιορισμούς.
  • ξ — τυχαίο διάνυσμα που αναπαριστά τις αβέβαιες παραμέτρους του προβλήματος (π.χ. ζήτηση, τιμές, καιρικές συνθήκες).
  • f(x,ξ) — αντικειμενική συνάρτηση, η τιμή της οποίας εξαρτάται τόσο από την ληφθείσα απόφαση x όσο και από την πραγματοποίηση του τυχαίου διανύσματος ξ.
  • 𝔼[] — τελεστής μαθηματικής προσδοκίας, υπολογιζόμενος ως προς την κατανομή πιθανότητας του διανύσματος ξ.

Ενα θεμελιώδης αρχή που διέπει τα πολυσταδιακά στοχαστικά μοντέλα είναι η αρχή μη-προδρομικότητας (αγγλ. non-anticipativity principle). Αυτή ορίζει ότι οι αποφάσεις που λαμβάνονται σε οποιοδήποτε στάδιο μπορούν να εξαρτώνται μόνο από πληροφορίες διαθέσιμες μέχρι εκείνη τη στιγμή και δεν μπορούν να «κοιτούν στο μέλλον»[2].

Δισταδιακό πρόβλημα με δικαίωμα αντιστάθμισης

Το πλέον διαδεδομένο μοντέλο είναι το δισταδιακό πρόβλημα με δικαίωμα αντιστάθμισης (αγγλ. two-stage stochastic program with recourse)[1]. Η διαδικασία λήψης αποφάσεων χωρίζεται σε δύο στάδια:

  1. Πρώτο στάδιο: Λαμβάνεται η απόφαση «εδώ και τώρα» (here-and-now) — καθορίζεται το διάνυσμα x. Αυτή η απόφαση πρέπει να ληφθεί πριν καταστεί γνωστή η συγκεκριμένη πραγματοποίηση του τυχαίου διανύσματος ξ.
  2. Δεύτερο στάδιο: Αφού το τυχαίο γεγονός έχει συμβεί, λαμβάνεται μια διορθωτική ή αντισταθμιστική απόφαση (recourse decision) — το διάνυσμα y(ξ), με στόχο την ελαχιστοποίηση των αρνητικών συνεπειών ή την αξιοποίηση ευνοϊκών ευκαιριών που προέκυψαν από τον συνδυασμό της απόφασης του πρώτου σταδίου x και του αποτελέσματος ξ.

Μαθηματικά, το δισταδιακό πρόβλημα στοχαστικού γραμμικού προγραμματισμού διατυπώνεται ως εξής: minxn1{cTx+𝔼ξ[Q(x,ξ)]} υπό τους περιορισμούς του πρώτου σταδίου: Ax=b,x0. Εδώ Q(x,ξ) — η συνάρτηση αντιστάθμισης (recourse function), που αναπαριστά την βέλτιστη τιμή του προβλήματος δευτέρου σταδίου: Q(x,ξ)=minyn2{q(ξ)TyT(ξ)x+Wy=h(ξ),y0} όπου ξ — τυχαίο διάνυσμα που περιλαμβάνει τις παραμέτρους q(ξ),T(ξ) και h(ξ)· ενώ c,A,b και W — ντετερμινιστικές παράμετροι[2].

Βασικές ιδιότητες και θεωρήματα

  • Κυρτότητα: Ένα από τα θεμελιώδη αποτελέσματα της θεωρίας είναι ότι για το δισταδιακό πρόβλημα στοχαστικού γραμμικού προγραμματισμού η αναμενόμενη συνάρτηση αντιστάθμισης Q(x)=𝔼ξ[Q(x,ξ)] είναι κυρτή συνάρτηση. Αυτή η ιδιότητα έχει τεράστια σημασία, καθώς εγγυάται ότι το συνολικό πρόβλημα πρώτου σταδίου είναι πρόβλημα κυρτού προγραμματισμού, για το οποίο υπάρχουν αποτελεσματικές μέθοδοι επίλυσης και το ολικό βέλτιστο συμπίπτει με το τοπικό[1].
  • Ντετερμινιστικό ισοδύναμο: Εάν το τυχαίο διάνυσμα ξ έχει πεπερασμένο αριθμό πιθανών πραγματοποιήσεων (σεναρίων) ξ1,,ξK με πιθανότητες p1,,pK, τότε το πρόβλημα στοχαστικού προγραμματισμού μπορεί να αναδιατυπωθεί ως ένα ενιαίο μεγάλο ντετερμινιστικό πρόβλημα βελτιστοποίησης. Στην περίπτωση αυτή η μαθηματική προσδοκία αντικαθίσταται από σταθμισμένο άθροισμα επί όλων των σεναρίων. Ωστόσο, το μέγεθος αυτού του προβλήματος αυξάνεται γραμμικά με τον αριθμό των σεναρίων, γεγονός που οδηγεί στην «κατάρα της διαστατικότητας» και καθιστά αυτήν την προσέγγιση υπολογιστικά ανεπίλυτη για μεγάλο αριθμό σεναρίων[2].

Σύγκριση με τη ρωμαλέα βελτιστοποίηση

Ο στοχαστικός προγραμματισμός είναι μία από τις αρκετές προσεγγίσεις βελτιστοποίησης υπό συνθήκες αβεβαιότητας. Η βασική διαφορά του από τη ρωμαλέα βελτιστοποίηση έγκειται στον τρόπο μοντελοποίησης της αβεβαιότητας και στο κριτήριο βελτιστότητας[4].

Σύγκριση προσεγγίσεων βελτιστοποίησης υπό συνθήκες αβεβαιότητας
Κριτήριο Στοχαστική βελτιστοποίηση Ρωμαλέα βελτιστοποίηση
Αναπαράσταση αβεβαιότητας Οι παράμετροι είναι τυχαίες μεταβλητές με γνωστή κατανομή πιθανότητας Οι παράμετροι ανήκουν σε ένα δεδομένο σύνολο αβεβαιότητας, κατανομή δεν απαιτείται
Κριτήριο βελτιστότητας Βελτιστοποίηση της μαθηματικής προσδοκίας της αντικειμενικής συνάρτησης Βελτιστοποίηση στο χειρότερο σενάριο (minimax)
Χαρακτήρας λύσης Πολιτική βέλτιστη «κατά μέσο όρο», μπορεί να είναι μη αποδεκτή για σπάνια σενάρια Λύση εγγυημένα αποδεκτή για όλες τις πραγματοποιήσεις· μπορεί να είναι συντηρητική

Παραδείγματα

  • Πρόβλημα εφημεριδοπώλη (αγγλ. newsvendor problem): Κλασικό πρόβλημα διαχείρισης αποθεμάτων, όπου ο πωλητής πρέπει να αποφασίσει πόσες μονάδες εμπορεύματος να αγοράσει, χωρίς να γνωρίζει τη μελλοντική ζήτηση με ακρίβεια. Η λύση εξισορροπεί τον κίνδυνο ζημίας από πλεονάσματα και τον κίνδυνο χαμένων κερδών από ελλείψεις.
  • Πρόβλημα αγρότη: Ένας αγρότης αποφασίζει πόσα στρέμματα γης να διαθέσει για διαφορετικές καλλιέργειες στη συνολική έκταση, χωρίς να γνωρίζει τις μελλοντικές καιρικές συνθήκες, οι οποίες επηρεάζουν τις αποδόσεις. Αφού οι καιρικές συνθήκες καταστούν γνωστές, ο αγρότης μπορεί να λάβει διορθωτικές ενέργειες (π.χ. να πουλήσει πλεονάσματα ή να αγοράσει την ελλείπουσα παραγωγή στην αγορά)[5].

Δείτε επίσης

  • Μαθηματικός προγραμματισμός
  • Επιχειρησιακή έρευνα
  • Ρωμαλέα βελτιστοποίηση
  • Δυναμικός προγραμματισμός
  • Θεωρία ελέγχου

Παραπομπές

[1] [2] [3] [4] [5] </references>

  1. 1.0 1.1 1.2 1.3 1.4 Shapiro, A., Dentcheva, D., & Ruszczyński, A. (2009). Lectures on Stochastic Programming: Modeling and Theory. Society for Industrial and Applied Mathematics (SIAM).
  2. 2.0 2.1 2.2 2.3 2.4 Birge, J. R., & Louveaux, F. (2011). Introduction to Stochastic Programming (2nd ed.). Springer Science+Business Media.
  3. 3.0 3.1 "Стохастическое программирование". Википедия. [1]
  4. 4.0 4.1 Gorissen, B. L., Yanıkoğlu, İ., & den Hertog, D. (2015). A practical guide to robust optimization. Omega, 53, 124-137.
  5. 5.0 5.1 Bricker, D. L. SLPwR: Farmer Problem. University of Iowa. [2]