楕円積分

積分学において楕円積分は、特定の積分の値として定義されるいくつかの関連関数の一つであり、ジュリオ・ファニャーノレオンハルト・オイラー 1750年頃)によって初めて研究されました。その名称は、楕円弧長を求める問題との関連性に由来しています

現代数学では、「楕円積分」を、次の式で表される関数 fと定義しています。

ここで、R2 つの引数の有理関数、 P重根を持たない 3 次または 4 次多項式、 cは定数です。

一般に、この形式の積分は初等関数で表現することはできない。この一般規則の例外は、P が重根を持つ場合、R ( x , y )がyの奇数乗を含まない場合、および積分が擬楕円型である場合である。しかし、適切な簡約公式を用いることで、すべての楕円積分は、有理関数上の積分と、第一種、第二種、第三種の楕円積分とも呼ばれる三つのルジャンドル標準形を含む形式に変換することができる

楕円積分は、以下に示すルジャンドル形式に加えて、カールソン対称形式でも表すことができます。楕円積分の理論へのさらなる理解は、シュワルツ・クリストッフェル写像の研究を通して得られます。歴史的に、楕円関数は楕円積分の逆関数として発見されました。

引数表記

不完全楕円積分は2つの引数を持つ関数であり、完全楕円積分は1つの引数を持つ関数です。これらの引数は同じ楕円積分を与えるため、様々な異なる方法で表現されますが、いずれも等価です。ほとんどの教科書では、以下の命名規則を用いた標準的な命名規則が採用されています。

1つの引数を表す場合:

上記の3つの量は、それぞれ他の量によって完全に決定されます(ただし、それらが非負値である場合)。したがって、これらは互換的に使用できます。

他の引数も同様にφ(振幅)またはxまたはu ( x = sin φ = sn usnヤコビ楕円関数の 1 つ)として表すことができます

これらの量のいずれかの値を指定すると、他の値も決定されます。umにも依存することに注意してください。u関連するその他の関係には以下が含まれます 。

後者はデルタ振幅と呼ばれることもあり、 Δ( φ ) = dn uと記されます。文献によっては、相補パラメータ相補モジュラス、または相補モジュラー角と呼ばれることもあります。これらは、四分の一周期に関する記事でさらに定義されています

この表記法では、縦棒を区切り記号として使用すると、それに続く引数が(上記で定義した)「パラメータ」であることを示し、バックスラッシュはモジュラー角であることを示します。セミコロンの使用は、それに続く引数が振幅の正弦であることを意味します。このように異なる引数区切り記号を使用するのは混乱を招く可能性があり、楕円積分では伝統的であり、表記法の多くは、AbramowitzとStegunの参考書やGradshteynとRyzhikの積分表で使用されている表記法と互換性があります

文献では、楕円積分の表記法として他の慣例も用いられています。引数を交換した表記法F ( k , φ )がよく用いられ、同様に第二種積分についてはE ( k , φ )が用いられます。Abramowitzと Stegun は第二種と第三種の積分の定義において、引数 φ の後に縦棒が続かない限り、第一種積分F ( φ , k )を引数φの代わりに用いています。つまり、 E ( φ | k 2 )ではなくE ( F ( φ , k ) | k 2 )としています。さらに、彼らの完全積分では係数kの代わりに引数としてパラメータ k 2を用いています。つまりK ( k )ではなくK ( k 2 )としています。またGradshteynとRyzhikによって定義された第3種の積分Π( φ , n , k )では、振幅φが優先され、「特性」nは優先されません。

したがって、これらの関数を使用する際には表記法に注意する必要があります。なぜなら、様々な信頼できる参考文献やソフトウェアパッケージでは、楕円関数の定義において異なる表記法が用いられているからです。例えば、Wolfram社のMathematicaソフトウェアとWolfram Alphaは、第一種完全楕円積分を楕円係数kではなく、パラメータmで定義しています。

第一種不完全楕円積分

第一種不完全楕円積分 F次のように定義される。

これはルジャンドルの楕円積分の三角関数形式である。t = sin θおよびx = sin φを代入すると、ヤコビの代数形式が得られる。

同様に、振幅とモジュラー角に関しては次のようになります。

x = sn( u , k )の場合、次の式が得られます。このヤコビ楕円関数は、第一種不完全楕円積分の単純な逆関数であることが証明されます。

第一種不完全楕円積分には次の加法定理がある[要出典]

楕円係数は次のように変換できます。

第二種不完全楕円積分

ルジャンドルの三角関数形式における第二種不完全楕円積分 E

t = sin θおよびx = sin φを代入すると、ヤコビの代数形式が得られます。

同様に、振幅とモジュラー角の観点からは、次のようになります。

ヤコビ楕円関数との関係としては、

赤道から緯度φまでの子午線弧の長さは、 Eで表されますここで、 a長半径eは離心率です

第二種不完全楕円積分には次の加法定理がある[要出典]

楕円係数は次のように変換できます。

第三種不完全楕円積分

第三種不完全楕円 積分Π

または

nは特性値と呼ばれ、他の引数とは独立して任意の値を取ることができます。ただし、Π(1; π/2 | m )は、任意のmに対して無限大です

