Quantcast
Channel: 円周率近似値の日に生まれて理系じゃないわけないだろ! - knifeのblog
Viewing all articles
Browse latest Browse all 5376

三角関数の積の不思議

$
0
0

最近、FBの数学のコミュで、三角関数がたくさん出てくる式をよく見かけるようになった。

問題

sin(π/13)sin(3π/13)sin(4π/13)
sin(2π/13)sin(5π/13)sin(6π/13)

の厳密解を求めよ。

sin(θ)=sin(π-θ)
より、
分子のsinの中の分子のπの係数が奇数のものを偶数にすると

 

与式=
sin(12π/13)sin(10π/13)sin(4π/13)
sin(2π/13)sin(5π/13)sin(6π/13)

と変形出来る。
分母のsinの中の分子がすべて偶数になったので、
sin(2θ)=2sin(θ)cos(θ)
と変形でき、
 

=
{2sin(6π/13)cos(6π/13)}{2sin(5π/13)cos(5π/13)}{2sin(2π/13)cos(2π/13)}
sin(2π/13)sin(5π/13)sin(6π/13)
=8cos(6π/13)cos(5π/13)cos(2π/13)

と、分数も消え、三角関数の数が半減しました。

ここまでは、結構容易に求まったのですが、これが答えではなくて、ここから厳密解を求めなければならない。

オイラーの公式
e=cos(θ)+i・sin(θ)
から、ド・モアブルの定理
(cos(θ)+isin(θ))n=cos(nθ)+i・sin(nθ)
を導出するか、ド・モアブルの定理を知っているものとして使う。

z=cos(2π/13)+isin(2π/13)
とおくと、
z13=cos(2π)+isin(2π)=1
と表すことが出来る。

cos(2π/13)=z+1/z
cos(5π/13)=-cos(8π/13)=-cos(4・2π/13)=-(z4+1/z4)
cos(6π/13)=cos(3・2π/13)=z3+1/z3

なので、

与式=-8cos(6π/13)cos(8π/13)cos(2π/13)
=-8(z4+1/z4)(z3+1/z3)(z+1/z)
=-8(z7+z+1/z+1/z7)(z+1/z)
=-8(z8+z2+1+1/z6+z6+1+1/z2+1/z8)
=-8(z8+z6+z2+2+1/z2+1/z6+1/z8)

z13=1より、
1/z2=z11
1/z6=z7
1/z8=z5

となり、

与式=-8(z11+z8+z7+z6+z5+z2+2) …(1)

また、分母のsinの中の分子のπの係数を奇数から偶数に変換した
 

与式=
sin(π/13)sin(3π/13)sin(4π/13)
sin(2π/13)sin(8π/13)sin(6π/13)
=
sin(π/13)sin(3π/13)sin(4π/13)
{2sin(π/13)cos(π/13)}{2sin(4π/13)cos(4π/13)}{2sin(3π/13)cos(3π/13)}
=
1
8cos(π/13)cos(4π/13)cos(3π/13)

と逆数を作る。
 

cos(π/13)=cos(12π/13)=cos(6・2π/13)=z6+1/z6
cos(4π/13)=cos(2・2π/13)=z2+1/z2
cos(3π/13)=cos(10π/13)=cos(5・2π/13)=z5+1/z5

なので、

与式=1/(8(z6+1/z6)(z2+1/z2)(z5+1/z5))
=1/(8(z8+z4+1/z4+1/z8)(z5+1/z5))
=1/(8(z13+z9+z+1/z3+z3+1/z+1/z9+1/z13))

z13=1より、
1/z3=z10
1/z=z12
1/z9=z4
1/z13=1

となり、

与式=1/(8(z12+z10+z9+z4+z3+z+2)) …(2)

(1)式を8で割った値を
u=-(z11+z8+z7+z6+z5+z2+2)
とおくと、
1/u=(z12+z10+z9+z4+z3+z+2)
となり、
-u+1/u=z12+z11+z10+z9+z8+z7+z6+z5+z4+z3+z2+z+4
となる。
ここで、
z13=1
z13-1=0
より、
(z-1)(z12+z11+z10+z9+z8+z7+z6+z5+z4+z3+z2+z+1)=0
となって、
z12+z11+z10+z9+z8+z7+z6+z5+z4+z3+z2+z=-1
となり、
-u+1/u=3

