logo

푸아송 방정식의 기본해 📂편미분방정식

푸아송 방정식의 기본해

빌드업1

라플라스 방정식의 기본해

xRnx \in \mathbb{R}^{n}이고, x0x \ne 0에 대해 아래의 함수 Φ\Phi를 라플라스 방정식의 기본해라고 정의한다.

Φ(x):={12πlogxn=21n(n2)α(n)1xn2n3 \Phi (x) := \begin{cases} -\frac{1}{2\pi}\log |x| & n=2 \\ \frac{1}{n(n-2)\alpha (n)} \frac{1}{|x|^{n-2}} & n \ge 3 \end{cases}

xΦ(x)x \mapsto \Phi (x)와 같이 매핑하는 함수를 생각해보자. 이는 x0x \ne 0인 곳에서 하모닉이다. 원점을 00에서 yRny\in \mathbb{R}^{n}로 대칭이동했다고 하자. 그러면 함수 xΦ(xy)x \mapsto \Phi (x-y)xyx\ne y인 곳에서 하모닉이다. 이제 임의의 함수 f:RnRf : \mathbb{R}^{n} \to \mathbb{R}가 주어졌다고 하자. 그러면 다음의 함수는, ffyy에 대한 함수이므로, 여전히 변수 xx에 대해서 하모닉이다.

xΦ(xy)f(y) x \mapsto \Phi (x-y)f(y)

그러면 각각의 ykRny_{k}\in \mathbb{R}^{n}에 대해서 위의 함수가 하모닉이므로 이를 모두 더해도 여전히 하모닉이다.

xk=1NΦ(xyk)f(yk) is harmonic in Rn{y1,,yN} x \mapsto \sum _{k=1}^{N}\Phi (x-y_{k})f(y_{k})\text{ is harmonic in } \mathbb{R}^{n}\setminus \left\{ y_{1},\dots,y_{N} \right\}

여기에서 NN을 늘려가는 센스로 함수 uu를 다음과 같이 정의하자.

정의

Φ\Phi를 라플라스 방정식의 기본해라고 하자. 그러면 다음과 같이 정의되는 uu푸아송 방정식의 기본해 라고 한다.

u(x)=RnΦ(xy)f(y)dy=Φf(x)={12πR2log(xy)f(y)dy(n=2)1n(n2)α(n)Rnf(y)xyn2dy(n3) \begin{equation} \begin{aligned} u(x) &= \int_{\mathbb{R}^n} \Phi (x-y) f(y)dy = \Phi \ast f (x) \\ &= \begin{cases} \displaystyle -\dfrac{1}{2\pi} \int_{\mathbb{R}^2} \log (|x-y|) f(y) dy & (n=2) \\ \displaystyle \dfrac{1}{n(n-2)\alpha (n) }\int_{\mathbb{R}^n} \dfrac{f(y)}{|x-y|^{n-2}}dy & (n \ge 3) \end{cases} \end{aligned} \end{equation}

이때 \ast합성곱 을 의미한다.

설명

Δu=f \begin{equation} \Delta u = f \end{equation}

이제 우리는 uu푸아송 방정식 (2)(2)을 만족하기를 기대한다. ff에 적절한 조건이 주어지면, uu가 잘 정의되고 푸아송 방정식도 만족함을 알 수 있다. 그것 조건은 ff컴팩트 서포트를 가지면서, 2번 연속적으로 미분가능한 것이다.

fCc2 f \in C^{2}_{c}

  • 잘 정의됨

    fCc(Rn)f \in C_{c}(\mathbb{R}^n)이라고 하자. 그러면 다음을 만족하는 열린 볼 B(x,rx)B(x,r_{x})가 존재한다.

    suppfB(x,rx),rx>0 \text{supp}f \subset B(x, r_{x}),\quad r_{x}>0

    그러면 다음의 계산으로 uu가 잘 정의됨을 알 수 있다.

    u(x)RnΦ(xy)f(y)dy=B(x,rx)Φ(xy)f(y)dymaxfB(x,rx)Φ(xy)dy=maxfB(0,rx)Φ(y)dy< \begin{align*} \left| u(x) \right| &\le \int_{ \mathbb{R}^{n} } \left| \Phi (x-y) \right| \left| f(y) \right| dy \\ &= \int_{ B(x, r_{x}) } \left| \Phi (x-y) \right| \left| f(y) \right| dy \\ &\le \max \left| f \right| \int_{ B(x, r_{x}) } \left| \Phi (x-y) \right| dy \\ &= \max \left| f \right| \int_{ B(0, r_{x}) } \left| \Phi (y) \right| dy < \infty \end{align*}

정리

fCc2(Rn)f \in C^{2}_{c}(\mathbb{R}^{n})라고 하자. uu(1)(1)에서와 같다고 하자. 그러면 다음이 성립한다.

(i)\text{(i)} uC2(Rn)u\in C^2 (\mathbb{R}^n)

(ii)\text{(ii)} Δu=f in Rn-\Delta u=f\quad \text{ in } \mathbb{R}^n

증명

(i)\text{(i)}

고정된 xRnx \in \mathbb{R}^n이 주어져있고 0hR0 \ne h \in \mathbb{R}, i{1,,n}i\in \left\{ 1,\cdots,n\right\}라고 하자. 그러면 다음이 성립한다.

u(x+hei)u(x)h=RnΦ(y)f(x+heiy)f(xy)hdy \dfrac{u(x+he_{i})-u(x) }{h} =\int_{\mathbb{R}^n} \Phi (y)\dfrac{f(x+he_{i}-y) -f(x-y)}{h}dy

이때 ei=(0,,1,,0)e_{i}=(0,\cdots ,1,\cdots, 0)ii번째 성분만 11이고 나머지 성분은 00인 벡터이다. 그러면, ff가 미분가능하므로, 평균값 정리(MVT)에 의해 임의의 yRny \in \mathbb{R}^nh(0,h)h^{\prime} \in (0,h)에 대해서 다음이 성립한다.

f(x+heiy)f(xy)h=fxi(x+heiy) \dfrac{f(x+he_{i}-y) - f(x-y)}{h}=f_{x_{i}}(x+h^{\prime}e_{i}-y)

가정에 의해 fxiCc1(Rn)f_{x_{i}} \in C_{c}^1 (\mathbb{R}^n)이고, 컴팩트 공간에서 연속인 함수는 균등연속이므로 fxif_{x_{i}}Rn\mathbb{R}^n에서 균등 연속이다. 따라서 주어진 ϵ>0\epsilon >0에 대해서 zw<δ    fxi(z)fxi(w)<ϵ|z-w|<\delta \implies |f_{x_{i}}(z)-f_{x_{i}}(w)|<\epsilon을 만족하는 δ>0\delta >0가 존재한다. 만약 0<h<δ0<|h|<\delta이면 모든 yRny\in \mathbb{R}^n에 대해서 다음이 성립한다.

(x+heiy)(xy)=h<h<δ |(x+h^{\prime}e_{i}-y)-(x-y)|=|h^{\prime}|<|h|<\delta

따라서 다음이 성립한다.

(x+heiy)(xy)=h<h<δ    fxi(x+heiy)fxi(xy)<ϵ    supyRnf(x+heiy)f(xy)hfxi(xy)<ϵ \begin{align*} &&|(x+h^{\prime}e_{i}-y)-(x-y)|=|h^{\prime}|<|h| &< \delta \\ \implies && |f_{x_{i}}(x+h^{\prime}e_{i}-y)-f_{x_{i}}(x-y)| &< \epsilon \\ \implies && \sup \limits_{y\in \mathbb{R}^n} \left| \dfrac{ f(x+he_{i}-y)-f(x-y)}{h}-f_{x_{i}}(x-y) \right| &< \epsilon \end{align*}

그러므로 다음이 성립한다.

f(x+heiy)f(xy)hfxi(xy)ash0 \dfrac{f(x+he_{i}-y) - f(x-y)}{h} \rightrightarrows f_{x_{i}}(x-y)\quad \mathrm{as}\quad h\rightarrow 0

따라서 다음을 얻는다.

uxi(x)=limh0u(x+hei)u(x)h=RnΦ(x)fxi(xy)dy(i=1,,) \begin{align*} u_{x_{i}}(x) &= \lim \limits_{h \rightarrow 0}\dfrac{ u(x+he_{i})-u(x)}{h} \\ &= \int_{\mathbb{R}^n} \Phi (x) f_{x_{i}}(x-y)dy \quad (i=1,\cdots, ) \end{align*}

같은 방식으로 다음이 성립한다.

uxixj(x)=limh0uxi(x+hej)uxi(x)h=RnΦ(x)fxixj(xy)dy(i,j=1,,) \begin{align*} u_{x_{i}x_{j}}(x) &= \lim \limits_{h \rightarrow 0}\dfrac{ u_{x_{i}}(x+he_{j})-u_{x_{i}}(x)}{h} \\ &= \int_{\mathbb{R}^n} \Phi (x) f_{x_{i}x_{j}}(x-y)dy \quad (i,j=1,\cdots, ) \end{align*}


  1. Lawrence C. Evans, Partial Differential Equations (2nd Edition, 2010), p 22-23 ↩︎