J u m p t o c o n t e n t
M a i n m e n u
M a i n m e n u
N a v i g a t i o n
● M a i n p a g e
● C o n t e n t s
● C u r r e n t e v e n t s
● R a n d o m a r t i c l e
● A b o u t W i k i p e d i a
● C o n t a c t u s
● D o n a t e
C o n t r i b u t e
● H e l p
● L e a r n t o e d i t
● C o m m u n i t y p o r t a l
● R e c e n t c h a n g e s
● U p l o a d f i l e
S e a r c h
Search
A p p e a r a n c e
● C r e a t e a c c o u n t
● L o g i n
P e r s o n a l t o o l s
● C r e a t e a c c o u n t
● L o g i n
P a g e s f o r l o g g e d o u t e d i t o r s l e a r n m o r e
● C o n t r i b u t i o n s
● T a l k
( T o p )
1
I n t r o d u c t i o n : K a k u t a n i ' s s o l u t i o n t o t h e c l a s s i c a l D i r i c h l e t p r o b l e m
2
T h e D i r i c h l e t – P o i s s o n p r o b l e m
3
R e f e r e n c e s
T o g g l e t h e t a b l e o f c o n t e n t s
S t o c h a s t i c p r o c e s s e s a n d b o u n d a r y v a l u e p r o b l e m s
A d d l a n g u a g e s
A d d l i n k s
● A r t i c l e
● T a l k
E n g l i s h
● R e a d
● E d i t
● V i e w h i s t o r y
T o o l s
T o o l s
A c t i o n s
● R e a d
● E d i t
● V i e w h i s t o r y
G e n e r a l
● W h a t l i n k s h e r e
● R e l a t e d c h a n g e s
● U p l o a d f i l e
● S p e c i a l p a g e s
● P e r m a n e n t l i n k
● P a g e i n f o r m a t i o n
● C i t e t h i s p a g e
● G e t s h o r t e n e d U R L
● D o w n l o a d Q R c o d e
● W i k i d a t a i t e m
P r i n t / e x p o r t
● D o w n l o a d a s P D F
● P r i n t a b l e v e r s i o n
A p p e a r a n c e
F r o m W i k i p e d i a , t h e f r e e e n c y c l o p e d i a
Let
D
{\displaystyle D}
be a domain (an open and connected set ) in
R
n
{\textstyle \mathbb {R} ^{n}}
. Let
Δ
{\displaystyle \Delta }
be the Laplace operator , let
g
{\displaystyle g}
be a bounded function on the boundary
∂
D
{\displaystyle \partial D}
, and consider the problem:
{
−
Δ
u
(
x
)
=
0
,
x
∈
D
lim
y
→
x
u
(
y
)
=
g
(
x
)
,
x
∈
∂
D
{\displaystyle {\begin{cases}-\Delta u(x )=0,&x\in D\\\displaystyle {\lim _{y\to x}u(y )}=g(x ),&x\in \partial D\end{cases}}}
It can be shown that if a solution
u
{\displaystyle u}
exists, then
u
(
x
)
{\displaystyle u(x )}
is the expected value of
g
(
x
)
{\displaystyle g(x )}
at the (random) first exit point from
D
{\displaystyle D}
for a canonical Brownian motion starting at
x
{\displaystyle x}
. See theorem 3 in Kakutani 1944, p. 710.
The Dirichlet–Poisson problem [ edit ]
Let
D
{\displaystyle D}
be a domain in
R
n
{\textstyle \mathbb {R} ^{n}}
and let
L
{\displaystyle L}
be a semi-elliptic differential operator on
C
2
(
R
n
;
R
)
{\textstyle C^{2}(\mathbb {R} ^{n};\mathbb {R} )}
of the form:
L
=
∑
i
=
1
n
b
i
(
x
)
∂
∂
x
i
+
∑
i
,
j
=
1
n
a
i
j
(
x
)
∂
2
∂
x
i
∂
x
j
{\displaystyle L=\sum _{i=1}^{n}b_{i}(x ){\frac {\partial }{\partial x_{i}}}+\sum _{i,j=1}^{n}a_{ij}(x ){\frac {\partial ^{2}}{\partial x_{i}\,\partial x_{j}}}}
where the coefficients
b
i
{\displaystyle b_{i}}
and
a
i
j
{\displaystyle a_{ij}}
are continuous functions and all the eigenvalues of the matrix
α
(
x
)
=
a
i
j
(
x
)
{\displaystyle \alpha (x )=a_{ij}(x )}
are non-negative. Let
f
∈
C
(
D
;
R
)
{\textstyle f\in C(D;\mathbb {R} )}
and
g
∈
C
(
∂
D
;
R
)
{\textstyle g\in C(\partial D;\mathbb {R} )}
. Consider the Poisson problem :
{
−
L
u
(
x
)
=
f
(
x
)
,
x
∈
D
lim
y
→
x
u
(
y
)
=
g
(
x
)
,
x
∈
∂
D
(P1)
{\displaystyle {\begin{cases}-Lu(x )=f(x ),&x\in D\\\displaystyle {\lim _{y\to x}u(y )}=g(x ),&x\in \partial D\end{cases}}\quad {\mbox{(P1)}}}
The idea of the stochastic method for solving this problem is as follows. First, one finds an Itō diffusion
X
{\displaystyle X}
whose infinitesimal generator
A
{\displaystyle A}
coincides with
L
{\displaystyle L}
on compactly-supported
C
2
{\displaystyle C^{2}}
functions
f
:
R
n
→
R
{\displaystyle f:\mathbb {R} ^{n}\rightarrow \mathbb {R} }
. For example,
X
{\displaystyle X}
can be taken to be the solution to the stochastic differential equation:
d
X
t
=
b
(
X
t
)
d
t
+
σ
(
X
t
)
d
B
t
{\displaystyle \mathrm {d} X_{t}=b(X_{t})\,\mathrm {d} t+\sigma (X_{t})\,\mathrm {d} B_{t}}
where
B
{\displaystyle B}
is n -dimensional Brownian motion,
b
{\displaystyle b}
has components
b
i
{\displaystyle b_{i}}
as above, and the matrix field
σ
{\displaystyle \sigma }
is chosen so that:
1
2
σ
(
x
)
σ
(
x
)
⊤
=
a
(
x
)
,
∀
x
∈
R
n
{\displaystyle {\frac {1}{2}}\sigma (x )\sigma (x )^{\top }=a(x ),\quad \forall x\in \mathbb {R} ^{n}}
For a point
x
∈
R
n
{\displaystyle x\in \mathbb {R} ^{n}}
, let
P
x
{\displaystyle \mathbb {P} ^{x}}
denote the law of
X
{\displaystyle X}
given initial datum
X
0
=
x
{\displaystyle X_{0}=x}
, and let
E
x
{\displaystyle \mathbb {E} ^{x}}
denote expectation with respect to
P
x
{\displaystyle \mathbb {P} ^{x}}
. Let
τ
D
{\displaystyle \tau _{D}}
denote the first exit time of
X
{\displaystyle X}
from
D
{\displaystyle D}
.
In this notation, the candidate solution for (P1) is:
u
(
x
)
=
E
x
[
g
(
X
τ
D
)
⋅
χ
{
τ
D
<
+
∞
}
]
+
E
x
[
∫
0
τ
D
f
(
X
t
)
d
t
]
{\displaystyle u(x )=\mathbb {E} ^{x}\left[g{\big (}X_{\tau _{D}}{\big )}\cdot \chi _{\{\tau _{D}<+\infty \}}\right]+\mathbb {E} ^{x}\left[\int _{0}^{\tau _{D}}f(X_{t})\,\mathrm {d} t\right]}
provided that
g
{\displaystyle g}
is a bounded function and that:
E
x
[
∫
0
τ
D
|
f
(
X
t
)
|
d
t
]
<
+
∞
{\displaystyle \mathbb {E} ^{x}\left[\int _{0}^{\tau _{D}}{\big |}f(X_{t}){\big |}\,\mathrm {d} t\right]<+\infty }
It turns out that one further condition is required:
P
x
(
τ
D
<
∞
)
=
1
,
∀
x
∈
D
{\displaystyle \mathbb {P} ^{x}{\big (}\tau _{D}<\infty {\big )}=1,\quad \forall x\in D}
For all
x
{\displaystyle x}
, the process
X
{\displaystyle X}
starting at
x
{\displaystyle x}
almost surely leaves
D
{\displaystyle D}
in finite time. Under this assumption, the candidate solution above reduces to:
u
(
x
)
=
E
x
[
g
(
X
τ
D
)
]
+
E
x
[
∫
0
τ
D
f
(
X
t
)
d
t
]
{\displaystyle u(x )=\mathbb {E} ^{x}\left[g{\big (}X_{\tau _{D}}{\big )}\right]+\mathbb {E} ^{x}\left[\int _{0}^{\tau _{D}}f(X_{t})\,\mathrm {d} t\right]}
and solves (P1) in the sense that if
A
{\displaystyle {\mathcal {A}}}
denotes the characteristic operator for
X
{\displaystyle X}
(which agrees with
A
{\displaystyle A}
on
C
2
{\displaystyle C^{2}}
functions), then:
{
−
A
u
(
x
)
=
f
(
x
)
,
x
∈
D
lim
t
↑
τ
D
u
(
X
t
)
=
g
(
X
τ
D
)
,
P
x
-a.s.,
∀
x
∈
D
(P2)
{\displaystyle {\begin{cases}-{\mathcal {A}}u(x )=f(x ),&x\in D\\\displaystyle {\lim _{t\uparrow \tau _{D}}u(X_{t})}=g{\big (}X_{\tau _{D}}{\big )},&\mathbb {P} ^{x}{\mbox{-a.s.,}}\;\forall x\in D\end{cases}}\quad {\mbox{(P2)}}}
Moreover, if
v
∈
C
2
(
D
;
R
)
{\textstyle v\in C^{2}(D;\mathbb {R} )}
satisfies (P2) and there exists a constant
C
{\displaystyle C}
such that, for all
x
∈
D
{\displaystyle x\in D}
:
|
v
(
x
)
|
≤
C
(
1
+
E
x
[
∫
0
τ
D
|
g
(
X
s
)
|
d
s
]
)
{\displaystyle |v(x )|\leq C\left(1+\mathbb {E} ^{x}\left[\int _{0}^{\tau _{D}}{\big |}g(X_{s}){\big |}\,\mathrm {d} s\right]\right)}
then
v
=
u
{\displaystyle v=u}
.
References [ edit ]
R e t r i e v e d f r o m " https://en.wikipedia.org/w/index.php?title=Stochastic_processes_and_boundary_value_problems&oldid=966674690 "
C a t e g o r i e s :
● B o u n d a r y v a l u e p r o b l e m s
● P a r t i a l d i f f e r e n t i a l e q u a t i o n s
● S t o c h a s t i c d i f f e r e n t i a l e q u a t i o n s
● T h i s p a g e w a s l a s t e d i t e d o n 8 J u l y 2 0 2 0 , a t 1 3 : 5 5 ( U T C ) .
● T e x t i s a v a i l a b l e u n d e r t h e C r e a t i v e C o m m o n s A t t r i b u t i o n - S h a r e A l i k e L i c e n s e 4 . 0 ;
a d d i t i o n a l t e r m s m a y a p p l y . B y u s i n g t h i s s i t e , y o u a g r e e t o t h e T e r m s o f U s e a n d P r i v a c y P o l i c y . W i k i p e d i a ® i s a r e g i s t e r e d t r a d e m a r k o f t h e W i k i m e d i a F o u n d a t i o n , I n c . , a n o n - p r o f i t o r g a n i z a t i o n .
● P r i v a c y p o l i c y
● A b o u t W i k i p e d i a
● D i s c l a i m e r s
● C o n t a c t W i k i p e d i a
● C o d e o f C o n d u c t
● D e v e l o p e r s
● S t a t i s t i c s
● C o o k i e s t a t e m e n t
● M o b i l e v i e w