Μετάβαση στο κύριο περιεχόμενο

Γραμμικός Προγραμματισμός: Κεφάλαιο 2.0 - λίγη γραμμική άλγεβρα & αναλυτική γεωμετρία - η γεωμετρική φύση του προβλήματος

Ως εδώ είχαμε προχωρήσει χωρίς να εισάγουμε βαρύ μαθηματικό συμβολισμό. Δυστυχώς όμως για να περάσουμε από τις 2 διαστάσεις στις n δεν μπορούμε να το κάνουμε χωρίς να βασιστούμε σε στοιχειώδη γραμμική άλγεβρα και αναλυτική γεωμετρία. Η αναγνώστρια/ης που δεν διαθέτει αυτές τις γνώσεις παραπέμπεται στη σειρά βίντεο του 3blue1brown (στα αγγλικά).

Ας σταθούμε λίγο στην παρατήρηση που είχαμε κάνει στο προηγούμενο κεφάλαιο:
Καθώς η συνάρτηση στόχος σαρώνει την εφικτή περιοχή, θα σταματήσει σε κάποια από τις κορυφές του πολυγώνου.
Αυτό το βασικό θεώρημα του γραμμικού προγραμματισμού, μας επιτρέπει να ψάχνουμε μόνο στις κορυφές του πολυτόπου, για την βέλτιστη του ΓΠ. Η κεντρική ιδέα της μέθοδου Simplex είναι η ακόλουθη:
  1. Έστω v μια κορυφή της εφικτής περιοχής 
  2. Όσο υπάρχει ένα γείτονας v' του v με καλύτερη τιμή της συνάρτησης στόχου: 
  3.     Θέσε v := v'
Το παραπάνω απαιτεί κάθε φορά να βλέπουμε το πολύτοπο από την οπτική γωνία μιας κορυφής. Θα δούμε στο επόμενο κεφάλαιο, πως αυτό είναι εύκολο να το κάνουμε για την αρχή των αξόνων. Αυτό που θα κάνουμε είναι να μετασχήματίζουμε τον χώρο, ώστε κάθε φορά η κορυφή του πολυτόπου που βρισκόμαστε να είναι η αρχή των αξόνων. Με αυτό το τελευταίο πρόβλημα θα καταπιαστούμε σε αυτό το κεφάλαιο. Δεδομένου ενός συστήματος αξόνων S και ενός σημείου \(\vec{y}\) πως μπορούμε να πάμε σε ένα σύστημα S' ώστε το \(\vec{y}\) να είναι το \(\vec{0}\).

