Wie zeichnet man einen Stern (und mehr) in Polarkoordinaten

Die Frage nach der Formel fĂŒr ein Polygon in Polarkoordinaten stellt sich regelmĂ€ĂŸig nach thematischen Ressourcen - und bleibt ebenso regelmĂ€ĂŸig ohne klare Antwort. Im besten Fall stĂ¶ĂŸt man durch die Funktion des Restes der Division auf eine Lösung - die aus mathematischer Sicht nicht "sauber" ist, da keine analytischen Transformationen fĂŒr die Funktion möglich sind. Anscheinend sind echte Mathematiker zu beschĂ€ftigt, Probleme des Jahrtausends zu lösen und nach einem einfachen Beweis fĂŒr Fermats Theorem zu suchen, um auf solche banalen Probleme zu achten. GlĂŒcklicherweise ist in dieser Angelegenheit die Vorstellungskraft wichtiger als das Wissen, und um dieses Problem zu lösen, mĂŒssen Sie kein Professor fĂŒr topologische Wissenschaften sein - Wissen auf Schulebene reicht aus.



Die Formel fĂŒr ein gleichseitiges Polygon in Polarkoordinaten sieht sehr einfach aus

ρ=cos⁥(2sin−1⁥(k)+πm2n)cos⁥(2sin−1⁥(kcos⁥(nϕ))+πm2n)



und hat folgende Parameter:



ϕ ist der Winkel;

n ist die Anzahl der konvexen Eckpunkte;

m - bestimmt, wie viele Eckpunkte die Seiten auf einer geraden Linie liegen. Negative Werte sind ebenfalls zulÀssig - das Vorzeichen bestimmt, in welche Richtung sich der Stern biegt;

k - Steifheit - beik=0 Wir erhalten einen Kreis unabhÀngig von anderen Parametern, zk=1 - Polygon mit geraden Linien bei Zwischenwerten von0 bis1 - Zwischenfiguren zwischen einem Kreis und einem Polygon.



Mit dieser Formel können Sie einen Stern auf zwei Arten zeichnen:



1)n=5,m=3





2) n=5/4,m=0 . In diesem Fall mĂŒssen Sie zwei Runden anstelle von einer machen:





Parameter m wirkt sich wie folgt auf das Polygon aus (hier Àndert es sich von -1 auf 5):





Parameter k in der Animation:





Komplexe Form und Modifikationen



Sie können die ursprĂŒngliche Formel in einer komplexen Form umschreiben, und trotz des Vorhandenseins imaginĂ€rer Einheiten bleibt der Wert des Radius weiterhin gĂŒltig:



ρ=41/n(1−k2+ik)−1/n(1+(1−k2+ik)2/neiπmn)(1−k2cos2⁥(nϕ)+ikcos⁥(nϕ))1/n41/n+eiπmn(21−k2cos2⁥(nϕ)+2ikcos⁥(nϕ))2/n



Auf den ersten Blick mag dies sinnlos erscheinen, da die Formel etwas umstĂ€ndlicher geworden ist - aber keine Schlussfolgerungen ziehen. Erstens enthĂ€lt es keinen Arkussinus, der die mathematische Bedeutung der Formel vollstĂ€ndig Ă€ndert und es Ihnen ermöglicht, die Konstruktion eines sternförmigen Polygons anders zu betrachten. Zweitens können damit auch kompakte Formeln fĂŒr SonderfĂ€lle erhalten werden . Drittens (und am interessantesten) kann es kreativ modifiziert und andere, unerwartete Formen erhalten werden. Damit das Auftreten einer möglichen imaginĂ€ren Komponente im Radius keine Mehrdeutigkeit bei der Berechnung verursacht, kann sie durch Multiplikation mit sofort auf kartesische Koordinaten reduziert werdenit(int)1/n1+(int)2/n . Hier sind Beispiele fĂŒr einige der Modifikationen:eiϕ



(−1)2/3eiϕ(icos⁥(3ϕ)+sin2⁥(3ϕ))1/3(−1)2/3+22/3(icos⁥(3ϕ)+sin2⁥(3ϕ))2/3





eiϕicos⁥(2ϕ)+sin2⁥(2ϕ)ei/2+2(icos⁥(2ϕ)+sin2⁥(2ϕ))3/2





e14i(4ϕ+π)2(−1)1/4cos⁥(2ϕ)+sin2⁥(2ϕ)−1−i



