Αριθμητική Ανάλυση - Μεταπτυχιακού

Ακαδημαϊκού Έτους 2008 - 2009

Τετάρτη , 29 Οκτωβρίου 2008

2η Εβδομάδα

ΠΑΝΕΠΙΣΤΗΜΙΟ ΠΑΤΡΩΝΠΑΝΕΠΙΣΤΗΜΙΟ ΠΑΤΡΩΝΤΜΗΜΑ ΜΑΘΗΜΑΤΙΚΩΝΤΜΗΜΑ ΜΑΘΗΜΑΤΙΚΩΝΤΟΜΕΑΣ ΤΟΜΕΑΣ ΥΠΟΛΟΓΙΣΤΙΚΩΝΥΠΟΛΟΓΙΣΤΙΚΩΝ ΜΑΘΗΜΑΤΙΚΩΝ ΜΑΘΗΜΑΤΙΚΩΝΚΑΙ ΠΛΗΡΟΦΟΡΙΚΗΣΚΑΙ ΠΛΗΡΟΦΟΡΙΚΗΣ

ΑΡΙΘΜΗΤΙΚΗ ΕΠΙΛΥΣΗ ΕΞΙΣΩΣΕΩΝ (συν.)

Θέμα: Τάξη Σύγκλισης μιάς Ακολουθίας

και Επιτάχυνση Σύγκλισης μιάς

Βραδέως Συγκλίνουσας Ακολουθίας

ΑΛΓΕΒΡΙΚΕΣ ΕΞΙΣΩΣΕΙΣ

Ως γνωστόν, η τάξη σύγκλισης (order of convergence) μιάς συγκλίνουσας ακολουθίας, ορίζεται ως εξής: Η ακολουθία x0, x1, x2, ... xν,… συγκλίνει στο όριο x*

Ορίζουμε ως σφάλμα του τρέχοντα όρου xk : ek = xk – x*. Ορίζουμε την τάξη σύγκλισης P(>1), εάν υπάρχει k>0 και ισχύει:

Στην ειδική περίπτωση του p=1, τότε το k πρέπει να είναι και μικρότερο του 1 δηλαδή 0<k<1, οπότε έχουμε τάξη σύγκλισης 1 (για να υπάρχει μείωση σφάλματος).

Εξ άλλου, η σύγκλιση είναι γραμμική (linear )εάν ισχύει:

Τέλος , εάν ισχύει:

eν+1 = keν, ν=0,1,2,… με 0 < |k |<1, τότε η σύγκλιση καλείται γεωμετρική.ΑΣΚΗΣΗ: Στις γνωστές μεθόδους : διχοτόμησης, εσφαλμένης θέσης,Χορδής και Newton-Raphson, εύρατε τις τάξεις σύγκλισης αυτών.

1. Τάξη Σύγκλισης - Ορισμοί

1| |lim .

| |pe

ke

1lim και 0 | | 1.e

k ke

(α) Γεωμετρική σύγκλιση και σχήμα Aitken Εάν πάρουμε μια συγκλίνουσα ακολουθία με γεωμετρική σύγκλιση, τότε βάσει

του ορισμού τα διάφορα σφάλματα ek θα έχουν την ακόλουθη συμπεριφορά για k = ν-1, ν και ν+1:

Διά διαιρέσεως των (1) κατά μέλη λαμβάνουμε:

που επιλυόμενη ως προς x* μας δίδει:

Δηλαδή, στις ακολουθίες με γεωμετρική σύγκλιση, με τρεις όρους ( τους xν-1,

xν και xν+1 ) μπορούμε να υπολογίσουμε το όριο!

Εφαρμογή: Δείξτε ότι η ακολουθία : , ν=0, 1, 2, …,συγκλίνει γεωμετρικά, αφού βρείτε πρώτα το όριό της.

2. Επιτάχυνση Σύγκλισης

3 1 0

0.2 0.5 0

x x

x x

*1

*1

, με 0,1,2,...,(1)

και 0 | | 1, όπου το όριο.k ke e e x x

