Fixpunkte

Evergreens bei der Analyse von Abbildungen sind die Suche nach Fixpunkten oder Fixgeraden.

Definition 1: Fixpunkt

Ein Punkt xx^* heisst Fixpunkt (der Periode 11) von ff, wenn f(x)=xf(x^*) = x^* gilt.

Note 1

Ist xx^* ein Fixpunkt der Periode 11 von ff, dann gilt natürlich auch fn(x)=xf^n(x^*) = x^* für jedes nNn\in\mathbb{N}. Ferner: Seien nNn\in\mathbb{N} und xx^* ein Fixpunkt der Periode kNk\in\mathbb{N} mit nmodk0n\mod{k}\equiv0, dann ist xx^* auch ein Fixpunkt der Periode nn.

Example 1
  • f(x)=xf(x) = \sqrt{x} hat die Fixpunkte 00 und 11.
  • g(x)=2x+1g(x) = 2x+1 hat den Fixpunkt x=1x = -1.
  • h(x)=x2h(x) = x-2 hat keinen Fixpunkt.
Definition 2: Fixpunkt-Ordnung

Ein Punkt xx^* heisst Fixpunkt der Ordnung kk von ff, wenn

fk(x)=xf^k(x^*) = x^*

gilt.

Example 2

Wir betrachten

f(x)=x+1.f(x) = -x+1.

Dann können wir leicht nachrechnen, dass x=0.25x = 0.25 ein Fixpunkt der Periode 22 von ff ist.

Exercise 1: Fixpunkte

Gibt es noch weitere Fixpunkte?

Solution

Der Fixpunkt der Periode 1 ist 12\frac{1}{2}. Für einen Fixpunkt der Periode 2 gilt: (x+1)+1=!x-(-x+1)+1\stackrel{!}{ = }x, also x=xx = x. Das heisst, es gibt unendlich viele Fixpunkte der Periode 2 für ff; jedes xRx\in\mathbb{R}.

Bleibt noch die Frage, wie man Fixpunkte höherer Ordnung findet? Dazu illustrativ ein

Example 3

Wir betrachten

f(x)=2x1.f(x) = 2x-1.

Einen Fixpunkt der Periode 11 finden wir rasch via die Bedingung f(x)=xf(x) = x, also

2x1=x2x-1 = x

was unmittelbar x=1x = 1 liefert.

Um nun Fixpunkte der Periode 22 zu finden berechnen wir

f2(x)=2(2x1)1=4x3f^2(x) = 2(2x-1)-1 = 4x-3

und erhalten

4x3=x.4x-3 = x.

Wir sind enttäuscht, weil x=1x = 1 bereits wieder Lösung ist. Scheinbar gibt es hier keine reinen Fixpunkte der Periode 22. "Keine reinen" deshalb, weil ja ein Fixpunkt erster Periode ebenfalls Fixpunkt nn-ter Periode ist.

Wir vermuten anhand der obigen Rechnung, dass eine affine Funktion keine echten Fixpunkte höherer Periode hat.

Exercise 2: Fixpunkte linearer Funktion

Beweise obige Vermutung.

Solution

Wir rechnen alle Fixpunkte nach.

mx+q=!xx(m1)=qx=q1m\begin{align*} mx+q &\stackrel{!}{ = } x\\ x(m-1) & = -q\\ x & = \frac{q}{1-m} \end{align*}

falls m1m\neq1; für m=1m = 1 ist q=0q = 0 und man hat mit f(x)=xf(x) = x eine Fixpunktgerade. Für 2-Zyklen ist

m(mx+q)+q=!xx(1m2)=q(m+1)x=q(m+1)1m2x=q1m\begin{align*} m(mx+q)+q &\stackrel{!}{ = } x\\ x(1-m^2) & = q(m+1)\\ x & = \frac{q(m+1)}{1-m^2}\\ x & = \frac{q}{1-m} \end{align*}

