Σύνολο Julia

μικρογραφία|352x352εσ|Σύνολο Julia για Zn+1 = Zn2 + C C = -0,8 + i0,156. Η εικόνα είναι κεντραρισμένη στην αρχή (0,0) χωρίς ζουμ.

Στο πλαίσιο της σύνθετης δυναμικής, ενός κλάδου των μαθηματικών, το **σύνολο Julia** (Τζούλια) και το σύνολο Φατού είναι δύο συμπληρωματικά σύνολα (τα "κορδόνια" του Julia και η "σκόνη" του Φατού) που ορίζονται από μια συνάρτηση. Ανεπίσημα, το σύνολο Φατού της συνάρτησης αποτελείται από τιμές με την ιδιότητα ότι όλες οι κοντινές τιμές συμπεριφέρονται παρόμοια όταν η συνάρτηση επαναλαμβάνεται επανειλημμένα, ενώ το σύνολο Τζούλια αποτελείται από τιμές τέτοιες ώστε μια αυθαίρετα μικρή διαταραχή μπορεί να προκαλέσει ριζικές αλλαγές στην ακολουθία των επαναλαμβανόμενων τιμών της συνάρτησης. Έτσι, η συμπεριφορά της συνάρτησης στο σύνολο Φατού είναι "κανονική", ενώ στο σύνολο Τζούλια η συμπεριφορά του είναι "χαοτική".

Το σύνολο Julia μιας συνάρτησης  συμβολίζεται συνήθως ως \operatorname{J}(f), και το σύνολο Fatou συμβολίζεται ως  \operatorname{F}(f). [α] Τα σύνολα αυτά πήραν το όνομά τους από τους Γάλλους μαθηματικούς Γκαστόν Τζούλια  και Πιέρ Φατού, το έργο των οποίων εγκαινίασε τη μελέτη της μιγαδικής δυναμικής στις αρχές του 20ού αιώνα.

## Επίσημος ορισμός
Έστω f(z) μια μη σταθερή ολομορφική συνάρτηση της σφαίρας Ρίμαν στον ίδιο της τον εαυτό. Τέτοιες συναρτήσεις f(z) είναι ακριβώς οι μη σταθερές μιγαδικές ορθολογικές συναρτήσεις, δηλαδή f(z) = p(z)/q(z) όπου p(z) and q(z) είναι μιγαδικά πολυώνυμα. Ας υποθέσουμε ότι τα *p* και *q* δεν έχουν κοινές ρίζες και ότι τουλάχιστον ένα από αυτά έχει βαθμό μεγαλύτερο από 1. Τότε υπάρχει πεπερασμένος αριθμός ανοικτών συνόλων  F_1, ..., F_r τα οποία είναι αναλλοίωτα από το f(z) και είναι τέτοια ώστε :

#Η ένωση των συνόλων F_i είναι πυκνή στο επίπεδο και
#η f(z) συμπεριφέρεται κανονικά και εξίσου σε κάθε ένα από τα σύνολα F_i.

Η τελευταία πρόταση σημαίνει ότι τα άκρα των ακολουθιών των επαναλήψεων που δημιουργούνται από τα σημεία του F_i είναι είτε ακριβώς το ίδιο σύνολο, το οποίο είναι τότε ένας πεπερασμένος κύκλος, είτε κύκλοι συνόλων κυκλικής μορφής ή δακτυλιοειδών συνόλων που βρίσκονται συγκεντρωτικά. Στην πρώτη περίπτωση, ο κύκλος είναι "ελκυστικός", στη δεύτερη περίπτωση είναι "ουδέτερος".

Αυτά τα σύνολα F_i είναι τα πεδία Φατού της f(z), και η ένωσή τους είναι το σύνολο Φατού \operatorname{F}(f) της f(z). Καθένα από τα πεδία Φατού περιέχει τουλάχιστον ένα κρίσιμο σημείο της f(z), δηλαδή ένα (πεπερασμένο) σημείο *z* που ικανοποιεί το f'(z) = 0, ή f(z) = \infty  αν ο βαθμός του αριθμητή p ( z ) p(z) είναι τουλάχιστον διπλάσιος από τον βαθμό του παρονομαστή q ( z ) q(z), ή αν f(z) = 1/g(z) + c για κάποιο *c* και μια ορθολογική συνάρτηση g(z) που ικανοποιεί αυτή τη συνθήκη.

