File:Newton-lplane-Mandelbrot-smooth.jpg
Originaldatei (6.000 × 4.800 Pixel, Dateigröße: 14,19 MB, MIME-Typ: image/jpeg)
Bildtexte
|
BeschreibungNewton-lplane-Mandelbrot-smooth.jpg |
English: Computergraphical study of the critical point 0 of Newton's method for a family of cubic polynomials in the complex -plane. For details see below. Mandelbrot set occurring in the analysis of Newton's method (detail |
|||
Datum | ||||
Quelle | Eigenes Werk | |||
Urheber | Georg-Johann Lay | |||
Genehmigung (Weiternutzung dieser Datei) |
|
Beschreibung[Bearbeiten]
in the complex plane. denotes the Newton operator
For in the black part of the plane, the critical point 0 of does not converge to a zero of . This means that the set of start values for which Newton's method does not converge to a zero of is a set of full measure. The black set is
Coloring[Bearbeiten]
Note: I added the relevant part of the C-source to document what went on.
- getLambdaColor
- gets the color for one pixel
- Nf
- perform one step of Newton's method, returns the next z als well as the value of the hard coded ƒλ(z)
- cpolar
- transform from cartesian coordinates to polar coordinates
- hsv2rgb
- map HSV to RGB color space
- argd[]
- some values that can be passed via command line to fine trim the coloring function.
The coloring function itself cannot be derived or explained at this point. It is based on trial and error, observation, intuition and experience to get a function that yields appealing results.
Resolution (both arithmetic and graphical) and graphics are taken care of by higher level procedures which do not contribute to the basic understanding.
Color getLambdaColor (double x, double y)
{
cplx z = {0,0};
lambda = (cplx) {x+argd[2], y+argd[3]};
int i;
cplx f;
double eps = 1e-8;
double le = 1./log(eps);
for (i=0; i < argd[1]; i++)
{
double v, s, h, b2;
z = Nf(z, &f);
if (isinf (z.x) || isinf (z.y))
return Cwhite;
if (isinf (f.x) || isinf (f.y))
return Cwhite;
b2 = f.x*f.x + f.y*f.y;
if (isinf(b2)) exit(4);
if (b2 < eps*eps)
{
double b = 0.5*log(b2)*le;
if (isinf(b)) b = 2;
b = i-b;
z = cpolar(z);
h = z.y/2/M_PI-.09;
v = b / argd[4];
s = 0.9-0.7*pow(v, 1.5);
if (v >= 1)
{
double q = 1.-log (b-argd[4])/log(argd[1]-argd[4]);
s = 4*q*(1-q);
if (s > 1) s = 1;
if (s< 0) s = 0;
s = 0.2+0.6*pow(s, 10);
v = 1;
}
return hsv2rgb (h, s, v);
}
}
return Cblack;
}
cplx Nf (cplx z, cplx *f)
{
cplx z2 = cprod (z,z);
*f = csum (z2, lambda);
*f = ccsum (-1, *f);
*f = cprod (*f, z);
*f = cdiff (*f, lambda);
cplx N = ccprod (3., z2);
N = csum (N, lambda);
N = ccsum (-1, N);
cplx Z = cprod (z, z2);
Z = csum (Z, Z);
Z = csum (Z, lambda);
return cquot (Z, N);
}
Dateiversionen
Klicke auf einen Zeitpunkt, um diese Version zu laden.
Version vom | Vorschaubild | Maße | Benutzer | Kommentar | |
---|---|---|---|---|---|
aktuell | 23:44, 17. Okt. 2008 | 6.000 × 4.800 (14,19 MB) | Georg-Johann (Diskussion | Beiträge) | {{Information |Description= |Source= |Date= |Author= |Permission= |other_versions= }} | |
20:06, 11. Apr. 2008 | 6.000 × 4.800 (1,96 MB) | Georg-Johann (Diskussion | Beiträge) | {{PD-self}} |
Du kannst diese Datei nicht überschreiben.
Dateiverwendung
Die folgenden 4 Seiten verwenden diese Datei:
Globale Dateiverwendung
Die nachfolgenden anderen Wikis verwenden diese Datei:
- Verwendung auf de.wikipedia.org
- Mandelbrot-Menge
- Wikipedia:Exzellente Bilder
- Wikipedia:Exzellente Bilder/Naturwissenschaften
- Benutzer:Georg-Johann/Bilder
- Wikipedia:Kandidaten für exzellente Bilder/Archiv2008/24
- Datei:Newton-lplane-Mandelbrot-smooth.jpg
- Wikipedia:Eine-Million-Artikel-Seite
- Wikipedia:Eine-Million-Artikel-Seite/Dateien
- Benutzer:Flawed reality
- Verwendung auf en.wikipedia.org
- Verwendung auf fr.wikipedia.org
Metadaten
Diese Datei enthält weitere Informationen, die in der Regel von der Digitalkamera oder dem verwendeten Scanner stammen. Durch nachträgliche Bearbeitung der Originaldatei können einige Details verändert worden sein.
_error | 0 |
---|