Math Equations

Tuesday, April 22, 2014

A parametric curve


In the frame above you should see a figure being generated for you. If you look at it when it's done you'll see that it looks a lot like a pokeball with symmetry added to it, and it's actually a fractal which can be described as follows: there is the unit circle, and each circle contains four smaller circles which are tangent to themselves and the bigger circle.

Are the points of the fractal the points of a parametric curve? At first I was convinced that they are not because the curve must be continuous. However, because all points are connected (by tangency) that's not a problem. Then came the constraint that the particle, which we imagine moving along the curve, must not come to rest, and after a lot of thinking I came up with the motion shown in the figure. The speed is inevitably infinite from time to time because the length of the curve is, and we just have to accept that, but it shouldn't ever be zero because it gets greater as the circles get smaller (by a factor of \(\frac{8}{1 + \sqrt2}\) which you'll be able to verify once I've shown you my definition of the curve).

It should be noted that the curve that I'm going to present does not have a velocity at the points of internal tangency (e.g. \(t = 0\)). Because the direction exists even at those points, I'm sure it is possible to remedy this by using a substitution \(u = g(t)\) where \(g\) is some increasing function that is piecewise differentiable.

OK, let's start. We'll be using complex numbers rather than vectors as values of \(\bar{r}\) because complex numbers are great for doing rotations. Furthermore we choose \(t \in I = [0, 1)\). Define the first iteration of the fractal parametric curve to be the unit circle \(\bar{r}_0(t) = \exp{(2\pi i \cdot t)}\). Then suppose we know how to draw the nth iteration and define the next iteration to be four copies of the current one scaled down to fit inside a new unit circle (which is also part of the new iteration):
\[
\bar{r}_{n+1}(t) =
\begin{cases}
\frac{\bar{r}_n(8t) + \sqrt2}{1 + \sqrt2} \enspace\mathrm{if}\, 0 \leq t < \frac{1}{8} \\
\exp{(2\pi i \cdot 2(t - \frac{1}{8}))} \enspace\mathrm{if}\!\,\, \frac{1}{8} \leq t < \frac{1}{4} \\
i \bar{r}_{n+1}(t-\frac{1}{4}) \enspace\mathrm{if}\!\,\, \frac{1}{4} \leq t < 1
\end{cases}
\]
(The \(1 + \sqrt2\) comes from the ratio between two successive radii and can be found trigonometrically.)

By induction, we know what all iterations look like. To get \(\bar{r}(t)\), we take the limit \(\bar{r}_n(t)\) as \(n \to \infty\). The limit exists because I've set it up so that in most cases \(\bar{r}_n(t) = \bar{r}_{n+1}(t)\) ultimately. In the cases where that's not true, the value will be enclosed by \(n+1\) circles with exponentially decreasing radii, the greatest of which is the unit circle. Because every circle but the unit circle is inside of the previous one, the circles will converge around a single point which will be the limit.

That's it basically. Some interesting points are those which lie on infinitely small circles. For these points the direction of the curve might not exist. One such point is \(t = \frac{4}{7}, \bar{r}(t) = 1 - \sqrt2\). Can you catch 'em all?

Edit: another interesting thing is that \(\bar{r}_n(\frac{4}{7})\) seems to lie on a line for all \(n\). Are there other \(t\) with this property?

No comments:

Post a Comment