next up previous contents index
Next: ΚΡΙΤΙΚΗ Up: ΡΙΖΕΣ ΜΗ-ΓΡΑΜΜΙΚΩΝ ΕΞΙΣΩΣΕΩΝ Previous: ΠΑΡΑΔΕΙΓΜΑ   Contents   Index


ΜΕΘΟΔΟΣ ΔΙΧΟΤΟΜΗΣΗΣ (BOLZANO)

Στο προηγούμενο παράδειγμα, ας προσπαθήσουμε να εφαρμόσουμε μια απλοϊκή ίσως διαδικασία. Υποθέτουμε ότι το κεφάλαιο που κατατίθεται κάθε μήνα είναι 1000 και το ποσό που έχουμε στόχο να συλλέξουμε σε διάστημα 50 μηνών είναι 65000.

Αν δοκιμάσουμε λοιπόν για δυο τυχαίες (αλλά λογικές) τιμές του επιτοκίου, έστω $x_{1} = 0.10$ και $x_{2} = 0.15$, θα βρούμε:

$\displaystyle f\left( {0.10} \right)$ $\textstyle =$ $\displaystyle 3286.4$  
$\displaystyle f\left( {0.15} \right)$ $\textstyle =$ $\displaystyle - 3881.8$  

που σημαίνει ότι με επιτόκιο 0.10 (δηλαδή 10% ετησίως), υπολειπόμεθα του στόχου μας κατά 3286.4, ενώ με επιτόκιο 0.15 τον ξεπερνούμε κατά 3881.8. Επομένως, η φυσική μας επιλογή θα είναι να δοκιμάσουμε για επιτόκιο 0.125, δηλαδή μια τιμή στο μέσο του διαστήματος $(0.1, 0.15)$. Για επιτόκιο $x_{3} = 0.125$ βρίσκουμε ότι:

\begin{displaymath}
f(0.125) = - 174.5
\end{displaymath}

άρα απέχουμε από το στόχο μας μόνο κατά 174.5 ευρώ, ενώ φαίνεται πλέον πως η σωστή τιμή του επιτοκίου βρίσκεται στο διάστημα $(0.10, 0.125)$. Επομένως, επιλέγοντας μια τιμή στο μέσο του διαστήματος, έστω $x_{4} = 0.1125$, βρίσκουμε ότι

\begin{displaymath}
f\left( {0.1125} \right) = 1585.6 \, .
\end{displaymath}

Επομένως η σωστή τιμή βρίσκεται στο διάστημα $(0.1125, 0.125)$, άρα μπορούμε να επιλέξουμε ένα νέο επιτόκιο, που θα είναι $x_{5} =
(0.1125 + 0.125)/2$, και να επαναλάβουμε την παραπάνω διαδικασία. Μετά από μερικές ακόμη προσπάθειες, θα καταλήξουμε στην τιμή του κεφαλαίου που επιθυμούμε να συγκεντρώσουμε μετά από καταθέσεις 50 μηνών.

Η `` απλοϊκή '' μέθοδος που χρησιμοποιήσαμε στο προηγούμενο παράδειγμα αποτελεί μια από τις κλασικές μεθόδους εύρεσης ριζών εξισώσεων και ονομάζεται μέθοδος διχοτόμησης ή μέθοδος του Bolzano. Η παραπάνω διαδικασία μπορεί να κωδικοποιηθεί ως ακολούθως :

`Εστω ότι μια ρίζα βρίσκεται στο διάστημα $[a_0,b_0]$, τότε $f(a_0)\cdot f(b_0)< 0$. Αν $\mu _{0}=(a_0+b_0)/2$, τότε:

Αν ισχύει η (III), τότε έχει υπολογισθεί η ρίζα και σταματά η διαδικασία, αλλιώς ορίζω νέο διάστημα
\begin{displaymath}
\left[a_1,b_1 \right] = \left\{ {\begin{array}{l}
\left[\...
...0\right] \quad \mbox{αν} \quad (I) \\
\end{array}} \right.
\end{displaymath} (7)

Figure: Γραφική απεικόνιση της διαδικασίας που αναπτύχθηκε για τη μέθοδο διχοτόμησης
\includegraphics[width=8.5cm,height=7.5cm]{Fig_1.1.ps}


Table: Σχέδιο προγράμματος για τη εύρεση ρίζας μη-γραμμικής εξίσωσης με τη μέθοδο της διχοτόμησης.
REPEAT
SET $x_3=(x_1+x_2)/2$
IF $f(x_3) \cdot f(x_1) < 0$
SET$x_{2}=x_{3}$
ELSE
SET$x_{1}=x_{3}$
ENDIF
UNTIL $(\vert x_1-x_2\vert < E )$ OR $f(x_3)=0$




Subsections
next up previous contents index
Next: ΚΡΙΤΙΚΗ Up: ΡΙΖΕΣ ΜΗ-ΓΡΑΜΜΙΚΩΝ ΕΞΙΣΩΣΕΩΝ Previous: ΠΑΡΑΔΕΙΓΜΑ   Contents   Index
Kostas Kokkotas 2005-06-13