7

円の内接n角形の平均面積

34
0

nつの弦の長さの平均値

問題

x2+y2=1上にn3つの点{P0=(cosθ0,sinθ0)=(1,0),P1=(cosθ1,sinθ1),,Pn1=(cosθn1,sinθn1)}
があるとする(0=θ0<θ1<θ2<<θn1)。このときn角形P0P1Pnの面積の平均値を求める。

この記事は結果こそ簡潔なものの、煩雑な計算のみしか掲載していない。

個々の計算は筆者が深夜テンションでやったものなので、間違いの見つかる可能性が大いにあり。

立式

1inに対して
ui=θiθi1.
と置く。ただしθn=2πとし、un=2πi=1n1uiui(1in1)の関数とする。
すると平均値は
02π02πu102π(u1++un2)i=1n12sinuidun1dun2du102π02πu102π(u1++un2)dun1dun2du1.(1)

となる。覚えていたら追記する。

計算

 まず分母から計算をする。
02π02πu102π(u1++un2)dun1du1=02π02πu102π(u1++un3)(2π(u1++un2))dun2du1=02π02πu102π(u1++un4)(2π(u1++un3))22dun3du1==02π02πu102π(u1++un1(k+1))(2π(u1++un1k))kk!dun1kdu1==(2π)n1(n1)!.(2)
 次に分子を計算する。02π02πu102π(u1++un2)i=1n12sinuidun1dun2du1は総和のiに対して場合分けをする。

i=nのとき

02π02πu102π(u1++un2)sinundun1du1=02π02πu102π(u1++un2)sin(2π(u1+un1))dun1du1=02π02πu102π(u1++un3){sin(2π(u1+un2)π2)+sinπ2}dun2du1=02π02πu102π(u1++un3)sin(2π(u1+un2)π2)dun2du1+sinπ202π02πu102π(u1++un3)dun2du1
ここで02π02πu102π(u1++un3)dun2du1は式(2)と全く同じ計算をすることで得られる。

=02π02πu102π(u1++un3)sin(2π(u1+un2)π2)dun2du1+sinπ202π02πu102π(u1++un3)dun2du1=02π02πu102π(u1++un3)sin(2π(u1+un2)π2)dun2du1+(2π)n2(n2)!sinπ2=02π02πu102π(u1++un4)sin(2π(u1+un3)2π2)dun3du1+(2π)n2(n2)!sinπ2+sin2π202π02πu102π(u1++un4)dun3du1=02π02πu102π(u1++un4)sin(2π(u1+un3)2π2)dun3du1+(2π)n2(n2)!sinπ2+(2π)n3(n3)!sin2π2==02π02πu102π(u1++un1(k+1))sin(2π(u1+un1k)kπ2)dun1kdu1+(2π)n2(n2)!sinπ2++(2π)n1k(n1k)!sinkπ2==02πsin(2πu1(n2)π2)du1+k=1n2(2π)n1k(n1k)!sinkπ2=k=1n2(2π)kk!sin(n1k)π2.(3)
ただし最後の行においては総和のkの向きを変えた。

1in1の場合

02π02πu102π(u1++un2)sinuidun1dun2du1=02π02πu102π(u1++ui1)Isinuiduidui1du1
ただし
I=02π(u1++ui)02π(u1++un2)dun1dui+1
と置いた。また、i=1のときには2π(u1++un2)2πとする。このIは式(1)とほぼ同じように計算できて、
I=(2π(u1++ui))ni1(ni1)!.
となる。よって、これを代入して部分積分をしていけばよい。一つ目の積分は
02π(u1++ui1)(2π(u1++ui))ni1(ni1)!sinuidui=[(2π(u1++ui))ni1(ni1)!sin(uiπ2)]02π(u1++ui1)+02π(u1++ui1)(2π(u1++ui))ni2(ni2)!sin(uiπ2)dui=(2π(u1++ui1))ni1(ni1)!sinπ2+[(2π(u1++ui))ni2(ni2)!sin(ui2π2)]02π(u1++ui1)+02π(u1++ui1)(2π(u1++ui))ni3(ni3)!sin(ui2π2)dui==(2π(u1++ui1))ni1(ni1)!sinπ2++(2π(u1++ui1))nik(nik)!sinkπ2+02π(u1++ui1)(2π(u1++ui))ni(k+1)(ni(k+1))!sin(uikπ2)dui=k=1ni1(2π(u1++ui1))nik(nik)!sinkπ2+02π(u1++ui1)sin(ui(ni1)π2)dui=k=1ni(2π(u1++ui1))nik(nik)!sinkπ2+sin(2π(u1++ui1)(ni)π2).
したがって、元の積分は
02π02πu102π(u1++un2)sinuidun1dun2du1=02π02πu102π(u1++ui2)k=1ni(2π(u1++ui1))nik(nik)!sinkπ2dui1dui2du1+02π02πu102π(u1++ui2)sin(2π(u1++ui1)(ni)π2)dui1dui2du1.
となる。ただしi=1では積分がないので、
02π02πu102π(u1++un2)sinu1dun1dun2du1=k=1n1(2π)n1k(n1k)!sinkπ2+sin(2π(n1)π2)=k=1n2(2π)kk!sin(n1k)π2.(4)
となる。ここでi=1の最後の変形は第二項を第一項のk=n1と打ち消し、そのあとkの順を逆にした。
 i1の場合の積分を続ける。(2)式とほぼ同様の計算により、
02π02πu102π(u1++ui2)k=1ni(2π(u1++ui1))nik(nik)!sinkπ2dui1dui2du1=k=1nisinkπ202π02πu102π(u1++ui2)(2π(u1++ui1))nik(nik)!dui1dui2du1=k=1ni(2π)n1k(n1k)!sinkπ2=k=i1n2(2π)kk!sin(n1k)π2.(5)
と計算される。また、(3)の計算とほぼ同じようにして(sinの位相に定数加えるだけだが、定積分で位相がマイナスされていることに注意)
02π02πu102π(u1++ui2)sin(2π(u1++ui1)(ni)π2)dui1dui2du1=k=1i2(2π)kk!sin(π2(ik1)+π2(ni))=k=1i2(2π)kk!sin(n1k)π2.(6)
と計算される。ただしi=2ならば0とする。

分子の値まとめ

 (3),(4),(5),(6)を用いると、分子の値は以下のようになる。
12[(3)+i=2n1{(4)+(5)}+(6)].
具体的値を求める。その際z(k):=(2π)kk!sin(n1k)π2とする。また、総和の範囲がi=2の時に空となるが、その時の和は0と定義する(このほうがきれいだからという理由)。具体的値は次のように求まる。
12[k=1n2z(k)+i=2n1{k=i1n2z(k)+k=1i2z(k)}+k=1n2z(k)]=12{2k=1n2z(k)+i=2n1k=1n2z(k)}=n2k=1n2z(k).
あまりにも簡潔すぎて、今までの計算の冗長性を疑ってしまう(実際漸化式を立てたほうが楽だったと今更気づいた)。

答え

 (1)にすべて代入すると
n2k=1n2(2π)kk!sin(n1k)π2(2π)n1/(n1)!=πn!(2π)nk=1n2(2π)kk!sin(n1k)π2.
となる。覚えていたらもうちょっと追記します。

投稿日:20201026
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

epidemic
epidemic
42
3010
ネタ切れ中; TeXの空白やピリオドの様式がよく分からん; 日本語記事の少ない話題を主に書く;

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. nつの弦の長さの平均値
  2. 問題
  3. 立式
  4. 計算
  5. 答え