e e A x

* *1

* *1

,x x x x

x x x x

2

1*1

1 1

(2) .2

x xx x

x x x

11

2x

Ο Aitken εφήρμοσε την παραπάνω σχέση (2) σε βραδέως συγκλίνουσες

ακολουθίες, όπως η γραμμική σύγκλιση ( π.χ. στην γενική επαναληπτική ),

όπου τότε , λόγω ορισμού , ισχύουν για τα σφάλματα οι σχέσεις :

και με τον παραπάνω τρόπο από κάθε τρεις όρους της ακολουθίας x0, x1,x2,...

xν,… ( βραδέως συγκλίνουσας ), υπολόγισε έναν όρο yν μιάς νέας ακολουθίας

της : yν, ν=1,2,… με γενικό όρο:

Φυσικά, λόγω της (3),που δεν ικανοποιείται ως ισότητα, η (4) δεν δίδει το

όριο, αλλά έναν όρο που είναι κοντά σ’αυτό.

Έχει αποδειχθεί ότι η ακολουθία Aitken (4) έχει ταχύτερη σύγκλιση της

αρχικής

(β) Ακολουθία Aitken

1(3) e e με 0 | | 1,k k

2 2

1ν 1 1 2

1 1

(4) y .2

x x xx x

x x x x

..

, 1, 2,...

εφx – 1.5x + 0.1 = 0 ,

με αρχική τιμή x0=0 , αξιοποιούμε την ακολουθία της Γενικής Επαναλη-

πτικής Μεθόδου (Γ.Ε.Μ.) και υπολογίζουμε τους όρους xν , από τον τύπο :

ενώ παράλληλα κάνουμε χρήση της ακολουθίας Aitken και υπολογίζουμε τους

όρους yv, ν=1,2,…, με αποτελέσματα που δίδονται από τον πίνακα 1 : Πίνακας 1.

(γ) Παράδειγμα : Για την εύρεση μιάς ρίζας της εξίσωσης :

1 (0.1 ) /1.5 , 0,1,2,...x x

ν Xv (Γ.Ε.Μ.) Yv (Aitken)

0 0.5

1 0.066667 0.200595

2 0.111177 0.202402

3 0.141092 0.203918

4 0.161357 0.204854

... ... ...

10 0.200994 0.205905

20 0.205791

25 0.205900

Στην (4), ο όρος της Aitken που προκύπτει, θα μπορούσε να ληφθεί ως αρχικόςτης Γ.Ε.Μ. απ’ όπου έχει παραχθεί, κι έτσι οι 2 επόμενοι όροι της Γ.Ε.Μ. να πα- ραχθούν με αρχικό εκείνον της Aitken , που βρέθηκε από την (4). Αυτή η τροποποίηση δημιουργεί επιτάχυνση, χωρίς την συνέχεια της Aitken. Δηλαδή , ξεκινούμε με τρεις όρους της Γ.Ε.Μ. και υπολογίζουμε έναν ό- ρο της Aitken. Με βάση αυτόν και με Γ.Ε.Μ. υπολογίζουμε 2 επί πλέον όρους της Γ.Ε.Μ. και με τον αρχικό (Aitken) υπολογίζουμε έναν νέο όρο Aitken, κλπ. Αυτή είναι η διαδικασία τηςεπαναληπτικής Steffensen, που ενώ είναι απλού- στερη της Aitken δίδει παρομοίως επιταχύνσεις συγκλίσεων .

Παράδειγμα: Για την προηγούμενη εξίσωση και με τα ίδια δεδομένα του προβλή- ματος η ακολουθία Steffensen δίδει:

z1=0.200595, z2=0.205912, z3=0.205922,

όπου ο τρίτος όρος z3 είναι πολύ πιό κοντά στη ρίζα x*=0.2059217, από ότι ο

10ος όρος της Aitken (0.205905) και ο 25ος της Γ.Ε.Μ. (0.205900), όπωςφαίνεται από τον προηγούμενο πίνακα 1.

(δ) Παραλλαγή της Aitken - Μέθοδος Steffensen