Wie Sie wahrscheinlich bemerkt haben, ist die Rotation des Vektors nicht mehr gleichmĂ€ĂŸig - und zwar genau aufgrund des Auftretens der imaginĂ€ren Komponente im Radius.



Quads und so



Unsere Formel hat einen wunderbaren Sonderfall - ein Quadrat, fĂŒr das die Formel geschrieben werden kann

ρ=22+2+2cos⁥(4ϕ)



oder

ρ=21+1−sin2⁥(2ϕ)



(WÀhlen Sie aus, welche Ihnen am besten gefÀllt).



In einem etwas weiter entwickelten Fall können Sie Zwischenfiguren zwischen einem Kreis und einem Quadrat durch einen Punkt definieren(k,k) auf der OberflÀche

ρ=21+1−(2k2−1)sin2⁥(2ϕ)k4





Sie können diesen Formen auch Variationen hinzufĂŒgen, wĂ€hrend Sie den Zustand beibehalten, in dem sie den Punkt durchlaufen (k,k) - den Parameter selbst direkt modulieren kabhĂ€ngig vom Winkel, so dass beim Durchlaufen der Diagonale sein Multiplikator gleich eins ist. Zum Beispiel stattdessen ersetzenk Funktion k1−zcos2⁥(2ϕ)erhalten wir einen zusĂ€tzlichen Parameter zmit denen zusĂ€tzliche Biegungen eingestellt werden können. Insbesondere fĂŒrz=1/4 Sie erhalten Folgendes:







In einem noch erweiterten Fall können Sie nicht nur ein Quadrat, sondern auch ein Rechteck und immer noch in Polarkoordinaten definieren:

ρ=4a2b2((b2−a2)cos⁥(2ϕ)+a2+b2)(1−4a2b2ksin2⁥(2ϕ)((b2−a2)cos⁥(2ϕ)+a2+b2)2+1)



Und berechnen Sie sogar seine FlĂ€che (ĂŒber elliptische Integrale):

S=4ab(k−1)K(k)+E(k)k

Hinweis
k (0 1) , πab 4ab.


Auf diese Weise können Profile mit einem Übergang von einem Kreis zu einem Rechteck mit einer kontrollierten SchnittflĂ€che erstellt werden. Hier ist die FlĂ€che konstant:





Und hier erweitert sich der Bereich exponentiell:





Gehen Sie zu den kartesischen Koordinaten



Jede Formel in Polarkoordinaten kann durch eine Gleichung in kartesischen Koordinaten und auf mindestens zwei Arten ausgedrĂŒckt werden - je nachdem, welche Form des Gradienten am Rand der Figur sich Ă€ndert. Dazu reicht es aus, den Winkel durch den Arkustangens der Koordinaten zu berechnen und die Formel durch Subtrahieren durch den Radiusvektor auf eine Konstante zu bringen

0=x2+y2−cos⁥(2sin−1⁥(k)+πm2n)cos⁥(2sin−1⁥(kcos⁥(ntan−1⁥(x,y)))+πm2n)



oder Teilung

1=x2+y2cos⁥(2sin−1⁥(kcos⁥(ntan−1⁥(x,y)))+πm2n)cos⁥(2sin−1⁥(k)+πm2n)



Die zweite Option ist vorzuziehen, da sie gerade VerlÀufe entlang der Seiten des Polygons ergibt.



Hinweis
, (0,0) - — , , ( −cos⁥(2sin−1⁥(k)+πm2n)sec⁥(2sin−1⁥(kcos⁥(πn2))+πm2n) ).



cos⁡(ntan−1⁡(x,y)) (x+iy)n+(x−iy)n2(x2+y2)n/2, A034839.


Der Wert der Formel von der rechten Seite der Gleichung (im 2. Fall) Ă€ndert sich von 0 Vor 1 wenn Punkt (x,y) kommt in die Figur und von 1ad infinitum - wenn nicht. Durch Auswahl verschiedener Funktionen zum Konvertieren in Helligkeit erhalten Sie verschiedene Rasteroptionen. FĂŒr den Exponenten (e−x−1 fĂŒr die erste und e−x fĂŒr die zweite Option) bekommen wir

oder wenn mit SĂ€ttigung



Es kann ein klassischer Tiefpassfilter verwendet werden 11+xp, indem p- die Reihenfolge des Filters, die den DĂ€mpfungsgrad bestimmt.



FĂŒr die erste Option:



Und zum zweiten:



Sie können auch eine stĂŒckweise kontinuierliche Funktion verwenden, indem Sie die Interpolationsgrenzen explizit festlegen.