wobei m±1m\neq\pm1 galt. Dies ist der Fixpunkt der Periode 1 und wir betrachten daher den Fall m=1m = -1 separat. Daraus folgt unter der Fixpunktbedingung sofort x=xx = x, so dass jeder Punkt xx für f(x)=x+qf(x) = -x+q ein 2-Zyklus ist.

Für 3-Zyklen:

m(m(mx+q)+q)+q=!xm(m2x+mq+q)+q=xq(m2+m+1)=x(1m3)x=qm2+m+11m3\begin{align*} m(m(mx+q)+q)+q &\stackrel{!}{ = } x\\ m(m^2x+mq+q)+q & = x\\ q(m^2+m+1) & = x(1-m^3)\\ x & = q\frac{m^2+m+1}{1-m^3} \end{align*}

Berechne (m31)÷(m2+m+1)(m^3-1)\div(m^2+m+1) mit Polynomdivision und sieh dann, dass x=qm2+m+11m3=q1mx = q\frac{m^2+m+1}{1-m^3} = \frac{q}{1-m}. Das heisst, es gibt keine weiteren Fixpunkte.

Die Vermutung war also nicht ganz korrekt: Fixpunkte der Periode 22 liegen auf den Senkrechten zur Winkelhalbierenden.

Definition 3: Attraktor

Ein Fixpunkt xx heisst anziehend oder stabil, wenn für alle Punkte x~\tilde{x} in einer Umgebung von xx gilt, dass

limnfn(y)=x.\lim_{n\to\infty}f^n(y) = x.

Man nennt xx auch etwa Attraktor.

Example 4

Beispielsweise ist für

f(x)=xf(x) = \sqrt{x}

der Punkt x=1x = 1 ein Attraktor.

Exercise 3: Attraktor?

Starte mit den Werten 33 und 0.20.2 und betrachte deren Verhalten unter

limnxn.\lim_{n\to\infty}\sqrt[n]{x}.

Wie könnte man allgemein zeigen, dass 11 ein Attraktor für x\sqrt{x} ist?

Solution

Man wählt eine kleine Umgebung von 11, ±ε\pm\varepsilon, und beobachtet das Geschehen. Einmal iteriert ist f(1+ε)=1εf(1+\varepsilon) = \sqrt{1-\varepsilon} und kk Mal iteriert ergibt 1+ε2k=(1+ε)12k\sqrt[2k]{1+\varepsilon} = (1+\varepsilon)^{\frac{1}{2k}}, was für kk\to\infty gegen den Grenzwert (1+ε)0=1(1+\varepsilon)^{0} = 1 konvergiert. (Limit kommentiert)

Definition 4: Repeller

Ein Fixpunkt xx heisst abstossend oder instabil, wenn für alle Punkte x~x\tilde{x}\neq x in einer Umgebung von xx gilt, dass

limnfn(y)x.\lim_{n\to\infty}f^n(y)\neq x.

Man nennt xx auch Repeller.

Example 5

Für f(x)=x2f(x) = x^2 ist der Fixpunkt x=1x = 1 abstossend.

Für einige ausgewählte Übungen sind weiter unten vorgezeichnete Graphen auf dem Silbertablett serviert.

Exercise 4: Repeller?

Betrachte die Startwerte x=2x = 2 bzw. x=0.5x = 0.5 und zeige allgemein, dass 11 ein Repellor für x2x^2 ist.

Solution

Orbits: {2,4,16,256,}\{2,4,16,256,\dots\} und {0.5,0.25.0.125,}\{0.5,0.25.0.125,\dots\}. Nun betrachten wir einen Punkt in der Nähe von 11 wobei ε>0\varepsilon>0, 1+ε1+\varepsilon. f(1+ε)=(1+ε)2=1+2ε+ε2f(1+\varepsilon) = (1+\varepsilon)^2 = 1+2\varepsilon+\varepsilon^2. Dieser Wert wird pro Iteration zunehmen, also ist 11 abstossend nach oben. Allerdings ist f(1ε)=(1ε)2=12ε+ε2f(1-\varepsilon) = (1-\varepsilon)^2 = 1-2\varepsilon+\varepsilon^2, das gegen 00 tendiert.

