極座標で星(およびそれ以上)を描く方法

極座標のポリゴンの公式の問題は、主題のリソースで定期的に発生します-そして、明確な答えがないまま定期的に残ります。最良の場合、分割の残りの部分の関数を介して解決策に出くわします。これは、関数に対して分析変換を実行できないため、数学的な観点からは「クリーン」ではありません。どうやら、本当の数学者は千年紀の問題を解決するのに忙しくて、そのような平凡な問題に注意を払うにはフェルマットの定理の簡単な証拠を探しています。幸いなことに、この問題では、想像力が知識よりも重要であり、この問題を解決するには、トポロジー科学の教授である必要はありません。学校レベルの知識で十分です。



極座標での等辺ポリゴンの式は非常に単純に見えます

ρ=cos(2sin1(k)+πm2n)cos(2sin1(kcos(nϕ))+πm2n)



次のパラメータがあります。



ϕ -角度;

n -凸状の頂点の数。

m-側面が1つの直線上にある頂点の数を決定します。負の値も許容されます-記号は星が曲がる方向を決定します;

k -剛性-で k=0 他のパラメータに関係なく、円を取得します。 k=1 -からの中間値での直線を持つポリゴン 01-円とポリゴンの中間の図形。



この式を使用すると、2つの方法で星を描くことができます:



1)n=5,m=3





2) n=5/4,m=0..。この場合、1ターンではなく2ターンする必要があります。





パラメータ m ポリゴンに次のように影響します(ここでは-1から5に変更されます)。





パラメータ k アニメーションで:





複雑なフォームと変更



元の式を複雑な形式で書き直すことができます。仮想単位が含まれていても、半径の値は引き続き有効です。



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



公式がもう少し面倒になっているので、一見、これは無意味に思えるかもしれませんが、結論にジャンプしないでください。まず、アークサインが含まれていません。これにより、式の数学的な意味が完全に変わり、星型ポリゴンの構築を別の方法で見ることができます。第二に、それはまた、例えば、特別な場合のためのコンパクトな式を得るために使用することができますit(int)1/n1+(int)2/n..。第三に(そして最も興味深いことに)、それは創造的に修正され、他の予期しない形を得ることができます。半径内の可能な虚数成分の出現が計算のあいまいさを引き起こさないようにするために、を掛けることによってそれをすぐにカルテシアン座標に減らすことができます。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ϕ)1i



お気づきかもしれませんが、ベクトルの回転は均一ではなくなりました。これは、半径に虚数成分が出現したためです。



クワッドなど



私たちの式には素晴らしい特別なケースがあります-正方形、その式は次のように書くことができます

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



または

ρ=21+1sin2(2ϕ)



(あなたが一番好きなものを選んでください)。



もう少し発展したケースでは、点を通る円と正方形の中間の図形を定義できます(k,k) 表面に

ρ=21+1(2k21)sin2(2ϕ)k4





ポイントを通過する状態を維持しながら、これらの形状にバリエーションを追加することもできます (k,k) -パラメータ自体を直接変調する k角度に応じて、対角線を通過するときにその乗数が1に等しくなります。たとえば、代わりにk 関数 k1zcos2(2ϕ)、追加のパラメータを取得します z追加の曲げを調整できます。特に、z=1/4 次のようになります。







さらに拡張されたケースでは、正方形だけでなく長方形を定義しても、極座標で定義できます。

ρ=4a2b2((b2a2)cos(2ϕ)+a2+b2)(14a2b2ksin2(2ϕ)((b2a2)cos(2ϕ)+a2+b2)2+1)



そして、その面積を(楕円積分を介して)計算することさえできます:

S=4ab(k1)K(k)+E(k)k

注意
k (0 1) , πab 4ab.


これにより、断面積が制御された円から長方形に遷移するプロファイルを作成できます。ここで、面積は一定です。





そしてここで、領域は指数関数的に拡大します。





カルテシアン座標に移動



極座標の式は、カルテシアン座標の式を介して、少なくとも2つの方法で表すことができます。図の境界の勾配の形式がどのように変化するかに応じて異なります。これを行うには、座標のアークタンジェントを通る角度を計算し、減算することによって半径ベクトルを介して式を定数にするだけで十分です。

0=x2+y2cos(2sin1(k)+πm2n)cos(2sin1(kcos(ntan1(x,y)))+πm2n)



または分割

1=x2+y2cos(2sin1(kcos(ntan1(x,y)))+πm2n)cos(2sin1(k)+πm2n)



2番目のオプションは、ポリゴンの側面に沿って直線のグラデーションを与えるため、推奨されます。



注意
, (0,0) - — , , ( cos(2sin1(k)+πm2n)sec(2sin1(kcos(πn2))+πm2n) ).



cos(ntan1(x,y)) (x+iy)n+(xiy)n2(x2+y2)n/2, A034839.


式の右辺からの式の値(2番目の場合)は、 01 ポイントの場合 (x,y) 図の中に入り、から 1adinfinitum-そうでない場合。さまざまな関数を選択して明るさに変換することで、さまざまなラスター化オプションを取得できます。指数の場合(ex1 最初と ex 2番目のオプションの場合)