Ως γνωστόν , η ακολουθία N-R για την εξίσωση φ(x)=0 είναι:

και x0 δοσμένο.

Έστω ότι αυτή συγκλίνει και το όριο της είναι x*, οπότε εάν γράψουμε:

( 6 ) x* = xν+ (x*-xν).τότε, το ανάπτυγμα Taylor της φ(x) (με 3 μόνο όρους) γύρω από το όριο x* δί- δει:

,

με x*<ξ<xν , ή , με ανακατάταξη των όρων και διαίρεση διά του

,Ή, ,που εάν:

δηλαδή, τετραγωνική σύγκλιση ( βάση του ορισμού της &1 ).

3. Σύγκλιση Ανωτέρας (Δευτέρας) Τάξεως - Μέθοδος Newton

1

( )(5) , 0,1,2,...

'( )

xx x

x

* * 2* * ( )

0=φ( ) ( [ ]) ( ) '( ) ''( ) ,1! 2!

x x x xx x x x x x

* * 2( ) 1 ''( )( ) 0

'( ) 2 '( )

xx x x x

x x

* * 2

1

1 ''( )( )

2 '( )x x x x

x

'( ) 0 :x

''( ) και '( ) και πάρουμε απόλυτες τιμές:M x x 2

1

1,

2e e

(α) Η επεκτεταμένη μέθοδος N-R έχει ακολουθία που δίδεται από τον τύπο:

και προκύπτει από το ανάπτυγμα Taylor της φ(x) γύρω από το xν+1=xν + h :

εάν προσδιορίσουμε το h έτσι ώστε το xν+1 να είναι ρίζα, δηλαδή φ(xν+1)=0,

πλην όμως στην απείρου βαθμού αλγεβρική εξίσωση ως προς h που προκύ-

πτει μπορούμε να την αποκόψουμε από τον τριτοβάθμιο και πέρα όρο,

οπότε όμως θα πάρουμε μιά προσέγγιση του h , από την λύση της εξίσωσης:

στην οποία επιπλέον αντικαθιστούμε το , από το γνωστό N.-R.

4. Σύγκλιση Τρίτης και Άνω Τάξης (Επεκτεταμένη Ν.-R.)

2

1 03

( ) ( ) ''( )1(7) , 0,1,2,..., δοσμένο,

'( ) 2 '( )

x x xx x v x

x x

2 3

1(8) ( ) ( ) ( ) '( ) ''( ) '''( ) ...1 2 6

h h hx x h x x x x

2

(9) ( ) '( ) ''( ) 0,1 2

h hx x x

2

2 ( )

'( )

xh

x

της κλασσικής μεθόδου ( για να αποφύγουμε την εξαγωγή ρίζας ), μιά που ο

στόχος μας είναι μια προσέγγιση του h να υπολογίσουμε.

Έτσι, έχουμε για το h τελικά την τιμή:

οπότε είναι πλέον προφανής ο υπολογιστικός τύπος της ακολουθίας (7).

(β) Για την τάξη σύγκλισης της επεκτεταμένης N.-R. θα κάνουμε χρήση

του αναπτύγματος Taylor γύρω από τη ρίζα της φ(x), οπότε έχουμε:

Οπότε, εάν διαιρέσουμε διά και θέσουμε όπως και προηγούμενα:

παίρνουμε τελικά:

2

* * 3( ) ''( ) ( )1 1 '''( )( ) 0.

'( ) 2 '( ) '( ) 6 '( )vx x x

x x x xx x x x

2( ) ''( ) ( )1

,'( ) 2 '( ) '( )

vx x xh

x x x

* * * * 2 * 31 10 ( ) ( ) ( ) ( ) '( ) ( ) ''( ) ( ) '''( ).

2 6vx x x x x x x x x x x x x

2

* ( ),

'( )

xh x x

x

'( ) 0x

Ή :

που προφανώς πληροί τον ορισμό της κυβικής σύγκλισης. Άρα, η επεκτεταμένη N-R μέθοδος έχει σύγκλιση κυβική.

