logo

モリファイアー 📂偏微分方程式

モリファイアー

定義1

関数ηC(Rn)\eta \in C^{\infty}(\mathbb{R}^{n})を以下のように定義しよう。

η(x):={Cexp(1x21)x<10x1 \begin{equation} \eta (x) := \begin{cases} C \exp \left( \dfrac{1}{|x|^2-1} \right) & |x|<1 \\ 0 & |x| \ge 1\end{cases} \label{1} \end{equation}

このようなη\etaモリファイヤーmollifierと言う。特にC>0C>0Rnηdx=1\displaystyle \int_{\mathbb{R}^{n}} \eta dx=1を満たす定数の時、η\etaスタンダードモリファイヤーstandard mollifierと言う。

ϵ>0\epsilon>0に対してηϵ\eta_{\epsilon}を以下のように定義しよう。

ηϵ(x):=1ϵnη(xϵ) \eta_\epsilon (x) := \dfrac{1}{\epsilon^n}\eta\left( \dfrac{x}{\epsilon} \right)

η\etaがスタンダードモリファイヤーの場合、Rnηϵdx=1\int_{\mathbb{R^n}} \eta_{\epsilon} dx=1が成り立つ。これは変数変換で簡単に確認できる。

説明

モリファイヤーは、畳み込みを通じて滑らかでない関数を滑らかにする役割を持っている。mollifyは「なだめる」という意味があり、これは、η\etaが微分不可能な関数を微分可能に近似させることをモリファイすると表現したものである。

一方で、モリファイヤーの役割を考えると、(def)\eqref{def}のような形である必要は必ずしもない。

一般化

φC(Rn)\varphi \in C^{\infty}(\mathbb{R}^{n})が以下の三つの条件を満たす場合、φ\varphiをモリファイヤーと言う。

  • コンパクトサポートを持つ。
  • Rnφdx=1\displaystyle \int_{\mathbb{R^n}} \varphi dx =1
  • limϵ0φϵ=limϵ01ϵnφ(xϵ)=δ(x)\lim \limits_{\epsilon \rightarrow 0}\varphi_\epsilon = \lim \limits_{\epsilon \rightarrow 0} \dfrac{1}{\epsilon^n}\varphi \left( \dfrac{x}{\epsilon} \right) =\delta (x) この時δ(x)\delta (x)ディラックのデルタ関数である。

以下の条件を満たす場合、ポジティブモリファイヤーpositive mollifierと言う。

  • φ(x)0,xRn\varphi (x) \ge 0, \quad \forall x \in \mathbb{R^n}

以下の条件を満たす場合、シメトリックモリファイヤーsymmetric mollifierと言う。

  • μC\mu \in C^{\infty}に対して、φ(x)=μ(x)\varphi (x)=\mu (|x|)

一方で、(1)\eqref{1}で定義されたη\etaが実際にスムース関数であるかは、以下のように証明できる。

証明

まずffggを以下のように定義しよう。

f(s):={Cexp(1s1)s<10s1g(x):=x2=x12+x22++xn2,xRn \begin{align*} f(s) &:= \begin{cases} C \exp\left( \dfrac{1}{s-1} \right) & s<1 \\ 0 & s \ge 1 \end{cases} \\g(x) &:= |x|^2={x_{1}}^2+{x_{2}}^2+\cdots+{x_{n}}^2 , \quad x\in \mathbb{R}^{n} \end{align*}

そしてfCf\in C^\inftygCg\in C^\inftyであることを示し、最終的にη=fgC\eta=f \circ g \in C^\inftyであることを示そうとする。


  • Part 1. fCf \in C^\infty

    最初に、(,1)(-\infty, 1)ffの導関数を計算すると、以下のようになる。

    f(s)= Ce1s11(s1)2f(s)= Ce1s11(s1)4+Ce1s12(s1)3f(s)= Ce1s11(s1)6+Ce1s12(s1)5+Ce1s16(s1)4 \begin{align*} f^{\prime}(s) =&\ Ce ^{ \frac{1}{s-1}} \frac{-1}{(s-1)^2} \\ f^{\prime \prime}(s) =&\ Ce ^{ \frac{1}{s-1}} \frac{-1}{(s-1)^4} + Ce ^{ \frac{1}{s-1}} \frac{2}{(s-1)^3} \\ f^{\prime \prime \prime}(s) =&\ Ce ^{ \frac{1}{s-1}} \frac{-1}{(s-1)^6} +Ce ^{ \frac{1}{s-1}} \frac{2}{(s-1)^5}+Ce ^{ \frac{1}{s-1}} \frac{-6}{(s-1)^4} \end{align*}

    したがって、すべてのk1k \ge 1に対して、以下の式を満たす定数a0(k), a1(k), , a2k(k)a^{(k)}_{0},\ a^{(k)}_{1},\ \cdots,\ a^{(k)}_{2k}が存在することを数学的帰納法を用いて示そうとする。

    f(k)(s)={e1s1(a0(k)+a1(k)s1++a2k(k)(s1)2k)s<10s1 f^{(k)}(s) = \begin{cases} e^{\frac{1}{s-1}} \left( a^{(k)}_{0} + \frac{a^{(k)}_{1}}{s-1}+\cdots+\frac{a^{(k)}_{2k}}{(s-1)^{2k}} \right) & s<1 \\ 0 & s \ge 1 \end{cases}

    • Part 1-1. k=1k=1

      f(s)={e1s1C(s1)2s<10s>1 f^{\prime}(s) = \begin{cases} e^{\frac{1}{s-1} }\frac{-C}{(s-1)^2} & s<1 \\ 0 & s >1\end{cases}

      だからa0(1)=a1(1)=0a^{(1)}_{0} =a^{(1)}_{1}=0, a2(1)=Ca^{(1)}_2=-Cである。そして、次が成り立つ。

      limh0+f(1+h)f(1)h=limh0+Ce1hh=limh0+C1he1h=0 \lim \limits_{h \rightarrow 0^+} \dfrac{ f(1+h)-f(1) }{h}=\lim \limits_{h \rightarrow 0^+} \dfrac{Ce^{ \frac{1}{h} } }{h}=\lim \limits_{h \rightarrow 0^+} \dfrac{C\frac{1}{h}}{e^{-\frac{1}{h} } }=0

      limh0f(1+h)f(1)h=limh0Ce1hh=limh0C1he1h=0 \lim \limits_{h \rightarrow 0^-} \dfrac{f(1+h)-f(1) }{h} = \lim \limits_{h \rightarrow 0^-} \dfrac{Ce^{ \frac{1}{h} } }{h}=\lim \limits_{h \rightarrow 0^-} \dfrac{C\frac{1}{h}}{e^{-\frac{1}{h} } }=0

      極値00が得られるのはロピタルの定理を使えば容易に得られる。これによりf(1)=0f^{\prime}(1)=0。従って、k=1k=1の時に成立する。

    • Part 1-2. 任意のkkの時に成り立つのであれば、k+1k+1の時も成り立つ。

      任意のkkに対して成り立つと仮定する。すると、以下の式を満たす定数a0(k), a1(k), , a2k(k)a^{(k)}_{0},\ a^{(k)}_{1},\ \cdots,\ a^{(k)}_{2k}が存在する。

      f(k)(s)={e1s1(a0(k)+a1(k)s1++a2k(k)(s1)2k)s<10s1 f^{(k)}(s) = \begin{cases} e^{\frac{1}{s-1}} \left( a^{(k)}_{0} + \frac{a^{(k)}_{1}}{s-1}+\cdots+\frac{a^{(k)}_{2k}}{(s-1)^{2k}} \right) & s<1 \\ 0 & s \ge 1 \end{cases}

      (,1)(-\infty, 1)f(k+1)(s)f^{(k+1)}(s)を計算すると、以下のようになる。

      f(k+1)(s)= e1s11(s1)2(a0(k)+a1(k)s1++a2k(k)(s1)2k)+e1s1(a1(k)(s1)2++2ka2k(k)(s1)2k+1)= e1s1(a0(k+1)+a1(k+1)s1++a2k+2(k+1)(s1)2k+2) \begin{align*} f^{(k+1)}(s) =&\ e^{\frac{1}{s-1}} \dfrac{-1}{(s-1)^2}\left( a^{(k)}_{0} + \frac{a^{(k)}_{1}}{s-1}+\cdots+\frac{a^{(k)}_{2k}}{(s-1)^{2k}} \right) + e^{\frac{1}{s-1}}\left( \frac{-a^{(k)}_{1}}{(s-1)^2}+\cdots+\frac{-2k a^{(k)}_{2k}}{(s-1)^{2k+1}} \right) \\ =&\ e^{\frac{1}{s-1}} \left( a^{(k+1)}_{0} + \frac{a^{(k+1)}_{1}}{s-1}+\cdots+\frac{a^{(k+1)}_{2k+2}}{(s-1)^{2k+2}} \right) \end{align*}

      この時、それぞれのa0(k+1), , a2k+2(k+1)a^{(k+1)}_{0},\ \cdots,\ a^{(k+1)}_{2k+2}a0(k), , a2k(k)a^{(k)}_{0},\ \cdots,\ a^{(k)}_{2k}で決定される。また、(1,)(1, \infty)ではf(k+1)(s)=0f^{(k+1)}(s)=0である。そして、以下の式が満たされるので、f(k+1)=0f^{(k+1)}=0である。

      limh0+f(1+h)f(1)h= limh0+e1s1(a0(k)+a1(k)h++a2k(k)(h)2k)h= limh0+(a0(k)h+a1(k)h2++a2k(k)(h2k+1)2k)e1h=0 \begin{align*} \lim \limits_{h \rightarrow 0^+} \dfrac{ f(1+h)-f(1)}{h}=&\ \lim \limits_{h \rightarrow 0^+} \dfrac{e^{\frac{1}{s-1}} \left( a^{(k)}_{0} + \frac{a^{(k)}_{1}}{h}+\cdots+\frac{a^{(k)}_{2k}}{(h)^{2k}} \right) }{h} \\ =&\ \lim \limits_{h \rightarrow 0^+} \dfrac{ \left( \frac{a^{(k)}_{0}}{h} + \frac{a^{(k)}_{1}}{h^2}+\cdots+\frac{a^{(k)}_{2k}}{(h^{2k+1})^{2k}} \right) }{e^{-\frac{1}{h}} } =0 \end{align*}

      limh0f(1+h)f(1)h= limh0e1s1(a0(k)+a1(k)h++a2k(k)(h)2k)h= limh0(a0(k)h+a1(k)h2++a2k(k)(h2k+1)2k)e1h=0 \begin{align*} \lim \limits_{h \rightarrow 0^-} \dfrac{ f(1+h)-f(1)}{h}=&\ \lim \limits_{h \rightarrow 0^-} \dfrac{e^{\frac{1}{s-1}} \left( a^{(k)}_{0} + \frac{a^{(k)}_{1}}{h}+\cdots+\frac{a^{(k)}_{2k}}{(h)^{2k}} \right) }{h} \\ =&\ \lim \limits_{h \rightarrow 0^-} \dfrac{ \left( \frac{a^{(k)}_{0}}{h} + \frac{a^{(k)}_{1}}{h^2}+\cdots+\frac{a^{(k)}_{2k}}{(h^{2k+1})^{2k}} \right) }{e^{-\frac{1}{h}} }=0 \end{align*}

      ロピタルの定理を繰り返し使うことで、最終的に00を得られる。従って、数学的帰納法によりfC(R)f \in C^\infty(\mathbb{R})が証明される。

  • Part 2. gC(Rn)g \in C^\infty( \mathbb{R^n})

    これはggの定義によって自明である。

従って、Part 1., Part 2. により、次が成り立つ。

η=fgC(Rn) \eta = f \circ g \in C^\infty ( \mathbb{R^n} )


  1. Lawrence C. Evans, Partial Differential Equations (第2版, 2010年), p713-714 ↩︎