または、飽和している場合



クラシックなローパスフィルターを使用できます 11+xp、 その中で p-減衰の程度を決定するフィルターの次数。



最初のオプションの場合:



そして2番目に:



補間境界を明示的に設定することにより、ピースワイズ連続関数を使用することもできます。



ラスター化自体に加えて、変形を設定することもできます。たとえば、チェッカーボードを円に押し込みます。





または、球の上に引っ張ることもできます。

x=u21+|u2v2u2+v2|



y=v1+2|u2v2u2+v2|



z=112u2(1+|u2v2u2+v2|)12v2(1+|u2v2u2+v2|)





付録:式がどのように導き出されたか



数学的なテキストでの物語の古典的なスタイルは、レンマ/定理とその証明の交代で構成されています-証明可能なステートメントが上からの啓示によって著者の頭に現れたかのように。そして、これにはいくつかの真実がありますが、多くの場合、式の出現の前にいくつかの研究が行われ、その説明は正式な証明よりもその意味をより深く理解することができます。そして、ステートメントの忠実度は、ステートメントにつながったステップの忠実度を通じて追跡できます。



したがって、ここでも-記事が複雑な形式の式で始まっている場合、その外観は自明ではなく直感に反し、宣言されたプロパティには追加の証拠が必要になります。しかし、三角測量形式の記録では、その出現の履歴を追跡することは非常に可能です。



1)最も単純なケース、つまり極座標で直線を描くタスクから始めます。これを行うには、方程式を解く必要がありますrcos(ϕ)=1その解決策は明らかです rsec(ϕ)..。





2)さらに、直線でねじれを提供するには、secant引数を「ループ」する必要があります。他のソリューションが分割残余の形で「ダーティハック」を使用するのはこの段階です。また、直接および逆正弦関数の順次取得も使用します-sin1(sin(ϕ))





このアプローチにより、結果の式に対して標準的な数学演算を実行できます。
例えば
:

sin1(sin(ϕ))ϕ=cos(ϕ)1sin2(ϕ)







同じ表記法のおかげで、三角関数の複雑な表現を使用して、極座標の二乗関数をより美的な外観に単純化することができます。Wolfram Mathematicaでは、これはTrigToExp関数とExpToTrig関数を使用して実行できます。



コード
Sec[1/2 ArcSin[k Sin[2 \[Phi]]]]^2//TrigToExp//ExpToTrig//Sqrt[#]&//FullSimplify

22+21k2sin2(2ϕ)



同じ記録のおかげで、追加の乗数を使用して、円と正方形の間の滑らかな中間図を取得できます kそのため、アークサインの議論は1つに満たない- sin1(ksin(ϕ))





そして、関数が特定のポイントと交差するためには、方程式を作成して再計算する必要があります k

21+1ksin2(π2)=k



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

k4(k21)k4





3)パラメータ n そして m創造的な方法で追加されたばかりで、実際、それらの効果は実験的に調査されました。



4)パラメトリックビューに移動し、軸を「ストレッチ」することで、長方形を簡単に取得できます。

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



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



しかしその後 t 角度を意味しなくなりました t座標軸への投影を通じてベクトルを記述する単なるパラメータです。極座標に戻るには、ベクトルの長さ(二乗和の根を通る)、角度(比率の逆正接を通る)を見つけ、この角度を次のように表す必要があります。ϕ 代わりに、結果の式を置き換えます t..。



コード
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,tan1(btan(t)a)}





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



ttan1(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(2tan1(atan(ϕ)b))+1)



このような式を単純化することはすでにより困難であり、これにはいくつかの段階が必要になります。



  1. 置き換えてカルテシアン座標に移動します ϕtan1(x,y);
  2. 指数形式に移動します。
  3. 簡素化する;
  4. 逆の交換をする xcos(ϕ) そして ysin(ϕ);
  5. 指数形式に戻ります。
  6. 簡素化する。


その結果、次の式が得られます。



コード
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((b2a2)cos(2ϕ)+a2+b2)(1+((a2+b2)cos(2ϕ)a2+b2)2((b2a2)cos(2ϕ)+a2+b2)2)









結論



ご覧のとおり、ポリゴンのような単純でありふれたものでも、何か新しいものを見つけて思いつくことができます。そして、話はそれだけではありません。一般的な場合の面積式は不明のままであり、通常のポリゴンだけでなく任意のポリゴンの式も不明のままであり、パワーおよび三角測量系列への拡張は考慮せずに残されました。3次元の場合にも同様の式が存在する可能性があります。



したがって、数学のすべてがすでに発明されており、普通の人の理解を超えた問題しかないという場合は、信じないでください。実際の数学者が気付いていない、純粋に実用的な問題がたくさんあります。または、周囲に十分な誇大宣伝がないため、または解決のためにそれらを達成する方法の大まかな考えをすでに持っているため、解決に興味がありません。ウィキペディアで解決策が入手できない問題に取り組むことを恐れないでください。解決策を公開することを恐れないでください。また、すべての役に立たないことについての記事の下にあるコメントを読むことを恐れないでください。



PSは、Mathematicaの元のドキュメントをここからダウンロードします



All Articles