Theorem 1: Attraktor-Bedingung

Gilt für einen Fixpunkt xx, dass f(x)>1|f'(x)>1|, dann ist xx ein Repellor. Ist f(x)<1|f'(x)|<1, so handelt es sich um einen Attraktor.

Proof

Übung, falls man zum Beispiel Lipschitz-Stetigkeit und den Mittelwertsatz googlen möchte.

Sei ff Lipschitz-stetig. Der Abstand zwischen xn+1x_{n+1} und xx^* ist xn+1x=f(xn)f(x)Lxnx|x_{n+1}-x^*| = |f(x_n)-f(x^*)|\leq L|x_n-x^*|. Nach Iteration hat man f(f(xn))xLf(xn+1)xL2xnx|f(f(x_n))-x^*|\leq L|f(x_{n+1})-x^*|\leq L^2|x_n-x^*|. Also ist xkxLkx0x|x_k-x^*|\leq L^k|x_0-x^*|. Ist L<1|L|<1 dann gilt Lk0L^k\to0 für kk\to\infty. Ist ff nicht Lipschitz-stetig, so kann man den Mittelwertsatz anwenden, wenn ff differenzierbar und stetig ist. Es ist xk+1x=f(xk)f(x)=f(ζ)(xkx)x_{k+1}-x^* = f(x_k)-f(x^*) = f'(\zeta)\cdot(x_k-x^*) mit ζ(xk,x)\zeta\in(x_k,x^*). Also betragsmässig f(xk)x=f(ζ)xkx|f(x_k)-x^*| = |f'(\zeta)|\cdot|x_k-x^*|, was für f(ζ)<1|f'(\zeta)|<1 den Abstand pro Iterationsschritt verkürzt.

Exercise 5: Repeller, Attraktor, ...

Gib jeweils ein Beispiel einer Funktion mit

a) einem Attraktor bei x=3x = 3.

b) einem Repellor bei x=2x = -2.

c) zwei Fixpunkten, bei x=1x = 1 bzw. x=4x = 4.

d) einem Fixpunkt der Periode 22 in x=7x = 7, für den f(7)=19f(7) = 19 ist.

Solution

a) Ich nehme eine Funktion mit Ableitung an der Stelle 33 betragsmässig kleiner 11. Zum Beispiel f(x)=23x+1f(x) = \frac{2}{3}x+1, denn f(x)=23f'(x) = \frac{2}{3}. f(3)=3f(3) = 3 ist Fixpunkt und f(3)=23|f'(3)| = \frac{2}{3}.

b) Wiederum bietet sich eine lineare Funktion mit Steigung grösser 11 an: g(x)=2x+2g(x) = 2x+2 erfüllt die Anforderungen.

c) Die beiden Werte sind Lösung der quadratischen Gleichung 0=(x1)(x4)=x25x+40 = (x-1)(x-4) = x^2-5x+4 die wir als Fixpunktbedingung umformulieren: x24x+4=xx^2-4x+4 = x. Also erfüllt h(x)=(x2)2h(x) = (x-2)^2 die Bedingung.

d) Aufgrund der vorhergehenden Ausführungen wissen wir, dass für jede lineare Funktion mit Steigung 1-1 jedes xRx\in\mathbb{R} ein 22-Zyklus ist. Für k(x)=x+26k(x) = -x+26 ist auch k(7)=19k(7) = 19 erfüllt.

Exercise 6: Einzugsgebiet

Welche Aussagen kannst du über die Fixpunkte folgender Funktionen machen? Bestimme gegebenenfalls das "Einzugsgebiet", falls es sich um einen Attraktor handelt.

a) cos(x)\cos(x)

b) 2sin(x)2\sin(x)

