Ziegenproblem (Geometrie)

Ziegenproblem (Geometrie)
Ziegenproblem (Geometrie)

Das Ziegenproblem – auch Die grasende Ziege genannt[1] – ist ein seit dem 18. Jahrhundert bekanntes Problem der Unterhaltungsmathematik. Die erste Veröffentlichung erfolgte 1748 in dem in England einmal jährlich erscheinenden The Ladies Diary: or, the Woman’s Almanack.

Problemstellung

Wie groß muss bei der gezeigten Abbildung r {\displaystyle r} sein, damit die rote Fläche die Hälfte der Kreisfläche ist? Konkrete Motivation: Am Punkt Q {\displaystyle Q} sei eine Ziege (oder ein anderes Tier) angebunden. Wie lang muss die Leine sein, damit das Tier auf genau der Hälfte der Kreisfläche grasen kann?

Lösung mit Berechnung der Linsenfläche

Die von der Ziege erreichbare Fläche hat die Form einer asymmetrischen Linse (siehe Berechnungsskizze), die von zwei Kreisbögen begrenzt wird.[2]

Um den Flächeninhalt A {\displaystyle A} der durch die zwei Kreisbögen begrenzten Fläche zu bestimmen, kann man diese in zwei Kreissegmente zerlegen, wobei die Trennungslinie a {\displaystyle a} in den beiden Schnittpunkten D {\displaystyle D} und E {\displaystyle E} der Kreisbögen endet. Mit R {\displaystyle R} wird der Radius des Kreises, der die Wiese darstellt, und mit r {\displaystyle r} derjenige des Kreises, dessen Mittelpunkt Q {\displaystyle Q} auf dem Kreisrand des anderen liegt, und mit d {\displaystyle d} wird der Abstand zwischen den zwei Kreismittelpunkten P {\displaystyle P} und Q {\displaystyle Q} bezeichnet.

Höhen d1 und d2 der rechtwinkligen Dreiecke MDP bzw. DMQ

Berechnungsskizze Die grasende Ziege

Gegeben: R {\displaystyle R} sowie d = d 1 + d 2 {\displaystyle d=d_{1}+d_{2}}

Nach dem Satz des Pythagoras gilt:

( 1 ) R 2 = d 1 2 + ( a 2 ) 2 ( 1.1 ) ( a 2 ) 2 = R 2 d 1 2 ( 2 ) r 2 = ( d d 1 ) 2 + ( a 2 ) 2 , {\displaystyle {\begin{aligned}(1)\;\;R^{2}={}&d_{1}^{2}+\left({\frac {a}{2}}\right)^{2}\;\;\Rightarrow \;\;(1.1)\;\left({\frac {a}{2}}\right)^{2}=R^{2}-d_{1}^{2}\\(2)\;\;\;r^{2}={}&\left(d-d_{1}\right)^{2}+\left({\frac {a}{2}}\right)^{2},\end{aligned}}}

( a 2 ) 2 {\displaystyle \left({\frac {a}{2}}\right)^{2}} von ( 1.1 ) {\displaystyle \left(1.1\right)} in ( 2 ) {\displaystyle \left(2\right)} eingesetzt ergibt

( 3 ) r 2 = ( d d 1 ) 2 + R 2 d 1 2 , {\displaystyle \;\left(3\right)\;\;r^{2}=\left(d-d_{1}\right)^{2}+R^{2}-d_{1}^{2},}

ausmultipliziert und umgeordnet ergibt

( 4 ) r 2 R 2 = d 2 2 d d 1 , {\displaystyle \;\left(4\right)\;\;r^{2}-R^{2}=d^{2}-2d\cdot d_{1},}

daraus d 1 {\displaystyle d_{1}} ergibt

( 5 ) d 1 = d 2 r 2 + R 2 2 d , {\displaystyle \;\left(5\right)\;\;d_{1}={\frac {d^{2}-r^{2}+R^{2}}{2d}},} [3]

