푸아송 방정식의 기본해
📂편미분방정식 푸아송 방정식의 기본해 빌드업 라플라스 방정식의 기본해
x ∈ R n x \in \mathbb{R}^{n} x ∈ R n 이고, x ≠ 0 x \ne 0 x = 0 에 대해 아래의 함수 Φ \Phi Φ 를 라플라스 방정식의 기본해라고 정의한다.
Φ ( x ) : = { − 1 2 π log ∣ x ∣ n = 2 1 n ( n − 2 ) α ( n ) 1 ∣ x ∣ n − 2 n ≥ 3
\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 ) := { − 2 π 1 log ∣ x ∣ n ( n − 2 ) α ( n ) 1 ∣ x ∣ n − 2 1 n = 2 n ≥ 3
x ↦ Φ ( x ) x \mapsto \Phi (x) x ↦ Φ ( x ) 와 같이 매핑하는 함수를 생각해보자. 이는 x ≠ 0 x \ne 0 x = 0 인 곳에서 하모닉 이다. 원점을 0 0 0 에서 y ∈ R n y\in \mathbb{R}^{n} y ∈ R n 로 대칭이동했다고 하자. 그러면 함수 x ↦ Φ ( x − y ) x \mapsto \Phi (x-y) x ↦ Φ ( x − y ) 는 x ≠ y x\ne y x = y 인 곳에서 하모닉이다. 이제 임의의 함수 f : R n → R f : \mathbb{R}^{n} \to \mathbb{R} f : R n → R 가 주어졌다고 하자. 그러면 다음의 함수는, f f f 가 y y y 에 대한 함수이므로, 여전히 변수 x x x 에 대해서 하모닉이다.
x ↦ Φ ( x − y ) f ( y )
x \mapsto \Phi (x-y)f(y)
x ↦ Φ ( x − y ) f ( y )
그러면 각각의 y k ∈ R n y_{k}\in \mathbb{R}^{n} y k ∈ R n 에 대해서 위의 함수가 하모닉이므로 이를 모두 더해도 여전히 하모닉이다.
x ↦ ∑ k = 1 N Φ ( x − y k ) f ( y k ) is harmonic in R n ∖ { y 1 , … , y N }
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\}
x ↦ k = 1 ∑ N Φ ( x − y k ) f ( y k ) is harmonic in R n ∖ { y 1 , … , y N }
여기에서 N N N 을 늘려가는 센스로 함수 u u u 를 다음과 같이 정의하자.
정의 Φ \Phi Φ 를 라플라스 방정식의 기본해라고 하자. 그러면 다음과 같이 정의되는 u u u 를 푸아송 방정식의 기본해 라고 한다.
u ( x ) = ∫ R n Φ ( x − y ) f ( y ) d y = Φ ∗ f ( x ) = { − 1 2 π ∫ R 2 log ( ∣ x − y ∣ ) f ( y ) d y ( n = 2 ) 1 n ( n − 2 ) α ( n ) ∫ R n f ( y ) ∣ x − y ∣ n − 2 d y ( n ≥ 3 )
\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}
u ( x ) = ∫ R n Φ ( x − y ) f ( y ) d y = Φ ∗ f ( x ) = ⎩ ⎨ ⎧ − 2 π 1 ∫ R 2 log ( ∣ x − y ∣ ) f ( y ) d y n ( n − 2 ) α ( n ) 1 ∫ R n ∣ x − y ∣ n − 2 f ( y ) d y ( n = 2 ) ( n ≥ 3 )
이때 ∗ \ast ∗ 는 합성곱 을 의미한다.
설명 Δ u = f
\begin{equation}
\Delta u = f
\end{equation}
Δ u = f
이제 우리는 u u u 가 푸아송 방정식 ( 2 ) (2) ( 2 ) 을 만족하기를 기대한다. f f f 에 적절한 조건이 주어지면, u u u 가 잘 정의되고 푸아송 방정식도 만족함을 알 수 있다. 그것 조건은 f f f 가 컴팩트 서포트 를 가지면서, 2번 연속적으로 미분가능 한 것이다.
f ∈ C c 2
f \in C^{2}_{c}
f ∈ C c 2
잘 정의됨
f ∈ C c ( R n ) f \in C_{c}(\mathbb{R}^n) f ∈ C c ( R n ) 이라고 하자. 그러면 다음을 만족하는 열린 볼 B ( x , r x ) B(x,r_{x}) B ( x , r x ) 가 존재한다.
supp f ⊂ B ( x , r x ) , r x > 0
\text{supp}f \subset B(x, r_{x}),\quad r_{x}>0
supp f ⊂ B ( x , r x ) , r x > 0
그러면 다음의 계산으로 u u u 가 잘 정의됨을 알 수 있다.
∣ u ( x ) ∣ ≤ ∫ R n ∣ Φ ( x − y ) ∣ ∣ f ( y ) ∣ d y = ∫ B ( x , r x ) ∣ Φ ( x − y ) ∣ ∣ f ( y ) ∣ d y ≤ max ∣ f ∣ ∫ B ( x , r x ) ∣ Φ ( x − y ) ∣ d y = max ∣ f ∣ ∫ B ( 0 , r x ) ∣ Φ ( y ) ∣ d y < ∞
\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*}
∣ u ( x ) ∣ ≤ ∫ R n ∣ Φ ( x − y ) ∣ ∣ f ( y ) ∣ d y = ∫ B ( x , r x ) ∣ Φ ( x − y ) ∣ ∣ f ( y ) ∣ d y ≤ max ∣ f ∣ ∫ B ( x , r x ) ∣ Φ ( x − y ) ∣ d y = max ∣ f ∣ ∫ B ( 0 , r x ) ∣ Φ ( y ) ∣ d y < ∞
정리 f ∈ C c 2 ( R n ) f \in C^{2}_{c}(\mathbb{R}^{n}) f ∈ C c 2 ( R n ) 라고 하자. u u u 가 ( 1 ) (1) ( 1 ) 에서와 같다고 하자. 그러면 다음이 성립한다.
(i) \text{(i)} (i) u ∈ C 2 ( R n ) u\in C^2 (\mathbb{R}^n) u ∈ C 2 ( R n )
(ii) \text{(ii)} (ii) − Δ u = f in R n -\Delta u=f\quad \text{ in } \mathbb{R}^n − Δ u = f in R n
증명 (i) \text{(i)} (i) 고정된 x ∈ R n x \in \mathbb{R}^n x ∈ R n 이 주어져있고 0 ≠ h ∈ R 0 \ne h \in \mathbb{R} 0 = h ∈ R , i ∈ { 1 , ⋯ , n } i\in \left\{ 1,\cdots,n\right\} i ∈ { 1 , ⋯ , n } 라고 하자. 그러면 다음이 성립한다.
u ( x + h e i ) − u ( x ) h = ∫ R n Φ ( y ) f ( x + h e i − y ) − f ( x − y ) h d y
\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
h u ( x + h e i ) − u ( x ) = ∫ R n Φ ( y ) h f ( x + h e i − y ) − f ( x − y ) d y
이때 e i = ( 0 , ⋯ , 1 , ⋯ , 0 ) e_{i}=(0,\cdots ,1,\cdots, 0) e i = ( 0 , ⋯ , 1 , ⋯ , 0 ) 는 i i i 번째 성분만 1 1 1 이고 나머지 성분은 0 0 0 인 벡터이다. 그러면, f f f 가 미분가능하므로, 평균값 정리(MVT) 에 의해 임의의 y ∈ R n y \in \mathbb{R}^n y ∈ R n 와 h ′ ∈ ( 0 , h ) h^{\prime} \in (0,h) h ′ ∈ ( 0 , h ) 에 대해서 다음이 성립한다.
f ( x + h e i − y ) − f ( x − y ) h = f x i ( x + h ′ e i − y )
\dfrac{f(x+he_{i}-y) - f(x-y)}{h}=f_{x_{i}}(x+h^{\prime}e_{i}-y)
h f ( x + h e i − y ) − f ( x − y ) = f x i ( x + h ′ e i − y )
가정에 의해 f x i ∈ C c 1 ( R n ) f_{x_{i}} \in C_{c}^1 (\mathbb{R}^n) f x i ∈ C c 1 ( R n ) 이고, 컴팩트 공간에서 연속인 함수는 균등연속 이므로 f x i f_{x_{i}} f x i 는 R n \mathbb{R}^n R n 에서 균등 연속 이다. 따라서 주어진 ϵ > 0 \epsilon >0 ϵ > 0 에 대해서 ∣ z − w ∣ < δ ⟹ ∣ f x i ( z ) − f x i ( w ) ∣ < ϵ |z-w|<\delta \implies |f_{x_{i}}(z)-f_{x_{i}}(w)|<\epsilon ∣ z − w ∣ < δ ⟹ ∣ f x i ( z ) − f x i ( w ) ∣ < ϵ 을 만족하는 δ > 0 \delta >0 δ > 0 가 존재한다. 만약 0 < ∣ h ∣ < δ 0<|h|<\delta 0 < ∣ h ∣ < δ 이면 모든 y ∈ R n y\in \mathbb{R}^n y ∈ R n 에 대해서 다음이 성립한다.
∣ ( x + h ′ e i − y ) − ( x − y ) ∣ = ∣ h ′ ∣ < ∣ h ∣ < δ
|(x+h^{\prime}e_{i}-y)-(x-y)|=|h^{\prime}|<|h|<\delta
∣ ( x + h ′ e i − y ) − ( x − y ) ∣ = ∣ h ′ ∣ < ∣ h ∣ < δ
따라서 다음이 성립한다.
∣ ( x + h ′ e i − y ) − ( x − y ) ∣ = ∣ h ′ ∣ < ∣ h ∣ < δ ⟹ ∣ f x i ( x + h ′ e i − y ) − f x i ( x − y ) ∣ < ϵ ⟹ sup y ∈ R n ∣ f ( x + h e i − y ) − f ( x − y ) h − f x i ( x − y ) ∣ < ϵ
\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*}
⟹ ⟹ ∣ ( x + h ′ e i − y ) − ( x − y ) ∣ = ∣ h ′ ∣ < ∣ h ∣ ∣ f x i ( x + h ′ e i − y ) − f x i ( x − y ) ∣ y ∈ R n sup h f ( x + h e i − y ) − f ( x − y ) − f x i ( x − y ) < δ < ϵ < ϵ
그러므로 다음이 성립한다.
f ( x + h e i − y ) − f ( x − y ) h ⇉ f x i ( x − y ) a s h → 0
\dfrac{f(x+he_{i}-y) - f(x-y)}{h} \rightrightarrows f_{x_{i}}(x-y)\quad \mathrm{as}\quad h\rightarrow 0
h f ( x + h e i − y ) − f ( x − y ) ⇉ f x i ( x − y ) as h → 0
따라서 다음을 얻는다.
u x i ( x ) = lim h → 0 u ( x + h e i ) − u ( x ) h = ∫ R n Φ ( x ) f x i ( x − y ) d y ( 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*}
u x i ( x ) = h → 0 lim h u ( x + h e i ) − u ( x ) = ∫ R n Φ ( x ) f x i ( x − y ) d y ( i = 1 , ⋯ , )
같은 방식으로 다음이 성립한다.
u x i x j ( x ) = lim h → 0 u x i ( x + h e j ) − u x i ( x ) h = ∫ R n Φ ( x ) f x i x j ( x − y ) d y ( 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*}
u x i x j ( x ) = h → 0 lim h u x i ( x + h e j ) − u x i ( x ) = ∫ R n Φ ( x ) f x i x j ( x − y ) d y ( i , j = 1 , ⋯ , )
■