はじめまして。楕円関数論に興味を持って勉強しております。
ガウスが発見した算術幾何平均の展開式、

\frac{1}{M(1+x,1-x)} = 1 + frac{1}{2}^2 x^2 + frac{3}{8}^2 x^4 + ...

の導出法をご存知の方教えてください。または、載っている本を教えて下さい。
よろしくお願い致します。

このQ&Aに関連する最新のQ&A

A 回答 (3件)

siegmund です.



前の回答は \frac の \ が抜けたり,\pi のはずが \pai になっていたり,
大分書き損なっていました.
最初はπと書いていたので,つい日本語変換に引きずられて \pai に
なってしまいました.
まあ,実害はありませんでしたが......

さて,
(1)  M(1+x,1-x) = \frac{1}{1+t^2} M(1+t^2,1-t^2)
(2)  x = \frac{2t}{1+t^2}
から出発するのでしたら(string さんのkの代わりにxと書いています),
(3)  \frac{1}{M(1+x,1-x)} = 1 + a_2 x^2 + a_4 x^4 + a_6 x^6 + ...
と書いておきます.
(1)の右辺の逆数にこれを適用して
(4)  \frac{1+t^2}{M(1+t^2,1-t^2)}
    = (1+t^2) (1 + a_2 t^2 + a_4 t^4 + a_6 t^6 + ...)
    = 1 + t^2 + a_2 t^4 + a_2 t^6 + ...
が得られます.
一方,(3)に(2)の展開形
(5)  x^2 = \frac{4t^2}{(1+t^2)^2}
      = 4t^2 - 8t^4 + 12t^6 - 16t^8 + ...
      = 4t^2 \sum_{j=1}^\infty (-1)^j j t^{2(j-1)}
を代入整理して
(6)  \frac{1}{M(1+x,1-x)}
    = 1 + 4a_2 t^2 + (-8a_2 + 16a_4)t^4 + (12a_2 - 64a_4 + 64a_6) t^6 + ...
になります.
(1)によって,(4)と(6)は等しいのですから,t^n の係数を等しいと置いて
(7)  4a_2 = 1
(8)  -8a_2 + 16a_4 = a_2
(9)  12a_2 - 64a_4 + 64a_6 = a_2
から,順次
(10)  a_2 = 1/4 = (1/2)^2
(11)  a_4 = 9/64 = (\frac{1・3}{2・4})^2
(12)  a_6 = 25/256 = (\frac{1・3・5}{2・4・6})^2
が得られます.
すなわち,
(13)  \frac{1}{M(1+x,1-x)}
     = 1 + (\frac{1}{2})^2 x^2 + (\frac{1・3}{2・4})^2 x^4
      + (\frac{1・3・5}{2・4・6})^2 x^6 + ...
(5)は展開の形が素直ですから,少し頑張れば(13)の一般項も出せそうです.

「知っていることを先に使って」いますが,
(13)はちょうど Gauss の超幾何関数
(14)  F(α,β,γ;z)
    = 1 + \frac{αβ}{1・γ}z
     + \frac{α(α+1)β(β+1)}{1・2γ(γ+1)}z^2
     + \frac{α(α+1)(α+2)β(β+1)(β+2)}{1・2・3γ(γ+1)(γ+2)}z^3 + ...
の,z = x^2 の場合になっています.
第1種完全楕円積分 K(k) は
(15)  K(k) = (π/2)F(1/2,1/2,1;k) ですから,ここらへんは前の私の回答と
つながっています.

ここからあとは,「こうなりそう」です.
確かめていませんので,間違っていたら(うまく行かなかったら)ご容赦を.
(1)(2)は関数方程式になっていますから,
微分して 1/M(1+x,1-x) に対する微分方程式が導けそうです.
全部tで書いておくと,前の因子が 1+t^2 ですから,
tの2階微分で微分方程式が作れそうです.
1/M(1+x,1-x) = Q(x^2) = Q(z) とでも置くと(z=x^2),
上の方で書いたこととのつながりから,
Q(z) の微分方程式は Gauss の超幾何微分方程式
(16)  z(1-z)Q''(z) + {γ-(α+β+1)z}Q'(z) - αβz = 0
のα=β=1/2,γ=1 になるはずと思います.
(16)の解を級数展開で求めたものが(14)です.
    • good
    • 0