wegen d = d 1 + d 2 {\displaystyle d=d_{1}+d_{2}} wird ( 5 ) {\displaystyle \left(5\right)} entsprechend ergänzt

( 6 ) d = d 2 r 2 + R 2 2 d + d 2 + r 2 R 2 2 d , {\displaystyle \;\left(6\right)\;\;d={\frac {d^{2}-r^{2}+R^{2}}{2d}}+{\frac {d^{2}+r^{2}-R^{2}}{2d}},}

wegen d 2 = d d 1 {\displaystyle d_{2}=d-d_{1}} ergibt sich schließlich

( 7 ) d 2 = d 2 + r 2 R 2 2 d . {\displaystyle \;\left(7\right)\;\;d_{2}={\frac {d^{2}+r^{2}-R^{2}}{2d}}.} [3]

Radius r

Mittels zweimaliger Anwendung der Formel für den Flächeninhalt eines Kreissegments mit d {\displaystyle d} als Abstand des Kreismittelpunktes P {\displaystyle P} bzw. Q {\displaystyle Q} bis zur Kreissehne a {\displaystyle a}

( 8 ) A = R 2 cos 1 ( d R ) d ( R d ) ( R + d ) {\displaystyle \;\left(8\right)\;\;A=R^{2}\cos ^{-1}\left({\frac {d}{R}}\right)-d{\sqrt {\left(R-d\right)\left(R+d\right)}}} [4]

und den darin eingesetzten Termen d 1 {\displaystyle d_{1}} von ( 5 ) {\displaystyle (5)} und d 2 {\displaystyle d_{2}} von ( 7 ) , {\displaystyle (7),} erhält man nach entsprechender Umformulierung die Formel für den Flächeninhalt der asymmetrischen Linse:[3]

( 9 ) A = A 1 ( R , d 1 ) + A 2 ( r , d 2 ) = R 2 cos 1 ( d 2 r 2 + R 2 2 d R ) + r 2 cos 1 ( d 2 + r 2 R 2 2 d r ) 1 2 ( d r + R ) ( d + r + R ) ( d + r R ) ( d + r + R ) . {\displaystyle {\begin{aligned}(9)\;\;A={}&A_{1}\left(R,d_{1}\right)+A_{2}\left(r,d_{2}\right)\\={}&R^{2}\cos ^{-1}\left({\frac {d^{2}-r^{2}+R^{2}}{2dR}}\right)+r^{2}\cos ^{-1}\left({\frac {d^{2}+r^{2}-R^{2}}{2dr}}\right)\;-\\&{}{\frac {1}{2}}{\sqrt {(d-r+R)(d+r+R)(d+r-R)(-d+r+R)}}\,.\end{aligned}}}

Für R = d = 1 {\displaystyle R=d=1} und halber Kreisfläche vereinfacht sich dies zu

( 10 ) 1 2 π = cos 1 ( 1 1 2 r 2 ) + r 2 cos 1 ( 1 2 r ) 1 2 r 4 r 2 . {\displaystyle \;\left(10\right)\;{\frac {1}{2}}\pi =\cos ^{-1}\left(1-{\frac {1}{2}}r^{2}\right)+r^{2}\cos ^{-1}\left({\frac {1}{2}}r\right)-{\frac {1}{2}}r{\sqrt {4-r^{2}}}.} [2]

Diese Gleichung kann nur numerisch gelöst werden und ergibt r = 1,158 7284 {\displaystyle r=1{,}1587284\ldots } (Folge A133731 in OEIS).

Lösung mit Integration

Aus der Integration über die rechte Hälfte der Linsenfläche mit

1 4 π = 0 r 2 r 4 4 ( r 2 t 2 + 1 t 2 1 ) d t {\displaystyle {\frac {1}{4}}\pi =\int _{0}^{\sqrt {r^{2}-{\frac {r^{4}}{4}}}}\left({\sqrt {r^{2}-t^{2}}}+{\sqrt {1-t^{2}}}-1\right)\,\mathrm {d} t}

