2019/8/5 HighSchoolMath 高校数学の話題 黒木玄 (Gen Kuroki) ~ このノートでは高校の数学の教科書にあるような話題を扱い, その数学的背景について解説する. タイポや自明な誤りは自分で訂正して読むこと. 本質的な誤りがあれば
63
0
0
全文
(2) 2019/8/5. HighSchoolMath. 9 ベータ函数の応用 9.1 1/6公式 9.2 sinのべきの定積分 9.3 ガンマ函数とベータ函数の関係 9.4 ベータ函数の極限によるガンマ函数の表示とWallisの公式 9.5 Gaussの超幾何函数への一般化 9.6 Kummerの超幾何函数 10 Taylor展開 10.1 Taylorの公式の証明 10.2 Taylorの公式の剰余項の評価 10.3 有名な交代級数の例 . In [1]:. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38. using Printf const e = ℯ linspace(start, stop, length) = range(start, stop=stop, length=length) using Base64 using Plots pyplot() #gr(); ENV["PLOTS_TEST"] = "true" #clibrary(:colorcet) #clibrary(:misc) function pngplot(P...; kwargs...) sleep(0.1) pngfile = tempname() * ".png" savefig(plot(P...; kwargs...), pngfile) showimg("image/png", pngfile) end pngplot(; kwargs...) = pngplot(plot!(; kwargs...)) showimg(mime, fn; scale="") = open(fn) do f base64 = base64encode(f) if scale == "" display("text/html", """<img src="data:$mime;base64,$base64">""") else display("text/html", """<img src="data:$mime;base64,$base64" width="$scale">""") end end using SymPy #sympy.init_printing(order="lex") # default #sympy.init_printing(order="rev-lex") const latex = sympy.latex using LaTeXStrings using SpecialFunctions using QuadGK using Elliptic.Jacobi: cd, sn. 1 三角函数の加法定理 三角函数の加法定理:. cos(x + y) = cosx cosy − sinx siny, sin(x + y) = cosx siny + sinx cosy.. □. この公式を認めて使えば, 三函数に関する他の多くの公式が導かれることは知っているだろう. だからよく 三角函数の加法定理だけは覚えておいて, 他の公式はそれから導けばよい. というような教え方がされている場合がある. しかし, この教え方は数学の理解という観点からはひどく中途半端である. なぜならば, 三角函数の加法定理自体がそう難しくない結果だからである. しかもその本質は中学校レベルの幾何の問題に過ぎない. 三角函数の加法定理は簡単に導けるので, 三角函数の加法定理さえ覚える必要はない.. 1.1 三角函数の加法定理の導出は易しい localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 2/63.
(3) 2019/8/5. HighSchoolMath. 以下の図を見て欲しい. . In [2]:. 1. showimg("image/jpeg", "images/trigonometric1.jpg", scale="40%"). 縦の青の点線が黒線と重なっていることを嫌うなら次のように図を描けばよい. . In [3]:. 1. showimg("image/jpeg", "images/trigonometric2.jpg", scale="40%"). 1.2 三角函数の加法定理の導出は中学校レベル. cosθ sinθ. 三角函数 , の正式な定義のためには, まず弧度法の意味での角度を定義し, 弧度法の意味ので角度の函数としてそれらを定 義しなければいけなくなる. 角度の測り方を固定しない場合には , のように角度の測り方の不定性によって角度に定 数倍 ( 倍)の違いが生じる.. a. しかし,. cos(aθ) sin(aθ). cosθ, sinθ の代わりに, c(θ) = cos(aθ), s(θ) = sin(aθ) を使っても, 三角函数の加法定理の形は変わらない: c(x + y) = c(x)c(y) − s(x)s(y), s(x + y) = c(x)s(y) + s(x)c(y).. このことから, 三角函数の加法定理を理解するためには弧度法による角度の定義を知っている必要がないことがわかる.. cos sin. 以下の図の問題を見て欲しい. それは実質的に三角函数の加法定理を示せという内容の問題である. そのことから, , という記号 を使わずに, 三角函数の加法定理と同等のことを述べることができることがわかる. そして, その問題の解答は完全に中学校数学の範 囲内の議論で可能である. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 3/63.
(4) 2019/8/5. . In [4]:. HighSchoolMath. 1. showimg("image/jpeg", "images/trigonometric3.jpg", scale="80%"). 1.3 三角函数の加法定理は複数の方法で得られる 自力で何も知らない状態から三角函数の加法定理を証明しようとすれば, 図の描き方に複数の選択肢があることに気付く. 実際にやっ てみればわかるようにどのように図を描いても, 結果的に三角函数の加法定理が得られる. 要するに, 三角函数の加法定理の証明のた めには, 知らなければできそうもないテクニカルな議論をする必要はなく, どのように図を描いても証明できる. ああやっても証明で きるし, こうやっても証明できる. そのようなことに気付けば, 三角函数の加法定理は真に易しい結果であることを納得できるはずで ある.. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 4/63.
(5) 2019/8/5. . In [5]:. HighSchoolMath. 1. showimg("image/jpeg", "images/trigonometric4.jpg", scale="80%"). sin の倍角の公式はこれの右上の図を使うと容易に証明可能である. 右上の図で α = β のとき a = 1 となるので, sin(2α) = cosα(sinα + 1sinα) = 2cosαsinα. 1.4 三角函数の加法定理と内積の関係 三角函数の加法定理:. cos(α + β) = cosα cosβ − sinα sinβ, { cos(α − β) = cosα cosβ + sinα sinβ, sin(α + β) = sinα cosβ + cosα sinβ, { sin(α − β) = sinα cosβ − cosα sinβ.. これらを成分がを極座標された2つの2次元ベクトル. ⃗ ⃗ b = ‖b‖cosβ ⃗ a⃗ = [ ac ] = [ ‖‖aa‖cosα , b = [ d ] [ ‖b‖sinβ ⃗ ] ⃗ ] ‖sinα に適用してみよう. これらの対応する成分の積の和を計算すると cos の加法定理より, ⃗ ⃗ ab + cd = ‖a‖‖⃗ b‖(cosα cosβ + sinα sinβ) = ‖a‖‖⃗ b‖cos(β − α) ⃗ b‖⃗ 倍になることがわかった. となる. これで, 2つの2次元ベクトルの対応する成分の和は2つのベクトルのあいだの角度の cos の ‖a‖‖ 我々はこれを「内積」と呼ぶのであった. 次に. ad − bc を計算してみよう:. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 5/63.
(6) 2019/8/5. HighSchoolMath. ⃗ ⃗ − α). ad − bc = ‖a‖‖⃗ b‖(cosα sinβ − sinα cosβ) = ‖a‖‖⃗ b‖sin(β これの絶対値は2つのベクトルを辺とする平行四辺形の面積である. 我々は ad − bc を行列 a b [c d] の行列式と呼んでいるのであった. 2 × 2 の行列式は平行四辺形の面積の ±1 倍という幾何学的意味を持っている. 2つの2次元ベクト ルの外積を ad − bc で定義することもできる. このように, 三角函数の加法定理は2次元ベクトルの内積や 2 × 2 の行列式(もしくは2次元ベクトルの外積)の幾何学的意味を記述して いる公式ともみなされる.. 2 3次方程式と4次方程式の解法 高校数学レベルでの3次方程式と4次方程式の解法を解説する.. 2.1 ある3次式の因数分解から3次方程式の解法へ 高校で次の因数分解の公式を習う:. x3 + y3 + z3 − 3xyz = (x + y + z)(x2 + y2 + z2 − xy − xz − yz). 1の原始3乗根を ω と書く: ω 2 + ω + 1 = 0 , ⎯⎯3 i −1 ± √ ±2πi/3 ω=e = 2 . 以下では, ω 2 + ω + 1 = 0 とそれから導かれる ω 3 = 1, ω ≠ 1 のみを使う. 1の原始3乗根 ω を使うと上の因数分解の公式は x3 + y3 + z3 − 3xyz = (x + y + z)(x + ωy + ω2 z)(x + ω2 y + ωz) と書き直される. 最初の因数分解の公式の右辺の −xy , −xz , −yz の係数 −1 は ω 2 + ω = −1 によって再現される. さらに, p = yz, q = y3 + z 3 とおくと, 上の因数分解の公式は x3 − 3px + q = (x + y + z)(x + ωy + ω2 z)(x + ω2 y + ωz) と書き直される. この公式を使うと3次方程式の解の公式を作れる.. p q に対して, p = yz, q = y3 + z3 を満たす y,z は以下のようにして求めることができる. y3 z3 = p3 と. 任意に与えられた , より. y3 + z3 = q. λ2 − qλ + p3 = (λ − y3 )(λ − z3 ).. y,z の立場を交換すれば ⎯q⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯q⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 − 4p3⎯ 2 − 4p3⎯ q + q − √ √ 3 3 y= 2 , z= 2 が成立している. 右辺の3乗根を取れば y,z も求まる: ⎯q⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯q⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯q⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯q⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 − 4p3⎯⎯ 2 − 4p3⎯⎯ + − √ √ 3 3 y=√ , z=√ . 2 2 上の x 3 − 3px + q の因数分解の公式より, x 3 − 3px + q = 0 の解はこの y,z を使って, x = −y − z, − ωy − ω2 z, − ω2 y − ωz. ゆえに, 必要ならば. と表わされる. これは本質的に所謂Cardanoの公式(カルダノの公式)である. 問題: 以上の計算を確認し, さらに解の公式を実際に作ってみよ.. □. 解答略.. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 6/63.
(7) 2019/8/5. . In [6]:. Out[6]:. . In [7]:. Out[7]: . In [8]:. Out[8]: . In [9]:. HighSchoolMath. 1 2 3 4 5. # 1の原始3乗根 x, y, z = symbols("x y z") ω₃ = factor((-1+√Sym(-3))/2) latexstring(raw"\ds ω =", latex(ω₃)). ⎯⎯3i −1 + √ ω= 2 1 2 3 4 5 6. # 因数分解の公式の確認 ω = symbols("ω") f3_factored = (x+y+z)*(x+ω*y+ω^2*z)*(x+ω^2*y+ω*z) f3 = simplify(f3_factored(ω=>ω₃)) latexstring(latex(f3_factored), "=", latex(f3)). (x + y + z) (x + yω + zω2 )(x + yω2 + zω) = x3 − 3xyz + y3 + z3 1 2 3 4 5 6 7. pp = y*z qq = y^3+z^3 latexstring("p=", latex(pp)) |> display latexstring("q=", latex(qq)) |> display p, q = symbols("p q") latexstring(latex(x^3 - 3p*x + q), "=", latex((x^3 - 3pp*x + qq))). p = yz q = y3 + z3 −3px + q + x3 = x3 − 3xyz + y3 + z3 1 2 3 4 5 6 7 8 ▾. # (m-y^3)(m-z^3) の係数 m = symbols("m") f2_org = (m-y^3)*(m-z^3) f2 = expand((m-y^3)*(m-z^3)) latexstring(latex(f2_org), "=", latex(f2)) |> display c2 = coeffs(f2, m) latexstring(raw"\text{coefficients:}\ ", latex(c2)). (m − y3 )(m − z3 ) = m2 − my3 − mz3 + y3 z3 UndefVarError: coeffs not defined Stacktrace: [1] top-level scope at In[9]:7. In [10]:. Out[10]:. 1 2 3 4 5 6 7 ▾. # 2次方程式の解 p, q = symbols("p q") equ = m^2-q*m+p^3 sol = factor.(solve(equ, m)) latexstring(latex(equ), "=0") |> display latexstring(raw"\ds m = ", latex(sol[1]), ",", latex(sol[2])). m2 − mq + p3 = 0 ⎯−4⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 3 + q2⎯ q + √⎯−4 3 + q2⎯ q − p p √ m= , 2 2. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 7/63.
(8) 2019/8/5. In [11]:. HighSchoolMath. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17. # 3次方程式の解が得られていることの確認 ▾ y = sol[2]^(Sym(1)/3) z = p/y ▾ X = -[ y+z ω*y+ω^2*z ω^2*y+ω*z ] equ = @. X^3 - 3p*X + q res = @.(simplify((f->f(ω=>ω₃))(equ))) for i in 1:3 ▾ latexstring("x_{$i}^3-3px_{$i}+q=", latex(res[i])) |> display end ▾ latexstring(raw"\ds x_1 = ", latex(X[1])) |> display ▾ latexstring(raw"\ds x_2 = ", latex(X[2])) |> display ▾ latexstring(raw"\ds x_3 = ", latex(X[3])) |> display. x13 − 3px1 + q = 0 x23 − 3px2 + q = 0 x33 − 3px3 + q = 0 3 ⎯⎯2p 2 23 √3 ⎯q⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ + √⎯−4⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ p3 + q2⎯ √ x1 = − 3 ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯− ⎯ ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 3 2 q + −4 p + q √ √ 3 ⎯⎯2pω2 2 23 ω√3 ⎯q⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ + √⎯−4⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ p3 + q2⎯ √ x2 = − 3 ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯− ⎯ ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 3 2 q + −4 p + q √ √ 3 ⎯⎯2pω 2 23 ω2 √3 ⎯q⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ + √⎯−4⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ p3 + q2⎯ √ x3 = − 3 ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯− ⎯ ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 3 2 q + −4 + √ p q √ In [12]:. 1 # 3次方程式の解 (SymPyのsolve函数による直接計算) 2 3 x, p, q = symbols("x p q") 4 equ = x^3-3p*x+q 5 sol = solve(equ, x) 6 latexstring("x^3-3px+q=0") |> display 7 display(L"x = x_1,x_2,x_3") 8 ▾ latexstring(raw"\ds x_1 = ", latex(sol[1])) |> display 9 ▾ latexstring(raw"\ds x_2 = ", latex(sol[2])) |> display 10 ▾ latexstring(raw"\ds x_3 = ", latex(sol[3])) |> display. x3 − 3px + q = 0 x = x1 , x2 , x3. ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 3+729q2⎯ −2916 √ p √2 + 2 3p x1 = − ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ − ⎯ 3 3 27q + √−2916p3+729q2 √2 2 ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 1 − √3i 3 27q + √−2916p3+729q2 − 2 ( 2 2 )√ 2 3p x2 = − − ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 3 1 − √3i 3 27q + √−2916p3+729q2 − √ 2 2 2 2 ( ) ⎯27q ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 3+729q2⎯ 3 i −2916 p √ √ 3 1 −2 + 2 )√ 2 + 2 ( 3p = − − x3 ⎯27q ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 3 3+729q2⎯ 3 i −2916 p √ √ 3 1 − + + 2 ( 2 2 )√ 2 3 27q. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 8/63.
(9) 2019/8/5. HighSchoolMath. 2.1.1 巡回行列式 1の原始3乗根. ω を用いた因数分解の公式 x3 + y3 + z3 − 3xyz = (x + y + z)(x + ωy + ω2 z)(x + ω2 y + ωz). は行列式を用いて次のように書き直される:. この公式は. 1 の原始 n 乗根 ζ を用いた公式 || x0 x1 x2 || x x x || n−1 0 1 || ⋱ xn−1 x0 || x2 ⋱ ⋱ || | x1 x2 ⋱. || x y z || 2 || z x y || = ∏(x + ωky + ω2kz). || y z x || k=0 ⋱ ⋱ ⋱ ⋱ xn−1. xn−1 ||| ⋱ ||| n−1 (x0 + ζ k x1 + ζ 2k x2 + ⋯ + ζ (n−1)k xn−1 ). x2 || = ∏ x1 ||| k=0 x0 ||. に一般化される. この公式の証明は例えば 佐武一郎, 線型代数学, 数学選書1, 裳華房 の第II章の研究課題に書いてある.. 2.2 ある4次式の展開公式から4次方程式の解法へ. f を次のように定める: このとき,. f = (w + x + y + z)(w + x − y − z)(w − x + y − z)(w − x − y + z). p = x2 + y2 + z2 , q = xyz, r = x2 y2 + x2 z2 + y2 z2. (∗). とおくと,. f = w4 − 2pw2 + 8qw + p2 − 4r. ゆえに, もしも与えられた p,q,r に対して, 条件(∗ )を満たす x,y,z を求めることができたならば, w に関する4次方程式 f = 0 は次の ように解ける:. w = −x − y − z, − x + y + z, x − y + z, x + y − z. 与えられた p,q,r に対して条件 (∗ ) を満たす x,y,z を求めるためには, 条件 x2 + y2 + z2 = p, x2 y2 + x2 z2 + y2 z2 = r, x2 y2 z2 = q2 (∗∗) を満たす x 2 , y2 , z 2 を求め, それらの平方根を取ればよい. ただし, 平方根の取り方には ±1 倍の不定性があることに注意せよ. 条件 xyz = q より, x,y,z のうち2つが決まれば残りは一意的に決まる. 条件 (∗∗ ) を満たす x 2 , y2 , z 2 は3次方程式の解と係数の関係より, 次の3次方程式の解である: λ3 − pλ2 + rλ − q2 = 0. この3次方程式は λ = Λ + p/3 とおけば次の形になる: 2 3 Λ3 − PΛ + Q = 0, P = p3 − r, Q = − 2p27 + pr3 − q2 . この形の3次方程式は前節の結果を用いれば解ける. 以上の方法はEulerの方法と呼ばれているらしい(https://en.wikipedia.org/wiki/Quartic_function) (https://en.wikipedia.org/wiki/Quartic_function)). 問題: 以上の計算を確認し, さらに解の公式を実際に作ってみよ.. □. 解答略. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 9/63.
(10) 2019/8/5. In [13]:. HighSchoolMath. 1 2 3 4 5 6 7. # 展開結果の確認 w, x, y, z = symbols("w x y z") f4_org = (w+x+y+z)*(w+x-y-z)*(w-x+y-z)*(w-x-y+z) f4 = expand(f4_org) latexstring(latex(f4_org)) |> display latexstring("=", latex(f4)) |> display. (w − x − y + z)(w − x + y − z)(w + x − y − z)(w + x + y + z) = w4 − 2w2 x2 − 2w2 y2 − 2w2 z2 + 8wxyz + x4 − 2x2 y2 − 2x2 z2 + y4 − 2y2 z2 + z4 In [14]:. 1 # 展開結果の 0,1,2,3,4 次の係数の確認. 2 # 3次の係数は0になっている. 3 4 ▾ display(L"\text{coefficients:}\ ") 5 c4 = coeffs(f4, w). coefficients: UndefVarError: coeffs not defined Stacktrace: [1] top-level scope at In[14]:5. In [15]:. 1 ▾ p = c4[2+1]/(-2) 2 latexstring("p=", latex(p)) UndefVarError: c4 not defined Stacktrace: [1] top-level scope at In[15]:1. In [16]:. 1 ▾ q = c4[2]/8 2 latexstring("q=", latex(q)) UndefVarError: c4 not defined Stacktrace: [1] top-level scope at In[16]:1. In [17]:. 1 r = x^2*y^2 + x^2*z^2 + y^2*z^2 2 ▾ @show simplify(c4[1] - (p^2-4r)) 3 latexstring("r=", latex(r)) UndefVarError: c4 not defined Stacktrace: [1] top-level scope at show.jl:555 [2] top-level scope at In[17]:2. In [18]:. 1 2 3 4 5 6. # p,q,rを使った公式の確認 ff4 = w^4 - 2p*w^2 + 8q*w + (p^2-4r) simplify(f4 - ff4) |> display latexstring("w^4-2pw^2+8qw+(p^2-4r)") |> display latexstring("=", latex(f4)) |> display. −p2 + 2pw2 − 8qw − 2w2 x2 − 2w2 y2 − 2w2 z2 + 8wxyz + x4 + 2x2 y2 + 2x2 z2 + y4 + 2y2 z2 + z4 w4 − 2pw2 + 8qw + (p2 − 4r) = w4 − 2w2 x2 − 2w2 y2 − 2w2 z2 + 8wxyz + x4 − 2x2 y2 − 2x2 z2 + y4 − 2y2 z2 + z4 localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 10/63.
(11) 2019/8/5. In [19]:. Out[19]: In [20]:. HighSchoolMath. 1 2 3 4. m = symbols("m") equ = m^3 - p*m^2 + r*m - q^2 sol = simplify(factor(equ)) latexstring(latex(equ), "=", latex(sol)). m3 − m2 p + m (x2 y2 + x2 z2 + y2 z2 ) − q2 = m3 − m2 p + mx2 y2 + mx2 z2 + my2 z2 − q2 1 2 3 4 5 6. # 補助的な3次方程式の形の確認 p, q, (f3 = (g3 = (c3 =. r, m, M = symbols("p q r m M") m^3 - p*m^2 + r*m - q^2) |> display simplify(expand(f3(m => M+p/3)))) |> display coeffs(g3, M)). m3 − m2 p + mr − q2 2 3 M 3 − Mp3 + Mr − 2p27 + pr3 − q2 UndefVarError: coeffs not defined Stacktrace: [1] top-level scope at In[20]:6. In [21]:. 1 2. # w, p, q, r = symbols("w p q r") # factor.(solve(w^4 - 2p*w^2 + 8q*w + (p^2-4r), w)). 3 べき乗和とベルヌイ多項式 高校数学ではべき乗和. 1 + 2 + ⋯ + n = n(n 2+ 1) , + 1) , 12 + 22 + ⋯ + n2 = n(n + 1)(2n 6 2 2 13 + 23 + ⋯ + n3 = n (n4+ 1). について習う. これらの公式の背景にベルヌイ多項式が控えていることを解説したい. ベルヌイ多項式については以下のツイッターにおける以下のスレッドも参照せよ: https://twitter.com/genkuroki/status/1111938896844095488 (https://twitter.com/genkuroki/status/1111938896844095488). 3.1 Bernoulli多項式 Bernoulli多項式(ベルヌイ多項式). Bk (x) が次のTaylor展開で定義される: zezx = ∞ Bk (x) zk. ez − 1 ∑ k! k=0. これの左辺をBernoulli多項式の母函数と呼ぶ.. Bk = Bk (0) はBernoulli数と呼ばれている. zezx = z exz = ∞ Bi zi ∞ xj zj ez − 1 ez − 1 ∑ i=0 i! ∑ j=0 j! ∞ + j)! ∞ k i+j k = ∑ (i i!j! Bi xj (i z+ j)! = ∑ ∑ (kj)Bk−j xj zk!. Bernoulli数: Bernoulli多項式の定数項. i,j=0. k=0 j=0. より, Bernoulli多項式はBernoulli数によって次のように表わされることがわかる:. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 11/63.
(12) 2019/8/5. HighSchoolMath. k Bk (x) = ∑ (kj)Bk−j xj . j=0. x で偏微分すると z2 ezx = ∞ B′k (x) zk ez − 1 ∑ k=0 k!. Bernoulli多項式の微分: Benoulli多項式の定義式の両辺を. となり, これの左辺は. z2 ezx = ∞ Bm (x) zm+1 = ∞ Bk−1 (x) zk ∑ m! ∑ ez − 1 m=0 k=1 (k − 1)!. と書けるので, 上と比較して. B′k (x) = kBk−1 (x) が得られる. B0 (x) = 1 なので Bk (x) は最高次の係数が 1 である k 次の多項式になる. Bernoulli多項式の積分: 上の結果より,. d Bk+1 (x) = Bk (x) dx k + 1. であるから,. Bk+1 (b) − Bk+1 (a) . ∫a Bk (x)dx = k+1 b. これの母函数表示は次の通り:. zk = b zezx dx = ezx x=b = ezb − eza . B (x)dx k ∑ k! ∫a ez − 1 [ ez − 1 ]x=a ez − 1 k=0 ∫a 特に a = 0 , b = 1 のとき, 右辺は 1 になるので, 1 ∫0 Bk (x)dx = δk0 となることもわかる. a = x , b = x + 1 の場合には ∞ x+1 k ez(x+1) − ezx zx ∞ k zk z (y)dy B ∑ ∫x k k! = ez − 1 = e = ∑ x k! ∞. b. k=0. k=0. となるので,. x+1. k ∫x Bk (y)dy = x . 3.2 Bernoulli多項式とべき乗和の関係 べき乗和. Sk (n) が次のように定義される:. n. Sk (n) = ∑ j k = 1k + 2k + ⋯ + nk . j=1. べき乗和の母函数は次のように計算される:. Sk (n) zk = ∞ n j k zk = n ∞ j k zk ∑ ∑ ∑ k=0 k! k=0 ∑ j=1 k! j=1 ∑ k=0 k! n (n+1)z z = ∑ ejz = e ez −−1 e . ∞. j=1. Bernoulli多項式の積分に関する結果より, localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 12/63.
(13) 2019/8/5. e(n+1)z − ez ez − 1. HighSchoolMath. n+1 = ∑ ∫1 Bk (x)dx zk!k k=0 ∞ = ∑ Bk+1 (n +k 1)+ −1 Bk+1 (1) zk!k . k=0 ∞. したがって,. Sk (n) = ∫1 Bk (x)dx = Bk+1 (n +k 1)+ −1 Bk+1 (1) . 特に Sk (n) は n について k + 1 次の多項式になり, 最高次の係数は 1/(k + 1) になることがわかる. n+1. 以上では母函数表示を経由して計算したが,. Bk+1 (x) k ∫x Bk (y)dy = x , ∫ Bk (x)dx = k + 1 x+1. を使えば. n. n. j+1. Sk (n) = ∑ j k = ∑ ∫j Bk (x)dx j=1 j=1 n+1 = ∫1 Bk (x)dx = Bk+1 (n +k 1)+ −1 Bk+1 (1) と同じ公式がより平易に得られる. In [22]:. 1 ▾ # ベルヌイ多項式のリスト k = 0,1,2,…,12 2 3 B(k,x) = sympy.bernoulli(k,x) 4 SS(k,x) = (B(k+1,x+1) - B(k+1,Sym(1)))/(k+1) 5 x = symbols("x") 6 for k in 0:12 7 latexstring("B_{$k}(x)=", latex(B(k,x))) |> display 8 end. B0 (x) = 1 B1 (x) = x − 12 B2 (x) = x2 − x + 16 B3 (x) = x3 − 3x22 + x2 B4 (x) = x4 − 2x3 + x2 − 301 B5 (x) = x5 − 5x24 + 5x33 − x6 B6 (x) = x6 − 3x5 + 5x24 − x22 + 421 B7 (x) = x7 − 7x26 + 7x25 − 7x63 + x6 B8 (x) = x8 − 4x7 + 14x3 6 − 7x34 + 2x32 − 301 B9 (x) = x9 − 9x28 + 6x7 − 21x5 5 + 2x3 − 3x10 B10 (x) = x10 − 5x9 + 15x2 8 − 7x6 + 5x4 − 3x22 + 665 B11 (x) = x11 − 11x210 + 55x6 9 − 11x7 + 11x5 − 11x2 3 + 5x6 691 B12 (x) = x12 − 6x11 + 11x10 − 33x2 8 + 22x6 − 33x2 4 + 5x2 − 2730. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 13/63.
(14) 2019/8/5. In [23]:. HighSchoolMath. 1 # べき乗和の公式のリスト k = 0,1,2,…,12 2 3 x = symbols("x") 4 for k in 0:12 5 ▾ latexstring("S_{$k}(x)=", latex(factor(SS(k,x)))) |> display 6 end. S0 (x) = x S1 (x) = x(x+1) 2 x(x+1)(2x+1) S2 (x) = 6 2 S3 (x) = x2(x+1) 4 S4 (x) = x(x+1)(2x+1)30(3x2+3x−1) S5 (x) = x2(x+1)2(122x2+2x−1) S6 (x) = x(x+1)(2x+1)(423x4+6x3−3x+1) S7 (x) = x2(x+1)2(3x4+624x3−x2−4x+2) S8 (x) = x(x+1)(2x+1)(5x6+15x905+5x4−15x3−x2+9x−3) S9 (x) = x2(x+1)2(x2+x−1)(202x4+4x3−x2−3x+3) S10 (x) = x(x+1)(2x+1)(x2+x−1)(3x6+966x5+2x4−11x3+3x2+10x−5) S11 (x) = x2(x+1)2(2x8+8x7+4x6−1624x5−5x4+26x3−3x2−20x+10) S12 (x) = x(x+1)(2x+1)(105x10+525x9+525x8−1050x7−11902730x6+2310x5+1420x4−3285x3−287x2+2073x−691) 問題: k が3以上の整数のとき Bk (0) = 0 となることを示せ. 証明: Bk (x) の母函数で x = 0 とおくと, ∞ B (0) k zk + z = z + z = z ez + 1 = z ez/2 + e−z/2 . ∑ 2 ez − 1 2 2 ez − 1 2 ez/2 − e−z/2 k=0 k! これは z の偶函数なので, 左辺のべき級数の奇数次の係数はすべて消える. ゆえに k が3以上の奇数ならば Bk (0) = 0 となる. □ 問題: 以下を示せ:. Bk (1 − x) = (−1)kBk (x). (2) k が奇数のとき Bk (1/2) = 0 . (3) k が 3 以上の奇数のとき Bk (1) = Bk (0) = 0 . (4) k ≧ 2 ならば Bk (0) = Bk (1). 証明: (1) Bk (x) の母函数で x に 1 − x を代入すると, ∞ B (1 − x) k k = zez(1−x) = ze−zx z ∑ ez − 1 1 − e−z k=0 k! (−z)x ∞ Bk (x) k = (−z)e e−z − 1 = ∑ k! (−z) . (1). k=0. Bk (1 − x) = (−1)kBk (x) となることがわかる. (2) x = 1/2 のとき 1 − x = x であり, k が奇数のとき Bk (1 − x) = −Bk (x) なので Bk (1/2) = 0 となることがわかる. (3) 1つ前の問題より, k が3以上の奇数のとき Bk (0) = 0 なので Bk (1) = Bk (1 − 0) = −Bk (0) = 0 となる. (4) k ≧ 2 であるとする. k が奇数ならば Bk (0) = 0 = Bk (1) となり, k が偶数ならば Bk (0) = Bk (1 − 1) = Bk (1) となる. □ なので両辺を比較すると,. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 14/63.
(15) 2019/8/5. HighSchoolMath. 問題: 以下を示せ:. Sk (x) は x で割り切れる. (2) k ≧ 1 のとき Sk (x) は x + 1 で割り切れる. (3) k が2以上の偶数ならば Sk (x) は x(x + 1)(2x + 1) で割り切れる. (4) k が3以上の奇数ならば Sk (x) は x 2 (x + 1) 2 で割り切れる. 証明: (1) Sk (0) = 0 を示せばよいが, Sk (0) = Bk+1 (1)k +− 1Bk+1 (1) = 0. (2) k ≧ 1 のとき, 1つ前の問題の(4)より Bk+1 (0) = Bk+1 (1) が成立するので, Sk (−1) = Bk+1 (0)m −+ B1k+1 (1) = 0. ゆえに Sk (x) は k + 1 で割り切れる. (3) k が2以上の偶数のとき, 1つ前の問題の(2),(3)より Bk+1 (1/2) = 0 , Bk+1 (1) = 0 が成立するので Sk (−1/2) = Bk+1 (1/2)k +−1Bk+1 (1) = 0. ゆえに Sk (x) は 2x + 1 で割り切れる. 上の(1),(2)より, Sk (x) は x と x + 1 でも割り切れるので, x(x + 1)(2x + 1) で割り切れること (1). がわかる. (4). x+1. x. Sk (x) = ∫1 Bk (t)dt = ∫0 Bk (t + 1)dt より,. Sk′ (x) = Bk (x + 1) なので, 1つ前の問題の(3)より, k が3以上の奇数ならば Sk′ (0) = Bk (1) = 0 , Sk′ (−1) = Bk (0) = 0 となる. ゆえに Sk′ (x) は x と x + 1 で割り切れる. 上の(1),(2)より, Sk (x) は x と x + 1 で割り切れる. これらより, Sk (x) が x 2 (x + 1) 2 で割り切れることがわかる. □ x(x + 1) となることを示せ. 問題: S1 (x) が x(x + 1) で割り切れることを用いて, S1 (x) = 2 x(x + 1) となる 証明: S1 (x) は最高次の係数が 1/2 の2次の多項式になるので, それが x(x + 1) で割り切れることを使えば, S1 (x) = 2 ことがただちに導かれる. □ x(x + 1)(2x + 1) となることを示せ. 問題: S2 (x) が x(x + 1)(2x + 1) で割り切れることを用いて, S2 (x) = 6 証明: S2 (x) は最高次の係数が 1/3 の3次の多項式になるので, それが x(x + 1)(2x + 1) で割り切れることを使えば, + 1) となることがただちに導かれる. □ S2 (x) = x(x + 1)(2x 6 x2 (x + 1)2 となることを示せ. 問題: S3 (x) が x 2 (x + 1) 2 で割り切れることを用いて, S3 (x) = 4 x2 (x + 1)2 と 証明: S3 (x) は最高次の係数が 1/4 の4次の多項式になるので, それが x 2 (x + 1) 2 で割り切れることを使えば, S3 (x) = 4 なることがただちに導かれる. □ 問題: S4 (x) が x(x + 1)(2x + 1) で割り切れることを用いて, S4 (x) を求めよ. 解答例: S4 (x) は最高次の係数が 1/5 の5次の多項式になり, x(x + 1)(2x + 1) で割り切れるので, S4 (x) = 101 x(x + 1)(2x + 1)(x2 + ax + b) と書ける. S4 (1) = 1 , S4 (2) = 17 より, 3 (1 + a + b) = 1, 3(4 + 2a + b) = 17. 5 localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 15/63.
(16) 2019/8/5. HighSchoolMath. 3(3 + a) = 12 となるので, a = 1 が得られ, それを左の等式に代入すると 3(2 + b) = 5 とな. 左の等式の5倍を→の等式から引くと り, が得られる. したがって,. b = −1/3. S4 (x) = 101 x(x + 1)(2x + 1)(x2 + x − 1/3) = 301 x(x + 1)(2x + 1)(3x2 + 3x − 1). □ 問題: S5 (x) が x 2 (x + 1) 2 で割り切れることを用いて, S5 (x) を求めよ. 解答例: S5 (x) は最高次の係数が 1/6 の6次多項式になり, x 2 (x + 1) 2 で割り切れるので, S5 (x) = 16 x2 (x + 1)2 (x2 + ax + b) と書ける. S5 (1) = 1 , S5 (2) = 33 より, 2 (1 + a + b) = 1, 6(4 + 2a + b) = 33. 3 前者の9倍を後者から引くと 6(3 + a) = 24 となるので, a = 1 が得られ, それを前者に代入すると, 2(2 + b) = 3 となるので, b = −1/2 が得られる. ゆえに S5 (x) = 16 x2 (x + 1)2 (x2 + x − 1/2) = 121 x2 (x + 1)2 (2x2 + 2x − 1). □ であることがわかる.. 3.3 べき乗和の直接的な取り扱い ツイッターでぴよぴよさん (https://twitter.com/piyopiyo1229)にBernoulli多項式に頼らない直接的で平易なべき乗和の取り扱い (https://twitter.com/piyopiyo1229/status/1052534499169394688)について教わったのでその内容を以下で説明する.. Sk (n) は n に関する k + 1 次の多項式になる. 証明: k + 1 次の多項式 f(x) で f(0) = 0 , f(x) = f(x − 1) + x k を満たすものが存在することを示せば十分である. (そのような f(x) k+1 について Sk (n) = f(n) となる.) f(x) = am xm とおき, f(0) = 0, f(x) = f(x − 1) + xk を満たす am 達が一意に定まることを示そ ∑ m=0 う. f(x − 1) は k+1 k+1 k+1 f(x − 1) = ∑ ai ∑ (mi)(−1)m−ixi = ∑ (∑(−1)m−i(mi)am ) xi 定理: べき乗和. m=0 i=0. i=0 m=i. f(0) = 0, f(x) = f(x − 1) + xk は係数 am 達に関する以下の連立一次方程式に書き直される: − (k + 1)ak+1 + 1 = 0, k+1 (−1)m−i(mi)am = 0 (i = k − 1,k − 2,…,1) ∑ m=i+1 a0 = 0. 1 から順番に ak , ak−1 ,…, a1 , a0 が決まり, この連立一次方程式の解 (a0 , a1 ,…, am ) が唯一つ存在すること これより, ak+1 = k+1 がわかる. □ 1 になることがわかる. さらに, ak に関する方程式は i = k − 1 の場合から 注意: 上の証明より, Sk (n) の最高次の係数は ak+1 = k+1 と書けることから, 条件. 得られ,. −kak + (k +2 1)k ak+1 = 0. ak = 1/2 となることもわかる. □ 注意: f(x) の次数が k + 1 次より大きいならば f(x) − f(x − 1) の次数は k + 1 次以上になるので, f(x) − f(x − 1) = x k が成立する ことはありえない. このことから Sk (n) は n について k + 2 次以上の多項式になることはありえないことがわかる. □ 定理: べき乗和を表す多項式 Sk (x) は以下を満たしている. (1) Sk (x) は Sk (0) = 0 , Sk (x) = Sk (x − 1) + x k という条件で一意的に特徴付けられ, その次数は k + 1 であり, その最高次の係数は 1/(k + 1) になる. になるので,. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 16/63.
(17) 2019/8/5. HighSchoolMath. k ≧ 1 ならば Sk (−1) = 0 となり, Sk (x) は x(x + 1) で割り切れる. (3) k ≧ 1 ならば Sk (x) = (−1) k+1 Sk (−1 − x). (4) k が2以上の偶数ならば Sk (−1/2) = 0 となり, Sk (x) は x(x + 1)(2x + 1) で割り切れる. (5) Sk′ (x) = (−1) k Sk′ (−1 − x) . (6) k が3以上の奇数ならば Sk′ (0) = Sk′ (−1) = 0 となり, Sk (x) は x 2 (x + 1) 2 で割り切れる. (2). 証明: (1)はすでに証明されている.. k ≧ 1 のとき, Sk (x) = Sk (x − 1) + xk で x = 0 とおくと, 0 = Sk (−1) が得られる. そのとき, Sk (x) は x + 1 で割り切れ, Sk (0) = 0 より x でも割り切れる. (3) Sk (x) = Sk (x − 1) + x k の x に −x を代入すると, Sk (−x) = Sk (−x − 1) + (−1) k x k . これは −Sk (−x − 1) = −Sk (−1 − (x − 1)) + (−1)kxk , (−1)k+1Sk (−x − 1) = (−1)k+1Sk (−1 − (x − 1)) + xk と書き直される. k ≧ 1 のと k+1 き, Sk (−1 − x) で x = 0 とおくと Sk (−1 − 0) = Sk (−1) = 0 となる. ゆえに, (−1) Sk (−1 − x) は Sk (x) を一意的に特徴付ける条 件を満たしているので, (−1) k+1 Sk (−1 − x) = Sk (x) となることがわかる. (4) k は2以上の偶数であると仮定する. このとき, (3)より Sk (x) = −Sk (−1 − x) となる. ゆえに x = −1/2 とおくと Sk (−1/2) = 0 が 得られる. そのとき Sk (x) は 2x + 1 で割り切れ, (2)より x(x + 1) でも割り切れる. (2). (5)は(3)からただちに得られる.. k は3以上の奇数であると仮定する. このとき, (5)より Sk′ (x) = −Sk′ (−1 − x) となり, 特に Sk′ (0) = −Sk′ (−1) が得られる. k Sk′ (x) = Sk′ (x − 1) + x の両辺を微分すると Sk′ (x) = Sk′ (x − 1) + kxk−1 なので,2 特に Sk′ (0) 2= Sk′ (−1) が得らえる. それらより, Sk (0) = Sk (−1) = 0 が得られる. Sk (0) = Sk (−1) = 0 と合わせると, Sk (x) は x と (x + 1) で割り切れることがわかる. □ (6). 3.4 第2種Stirling数とべき乗和 以下の内容については https://twitter.com/genkuroki/status/1052837557732433921 (https://twitter.com/genkuroki/status/1052837557732433921) も参照せよ.. {1,2,…,n} を空でない k 個の部分集合への分割(k 分割)の個数を第2種Stirling数と呼び, {nk} と書くことにする. 集合 {1,2,…,n,n + 1} の k 分割は, {1,2,…,n} の k − 1 分割と {n + 1} で構成された分割と {1,2,…,n} の k 分割中の k 個の部分 集合のどれかに n + 1 を付け加えてできる分割のどちらかになるので, n+1 = n +k n { k } { k − 1 } {k } 0 を満たしている. この漸化式と {k} = δk0 によって第2種Stirling数は一意的に特徴付けられる. 定理: ∂ = d/dx とおく. 微分作用素 x∂ の n 乗は以下のように表わされる: n (x∂)n = ∑ {nk}xk ∂k . (∗) 集合. k=0. 証明: き. n に関する数学的帰納法. n = 1 のとき {1k} = δk1 より(∗)は成立している. n について(∗)が成立していると仮定する. このと n n (x∂)n+1 = ∑ {nk}x∂xk ∂k = ∑ {nk}(xk+1 ∂k+1 + kxk ∂k ) k=0 k=0 n+1 n n = ∑ {l − 1}xl ∂l + ∑ k{nk}xk ∂k l=1 k=0 n+1 n+1 = ∑ ({k −n 1} + k{nk}) xk ∂k = ∑ {n +k 1}xk ∂k . k=0. k=0. n が n + 1 の場合にも(∗)が成立することがわかった. □ 注意: 上の定理の公式は帰納法によらずに「項と k 分割の一対一対応」を構成することによっても証明可能である. これで. https://twitter.com/genkuroki/status/1052837562342027264 (https://twitter.com/genkuroki/status/1052837562342027264) localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 17/63.
(18) 2019/8/5. に簡単な説明がある.. HighSchoolMath. □. 系: 次が成立している:. n an = ∑ {nk}a(a − 1)⋯(a − k + 1).. k=0 a 証明: x に上の定理の公式の両辺を作用させ, さらに両辺を x a で割ればこの公式が得られる.. a(a − 1)⋯(a − r) (r + 1 個の因子の積)と a(a − 1)⋯(a − r + 1) (r 個の因子の積)ついて (a + 1)a⋯(a − r + 1) − a(a − 1)⋯(a − r) = (r + 1)a(a − 1)⋯(a − r + 1). これの両辺を a = n,n − 1,…,1,0 について足し上げて, 全体を r + 1 で割ると, (n + 1)n(n − 1)⋯(n − r + 1) = n a(a − 1)⋯(a − r + 1). ∑ r+1 a=0. したがって, 上の系の公式を書き直した. k ak = ∑ {kr}a(a − 1)⋯(a − r + 1) r=0. a = 0,1,2,…,n について足し上げると, n k k (n + 1)n(n − 1)⋯(n − r + 1) k= a . ∑ ∑ { r} r+1. を. a=0. 左辺は. r=0. k = 0 のとき 1 + S0 (n) = n + 1 に一致し, k ≧ 1 のとき Sk (n) に一致する. {k0} = δk0 であることに注意せよ. □. 注意: 第2種Stirling数は次の母函数表示を持つ:. ∞ n exp(x(et − 1)) = ∑ ∑ {nk}xk tn!n .. このことは. exp(x) に exp(tx∂) を作用させると,. n=0 k=0. exp(tx∂)exp(x) = exp(xet ). となり, 上の定理の公式より,. となることからわかる. In [24]:. 1 2 3 4 5 6 7 8 9 10. □. ∞ n exp(tx∂)exp(x) = ∑(x∂)n exp(x) tn! n=0 ∞ n ∞ n n n = ∑ ∑ {nk}xk ∂k exp(x) tn! = exp(x) ∑ ∑ {nk}xk tn! n=0 k=0 n=0 k=0. # a(a-1)…(a-r+1) の和 a = symbols("a", integer=true) n = symbols("n", integer=true, positive=true) ff(a,r) = iszero(r) ? typeof(a)(1) : prod(a-i for i in 0:r-1) for r in 0:4 s = Sum(ff(a,r), (a,0,n)) latexstring(latex(s), "=", latex(factor(doit(s)))) |> display end. UndefVarError: Sum not defined Stacktrace: [1] top-level scope at .\In[24]:8. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 18/63.
(19) 2019/8/5. In [25]:. HighSchoolMath. 1 ▾ # Σ_{a=0}^n a^k を a(a-1)…(a-r+1) の和に帰着する方法で計算 2 3 stirlingsecond(n,k) = sympy.functions.combinatorial.numbers.stirling(n, k, kind=2) 4 n = symbols("n", integer=true, positive=true) 5 SSS(k,n) = expand(sum(stirlingsecond(k,r)*ff(n+1, r+1)/(r+1) for r in 0:k)) 6 7 for k in 0:12 8 latexstring(raw"\sum_{j=0}^n", "j^{$k} =", latex(factor(SSS(k,n)))) |> display 9 end. ∑nj=0 j 0 = n + 1 ∑nj=0 j 1 = n(n+1) 2 n(n+1)(2n+1) ∑nj=0 j 2 = 6 2 ∑nj=0 j 3 = n2(n+1) 4 ∑nj=0 j 4 = n(n+1)(2n+1)30(3n2+3n−1) ∑nj=0 j 5 = n2(n+1)2(122n2+2n−1) ∑nj=0 j 6 = n(n+1)(2n+1)(423n4+6n3−3n+1) ∑nj=0 j 7 = n2(n+1)2(3n4+624n3−n2−4n+2) ∑nj=0 j 8 = n(n+1)(2n+1)(5n6+15n905+5n4−15n3−n2+9n−3) ∑nj=0 j 9 = n2(n+1)2(n2+n−1)(202n4+4n3−n2−3n+3) ∑nj=0 j 10 = n(n+1)(2n+1)(n2+n−1)(3n666+9n5+2n4−11n3+3n2+10n−5) ∑nj=0 j 11 = n2(n+1)2(2n8+8n7+4n6−1624n5−5n4+26n3−3n2−20n+10) ∑nj=0 j 12 = n(n+1)(2n+1)(105n10+525n9+525n8−1050n7−11902730n6+2310n5+1420n4−3285n3−287n2+2073n−691) 3.5 第2種Stirling数とBernoulli数の関係 Bernoulli数(正確にはBernoulli多項式. Bn (x) の x = 1 での値)を第2種Stirling数で表す公式が知られている: ∞ Bn (1) = ∑ (−1)mm!{mn ++ 11} m 1+ 1 . m=0. そしてさらに, 二項係数がPascalの三角形で計算できるのと同様に, Akiyama-Tanigawa法によって, Bernoulli数を計算できることを示 せる. そのことは, ( ) から ( ) を帰納的に. a0m m = 0,1,2,…. anm n ≧ 1 an+1,m = (m + 1)(anm − an,m+1 ). によって定めると,. ∞ an0 = ∑ (−1)mm!{mn ++ 11}a0m m=0 が成立することからわかる. a0m = 1/(m + 1) のとき an0 = Bn (1) となる. 以上の結果の証明については次の論文を参照せよ. Masanobu Kaneko. The Akiyama-Tanigawa algorithm for Bernoulli numbers. Journal of Integer Sequences, Vol. 3 (2000), Article 00.2.9. Google Scholar (https://scholar.google.com/scholar?cluster=247561690595205390). localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 19/63.
(20) 2019/8/5. In [26]:. HighSchoolMath. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36. ▾ struct AkiyamaTanigawa{T,S<:Integer} ▾ a::Array{T,2} N::S end Base.length(AT::AkiyamaTanigawa) = AT.N ▾ function AkiyamaTanigawa(a0::AbstractArray{T,1}) where T N = size(a0,1) a = zeros(eltype(a0), N, N) ▾ a[0+1,:] = a0 for n in 0:N-2 for m in 0:N-n-2 ▾ a[(n+1)+1,m+1] = (m+1) * (a[n+1, m+1] - a[n+1, (m+1)+1]) end end AkiyamaTanigawa(a, N) end function displayAT(AT::AkiyamaTanigawa) N = length(AT) ▾ s = raw"\begin{matrix}" * "\n" for n in 0:N-1 ▾ s *= prod(x * raw"&" for x in @.(latex(Sym(AT.a[n+1,1:N-n])))) s = replace(s, r"&$"=>"") s *= "& "^n * raw"\\\\" * "\n" end ▾ s *= raw"\end{matrix}" * "\n" l = latexstring(raw"\displaystyle", s) display(l) end L = 10 a0 = 1 .// collect(1:L+1) AT = AkiyamaTanigawa(a0) displayAT(AT). 1. 1 2 1 1 2 3 1 1 6 6 1 0 30 − 301 − 301 0 − 421 1 1 42 42 0 301 − 301 − 301 0 − 665 5 66. 1 3 1 4 3 20 1 20 − 1403 − 281 1 140 1 20 1 220. 1 4 1 5 2 15 2 35 − 1051 − 1054 − 1051 8 165. 1 5 1 6 5 42 5 84. 0 − 281 − 2315. 1 6 1 7 3 28 5 84 1 140 − 92429. 1 7 1 8 7 72 7 120 49 3960. n+1 an0 , an1 ,…, a12 = 1/4 a13 = 1/5 1/20. 1 8 1 9 4 45 28 495. 1 9 1 10 9 110. 1 10 1 11. 1 11. anm an,m+1 の差を取り, m + 1 倍したものが an+1,m になっている. 例え 3/20 になっている. そのように計算した結果の左端の an0 が Bn (1). 以上の第 行目が である. 隣り合った と ば, 2行目の と の差 の3倍が3行目の に一致していることを, 以下のセルの計算結果と比較するとわかる.. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 20/63.
(21) 2019/8/5. In [27]: Out[27]:. HighSchoolMath. 1 ▾ [B(n,1) for n in 0:L]. 1 1 216 0 − 301 10 420 − 301 50 66 3.6 べき乗和とHurwitzのゼータ函数の関係. s > 1 , x ≠ 0,−1,−2,… に対して, Hurwitz(フルヴィッツ)のゼータ函数 ζ(s,x) が ∞ ζ(s,x) = ∑ (x +1 k)s = x1s + (x +1 1)s + (x +1 2)s + ⋯ k=0. m = 0,1,2,… に対して ζ(−m,x) − ζ(−m,x + n) = (xmm + ⋯ + (xm + n)m + (x + nm+ 1)m + ⋯) − ((x + n)m + (x + n + 1)m + ⋯) = x + (x + 1) + ⋯ + (x + n). によって定義される. 形式的には. なので, 特に. ζ(−m,1) − ζ(−m,n + 1) = 1m + 2m + ⋯ + nm = Sm (n).. この公式はHurwitzのゼータ函数の解析接続によって論理的に正当化される. 一方, ガンマ函数の応用としてよく使われる公式. 1 = 1 ∞ e−at t s−1 dt as Γ(s) ∫0. a = x,x + 1,x + 2,… の場合をHurwitzのゼータ函数の定義式に代入して, 無限和と積分の順序を交換して, 等比級数の和の公式を. の 使うと,. 1 ∞ ∞ e−(x+k)t t s−1 dt = 1 ∞ e−xt t s−1 dt = 1 ∞ te(1−x)t t s−2 dt. ζ(s,x) = Γ(s) ∑ ∫0 Γ(s) ∫0 1 − e−t Γ(s) ∫0 et − 1 k=0. このようにして, 自然にBernoulli多項式に. 1 − x を代入したものの母函数 te(1−x)t = ∞ Bk (1 − x) t k et − 1 ∑ k! k=0. が出て来る. この結果はべき乗和を無限和に拡張して得られるHurwitzのゼータ函数の中に自然にBernoulli多項式の母函数が現われる ことを意味している. さらに, から までの積分を から までの積分と から までの積分の和に分解し, から までの積. 0 ∞ 0 1 1 ∞ 0 1 N (1−x)t k te t 分の中の Bk (1 − x) の母函数 Bk (1 − x) k! を引いて足したもので置き換え, 足した分から得らえる項を 0 か et − 1 をそれから ∑ k=0 ら 1 まで積分することによって, 次が得られる: 1 ∞ te(1−x)t t s−2 dt + 1 te(1−x)t − N Bk (1 − x) t k t s−2 dt + N Bk (1 − x) 1 . ζ(s,x) = Γ(s) ∑ k! s + k − 1 ] ∫0 ( et − 1 ∑ k! ) [∫1 et − 1 k=0. この公式の右辺は によって,. k=0. Ims > −N で意味を持ち, そこへの ζ(s,x) の解析接続を与える. m ∈ ℤ, 0 ≦ m < N のとき s → −m とすること m ζ(−m,x) = (−1) Bmm+1+ 1(1 − x) = − Bmm+1+(x)1 .. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 21/63.
(22) 2019/8/5. HighSchoolMath. ここで Bk (1 − x) = (−1) k Bk (x) となることを使った. この結果を使っても, べき乗和をBernoulli多項式で表す公式. Sm (n) = ζ(−m,1) − ζ(−m,n + 1) = Bm+1 (n +m1)+−1 Bm+1 (1). が得られる. 以上の経路でのこの公式の証明はHurwitzのゼータ函数という解析学の対象を用いた分だけ難しくなっているが, Bernoulli多項式の母 函数がどのような形で自然に現われるかがよくわかる証明になっている. べき乗和を真に理解するためにはHurwitzのゼータ函数のような解析学の対象にまで視界を広げる必要がある.. 4 平面上の点と直線の距離. a,b,c は実数であり, (a,b) ≠ (0,0) であると仮定する. 高校の数学の教科書には, xy 平面上の直線 ax + by + c = 0 と xy 平面上の点 (X,Y) の距離 d が bY + c| d = |aX +⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 √a + b2⎯ と表わされることが書いてある. これは, そこに登場する様々な量の幾何学的な意味を理解していれば自明な公式に過ぎないことを以 下で説明したい.. xyz 空間内の傾いた平面 z = ax + by + c を考えよう. この平面の傾きは (a,b) で決まっている. 定理: xyz 空間内における z = ax + by + c のグラフはベクトル (a,b) の方向が登り方向の傾いた平面になり, その方向の傾きの大き ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ (a,b) ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ さは √ a2 + b2 になる. すなわち, 単位ベクトル の分だけ (x,y) をずらすと高さが √ a2 + b2 だけ増す. ⎯ ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯ 2 2 √a + b 証明: z = ax + by + c は (x,y) を (Δx,Δy) だけずらすと, aΔx + bΔy の分だけ変化する. Cauchy-Schwartzの不等式より, ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 + b2⎯ ⎯(Δx − √⎯a⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ √ )2 + (Δy)2 ≦ aΔx + bΔy ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 + b2⎯ (Δx)2 + (Δy)2 ≦ √⎯a⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ √ (a,b) ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ が成立している. ゆえに √ (Δx) 2 + (Δy) 2 = 1 という条件のもとでの aΔx + bΔy の最大値は (Δx,Δy) = のときの ⎯a⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 + b2⎯ √ 2 + b2⎯ である. □ √⎯a⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ このように z = ax + by + c における (a,b) は平面の傾きの方向と大きさを記述している. xy 平面上の点 (X,Y) から直線 ax + by + c = 0 への距離を d と書き, 点 (X,Y) から直線 ax + by + c = 0 におろした垂線の足を点 (X0 , Y0 ) と書くことにする. このとき, 点 (X,Y) は (X0 , Y0 ) からベクトル ±(a,b) と同じ方向に距離 d の位置にある. したがって, z = ax + by + c の点 (X,Y) における値は 2 + b2⎯ d aX + bY + c = ±√⎯a⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ になる. (X,Y) が (X0 , Y0 ) から見てベクトル (a,b) と同じ方向にあれば符号は + になり, その反対側にあれば符号は − になる. これ より,. bY + c| d = |aX√+⎯a⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 + b2⎯ .. z = ax + by + c の傾き方を調べることは2つのベクトル (a,b) と (Δx,Δy) の内積 ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 + b2⎯ ⎯(Δx aΔx + bΔy = √⎯a⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ √ )2 + (Δy)2 cosθ ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ を調べることに他ならない. ここで θ はそれら2つのベクトルのなす角度である. この公式から, √ (Δx) 2 + (Δy) 2 = 1 という条件の ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ もとで aΔx + bΔy が最大になるのは cosθ = 1 のときであり, 最大値は √ a2 + b2 であることがわかる. さらに, cosθ = 1 は 2つ ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ のベクトル (a,b) と (Δx,Δy) が同じ方向を向いていることを意味するので, そのことから, √ (Δx) 2 + (Δy) 2 = 1 という条件のもと (a,b) で aΔx + bΔy が最大になるのは, (Δx,Δy) = 2 + b2⎯ のときであることもわかる. √⎯a⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 補足: 平面. 高校での授業で「点と直線の距離の公式は内積を使えば容易に導ける」と習った人がいるかもしれないが, 内積の使用は実質的に「平 面 の傾き方」を調べていることに他ならない.. z = ax + by + c. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. □. 22/63.
(23) 2019/8/5. In [28]:. HighSchoolMath. 1 2 3 4 5 6. a, b, c = 0.8, 0.5, 1.0 f(x,y) = a*x + b*y + c x = -10:0.1:5 y = -10:0.1:15 surface(x, y, abs.(f.(x',y)), colorbar=false, size=(400,300)) plot!(title="\$z = |$a x + $b y + $c|\$"). Out[28]:. 5 Jensenの不等式と相加相乗調和平均 相加相乗平均の不等式はそれより圧倒的に一般的なJensenの不等式の特別な場合になっていることを解説する. さらに, 相加相乗調和 平均の一般化になっている 乗平均についても解説する.. p. 5.1 Jensenの不等式. I は実数の区間であるとする. 例えば I = ℝ , I = (0,∞) のような場合を考える. a1 ,…, an ∈ I であるとし, p1 ,…, pn ≧ 0, p1 + ⋯ + pn = 1 と仮定する. 区間 I 上の実数値函数 f(x) に対して, 実数 E[f(x)] を対応させる函数(汎函数) E[ ] を E[f(x)] = p1 f(a1 ) + ⋯ + pn f(an ) と定めると, I 上の実数値函数 f(x),g(x) と実数 α , β に対して以下が成立している. (1) 線形性: E[αf(x) + βg(x)] = αE[f(x)] + βE[g(x)] . (2) 単調性: I 全体上で f(x) ≦ g(x) ならば E[f(x)] ≦ E[g(x)] . (3) 規格化条件: I 上の定数値函数 α に対して, E[α] = α . 区間 I 上の実数値函数 f(x) が上に凸な函数であるとは, a,b ∈ I, 0 < t < 1 ⟹ (1 − t)f(a) + tf(b) ≦ f((1 − t)a + tb) が成立することだと定める. 下に凸な函数は不等式の向きを逆転することによって定義される.. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 23/63.
(24) 2019/8/5. In [29]:. HighSchoolMath. 1 # log x は上に凸な函数 2 3 x = 0:0.01:2.0 4 a, b = 0.3, 1.5 5 f(x) = log(x) 6 t = 0:0.01:1.0 7 g(a,b,t) = (1-t)*f(a) + t*f(b) 8 h(a,b,t) = (1-t)*a + t*b 9 plot(size=(400,250), legend=:topleft, xlims=(0,2.0), ylims=(-2.0, 0.8)) 10 plot!(x, f.(x), label=L"y = \log\,x") 11 plot!(h.(a,b,t), g.(a,b,t), label="") 12 ▾ plot!([a,a], [-10.0, f(a)], label=L"x = a", ls=:dash) 13 ▾ plot!([b,b], [-10.0, f(b)], label=L"x = b", ls=:dashdot). Out[29]:. In [30]:. 1 # e^x は下に凸な函数 2 3 x = -1:0.01:2 4 a, b = -0.3, 1.5 5 f(x) = e^x 6 t = 0:0.01:1.0 7 g(a,b,t) = (1-t)*f(a) + t*f(b) 8 h(a,b,t) = (1-t)*a + t*b 9 plot(size=(400,250), legend=:topleft, xlims=(-1,2), ylims=(0,8)) 10 plot!(x, f.(x), label=L"y = e^x") 11 plot!(h.(a,b,t), g.(a,b,t), label="") 12 ▾ plot!([a,a], [-0.0, f(a)], label=L"x = a", ls=:dash) 13 ▾ plot!([b,b], [-0.0, f(b)], label=L"x = b", ls=:dashdot). Out[30]:. Jensenの不等式:. f(x) が区間 I 上の上に凸な函数ならば E[f(x)] ≦ f(E[x]). (下に凸ならば不等式の向きが逆になる.). f(x) が C 1 級の場合に限定して証明する. そのように仮定しない場合には接線の存在を微分に頼らずに直接示す必要が出て来る. f(x) は上に凸であると仮定し, μ = E[x] = p1 a1 + ⋯ + pn an とおく. E[f(x)] ≦ f(μ) を示したい. x = μ における y = f(x) の接線 を y = a(x − μ) + f(μ) と書く. f(x) が上に凸であることより, I 全体上で f(x) ≦ a(x − μ) + f(μ) が成立する. ゆえに E[ ] の性質より, E[f(x)] ≦ E[a(x − μ) + f(μ)] = a(E[x] − μ) + f(μ) = f(μ). 最初の等号で E[ ] の単調性を使い, 2つ目の等号で E[ ] の線形性と規格化条件を使い, 3つ目の等号で E[x] = μ を使った. □ 証明:. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 24/63.
(25) 2019/8/5. HighSchoolMath. n. E[ ] p(x) E[f(x)] = ∫I f(x)p(x)dx がある. これは確率密度函数 p(x) を持つ確率分布における f(x) の期待値である. E[ ] は期待値汎函数(expected value functional)と呼 ばれる. □. 注意: 以上の証明法ならば に関する数学的帰納法を使わずに, しかも の定義に直接触れずに, その基本性質だけを使って証明を できた. と同じ性質を持つものの例として, 確率密度函数 に対する. E[ ]. In [31]:. 1 # 上に凸な函数 f(x) = log x の接線 2 3 x = 0:0.01:2.0 4 μ = 0.7 5 f(x) = log(x) 6 g(μ,x) = (1/μ)*(x-μ) + f(μ) 7 plot(size=(500,350), legend=:topleft, xlims=(0,1.7), ylims=(-2.2, 0.8)) 8 plot!(x, f.(x), label=L"y = \log\,x") 9 plot!(x, g.(μ,x), label=L"y = a(x-\mu)+f(\mu)", ls=:dashdot) 10 ▾ plot!([μ, μ], [-3.0, f(μ)], label=L"x=\mu", ls=:dash) 11 plot!(size=(400,250)). Out[31]:. 5.2 相加相乗調和平均の不等式 一般の相加相乗平均の不等式は き, に対して,. a1 ,…, an > 0. p1 = p2 = ⋯ = pn = 1/n と f(x) = logx の場合にJensenの不等式からただちに得られる. そのと E[logx] = log a1 + ⋯n + log an = log(a1 ⋯ an )1/n, E[x] = a1 + ⋯n + an. なので, Jensenの不等式より,. なので,. logx が単調増加函数であることより,. この不等式の. log(a1 ⋯ an )1/n = E[logx], logE[x] = log a1 + ⋯n + an (a1 ⋯ an )1/n ≦ a1 + ⋯n + an .. ai 達をそれらの逆数で置き換えて, 全体の分子分母を交換することによって, n 1/n 1 + ⋯ + 1 ≦ (a1 ,…, an ) a1 an. も得られる.. 5.3 p乗平均. x1 ,…, xn > 0 の p 乗平均 Mp (x1 ,…, xn ) を, p ≠ 0 に対して localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 25/63.
(26) 2019/8/5. HighSchoolMath. と定め,. p = 0 に対して. と定める.. M0 は相乗平均である. そして,. 1/p n 1 p Mp (x1 ,…, xn ) = ( n ∑ xi ) i=1. M0 (x1 ,…, xn ) = (x1 ,…, xn )1/n M1 (x1 ,…, xn ) = x1 + ⋯n + xn , M−1(x1 ,…, xn ) = 1 n 1 . x1 + ⋯ + xn. M1 は加法平均で, M−1 は調和平均である. p → 0 における p 乗平均の挙動を調べよう. 1 n xip = 1 n eplog xi n∑ n∑ i=1 i=1 n = 1n ∑(1 + plog xi + O(p2 )) i=1 = 1 + plog(x1 ⋯ xn )1/n + O(p2 ) なので, log(1 + X) = X + O(X 2 ) を使うと, n log Mp (x1 ,…, xn ) = 1p log 1n ∑ xip i=1 1 = p log(1 + plog(x1 ⋯ xn )1/n + O(p2 )) = log(x1 ⋯ xn )1/n + O(p) = M0 (x1 ,…, xn ) + O(p) これより, Mp (x1 ,…, xn ) は p = 0 でも解析的であることがわかる. 以下, Mp = Mp (x1 ,…, xn ) とおき, Mp の p に関する依存性について調べたい. d log Mp = d 1 log 1 n xip dp dp p n ∑ i=1 n ∑ni=1 xip log xi = − p12 log 1n ∑ xip + 1p (1/n)(1/n)∑ n xp i=1 i i=1 1 = p2 (1/n) ∑n xp i=1 i n n n 1 × ( n ∑ xip log xip − 1n ∑ xip log 1n ∑ xip ) . i=1 i=1 i=1 f(x) = xlogx とおくと, f ′ (x) = logx + 1, f ″ (x) = 1/x > 0 なので f(x) は下に凸な函数である. Jensenの不等式を 1 E[f(x)] = n ∑ni=1 f(xip ) と下に凸な函数 f(x) = xlogx に適用すると, 1 n xip log xip = E[xlogx] ≧ E[x]logE[x] n∑ i=1 n n = 1n ∑ xip log 1n ∑ xip . i=1 i=1 d これで dp log Mp ≧ 0 であることがわかった. Mp は p について単調増加函数になる: p ≦ q ⟹ Mp (x1 ,…, xn ) ≦ Mq (x1 ,…, xn ). この不等式は相加相乗平均の不等式の大幅な一般化になっている. 例えば M−1 ≦ M0 ≦ M1 より相加相乗調和平均の不等式 なので,. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 26/63.
(27) 2019/8/5. HighSchoolMath. n 1/n x1 + ⋯ + xn 1 + ⋯ + 1 ≦ (x1 ,…, xn ) ≦ n x1 xn が得られる.. p → ∞ のとき Mp (x1 ,…, xn ) → max{x1 ,…, xn } となることを示せ. 解答例: x1 = ⋯ = xk > xk+1 ≧ ⋯ ≧ xn と仮定してよい. n p x p = kx p 1 + i x ∑ i 1 ( ∑ ( x1 ) ) 問題:. i=1. i>k. なので,. n log Mp (x1 ,…, xn ) = 1p log 1n ∑ xip i=1 1 = − p logn + 1p logk + log xi p + 1p log(1 + ∑ ( xx1i ) ). i>k であり, i > k のとき 0 < xi /x1 < 1 なので, これは p → ∞ のとき logx1 に収束し, Mp (x1 ,…, xn ) → x1 = max{x1 ,…, xn } とな ることがわかる. □ 注意: 全く同様にして, p → −∞ のとき Mp (x1 ,…, xn ) → min{x1 ,…, xn } となることを示せる. □ xp + yp 1/p なので Mp (x,y) = 1 と y = (2 − xp )1/p と同値であり, M0 (x,y) = √⎯xy⎯⎯⎯ なの プロット: p ≠ 0 のとき, Mp (x,y) = ( 2 ) で M0 (x,y) = 1 と y = 1/x は同値である. Mp (x,y) = 1 のグラフをプロットしよう. In [32]:. 1 # M_p(x,y) = 1 のプロット 2 3 f(p,x) = iszero(p) ? 1/x : (2 - x^p)^(1/p) 4 P = plot(size=(500,500)) 5 ▾ ps = [-10, -2, -1, 0, 1, 2, 10] 6 for p in ps 7 a = 2^(1/p) 8 if p > 0 9 Δx = 2^(1/p)/1000 10 x = Δx:Δx:a 11 else 12 Δx = 3/1000 13 x = a+eps():Δx:3 14 end 15 plot!(x, f.(p,x), label="p = $p", ls=:auto) 16 end 17 plot(P, xlim=(0,3), ylim=(0,3), size=(300,300)). Out[32]:. 単位 内接 多角 積 localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 大値. 27/63.
(28) 2019/8/5. HighSchoolMath. 5.4 単位円に内接する多角形の周長と面積の最大値. sinα 0 ≦ α ≦ π. 正弦函数 が において上に凸な函数であることにJensenの不等式を適用することによって, 単位円に内接する の周長と面積が最大になるのは正 角形の場合であることを示そう. 一般に上に凸な函数 についてJensenの不等式より,. n. f(x). n 角形. 1 n f(xi ) ≦ f 1 n xi . n∑ (n ∑ i=1 i=1 ) が成立している. さらに, f(x) が強い意味で上に凸ならば(f(x) のグラフに局所的に直線になっている部分が存在しなければ), 等号が 成立するための必要十分条件は x1 = ⋯ = xn となることである. θ1 < ⋯ < θn < θn+1 = θ1 + 2π と仮定し, Ai = (cos θi ,sin θi ) (i = 1,…,n) とおき, 単位円に内接する n 角形 A1 ⋯ An を考え る.. αi = θi+1 − θi とおく. もしも′ αi′ > π となる i が存在するならば, 直線 Ai Ai+1 をそれに平行な原点を通る直線で線対称変換して得ら れる直線と単位円の交点を Ai , Ai+1 とし, Ai , Ai+1 のぞれぞれを A′i , A′i+1 で置き換えて得られる単位円に内接する n 角形を考えるこ とによって, 単位円に内接する n 角形の周長と面積を真に大きくすることができる. ゆえに, 単位円に内接する n 角形で周長の面積の 最大化に興味があるならば, すべての i = 1,…,n について αi ≦ π であると仮定してよい. 以下ではそのように仮定する. ′ = θi+1 − (αi − π) = θi + π , A′i = (cos θi′ ,sin θi′ ), 補足: αi = θi+1 − θi > π のとき, θi′ = θi + (αi − π) = θi+1 − π , θi+1 ⎯ ⎯⎯⎯⎯ ⎯ ⎯⎯⎯⎯ ⎯ ⎯⎯⎯ ⎯ ⎯ ⎯ ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ ⎯ ⎯ A′i+1 = (cos θi+1′ ,sin θi+1′ ) とおくと, 線分 Ai Ai+1 と線分 A′i A′i+1 は平行で同じ長さになり, n 角形 A1 ⋯ A′i A′i+1 ⋯ An の周長と面積 は n 角形 A1 ⋯ An のそれらよりも真に大きくなる. 図を描いてみよ! □ 以上の設定のもとで, n 角形 A1 ⋯ An の周長 L と面積 S について以下が成立している. ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 線分 Ai Ai+1 の長さは 2sin(αi /2) に等しいので, n n n L = 2 ∑ sin α2i = 2n 1n ∑ sin α2i ≦ 2nsin( 1n ∑ α2i ) = 2nsin πn . i=1. i=1. i=1. sinα 0 ≦ α ≦ π/2 sinα 0 ≦ α ≦ π/2 L A1 … An n 三角形 △Ai OAi+1 の面積は (1/2)sinαi に等しいので, n n n S = 12 ∑ sin αi = n2 1n ∑ sin αi ≦ n2 sin( 1n ∑ αi ) = n2 sin 2πn . i=1 i=1 i=1 この計算中の不等号は sinα が 0 ≦ α ≦ π で上に凸であることとJensenの不等式から従う. この不等式の最右辺は単位円に内接する 正 n 角形の面積に等しい. sinα が 0 ≦ α ≦ π で強い意味で凸であることを使えば, 逆に面積 S が最大になるのは単位円に内接する n 角形 A1 … An が正 n 角形になるときであることもわかる. この計算中の不等号は が で上に凸であることとJensenの不等式から従う. この不等式の最右辺は単位円に内接す る正 角形の周長に等しい. が で強い意味で凸であることを使えば, 逆に周長 が最大になるのは単位円に内接 する 角形 が正 角形になるときであることもわかる.. n n. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 28/63.
(29) 2019/8/5. In [33]:. HighSchoolMath. 1 # α_i > π ならば周長と面積をより大きくできること 2 3 t = range(0, 2π, length=400) 4 ▾ theta = 2π * [0, 3/20, 14/20, 16/20, 18/20, 1] 5 phi = copy(theta) 6 ▾ phi[2] = theta[3] - π 7 ▾ phi[3] = theta[2] + π 8 plot(size=(300,300), aspect_ratio=1, legend=false) 9 plot!(cos.(t), sin.(t), lw=0.5, color=:black) 10 plot!(cos.(theta), sin.(theta), color=:blue) 11 plot!(cos.(phi), sin.(phi), ls=:dash, color=:red). Out[33]:. 6 三角函数の微積分 高校数学の範囲内で三角函数の微分積分学を再構成してみせる. その結果は直接楕円函数論に一般化可能である.. 6.1 高校の数学の教科書の方針 高校の数学の教科書では以下のような筋道で. sinx の導函数を求めている.. 0 < x < π/2 のとき, 「面積の大小関係」によって が得られ, 全体に. 2 をかけて, 逆数を取って, sinx をかけると, 1 > sinx x > cosx. となることから, 挟み撃ちによって. を示す. そのとき. 1 sinx < 1 x < 1 sinx 2 2 2 cosx. sinx = 1 lim x→0 x. sin(−x) = sinx であることに注意せよ. これより −x x cosx − 1 = cos2 x − 1 = sin2 x → 1 (x → 0) x2 x2 (cosx + 1) x2 (cosx + 1) 2. も得られる. そして, 三角函数の加法公式. を使って,. h → 0 のとき. cos(x + y) = cosx cosy − sinx siny, sin(x + y) = cosx siny + sinx cosy.. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 29/63.
(30) 2019/8/5. HighSchoolMath. cos(x + h) − cosx = cosx cosh − sinx sinh − cosx h h cosx(cosh − 1) sinx sinh = h h2 − h → −sinx, sin(x + h) − sinx = cosx sinh + sinx cosh − sinx h h sinx(cosh − 1) cosx sinh = h h2 + h → cosx となることを示す. これで. (cosx)′ = −sinx, (sinx)′ = cosx. であることが示された. しかし, 以上の方針は次の節の方針と比較すると, 非常に遠回りになっており, 弧度法の意味での角度の定義(単位円弧の長さで角度を 定義すること)が不明瞭になっているという問題がある. In [34]:. 1. showimg("image/jpeg", "images/Jikkyo20140125limitsinc.jpg", scale="60%"). 6.2 曲線の長さが速さの積分になることの応用 高校数学IIIの教科書には. (x(t),y(t)), a ≦ t ≦ b の軌跡の長さ(曲線の長さ) L が. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. b. 30/63.
(31) 2019/8/5. HighSchoolMath. ′ (t)2 + y′ (t)2⎯ dt L = ∫a √⎯x⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ と表せることが説明されている. t を時間変数とみなすとき, 点 (x(t),y(t)) の運動の時刻 t における速度ベクトルは (x ′ (t), y′ (t)) にな ⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ り, 速さは √ x ′ (t) 2 + y′ (t) 2 と書ける. 上の公式は曲線の長さを速さの積分で表せることを意味している. b. In [35]:. 1. showimg("image/jpeg", "images/Jikkyo20140125ArcLength1.jpg", scale="60%"). In [36]:. 1. showimg("image/jpeg", "images/Jikkyo20140125ArcLength2.jpg", scale="60%"). これを使えば(曲線の長さを上の公式で定義すれば), 三角函数の微分の導出を非常に簡潔な議論で行うことができる. そのことを以下 で説明しよう.. (x(t),y(t)) = (√⎯1⎯⎯⎯⎯⎯⎯⎯⎯ − t 2⎯,t), −1 < t < 1 は単位円 x2 + y2 = 1 の右半分の上を動き, t は単位円上の点の y 座標になる. このとき, localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 31/63.
(32) 2019/8/5. HighSchoolMath. t , y′ (t) = 1, x′ (t)2 + y′ (t)2 = 1 x′ (t) = − √⎯1⎯⎯⎯⎯⎯⎯⎯⎯ 1 − t2 − t 2⎯ なので, 単位円 x 2 + y2 = 1 の右半分上の点 (x,y) に対応する弧度法の意味での角度(単位円弧の長さで定義された角度, 左回りに正, 右回りに負とみなす) θ = F(y) は y dt θ = F(y) = ∫0 ⎯⎯⎯⎯⎯⎯⎯⎯⎯ √1 − t 2⎯ と表わされる. この公式を使えば, 高校の数学の教科書の範囲内では不明瞭だった弧度法の意味での角度が積分論を使って明瞭に定義 される.. d y f(t)dt = f(y) が成立するので, dy ∫0. 一般に. dθ = F ′(y) = ⎯⎯⎯⎯⎯⎯⎯⎯⎯ 1 ⎯. dy √1 − y2. θ = F(y) は単位円上の点の y 座標に弧度法の意味での角度を対応させる函数であり, y = sinθ の定義は弧度法の意味ので角度 θ に 単位円上の点の y 座標を対応させる函数だったので, y = sinθ の定義は θ = F(y) の逆函数である. ゆえに, 逆函数の微分によって,. が得られる.. d sinθ = dy = √⎯1⎯⎯⎯⎯⎯⎯⎯⎯ 2⎯ = √⎯1⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯ 2 θ⎯ = cosθ − y − sin dθ dθ y. 要するに, 弧度法の意味での角度を単位円上の点の 座標を用いた積分で表わせば, 単に逆函数の微分として になることがわかる. この議論は非常にシンプルである.. sinθ の導函数が cosθ. 単位円の右側に制限した議論を単位円全体に拡張する作業は読者にまかせる. 現代数学的には円を多様体とみなして議論するのがよ い. 三角函数論を完璧に理解するためには「円を多様体とみなす」というような数学的教養が必要になる. In [37]:. Out[37]:. 1 2 3 4 5 6 7 ▾. t = symbols("t") x = √(1-t^2) y = t equ = diff(x, t)^2 + diff(y, t)^2 sol = simplify(equ) latexstring("x=", latex(x), raw",\quad y=", latex(y)) |> display latexstring(raw"\ds\left(\frac{dx}{dt}\right)^2+\left(\frac{dy}{dt}\right)^2=", latex(sol)). x = √⎯1⎯⎯⎯⎯⎯⎯⎯⎯ − t 2⎯, y = t dx 2 + dy 2 = − 1 ( dt ) ( dt ) t 2 − 1 問題: 直線 y = tx と単位円 x 2 + y2 = 1 の右半分の交点は 1 , t . (x(t),y(t)) = ( ⎯⎯⎯⎯⎯⎯⎯⎯⎯ √1 + t 2⎯ √⎯1⎯⎯⎯⎯⎯⎯⎯⎯ + t 2⎯ ) 原点を通る直線の傾き a をそれに対応する弧度法の意味での角度 θ に対応させる函数 θ = G(a) が a θ = G(a) = ∫0 1 +dtt 2 と書けることを確認せよ. これと逆に角度 θ を直線の傾き a に対応させる函数が tan の定義なので, a = tanθ の定義は θ = G(a) の 逆函数である. □ 解答略.. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 32/63.
(33) 2019/8/5. In [38]:. HighSchoolMath. 1 2 3 4 5 6. ▾ ▾ ▾ ▾ ▾. t, x, y = symbols("t x y") equ = [y-t*x, x^2+y^2-1] s = solve(equ, [x,y]) latexstring(raw"\text{euqation:}\ ", latex(equ[1]), raw"=0,\ ", latex(equ[2]), "=0") |> display latexstring(raw"\text{solution:}\ ") |> display s[2]. euqation: − tx + y = 0, x2 + y2 − 1 = 0 solution: Out[38]: (sqrt(1/(t^2 + 1)), t*sqrt(1/(t^2 + 1))) In [39]:. 1 ▾ X, Y= s[2][x], s[2][y] 2 sol = simplify(diff(X,t)^2 + diff(Y,t)^2) 3 latexstring(raw"\ds X=", latex(X), raw",\quad Y=", latex(Y)) |> display 4 ▾ latexstring(raw"\ds\left(\frac{dX}{dt}\right)^2+\left(\frac{dY}{dt}\right)^2=", latex(sol)) MethodError: no method matching getindex(::Tuple{Sym,Sym}, ::Sym) Closest candidates are: getindex(::Tuple, !Matched::Int64) at tuple.jl:24 getindex(::Tuple, !Matched::Real) at tuple.jl:25 getindex(::Tuple, !Matched::AbstractUnitRange{#s72} where #s72<:Real) at range.jl:287 ... Stacktrace: [1] top-level scope at In[39]:1. 問題: 以下のセルの画像を解読して, 双曲線函数の微積分の理論について整理せよ.. □. 解答略.. localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 33/63.
(34) 2019/8/5. In [40]:. HighSchoolMath. 1. showimg("image/jpeg", "images/sin-sinh.jpg", scale="80%"). localhost:8888/notebooks/HighSchoolMath/HighSchoolMath.ipynb#. 34/63.
関連したドキュメント
前章 / 節からの流れで、計算可能な関数のもつ性質を抽象的に捉えることから始めよう。話を 単純にするために、以下では次のような型のプログラム を考える。 は部分関数 (
これはつまり十進法ではなく、一進法を用いて自然数を表記するということである。とは いえ数が大きくなると見にくくなるので、.. 0, 1,
このような情念の側面を取り扱わないことには それなりの理由がある。しかし、リードもまた
としても極少数である︒そしてこのような区分は困難で相対的かつ不明確な区分となりがちである︒したがってその
大村 その場合に、なぜ成り立たなくなったのか ということ、つまりあの図式でいうと基本的には S１ という 場
神はこのように隠れておられるので、神は隠 れていると言わない宗教はどれも正しくな
自分ではおかしいと思って も、「自分の体は汚れてい るのではないか」「ひどい ことを周りの人にしたので