ヤコビ楕円関数との関係は

赤道から緯度φまでの子午線弧の長さもΠの特殊なケースと関連しています

第一種完全楕円積分

第一種完全楕円積分K ( k )のプロット

楕円積分は、振幅φ = ⁠のとき「完全」であると言われる。π/2であり、したがってx = 1 である。第一種完全楕円積分 K は、このように定義される。あるいは、第一種不完全積分を用いてより簡潔に

これはべき級数として表すことができる

ここでP nはルジャンドル多項式であり、これは次式と等価である。

ここでn !!は二重階乗を表す。ガウスの超幾何関数を用いると、第一種完全楕円積分は次のように表される。

第一種完全楕円積分は、四分の一周期と呼ばれることもある。これは算術幾何平均を用いて非常に効率的に計算することができる[1]

したがって、係数は次のように変換できます。

この式は0 ≤ k ≤ 1のすべてに対して有効です

ガンマ関数との関係

k 2 = λ ( i r )かつλモジュラーラムダ関数ならば、 K ( k )はガンマ関数を用いて閉じた形で表現できる[2]例えば、r = 2r = 3r = 7の場合には、それぞれ[3]となる。

そして

そして

より一般的には、虚二次体[注1]に含まれるという条件で十分である。[4] [5]例えば、k = e 5 πi /6のとき、iK /K = e 2 πi /3および[6]

漸近表現

この近似値は、3 × 10 −4k < ⁠)1/2最初の2項のみを保持すると、 k < の場合には0.01の精度で正しい1/2 . [要出典]

微分方程式

第一種楕円積分の微分方程式は

この方程式の2番目の解は である。この解は次の関係を満たす。

連分数

連分数展開は次の通りある: [7] ここでnomeはその定義の中にある

期間比の反転

ここでは、代わりにパラメータ を 用いた第一種完全楕円積分を用いる。これは、複素平面上で逆関数を求める際に平方関数が問題を引き起こすためである。そこで、

そして

シータ関数となります

方程式

(解が存在すると仮定して) 次のように解くことができる。

これは実際にはモジュラーラムダ関数です。

計算の目的のために、誤差解析は[8]で与えられる。

ここで、および

また

どこ

第二種完全楕円積分

第二種完全楕円積分E ( k )のプロット

第二種完全楕円積分 E次のように定義される。

あるいは、第二種不完全積分​​E ( φ , k )を用いてより簡潔に表すと、

長半径aと短半径bを持ち、離心率e = 1 − b 2 / a 2である楕円の場合、第二種完全楕円積分E ( e )は、長半径aを単位として測った楕円の円周 Cの4分の1に等しい。言い換えれば、

第二種完全楕円積分は冪級数として表すことができる[9]

これは次の式と同等である。

ガウスの超幾何関数を用いると、第二種完全楕円積分は次のように表される。

係数は次のように変換できます。

計算

第一種積分と同様に、第二種完全楕円積分は算術幾何平均を用いて非常に効率的に計算することができる。[1]

数列a ng nを定義する。ここでa 0 = 1g 0 = 1 − k 2 = k であり、再帰関係a n + 1 = a n + g n/2 , g n + 1 = a n g nが成り立つ。さらに、定義する。

定義上、

また

それから

実際には、算術幾何平均はある限界まで単純に計算されます。この式は、| k | ≤ 1のすべてにおいて二次収束します。計算をさらに高速化するために、c n + 1 = ⁠という関係式が用いられます。c n 2/4 a n + 1使えます。

さらに、k 2 = λ ( i r )かつλモジュラーラムダ関数)とすれば、E ( k )は を用いて閉じた形で表現できるため、無限和項を必要とせずに計算できる。例えば、r = 1r = 3r = 7の場合には、それぞれ次式が得られる[10] 。

そして

そして

微分方程式と微分方程式

この方程式の2番目の解はE ( √1 k2 ) −K ( √1 k2 )です

第三種完全楕円積分

第三種完全楕円積分Π( n , k )のプロット。nの値はいくつか固定されている

第三種完全楕円積分 Π次のように定義される。

第三種楕円積分は、特性 nの符号が逆になる形で定義されることもあるので注意が必要である。

第一種および第二種の完全楕円積分と同様に、第三種の完全楕円積分は算術幾何平均を用いて非常に効率的に計算することができる。[1]

偏微分

ヤコビゼータ関数

1829年、ヤコビはヤコビゼータ関数を定義しました。これは において周期的であり、最小周期です。これはヤコビzn関数との関係にあります。文献(例えばWhittaker and Watson (1927))では、 はWikipediaの を意味する場合があります。また、一部の著者(例えばKing (1924))はWikipediaの と両方に を使用しています

ルジャンドルの関係

ルジャンドルの関係式またはルジャンドル恒等式は、楕円係数の積分KとEとその反関係積分[11] [12]の関係を2次積分方程式で表したものである。

互いにピタゴラス対応である 2 つのモジュールの場合、次の関係は有効です。

例えば:

また、互いに接線方向に対応する 2 つのモジュールの場合、次の関係が有効です。