siegmund です.


またミスプリしちゃいました.
下記のように訂正します.

-------------

(5)  x^2 = \frac{4t^2}{(1+t^2)^2}
      = 4t^2 - 8t^4 + 12t^6 - 16t^8 + ...
      = 4t^2 \sum_{j=1}^\infty (-1)^{j+1} j t^{2(j-1)}

-------------

(13)はちょうど Gauss の超幾何関数
(14)  F(α,β,γ;z)
    = 1 + \frac{αβ}{1・γ}z
     + \frac{α(α+1)β(β+1)}{1・2γ(γ+1)}z^2
     + \frac{α(α+1)(α+2)β(β+1)(β+2)}{1・2・3γ(γ+1)(γ+2)}z^3 + ...
の,α=β=1/2,γ=1 ,z = x^2 の場合になっています.
    • good
    • 0
この回答へのお礼

siegmundさん、再度ご回答いただき誠にありがとうございます。

展開して係数を決めるというやり方でやればよいのですね。

私は

\frac{1}{M(1+x,1-x)} = (1+t^2) \frac{1}{M(1+t^2,1-t-2)}

を繰り返し適用するというやり方でやろうとしていました。

超幾何関数と対応関係もよくわかりました。
本当にどうもありがとうございました。

お礼日時:2001/03/09 05:44

a,b の算術幾何平均 M(a,b) は


(1)   frac{1}{M(a,b)}
     = \frac{2}{\pai} \int_0^{\pai/2}
      frac{d\theta}{\sqrt{a^2 \cos + b^2 \sin^2\theta}}
で表されます.
(1)の分母でaを外に引っ張り出すなどしますと
(2)   frac{1}{M(a,b)}
     = frac{2 K'(b/a)}{\pai a}
になります.
K'(k) = K(k') = K(\sqrt{1-k^2}) で,K は第一種完全楕円積分
(3)   K(k) = \int_0^{\pai/2} frac{d\theta}{\sqrt{1 - k^2 \sin^2 \theta}}
です.
あとは,a = 1+x,b = 1-x と,K(k) の展開式
(4)   K(k) = \frac{\pai}{2} \sum_{r=0}^\infty
        [ \frac{(2r-1)!!}{(2r)!!} ]^2 k^{2r}
を整理すれば,string さんの式が出てきます.

でも,(1)はどうやるんでしたっけ?
楕円 \vartheta_j 関数(j=0~3)の \vartheta_j(0|2\tau) を
\vartheta_n(0|\tau) で表す公式のどれかが算術平均と幾何平均の形になっていて,
そこがキーポイントになっていたと思います.
ちょっとよく思い出せません(というか,私の知識はそんなもの).

どなたか楕円関数自由自在という方のフォローがあるといいんですが....
    • good
    • 0
この回答へのお礼

ご回答を頂き、誠にありがとうございます。

算術幾何平均と楕円積分の対応を認めれば、なるほど簡単な議論で導出されるのですね。

ガウスは、算術幾何平均の性質、

M(1+k,1-k) = \frac{1}{1+t^2} M(1+t^2,1-t^2)
ただし、 k = \frac{2t}{1+t^2}

から、展開式を導いたようです。それと、楕円積分を二項展開して積分した式が一致することから、(レムニスケート積分の時に発見した)楕円積分と算術幾何平均との対応を改めて確認していったようです。

算術幾何平均の性質から導出したいと知恵をしぼっております。
また、お助けください。

お礼日時:2001/03/07 23:21

このQ&Aに関連する人気のQ&A

お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!

このQ&Aを見た人が検索しているワード

このQ&Aと関連する良く見られている質問

Qexp{L[1]x+L[2]x^2/2+L[3]x^3/3+…}=F[1]+F[2]x+F[3]x^2+…