Το συμπλήρωμα του \operatorname{F}(f) είναι το σύνολο Julia \operatorname{J}(f) του f(z). Αν όλα τα κρίσιμα σημεία είναι προπεριοδικά, δηλαδή δεν είναι περιοδικά αλλά τελικά προσγειώνονται σε έναν περιοδικό κύκλο, τότε το \operatorname{J}(f) είναι όλη η σφαίρα. Διαφορετικά, το \operatorname{J}(f) είναι ένα πουθενά πυκνό σύνολο (είναι χωρίς εσωτερικά σημεία) και ένα μη μετρήσιμο σύνολο (με την ίδια πληθικότητα με τους πραγματικούς αριθμούς). Όπως και το \operatorname{F}(f), το \operatorname{J}(f) παραμένει αναλλοίωτο από το f(z), και σε αυτό το σύνολο η επανάληψη είναι απωθητική, που σημαίνει ότι |f(z) - f(w)| > |z - w|  για όλα τα *w* σε μια γειτονιά του *z* (μέσα στο \operatorname{J}(f)). Αυτό σημαίνει ότι η f(z) συμπεριφέρεται χαοτικά στο σύνολο Julia. Αν και υπάρχουν σημεία στο σύνολο Julia των οποίων η ακολουθία των επαναλήψεων είναι πεπερασμένη, υπάρχει μόνο ένας μετρήσιμος αριθμός τέτοιων σημείων (και αποτελούν ένα απειροελάχιστο μέρος του συνόλου Julia). Οι ακολουθίες που παράγονται από σημεία εκτός αυτού του συνόλου συμπεριφέρονται χαοτικά, ένα φαινόμενο που ονομάζεται *προσδιοριστικό χάος*.

Υπήρξε εκτεταμένη έρευνα σχετικά με το σύνολο Φατού και το σύνολο Τζούλια των επαναλαμβανόμενων ορθολογικών συναρτήσεων, γνωστών ως ορθολογικοί χάρτες. Για παράδειγμα, είναι γνωστό ότι το σύνολο Φατού ενός ορθολογικού χάρτη έχει είτε 0, 1, 2 είτε άπειρες συνιστώσες. Κάθε συνιστώσα του συνόλου Φατού ενός ορθολογικού χάρτη μπορεί να ταξινομηθεί σε μία από τέσσερις διαφορετικές κατηγορίες.

## Αντίστοιχες περιγραφές του συνόλου Julia
* \operatorname{J}(f) είναι το μικρότερο κλειστό σύνολο που περιέχει τουλάχιστον τρία σημεία και είναι πλήρως αναλλοίωτο υπό *f*.
*  \operatorname{J}(f) είναι το κλειστό του συνόλου των απωθούμενων περιοδικών σημείων.
* Για όλα, εκτός από το πολύ δύο σημεία \;z \in X\;, το σύνολο Julia είναι το σύνολο των οριακών σημείων της πλήρους αντίστροφης τροχιάς \bigcup_n f^{-n}(z). (Αυτό υποδηλώνει έναν απλό αλγόριθμο για την απεικόνιση των συνόλων Julia, βλ. παρακάτω).
* Αν η *f* είναι μια ολόκληρη συνάρτηση, τότε \operatorname{J}(f) είναι το όριο του συνόλου των σημείων που συγκλίνουν στο άπειρο υπό επανάληψη.
* Αν το *f* είναι ένα πολυώνυμο, τότε το \operatorname{J}(f) είναι το όριο του γεμάτου συνόλου Τζούλια, δηλαδή τα σημεία των οποίων οι τροχιές κάτω από επαναλήψεις του *f* παραμένουν περιορισμένες.

## Ιδιότητες του συνόλου Julia και του συνόλου Φατού
Το σύνολο Julia και το σύνολο Φατού της *f* είναι και τα δύο πλήρως αναλλοίωτο υπό τις επαναλήψεις της ολομορφικής συνάρτησης *f*:
: f^{-1}(\operatorname{J}(f)) = f(\operatorname{J}(f)) = \operatorname{J}(f),
: f^{-1}(\operatorname{F}(f)) = f(\operatorname{F}(f)) = \operatorname{F}(f).

## Παραδείγματα
Για f(z) = z^{2} το σύνολο Julia είναι ο μοναδιαίος κύκλος και σε αυτό η επανάληψη δίνεται από τον διπλασιασμό των γωνιών (μια πράξη που είναι χαοτική στα σημεία των οποίων το όρισμα δεν είναι ορθολογικό κλάσμα του 2\pi). Υπάρχουν δύο πεδία Φατού: το εσωτερικό και το εξωτερικό του κύκλου, με επανάληψη προς το 0 και το ∞, αντίστοιχα.

Για g(z) = z^{2} - 2 το σύνολο Julia είναι το ευθύγραμμο τμήμα μεταξύ -2 και 2.  Υπάρχει μία περιοχή Φατού τα σημεία που δεν βρίσκονται στο ευθύγραμμο τμήμα επαναλαμβάνονται προς το ∞. (Εκτός από μια μετατόπιση και κλιμάκωση του πεδίου, αυτή η επανάληψη είναι ισοδύναμη με x \to 4(x - \tfrac{1}{2})^{2} στο μοναδιαίο διάστημα, το οποίο χρησιμοποιείται συνήθως ως παράδειγμα χαοτικού συστήματος).

Οι συναρτήσεις *f* και *g* είναι της μορφής z^2 + c, όπου *c* είναι ένας μιγαδικός αριθμός. Για μια τέτοια επανάληψη το σύνολο Julia δεν είναι γενικά μια απλή καμπύλη, αλλά είναι ένα φράκταλ, και για ορισμένες τιμές του *c* μπορεί να πάρει εκπληκτικές μορφές. Δείτε τις παρακάτω εικόνες.

