平均值定理在调和函数中的应用

题目

问题 2. 利用平均值定理的证明(参见第 7.2.4 小节),证明如果在球 B(y,r) B(y, r) B(y,r)Δu≥0 \Delta u \geq 0 Δu0,则

(a) u(y) u(y) u(y) 不超过 u u u 在该球边界球面 S(y,r) S(y, r) S(y,r) 上的平均值:

u(y)≤1σnrn−1∫S(y,r)u dS.u(y) \leq \frac{1}{\sigma_n r^{n-1}} \int_{S(y, r)} u \, dS.u(y)σnrn11S(y,r)udS.

(b) u(y) u(y) u(y) 不超过 u u u 在该球 B(y,r) B(y, r) B(y,r) 上的平均值:

u(y)≤1ωnrn∫B(y,r)u dV.u(y) \leq \frac{1}{\omega_n r^n} \int_{B(y, r)} u \, dV.u(y)ωnrn1B(y,r)udV.

© 对于满足 Δu≤0 \Delta u \leq 0 Δu0 的函数,给出类似的陈述(在下一个问题中,我们将它们称为 (a)’ 和 (b)’’)。

其中:

  • σn \sigma_n σnn n n 维单位球面的表面积(即 Rn \mathbb{R}^n Rn 中单位球面的 (n−1) (n-1) (n1) 维测度),
  • ωn \omega_n ωnn n n 维单位球的体积(即 Rn \mathbb{R}^n Rn 中单位球的 n n n 维测度),
  • Δu=∑i=1n∂2u∂xi2 \Delta u = \sum_{i=1}^n \frac{\partial^2 u}{\partial x_i^2} Δu=i=1nxi22u 是 Laplace 算子,
  • B(y,r)={x∈Rn:∣x−y∣B(y,r)={xRn:xy<r} 是以 y y y 为中心、r r r 为半径的开球,
  • S(y,r)={x∈Rn:∣x−y∣=r} S(y, r) = \{ x \in \mathbb{R}^n : |x - y| = r \} S(y,r)={xRn:xy=r} 是球 B(y,r) B(y, r) B(y,r) 的边界球面。

解决题目

我们利用调和函数平均值定理的证明方法(见第 7.2.4 小节)。假设 u u u 是足够光滑的函数(例如,二阶连续可微),以确保所有导数和积分操作有效。给定条件 Δu≥0 \Delta u \geq 0 Δu0B(y,r) B(y, r) B(y,r) 中成立,表明 u u u 是次调和函数(subharmonic function)。

(a) 证明 u(y)≤1σnrn−1∫S(y,r)u dS u(y) \leq \frac{1}{\sigma_n r^{n-1}} \int_{S(y, r)} u \, dS u(y)σnrn11S(y,r)udS

定义球面平均值函数:
ϕ(r)=1σnrn−1∫S(y,r)u(x)dS(x). \phi(r) = \frac{1}{\sigma_n r^{n-1}} \int_{S(y, r)} u(x) dS(x). ϕ(r)=σnrn11S(y,r)u(x)dS(x).
通过坐标变换 x=y+rξ x = y + r \xi x=y+rξ(其中 ∣ξ∣=1 |\xi| = 1 ξ=1,球面面积元素变换为 dS(x)=rn−1dS(ξ) dS(x) = r^{n-1} dS(\xi) dS(x)=rn1dS(ξ),且单位球面 ∣ξ∣=1 |\xi| = 1 ξ=1 的面积为 σn \sigma_n σn,因此:
ϕ(r)=1σn∫∣ξ∣=1u(y+rξ)dS(ξ). \phi(r) = \frac{1}{\sigma_n} \int_{|\xi|=1} u(y + r \xi) dS(\xi). ϕ(r)=σn1ξ=1u(y+rξ)dS(ξ).
r r r 求导:
ϕ′(r)=1σn∫∣ξ∣=1∇u(y+rξ)⋅ξdS(ξ). \phi'(r) = \frac{1}{\sigma_n} \int_{|\xi|=1} \nabla u(y + r \xi) \cdot \xi dS(\xi). ϕ(r)=σn1ξ=1u(y+rξ)ξdS(ξ).
在球面 S(y,r) S(y, r) S(y,r) 上,单位外法向量为 n=x−yr=ξ \mathbf{n} = \frac{x - y}{r} = \xi n=rxy=ξ,因此:
ϕ′(r)=1σnrn−1∫S(y,r)∇u(x)⋅ndS(x). \phi'(r) = \frac{1}{\sigma_n r^{n-1}} \int_{S(y, r)} \nabla u(x) \cdot \mathbf{n} dS(x). ϕ(r)=σnrn11S(y,r)u(x)ndS(x).
由散度定理(Gauss 定理):
∫S(y,r)∇u⋅ndS=∫B(y,r)ΔudV. \int_{S(y, r)} \nabla u \cdot \mathbf{n} dS = \int_{B(y, r)} \Delta u dV. S(y,r)undS=B(y,r)ΔudV.
代入得:
ϕ′(r)=1σnrn−1∫B(y,r)ΔudV. \phi'(r) = \frac{1}{\sigma_n r^{n-1}} \int_{B(y, r)} \Delta u dV. ϕ(r)=σnrn11B(y,r)ΔudV.
给定 Δu≥0 \Delta u \geq 0 Δu0B(y,r) B(y, r) B(y,r) 中,因此:
∫B(y,r)ΔudV≥0  ⟹  ϕ′(r)≥0. \int_{B(y, r)} \Delta u dV \geq 0 \implies \phi'(r) \geq 0. B(y,r)ΔudV0ϕ(r)0.
这表明 ϕ(r) \phi(r) ϕ(r)r r r 的非减函数。当 r→0+ r \to 0^+ r0+ 时,球面收缩到点 y y y,且 u u u 连续,故:
lim⁡r→0+ϕ(r)=u(y). \lim_{r \to 0^+} \phi(r) = u(y). r0+limϕ(r)=u(y).
由于 ϕ(r) \phi(r) ϕ(r) 非减,对任意 r>0 r > 0 r>0,有:
ϕ(r)≥lim⁡ρ→0+ϕ(ρ)=u(y). \phi(r) \geq \lim_{\rho \to 0^+} \phi(\rho) = u(y). ϕ(r)ρ0+limϕ(ρ)=u(y).
即:
u(y)≤1σnrn−1∫S(y,r)udS. u(y) \leq \frac{1}{\sigma_n r^{n-1}} \int_{S(y, r)} u dS. u(y)σnrn11S(y,r)udS.

(b) 证明 u(y)≤1ωnrn∫B(y,r)udV u(y) \leq \frac{1}{\omega_n r^n} \int_{B(y, r)} u dV u(y)ωnrn1B(y,r)udV

定义体积平均值函数:
ψ(r)=1ωnrn∫B(y,r)u(x)dV(x). \psi(r) = \frac{1}{\omega_n r^n} \int_{B(y, r)} u(x) dV(x). ψ(r)=ωnrn1B(y,r)u(x)dV(x).
球体积积分可分解为球面积分:
∫B(y,r)udV=∫0r∫S(y,s)udSds. \int_{B(y, r)} u dV = \int_0^r \int_{S(y, s)} u dS ds. B(y,r)udV=0rS(y,s)udSds.
由 (a) 中的 ϕ(s) \phi(s) ϕ(s),有 ∫S(y,s)udS=ϕ(s)σnsn−1 \int_{S(y, s)} u dS = \phi(s) \sigma_n s^{n-1} S(y,s)udS=ϕ(s)σnsn1,且 σn=nωn \sigma_n = n \omega_n σn=nωn(因为单位球体积 ωn=1nσn \omega_n = \frac{1}{n} \sigma_n ωn=n1σn),因此:
∫B(y,r)udV=∫0rϕ(s)σnsn−1ds=∫0rϕ(s)nωnsn−1ds. \int_{B(y, r)} u dV = \int_0^r \phi(s) \sigma_n s^{n-1} ds = \int_0^r \phi(s) n \omega_n s^{n-1} ds. B(y,r)udV=0rϕ(s)σnsn1ds=0rϕ(s)nωnsn1ds.
代入 ψ(r) \psi(r) ψ(r)
ψ(r)=1ωnrn∫0rϕ(s)nωnsn−1ds=nrn∫0rϕ(s)sn−1ds. \psi(r) = \frac{1}{\omega_n r^n} \int_0^r \phi(s) n \omega_n s^{n-1} ds = \frac{n}{r^n} \int_0^r \phi(s) s^{n-1} ds. ψ(r)=ωnrn10rϕ(s)nωnsn1ds=rnn0rϕ(s)sn1ds.
由 (a),对任意 s∈(0,r) s \in (0, r) s(0,r),球 B(y,s)⊂B(y,r) B(y, s) \subset B(y, r) B(y,s)B(y,r)Δu≥0 \Delta u \geq 0 Δu0B(y,r) B(y, r) B(y,r) 中,故在 B(y,s) B(y, s) B(y,s)Δu≥0 \Delta u \geq 0 Δu0,应用 (a) 得 ϕ(s)≥u(y) \phi(s) \geq u(y) ϕ(s)u(y)。因此:
ψ(r)≥nrn∫0ru(y)sn−1ds=u(y)nrn⋅rnn=u(y). \psi(r) \geq \frac{n}{r^n} \int_0^r u(y) s^{n-1} ds = u(y) \frac{n}{r^n} \cdot \frac{r^n}{n} = u(y). ψ(r)rnn0ru(y)sn1ds=u(y)rnnnrn=u(y).
即:
u(y)≤1ωnrn∫B(y,r)udV. u(y) \leq \frac{1}{\omega_n r^n} \int_{B(y, r)} u dV. u(y)ωnrn1B(y,r)udV.

© 对 Δu≤0 \Delta u \leq 0 Δu0 的函数的类似陈述

如果 Δu≤0 \Delta u \leq 0 Δu0B(y,r) B(y, r) B(y,r) 中(即 u u u 是超调和函数,superharmonic),则类似推导可得反向不等式:

  • (a)’ u(y) u(y) u(y) 不小于 u u u 在球面 S(y,r) S(y, r) S(y,r) 上的平均值:
    u(y)≥1σnrn−1∫S(y,r)udS. u(y) \geq \frac{1}{\sigma_n r^{n-1}} \int_{S(y, r)} u dS. u(y)σnrn11S(y,r)udS.
    证明概要: 类似 (a),ϕ′(r)=1σnrn−1∫B(y,r)ΔudV≤0 \phi'(r) = \frac{1}{\sigma_n r^{n-1}} \int_{B(y, r)} \Delta u dV \leq 0 ϕ(r)=σnrn11B(y,r)ΔudV0,故 ϕ(r) \phi(r) ϕ(r) 非增,且 lim⁡r→0+ϕ(r)=u(y) \lim_{r \to 0^+} \phi(r) = u(y) r0+limϕ(r)=u(y),因此 ϕ(r)≤u(y) \phi(r) \leq u(y) ϕ(r)u(y).

  • (b)’’ u(y) u(y) u(y) 不小于 u u u 在球 B(y,r) B(y, r) B(y,r) 上的平均值:
    u(y)≥1ωnrn∫B(y,r)udV. u(y) \geq \frac{1}{\omega_n r^n} \int_{B(y, r)} u dV. u(y)ωnrn1B(y,r)udV.
    证明概要: 类似 (b),ψ(r)=nrn∫0rϕ(s)sn−1ds \psi(r) = \frac{n}{r^n} \int_0^r \phi(s) s^{n-1} ds ψ(r)=rnn0rϕ(s)sn1ds,由 (a)’ 有 ϕ(s)≤u(y) \phi(s) \leq u(y) ϕ(s)u(y),因此 ψ(r)≤u(y) \psi(r) \leq u(y) ψ(r)u(y).

这些陈述在下一个问题中称为 (a)’ 和 (b)’’。

你可能感兴趣的:(算法)