ergibt sich die ebenfalls transzendente Gleichung

r = π π r 4 r 2 + ( 4 r 2 r ) sin 1 ( 1 2 r ) {\displaystyle r={\frac {\pi }{\pi r-{\sqrt {4-r^{2}}}+\left({\frac {4}{r}}-2r\right)\sin ^{-1}\left({\frac {1}{2}}r\right)}}}

mit der gleichen Lösung.

Geometrische Näherungslösung

Die grasende Ziege,
mit Approximation der halben Wiesenfläche (grün).
{\displaystyle } | Q C ¯ | = r {\displaystyle |{\overline {QC}}|=r} die Länge der Leine

Zwei sich schneidende Kreise und deren Schnittpunkt liefern den gesuchten Radius, der die kreisförmige Wiesenfläche nahezu halbiert.

Konstruktion

Es beginnt mit dem Einheitskreis um Punkt P {\displaystyle P} und dem Einzeichnen von zwei zueinander senkrecht stehenden Radien; dabei ergeben sich die Schnittpunkte A {\displaystyle A} und Q . {\displaystyle Q.} Es folgt der Kreisbogen k b {\displaystyle kb} um A {\displaystyle A} mit Radius | A B ¯ | = 1 3 | A P ¯ | . {\displaystyle |{\overline {AB}}|={\tfrac {1}{3}}\cdot |{\overline {AP}}|.} Er schneidet den Kreis in C {\displaystyle C} und bringt damit den gesuchten Radius r {\displaystyle r} als Länge der Strecke | Q C ¯ | . {\displaystyle |{\overline {QC}}|.} Der abschließende Kreisbogen um Punkt Q {\displaystyle Q} mit Radius r {\displaystyle r} ab C {\displaystyle C} schneidet den Kreis in D {\displaystyle D} und liefert nahezu eine Halbierung der Kreisfläche.

Nachrechnung

Berechnungsskizze

Aus der nebenstehenden Berechnungsskizze bzw. aus der obigen Konstruktionsbeschreibung ist zu entnehmen:

  • k 1 {\displaystyle k_{1}} ist der Einheitskreis mit der Gleichung
k 1 : x 2 + y 2 = 1 {\displaystyle k_{1}:x^{2}+y^{2}=1}
  • k b {\displaystyle kb} ist ein Teil des Kreises k 2 {\displaystyle k_{2}} mit dem Radius 1 3 | A P ¯ | {\displaystyle {\frac {1}{3}}\cdot |{\overline {AP}}|} und der Gleichung
k 2 : ( x 1 ) 2 + y 2 = 0 , 1 ¯ = 1 9 {\displaystyle k_{2}:\left(x-1\right)^{2}+y^{2}=0,{\bar {1}}={\frac {1}{9}}}
  • Punkt C {\displaystyle C} ist der Schnittpunkt des Kreises k 1 {\displaystyle k_{1}} mit dem Kreis k 2 . {\displaystyle k_{2}.} Die Abstände des Punktes C {\displaystyle C} in einem kartesischen Koordinatensystem sind: x = | F C ¯ | {\displaystyle x=|{\overline {FC}}|} und y = | P F ¯ | . {\displaystyle y=|{\overline {PF}}|.}

Die Länge der Strecke | F C ¯ | = x {\displaystyle |{\overline {FC}}|=x} erhält man durch Subtraktion der beiden Kreisgleichungen:[5]

k 2 : ( x 1 ) 2 + y 2 = 1 9 x 2 2 x + 1 + y 2 = 1 9 k 1 : x 2 + y 2 = 1 | k 2 k 1 1 2 x = 8 9 x = 17 18 . {\displaystyle {\begin{aligned}\;k_{2}:\;\;\left(x-1\right)^{2}\;+\;\;y^{2}&={\frac {1}{9}}\\x^{2}-2x+1+y^{2}&={\frac {1}{9}}\\k_{1}:\;\;\;\;\;\;\;\;x^{2}\;\;\;+\;\;\;y^{2}&=1|\;\;k_{2}-k_{1}\\1-2x&=-{\frac {8}{9}}\Rightarrow \\x&={\frac {17}{18}}.\end{aligned}}}