ZusĂ€tzlich zur Rasterung als solche können Sie auch Verformungen festlegen - z. B. ein Schachbrett in einen Kreis drĂŒcken:





Oder ziehen Sie es sogar ĂŒber eine Kugel:

Formel

x=u21+|u2−v2u2+v2|



y=v1+2|u2−v2u2+v2|



z=1−12u2(1+|u2−v2u2+v2|)−12v2(1+|u2−v2u2+v2|)





Anhang: Wie die Formel abgeleitet wurde



Der klassische Stil des GeschichtenerzĂ€hlens in mathematischen Texten besteht aus dem Wechsel von Deckspelzen / Theoremen und ihren Beweisen - als ob die nachweisbaren Aussagen in den Köpfen der Autoren durch eine Offenbarung von oben erschienen wĂ€ren. Und obwohl dies eine gewisse Wahrheit ist, gehen dem Erscheinen von Formeln hĂ€ufiger einige Forschungsarbeiten voraus, deren Beschreibung ein besseres VerstĂ€ndnis ihrer Bedeutung vermitteln kann als formale Beweise; und die Treue der Aussagen kann wiederum durch die Treue der Schritte verfolgt werden, die zu ihnen gefĂŒhrt haben.



Auch hier - wenn der Artikel mit einer Formel in komplexer Form beginnen wĂŒrde, wĂ€re sein Erscheinungsbild nicht offensichtlich und nicht intuitiv, und die deklarierten Eigenschaften wĂŒrden zusĂ€tzliche Beweise erfordern. In der trigonometrischen Form der Aufzeichnung ist es jedoch durchaus möglich, die Geschichte ihres Auftretens zu verfolgen.



1) Wir beginnen mit dem einfachsten Fall - der Aufgabe, eine gerade Linie in Polarkoordinaten zu zeichnen. Dazu mĂŒssen Sie die Gleichung lösenrcos⁥(ϕ)=1deren Lösung ist offensichtlich r→sec⁥(ϕ)...





2) Ferner muss das Sekantenargument "geloopt" werden, um Knicke in einer geraden Linie bereitzustellen. In diesem Stadium verwenden andere Lösungen einen "Dirty Hack" in Form eines Divisionsrestes. Es wird auch die sequentielle Aufnahme der direkten und inversen Sinusfunktion verwendet -sin−1⁥(sin⁥(ϕ))





Mit diesem Ansatz können Sie mathematische Standardoperationen fĂŒr die resultierende Formel ausfĂŒhren.
z.B
:

∂sin−1⁥(sin⁥(ϕ))∂ϕ=cos⁥(ϕ)1−sin2⁥(ϕ)







Dank der gleichen Notation können Sie die Quadratfunktion in Polarkoordinaten mithilfe der komplexen Darstellung trigonometrischer Funktionen zu einem Àsthetischeren Erscheinungsbild vereinfachen. In Wolfram Mathematica kann dies mit den Funktionen TrigToExp und ExpToTrig erfolgen:



Der Code
Sec[1/2 ArcSin[k Sin[2 \[Phi]]]]^2//TrigToExp//ExpToTrig//Sqrt[#]&//FullSimplify

↓

22+21−k2sin2⁥(2ϕ)



Dank derselben Aufnahme können Sie mit einem zusĂ€tzlichen Multiplikator glatte Zwischenfiguren zwischen einem Kreis und einem Quadrat erhalten kaufgrund dessen das arcsine-Argument eins verfehlt - sin−1⁥(ksin⁥(ϕ))::





Und damit die Funktion einen bestimmten Punkt schneidet, mĂŒssen Sie nur eine Gleichung erstellen und neu berechnen k::

21+1−kâ€Čsin2⁥(π2)=k



Der Code
Solve[(Sqrt[2/(1+Sqrt[1-k Sin[2 \[Phi]]^2])] /. \[Phi]->Pi/4)==x, k] /. x->k

↓

kâ€Č→4(k2−1)k4





3) Parameter n und mwurden nur auf kreative Weise hinzugefĂŒgt und ihre Auswirkungen wurden tatsĂ€chlich experimentell untersucht.



4) Das Rechteck ist leicht zu erhalten, indem Sie zur parametrischen Ansicht gehen und die Achsen "strecken"

x=acos⁡(t)21+1−sin2⁡(2t)



y=bsin⁡(t)21+1−sin2⁡(2t)