(γ) Γενικότερη Μελέτη της τάξης σύγκλισης Τέλος, γενικότερα, η τάξη σύγκλισης μιάς συγκλίνουσας ακολουθίας τηςμορφής:

μπορεί να μελετηθεί, ομοίως με τη βοήθεια του αναπτύγματος Taylor, γύρω

από το όριο x* της φ(xv) ως εξής :

Εάν υποθέσουμε ότι όλες οι παράγωγοι μέχρι και την (k-1) παράγωγο της

φ(x) μηδενίζονται στο x* , ενώ φk (x*) ≠ 0, τότε από την (11), εύκολα έχουμε:

* * 3 * 31 1

1 '''( ) 1 '''( )( ) ( ) ,

6 '( ) 6 '( )x x x x ή e x x e

x x

1 0(10) ( ), 0,1,2,3,..., με δοσμένο,x x v x

2* * * * *

1 1

( ) * *

( )( ) ( ) ( ) ( ) '( ) '( )

2!(11) ( )

... ( ) ( ) .!

vv v

kkv

v

xe x x x x x x x x

xx x x

k

Προφανώς, από την (12) , η τάξη σύγκλισης της (10), βάσει ορισμού είναι k. (δ) Παράδειγμα: Τα παραπάνω εύκολα επαληθεύονται στην περίπτωση της N-R , που γράφεται:

Προφανώς, έχουμε:

οπότε :

Επιπλέον ισχύει (να αποδειχθεί) : Φ''(x*)≠0.

( )( )

1

( ) ( )(12) ( ) .

! !

k kk kv

v

xe e

k k

* *1

( )( ), με : ( ) 0 την ρίζα της συνάρτησης.

'( ) v

xx x x x x

x

2

2 2

2 2

( ) '( ) '( ) ( ) ''( )( ) , και '( ) 1

'( ) [ '( )]

[ '( )] [ '( )] ( ) ''( ) ( ) ''( )

[ '( )] [ '( )]

x x x x xx x x

x x

x x x x x x

x x

* **

* 2

( ) ''( )'( ) 0.