Για ορισμένες συναρτήσεις *f*(*z*) γνωρίζουμε εκ των προτέρων ότι το σύνολο Julia είναι φράκταλ και όχι απλή καμπύλη. Αυτό συμβαίνει λόγω του ακόλουθου αποτελέσματος για τις επαναλήψεις μιας λογικής συνάρτησης:

**Θεώρημα** - Κάθε περιοχή Φατού έχει το ίδιο όριο, το οποίο κατά συνέπεια είναι το σύνολο Julia.

Αυτό σημαίνει ότι κάθε σημείο του συνόλου Julia είναι ένα σημείο συγκέντρωσης για κάθε έναν από τους τομείς Φατού. Επομένως, αν υπάρχουν περισσότεροι από δύο τομείς Φατού, *κάθε* σημείο του συνόλου Julia πρέπει να έχει σημεία περισσότερων από δύο διαφορετικών ανοικτών συνόλων απείρως κοντά, πράγμα που σημαίνει ότι το σύνολο Julia δεν μπορεί να είναι μια απλή καμπύλη. Ένα τέτοιο φαινόμενο συμβαίνει, για παράδειγμα, όταν η *f*(*z*) είναι η επανάληψη Νεύτων για την επίλυση της εξίσωσης \;P(z) := z^n - 1 = 0 ~ : ~ n > 2\;:

:f(z) = z - \frac{P(z)}{P'(z)} = \frac{\;1 + (n-1)z^n\;}{nz^{n-1}} ~ .

Η εικόνα στα δεξιά δείχνει την περίπτωση *n* = 3.

## Τετραγωνικά πολυώνυμα
Ένα πολύ δημοφιλές μιγαδικό δυναμικό σύστημα δίνεται από την οικογένεια των μιγαδικών τετραγωνικών πολυωνύμων, μια ειδική περίπτωση των ορθολογικών χαρτών. Τέτοια τετραγωνικά πολυώνυμα μπορούν να διατυπωθούν ως εξής 
: f_c(z) = z^2 + c ~,
όπου *c* είναι μια μιγαδική παράμετρος. Ορίστε κάποιο R > 0 αρκετά μεγάλο ώστε R^2 - R \ge |c|. (Για παράδειγμα, αν } είναι στο σύνολο Μάντελμπροτ, τότε |c| \le 2, οπότε μπορούμε απλά να αφήσουμε R = 2~.) Τότε το συμπληρωμένο σύνολο Julia για αυτό το σύστημα είναι το υποσύνολο του μιγαδικού επιπέδου που δίνεται από τη σχέση
: K(f_c) = \left\{z \in \mathbb C : \forall n \in \mathbb N, |f_c^n(z)| \le R \right\}~,
όπου f_c^n(z) είναι η *n*η επανάληψη της f_c(z). Το σύνολο Julia J(f_c) αυτής της συνάρτησης είναι το όριο της K(f_c).


File:JSr07885.gif|Σύνολα Julia για  z^2 + 0.7885\,e^{ia} ,  όπου το *a* κυμαίνεται από 0 έως 2\pi
File:Julia_circling.ogg|Το βιντεο των συνόλων Julia
Image:Julia-Menge.png|Πλήρες σύνολο Julia για *f<sub>c</sub>*, , όπου *φ* είναι χρυσή τομή
Image:Julia 0.4 0.6.png|Σύνολο Julia για *f<sub>c</sub>*, 
Image:Julia 0.285 0.png|Σύνολο Julia για *f<sub>c</sub>*, 
Image:Julia 0.285 0.01.png|Σύνολο Julia για *f<sub>c</sub>*, 
Image:Julia 0.45 0.1428.png|Σύνολο Julia για *f<sub>c</sub>*, 
Image:Julia -0.70176 -0.3842.png|Σύνολο Julia για *f<sub>c</sub>*, 
Image:Julia -0.835 -0.2321.png|Σύνολο Julia για *f<sub>c</sub>*, 
Image:Julia -0.8 0.156.png|Σύνολο Julia για *f<sub>c</sub>*, 
Image:Julia -0.7269 0.1889.png|Σύνολο Julia για *f<sub>c</sub>*, 
Image:Julia 0 0.8.png|Σύνολο Julia για *f<sub>c</sub>*, 
File:Julia Mandelbrot Relationship.png|Συλλογή συνόλων Julia τοποθετημένων σε ένα πλέγμα 100 × 100 έτσι ώστε το κέντρο κάθε εικόνας να αντιστοιχεί στην ίδια θέση στο μιγαδικό επίπεδο με την τιμή του συνόλου. Όταν τοποθετείται έτσι, η συνολική εικόνα μοιάζει με ένα φωτογραφικό μωσαϊκό που απεικονίζει ένα σύνολο Μάντελμπροτ.


Το επίπεδο παραμέτρων των τετραγωνικών πολυωνύμων&nbsp;- δηλαδή, το επίπεδο των πιθανών τιμών του *c*&nbsp;- δίνει το περίφημο σύνολο Μάντελμπροτ. Πράγματι, το σύνολο Μάντελμπροτ ορίζεται ως το σύνολο όλων των *c* ώστε το J(f_c) να είναι συνδεδεμένο. Για παραμέτρους εκτός του συνόλου Μάντελμπροτ, το σύνολο Τζούλια είναι ένας χώρος Κάντορ: σε αυτή την περίπτωση αναφέρεται μερικές φορές ως **σκόνη Φατού**.