Die Länge der Strecke | P F ¯ | {\displaystyle |{\overline {PF}}|} erhält man durch Einsetzen des x-Wertes in die Kreisgleichung des Kreises k 1 : {\displaystyle k_{1}:}

k 1 : ( 17 18 ) 2 + y 2 = 1 y 2 = 1 289 324 y 2 = 35 324 y = 35 18 {\displaystyle {\begin{aligned}\;k_{1}:\left({\frac {17}{18}}\right)^{2}+y^{2}&=1\\y^{2}&=1-{\frac {289}{324}}\\y^{2}&={\frac {35}{324}}\\y&={\frac {\sqrt {35}}{18}}\end{aligned}}}

Somit ist die Länge der Strecke | F Q ¯ | = 1 35 18 . {\displaystyle |{\overline {FQ}}|=1-{\frac {\sqrt {35}}{18}}.}

Für die Hypotenuse | Q C ¯ | = r {\displaystyle |{\overline {QC}}|=r} des rechtwinkligen Dreiecks Q C F {\displaystyle QCF} gilt nach dem Satz des Pythagoras:

r = ( | F Q ¯ | ) 2 + ( | F C ¯ | ) 2 = ( 1 35 18 ) 2 + ( 17 18 ) 2 = 18 35 3 r = 1,158 731117161276 [ L E ] {\displaystyle {\begin{aligned}\;r&={\sqrt {\left(|{\overline {FQ}}|\right)^{2}+\left(|{\overline {FC}}|\right)^{2}}}\\&={\sqrt {\left(1-{\frac {\sqrt {35}}{18}}\right)^{2}+\left({\frac {17}{18}}\right)^{2}}}\Rightarrow \\&={\frac {\sqrt {18-{\sqrt {35}}}}{3}}\Rightarrow \\r&=1{,}158731117161276\ldots [LE]\end{aligned}}}

Absoluter Fehler der konstruierten Länge r {\displaystyle r} der Leine; darin entspricht Radius r n {\displaystyle r_{n}} dem numerisch gelösten r {\displaystyle r} (s. oben):

F r = r r n = 1,158 731117 1,158 728473 = 0,000 002644 {\displaystyle F_{r}=r-r_{n}=1{,}158731117-1{,}158728473=0{,}000002644\ldots } [LE]

Für den relativen Fehler des konstruierten Radius r {\displaystyle r} gilt:

f r = ( r r n 1 ) 100 % , {\displaystyle f_{r}=\left({\frac {r}{r_{n}}}-1\right)\cdot 100\,\%,}

mit den eingesetzten Werten ergibt sich

f r = ( 1,158 731117 1,158 728473 1 ) 100 % 0,000 228 % . {\displaystyle f_{r}=\left({\frac {1{,}158731117}{1{,}158728473}}-1\right)\cdot 100\,\%\approx 0{,}000228\%.}

Den Radius r {\displaystyle r} eingesetzt in die vereinfachte Formel der Linsenfläche für den Einheitskreis (mit R = d = 1 {\displaystyle R=d=1} ), oben in Lösung mit Berechnung der Linsenfläche beschrieben, ergibt näherungsweise die konstruierte halbe, im Bild grüne, Wiesenfläche

A k = 1,158 731117 2 cos 1 ( 1 2 1,158 731117 ) + cos 1 ( 1 1 2 1,158 731117 2 ) 1 2 1,158 731117 4 1,158 731117 2 = 1,570 802165 {\displaystyle {A_{k}=1{,}158731117^{2}\cdot \cos ^{-1}\left({\frac {1}{2}}\cdot 1{,}158731117\right)+\cos ^{-1}\left(1-{\frac {1}{2}}\cdot 1{,}158731117^{2}\right)-{\frac {1}{2}}\cdot 1{,}158731117\cdot {\sqrt {4-1{,}158731117^{2}}}=1{,}570802165\ldots }} [FE].