フィボナッチ数列F[n]は、
F[1]=1,F[2]=1,F[n+2]=F[n+1]+F[n]
で定義され、リュカ数列L[n]は、
L[1]=1,L[2]=3,L[n+2]=L[n+1]+L[n]
で定義されます。このとき、

exp{L[1]x+L[2]x^2/2+L[3]x^3/3+…}=F[1]+F[2]x+F[3]x^2+…

が成り立つそうなのですが、どうしてなのですか?

右辺は、フィボナッチ数列の母関数と似ていてなんとか求められるのですが、左辺をどうして求めていいかわかりません。

なお、式は
http://mathworld.wolfram.com/FibonacciNumber.html
の(68)を参照しました。

Aベストアンサー

↓ここに証明がありますね。
http://maths.dur.ac.uk/~dma0rcj/PED/fib.pdf
(2.7 A surprising sum を見てください。)

参考URL:http://maths.dur.ac.uk/~dma0rcj/PED/fib.pdf

Q∫【1→2】{(x^2-x+4)/x(x^3+1)}dx

∫【1→2】{(x^2-x+4)/x(x^3+1)}dxという定積分の求め方がわかりません。
私はまず部分分数に分けて、
(x^2-x+4)/x(x^3+1)
=4/x-(4x^2-x+1)/(x^3+1)として、
∫【1→2】{(x^2-x+4)/x(x^3+1)}dx
=(16/3)*log2-(8/3)*log3+【1→2】∫(x-1)/(x^3+1)dx
というところまで求めたのですが、最後の定積分が求められず、ここで手が止ま
ってしまいました。
ちなみに最終的な答えは3*log(4/3)となるそうです。問題集には答えしか書か
れてないので困っています(^_^;)

Aベストアンサー

x^3+1=(x+1)(x^2-x+1)
でさらに部分分数分解を行ってください

Q∫∫【D】2x|y|dxdy, D={x^2+y^2≦1,x^2+y^2≦2x}

∫∫【D】2x|y|dxdy, D={x^2+y^2≦1,x^2+y^2≦2x}
という重積分について質問です。∫∫【D】2x|y|dxdyと∫∫【D】2xydxdyってどう違いますか?

この場合では、領域がx軸に関して対称だから、前者の場合も後者の場合もたまたま答えが同じになるけれど、理屈としては、y座標が負になっている部分をx軸に関して折り曲げた結果として、図形がx軸に関して対称だったために、y座標が正の部分を2倍することになったと考えればよいのでしょうか?
言葉が下手で、伝わりにくい文章ですみません。

Aベストアンサー

>この場合では、領域がx軸に関して対称だから、前者の場合も後者の場合もたまたま答えが同じになるけれど

本当にそうなります?
2xyはyについて奇関数、2x|y|はyについて偶関数です。
前者をx軸について対称な領域で積分すると"0"に、後者を同じ領域で積分するとx軸よりも上側の領域での積分の2倍になります。

QF_n=(a+b+c)^(2n+1)-{a^(2n+1)+b^(2n+1)+c^(2n+1)} の因数分解

F_n=(a+b+c)^(2n+1)-{a^(2n+1)+b^(2n+1)+c^(2n+1)} 
(n=1,2,3,4,5)
を因数分解せよ、という問題なのですが、どすればよいのでしょうか?

なお、答えは、

F_1=3(b+c)(c+a)(a+b)
F_2=5(b+c)(c+a)(a+b)(Σa^2+Σab)
F_3=7(b+c)(c+a)(a+b)(Σa^4+2Σa^3 b+3Σa^2 b^2+5Σa^2 bc)
F_4=3(b+c)(c+a)(a+b)(3Σa^6+9Σa^5 b+19Σa^4 b^2+35Σa^4 bc+23Σa^3 b^3+63Σa^3 b^2 c)
F_5=11(b+c)(c+a)(a+b)(Σa^8+4Σa^7 b+11Σa^6 b^2+21Σa^6 bc+9Σa^5 b^3+54Σa^5 b^2 c+23Σa^4 b^4+84Σa^4 b^3 c+123Σa^4 b^2 c^2+159Σa^3 b^3 c^2)