c) 0.5x2+0.50.5x^2+0.5

Solution

a) cos(x)=x\cos(x) = x ist analytisch nicht lösbar. Man kann beispielsweise gucken, wo die Graphen sich schneiden und Werte approximieren. Eine approximative Lösung mit einem Taylorpolynom wäre auch denkbar.

b) Dasselbe gilt wie in a) gilt hier auch.

c)

0.5x2+0.5=!xx22x+1=0(x1)2=0x=1\begin{align*} 0.5x^2+0.5 &\stackrel{!}{ = } x\\ x^2-2x+1 & = 0\\ (x-1)^2 & = 0\\ x & = 1 \end{align*}

Die Bedingung liefert x=1|x| = 1 und damit ist 11 sicher kein Attraktor.

Exercise 7: Fixpunkte?

Untersuche die Funktion

f(x)=x2+1f(x) = -x^2+1

auf Fixpunkte.

Solution

Aus der Fixpunktbedingung x2+1=!x-x^2+1\stackrel{!}{ = }x folgt 0=x2+x10 = x^2+x-1 deren Lösung die goldenen Schnitte φ\varphi und Φ-\Phi sind. Wegen f(φ)=2φ+1>1f'(\varphi) = 2\varphi+1>1 ist φ\varphi instabil und aus f(Φ)=2Φ+1<1f'(-\Phi) = -2\Phi+1<-1 erkennen wir, dass Φ-\Phi auch Repeller ist.

Die Fixpunkte der Periode 22 geben etwas mehr zu tun.

f(f(x))=(x2+1)2+1=x4+2x2=!x.f(f(x)) = -(-x^2+1)^2+1 = -x^4+2x^2\stackrel{!}{ = }x.

Darin stecken natürlich auch die Fixpunkte der Periode 11 (gegeben durch x2+x1=0x^2+x-1 = 0), die wir mit Hilfe einer Polynomdivision eliminieren, um die Polynomgleichung zu erhalten, welche die reinen Fixpunkte der Periode 2 liefert.

(x4+2x2x)÷(x2+x1)=x2+x(-x^4+2x^2-x)\div(x^2+x-1) = x^2+x

und damit x3=0x_3 = 0 und x4=1x_4 = 1. Es ist (f(f(0)))=403+40=0(f(f(0)))' = -4\cdot0^3+4\cdot0 = 0 und (f(f(1)))=0(f(f(1)))' = 0, also sind beide Attraktoren.

Exercise 8: Standard-Parabel

Untersuche die Standard-Parabel

f(x)=x2f(x) = x^2

auf Fixpunkte.

Solution

Aus der Fixpunktbedingung folgt unmittelbar x2x=!0x^2-x\stackrel{!}{ = }0 und damit x1=0,x2=1x_1 = 0, x_2 = 1. 00 ist Attraktor, 11 Repeller.

Für die Periode 2 Fixpunkte betrachten wir x4x=!0x^4-x\stackrel{!}{ = }0, also x3=1x^3 = 1 und haben immer noch die beiden gleichen Fixpunkte.

Example 6

Wie sieht es nun für eine verschobene Parabel

f(x)=x22f(x) = x^2-2

aus? Wir haben

x2x2=0x^2-x-2 = 0

mit den Lösungen x1=1x_1 = -1 und x2=2x_2 = 2.

Für Fixpunkte der Periode zwei muss

x44x2x+2=0x^4-4x^2-x+2 = 0

gelten. Um nun die echten 22-periodischen Fixpunkte zu erhalten, dividieren wir die Fixpunkte der Periode 11 aus:

(x44x2x+2)÷(x2x2)=x2+x1(x^4-4x^2-x+2)\div(x^2-x-2) = x^2+x-1

und setzen gleich 00. Die Lösung dieser schönen Gleichung kennen wir. Die Fixpunkte 22-ter Periode sind also

x=1±52.x = \frac{-1\pm\sqrt{5}}{2}.

(Parabel kommentiert)