Aber danach t wird jetzt nicht mehr Winkel bedeuten tIst nur ein Parameter, der einen Vektor durch seine Projektionen auf die Koordinatenachsen beschreibt. Um zu den Polarkoordinaten zurĂŒckzukehren, mĂŒssen Sie die LĂ€nge des Vektors (durch die Wurzel der Quadratsumme) und den Winkel (durch den Arkustangens des VerhĂ€ltnisses) ermitteln und diesen Winkel durchdrĂŒckenϕ und ersetzen Sie stattdessen den resultierenden Ausdruck t...



Der Code
With[{r = Sqrt[2/(1 + Sqrt[

1 - Sin[2 t]^2])]}, {Sqrt[(a r Cos[t])^2 + (b r Sin[t])^2],

ArcTan[(b r Sin[t])/(a r Cos[t])]}] // Simplify


↓

{2a2cos2⁡(t)+b2sin2⁡(t)cos2⁡(2t)+1,tan−1⁡(btan⁡(t)a)}





Solve[ArcTan[(b Tan[t])/a]==\[Phi], t]

↓

t→tan−1⁥(atan⁥(ϕ)b)





Sqrt[2] Sqrt[(a^2 Cos[t]^2 + b^2 Sin[t]^2)/(1 + Sqrt[Cos[2 t]^2])]

/. t -> ArcTan[(a Tan[\[Phi]])/b] // Simplify


↓



2a2b2sec2⁥(ϕ)(a2tan2⁥(ϕ)+b2)(cos2⁥(2tan−1⁥(atan⁥(ϕ)b))+1)



Die Vereinfachung einer solchen Formel ist bereits schwieriger, und dies erfordert mehrere Schritte:



  1. Gehen Sie zu kartesischen Koordinaten, indem Sie sie ersetzen ϕ→tan−1⁥(x,y);;
  2. gehe zur Exponentialform;
  3. vereinfachen;
  4. Machen Sie einen umgekehrten Austausch x→cos⁥(ϕ) und y→sin⁥(ϕ);;
  5. gehe zurĂŒck zur Exponentialform;
  6. vereinfachen.


Als Ergebnis erhalten wir die folgende Formel:



Der Code
Sqrt[2] Sqrt[(a^2 b^2 Sec[\[Phi]]^2) /

((1 + Sqrt[Cos[2 ArcTan[(a Tan[\[Phi]])/b]]^2])

(b^2 + a^2 Tan[\[Phi]]^2))] /. \[Phi] -> ArcTan[x, y]

// TrigToExp // Simplify

// # /. {x -> Cos[\[Phi]], y -> Sin[\[Phi]]} &

// TrigToExp // Simplify // FullSimplify


↓

2a2b2((b2−a2)cos⁥(2ϕ)+a2+b2)(1+((a2+b2)cos⁥(2ϕ)−a2+b2)2((b2−a2)cos⁥(2ϕ)+a2+b2)2)









Fazit



Wie Sie sehen können, können Sie selbst in einer so einfachen und alltĂ€glichen Sache wie einem Polygon etwas Neues finden und sich etwas einfallen lassen. Und die Geschichte endet nicht dort - die FlĂ€chenformel fĂŒr den allgemeinen Fall blieb unbekannt, die Formel fĂŒr ein beliebiges und nicht nur regulĂ€res Polygon blieb unbekannt, und die Erweiterung in Potenz- und trigonometrische Reihen wurde ohne BerĂŒcksichtigung gelassen. Wahrscheinlich gibt es auch eine Ă€hnliche Formel fĂŒr den dreidimensionalen Fall.



Wenn Ihnen also gesagt wird, dass alles in der Mathematik bereits erfunden wurde und es nur Probleme gibt, die fĂŒr einen normalen Menschen unverstĂ€ndlich sind, glauben Sie es nicht. Es gibt viele rein praktische Probleme, die echten Mathematikern nicht bekannt sind, oder sie sind nicht an ihrer Lösung interessiert, weil sie keinen ausreichenden Hype um sich haben oder weil sie bereits eine ungefĂ€hre Vorstellung davon haben, wie sie erreicht werden können. Haben Sie keine Angst, Probleme anzugehen, deren Lösungen in Wikipedia nicht verfĂŒgbar sind, haben Sie keine Angst, ihre Lösungen zu veröffentlichen, und haben Sie keine Angst, die Kommentare unter den Artikeln ĂŒber die Nutzlosigkeit von allem zu lesen.



PS das Originaldokument fĂŒr Mathematica herunterladen hier .



All Articles