のようなのですが、(b+c)(c+a)(a+b)を因数に持つことは分かりますが、残りの因数はどうやってもとめるのでしょうか?

一文字を変数と見て、地道に割り算するしかないのでしょうか?
効率的な計算方法はありますでしょうか?

F_n=(a+b+c)^(2n+1)-{a^(2n+1)+b^(2n+1)+c^(2n+1)} 
(n=1,2,3,4,5)
を因数分解せよ、という問題なのですが、どすればよいのでしょうか?

なお、答えは、

F_1=3(b+c)(c+a)(a+b)
F_2=5(b+c)(c+a)(a+b)(Σa^2+Σab)
F_3=7(b+c)(c+a)(a+b)(Σa^4+2Σa^3 b+3Σa^2 b^2+5Σa^2 bc)
F_4=3(b+c)(c+a)(a+b)(3Σa^6+9Σa^5 b+19Σa^4 b^2+35Σa^4 bc+23Σa^3 b^3+63Σa^3 b^2 c)
F_5=11(b+c)(c+a)(a+b)(Σa^8+4Σa^7 b+11Σa^6 b^2+21Σa^6 bc+9Σa^5 b^3+54Σa^5 b^2 c+23Σa^4 b^4+84Σa^4 b^3 c+123Σa^4 b^2 c^2+159Σa^3 b^3 c^...続きを読む

Aベストアンサー

最後までは計算していませんが、次の方法でできそうです。
F_n = (b+c)(c+a)(a+b)(Σ[ABC] k_ABC a^A b^B c^C) とおきます。
(ここで、A+B+C = 2n+1 です。)
展開すると、F_n = (a^2 b + 略 + 2abc)(Σ[ABC] k_ABC a^A b^B c^C) です。
そして、F_n を例えば、a で A+2 回偏微分、a で B+1 回偏微分、
a で C 回偏微分、した後、a,b,c に 0 を代入します。
F_n=(a+b+c)^(2n+1)-{a^(2n+1)+b^(2n+1)+c^(2n+1)} に対しても同じようにします。
このようにすると、例えば C > 0 であれば、
k_ABC (A+2)!(B+1)!(C)! = (2n+1)! となり、係数が得られます。

Q「(5x+3)^10でx^pとx^(p+1)の係数比が21:20になる時のpの値」と「x+y=1を満たす全x,yに対してax^2+2bxy+by^2

こんにちは。識者の皆様、宜しくお願い致します。

[問1] (5x+3)^10の展開式でx^pとx^(p+1)の係数比が21:20になる時のpの値を求めよ。
[問2]x+y=1を満たす全てのx,yに対して
ax^2+2bxy+by^2+cx+y+2=0が成立するように定数a,b,cの値を定めよ。

[1の解]
(5x+3)^10=10Σk=0[(10-k)Ck 5x^(10-k)3^k]なので
p=10-kの時(k=10-pの時)
p+1=10-kの時(k=9-pの時)より
a:b=pC(10-p) 5^p 3^(10-p):(1+p)C(9-p) 5^(1+p) 3^(9-p)
で 1/(10-p):(1+p)/(2p-8)/(2p-9)=7:4 から
23p^3-199p+218=0
となったのですがこれを解いてもp=6(予想される解)が出ません。
やり方が違うのでしょうか?

[2の解]
与式をx+yという対称式で表せばならないと思います(多分)。
どうすれば対称式で表せるのでしょうか?

Aベストアンサー

 (1)Cをばらして比を簡略化するところで計算間違いがありそうな気がします。その経過をもう少し詳しく書いてもらえませんか?
 (2)a,b,cを求めるにはまず、x+y=1 を満たすすべての(x,y)で成り立つのですから、x+y=1を満たす(x,y)をまず代入してみてはどうでしょうか。候補としては、(1,0)(0,1)(2,-1)など。
 それから計算されたa,b,c でx+y=1を満たすすべてのx,yで成り立つかどうかを確認するという手順でどうでしょうか?


人気Q&Aランキング

おすすめ情報