logo

調和関数のスムージング効果 📂偏微分方程式

調和関数のスムージング効果

定理

平均値の性質

u(x)= ⁣ ⁣ ⁣ ⁣ ⁣ ⁣B(x,r)udS= ⁣ ⁣ ⁣ ⁣ ⁣ ⁣B(x,r)udy \begin{align*} u(x) = -\!\!\!\!\!\! \int_{\partial B(x,r)} udS = -\!\!\!\!\!\! \int _{B(x,r)} udy \end{align*}

uC(Ω)u \in C(\Omega)がそれぞれの開いた球 B(x,r)ΩB(x,r)\subset \Omegaで平均値の性質を満たすとしよう。すると、次が成り立つ。

uC(Ω) u \in C^{\infty}(\Omega)

説明

ハーモニックであるなら、内部が滑らかだということだ。注意点は、境界 Ω\partial \Omegaでは滑らかさや連続性が保証されないということだ。

証明

ϵ>0\epsilon>0が与えられたとしよう。uuϵ\epsilon-モリフィケーションは次のようになる。

uϵ=ηϵuC(Ω>ϵ) u^\epsilon=\eta_\epsilon *u \in C^\infty(\Omega_{>\epsilon})

この時、Ω>ϵ:={xΩ:dist(x,Ω)>ϵ}\Omega_{>\epsilon} := \left\{ x \in \Omega : \mathrm{dist}(x, \partial \Omega) > \epsilon \right\}である。そして、xΩ>ϵx \in \Omega_{>\epsilon}としよう。すると、次が成り立つ。

uϵ(x)=Ωηϵ(xy)u(y)dy=1ϵnΩη(xyϵ)u(y)dy=1ϵnB(x,ϵ)η(xyϵ)u(y)dy \begin{align*} u^\epsilon (x) &= \int_\Omega \eta_\epsilon (x-y)u(y)dy \\ &= \dfrac{1}{\epsilon^n}\int_{\Omega} \eta \left( \frac{x-y}{\epsilon} \right) u(y)dy \\ &= \dfrac{1}{\epsilon^n}\int_{B(x,\epsilon)}\eta \left( \frac{|x-y|}{\epsilon} \right) u(y) dy \end{align*} 三つ目の等号は、モリファイアη\etaの定義により球の外では値が00であるため成り立つ。表面積分と半径に関する積分を分けると、次のようになる。

1ϵn0ϵη(rϵ)(B(x,r)udS)dr \dfrac{1}{\epsilon^n} \int_{0}^\epsilon \eta \left( \frac{r}{\epsilon} \right) \left( \int _{\partial B(x,r)} udS \right) dr

平均値の性質を使うと、次のようになる。

1ϵn0ϵη(rϵ)(B(x,r)udS)dr=1ϵn0ϵη(rϵ)(nα(n)rn1nα(n)rn1B(x,r)udS)dr=1ϵn0ϵη(rϵ)nα(n)rn1 ⁣ ⁣ ⁣ ⁣ ⁣ ⁣B(x,r)u(y)dS(y)dr=1ϵn0ϵη(rϵ)nα(n)rn1u(x)dr=1ϵnu(x)0ϵη(rϵ)nα(n)rn1dr \begin{align*} &\quad \dfrac{1}{\epsilon^n} \int_{0}^\epsilon \eta \left( \frac{r}{\epsilon} \right) \left( \int _{\partial B(x,r)} udS \right) dr \\ &= \dfrac{1}{\epsilon^n} \int_{0}^\epsilon \eta \left( \frac{r}{\epsilon} \right) \left( \dfrac{n\alpha (n)r^{n-1}}{n\alpha (n)r^{n-1}}\int _{\partial B(x,r)}u dS \right) dr \\ &= \dfrac{1}{\epsilon^n} \int_{0}^\epsilon \eta \left( \frac{r}{\epsilon} \right) n\alpha (n)r^{n-1} -\!\!\!\!\!\! \int _{\partial B(x,r)}u(y) dS(y)dr \\ &= \dfrac{1}{\epsilon^n} \int_{0}^\epsilon \eta \left( \frac{r}{\epsilon} \right) n\alpha (n)r^{n-1} u(x) dr \\ &= \dfrac{1}{\epsilon^n}u(x) \int_{0}^\epsilon \eta \left( \frac{r}{\epsilon} \right) n\alpha (n)r^{n-1}dr \end{align*}

ここでnα(n)rn1n\alpha (n)r^{n-1}は半径rrの球の表面積であるから、上の積分は以下のように書き換えることができる。

1ϵnu(x)B(0,ϵ)η(rϵ)dr=u(x)B(0,ϵ)ηϵ(y)dy=u(x) \dfrac{1}{\epsilon^n}u(x) \int_{B(0,\epsilon)}\eta \left( \frac{r}{\epsilon} \right) dr=u(x) \int_{B(0,\epsilon)} \eta_\epsilon (y)dy=u(x)

最後の等号は、ηϵ\eta_\epsilonの定義により、球B(r,ϵ)B(r,\epsilon)内の積分が11であるため成り立つ。したがって、すべてのϵ\epsilonに対してu=uϵ in Ω>ϵu=u^\epsilon \ \mathrm{in}\ \Omega_{>\epsilon}であり、uϵC(Ω)u^\epsilon \in C^\infty(\Omega)であるから、uC(Ω)u \in C^\infty (\Omega)