Σε πολλές περιπτώσεις, το σύνολο Julia του *c* μοιάζει με το σύνολο Μάντελμπροτ σε επαρκώς μικρές γειτονιές του *c*. Αυτό ισχύει, ειδικότερα, για τις λεγόμενες παραμέτρους Μισιούρεβιτς, δηλαδή για τις παραμέτρους *c* για τις οποίες το κρίσιμο σημείο είναι προ-περιοδικό. Για παράδειγμα:
* Στο *c* = *i*, το κοντύτερο, μπροστινό δάχτυλο του μπροστινού ποδιού, το σύνολο Julia μοιάζει με διακλαδισμένο κεραυνό.
* Στο *c* = -2, το άκρο της μακράς αγκαθωτής ουράς, το σύνολο Julia είναι ένα ευθύγραμμο τμήμα.

Με άλλα λόγια, τα σύνολα Julia J(f_c) είναι τοπικά παρόμοια γύρω από τα σημεία Μισιούρεβιτς.

## Γενικεύσεις
Ο ορισμός των συνόλων Τζούλια και Φατού μεταφέρεται εύκολα στην περίπτωση ορισμένων χαρτών των οποίων η εικόνα περιέχει το πεδίο ορισμού τους, κυρίως των μερομορφικών συναρτήσεων και των **χαρτών πεπερασμένου τύπου** του Άνταμ Επστάιν.

Τα σύνολα Julia ορίζονται επίσης συνήθως στη μελέτη της δυναμικής σε πολλές μιγαδικές μεταβλητές.

## Ψευδοκώδικας
Οι παρακάτω εφαρμογές ψευδοκώδικα κωδικοποιούν αυστηρά τις συναρτήσεις για κάθε φράκταλ. Λάβετε υπόψιν σας την υλοποίηση λειτουργιών μιγαδικών αριθμών, ώστε να είναι πιο δυναμικός και επαναχρησιμοποιήσιμος ο κώδικας.

### Ψευδοκώδικας για κανονικά σύνολα Julia
: f(z) = z^2 + c


```perl
R = escape radius  # choose R > 0 such that R**2 - R >= sqrt(cx**2 + cy**2)

for each pixel (x, y) on the screen, do:   
{
    zx = scaled x coordinate of pixel # (scale to be between -R and R)
       # zx represents the real part of z.
    zy = scaled y coordinate of pixel # (scale to be between -R and R)
       # zy represents the imaginary part of z.

    iteration = 0
    max_iteration = 1000
  
    while (zx * zx + zy * zy < R**2  AND  iteration < max_iteration) 
    {
        xtemp = zx * zx - zy * zy
        zy = 2 * zx * zy  + cy 
        zx = xtemp + cx
    
        iteration = iteration + 1 
    }
  
    if (iteration == max_iteration)
        return black;
    else
        return iteration;
}
```


### Ψευδοκώδικας για σύνολα πολλαπλών Julia
: f(z) = z^n + c

```perl
R = escape radius #  choose R > 0 such that R**n - R >= sqrt(cx**2 + cy**2)

for each pixel (x, y) on the screen, do:
{
    zx = scaled x coordinate of pixel # (scale to be between -R and R)
    zy = scaled y coordinate of pixel # (scale to be between -R and R)
  
    iteration = 0
    max_iteration = 1000
  
    while (zx * zx + zy * zy < R**2  AND  iteration < max_iteration) 
    {
        xtmp = (zx * zx + zy * zy) ^ (n / 2) * cos(n * atan2(zy, zx)) + cx;
	    zy = (zx * zx + zy * zy) ^ (n / 2) * sin(n * atan2(zy, zx)) + cy;
	    zx = xtmp;
    
        iteration = iteration + 1
    } 
    if (iteration == max_iteration)
        return black;
    else
        return iteration;
}
```


## Συνάρτηση δυναμικού και πραγματικός αριθμός επανάληψης
Το σύνολο Julia για την f(z) = z^{2} είναι ο μοναδιαίος κύκλος και στο εξωτερικό πεδίο Φατού, η *συνάρτηση δυναμικού* *φ*(*z*) ορίζεται από τη σχέση *φ*(*z*) = log|*z*|. Οι ισοδυναμικές γραμμές αυτής της συνάρτησης είναι ομόκεντροι κύκλοι. Καθώς |f(z)| = |z|^{2} έχουμε

:\varphi(z) = \lim_{k\to\infty} \frac{\log|z_k|}{2^k},

όπου z_k είναι η ακολουθία της επανάληψης που δημιουργείται από το *z*. Για τη γενικότερη επανάληψη f(z) = z^2 + c, έχει αποδειχθεί ότι αν το σύνολο Julia είναι συνδεδεμένο (δηλαδή αν το *c* ανήκει στο (συνηθισμένο) σύνολο Μάντελμπροτ), τότε υπάρχει ένας διχομορφικός χάρτης *ψ* μεταξύ του εξωτερικού πεδίου Φατού και του εξωτερικού του μοναδιαίου κύκλου έτσι ώστε |\psi(f(z))| = |\psi(z)|^{2}.  Αυτό σημαίνει ότι η συνάρτηση δυναμικού στην εξωτερική περιοχή Φατού που ορίζεται από αυτή την αντιστοιχία δίνεται από:

:\varphi(z) = \lim_{k\to\infty} \frac{\log|z_k|}{2^k}. 

Ο τύπος αυτός έχει νόημα και αν το σύνολο Julia δεν είναι συνδεδεμένο, οπότε για όλα τα *c* μπορούμε να ορίσουμε τη συνάρτηση δυναμικού στο πεδίο Φατού που περιέχει το ∞ με αυτόν τον τύπο. Για μια γενική ορθολογική συνάρτηση *f*(*z*) τέτοια ώστε το ∞ να είναι κρίσιμο σημείο και σταθερό σημείο, δηλαδή τέτοια ώστε ο βαθμός *m* του αριθμητή να είναι τουλάχιστον δύο φορές μεγαλύτερος από τον βαθμό *n* του παρονομαστή, ορίζουμε τη *δυνητική συνάρτηση* στο πεδίο Fatou που περιέχει το ∞ από:

:\varphi(z) = \lim_{k\to\infty} \frac{\log|z_k|}{d^k},

όπου *d* = *m* − *n* είναι ο βαθμός της ορθολογικής συνάρτησης. 

Αν το *N* είναι ένας πολύ μεγάλος αριθμός (π.χ. 10100), και αν το *k* είναι ο πρώτος αριθμός επανάληψης τέτοιος ώστε |z_k| > N, έχουμε αυτό

:\frac{\log|z_k|} {d^k} = \frac{\log(N)}{d^{\nu(z)}},

για κάποιο πραγματικό αριθμό \nu(z), ο οποίος θα πρέπει να θεωρηθεί ως ο *πραγματικός αριθμός επανάληψης*, και έχουμε το εξής:

:\nu(z) = k - \frac{\log(\log|z_k|/\log(N))}{\log(d)},

όπου ο τελευταίος αριθμός βρίσκεται στο διάστημα [0,&nbsp;1].

Για επανάληψη προς έναν πεπερασμένο κύκλο έλξης τάξης *r*, έχουμε ότι αν z^* είναι ένα σημείο του κύκλου, τότε f(f(...f(z^*))) = z^* (η *r*-πτυχωτή σύνθεση), και ο αριθμός

:\alpha = \frac{1}{\left |(d(f(f(\cdots f(z))))/dz)_{z=z^*} \right |} \qquad ( > 1)

είναι η έλξη του κύκλου. Αν *w* είναι ένα σημείο πολύ κοντά στο z^* και *w*&prime; είναι το *w* επαναλαμβανόμενο *r* φορές, έχουμε το εξής

:\alpha = \lim_{k\to\infty} \frac.

Επομένως ο αριθμός |z_{kr} - z^*|\alpha^{k} είναι σχεδόν ανεξάρτητος του *k*.. Ορίζουμε τη συνάρτηση δυναμικού στο πεδίο Φατού ως εξής:

:\varphi(z) = \lim_{k\to\infty} \frac{1}{(|z_{kr} - z^*|\alpha^{k})}.

Εάν ε είναι ένας πολύ μικρός αριθμός και *k* είναι ο πρώτος αριθμός επανάληψης τέτοιος ώστε |z_k - z^*|  , έχουμε το εξής

:\varphi(z) = \frac{1}{(\varepsilon \alpha^{\nu(z)})}

για κάποιον πραγματικό αριθμό\nu(z), ο οποίος θα πρέπει να θεωρηθεί ως ο πραγματικός αριθμός επανάληψης, και έχουμε το εξής:

:\nu(z) = k - \frac{\log(\varepsilon/|z_k - z^*|)}{\log(\alpha)}.

Αν η έλξη είναι ∞, δηλαδή ο κύκλος είναι υπερ-ελκυστικός, δηλαδή και πάλι ένα από τα σημεία του κύκλου είναι κρίσιμο σημείο, πρέπει να αντικαταστήσουμε το *α* με

