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

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

Στο προηγούμενο μέρος του κεφαλαίου δείξαμε πως μπορούμε να γράφουμε τα υπερεπίπεδα και τις ευθείες σε άλγεβρα και αποδείξαμε τις εξισώσεις: 
\(\vec{x} \cdot \vec{a} = -c \hspace{5mm} (1)\)
\(H_{\vec{a}, c} = \{ \vec{x} \in \mathbb{R}^n | \vec{a} \cdot  \vec{x} = - c\} \hspace{5mm}(2)\)
 Την (1) την ονομάσαμε Εσσιανή Κανονική Μορφή. Θα δούμε ότι είναι πολύ εύκολο να βρούμε την απόσταση ενός σημείου από ένα υπερεπίπεδο, αν μας δίνεται το υπερεπίπεδο στην μορφή της (2). Θυμόμαστε ότι το \(\vec{a} \perp Η_{\vec{a}, c}\) και \(\|a\| = 1\). Επίσης, ότι το υπερεπίπεδο απέχει \(|c|\) από την αρχή των αξόνων. Έστω Χ το σημείο του οποίου θέλουμε να βρούμε την απόσταση από το \(H_{\vec{a}, c}\), τότε, ξέρουμε  ήδη από το σχολείο, ότι όταν λέμε απόσταση σημείου από υπερεπίπεδο εννοούμε κάθετη απόσταση.



Σχήμα 1(α): Κάθετη απόσταση σημείου από υπερεπίπεδο
Σχήμα 1: Κάθετη απόσταση σημείου από υπερεπίπεδο





Ας προσπαθήσουμε να "πιάσουμε" το σημείο Χ. Θεωρούμε ένα σημείο \(X_0\) πάνω στο υπερεπίπεδο, τότε παρατηρούμε ότι το μέτρο του διανύσματος \(\vec{w} := \vec(OX) - \vec{Ox_0}\) παριστάνει κάποια έννοια απόστασης από το υπερεπίπεδο· όχι όμως την κάθετη!


Σχήμα 2: Το διάνυσμα w = OX - OX_0
Σχήμα 2: Το διάνυσμα w = OX - OX_0

Ξέρουμε ότι το \(\vec{a} \perp H_{\vec{a}, c}\), οπότε αρκεί να πάρουμε την συνιστώσα του \(\vec{w}\) που είναι παράλληλη στο \(\vec{a}\). Πιο μαθηματικά θα πάρουμε την προβολή του \(\vec{w}\) στο \(\vec{a}\). Δηλαδή καταλήγουμε στην (3).

\(d(H_{\vec{a}, c}, \vec{x}) = \|proj_{\vec{a}}(\vec{x} - \vec{x_0})\| \hspace{5mm} (3)\)

και γεωμετρικά:


Σχήμα 3: Γεωμετρική αναπαράσταση του τύπου (3)
Σχήμα 3: Γεωμετρική αναπαράσταση του τύπου (3)

Για το σχήμα, αρκεί να επισημάνουμε ότι το \(\vec{a'}\) είναι το \(\vec{a}\) απλά μετατοπισμένο παράλληλα στον εαυτό του, ώστε η αρχή του να είναι το κοντινότερο σημείο του \(Η_{\vec{a}, c}\) στο Χ. Ας δουλέψουμε λίγο τον τύπο (3):

\( d(H_{\vec{a}, c}, x) \) \(= \|proj_{\vec{a}}( \vec{x} - \vec{x_0}) \| \)

                        \( = \|\frac{\vec{a}(\vec{x} - \vec{x_0})}{\|\vec{a}\|^2}\vec{a}\| \)

                        \( = \frac{ | \vec{a} ( \vec{x} - \vec{x_0} ) | }{ \| \vec{a} \|^2 } \| \vec{a} \| \)
                       
                        \(\overset{\|\vec{a}\|=1}{=}|\vec{a}(\vec{x}-\vec{x_0})|\)

                        \(=|\vec{a}\cdot\vec{x}-\vec{a}\cdot\vec{x_0}|\)
                         \(\overset{\vec{x_0}\in H_{\vec{a},c}}{=}|\vec{a}\cdot\vec{x}+c|\) (4)

Τελικά (4) : \( d(H_{\vec{a}, c}, x) \) \( =  | \vec{a}\cdot\vec{x} + c | \), (συγκρίνετε την (4) με (1)). Τώρα έχουμε ολοκληρώσει τα μαθηματικά εργαλεία που χρειαζόμαστε για να προχωρήσουμε στη μέθοδο Simplex, στο επόμενο κεφάλαιο.

Σημειώσεις:
Αν δεν εμφανίζονται σωστά τα μαθηματικά σύμβολα, παρακαλώ πατήστε 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}...