dynamicfractalviewer
A Dynamic Fractal Viewer, Automatic Expression supported
view repo
A lot of formal and informal recreational study took place in the fields of Meromorphic Maps, since Mandelbrot popularized the map z <- z^2 + c. An immediate generalization of the Mandelbrot z <-z^n + c also known as the Multibrot family were also studied. In the current paper, general truncated polynomial maps of the form z <- ∑_p>=2 a_px^p +c are studied. Two fundamental properties of these polynomial maps are hereby presented. One of them is the existence of shape preserving transformations on fractal images, and another one is the existence of embedded Multibrot fractals inside a polynomial fractal. Any transform expression with transcendental terms also shows embedded Multibrot fractals, due to Taylor series expansion possible on the transcendental functions. We present a method by which existence of embedded fractals can be predicted. A gallery of images is presented alongside to showcase the findings.
READ FULL TEXT VIEW PDFA Dynamic Fractal Viewer, Automatic Expression supported
Mandelbrot popularized the map which bears his name ([1],[5]). The Mandelbrot set , comprise of points ‘’, such that , repeated infinite time, upon itself, remains bounded.
The iteration is defined as in [1] :-
Initialize current point and .
Use the Transform map to generate .
Repeat the procedure (2), to get until or is the maximum number of iteration.
If iterations expired, the original point .
Else .
Take another point in plane, and start again from step (1), till all the points in the plane are exhausted.
When is applied on itself once it can be said . Making things general, an ‘’ times repetition of upon itself would be termed as .
Then, Formally:-
(1.1) |
A plot of the set yields outstanding visual imagery, and patterns.
A generic way to produce any set of this kind would be then to have a transform function such as to:-
(1.2) |
Clearly then, is a very important parameter. However the removal of can be easily done with:-
(1.3) |
For in case of Mandelbrot set it can be easily shown that any point would escape to infinity, that means:-
so for , k is actually omitted.
The images generated by equation (1.3) are called “Escape Time Fractals”, also known as “orbits” fractals.
We have modified an open source ‘Mandelbrot Set Viewer’ computer program, and incorporated the ability to type in mathematical expressions for the function in equation (1.3). The source code is available here.
There is an expression handler in the program that can handle many standard mathematical functions. Given the expression, the program would evaluate the expression dynamically, to generate the ‘escape time fractal’ for the function . The program is capable of:-
zoom-in and zoom-out.
modifying the radii of convergence ‘’, to be expressed as a power of ‘’.
modifying number of iterations to perform as ‘infinity’.
saving the image as ‘png’ format image file, that is how illustrations of this paper were generated.
plotting the set in Colors, to be discussed more in next.
From a computational point of view, there is no way to perform the iteration of equation (1.3) up-to infinity. Hence, a maximum ‘’ number of iteration is performed. Using the custom software we can manipulate ‘’, . If a point exceeds the limit at iteration number , then a ratio is defined. Clearly means, the function remained bounded till the ‘’th iteration. Based upon the value of , the color of the point is chosen, among a fixed set of colors. Though a point might show , that does not imply that, it would remain so for an increased . Also, another point having current might get into the set, for an increased . The coloring can be formally defined as a mapping:-
(1.4) |
The custom application which produced the illustrations of present paper is capable of plotting in any integer number of colors, and can handle any mathematical expression and function. The color codes are dependent on the color map used. The program uses equation (1.5) for coloring.
(1.5) |
In equation (1.5) ‘’ is the total number of colors in the color map. The default value for color maps used is , and the default color is Gray Scale. The user of the software can choose from a list of color maps, or can create his own color map.
As the software can evaluate almost any mathematical expressions to generate the fractals of type equation (1.3), we noticed certain interesting findings.
The most interesting phenomenon we noticed is the existence of ‘embedded’ fractals while using Transcendental expressions, that is fractals, which are created by expressions with transcendental functions in it. The first observation was ‘’ fractal ( figure 8 on page 8 )contains ‘’ or Mandelbrot fractals inside it (figure 18 on page 18).
There are Mandelbrot fractals inside ‘’ fractal. Figure 19 on page 19 illustrates a fractal, and figure 20 on page 20 depicts the embedded Mandelbrot fractals inside it.
The second observation was, that the Mandelbrot set is capable of ‘Shape Preserving Transformations’, like rotation, scaling, and translation, just by multiplying the term by different (complex) factors. We have established a ‘generic’ way to transform polynomial fractals without a term.
Equation (1.1) can be further generalized into:-
(1.6) |
where is the set of Integers.
Further generalization can be done by:-
(1.7) |
We know that the only shape preserving transforms are:-
Translation
Rotation
Scaling
In this section we review the integer multibrot map (1.6), under the influence of shape preserving transformations. Our observation is that Standard shape preserving transformations works on the function family of ‘Multibrots’ with positive ‘n’ (equation (1.6) , ). We note that the family with does not preserve shape under rotation, simply because the transform cease to become a function :-
The expression is not a function, and therefore no 1-1 relationship between the old point, and mapped point can be established.
We show that the point which is now current transformed ‘’ , has the same trajectory of the earlier ‘’. So, if ‘’ was included in the set, so would be ‘’ with modification in ’s value, and, if ‘’ was not included in the set, so would not be ‘’. That would establish the transformation properties.
A variable transformation of the form translates a Multibrot Map.
Translation can be achieved by :-
(2.1) |
The result is the set image shifted co-ordinates to .
One point to be noted here is that as the radii of convergence is , if the point resides further than , then we need to increase the radii appropriately, for seeing the whole set’s image.
To prove this we note that all points which were mapped as according to equation (1.6), are now mapped as in the new mapping with equation (2.1). That would mean that if for equation (1.6), then for equation (2.1), and vice versa. We take a point and hence .
Clearly at 1st level the original transformation can be written as:-
(2.2) |
Compare with 1st level transformation using equation (2.1):-
where and signifies the transformed .
So we have established that:-
(2.3) |
Now, let it be true for some . Then we have:-
We have:-
(2.4) |
and
Replacing the value of we have:-
hence:-
(2.5) |
Transformation rotates a Multibrot Map.
A “clockwise” rotation of angle can be achieved by:-
(2.6) |
To prove that we note again:- and .
We have:-
Hence, using equation (2.2)
(2.7) |
Assume now that for some
Clearly we would have:-
Then, using equation (2.4) we get:-
(2.8) |
Using equation (2.7) as basis and having equation (2.8) implies that, using principle of mathematical induction, we can say all points which were mapped as according to equation (1.6), are now mapped as in the new mapping with equation (2.6). That establishes the rotational property.
∎
Note that for “Multibrots” with this formula holds true. As the term , this rotation would be in “anti-clock-wise” direction. For , this rotation is in “clockwise” direction. Also note that, for “Multibrots”with fractional , rotation leads to altogether different type of images. They are not invariant under rotation.
Transformation scales a Multibrot Map
A Scaling of factor can be achieved by :-
(2.9) |
The resulting figure is a scaled version of the original Map, by a factor of .
To show this, we start with and as initial points.
Clearly then:-
hence, using equation (2.2)
(2.10) |
Now, assume it is true for some
Clearly we can say:-
Then, using equation (2.4) we get:-
(2.11) |
Using equation (2.10) as basis and having equation (2.11) implies that, using principle of mathematical induction, we can say all points which were mapped as according to equation (1.6), are now mapped as in the new mapping with equation (2.9). That establishes the scaling property.
∎
Note that, “Multibrots”with fractional are invariant under scaling.
The most generic shape preserving transformation would be, then
(2.12) |
The transformation of equation (2.12) translates the origin by amount, then rotates the map by , and then scales down the map by .
Polynomials are important due to Taylors Series expansion available for any infinitely differentiable function. The generic formula for a polynomial transform is:-
(2.13) |
We note that while equation (2.1) can be immediately applied to a polynomial, due to variations in power transforms like rotation: equation (2.6) and scaling: equation (2.9) can not be used on equation (2.13). Hence, only with the translation part, we can redefine equation (2.12) for polynomials as:-
(2.14) |
We can not consider polynomial with a linear term because then the results from rotation (2.6) and scaling (2.9) becomes undefined with .
Left Truncated Polynomial.
A Left Truncated polynomial is defined as:-
(2.15) |
We define a truncated polynomial is a polynomial with at least no linear term, (and possibly constant term) to be defined formally in a later section, with equation (2.15). We can drop the requirement of not having a constant term, because in that case , and the equations (2.6) and (2.9) remains defined.
We have found out from equations (2.6) and (2.9) that they affect different powers of ‘z’ differently. But assume, we still want to find a way to rotate, and scale a polynomial fractal image appropriately. For notational convenience, we treat
as a vector, with components
.We note that, to generate an uniform scaling, .
(2.16) |
We note that, to generate an uniform rotation, .
(2.17) |
It is easy to see that the for the ‘effect’ of them all should be equal to some predefined value ‘’,
rewriting, we get:-
In the same way, to get a fixed angle of rotation ‘’ :-
Hence, a properly scaled and rotated original Polynomial fractal transform, with translation, would have components:-
(2.18) |
Using equation (2.18) we can modify and shift any ‘non-infinite’ ( ) truncated ( ) polynomial fractal, such a way that the original shape is preserved. So, the set image remains ‘invariant’ under equation (2.18).
Note that equation (2.18) can handle , but not the scaling part, as ‘’ needed to be adjusted.
We take the equation of :
(2.19) |
The image of equation (2.19) can be seen in figure 9 on page 9. The image of this can be rotated clockwise radian, by multiplying the by and by . Hence,
With the previous set of transforms in mind, we focus on what is observed self similarity in the Multibrot fractals.
The tool we seek is the magnify a small section of the fractal shape, and finding how it should look. We note that this has to deal with the transform scaling ( theorem 2.3 ). But that won’t be all.
The point with respect to which we need to magnify or zoom first needs to be made the origin of the co-ordinate, and then magnified , or scaled. That would make the transform looks like translation followed by scaling. But that is not so.
We established that the direction of the co-ordinate system matters for fractals, as in Multibrot with fractional powers does not show remain same under an induced co-ordinate rotation. Therefore, the transform should also take care of aligning the co ordinates properly, via rotation.
The Zoom Transform.
Zoom of a fractal image can be done using :-
(2.22) |
where is a positive constant factor, a complex number to have the co-ordinate origin shifted to, and denotes the rotation angle.
Existence of Self Similar Embedded Fractals.
A fractal of the type :-
has quasi self similar embedded fractals under magnification iff it is invariant under the equation (2.22) or the zoom transform.
It is to be noted here that all geometric transcendental functions can be expanded in Taylor series, and hence a transform of the form:
where is a geometric transcendental function, qualifies as an infinite polynomial map.
We have noticed an important property of such maps, that they have embedded Multibrot fractals inside small regions within them. Some of them we have already shared in the introduction section. We furnish some more examples here:-
These observations lead us to think that for those maps, the original transform function can probably be approximated by the Multibrot type function. This is the topic of discussion in this section.
What we informally seek is that in the neighborhood of a point, can we replace a function by another much simpler function, because, they behave approximately the same way. This concept is already available as asymptotic bounds proposed by Bachmann [6], and reintroduced by Landau [7], for real functions, with a parameter .
Big Theta : .
Let there exist functions and constants such that :-
then,
(3.1) |
The equation (3.1) is the Bachmann-Landau ‘’ notation, mostly used in computer science (Knuth [8]) for getting approximate bounds for functions (runtime analysis of algorithms). With simple modification, complex analytic functions can be included in notation.
Complex Generalization of Big Theta.
Iff for and constants at neighborhood of a point :-
(3.2) |
then, we define:-
(3.3) |
As both are analytic, then, there are infinite points in the neighborhood of where equation (3.2) will be true. On that neighborhood, will approximate with arbitrary accuracy.
There is a relationship for Complex Polynomials.
Let
Then, if this holds:-
then
In particular:-
and,
We rewrite the polynomial as:-
Taking mods we get:-
which immediately gives the result. ∎
The relation is true over iterations over polynomials.
Let be a polynomial function with
where Then,
We rewrite the polynomial as:-
and then the becomes:-
Taking mod, the result follows immediately.
∎
What we informally seek, is when we “zoom in” to a polynomial fractal, what “inner” fractal type we may find. If we find a fractal of type at a point ‘’ inside a transformation function , then we say that the function must have “dominated” the function at point .
Relationship can be used to find escape time fractals.
Let the transform function be:-
being a polynomial with . Let then:-
and
so that:-
We note that from theorem (3.2) :-
This would imply that:-
which would imply that:-
also
which would imply that:-
combining we get . ∎
Dominance : Relation guarantees a ‘sub-fractal’.
Lets . Then, there for the escape time fractal generated from the transform function:-
there will be a embedded sub fractal of the form:-
In general for , a truncated polynomial (definition 2.1), there would be sub fractals of type .
Theorem 3.3 makes proof of this theorem a triviality. However, another proof can be furnished from the perspective of scaling transform. We demonstrate the proof here.
We assume that:-
If that is the case, then we know that the fractal:-
would be bigger in size than the fractal:-
if the fractal exists due to the scaling theorem (2.9).
That would mean that the fractal is to be bound between the two fractals of the same shape, which is only possible, iff the shape of the fractal is similar to that of , which proves the theorem.
∎
The theorem (3.4) suggests, that there will be Multibrots found in the vicinity of truncated polynomials. In this section we elaborate what observations we have made, and how those observations fit into the property of dominance. We demonstrate simple truncated polynomial fractals.
To illustrate the result, we exhibit two transformation functions. The first transformation function:-
has Mandelbrot () fractals inside it. Figure 9 on page 9 shows the , and inner () in the next image.
The seccond transformation function:-
has Multibrot () fractals inside it. Figure 14 on page 14 shows the , and inner () in the next image.
These two observations shows validity of the dominance property.
The most important result of Dominance lies with the Taylor Series expansion. For a function which is capable of expansion by Taylor Series, “multibrot” type fractals can be generated by removing lower powers, or multiplying or dividing the function by .
We note down the expansion of :-
(3.4) |
and the expansion of :-
(3.5) |
In this subsection for brevity we have used the informal to denote .
Now clearly then due to dominance effect:-
And finally we explain, the initial fractals, which served as the starting point of this paper,in the next subsection.
Figure 17 on page 17 shows a fractal, and figure 18 on page 18 shows embedded Mandelbrots inside it.
(3.6) |
hence,
Now we proceed to explain the observation of ‘’. Figure 19 on page 19 shows a fractal, and figure 20 on page 20 shows embedded Mandelbrots inside it.
(3.7) |
We note down that equation (3.7) contains power of , and hence, we can not directly use the theorem of existence of embedded Multibrots.
Let us define:-
(3.8) |
and
(3.9) |
We state without proof:-
and
Existence of Multibrot inside Fractal of .
Fractal of equation family (3.9) would contain Multibrots inside.
We note that when , the equation (3.8) gets unbounded.
Assume is large enough so that:-
Now, when but and , we have:-
(3.10) |
In that case:-
That approximates into :-
Ignoring the higher powers of we get:-
(3.11) |
And hence equation (3.9) can be written as:-
(3.12) |
Hence we complete the explanations of the observations what we have stated in the beginning of the current section.
This paper started as a celebration of what Mandelbrot’s work has been. We started working on creating a software tool to produce ‘only’ Mandelbrot set in color, and ended up creating a software tool capable of evaluating any Mathematical expression, to generate escape time fractals, in any color. While experimenting with this tool, we came up with heather to unknown, unexplained phenomenon. The phenomenon of shape preserving transform, and the phenomenon of embedded inner fractal images, were so fascinating that we tried to explain them, so as to generate a preliminary theoretical basis. The usage of transformation vector to preserve shape and transform any finite truncated polynomial fractal is proven. The dominance phenomenon we borrowed, and reused from Computer Science, and shown that it is a very powerful mechanism to predict existence of embedded inner fractals. While dominance principle shows promise, the principle is fairly new at least in the fractal geometry domain. Further theoretical research work is needed to take the preliminary concepts discussed in this paper, into mainstream fractal geometry.
We would specially like to thank Carl Johansen, whose open sourced Mandelbrot viewer we borrowed, without which, not a single line of this paper could have been written.
Comments
There are no comments yet.