例えば:

接線モジュラー対応物に対するルジャンドルの関係は、ピタゴラス対応モジュラー上のランデンモジュラー変換を使用することにより、ピタゴラス対応物に対するルジャンドルの恒等式から直接得られます。

レムニスカティックケースの特別な恒等式

レムニスカティックの場合、楕円係数または比離心率 ε は2の平方根の半分に等しい。レムニスカティックの場合のルジャンドルの恒等式は次のように証明できる。

連鎖律によれば、これらの導関数は次のようになります。

微積分学の基本定理を使用すると、次の式を生成できます。

今述べた 2 つの積分の線形結合により、次の式が得られます

積の法則を使用して、ここに示した関数から x に関連する元の不定積分を形成すると、次の式が得られます。

この整数恒等式に値を挿入すると、次の恒等式が出現します。

ルジャンドルの等式からのこのレムニスカティックな抜粋は次のようになります。

全体的なケースの一般化

ここで、モジュラー一般の場合[13] [14] が解明される。この目的のために、完全楕円積分の導関数はモジュラスの後に導出され、それらを結合する。そして、ルジャンドルの恒等式バランスが決定される。

円関数の導関数は、同一の写像関数と円関数の逆数の負の積であるため、

これらは、この記事の上のセクションで示した K と E の導関数です。

円関数の導関数と組み合わせると、これらの導関数は有効になります。

ルジャンドルの恒等式は、任意の2つの完全楕円積分の積を含む。ルジャンドルの恒等式の等式スケールから関数の辺を導出するために、積の法則を以下のように適用する。

これら 3 つの方程式のうち、上の 2 つの方程式を加算し、下の方程式を引くと、次の結果が得られます。

方程式のバランスに関しては、常に値ゼロが返されます。

以前に決定された結果は、前のセクションで計算された係数のルジャンドル方程式と組み合わせる必要があります。

最後の 2 つの式を組み合わせると、次の結果が得られます。

連続関数の微分が常にゼロの値を取る場合、その関数は定数関数です。つまり、この関数は各横軸の値に対して同じ関数値を持ち、したがって関数グラフは水平な直線になります。

参照

参考文献

注記

  1. ^ K は複素平面解析的に拡張できる

参考文献

  1. ^ abc カールソン2010、19.8。
  2. ^ ボルウェイン、ジョナサン・M.; ボルウェイン、ピーター・B. (1987).円周率とAGM:解析的数論と計算複雑性の研究(初版). Wiley-Interscience. ISBN 0-471-83138-7296ページ
  3. ^ ボルウェイン、ジョナサン・M.; ボルウェイン、ピーター・B. (1987).円周率とAGM:解析的数論と計算複雑性の研究(初版). Wiley-Interscience. ISBN 0-471-83138-7298ページ
  4. ^ Chowla, S.; Selberg, A. (1949). 「エプスタインのゼータ関数について (I)」. Proceedings of the National Academy of Sciences . 35 (7): 373. Bibcode :1949PNAS...35..371C. doi :10.1073/PNAS.35.7.371. PMC 1063041. PMID 16588908.  S2CID 45071481  . 
  5. ^ チョウラ、S.セルバーグ、A. (1967)。 「エプスタインのゼータ関数について」。Reine und Angewandte Mathematik に関するジャーナル227 : 86-110 .
  6. ^ 「ルジャンドルの楕円積分(エントリ175b7a)」。
  7. ^ N.Bagis,L.Glasser.(2015)「ラマヌジャンの連分数の評価」Rend.Sem.Mat.Univ.Padova, Vol.133 pp 1-10
  8. ^ 「ヤコビ・シータ関数の近似値」『数学関数グリモア』フレドリック・ヨハンソン. 2024年8月29日閲覧
  9. ^ 「第2種の完全楕円積分:級数表現(式08.01.06.0002)」。
  10. ^ ボルウェイン、ジョナサン・M.; ボルウェイン、ピーター・B. (1987).円周率とAGM:解析的数論と計算複雑性の研究(初版). Wiley-Interscience. ISBN 0-471-83138-726ページ、161ページ
  11. ^ "Legendre-Relation" (ドイツ語) . 2022年11月29日閲覧
  12. ^ 「Legendre Relation」 . 2022年11月29日閲覧
  13. ^ 「積分 - 楕円曲線のルジャンドル関係式の証明」 。 2023年2月10日閲覧
  14. ^ インターネットアーカイブ (1991)、ポール・ハルモス数学生誕50周年記念、ニューヨーク:シュプリンガー・フェアラーク、ISBN 0-387-97509-82023年2月10日取得

出典

  • 「楕円積分」数学百科事典EMSプレス、2001 [1994]
  • Eric W. Weisstein、「楕円積分」(Mathworld)
  • 楕円積分評価のためのMATLABコード(楕円プロジェクト)
  • 完全楕円積分の有理近似(Exstrom Laboratories)
  • 楕円積分の加法定理の簡単な歴史
Retrieved from "https://en.wikipedia.org/w/index.php?title=Elliptic_integral&oldid=1314291870"