また、
-u×(1/u)=-1
なので、

2次方程式の解と係数の関係
-u+1/u=3
-u×(1/u)=-1
より、
x2-3x-1=0
の解となり、
x=(3±13)/2

-u=(3+13)/2
または、
-u=(3-13)/2

問題文の式のsin(kπ/13)は、k<(13-1)/2なので、0<kπ/13<πの範囲となり、
u>0
となり、
(-3+13)/2
のように、かなり面倒だが厳密解を求めることが出来た。

問題のように三角関数の角度の係数が歯抜けの等差数列になっているようなものの積が、厳密解となるのが面白いですよね。

このような性質を利用すると、
 

n
Π
k=1
sin(kπ/(2n+1))=2n+1/2n
n
Π
k=1
cos(kπ/(2n+1))=1/2n
n
Π
k=1
tan(kπ/(2n+1))=2n+1

といったような式になる。

これらは何を意味しているかというと、
単位円の円弧を(2n+1)等分して出来る点、y>0の範囲において、
cosはx座標、sinはy座標、のそれぞれの積である。

また、(0,1)から各点への線分の積なら
 

n-1
Π
k=1
(sin(2kπ/n))2+(1-cos(2kπ/n))2=n


のようなことも可能である。

では、今回の問題のような歯抜け状態のもので、厳密解が得られるものはどういったケースなのだろうか。

どうやら、ガロア群やガロア理論が関係しているようです。

例えば、今回の分母が13のケースでは、
13を法として、平方剰余と平方非剰余を考えます。

12≡1 (mod 13)
22≡4 (mod 13)
32≡9 (mod 13)
42≡3 (mod 13)
52≡12 (mod 13)
62≡10 (mod 13)
となって、
平方剰余の集合は、{1, 3, 4, 9, 10, 12}
平方非剰余の集合は、{2, 5, 6, 7, 8, 11}
となり、(1)式の指数が平方剰余、(2)式の指数が平方非剰余となっていることが解ります。

つまり、出題者は意図的に係数を決めたということです。

sin(π/13)sin(3π/13)sin(4π/13)sin(9π/13)sin(10π/13)sin(12π/13)=(13-313)/128
sin(2π/13)sin(5π/13)sin(6π/13)sin(7π/13)sin(8π/13)sin(11π/13)=(13+313)/128
cos(π/13)cos(3π/13)cos(4π/13)cos(9π/13)cos(10π/13)cos(12π/13)=(-11-313)/128
cos(2π/13)cos(5π/13)cos(6π/13)cos(7π/13)cos(8π/13)cos(11π/13)=(-11+313)/128
となります。


3≦2n+1≦21を表にまとめてみると、

 

平方剰余(sin)

 3579111315171921
10.86602540370.58778525220.43388373910.34202014330.28173255680.23931566420.20791169080.18374951780.16459459020.1490422661
2  0.7818314824    0.3612416661  
3    0.75574957430.6631226582    
4 0.58778525220.97492791210.98480775300.90963199530.82298386580.74314482540.67369564360.61421271260.5633200580
5    0.9898214418   0.7357239106 
6      0.9510565162 0.8371664782 
7   0.6427876096    0.91577332660.8660254037
8       0.9957341762  
9    0.54064081740.82298386580.95105651620.99573417620.99658449300.9749279121
10     0.66312265820.8660254037   
11        0.9694002659 
12     0.2393156642    
13       0.6736956436  
14          
15       0.3612416661 0.7818314824
16       0.18374951780.47594739300.6801727377
17        0.3246994692 
18         0.4338837391
19          
20          
0.86602540370.34549150280.33071891380.21650635090.10364452460.01705739190.12103072950.00198272450.00851347440.0163558660
3/2(5-5)/87/83/811/32(13-313)/12815/32(17-417)/25619/512(105-2121)/32768


平方非剰余(sin)

 3579111315171921
