よんログ

解析学

前提知識

  • 微分
  • 積分
  • 複素数

テイラー展開

テイラー展開 (Taylor expansion) は、関数のある一点での導関数たちの値から計算される項の無限和として関数を表したものである。

f(x)=f(a)+f(a)(xa)+12!f(x)(xa)2+13!f(a)(xa)3+=n=01n!f(n)(a)(xa)n\begin{aligned} f(x)&=f(a)+f{'}(a)(x-a)+\frac{1}{2!}f{''}(x)(x-a)^2+\frac{1}{3!}f{'''}(a)(x-a)^3+\dots\\ &=\sum_{n=0}^\infty{\frac1{n!}f^{(n)}(a)(x-a)^n} \end{aligned}

特に a=0a=0 のときのテイラー展開をマクローリン展開 (Maclaurin expansion) という。

f(x)=n=01n!f(n)(0)xnf(x)=\sum_{n=0}^\infty{\frac1{n!}f^{(n)}(0)x^n}

テイラー展開 (マクローリン展開) をすることによって、複雑な関数を多項式で表せる。

オイラーの公式

以下の関係式をオイラーの公式 (Euler's formula) という。

eiy=cosy+isinye^{iy}=\cos{y}+i\sin{y}

特に θ=π\theta=\pi のとき

eiπ+1=0e^{i\pi}+1=0

が導かれ、この関係式をオイラーの等式 (Euler's identity) という。

オイラー公式を用いれば、複素数の極座標表示は

z=r(cosθ+isinθ)=reiθz=r(\cos{\theta}+i\sin{\theta})=re^{i\theta}

と表される。

証明 (注: 厳密でない)

cosθ,sinθ\cos\theta,\sin\theta をそれぞれマクローリン展開すると

cosθ=112!θ2+14!θ416!θ6+sinθ=θ13!θ3+15!θ517!θ7+(1)\begin{aligned} \cos\theta&=1-\dfrac1{2!}\theta^2+\dfrac1{4!}\theta^4-\dfrac1{6!}\theta^6+\cdots\\ \sin\theta&=\theta-\dfrac1{3!}\theta^3+\dfrac1{5!}\theta^5-\dfrac1{7!}\theta^7+\cdots\tag{1} \end{aligned}

exe^x をマクローリン展開すると

ex=1+x+12!x2+13!x3+14!x4+15!x5+16!x6+17!x7+(2)e^x=1+x+\dfrac1{2!}x^2+\dfrac1{3!}x^3+\dfrac1{4!}x^4+\dfrac1{5!}x^5+\dfrac1{6!}x^6+\dfrac1{7!}x^7+\cdots\tag{2}

(2) に x=iθx=i\theta を形式的に代入すると

eiθ=1+iθ12!θ213!iθ3+14!θ4+15!iθ516!θ617!iθ7+=(112!θ2+14!θ416!θ6+)+i(θ13!iθ3+15!iθ517!iθ7+)=cosθ+isinθ\begin{aligned} e^{i\theta}&=1+i\theta-\dfrac1{2!}\theta^2-\dfrac1{3!}i\theta^3+\dfrac1{4!}\theta^4+\dfrac1{5!}i\theta^5-\dfrac1{6!}\theta^6-\dfrac1{7!}i\theta^7+\cdots\\ &=\left(1-\dfrac1{2!}\theta^2+\dfrac1{4!}\theta^4-\dfrac1{6!}\theta^6+\cdots\right)+i\left(\theta-\dfrac1{3!}i\theta^3+\dfrac1{5!}i\theta^5-\dfrac1{7!}i\theta^7+\cdots\right)\\ &=\cos\theta+i\sin\theta \end{aligned}

ε-δ論法

定義

関数 f(x)f(x)x=ax=a で連続であるとは

limxaf(x)=b\lim_{x\to a}f(x)=b

が成り立つことである。

これをε-δ論法で定義すると

ε>0,δ>0 s.t. xR,xa<δf(x)b<ε^\forall\varepsilon>0,^\exists\delta>0~\mathrm{s.t.}~^\forall x\in \R,|x-a|<\delta\Rightarrow|f(x)-b|<\varepsilon

となり、これは

任意の正の数 ε\varepsilon に対し、ある適当な正の数 δ\delta が存在して、δ<xa\delta<|x-a| を満たす全ての実数 xx に対し、f(x)b<ε|f(x)-b|<\varepsilon が成り立つという意味の条件である。

ガウス積分

定義

ex2dx=π\int_{-\infty}^\infty e^{-x^2}dx=\sqrt\pi
証明

I=ex2dxI=\int_{-\infty}^\infty e^{-x^2}dx とおく

I2=(ex2dx)(ex2dx)=e(x2+y2)dxdy\begin{aligned} I^2&=\left(\int_{-\infty}^\infty e^{-x^2}dx\right)\left(\int_{-\infty}^\infty e^{-x^2}dx\right)\\ &=\int_{-\infty}^\infty\int_{-\infty}^\infty e^{-(x^2+y^2)}dxdy \end{aligned}

x=rcosθ,y=rsinθx=r\cos\theta,y=r\sin\theta とおく

I2=02π0er2{π(r+dr)2×dθ2ππr2×dθ2π}=02π0er2(rdrdθ+12(dr)2dθ)=02π0er2rdrdθ=02πdθ0rer2dr=2π[12er2]0=π\begin{aligned} I^2&=\int_0^{2\pi}\int_0^\infty e^{-r^2}\{\pi(r+dr)^2\times\frac{d\theta}{2\pi}-\pi r^2\times\frac{d\theta}{2\pi}\}\\ &=\int_0^{2\pi}\int_0^\infty e^{-r^2}(rdrd\theta+\frac12(dr)^2d\theta)\\ &=\int_0^{2\pi}\int_0^\infty e^{-r^2}rdrd\theta\\ &=\int_0^{2\pi}d\theta\int_0^\infty re^{-r^2}dr\\ &=2\pi\left[-\frac12e^{-r^2}\right]_0^\infty\\ &=\pi \end{aligned}

I>0I>0 なので

I=πI=\sqrt\pi

双曲線関数

定義

sinh(x)=exex2cosh(x)=ex+ex2tanh(x)=exexex+ex\begin{aligned} \sinh(x)&=\frac{e^x-e^{-x}}{2}\\ \cosh(x)&=\frac{e^x+e^{-x}}{2}\\ \tanh(x)&=\frac{e^x-e^{-x}}{e^x+e^{-x}} \end{aligned}

性質

双曲線関数は三角関数に似た性質をもっている。

cosh2x+sinh2x=1tanhx=sinhxcoshx1tanh2x=1cosh2x(sinhx)=coshx(coshx)=sinhx(tanhx)=1cosh2xex2dx=π\begin{aligned} \cosh^2{x}+\sinh^2{x}&=1\\ \tanh{x}&=\frac{\sinh{x}}{\cosh{x}}\\ 1-\tanh^2{x}&=\frac{1}{\cosh^2{x}}\\ (\sinh{x})'&=\cosh{x}\\ (\cosh{x})'&=\sinh{x}\\ (\tanh{x})'&=\frac1{\cosh^2{x}}\\ \int_{-\infty}^\infty e^{-x^2}dx&=\sqrt\pi \end{aligned}
なぜ三角関数に性質が似るのか

オイラーの公式 eiθ=cosθ+isinθe^{i\theta}=\cos{\theta}+i\sin{\theta} を用いて、三角関数 cosθ,sinθ\cos{\theta},\sin{\theta} を表すと

cosθ=eiθ+eiθ2sinθ=eiθeiθ2i\begin{aligned} \cos{\theta}&=\frac{e^{i\theta}+e^{-i\theta}}{2}\\ \sin{\theta}&=\frac{e^{i\theta}-e^{-i\theta}}{2i} \end{aligned}

となり、双曲線関数と式の形が似ているため。

ロピタルの定理

定義

limxa\lim_{x\rightarrow a}00\dfrac00±±\dfrac{\pm\infty}{\pm\infty} の不定形であり、次の条件をすべて満たすとき

  1. limxaf(x)=limxag(x)\lim_{x\rightarrow a}f(x)=\lim_{x\rightarrow a}g(x)00 または ±\pm\infty
  2. aa を含むある開区間から aa を除くすべての点で g(x)0g'(x)\ne0
  3. 極限 limxaf(x)g(x)\lim_{x\rightarrow a}\dfrac{f'(x)}{g'(x)} が存在する

このとき、極限 limxaf(x)g(x)\lim_{x\rightarrow a}\dfrac{f(x)}{g(x)} も存在し

limxaf(x)g(x)=limxaf(x)g(x)\lim_{x\rightarrow a}\frac{f(x)}{g(x)}=\lim_{x\rightarrow a}\frac{f'(x)}{g'(x)}

が成り立つ。

例題

  1. limx0x+sinxx\lim_{x\rightarrow0}\dfrac{x+\sin{x}}x
  2. limxxex\lim_{x\rightarrow\infty}\dfrac{x}{e^x}
  3. limx0sin2xx+sinx\lim_{x\rightarrow0}\dfrac{\sin2x}{x+\sin{x}}
  4. limx0xsinxx3\lim_{x\rightarrow0}\dfrac{x-\sin{x}}{x^3}

フーリエ解析

フーリエ級数展開

ある関数 f(x)f(x) の性質が知りたい場合、より基本的な関数系 {φ0(x),φ1(x),φ2(x),}\{\varphi_0(x),\varphi_1(x),\varphi_2(x),\dots\} で級数展開する。

f(x)=n=0cnφn(x)f(x)=\sum_{n=0}^\infty{c_n\varphi_n(x)}
  1. f(x)f(x) が無限回微分可能な場合 → {1,x,x2,x3,}\{1,x,x^2,x^3,\dots\} で級数展開 (マクローリン展開)

    f(x)=n=0f(n)(0)n!xnf(x)=\sum_{n=0}^\infty\frac{f^{(n)}(0)}{n!}x^n
  2. f(x)f(x) が周期関数である場合 → {1,cosx,sinx,cos2x,sin2x,}\{1,\cos{x},\sin{x},\cos{2x},\sin{2x},\dots\} (三角関数系) で級数展開 (フーリエ級数展開)

    f(x)=a+n=1(ancosnx+bnsinnx)f(x)=a+\sum_{n=1}^\infty(a_n\cos{nx}+b_n\sin{nx})

級数展開に三角関数を使う利点は、直交性 (自身以外との内積が 0) であること。

関数の内積

f(x),g(x)f(x),g(x) を周期 2π2\pi の周期関数とすると、その内積は

ππf(x)g(x)dx\int_{-\pi}^\pi{f(x)g(x)dx}

で定義される。

三角関数系が直交性をもつことの証明
ππsinmxsinnx=ππcos(mn)xcos(m+n)x2dx={0(mn)π(m=n)\begin{aligned} &\int_{-\pi}^\pi{\sin{mx}\cdot\sin{nx}}\\ =&\int_{-\pi}^\pi\frac{\cos(m-n)x-\cos(m+n)x}2dx\\ =&\begin{cases} 0&(m\ne n)\\ \pi&(m=n) \end{cases}\\ \end{aligned}

同様に

ππcosmxcosnx={0(mn)π(m=n)ππcosmxsinnx=0\begin{aligned} \int_{-\pi}^\pi{\cos{mx}\cdot\cos{nx}}&=\begin{cases} 0\,(m\ne n)\\ \pi\,(m=n) \end{cases}\\ \int_{-\pi}^\pi{\cos{mx}\cdot\sin{nx}}&=0 \end{aligned}

また、三角関数と 11 との直交性は自明である。

フーリエ係数

a02+n=1(ancosnx+bnsinnx)an=1πππf(x)cosnxdxbn=1πππf(x)sinnxdx\begin{aligned} \frac{a_0}2+\sum_{n=1}^\infty(a_n\cos{nx}+b_n\sin{nx})\\ a_n=\frac1\pi\int_{-\pi}^\pi{f(x)\cos{nx}dx}\\ b_n=\frac1\pi\int_{-\pi}^\pi{f(x)\sin{nx}dx} \end{aligned}

このような形で表される級数を f(x)f(x)フーリエ級数といい、an,bna_n,b_nフーリエ係数という。

フーリエ係数の導出
f(x)=a+n=1(ancosnx+bnsinnx)(1)f(x)=a+\sum_{n=1}^\infty(a_n\cos{nx}+b_n\sin{nx})\tag1

(1) の両辺に cosmx\cos{mx} をかけて、[π,π][-\pi,\pi] で積分する。

ππf(x)cosmxdx=ππacosmxdx+n=1ππ(ancosnxcosmx+bnsinnxcosmx)dx=πamam=1πππf(x)cosmxdx\begin{aligned} \int_{-\pi}^\pi{f(x)\cos{mx}dx}&=\int_{-\pi}^\pi{a\cos{mx}d x}+\sum_{n=1}^\infty\int_{-\pi}^\pi{(a_n\cos{nx}\cos{mx}+b_n\sin{nx}\cos{mx})dx}\\ &=\pi a_m\\ a_m&=\frac1\pi\int_{-\pi}^\pi{f(x)\cos{mx}dx} \end{aligned}

同様に、(1) の両辺に sinmx\sin{mx} をかけて、[π,π][-\pi,\pi] で積分する。

ππf(x)sinmxdx=ππasinmxdx+n=1ππ(ancosnxsinmx+bnsinnxsinmx)dx=πbmbm=1πππf(x)sinmxdx\begin{aligned} \int_{-\pi}^\pi{f(x)\sin{mx}dx}&=\int_{-\pi}^\pi{a\sin{mx}d x}+\sum_{n=1}^\infty\int_{-\pi}^\pi{(a_n\cos{nx}\sin{mx}+b_n\sin{nx}\sin{mx})dx}\\ &=\pi b_m\\ b_m&=\frac1\pi\int_{-\pi}^\pi{f(x)\sin{mx}dx} \end{aligned}

同様に、(1) の両辺を積分すると

ππf(x)dx=ππadx+n=1ππ(ancosnx+bnsinnx)dx=2πaa=12πππf(x)dx\begin{aligned} \int_{-\pi}^\pi{f(x)dx}&=\int_{-\pi}^\pi{ad x}+\sum_{n=1}^\infty\int_{-\pi}^\pi{(a_n\cos{nx}+b_n\sin{nx})dx}\\ &=2\pi a\\ a&=\frac1{2\pi}\int_{-\pi}^\pi{f(x)dx} \end{aligned}

よって

an=1πππf(x)cosnxdxbn=1πππf(x)sinnxdxa=12πππf(x)dx=a02\begin{aligned} a_n&=\frac1\pi\int_{-\pi}^\pi{f(x)\cos{nx}dx}\\ b_n&=\frac1\pi\int_{-\pi}^\pi{f(x)\sin{nx}dx}\\ a&=\frac1{2\pi}\int_{-\pi}^\pi{f(x)dx}=\frac{a_0}2 \end{aligned}

複素フーリエ級数

n=cneinxcn=12πππf(x)einxdx\sum_{n=-\infty}^\infty{c_ne^{inx}}\\ c_n=\frac1{2\pi}\int_{-\pi}^\pi{f(x)e^{-inx}dx}

cnc_n の性質

f(x)f(x) が実関数であるとき、cn=cnc_{-n}=c_n^* が成り立つ。

複素フーリエ級数の導出
{einx=cosnx+isinnxeinx=cosnxisinnx\begin{aligned} \begin{cases} e^{inx}&=\cos{nx}+i\sin{nx}\\ e^{-inx}&=\cos{nx}-i\sin{nx} \end{cases} \end{aligned}

より

cosnx=einx+einx2sinnx=einx+inx2i\begin{aligned} \cos{nx}&=\frac{e^{inx}+e^{-inx}}2\\ \sin{nx}&=\frac{e^{inx}+-^{-inx}}{2i} \end{aligned}a02+n=1(ancosnx+bnsinnx)=a02+n=1(anibn2einx+an+ibn2einx)=n=cneinx\begin{aligned} &\frac{a_0}2+\sum_{n=1}^\infty(a_n\cos{nx}+b_n\sin{nx})\\ =&\frac{a_0}2+\sum_{n=1}^\infty(\frac{a_n-ib_n}2e^{inx}+\frac{a_n+ib_n}2e^{-inx})\\ =&\sum_{n=-\infty}^\infty{c_ne^{inx}} \end{aligned}

ライプニッツの法則

ライプニッツの法則は、f,gf,gnn 回微分可能な関数としたときの、それらの積 fgf\cdot gnn 階微分の一般化である。

定義

(fg)(n)=k=0nnCkf(nk)g(k)(f\cdot g)^{(n)}=\sum_{k=0}^n{_nC_kf^{(n-k)}g^{(k)}}

sup, inf

定義

AA を空でない実数 RR の集合とする。

max

maxA=αA の最大値は α\max{A}=\alpha\Rightarrow A~\text{の最大値は}~\alpha
  1. 任意の xAx\in A に対して xαx\le\alpha (α\alphaAA の上界)
  2. αA\alpha\in Aα\alphaAA に入っていけなければならない

sup

supA=αA の上限は α\sup{A}=\alpha\Rightarrow A~\text{の上限は}~\alpha
  1. 任意の xAx\in A に対して xαx \le\alpha (α\alphaAA の上界)
  2. α\alpha より小さい任意の実数 β\beta に対し β<x\beta\lt x なる xAx\in A が存在 → α\alpha 自体は AA に入ってなくてもよい

同様に

minA=αA の最小値は αinfA=αA の下限は α\begin{aligned} \min{A}=\alpha &\Rightarrow A~\text{の最小値は}~\alpha \\ \inf{A}=\alpha &\Rightarrow A~\text{の下限は}~\alpha \end{aligned}

性質

  1. supAAsupA=maxA\sup{A}\in A\Rightarrow\sup{A}=\max{A}
  2. A が上に有界supA が存在A\text{ が上に有界} \Rightarrow \sup{A}\text{ が存在}
  3. infAAinfA=minA\inf{A}\in A\Rightarrow\inf{A}=min{A}
  4. A が下に有界 infA が存在A\text{ が下に有界 }\Rightarrow\inf{A}\text{ が存在}

ガンマ関数

定義

Γ(s)=0xs1exp(x)dx(s>0)\Gamma(s)=\int_0^\infty{x^{s-1}\exp{(-x)}dx}\,(s>0)

性質

  1. s>1s>1 のとき

    Γ(s)=(s1)Γ(s1)\Gamma(s)=(s-1)\Gamma(s-1)
    証明
    Γ(s)=0xs1exdx=[xs1(ex)]0+0(s1)xs2exdx=(s1)Γ(s1)\begin{aligned} \Gamma(s)&=\int_0^\infty{x^{s-1}e^{-x}dx}\\ &=\left[x^{s-1}(-e^{-x})\right]_0^\infty+\int_0^\infty{(s-1)x^{s-2}e^{-x}dx}\\ &=(s-1)\Gamma(s-1) \end{aligned}
  2. ss が正の整数のとき

    Γ(s)=(s1)!\Gamma(s)=(s-1)!
    s が正の整数のとき Γ(s)=(s-1)! となる証明
    Γ(s)=(s1)Γ(s1)=(s1)(s2)Γ(s2)=(s1)(s2)(s3)1Γ(1)Γ(1)=0exdx=[ex]0=1Γ(s)=(s1)!\begin{aligned} \Gamma(s)&=(s-1)\Gamma(s-1)\\ &=(s-1)(s-2)\Gamma(s-2)\\ &=(s-1)(s-2)(s-3)\cdots1\cdot\Gamma(1)\\ \Gamma(1)&=\int_0^\infty{e^{-x}dx}\\ &=\left[-e^{-x}\right]_0^\infty\\ &=1\\ \therefore \Gamma(s)&=(s-1)! \end{aligned}
  3. Γ(12)=π\Gamma(\frac12)=\sqrt\pi
    Γ(1/2)=√π となる証明
    Γ(12)=0x12exdx\Gamma(\frac12)=\int_0^\infty{x^{-\frac12}e^{-x}dx}

    ここで x=t2x=t^2 とおくと

    dx=2tdtx0t0dx=2tdt\\ \begin{array}{c|c} x&0\to\infty\\ \hline t&0\to\infty \end{array}

    より

    Γ(12)=0t1et22tdt=20et2dt=et2dt\begin{aligned} \Gamma(\frac12)&=\int_0^\infty{t^{-1}e^{-t^2}\cdot2tdt}\\ &=2\int_0^\infty{e^{-t^2}dt}\\ &=\int_{-\infty}^\infty{e^{-t^2}dt} \end{aligned}

    となり、ガウス積分より

    Γ(12)=π\Gamma(\frac12)=\sqrt\pi

ラグランジュの未定乗数法

条件 g(x,y)=0g(x,y)=0 下での f(x,y)f(x,y) の極値を考える。

f(x,y),g(x,y)f(x,y),g(x,y)C1C^1 級であるとき、条件付き極値をとる点の候補は

  • g(x,y)=0g(x,y)=0\nabla{g(x,y)}=\bold{0}\land g(x,y)=0
  • f(x,y)=λg(x,y)g(x,y)=0\nabla{f(x,y)}=\lambda g(x,y)\land g(x,y)=0 (ラグランジュの未定乗数法, method of Lagrange multiplier)

を満たす点である。

以上の条件は、条件付き極値となるための必要条件であり、十分条件ではない。

参考文献