■レムニスケート積分(その6)
∫1/(1-x^2)^(1/2)dx
は円(2次曲線),
∫1/(1-x^4)^(1/2)dx
はレムニスケート(4次曲線)に対応していますが,周長が
∫1/(1-x^3)^(1/2)dx
で表される曲線はどのようなものになるでしょうか?
===================================
【1】変形バラ曲線
より一般に,「周長が
∫1/(1-x^n)^(1/2)dx
で表される曲線は
r^(n/2)=cos(n/2・θ)
である.」
(証)
{1+(rdθ/dr)^2}^(1/2)=1/(1-r^n)^(1/2)
rdθ/dr=(r^n/(1-r^n))^(1/2)
dθ/dr=(r^n-2/(1-r^n))^(1/2)
dr/dθ=((1-r^n)/r^n-2)^(1/2)
一方,
r^(n/2)=cos(n/2θ)
n/2・r^(n/2-1)dr/dθ=n/2・sin(n/2θ)
r^(n/2-1)dr/dθ=sin(n/2θ)
dr/dθ=sin(n/2θ)/r^(n/2-1)=((1-r^n)/r^n-2)^(1/2)
この曲線はバラ曲線(正葉曲線)の変形版になっている.円(n=2)もレムニスケート(n=4)もこの曲線族に属する.n=1の場合,曲線
r^(1/2)=cos(1/2・θ)
r=cos^2(1/2・θ)=(1+cosθ)/2
はカージオイドとなる.
なお,nが奇数とき4n次代数曲線,nが偶数のときn次代数曲線になることが示される.
===================================
【2】弧長積分
f(x)=1/(1-x^2)^(1/2)
のとき,
sin^(-1)z=∫(0,z)f(x)dx
であるから,
2∫(0,1)f(x)dx=3.141592・・・=π
∫(0,1)f(x)dx=π/2
となる.それでは,
f(x)=1/(1-x^4)^(1/2)
としたとき,
∫(0,1)f(x)dx=1.311028・・・=ω/2
は,どのようにすれば得られるのでしょうか?
ベータ関数において,a=m/n,b=1/2とおき,t=x^nと置換すると,
∫(0,1)x^(m-1)/(1-x^n)^(1/2)dx=Γ(m/n)√π/nΓ(m/n+1/2)
したがって,
(m,n)=(1,1)のとき,∫(0,1)1/(1-x^1)^(1/2)dx=2
(m,n)=(1,2)のとき,∫(0,1)1/(1-x^2)^(1/2)dx=π/2
(m,n)=(1,3)のとき,∫(0,1)1/(1-x^3)^(1/2)dx=Γ^3(1/3)/2^(4/3)3^(1/2)π
(m,n)=(1,4)のとき,∫(0,1)1/(1-x^4)^(1/2)dx=Γ^2(1/4)/2^(5/2)π^(1/2)
が得られます.
∫(0,1)1/(1-x^1)^(1/2)dx=2
∫(0,1)1/(1-x^2)^(1/2)dx=π/2
は初等的にも得ることができます.一方,
∫(0,1)1/(1-x^3)^(1/2)dx=Γ^3(1/3)/2^(4/3)3^(1/2)π
∫(0,1)1/(1-x^4)^(1/2)dx=Γ^2(1/4)/2^(5/2)π^(1/2)
は,特別な数と楕円積分を関係づけるものになっています.
なお,
∫(0,1)1/(1-x^3)^(1/2)dx=Γ^3(1/3)/2^(4/3)3^(1/2)π
を得るには,ガンマ関数の乗法公式(倍数公式)
Γ(x/2)Γ((x+1)/2)=π^(1/2)Γ(x)/2^(x-1)
と相反公式(相補公式)
Γ(x)Γ(1-x)=π/sinπx
また,
∫(0,1)1/(1-x^4)^(1/2)dx=Γ^2(1/4)/2^(5/2)π^(1/2)
を得るには乗法公式を用いています.
[1]n=1:Γ(1)√π/Γ(3/2)=2
[2]n=2:Γ(1/2)√π/2Γ(1)=π/2=1.5708
[3]n=3:Γ(1/3)√π/3Γ(5/6)=1.40218
[4]n=4:Γ(1/4)√π/4Γ(3/4)=1.31103
[5]n=5:Γ(1/5)√π/5Γ(7/10)=1.25373
[6]n=6:Γ(1/6)√π/6Γ(2/3)=1.21433
[7]n→∞のとき,
∫(0,1)1/(1-x^n)^(1/2)dx=Γ(1/n)Γ(1/2)/nΓ(1/n+1/2)→Γ(1/n)/n=Γ(1+1/n)→Γ(1)=1
この結果は図形的に考察すれば明らかである.
===================================
【3】第1種楕円積分の標準形への還元
積分
I=∫(c,1)dx/(1−x^3)^1/2
を考える.変換
y=(λx+μ)/(νx+ρ),x=(−ρy+μ)/(νy−λ)
を適当に選ぶことにより,常に
dy/(1+my^2+ny^4)^1/2
の形に還元できる.
λ=−1,μ=1−√3,ν=1,ρ=−1+√3
x=cに対するyの値を
y1=(c−1+√3)/(−c+1+√3)
とおくと
I=1/4√3∫(y1,1)2dy/((1−y^2)(2−√3+(2+√3)y^2)^1/2
ここで,y^2=1−z^2,z1^2=1−y1^2とおくと第1種楕円
I=1/4√3∫(0,z1)dz/((1−z^2)(1−k^2z^2)^1/2=1/4√3F(k,φ1)
となる.ただし,
k=(√2+√6)/4,k^2=(2+√3)/4
φ1=arccos((c−1+√3)/(−c+1+√3))=arcsin((1−y1^2)^1/2)
sinφ1=z1,z1=(1−y1^2)^1/2
c=1→z1=0
c=0→z1=(4√3−6)^1/2=α
ところで,第1種楕円積分の標準形
F(k,x)=∫(0,x)1/{(1-x^2)(1-k^2x^2)}^(1/2)f(x)dx
については加法公式が使える.これで,この曲線の弧長を求める問題は,加法公式
z={x(1+my^2+ny^4)^1/2+y(1+mx^2+nx^4)^1/2}/(1−nx^2y^2)
m=−(1+k^2),n=k^2,α=(4√3−6)^1/2
の問題に還元されたことになる.
倍角公式
x’=2xy/(1−nx^4)=2x(1+mx^2+nx^4)^1/2/(1−nx^4)
より,2等分点に対応する楕円曲線上の点
2x(1+mx^2+nx^4)^1/2/(1−nx^4)=α
の解をx=βとすると,
1−((c−1+√3)/(−c+1+√3))^2=β^2
c=√3+1−2√3/((1−β^2)^1/2+1)
がこの曲線の2等分点である.
2x(1+mx^2+nx^4)^1/2/(1−nx^4)=α
を求めたところ,解析解
x=−1+√3
が得られた.
したがって,2等分点は
c=√3+1−(2√3+3)(1−√(2√3−3)=0,671618(→作図可能).4等分点,8等分点も同様に計算することができる.
===================================