1          
20.86602540380.9510565163 0.64278760970.54064081750.46472317200.4067366431 0.32469946920.2947551744
3 0.95105651630.97492791220.8660254038  0.58778525230.52643216290.47594739300.4338837391
4          
5  0.78183148250.9848077530 0.93501624270.86602540380.7980172273 0.6801727378
6  0.43388373910.86602540380.98982144190.9927088741 0.8951632914 0.7818314825
7    0.90963199540.99270887410.99452189540.9618256432  
8   0.34202014330.75574957440.93501624270.9945218954 0.96940026590.9308737486
9          
10    0.2817325568  0.96182564320.99658449300.9972037972
11     0.46472317200.74314482550.8951632914 0.9972037972
12      0.58778525230.79801722730.91577332670.9749279122
13      0.4067366431 0.83716647830.9308737486
14      0.20791169080.52643216290.73572391070.8660254038
15        0.6142127127 
16          
17         0.5633200581
18        0.1645945903 
19         0.2947551744
20         0.1490422662
0.86602540380.90450849720.33071891390.16237976320.10364452470.18606760800.00756442060.13082977540.00851347450.0012244635
3/2(5+5)/87/833/3211/32(13+313)/12815/512(17+417)/25619/512(105+2121)/134217728


平方剰余(cos)

 3.00000000005.00000000007.00000000009.000000000011.000000000013.000000000015.000000000017.000000000019.000000000021.0000000000
10.50000000000.80901699440.90096886790.93969262080.95949297360.97094181740.97814760070.98297309970.98636130340.9888308262
2  0.6234898019    0.9324722294  
3    0.65486073390.7485107482    
4 -0.8090169944-0.22252093400.17364817770.41541501300.56806474670.66913060640.73900891720.78914050940.8262387743
5    0.1423148383   0.6772815716 
6      0.3090169944 0.5469481581 
7   -0.7660444431    0.40169542470.5000000000
8       0.0922683595  
9    -0.8412535328-0.5680647467-0.3090169944-0.09226835950.08257934550.2225209340
10     -0.7485107482-0.5000000000   
11        -0.2454854871 
12     -0.9709418174    
13       -0.7390089172  
14          
15       -0.9324722294 -0.6234898019
16       -0.9829730997-0.8794737512-0.7330518718
17        -0.9458172417 
18         -0.9009688679
19          
20          
0.5000000000-0.6545084972-0.1250000000-0.1250000000-0.0312500000-0.17044260800.03125000000.0039062500-0.0019531250-0.0374319363
1/2(-3-5)/8-1/8-1/8-1/32(-11-313)/1281/321/256-1/512(-5-21)/256


平方非剰余(cos)

 3579111315171921
1          
2-0.50000000000.3090169944 0.76604444310.84125353280.88545602570.9135454576 0.94581724170.9555728058
3 -0.30901699440.22252093400.5000000000  0.80901699440.85021713570.87947375120.9009688679
4          
5  -0.6234898019-0.1736481777 0.35460488700.50000000000.6026346364 0.7330518718
6  -0.9009688679-0.5000000000-0.14231483830.1205366803 0.4457383558 0.6234898019
7    -0.4154150130-0.12053668030.10452846330.2736629901  
8   -0.9396926208-0.6548607339-0.3546048870-0.1045284633 0.24548548710.3653410244
9          
10    -0.9594929736  -0.2736629901-0.08257934550.0747300936
11     -0.8854560257-0.6691306064-0.4457383558 -0.0747300936
12      -0.8090169944-0.6026346364-0.4016954247-0.2225209340
13      -0.9135454576 -0.5469481581-0.3653410244
14      -0.9781476007-0.8502171357-0.6772815716-0.5000000000
15        -0.7891405094 
16          
17         -0.8262387743
18        -0.9863613034 
19         -0.9555728058
20         -0.9888308262
-0.5000000000-0.09549150280.1250000000-0.03125000000.0312500000-0.0014323920-0.00195312500.00390625000.0019531250-0.0000254776
-1/2(-3+5)/81/8-1/321/32(-11+313)/128-1/5121/2561/512(-5-21)/13684


のように、二重根号程度の厳密解に収まっている。


ではでは


Viewing all articles
Browse latest Browse all 5376

Trending Articles