[ '( )]

x xx

x

Άρα η N-R συγκλίνει τετραγωνικά (έχει τάξη σύγκλισης 2), που επαλη-θεύει το συμπέρασμα της περίπτωσης (α).Άσκηση: Με παρόμοιο τρόπο δείξτε ότι η επεκτεταμένη N-R μέθοδος έχει

σύγκλιση κυβική, δηλαδή δείξτε ότι η επαναληπτική συνάρτηση του υπολογιστικού αλ- γορίθμου της επεκτεταμένης N.-R.:

πληροί τις σχέσεις:

με ξ ρίζα της : φ(x)=0. (ε) Περίπτωση πολλαπλής ρίζας της N.-R. και τροποποίηση της Η ακολουθία σφαλμάτων της μεθόδου N.-R. αποδεικνύεται (βλέπε Young - Gregory) ότι στην περίπτωση πολλαπλής ρίζας ικανοποιεί την:

Εάν τώρα η πολλαπλότητα k είναι 1 τότε προφανώς η σύγκλιση είναι τε- τραγωνική, αφού τότε ο γραμμικός όρος θα μηδενίζεται, άλλως καθίσταται

2( ) 1 ''( ) ( )

( ) ,'( ) 2 '( ) '( )

x x xx x

x x x

'( ) 0, ''( ) 0 και '''( ) 0,

,e ν=1,2,...

21

1(13) e e + ( ), όπου k η πολλαπλότητα της ρίζας.v

kO e

k

γεωμετρική , με συντελεστή ίσο με:

Το σημαντικό είναι ότι εφόσον γνωρίζουμε την πολλαπλότητα της ρίζας,τότεμπορούμε να αποκαταστήσουμε την τετραγωνική σύγκλιση εφαρμόζοντες τοακόλουθο τροποποιημένο σχήμα:

Το πρόβλημα φυσικά εμφανίζεται ως ακολουθία βραδέως συγκλίνουσα,προφανώς λόγω της πολλαπλής ρίζας, και ουδεμία γνώση έχουμε της πολ- λαπλότητος αυτής, πέρα των όρων της ακολουθίας και τηςγεωμετρικής συγκλίσεώς της. Στη περίπτωση αυτή αξιοποιούμε την παρατήρηση της παραγράφου 2α και υ-πολογίζουμε μιά προσέγγιση της ρίζας ξ από τον τύπο Aitken:

και στη συνέχεια υπολογίζουμε δύο νέους όρους της ακολουθίας N.-R.,οπότε η πολλαπλότητα k θα δίδεται από τον πλησιέστερο ακέραιο του :

1 (βλέπε γεωμετρική σύγκλιση, παρ. 1).

k

k

1

( )(14) .

'( )

xx x k

x

1 1, , ,...x x x

2

11

1 1

(15) ,2

x xx

x x x

2 3, ,x x

2

3 2

(16) .v

v v

xk

x x

Παράδειγμα: Στην εξίσωση :2

0,2

xx

εφαρμόζουμε:

(α) Το κλασσικό σχήμα Ν.-R.

(β) Το τροποποιημένο σχήμα (12), με κ=2,

οπότε εάν x0=π/2, τότε γιά τον υπολογισμό της διπλής της ρίζας έχουμε τα αποτελέσματα του πινακα 2 , που ακολουθεί :

Πίνακας 2.

s/n (α) (β) Παρατηρήσεις

X1 1.78540 2.00000 1. Ο υπολογισμός των σφαλμάτων στους 4 πρώτους όρους δίδει:

0.11009

0.05093

0.02466

0.01214

δηλαδή, το σφάλμα είναι σχεδόν το ήμισυ των προηγουμένων πράγμα που αναμένεται, αφού θα έχουμε:

ή

2. Προφανώς στην πολλαπλή ρίζα ρ θα έχουμε:

οπότε στη συνάρτηση Φ(χ) με :

η ρίζα ρ είναι απλή.

x2 1.84456 1.90100

x3 1.87083 1.89551

x4 1.88335 1.89549

x5 1.88946

X6 1.89249

x7 1.89399

X8 1.89475

x9 1.89512

x10 1.89531

x11 1.89540

x12 1.89545

x13 1.89547

x14 1.89548

x15 1.89549

1( ) ( ),x x x

1

1 1

( )( )( )

( ) ( ) ( )

x xxx

x x x x

21

1 1 ( ),ve e O e

1

1 11 .

2v v ve e e

6. Ταχύτητα σύγκλισης άλλων υπολογιστικών αλγορίθμων

(α) Μέθοδος Διχοτόμησης (Bisection method)

Με την προϋπόθεση ότι έχει εντοπισθεί μια ρίζα x* της φ(x) στο διάστημα (α,β),

τότε ισχύει η ακόλουθη σχέση φράγματος ,για τη ν-οστή επανάληψη:

οπότε η εξασφάλιση της ακρίβειας ε για τον υπολογισμό της x* απαιτεί

τουλάχιστον ν επαναλήψεις, με ν:

Σημαντικό:Κάθε κύκλος της μεθόδου βελτιώνει την ακρίβεια της ρίζας κατά 1

δυαδικό ψηφίο,οπότε για ένα δεκαδικό ψηφίο ακρίβειας απαιτούνται τρείς(3)

τουλάχιστον κύκλοι του αλγορίθμου της μεθόδου της διχοτόμησης.

*1

(17) ,2vx x

log( ) log(18) 1.

log 2

(β) Μέθοδος Εσφαλμένης Θέσεως (Ε.Θ.)

Με την προϋπόθεση του εντοπισμού της ρίζας x* της φ(x) στο διάστημα (α,β),

και επιπλέον την μονοτονία της συνάρτησης με τη συνθήκη :

οπότε η συμπεριφορά της αποδίδεται από το παρακάτω σχήμα

και τον ακόλουθο υπολογιστικό αλγόριθμο:

(Ο xv+1 είναι ο βεβαρυμένος μέσος

των xv-1 και xv με βάρη |φ(xv)| και |φ(xv-1 )|),

που στην προκειμένη περίπτωση γράφεται:

ή εάν προσθέσουμε κατά μέλη σ’ αυτή την τιμή 0 = φ(x*), έχουμε, εάν κάνουμε

χρήση του Θεωρήματος Μέσης Τιμής του Διαφορικού Λογισμού:`

1 11

1

( ) ( )(19) ,

( ) ( )v v v v

vv v

x x x xx

x x

''( ) 0,x ( , ),x

1

( ) ( )( ),v

v v vv

x ax x x

x a

που γράφεται, με τη χρήση του άνω και κάτω φράγματος

ή τελικά:

Το φράγμα σφάλματος της (20), στην προηγούμενη περίπτωση μπορεί να είναι«προβληματικό», λόγω σταθερότητος του άνω άκρου (δεξιού), των διαδοχικώνδιαστημάτων εγκλεισμού της ρίζας, πράγμα που συμβαίνει στο παρακάτωπαράδειγμα.(γ) Παράδειγμα :Στην συνάρτηση φ(x) η ακολουθία (19) της Ε.Θ. για δίδεται από την:

από την οποία μπορούμε να υπολογίσουμε το πλήθος ν των όρων που απαι- τούνται ώστε να έχουμε όρο πέραν του ½, που είναι:

και αν λάβουμε δ=10-6 τότε απαιτούνται 693.000 επαναλήψεις,

1

v

v v

x ax x

1'( )

vx a

* *

2( ) ( ) ' ,v vx x x x

'( ) :x M

* 1 21 1

2

'( ) '( ),

'( )v v vx x x x

*1 1(20) .v v v

Mx x x x

12vx

1

(21) 1 ,1

vx

1

2 vx

0.693log 2 / log(1 )v

πράγμα που δείχνει την προβληματικότητα του αλγορίθμου, στην προκειμένη

περίπτωση που παραμένει ένα άκρο του διαστήματος σταθερό, οπότε

ενδείκνυται η χρήση τροποποιημένων σχημάτων όπως είναι το Illinois, ή το

Pegasus, κλπ.

(δ) Τροποποιημένα σχήματα εσφαλμένης θέσεως

Η γενική φιλοσοφία των τροποποιημένων σχημάτων είναι η υπέρβαση της

αδυναμίας με κατάλληλη μείωση (π.χ. με υποδιπλασιασμό) της τιμής της

συναρτήσεως στο άκρο που παραμένει σταθερό σε επόμενο κύκλο. Π.χ. στο

σχήμα το άκρο 1 παραμένει σταθερό στον δεύτερο κύκλο, άρα τότε αντί για την

φ(1) παίρνουμε την τιμή kφ(α), με k<1,

οπότε τότε αντί για τον όρο x3 θα είχαμε

τον όρο x'3, που είναι πολύ πλησιέστερα

στη ρίζα Η συγκεκριμένη τιμή του k

είναι θέμα μεθόδου. Εάν τότε

έχουμε το σχήμα Illinois. Εάν

τότε έχουμε το σχήμα Pegasus

.x12 ,k

1 1( ) | ( ) ( )v v vk x x x

Τέλος, μια άλλη επιλογή είναι να πάρουμε τις διηρημένες διαφορές:

Αποδεικνύεται ότι η μέση ταχύτητα σύγκλισης της Illinois είναι

ενώ της Pegasus είναι και της τελευταίας είναι

(ε) Η μέθοδος της χορδής (Secant method)

Η μέθοδος της χορδής δεν απαιτεί εγκλεισμό της ρίζας μέσα σε ένα διάστημα

και δεν κάνει χρήση της παραγώγου, έχει «ικανοποιητική ταχύτητα» όταν

συγκλίνει: πλην όμως μπορεί να αστοχήσει! Πάντως η απλότητα

του υπολογιστικού της αλγορίθμου την κάνει ελκυστική:

και οδηγεί στην Newton Raphson όταν ξεκινήσουμε όχι από 2 σημεία της

συνάρτησης φ(x) (που καθορίζουν τη χορδή), αλλά από ένα σημείο και την

εφαπτομένη της συνάρτησης εκεί, που εξασφαλίζει την τετραγωνική

σύγκλιση της μεθόδου.

31.442 3Iu

1

1 2

, αν >0, ( ) ( )

με , όπου [ , ] .1,, αν 0

2

v v

v v

x x x yk x y

x x x y

41.642 7.275Pu . . 1.68.u

1.618,xu

11

1 1

1

( )( )(22) ,

( ) ( ) ( ) ( )v v vv

v v vv v v v

v v

x x xxx x x

x x x xx x

Ως εφαρμογή, ας πάρουμε την εξίσωση:Παράδειγμα:

0,2

xx

που εύκολα διαπιστώνεται ότι η ρίζα της είναι στο διάστημα (π/2, π). Στη συνέχεια θα εφαρμόσουμε τις μεθόδους της εσφαλμένης θέσεως, καθώς και τις τροποποιημένες μορφές της Illinois και Pegasus, με αρχικό διάστημα το (π/2, π), οπότε λαμβάνουμε τα παρακάτω αριθμητικά αποτελέσματα για x0=π/2 και x1=π:

Xκ/κ Εσφ. θέσης Illinois Pegasus

2 1.75960 1.75960 1.75960

3 1.84420 1.91904 1.99169

4 1.87701 1.89349 1.88772

5 1.88895 1.89547 1.89509

6 1.89320 1.89552 1.90226

7 1.89469 1.89549 1.89549

8 1.89521

9 1.89540

10 1.89546

11 1.89548

12 1.89549

7. Αλγεβρικές εξισώσεις (α) Μέθοδος Νewton - Raphson Στην περίπτωση των αλγεβρικών εξισώσεων η μέθοδος N.-R. μπορεί νααξιοποιηθεί για τον υπολογισμό κάθε ρίζας ( πραγματικής ή μιγαδικής ) τηςεξίσωσης, με βέλτιστο τρόπο υπολογισμού των τιμών της συνάρτησης καιτων παραγώγων αυτής πάσης τάξεως, με εφαρμογή σχημάτων Horner, που ωςγνωστό έχουν απλό αλγόριθμο. I. Περίπτωση πραγματικής ρίζας Ας υποθέσουμε την γενική μορφή της αλγεβρικής εξισώσεως κ - βαθμού:

και την εφαρμογή του κλασσικού σχήματος N.-R.:

Η γραφή του Pk στην μορφή Horner:

καθιστά τον υπολογισμό του Pk(x) βέλτιστο ,με k το πολύ πολλαπλασιασμούς και k το πολύ προσθέσεις.

1 2 21 2 2 1 0(23) ( ) ... 0,k k k

k k k kP x a x a x a x a x a x a

1 0'

( )(24) , 0,1,2,... και δοσμένο.

( )k v

v vk v

P xx x v x

P x

0 1 2 3 1(25) ( ) ... ... ,k k kP x a x a x a x a x a a x

Έτσι, για x = x0 η (25) μπορεί να υλοποιηθεί με τον αλγόριθμο:

Προφανώς στον (26) οι συντελεστές δίδουν τις εντόςτων παρενθέσεων εκφράσεις της (25) ενώ ο και ο Τέλος, εύκολα μπορούμε να αποδείξουμε ότι:

με:

αρκεί να εκτελεσθούν οι πράξεις στο δεξιό μέλος της (27) και να εξισωθούν οισυντελεστές των ομοίων δυνάμεων. Εξάλλου, με παραγώγιση της (27), συνάγεται εύκολα η σχέση:

που σημαίνει ότι στην (24) ο αριθμητής και ο παρονομαστής του κλάσματος

μπορούν να υπολογιστούν με δύο συνεχείς εφαρμογές του σχήματος Hornerγια τον υπολογισμό των τιμών πράγμα που απλοποιεί τηνόλη διαδικασία.

1 2 2 1 0

0 0 0 1 0 3 0 2 0 1(1)

1 2 2 1 0 0

( ) ...

(26) ...

( ) ... ( )k

k k k k

k k

k k k k

P x a a a a a a

x x x x x x x

P x P x

1 2 2 1, , ,..., και k k k

k ka 0 0( ).kP x

(1)0 0(27) ( ) ( ) ,

kkP x x x P x

(1) 1 21 2 1(28) ( ) ,

k

k kk kP x x x x

' (1)0 0(29) ( ) ( ),k kP x P x

'( ) και ( ),k v kP x P x

II. Περίπτωση μιγαδικού ζεύγους ριζών

Το σχήμα N.-R. (24) μπορεί και πάλι να εφαρμοσθεί, αρκεί η αρχική τιμή x0 να

είναι μιγαδικός αριθμός, π.χ.

οπότε παρουσιάζεται η ανάγκη υπολογισμού τιμών πολυωνύμων για μιγαδική

ανεξάρτητη μεταβλητή. Μια τροποποίηση του σχήματος Horner μπορεί να

αξιοποιηθεί για την απλοποίηση των πράξεων και την αποφυγή μιγαδικής

αριθμητικής. Η τροποποίηση συνίσταται στην χρήση του τριωνύμου που έχει

ως ρίζες τις δηλαδή υπολογίζουμε το άθροισμα και το γινόμενό τους :

οπότε το τριώνυμο θα είναι το:

ενώ το μπορεί να γραφεί ως εξής (π.χ. ταυτότητα της διαίρεσης ) :

με:

0(30) ,x a i

,a i 2 22 , ,k a i a i a a i a i a

2 ,x kx ( )kP x

2 (2)(31) ( ) ( ) ,kkP x x kx P x Ax B

(2) 2 31 3 2( ) ,k k

k k kP x B x B x B x B

ενώ οι συντελεστές : θα δίδονται από τις :

Δηλαδή θα έχουμε για τους συντελεστές τις σχέσεις προσδιορισμού , που λαμ-

βάνονται από την (31), με εξίσωση των συντελεστών των ομοιόβαθμων όρων :

Οπότε στο τέλος, για τον υπολογισμό της τιμής του πολυωνύμου θα έχουμε να

πολλαπλασιάσουμε την τιμή του A επί τον μιγαδικό α+βi και να προσθέσουμε

την τιμή του B, όπως φαίνεται σαφώς στην τελευταία σχέση.

:

1 2 3 2, , ,..., , , ,k k kB B B B B A B

1 2 2 1 0

1 23

34 2

(2)1 2 2

( ) ...

(32) ...

( ) ...k

k k k k

kk

k

k k k

P x a a a a a a

kB kBkBk kB

B BB B

P x B B B B B

1 21 1

1 1 2 3 0 22 2 1

, 3, 4,..., 2

,

( ) ( ) .

k ki i i i

k k k

k k k k

k

B aB a kB B i k k

B a kBA B kB B a B

B a kB B

P i A i B

Παράδειγμα :

Στο πολυώνυμο :

ζητείται να υπολογιστούν με βέλτιστο τρόπο οι τιμές :

Aπάντηση : Προφανώς, θα κάνουμε χρήση του σχήματος Horner για μιγαδικούς, οπότε κατά τα προαναφερθέντα βρίσκουμε το άθροισμα κ=2 και το γινόμενο -λ = 5, καθώς και την σχετική διάταξη του σχήματος Horner για μιγαδική τιμή :

Άρα : Α=-50 και Β= 34 – 2(-50)=134, οπότε Ρ(1+2i) =-50(1+2i)+134=84-100i, και για την συζυγή τιμή θα έχουμε Ρ(1-2i) =84+100i.

5 2( ) 2 3 6 1,

(1 2 ), (1 2 ) .i i

2

( ) : 2 0 0 3 6 1

2 : 4 8 4 54 100

5 : 10 20 10 135

_________________________________________

( ) : 2 4 2 27 50 34