基于分析Laplace方程“放射状”函数特解的基本解引入
作者：互联网
基于分析Laplace方程“放射状”函数特解的基本解引入
参考文献：【偏微分方程笔记(2)——Laplace(位势)方程的基本解】
1. 基本定义
关于函数
u
(
x
1
,
x
2
,
.
.
.
,
x
n
)
u(x_1,x_2,...,x_n)
u(x1,x2,...,xn)的
n
n
n维Laplace方程是:
Δ
u
=
u
x
1
x
1
+
u
x
2
x
2
+
.
.
.
+
u
x
n
x
n
=
0
(1)
\Delta u=u_{x_1x_1}+u_{x_2x_2}+...+u_{x_nx_n}=0 \tag 1
Δu=ux1x1+ux2x2+...+uxnxn=0(1)
它的解称为调和函数或势函数。
2. 寻找特解
由于Laplace方程作旋转以后是不变的，因此我们似乎可以先去找“放射状（radial）”的函数，也就是：
r
=
∣
x
∣
=
x
1
2
+
.
.
.
+
x
n
2
(2)
r=|x|=\sqrt{x_1^2+...+x_n^2} \tag 2
r=∣x∣=x12+...+xn2
(2)
首先尝试在
n
n
n维线性空间中寻找满足Laplace方程
(
1
)
(1)
(1)的一个解
u
u
u，具有形式：
u
(
x
)
=
v
(
r
)
(3)
u(x)=v(r) \tag 3
u(x)=v(r)(3)
首先注意到对于
i
=
1
,
2
,
.
.
.
,
n
i=1,2,...,n
i=1,2,...,n，有：
∂
r
∂
x
i
=
1
2
x
1
2
+
.
.
.
+
x
n
2
2
x
i
=
x
i
r
(
x
≠
0
)
(4)
\frac{\partial r}{\partial x_i} = \frac{1}{2\sqrt{x_1^2+...+x_n^2}} 2x_i = \frac{x_i}{r}(x \not =0) \tag 4
∂xi∂r=2x12+...+xn2
12xi=rxi(x=0)(4)
因此：
u
(
x
)
u(x)
u(x)的一阶导数有：
u
x
i
=
∂
v
∂
r
∂
r
∂
x
i
=
v
′
(
r
)
x
i
r
(5)
u_{x_i}=\frac{\partial v}{\partial r} \frac{\partial r}{\partial x_i} = v'(r) \frac{x_i}{r} \tag 5
uxi=∂r∂v∂xi∂r=v′(r)rxi(5)
u
(
x
)
u(x)
u(x)的二阶导数有：
u
x
i
x
i
=
∂
u
x
i
∂
x
i
=
∂
∂
x
i
(
∂
v
∂
r
∂
r
∂
x
i
)
=
∂
2
v
∂
x
i
∂
r
∂
r
∂
x
i
+
∂
2
r
∂
x
i
∂
x
i
∂
v
∂
r
=
[
∂
∂
r
∂
r
∂
x
i
(
v
′
(
r
)
)
]
⋅
x
i
r
+
[
∂
∂
x
i
(
x
i
r
)
]
⋅
v
′
(
r
)
(6)
u_{x_ix_i}=\frac{\partial u_{x_i}}{\partial x_i}=\frac{\partial}{\partial x_i} (\frac{\partial v}{\partial r} \frac{\partial r}{\partial x_i}) \\ = \frac{\partial^2 v}{\partial x_i \partial r} \frac{\partial r}{\partial x_i} + \frac{\partial^2 r}{\partial x_i \partial x_i} \frac{\partial v}{\partial r} \\ =[\frac{\partial}{\partial r} \frac{\partial r}{\partial x_i}(v'(r))]·\frac{x_i}{r} + [\frac{\partial}{\partial x_i} (\frac{x_i}{r})]·v'(r) \tag 6
uxixi=∂xi∂uxi=∂xi∂(∂r∂v∂xi∂r)=∂xi∂r∂2v∂xi∂r+∂xi∂xi∂2r∂r∂v=[∂r∂∂xi∂r(v′(r))]⋅rxi+[∂xi∂(rxi)]⋅v′(r)(6)
对于前式中括号中的式子：
∂
∂
r
(
x
i
v
′
(
r
)
r
)
=
x
i
[
v
′
′
(
r
)
r
−
v
′
(
r
)
]
r
2
(7)
\frac{\partial}{\partial r} (\frac{x_i v'(r)}{r}) = \frac{x_i[v^{''}(r)r-v'(r)]}{r^2} \tag 7
∂r∂(rxiv′(r))=r2xi[v′′(r)r−v′(r)](7)
对于后式中括号中的式子：
∂
∂
x
i
(
x
i
r
)
=
1
r
(8)
\frac{\partial}{\partial x_i} (\frac{x_i}{r}) = \frac{1}{r} \tag 8
∂xi∂(rxi)=r1(8)
故
u
(
x
)
u(x)
u(x)的二阶偏导数：
u
x
i
x
i
=
v
′
′
(
r
)
x
i
2
r
2
+
v
′
(
r
)
(
1
r
−
x
i
2
r
3
)
(9)
u_{x_ix_i}=v^{''}(r) \frac{x_i^2}{r^2} + v'(r) (\frac{1}{r}-\frac{x_i^2}{r^3}) \tag 9
uxixi=v′′(r)r2xi2+v′(r)(r1−r3xi2)(9)
因此：
Δ
u
=
v
′
′
(
r
)
+
n
−
1
r
v
′
(
r
)
(10)
\Delta u=v^{''}(r)+\frac{n-1}{r} v'(r) \tag {10}
Δu=v′′(r)+rn−1v′(r)(10)
因此
Δ
u
=
0
\Delta u=0
Δu=0当且仅当
v
′
′
+
n
−
1
r
v
′
=
0
(11)
v^{''}+\frac{n-1}{r} v'=0 \tag {11}
v′′+rn−1v′=0(11)
3. 解常微分方程
如果
v
′
≠
0
v' \not =0
v′=0，则：
l
n
(
∣
v
′
∣
)
′
=
v
′
′
v
′
=
1
−
n
r
(12)
ln(|v'|)'=\frac{v^{''}}{v'}=\frac{1-n}{r} \tag {12}
ln(∣v′∣)′=v′v′′=r1−n(12)
存在常数
a
a
a使得：
v
′
(
r
)
=
a
r
n
−
1
(13)
v'(r)=\frac{a}{r^{n-1}} \tag {13}
v′(r)=rn−1a(13)
因此如果
r
>
0
r>0
r>0，我们有：
当
n
=
1
n=1
n=1时：
v
(
r
)
=
a
r
(14)
v(r)=ar \tag {14}
v(r)=ar(14)
当
n
=
2
n=2
n=2时：
v
(
r
)
=
b
l
n
r
+
c
(15)
v(r)=b \ ln r+c \tag {15}
v(r)=b lnr+c(15)
当
n
≥
3
n \geq 3
n≥3时：
v
(
r
)
=
b
r
n
−
2
+
c
(16)
v(r)=\frac{b}{r^{n-2}}+c \tag {16}
v(r)=rn−2b+c(16)
这里
b
,
c
b,c
b,c均为常数。
4. 基本解
该方程组称为Laplace方程的基本解，这里
α
(
n
)
=
R
n
\alpha(n)=R^n
α(n)=Rn中单位球的体积=
π
n
τ
(
n
2
+
1
)
\frac{\sqrt{\pi^n}}{\tau (\frac{n}{2}+1)}
τ(2n+1)πn
。
标签：xi,partial,特解,放射状,tag,frac,+...+,Laplace 来源： https://blog.csdn.net/linjing_zyq/article/details/120356477