:\alpha = \lim_{k\to\infty} \frac{\log|w' - z^*|}{\log|w - z^*|},

όπου *w*&prime; είναι το *w* επαναλαμβανόμενο *r* φορές και ο τύπος για το *φ*(*z*) από:

:\varphi(z) = \lim_{k\to\infty} \frac{\log(1/|z_{kr} - z^*|)}{\alpha^k}. 

Και τώρα ο πραγματικός αριθμός επανάληψης δίνεται από:

:\nu(z) = k - \frac{\log(\log|z_k - z^*|/\log(\varepsilon))}{\log(\alpha)}.

Για τον χρωματισμό πρέπει να έχουμε μια κυκλική κλίμακα χρωμάτων (κατασκευασμένη μαθηματικά, παραδείγματος χάριν) που να περιέχει χρώματα *Η* αριθμημένα από το 0 έως το *Η*-1 (*Η* = 500, ενδεικτικά). Πολλαπλασιάζουμε τον πραγματικό αριθμό \nu(z) με έναν σταθερό πραγματικό αριθμό που καθορίζει την πυκνότητα των χρωμάτων στην εικόνα και παίρνουμε το ολοκληρωτικό μέρος αυτού του αριθμού modulo *H*.

Ο ορισμός της συνάρτησης δυναμικού και ο τρόπος χρωματισμού μας προϋποθέτουν ότι ο κύκλος είναι ελκτικός, δηλαδή όχι ουδέτερος. Εάν ο κύκλος είναι ουδέτερος, δεν μπορούμε να χρωματίσουμε την περιοχή Fatou με φυσικό τρόπο. Καθώς το τέρμα της επανάληψης είναι μια περιστροφική κίνηση, μπορούμε, λόγου χάρη, να χρωματίσουμε με βάση την ελάχιστη απόσταση από τον κύκλο που αφήνει σταθερό η επανάληψη.

## Γραμμές πεδίου

Σε κάθε πεδίο Φατού (που δεν είναι ουδέτερο) υπάρχουν δύο συστήματα γραμμών ορθογώνια μεταξύ τους: οι *ισοδυναμικές γραμμές* (για τη συνάρτηση δυναμικού ή τον πραγματικό αριθμό επανάληψης) και οι *γραμμές πεδίου*.

Αν χρωματίσουμε το πεδίο Φατού σύμφωνα με τον αριθμό επανάληψης (και *όχι* τον πραγματικό αριθμό επανάληψης \nu(z), όπως ορίστηκε στην προηγούμενη ενότητα), οι ζώνες επανάληψης δείχνουν την πορεία των ισοδυναμικών γραμμών. Αν η επανάληψη είναι προς το ∞ (όπως συμβαίνει με την εξωτερική περιοχή Φατού για τη συνήθη επανάληψη z^{2} + c), μπορούμε εύκολα να δείξουμε την πορεία των γραμμών πεδίου, δηλαδή αλλάζοντας το χρώμα ανάλογα με το αν το τελευταίο σημείο στην ακολουθία της επανάληψης είναι πάνω ή κάτω από τον άξονα *x* (πρώτη εικόνα), αλλά σε αυτή την περίπτωση (πιο συγκεκριμένα: όταν το πεδίο Fatou είναι υπερ-επικρατές) δεν μπορούμε να σχεδιάσουμε τις γραμμές πεδίου συνεκτικά - τουλάχιστον όχι με τη μέθοδο που περιγράφουμε εδώ. Στην περίπτωση αυτή μια γραμμή πεδίου ονομάζεται επίσης εξωτερική ακτίνα.

Έστω *z* ένα σημείο στο ελκτικό πεδίο Φατού. Αν επαναλάβουμε το *z* πολλές φορές, το τέρμα της ακολουθίας επανάληψης είναι ένας πεπερασμένος κύκλος *C*, και το πεδίο Φατού είναι (εξ ορισμού) το σύνολο των σημείων των οποίων η ακολουθία επανάληψης συγκλίνει προς το *C*. Οι γραμμές πεδίου ξεκινούν από τα σημεία του *C* και από το (άπειρο πλήθος) σημείων που επαναλαμβάνονται *προς* ένα σημείο του *C*. Και καταλήγουν στο σύνολο Julia σε σημεία που δεν είναι χαοτικά (δηλαδή δημιουργούν έναν πεπερασμένο κύκλο). Έστω *r* η τάξη του κύκλου *C* (ο αριθμός των σημείων του) και έστω z^* ένα σημείο του *C*. Έχουμε f(f(\dots f(z^*))) = z^* (η r-πολλαπλή σύνθεση) και ορίζουμε τον μιγαδικό αριθμό α ως εξής

:\alpha = (d(f(f(\dots f(z))))/dz)_{z=z^*}.

Αν τα σημεία του *C* είναι z_i, i = 1, \dots, r (z_1 = z^*), το α είναι το γινόμενο των αριθμών *r* f'(z_i). Ο πραγματικός αριθμός 1/|α| είναι η *έλξη* του κύκλου, και η παραδοχή μας ότι ο κύκλος δεν είναι ούτε ουδέτερος ούτε υπερ-ελκυστικός, σημαίνει ότι 1 z^* είναι σταθερό σημείο για το f(f(\dots f(z))), και κοντά σε αυτό το σημείο ο χάρτης f(f(\dots f(z))) έχει (σε σχέση με τις γραμμές πεδίου) χαρακτήρα περιστροφής με το επιχείρημα β του α (δηλαδή \alpha = |\alpha|e^{\beta i}).

Για να χρωματίσουμε το πεδίο Φατού, επιλέξαμε έναν μικρό αριθμό ε και θέσαμε τις ακολουθίες επανάληψης z_k (k = 0, 1, 2, \dots, z_0 = z) να σταματούν όταν |z_k - z^*| , και χρωματίζουμε το σημείο *z* σύμφωνα με τον αριθμό *k* (ή τον πραγματικό αριθμό επανάληψης, αν προτιμάμε έναν ομαλό χρωματισμό). Αν επιλέξουμε μια κατεύθυνση από το z^* που δίνεται από μια γωνία *θ*, η γραμμή πεδίου που εξέρχεται από το z^* σε αυτή τη διεύθυνση αποτελείται από τα σημεία *z* τέτοια ώστε το επιχείρημα *ψ* του αριθμού z_k - z^* να ικανοποιεί τη συνθήκη ότι

:\psi - k\beta = \theta \mod \pi. \, 

Πράγματι, αν περάσουμε μια ζώνη επανάληψης προς την κατεύθυνση των γραμμών πεδίου (και μακριά από τον κύκλο), ο αριθμός των επαναλήψεων *k* αυξάνεται κατά 1 και ο αριθμός ψ αυξάνεται κατά β, οπότε ο αριθμός \psi - k\beta \mod \pi είναι σταθερός κατά μήκος της γραμμής πεδίου.

Ένας χρωματισμός των γραμμών πεδίου του πεδίου Φατού σημαίνει ότι χρωματίζουμε τους χώρους μεταξύ ζευγών γραμμών πεδίου: επιλέγουμε έναν αριθμό κανονικά τοποθετημένων κατευθύνσεων που ξεκινούν από το z^*, και σε κάθε μία από αυτές τις κατευθύνσεις επιλέγουμε δύο κατευθύνσεις γύρω από αυτή τη διεύθυνση. Καθώς μπορεί να συμβεί οι δύο γραμμές πεδίου ενός ζεύγους να μην καταλήγουν στο ίδιο σημείο του συνόλου Julia, οι χρωματισμένες γραμμές πεδίου μας μπορούν να διακλαδίζονται (ατελείωτα) στην πορεία τους προς το σύνολο Julia. Μπορούμε να χρωματίσουμε με βάση την απόσταση από την κεντρική γραμμή της γραμμής πεδίου, και μπορούμε να αναμείξουμε αυτόν τον χρωματισμό με τον συνηθισμένο χρωματισμό. Τέτοιες εικόνες μπορούν να είναι πολύ διακοσμητικές (δεύτερη εικόνα).

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

## Σχεδιάζοντας το σύνολο Julia

Μέθοδοι : 
* Μέθοδος εκτίμησης απόστασης για το σύνολο Julia (DEM/J)
* Μέθοδος αντίστροφης επανάληψης (IIM)

### Χρησιμοποιώντας την αντίστροφη επανάληψη (IIM)
 
 

Όπως αναφέρθηκε παραπάνω, το σύνολο Julia μπορεί να βρεθεί ως το σύνολο των οριακών σημείων του συνόλου των προ-εικόνων (ουσιαστικά) οποιουδήποτε δεδομένου σημείου. Έτσι μπορούμε να προσπαθήσουμε να σχεδιάσουμε το σύνολο Julia μιας δεδομένης συνάρτησης ως εξής. Ξεκινάμε με οποιοδήποτε σημείο *z* που γνωρίζουμε ότι βρίσκεται στο σύνολο Julia, όπως ένα απωθητικό περιοδικό σημείο, και υπολογίζουμε όλες τις προ-εικόνες του *z* κάτω από κάποια υψηλή επανάληψη f^n της *f*.

Δυστυχώς, καθώς ο αριθμός των επαναλαμβανόμενων προ-εικόνων αυξάνεται εκθετικά, αυτό δεν είναι εφικτό υπολογιστικά. Ωστόσο, μπορούμε να προσαρμόσουμε αυτή τη μέθοδο, με παρόμοιο τρόπο όπως η μέθοδος του "τυχαίου παιγνίου" για επαναληπτικό σύστημα συναρτήσεων. Δηλαδή, σε κάθε βήμα, επιλέγουμε τυχαία μία από τις αντίστροφες εικόνες του *f*.

Για παράδειγμα, για το τετραγωνικό πολυώνυμο *f<sub>c</sub>*, η αντίστροφη επανάληψη περιγράφεται από τη σχέση

:z_{n-1} = \sqrt{z_n - c} .

Σε κάθε βήμα επιλέγεται τυχαία μία από τις δύο τετραγωνικές ρίζες.  

Ας σημειωθεί ότι ορισμένα τμήματα του συνόλου Julia είναι αρκετά δύσκολο να προσεγγιστούν με τον αντίστροφο αλγόριθμο Julia. Για το λόγο αυτό, πρέπει να τροποποιηθεί ο IIM/J ( ονομάζεται MIIM/J) ή να χρησιμοποιηθούν άλλες μέθοδοι για την παραγωγή καλύτερων εικόνων.

### Χρήση DEM/J

File:Demj.jpg| c=-0.74543+0.11301*i
File:Julia dem.png| c=-0.75+0.11*i 
File:Julia dem c=-0.1+0.651.png| c=-0.1+0.651*i
File:Juliasetsdkjulia1.jpg|Σύνολο Julia που σχεδιάστηκε με εκτίμηση της απόστασης, η επανάληψη είναι της μορφής f_c(z) = 1 - z^2 + \frac{z^5}{2 + 4z} + c
File:Juliasetsdkland1.jpg|Τρισδιάστατη απεικόνιση του συνόλου Julia με εκτίμηση της αποστάσεως


Καθώς ένα σύνολο Julia είναι απείρως λεπτό, δεν μπορούμε να το σχεδιάσουμε αποτελεσματικά με αντίστροφη επανάληψη από τα εικονοστοιχεία. Θα εμφανίζεται κατακερματισμένο λόγω της μη πρακτικής εξέτασης απείρως πολλών σημείων εκκίνησης. Δεδομένου ότι ο αριθμός των επαναλήψεων αλλάζει έντονα κοντά στο σύνολο Julia, μια μερική λύση είναι να υπονοήσουμε το περίγραμμα του συνόλου από τα πλησιέστερα χρωματικά περιγράμματα, αλλά το σύνολο θα τείνει να φαίνεται θολό.

Ένας καλύτερος τρόπος για να σχεδιάσετε το σύνολο Julia σε ασπρόμαυρο χρώμα είναι να εκτιμήσετε την απόσταση των εικονοστοιχείων (DEM) από το σύνολο και να χρωματίσετε κάθε εικονοστοιχείο του οποίου το κέντρο βρίσκεται κοντά στο σύνολο. Ο τύπος για την εκτίμηση της απόστασης προκύπτει από τον τύπο της συνάρτησης δυναμικού *φ*(*z*). Όταν οι ισοδυναμικές γραμμές για την *φ*(*z*) βρίσκονται κοντά, ο αριθμός |\varphi'(z)| είναι μεγάλος, και αντίστροφα, επομένως οι ισοδυναμικές γραμμές για τη συνάρτηση \delta(z) = \varphi(z)/|\varphi'(z)| πρέπει να βρίσκονται περίπου κανονικά. Έχει αποδειχθεί ότι η τιμή που βρίσκεται από αυτόν τον τύπο (μέχρι ένα σταθερό παράγοντα) συγκλίνει προς την πραγματική απόσταση για το z που συγκλίνει προς το σύνολο Julia.

Υποθέτουμε ότι το *f*(*z*) είναι ορθολογικό, δηλαδή f(z) = p(z)/q(z) όπου *p*(*z*) και *q*(*z*) είναι μιγαδικά πολυώνυμα βαθμού *m* και *n*, αντίστοιχα, και πρέπει να βρούμε την παράγωγο των παραπάνω εκφράσεων για το *φ*(*z*). Και επειδή μεταβάλλεται μόνο το z_k, πρέπει να υπολογίσουμε την παράγωγο z'_k του z_k ως προς *z*. Καθώς όμως z_k = f(f(\cdots f(z))) (η *k*-πολλαπλή σύνθεση), το z'_k είναι το γινόμενο των αριθμών f'(z_k), και η ακολουθία αυτή μπορεί να υπολογιστεί αναδρομικά με z'_{k+1} = f'(z_k)z'_k, ξεκινώντας με z'_0 = 1 (*πριν* τον υπολογισμό της επόμενης επανάληψης z_{k+1} = f(z_k)).

Για επανάληψη προς το ∞ (πιο συγκεκριμένα όταν , έτσι ώστε το ∞ να είναι ένα υπερ-ελκυστικό σταθερό σημείο), έχουμε

:|\varphi'(z)| = \lim_{k\to\infty} \frac{|z_k|d^k},

(*d* = *m* - *n*) και συνεπώς:

:\delta(z) = \varphi(z)/|\varphi'(z)| = \lim_{k\to\infty} \log|z_k||z_k|/|z'_k|. \, 

Για επανάληψη προς έναν πεπερασμένο ελκυστικό κύκλο (που δεν είναι υπερ-ελκυστικός) που περιέχει το σημείο z* και έχει τάξη *r*, έχουμε

:|\varphi'(z)| = \lim_{k\to\infty} |z'_{kr}|/(|z_{kr} - z^*|^{2}\alpha^{k}), \, 

και συνεπώς: 

:\delta(z) = \varphi(z)/|\varphi'(z)| = \lim_{k\to\infty} |z_{kr} - z^*|/|z'_{kr}|. \, 

Για έναν εξαιρετικά ελκυστικό κύκλο, ο τύπος έχει ως εξής:

:\delta(z) = \lim_{k\to\infty} \log|z_{kr} - z^*|^2/|z'_{kr}|. \, 

Υπολογίζουμε αυτόν τον αριθμό όταν σταματά η επανάληψη. Σημειώστε ότι η εκτίμηση της απόστασης είναι ανεξάρτητη από την έλξη του κύκλου. Αυτό σημαίνει ότι έχει νόημα για υπερβατικές συναρτήσεις "άπειρου βαθμού" (π.χ. sin(*z*) και tan(*z*)).

Εκτός από τη σχεδίαση του ορίου, η συνάρτηση απόστασης μπορεί να εισαχθεί ως 3η διάσταση για τη δημιουργία ενός στερεού φρακταλικού τοπίου.

## Δημοσιεύσεις
* 
*  &nbsp; 
* First appeared in as a  available as 
* 
* 
* 

## Δείτε επίσης
*Σύνολο Μάντελμπροτ
*Απολλώνιο πρόβλημα
*Οικοδομήσιμο Σύμπαν
*Σύνολο Κάντορ

## Παραπομπές


Κατηγορία:Φράκταλ
Κατηγορία:Μαθηματική λογική
Κατηγορία:Φιλοσοφία των μαθηματικών
Κατηγορία:Ιστορία των Μαθηματικών