Science Fair Projects Ideas - Möbius transformation

All Science Fair Projects

      

Science Fair Project Encyclopedia for Schools!

  Search    Browse    Forum  Coach    Links    Editor    Help    Tell-a-Friend    Encyclopedia    Dictionary     

Science Fair Project Encyclopedia

For information on any area of science that interests you,
enter a keyword (eg. scientific method, molecule, cloud, carbohydrate etc.).
Or else, you can start by choosing any of the categories below.

Möbius transformation

Möbius transformations should not be confused with the Möbius transform.
Contents

Geometry

In mathematics, a Möbius transformation, also called a homographic function, is a conformal mapping that is a bijection on the extended complex plane (that is, the complex plane augmented by the point at infinity, written ∞.) It is named in honor of August Ferdinand Möbius.

The general formula is given by

w = \frac{a z + b}{c z + d}

almost everywhere where a, b, c, d are any complex numbers satisfying ad - bc ≠ 0. There are two special cases not covered by the formula above:

  • the point z = -d/c\, is mapped to w=\infty
  • the point z=\infty is mapped to w = a/c\,

We can have Möbius transformations over the real numbers, as well as for the complex numbers. In both cases, we need to augment the domain with a point at infinity.

It can be shown that the inverse and composition of two Möbius transformations are similarly defined, and so the Möbius transformations form a group under composition - called the Möbius group.

The geometric interpretation of the Möbius group is that it is the group of automorphisms of the Riemann sphere. The bilinear transform is a special case of a Möbius transformation.

Any Möbius map can be composed from the elementary transformations - dilations, translations and inversions. If we define a line to be a circle passing through infinity, then it can be shown that a Möbius transformation maps circles to circles, by looking at each elementary transformation.

The Möbius transformation cross-ratio preservation theorem states that the cross-ratio

\frac{(w_1-w_2)(w_3-w_4)}{(w_1-w_4)(w_3-w_2)} = \frac{(z_1-z_2)(z_3-z_4)}{(z_1-z_4)(z_3-z_2)}

is invariant under a Möbius transformation that maps from z to w.

Equations

The transformation

w = \frac{a z + b}{c z + d}

can be usefully expressed as a matrix

\mathfrak{H} = \begin{pmatrix} a & b \\ c & d \end{pmatrix}

In this form, the matrix may be multiplied by any scalar λ and still represent the same transformation. This means that a Möbius transformation on C therefore has six real degrees of freedom.The matrix view of a Möbius transformation corresponds to a projectivity on the projective line over C.

Composition

Let \mathfrak{H}_1, \mathfrak{H}_2 be two Möbius transformations:

\mathfrak{H}_1 = \begin{pmatrix} a_1 & b_1 \\ c_1 & d_1 \end{pmatrix}, \;\;   \mathfrak{H}_2 = \begin{pmatrix} a_2 & b_2 \\ c_2 & d_2 \end{pmatrix}

If these transformations are carried out in succession, first z_1 = \mathfrak{H}_1(z) then Z = \mathfrak{H}_2(z_1) to obtain Z = \mathfrak{H}_2[\mathfrak{H}_1(z)], the result can be readily seen to be another Möbius transformation \mathfrak{H}_3 which appears as the product of the two matrices \mathfrak{H}_1, \mathfrak{H}_2

\mathfrak{H}_3 = \mathfrak{H}_2 \mathfrak{H}_1 =  \begin{pmatrix} a_2 a_1 + b_2 c_1 & a_2 b_1 + b_2 d_1 \\  c_2 a_1 + d_2 c_1 & c_2 b_1 + d_2 d_1 \end{pmatrix}

Thus, Möbius transformations form a group.

Inversion

The inverse of a Möbius transformation \mathfrak{H} can be derived as

z = \mathfrak{H}^{-1}(Z) = \frac{d Z - b}{-c Z + a}

and so

\mathfrak{H}^{-1} = \begin{pmatrix} \;\;d & -b \\ -c & \;\;a \end{pmatrix}

Fixed points, characteristic constant

Any Möbius transformation \mathfrak{H} which is not the identity mapping will have two fixed points γ12, invariant under transformation by \mathfrak{H}. Note that the fixed points are counted here "with multiplicity" so the fixed points may coincide. In this case the Möbius transformation is also called a parabolic transform. Either or both of these fixed points may be the point at infinity: this will happen when c = 0. If this is the case, then the transformation will be an affine transformation (some combination of rotation, dilation, and translation). If both fixed points are at infinity, then the transformation is a pure translation with parameters a = λ,b = λΔ,c = 0,d = λ i.e. the map \gamma \mapsto \gamma+\Delta.

Determination of fixed points

The fixed points can be derived from the fixed point equation

\gamma={{a\gamma +b}\over {c\gamma +d}}

as the two roots

\gamma = \frac{(a - d) \pm \sqrt{(a - d)^2 + 4 c b}}{2 c}

of the quadratic equation

c \gamma^2 - (a - d) \gamma - b = 0 \ ,