Το πρόβλημα θα ήταν τετρημένο, αν μας δίνεται το σημείο \(\vec{y}\) (απλά \(\forall \vec{x} \in \mathbb{R}^n\) στον χώρο S, θέτουμε \( \vec{x'} := \vec{x} - \vec{y}\) στον S'), κάτι τέτοιο όμως δεν θα συμβαίνει σε ένα ΓΠ. Στο ΓΠ αυτό που μας δίνεται είναι οι περιορισμοί· και όπως είδαμε μπορούμε να τους αντιστοιχήσουμε σε ευθείες στον δισδιάστατο χώρο, σε επίπεδα στον τρισδιάστατο χώρο ή υπερεπίπεδα στον n-διάστατο χώρο. Και όπως στον δισδιάστατο χώρο μπορούμε να ορίσουμε ένα σημείο ως τομή δύο ευθειών, έτσι και στον n-διάστατο χώρο μπορούμε να ορίσουμε ένα σημείο ως τομή n υπερεπιπέδων.

Πως μπορούμε να κάνουμε την μετάβαση που θέλουμε δεδομένων των υπερεπιπέδων που σχηματίζουν το σημείο \(\vec{y}\);. Αρκεί να εκφράσουμε τις συντεταγμένες κάθε σημείου \(\vec{x} \in \mathbb{R}^n\) ως αποστάσεις, όχι πια από τους άξονες, αλλά από τα υπερεπίπεδα που ορίζουν το \(\vec{y}\).

Στον S θα έχουμε:

Σχήμα 1(α): Οι συντεταγμένες ενός σημείου, ως αποστάσεις από τους άξονες (στο S)
Σχήμα 1(α): Οι συντεταγμένες ενός σημείου, ως αποστάσεις από τους άξονες (στο S)

Ενώ στο S' θα έχουμε:

Σχήμα 1(β): Οι συντεταγμένες ενός σημείου, ως αποστάσεις από τις ευθείες που ορίζουν το Υ (στο S')
Σχήμα 1(β): Οι συντεταγμένες ενός σημείου, ως αποστάσεις από τις ευθείες που ορίζουν το Υ (στο S')

Το Σχήμα 1(β) φαίνεται περίεργο, αλλά δεν έχει τίποτα λάθος. Θυμόμαστε από στην γραμμική άλγερβρα ότι οι άξονες ενός συστήματος αναφοράς δεν χρειάζεται να είναι κάθετοι (ορθοκανονικοί), αρκεί, τα διανύσματα που τους ορίζουν να είναι γραμμικός ανεξάρτητα.

Τέλος παρατηρούμε ότι εφόσον έχουμε καταφέρει την μετάβαση από το S στο S', δεν απαιτείται κάτι το διαφορετικό για να πάμε από το S' στο S''. Επαναλαμβάνουμε το ίδιο βήμα με διαφορετικό σημείο \(\vec{y'}\), αφού το \(\vec{y}\) θα είναι πια στο S' το \(\vec{0}\).

Στο Κεφάλαιο 2.1 θα μιλήσουμε για ευθείες και υπερεπίπεδα και την σύνδεση μεταξύ την γλώσσας της άλγεβρας (που θα χρησιμοποιήσουμε για τους υπολογισμούς μας) και της γεωμετρικής εποπτείας (η οποία μας δίνει καλύτερη διαίσθηση). Ενώ στο Κεφάλαιο 2.2 θα δούμε την απόσταση σημείου από υπερεπίπεδο, όπου και θα ολοκληρώσουμε την απόδειξη του τύπου, που υλοποιεί τον μετασχηματισμό, που περιγράψαμε πιο πάνω με λόγια.

Σημειώσεις:
Αν δεν εμφανίζονται σωστά τα μαθηματικά σύμβολα, παρακαλώ πατήστε Ctrl+F5. 
Τα σχήματα έγιναν με τον GeoGebra5. Το πρόγραμμα διατίθεται για Windows, Mac και Linux.
Για μαθηματικά σύμβολα στο blogger ακολούθησα αυτές εδώ τις οδηγίες

Πλοήγηση

Σχόλια

Δημοφιλείς αναρτήσεις από αυτό το ιστολόγιο

Γραμμικός Προργαμματισμός: Κεφάλαιο 3.1 - Η Μέθοδος Simplex

Τώρα έχουμε επιτέλους τα κατάλληλα εργαλεία (συγκεκριμένα την εξίσωση (4) ), ώστε να λύσουμε ένα ΓΠ. Συγκεκριμένα ας θεωρήσουμε το ακόλουθο ΓΠ.:                 \(max\) \(z = 2x_1 + 5x_2\) υ.π.:                 (1)    \(2x_1 - x_2 \leq 4\)                 (2)    \(x_1 + 2x_2 \leq 9\)                 (3)    \(-x_1 + x_2 \leq 3\)                 (4)    \(x_1, \geq 0\)                 (5)    \(x_2, \geq 0\) Ας ονομάσουμε την παραπάνω γραφή του ΓΠ (Π1). Δεν θα κάνουμε κάτι άλλο από το να υλοποιήσου...

Γραμμικός Προγραμματισμός: Περιεχόμενα

Ο Γραμμικός Προγραμματισμός είναι εξαιρετικό μαθηματικό εργαλείο, με τεράστια εκφραστική ισχύ· πάρα πολλά προβλήματα μπορούμε να τα γράψουμε ως Γραμμικά Προγράμματα. Ακόμη προσφέρει γόνιμο έδαφος για να παρατηρηθεί η σχέση μεταξύ άλγεβρας και γεωμετρίας, η οποία αποτυπώνεται στην γλώσσα της αναλυτικής γεωμετρίας. Η γεωμετρία θα μας προσφέρει μια διαίσθηση του προβλήματος, ενώ η άλγεβρα θα επιταχύνει τους υπολογισμούς μας και θα επιτρέψει να γραφτεί κάποιο πρόγραμμα (σε υπολογιστή) για να επιλύει ΓΠ. Αυτή η μίνι σειρά μαθημάτων αποσκοπεί να εξοικειώσει τον αναγνώστη με τις βασικές έννοιες του γραμμικού προγραμματισμού. Δεν απαιτούνται προηγούμενες γνώσεις· μόνο κάποια εξοικείωση κυρίως με τα διανύσματα και τις πράξεις τους. Μία μαθήτρια λυκείου που έχει εμπεδώσει την ύλη της αναλυτικής γεωμετρίας που διδάσκεται στην Β' Λυκείου, δεν θα πρέπει να έχει πρόβλημα να διαβάσει τα κεφάλαια 0 - 3. Περιεχόμενα Κεφάλαιο 0: Εισαγωγή Κεφάλαιο 1: Γραφική Μέθοδος Επίλυσης Κεφάλαιο 2.0: Λί...

Γραμμικός Προγραμματισμός: Κεφάλαιο 2.1 - λίγη γραμμική άλγεβρα & αναλυτική γεωμετρία - ευθείες και υπερεπίπεδα

Είδαμε στα προηγούμενα κεφάλαια την σημασία της γεωμετρίας των ευθειών και των υπερεπιπέδων για τον Γραμμικό Προγραμματισμό. Σε αυτό το κεφάλαιο θα πάμε πιο βαθιά σε αυτή την έννοια και θα δείξουμε πως να την "κωδικοποιούμε" στην γλώσσα της άλγεβρας. Ευθείες, Επίπεδα και Υπερεπίπεδα Τι ακριβώς είναι μια ευθεία ή ένα υπερεπίπεδο; Πως γράφουμε μια ευθεία σε αλγεβρική μορφή; Θα πρέπει να σκεφτόμαστε την αναλυτική γεωμετρία σαν ένα υπολογιστικό σύστημα, σαν μια γλώσσα προγραμματισμού. Στην ευκλείδια γεωμετρία που κάναμε στο σχολείο η ευθεία και το επίπεδο είναι κάτι τελείως χειροπιαστό, στην άλγεβρα όμως έχουμε μόνο πράξεις, μεταβλητές και εξισώσεις, τι άλλο θα μπορούσε να είναι λοιπόν μια ευθεία από μια εξίσωση, η εξίσωση θα ισχύει μόνο για τα σημεία εκείνα τα οποία ανήκουν στην ευθεία . Στο σχολείο γράφαμε το εξής: \(y = \frac{-a}{b}x + \frac{-c}{b}\) ή \(y = \lambda x+\beta\) θα δούμε πως η πρώτη μορφή είναι πολύ πιο χρήσιμη από την δεύτερη: $$  y = \frac{-a}{b}...