Flächeninhalt der halben Wiese (halber Einheitskreis)

A W = π 2 = 1,570 796326 {\displaystyle A_{W}={\frac {\pi }{2}}=1{,}570796326\ldots } [FE]

Absoluter Fehler der approximierten halben Wiesenfläche

F A k = A k A W = 1,570 802165 1,570 796326 = 0,000 005829 {\displaystyle F_{A_{k}}=A_{k}-A_{W}=1{,}570802165-1{,}570796326=0{,}000005829\ldots } [FE]

Relativer Fehler der approximierten halben Wiesenfläche (Formel siehe oben bei f r {\displaystyle f_{r}} )

f A k = ( 1,570 802165 1,570 796326 1 ) 100 % 0,000 371 % {\displaystyle f_{A_{k}}=\left({\frac {1{,}570802165}{1{,}570796326}}-1\right)\cdot 100\,\%\approx 0{,}000371\,\%}

Verdeutlichung der Approximation

Hätte z. B. die kreisförmige Wiese einen Radius gleich 100 m , {\displaystyle 100\;\mathrm {m} ,} dann wäre die Leine um ca. 0 , 3 m m {\displaystyle 0{,}3\;\mathrm {mm} } zu lang und die Ziege könnte – angebunden am Punkt Q {\displaystyle Q} an eine Leine mit der Länge 115 , 8731 m {\displaystyle 115,8731\;\mathrm {m} } – außer der für sie vorgesehenen Hälfte der Wiesenfläche (rund 15.707 , 963 m 2 {\displaystyle 15.707,963\;\mathrm {m} ^{2}} ), noch zusätzlich 5 , 8 d m 2 {\displaystyle 5{,}8\;\mathrm {dm} ^{2}} abgrasen, das wären etwa 40 m m 2 {\displaystyle 40\;\mathrm {mm} ^{2}} weniger als ein DIN-A4-Blatt.

Geschlossene Lösung

Mit Methoden der komplexen Geometrie fand Ingo Ullisch im Jahr 2020 folgende geschlossene Lösung[6][7][8]

r = 2 cos ( 1 2 | z 3 π / 4 | = π / 4 z / ( sin z z cos z π / 2 ) d z | z 3 π / 4 | = π / 4 1 / ( sin z z cos z π / 2 ) d z ) {\displaystyle r=2\cos \left({\frac {1}{2}}{\frac {\oint _{|z-3\pi /4|=\pi /4}z/(\sin z-z\cos z-\pi /2)\,dz}{\oint _{|z-3\pi /4|=\pi /4}1/(\sin z-z\cos z-\pi /2)\,dz}}\right)}

Erweiterungen

Die Ziege im Weltall

Der dreidimensionale Fall mit Einheitskugel oben und Ziegenkugel unten

Im dreidimensionalen Fall befindet sich Punkt Q {\displaystyle Q} auf der Oberfläche einer Einheitskugel, und die Fragestellung ist, wie groß der Radius r {\displaystyle r} der zweiten Kugel sein muss, damit der Schnittkörper genau die Hälfte des Volumens der Einheitskugel hat.

Der vom Tier erreichbare Teil der Einheitskugel hat die Form einer dreidimensionalen Linse mit unterschiedlich gewölbten Seiten, die von den beiden Kugelkalotten begrenzt wird.

Das Volumen V {\displaystyle V} einer Linse bei zwei Kugeln mit Radien R , r {\displaystyle R,r} und Mittelpunktabstand d {\displaystyle d} ist:

V = π ( R + r d ) 2 ( d 2 + 2 d r 3 r 2 + 2 d R + 6 r R 3 R 2 ) 12 d {\displaystyle V={\frac {\pi (R+r-d)^{2}\left(d^{2}+2dr-3r^{2}+2dR+6rR-3R^{2}\right)}{12d}}}

was sich bei R = d = 1 {\displaystyle R=d=1} und halbem Kugelvolumen vereinfacht zu

1 2 4 3 π = 1 4 π r 4 + 2 3 π r 3 {\displaystyle {\frac {1}{2}}\cdot {\frac {4}{3}}\pi =-{\frac {1}{4}}\pi r^{4}+{\frac {2}{3}}\pi r^{3}}

woraus als Lösung folgt r = 1,228 5 {\displaystyle r=1{,}2285\ldots }

Es kann gezeigt werden, dass sich r {\displaystyle r} bei weiter zunehmender Dimensionalität dem Wert 2 {\displaystyle {\sqrt {2}}} annähert.

Die Ziege am Silo

Das Ziegenproblem mit Silo für Seillängen von 1 4 π r {\displaystyle {\tfrac {1}{4}}\cdot \pi r} , 1 2 π r {\displaystyle {\tfrac {1}{2}}\cdot \pi r} , 3 4 π r {\displaystyle {\tfrac {3}{4}}\cdot \pi r} , 1 π r {\displaystyle 1\cdot \pi r} und 1 r {\displaystyle 1\cdot r} (grün)
Die Ziege außerhalb des Kreises

Im zweidimensionalen Fall kann auch die Frage nach der Größe der erreichbaren Fläche außerhalb des roten Kreises gestellt werden. Das entspricht der Situation, dass das Tier an einem Silo festgebunden ist.

In diesem Fall besteht die Fläche aus einem Halbkreis (hellblau) mit Radius r {\displaystyle r} und zwei Flächen, die durch den roten Kreis und die Kreisevolvente begrenzt sind (dunkelblau). Aus der Sektorformel von Leibniz folgt der Inhalt einer der dunkelblauen Flächen. Die gesamte erreichbare Fläche (hell- und dunkelblau) beträgt dann

A = 1 2 π r 2 + 1 3 r 3 {\displaystyle A={\frac {1}{2}}\pi r^{2}+{\frac {1}{3}}r^{3}}

unter der Bedingung, dass r π {\displaystyle r\leq \pi } (andernfalls überschneiden sich die beiden dunkelblauen Flächen auf der Rückseite des Silos).

Siehe auch

Literatur

Einzelnachweise

  1. Heinrich Hemme: Die Hölle der Zahlen. 92 mathematische Rätsel mit ausführlichen Lösungen. Vandenhoeck & Ruprecht, Göttingen 2007, ISBN 978-3-525-40841-4, S. 32 und 102 f. 
  2. a b Eric Weisstein: Goat Problem. WolframMathworld, abgerufen am 17. November 2019. 
  3. a b c Eric Weisstein: Circle-Circle Intersection. WolframMathworld, abgerufen am 17. November 2019. 
  4. Eric Weisstein: Circular Segment. WolframMathworld, abgerufen am 17. November 2019. 
  5. Schnitt und Schnittwinkel zweier Kreise. (PDF) Nichtlineare analytische Geometrie. In: Dimensionen – Mathematik 7. Verlag E. Dorner, Wien, 2016, S. 1, abgerufen am 19. November 2020. 
  6. Ingo Ullisch: A Closed-Form Solution to the Geometric Goat Problem. In: The Mathematical Intelligencer. 18. Februar 2020, ISSN 0343-6993, doi:10.1007/s00283-020-09966-0 (englisch). 
  7. Steve Nadis: A After Centuries, a Seemingly Simple Math Problem Gets an Exact Solution. Quanta Magazin, 9. Dezember 2020, abgerufen am 23. Februar 2022. 
  8. Ingo Ullisch: Correction to: A Closed-Form Solution to the Geometric Goat Problem. Abgerufen am 3. März 2024.