which follows from the fixed point equation by multiplying both sides with the denominator cγ + d and collecting equal powers of γ. Note that the quadratic equation degenerates into a linear equation if c = 0, this corresponds to the situation that one of the fixed points is the point at infinity. In this case the second fixed point is finite if a-d \ne 0 otherwise the point at infinity is a fixed point "with multiplicity two" (the case of a pure translation).

Construction of Möbius transformations with prescribed fixed points

Non-parabolic case:

Let us first discuss the case where the transformation has two different fixed points which are finite.

A Möbius \mathfrak{H} transformation is uniquely defined by its two fixed points γ12 and by its so-called characteristic constant k:

\mathfrak{H} = \begin{pmatrix}   k \gamma_2 - \gamma_1 & (1 - k) \gamma_1     \gamma_2 \\   k          - 1        &         \gamma_2 - k \gamma_1 \end{pmatrix}

All transformations with the same characteristic constant are similar. The above transform can be written in a normalized form by multiplying each entry with the (non-zero) scalar c={{1}\over{\gamma_2-\gamma_1}}. The resulting matrix {\mathfrak{H}}' representing the same Möbius transformation has matrix trace equal to k+1 and determinant equal to k. This implies that the characteristic polynomial of the matrix {\mathfrak{H}}' has roots \lambda_1, \, \lambda_2 equal to \,\lambda_1=k, \, \lambda_2=1. Thus the characteristic constant k coincides with one of the two different ratios of eigenvalues \lambda_2 \over \lambda_1 of the matrix representing the transformation (note that each ratio is invariant under multiplication of matrices with non-zero scalars). This property shows that the characteristic constant is preserved under a similarity transformation (corresponding to a conjugation of the matrix representing the transformation).

Every transformation having two different fixed points is similar to some particular linear transformation having one fixed point at infinity and another at 0, i.e a map \gamma \mapsto \lambda \gamma. In the degenerate case k = 1 the transform defined above is the identity transform.

Parabolic case:

If there is only one (finite) fixed point γ then the Möbius transformation is of the form

\mathfrak{H} = \begin{pmatrix}    1+\gamma k    & - k \gamma^2  \\    k             &   1-\gamma k  \end{pmatrix}

where k plays the role of the characteristic constant.

Geometric interpretation of the characteristic constant

The following picture depicts the two fixed points of a Möbius transformation in the non-parabolic case:


The characteristic constant can be expressed in terms of its logarithm:

e^{\rho + \alpha i} = k \;

When expressed in this way, ρ becomes an expansion factor. It indicates how repulsive the fixed point γ1 is, and how attractive γ2 is. If ρ = 0, then the fixed points are neither attractive nor repulsive but indifferent, and the transformation is said to be elliptical. These transformations tend to move all points in circles around the two fixed points . If one of the fixed points is at infinity, the this is equivalent to doing an affine rotation around a point.



α is a rotation factor, indicating to what extent the transform rotates the plane anti-clockwise about γ1 and clockwise about γ2. If α is zero (or a multiple of ), then the transformation is said to be hyperbolic. These transformations tend to move points along circular paths from one fixed point toward the other.



If both ρ and α are nonzero, then the transformation is said to be loxodromic. These transformations tend to move all points in S-shaped paths from one fixed point to the other.


Iterating a transformation

If a transformation \mathfrak{H} has fixed points γ12, and expansion and rotation factors ρ and α, then \mathfrak{H}' = \mathfrak{H}^n will have γ1' = γ12' = γ2,ρ' = nρ,α' = nα.

This can be used to continuously iterate a transformation.

These images show three points (red, blue and black) continuously iterated under transformations with various characteristic constants



Poles of the transformation

The point

z_\infty = - \frac{d}{c}

is called the pole of \mathfrak{H}; it is that point which is transformed to the point at infinity under \mathfrak{H}.

The inverse pole

Z_\infty = \frac{a}{c}

Is that point to which the point at infinity is transformed. The point midway between the two poles is always the same as the point midway between the two fixed points:

\gamma_1 + \gamma_2 = z_\infty + Z_\infty

These four points are the vertices of a parallelogram which is sometimes called the characteristic parallelogram of the transformation.

A transform \mathfrak{H} can be specified with two fixed points γ12 and the pole z_\infty.

\mathfrak{H} = \begin{pmatrix}   Z_\infty & - \gamma_1 \gamma_2 \\   1        & - z_\infty \end{pmatrix}, \;\;   Z_\infty = \gamma_1 + \gamma_2 - z_\infty

This allows us to derive a formula for conversion between k and z_\infty given γ12:

z_\infty = \frac{k \gamma_1 - \gamma_2}{1 - k}
k  = \frac{\gamma_2 - z_\infty}{\gamma_1 - z_\infty} = \frac{Z_\infty - \gamma_1}{Z_\infty - \gamma_2} = \frac {a - c \gamma_1}{a - c \gamma_2}

Which reduces down to

k = \frac{(a + d) + \sqrt {(a - d)^2 + 4 b c}}{(a + d) - \sqrt {(a - d)^2 + 4 b c}}

The last expression coincides with one of the (mutually reciprocal) eigenvalue ratios \lambda_1\over \lambda_2 of the matrix

