■立方体と正八面体の断面(その9)

 今回のコラムでは

  ∫(0,∞)sin^k(ax)/x^kdx=a^(n-1)π/2-a^(n-1)π/(m-1)!2^(n-1)Σ(-1)^(r-1)(n-1,r-1)(n-2r)^(n-1)

に注目し,シンク関数の積分不等式

  1/π∫(-∞,∞)|sin(x)/x|^kdx≦√(2/k)

  1/π∫(0,∞)|sin(x)/x|^kdx≦1/√(2k)  (等号はk=2のときに限る)

を検証したい.

 すると,この証明がボールの不等式

  1≦Vn(a)≦√2

の証明につながっていくのである.

===================================

【1】シンク積分と超立方体の断面積

 シンク積分は,幾何学的には立方体の断面積として理解することができる.以下,その様子をみていくことにしたい.

 公式集

  "Table of Integrals,Series and Products"

  I.S.Gradshteyn and I.M.Ryzhik,6th Edition,Academic Press

p431-432には,

  ∫(0,∞)sin(ax)sin(bx)sin(cx)/x^ndx

とくにn=1,2,3の場合が掲げられていてる.

 それらはSinとSignがやたらとでてくる公式であって,例えば,a,b,cを正の定数0<a≦b≦cとし,

(1)a+b≦cならば,

  2/π∫(0,∞)sin(ax)sin(bx)sin(cx)/x^3dx=ab

(2)a+b>cならば,

  2/π∫(0,∞)sin(ax)sin(bx)sin(cx)/x^3dx=ab-1/4(a+b-c)^2

が成り立つ,などである.

===================================

【2】シンク関数の積分不等式

  1/π∫(0,∞)sin(x)/xdx=1/2<1/√2

  1/π∫(0,∞)sin^2(x)/x^2dx=1/2=1/2

  1/π∫(0,∞)sin^3(x)/x^3dx=3/8<1/√6

  1/π∫(0,∞)sin^4(x)/x^4dx=1/3<1/√8

  1/π∫(0,∞)sin^5(x)/x^5dx=115/384<1/√10

  1/π∫(0,∞)sin^6(x)/x^6dx=11/40<1/√12

  1/π∫(0,∞)sin^7(x)/x^7dx=5887/23040<1/√14

  1/π∫(0,∞)sin^8(x)/x^8dx=151/630<1/4

  1/π∫(0,∞)sin^9(x)/x^9dx=259732/1146880<1/√18

  1/π∫(0,∞)sin^10(x)/x^10=15619/72576<1/√20

  1/π∫(0,∞)sin^11(x)/x^11dx=381773117/1857945600<1/√22

  1/π∫(0,∞)sin^12(x)/x^12dx=655177/3326400<1/√24

  1/π∫(0,∞)sin^13(x)/x^13dx=20646903199/108999475200<1/√26

  1/π∫(0,∞)sin^14(x)/x^14dx=27085381/148262400<1/√28

  1/π∫(0,∞)sin^15(x)/x^15dx=467168310097/2645053931520<1/√30

 検証をするだけでは物足りないので,

  ∫(0,∞)sin^k(x)/x^kdx≦1/√2k

の証明のあらすじだけ紹介すると,

  sinx/x=1-x^2/6+x^4/120-・・・<1-x^2/6+x^4/120

  exp(-x^2/6)=1-x^2/6+x^4/72-x^6/1296+・・・>1-x^2/6+x^4/72-x^6/1296

  exp(-x^2/6)-sinx/x=x^4/72-x^6/1296-x^4/120>=0  (x^2≦36/5<π^2のとき)

などの不等式を利用すると証明することができる.

  0≦sinx/x≦exp(-x^2/6)  (x^2≦36/5<π^2のとき)

より

  ∫(0,∞)sin^k(x)/x^kdx≦∫(0,∞)exp(-kx^2/6)dx1/√2k=1/2√(6/πk)≦1/√2k

===================================