\mathfrak{H} = \begin{pmatrix}   a & b \\   c & d \end{pmatrix}

representing the transform (compare the discussion in the preceding section about the characteristic constant of a transformation). Its characteristic polynomial is equal to

\mbox{det} (\lambda I_2- \mathfrak{H}) =\lambda^2-\mbox{tr} \mathfrak{H}\,\lambda+ \mbox{det} \mathfrak{H} =\lambda^2-(a+d)\lambda+(ad-bc)

which has roots

\lambda_{i}=\frac{(a + d) \pm \sqrt {(a - d)^2 + 4 b c}}{2}=\frac{(a + d) \pm \sqrt {(a + d)^2 - 4(ad-b c)}}{2} \ .

Specifying a transformation by three points

Direct approach

Any set of three points

Z_1 = \mathfrak{H}(z_1), \;\;   Z_2 = \mathfrak{H}(z_2), \;\;   Z_3 = \mathfrak{H}(z_3)

uniquely defines a transformation \mathfrak{H}. To calculate this out, it is handy to make use of a transformation that is able to map three points onto (0,0), (1, 0) and the point at infinity.

\mathfrak{H}_1 = \begin{pmatrix}  \frac{z_2 - z_3}{z_2 - z_1} & -z_1 \frac{z_2 - z_3}{z_2 - z_1}  \\ 1 & -z_3  \end{pmatrix}, \;\; \mathfrak{H}_2 = \begin{pmatrix}  \frac{Z_2 - Z_3}{Z_2 - Z_1} & -Z_1 \frac{Z_2 - Z_3}{Z_2 - Z_1}  \\ 1 & -Z_3  \end{pmatrix}

One can get rid of the infinities by multiplying out by z2 - z1 and Z2 - Z1 as previously noted.

\mathfrak{H}_1 = \begin{pmatrix}  z_2 - z_3 & z_1 z_3 - z_1 z_2  \\ z_2 - z_1 & z_1 z_3 - z_3 z_2  \end{pmatrix} , \;\; \mathfrak{H}_2 = \begin{pmatrix}  Z_2 - Z_3 & Z_1 Z_3 - Z_1 Z_2  \\ Z_2 - Z_1 & Z_1 Z_3 - Z_3 Z_2  \end{pmatrix}

The matrix \mathfrak{H} to map z1,2,3 onto Z1,2,3 then becomes

\mathfrak{H} = \mathfrak{H}_2^{-1} \mathfrak{H}_1

You can multiply this out, if you want, but if you are writing code then it's easier to use temporary variables for the middle terms.

Alternate method using cross-ratios of point quadruples

This construction exploits the fact (mentioned in the first section) that the cross-ratio

\mbox{cr}(z_1,z_2,z_3,z_4)= {{(z_1-z_3)(z_2-z_4)}\over{(z_1-z_4)(z_2-z_3)}}

is invariant under a Möbius transformation mapping a quadruple (z1,z2,z3,z4) to (w1,w2,w3,w4) via w_i=\mathfrak{H}z_i. If \mathfrak{H} maps a triple (z1,z2,z3) of pairwise different zi to another triple (w1,w2,w3), then the Möbius transformation \mathfrak{H} is determined by the equation

\mbox{cr}(\mathfrak{H}(z),w_1,w_2,w_3)=\mbox{cr}(z,z_1,z_2,z_3)

or written out in concrete terms:

{{(\mathfrak{H}(z)-w_2)(w_1-w_3)} \over{(\mathfrak{H}(z)-w_3)(w_1-w_2)}} ={{(z-z_2)(z_1-z_3)}\over{(z-z_3)(z_1-z_2)}}\ .

The last equation can be transformed into

{{\mathfrak{H}(z)-w_2} \over{\mathfrak{H}(z)-w_3}} ={{(z-z_2)(w_1-w_2)(z_1-z_3)}\over{(z-z_3)(w_1-w_3)(z_1-z_2)}} \ .

Solving this equation for \mathfrak{H}(z) one obtains the sought transformation.

References

Not to be confused with:

See also

External link

A java applet allowing you to specify a transformation via its fixed points and so on may be found here.


This page contains material from this article and this article at PlanetMath, used under the GFDL by permission.

Last updated: 06-01-2005 22:11:37
10-26-2009 08:16:03
The contents of this article is licensed from www.wikipedia.org under the GNU Free Documentation License. Click here to see the transparent copy and copyright details
Science kits, science lessons, science toys, maths toys, hobby kits, science games and books - these are some of many products that can help give your kid an edge in their science fair projects, and develop a tremendous interest in the study of science. When shopping for a science kit or other supplies, make sure that you carefully review the features and quality of the products. Compare prices by going to several online stores. Read product reviews online or refer to magazines.

Start by looking for your science kit review or science toy review. Compare prices but remember, Price $ is not everything. Quality does matter.
Science Fair Coach
What do science fair judges look out for?
ScienceHound
Science Fair Projects for students of all ages
All Science Fair Projects.com Site
All Science Fair Projects Homepage
Search | Browse | Links | From-our-Editor | Books | Help | Contact | Privacy | Disclaimer | Copyright Notice