我们可以将上述紧算子的理论应用到 Dirichlet 问题之上 (我们在这个章节并不需要区域是光滑的, Ω⊂Rn 只要是有界的开集即可, 而我们只在 H01(Ω) 中研究问题而不再关心函数在边界 ∂Ω 上的限制) . 为此, 我们给出一个紧算子的例子. 从某种意义上说, 这是最重要的一类紧算子的例子:
任意给定有界开区域 Ω⊂Rn, 那么, 自然的嵌入映射ι:H01(Ω)→L2(Ω), u↦u,是紧算子.
证明. 我们将利用 Fourier 级数的理论来证明这个重要的定理. 首先, 我们不妨假设 Ω⊂⊂(0,2π)n (否则, 我们可以把 Ω 放到更大的一个 Rn 中的一个盒子中去, 从而用一个周期与 2π 不同的 Fourier 级数即可) .
我们可以把 C0∞(Ω) 中的函数在 Rn−Ω 中用 0 来延拓, 这给出连续映射Ext:C0∞(Ω)→H1(Rn).其中, 上述连续性之所以成立是因为我们可以在 Ω 上运用 Poincaré 不等式. 根据 C0∞(Ω) 在 H01(Ω) 中的稠密性, 我们就得到了等距嵌入ι:(H01(Ω),∥⋅∥H1)↪H1(Rn).由于 Ω⊂(0,2π)n 是相对紧的, 我们选取 χ∈C0∞((0,2π)n), 使得 χ∣∣Ω≡1,0⩽χ⩽1. 我们现在证明映射Tχ:H1(Rn)→L2(Rn), f(x)↦χ(x)f(x),是紧算子. 我们注意到, 作为 L2((0,2π)n) 中的函数, 利用 Fourier 级数的展开, 我们有χ(x)f(x)=χ(x)k∈Zn∑ckeik⋅x=χ(x)k∈Zn∑((2π)n1∫(0,2π)nf(y)e−ik⋅ydy)eik⋅x==Tχ,N(f)∣k∣<N∑ckχ(x)eik⋅x+=Rχ,N(f)χ(x)∣k∣⩾N∑ckeik⋅x.我们把 Tχ 分解成两个算子的和, 其中, Tχ,N:H1(Rn)→L2(Rn), f(x)↦∣k∣<N∑ckχ(x)eik⋅x,而Rχ,N:H1(Rn)→L2(Rn), f(x)↦χ(x)∣k∣⩾N∑ckeik⋅x.这里, 整数 N 是待定的. 由于 Tχ,N 的像完全落在 {eik⋅x}∣k∣<N 这有限多个函数所生成的有限维线性空间中, 所以它是有限秩的算子, 从而, Tχ,N 是紧算子.
为了研究 Rχ,N, 我们将利用 ∇f 也是 L2 的函数这个事实. 实际上, 我们有∥Rχ,N(f)∥L2(Rn)2⩽∥∥∣k∣⩾N∑ckeik⋅x∥∥L2((0,2π)n)2⩽N21∣k∣⩾N∑∣k∣2∣ck∣2⩽N21∥∇f∥L22.其中, 我们把导数转化为频率空间上衰减. 这表明, Tχ 这个算子可以用一列紧算子 Tχ,N 来逼近, 其中, N→∞, 所以 Tχ 为紧算子.
最后, 我们把嵌入
ι:H01(Ω)→L2(Ω)写成如下两个连续线性映射的复合:
H01(Ω)→H1(Rn)⟶TχL2(Ω)⊂L2(Rn).由于任何算子与紧算子复合之后是紧算子, 所以我们就证明了
ι 也是紧算子.
我们现在考虑如下的交换图表: H−1(Ω)H01(Ω)L2(Ω)L2(Ω)(−△)−1ιι(−△)−1据此, 通过复合, 我们就可以定义(−△)−1:L2(Ω)⟶L2(Ω).从此之后, 当我们谈论 (−△)−1 的时候, 我们总假设它的定义域和值域都是 L2(Ω).
由于右边竖列的箭头是紧算子, 复合之后的算子 (−△)−1 也是紧的. 我们现在验证 (−△)−1 是自伴算子: 对任意的 f1,f2∈L2(Ω), 令 u1=(−△)−1f1, u2=(−△)−1f2, 我们知道, u1,u2∈H01(Ω), 所以, ((−△)−1f1,f2)=−(u1,△u2)=(∇u1,∇u2),这里, 我们用到了 u1∈H01(Ω) 事实, 请参考引理 73.7. 类似地, 我们还有(f1,(−△)−1f2)=(∇u1,∇u2).这就证明了自伴性. 此时, 就可以对 (−△)−1 使用抽象的紧算子理论了.
另外, 我们之前还证明了 (−△)−1 是连续的正线性算子, 即对任意的 f∈L2(Ω), 我们有((−△)−1f,f)L2⩾0.这个不等式中等号成立当且仅当 f=0. 这表明, (−△)−1 的所有特征值都是正实数 (没有 0!) .
运用 Hilbert–Schmidt 定理, 那么, 我们可把 (−△)−1 的特征值的集合写成σ((−△)−1)={μk}k⩾1使得μ1⩾μ2⩾μ3⩾⋯, k→∞limμk=0.我们用 φk(x)∈L2(Ω) 表示与 μk 相对应的长度 (L2-范数) 为 1 的特征函数. 我们再令 λk=μk−1, 所以, −△φk=λkφk.由于 φk∈L2(Ω), 根据 Dirichlet 问题的解, 我们知道 φk∈H01(Ω). 综上所述, 我们证明了如下关于 (−△)−1 的谱分解定理:
假设 Ω⊂Rn 是有界的开区域, 那么, 存在单调上升的无界序列0<λ1⩽λ2⩽λ3⩽⋯, k→∞limλk=+∞,以及 L2(Ω) 的一组 Hilbert 基 {φk}k⩾1, 使得−△φk=λkφk.进一步, 对任意的 k⩾1, 我们有 φk∈H01(Ω).
我们现在考虑 Ω=(0,d)n 是一个正方体的情况, 其中 d>0 是常数. (我们注意到 Ω 的边界并不是光滑的. ) 我们要详细地计算 −△ 在 Ω 上的所有特征值和特征函数. 根据 Fourier 级数的理论, L2(Ω) 的一族 Hilbert 基可以取作{d−2ned2πik⋅x}k∈Zn.由于我们要求 φ∣∣∂Ω=0, 所以, 我们希望选取 sin 型的函数.
先研究 n=1 的情形, 我们选取 {sin(dkπx)}k∈Z⩾1 作为备选. 我们指出这里和 Fourier 级数的不同之处: 指数中的 2π 变成了 π. 很明显, 我们有
• | −△(sin(dkπx))=d2π2k2sin(dkπx); |
• | sin(dkπx)∣∣∂(0,d)=0, 所以, 根据一维的 H01 空间的描述, 我们有 sin(dkπx)∈H01((0,d)); |
• | 当 k=l 时, 我们有如下的正交性: (sin(dkπx),sin(dlπx))L2=∫0dsin(dkπx)sin(dlπx)dx=0. |
为了说明这给出了 −△=−dx2d2 的所有特征函数, 我们在区间 (0,d) 上解 (常) 微分方程−u′′=λ2u,其中 λ>0. 注意到, 我们要找的 u 满足 u∈H01((0,d)), 根据 Sobolev 嵌入定理, u∈C0((0,d)). 据此, 根据方程, u′′∈C0, 所以, u∈C2, 再代入方程, u′′∈C2, 所以, u∈C4. 如此迭代, 我们知道 u 是足够光滑的函数, 从而可以能用经典的常微分方程理论 (解存在唯一) . 所以, 这个方程的通解可写成u(x)=Aeiλx+Be−iλx,其中, A,B 是复数. 再利用 u(0)=u(d)=0, 我们有A+B=0, Aeiλd+Be−iλd=0,这说明B=−A, Asin(λd)=0.由于 A=0, 所以, u(x)=Csin(dπkx).这说明, 我们已经列出了所有的特征函数和特征值.
现在假设维数是 n, 我们记 k=(k1,⋯,kn)∈(Z⩾1)n, 那么{j=1∏nsin(dπkj⋅xj)}(k1,⋯,kn)∈(Z⩾1)n是所有的特征函数.
我们首先计算−△(j=1∏nsin(dπkj⋅xj))=d2π2∣k∣2j=1∏nsin(dπkj⋅xj).其次, 我们仍然有正交关系: 当 k=k′ 时, 我们就有∫(0,d)nj=1∏nsin(dπkj⋅xj)j=1∏nsin(dπkj′⋅xj)dx=0.对于 k∈(Z⩾1)n, 我们定义φk=j=1∏nsin(dπkj⋅xj),我们来证明 φk∈H01((0,d)n): 这显然是一个 H1((0,d)n) 中的函数, 下面用光滑的有紧支集的函数来逼近它: 由于对每个 j, sin(dπkj⋅x)∈H01((0,d)), 所以, 存在 {φj(p)(x)}p⩾1⊂C0∞((0,d)), 使得p→∞lim∥∥φj(p)(x)−sin(dπkj⋅x)∥∥H1=0.所以, ∥∥j=1∏nφj(p)(xj)−j=1∏nsin(dπkj⋅xj)∥∥H12=j=1∏n∥∥φj(p)−sin(dπkj⋅xj)∥∥H012→0.最终, 为了证明这是所有的特征函数, 我们来证明 {φk(x)}k∈(Z⩾1)n 构成 L2 的 Hilbert 基. 我们假设 n⩾2 (n=1 的情形已经完成) , 只要证明与 {φk}k∈(Z⩾1)n 都垂直的函数只有 0 即可 (这表明这些函数所张成的空间的闭包是整个 L2) . 这与证明高维的 Fourier 级数是一组基是完全一样的 (利用 Fubini 定理) , 请参考上学期 5 月 14 日的讲义.
我们现在研究 Ω=(0,d)n 上 Laplace 算子的特征值分布问题. 我们定义Ψ(λ)=∣{λk∣λk⩽λ}∣.按照上述计算, 我们有Ψ(λ)=∣∣{(k1,⋯,kn)∣∣ki⩾1,k12+⋯+kn2⩽π2d2λ}∣∣.这是在圆内的整点问题 (Gauss) : Ψ(λ) 是半径为 πdλ 的球在第一卦限中的整点的个数, 从而Ψ(λ)∼cnλ2n ⇔ λ→∞limcnλ2nΨ(λ)=1,其中 cn 是只依赖于维数的常数, 实际上, 通过计算半径为 πdλ 的球在第一卦限中的体积, 我们知道cn=(2π)n∣Bn(1)∣dn,其中 ∣Bn(1)∣=Γ(2n+1)π2n 是 Rn 中单位球的体积, 所以, Ψ(λ)∼(2π)nΓ(2n+1)dnλ2n.
我们考虑 λk 的渐近大小, 其中 k→∞. 我们在上式中取 λ=λk, 按照定义, Ψ(λk)=k, 从而k∼(2π)nΓ(2n+1)dnλk2n=(2π)nΓ(2n+1)∣Ω∣λk2n.所以, λk∼∣Bn(1)∣n2∣Ω∣n2(2π)2kn2, k→∞.我们将证明, 这个公式对于一般的区域 Ω 都成立, 这就是所谓的 Weyl 渐近公式.
如果 Ω 是连通的, 我们可以证明 λ1<λ2, 也就是说第一特征值的重数是 1.
对于上面的例子, 这一点很容易验证. 在这个例子中, 我们还看到, λ2=λ3=⋯=λn+1.
我们现在说明, 特征函数 φk(x) 是 Ω 上的光滑函数.
由于光滑性是局部性质, 所以, 我们只要证明对任意的 χ∈C0∞(Ω), 我们都有 χ(x)φk(x)∈C∞(Ω) 即可. 为了书写方便, 我们令 u=φk(x) 并且 λ=λk, 所以, −△u=λu.所以, −△(χ⋅u)=χλu−2∇χ⋅∇u−u△χ=∈L2(λχ−△χ)u−∈L22∇χ⋅∇u.按照定义, 由于 u∈H1(Ω), 所以右边都是 L2 的函数, 根据 △ 的正则性 (注意到, 由于 χ⋅u 的支集远离 ∂Ω, 我们不妨假设 ∂Ω 是光滑的即可) , 我们知道 χ⋅u∈H2. 所以, 根据 χ 选取的任意性, u 限制在任意的紧集 K⊂Ω 上都是 H2 的.
重复这个做法, 由于−△(χ⋅u)=∈H1(λχ−△χ)u−∈H12∇χ⋅∇u,所以, u 限制在任意的紧集 K⊂Ω 上都是 H3 的.
以此类推, 对任意的 χ∈C0∞(Ω), χ⋅u∈Hk(Rn). 根据 Sobolev 嵌入定理, 我们就知道 χ⋅u 是光滑的.
我们可以对 −△ 的特征值进行如下的变分表述:
假设 u∈H01(Ω) 并且 u=0, 我们令R(u)=∥u∥2⟨−△u,u⟩L2=形式上∥u∥2(−△u,u)L2.上面的定义之所以有意义是因为 −△u∈H−1(Ω) 而 u∈H01(Ω). 记 H=H01(Ω), 那么, 我们有
1) | 对于 λ1, 我们有λ1=u∈H,u=0infR(u). |
2) | 对于 k⩾2, 我们有λk=u=0u⊥φ1,⋯,u⊥φk−1,infR(u),其中, 上面的垂直关系是 (可以是) 用 L2-内积来定义的. |
3) | 令 Grk−1(H) 为 H 的所有 k−1 维线性子空间所构成的集合, 其中 k⩾1. 对于 P∈Grk−1(H), 令μ(P)=u⊥P,u=0infR(u).那么, λk=P∈Grk−1(H)supμ(P). |
4) | 对于 Q∈Grk(H), 定义ν(Q)=u∈Q,u=0supR(u).那么, λk=Q∈Grk(H)infν(Q). |
证明. 我们先证明 1) 和 2):
考虑 u=j⩾1∑αjφj∈L2(Ω). 由于 u⊥φ1,⋯,u⊥φk−1, 所以, u=j⩾k∑αjφj.从而, (根据定理之前的注记, ) 我们可以计算R(u)=∣αk∣2+∣αk+1∣2+⋯λk∣αk∣2+λk+1∣αk+1∣2+⋯⩾∣αk∣2+∣αk+1∣2+⋯λk∣αk∣2+λk∣αk+1∣2+⋯=λk.另外, R(φk)=λk, 这就给出了 1) 和 2).
为了证明 3), 我们先选取 Pk=φ1∧⋯∧φk−1, 这个符号代表的是由 φ1,⋯,φk−1 所张成的 k−1-维线性子空间. 根据上面的计算, 我们知道μ(Pk)=λk.所以, P∈Grk−1(H)supμ(P)⩾λk.为了证明反过来的不等式, 我们利用下面的观察:
• | 对每一个 P∈Grk−1(H), 总有不全为零的 α1,⋯,αk, 使得 α1φ1+⋯+αkφk⊥P. |
这是一道标准的线性代数习题: 为了让 α1φ1+⋯+αkφk⊥P, 我们假设 v1,⋯,vk−1 是 P 的一组基, 那么, 这个垂直的条件等价于⎩⎨⎧α1⋅(φ1,v1)+⋯+αk⋅(φk,v1)=0,⋯ ⋯,α1⋅(φ1,vk−1)+⋯+αk⋅(φk,vk−1)=0.这是 k 个未知数 k−1 个方程, 所以有解.
利用这个观察, 我们有μ(P)⩽R(α1φ1+⋯+αkφk)=∣α1∣2+⋯+∣αk∣2λ1∣α1∣2+⋯+λk∣αk∣2⩽∣α1∣2+⋯+∣αk∣2λk∣α1∣2+⋯+λk∣αk∣2=λk.所以, P∈Grk−1(H)supμ(P)⩽λk.这就证明了 3).
最后, 我们来证明 4). 通过选取 Qk=φ1∧⋯∧φk, 那么, 我们有ν(Qk)=α1,⋯,αk∈Csup(∑j⩽kαjφj,∑j⩽kαjφj)(−△(∑j⩽kαjφj),∑j⩽kαjφj)=α1,⋯,αk∈Csup∣α1∣2+⋯+∣αk∣2λ1∣α1∣2+⋯+λk∣αk∣2⩽α1,⋯,αk∈Csup∣α1∣2+⋯+∣αk∣2λk∣α1∣2+⋯+λk∣αk∣2=λk.又因为 φk 可以实现上面的最大值, 所以ν(Qk)=λk.从而, Q∈Grk(H)infν(Q)⩽λk.为了证明反过来的不等式, 我们利用下面的线性代数事实 (证明仿照前述) :
• | 对每一个 Q∈Grk(H), 总有不为零的 u∈Q, 使得 u⊥φ1,⋯,u⊥φk−1. |
据此以及 1) 和 2) 证明的过程, 我们有
μ(Q)⩾R(u)⩾λk.也就是说,
Q∈Grk(H)infν(Q)⩾λk.这就完成了证明.
根据 1), 我们知道对任意的 u∈H01(Ω), 有∥∇u∥L2(Ω)2⩾λ1∥u∥L2(Ω)2.这表明 λ1 是使得 Poincaré 不等式成立的最佳常数.
给定两个有界开区域 Ω1⊂Ω2⊂Rn, 对于每个 k⩾1, 我们都有λk(Ω1)⩾λk(Ω2).
证明. 通过将函数用零来延拓, 我们已经构造了自然的 (连续) 嵌入映射:
ι:H01(Ω1)→H01(Ω2).我们用
{φj}j⩾1 表示 Laplace 算子在小区域
Ω1 上的特征函数. 根据上面的定理, 我们有
λk(Ω1)=u=0u∈φ1∧⋯∧φksupR(u)=u=0u∈ι(φ1)∧⋯∧ι(φk)supR(u).其中, 第二个等号已经开始在
H01(Ω2) 中进行计算. 所以,
λk(Ω1)=u=0u∈ι(φ1)∧⋯∧ι(φk)supR(u)=ν(ι(φ1)∧⋯∧ι(φk))⩾λk(Ω2).最后一个等号利用的是
λk(Ω2) 在 4) 中的表述. 证明完毕.
对任意的有界开区域 Ω, 我们可以找到常数 c1 和 c2, 使得当 k→∞ 时, 我们有c1kn2⩽λk(Ω)⩽c2kn2, k→∞.
证明. 由于
Ω 是有界开区域, 所以我们总能找到
d,D>0, 使得
(0,d)n⊂Ω⊂(0,D)n (通过平行移动, 这不改变特征值) . 从而, 根据特征值的比较定理, 我们有
λk((0,D)n)⩽λk(Ω)⩽λk((0,d)n).根据我们之前的例子,
λk((0,D)n)∼∣Bn(1)∣n2D2(2π)2kn2, λk((0,d)n)∼∣Bn(1)∣n2d2(2π)2kn2,所以, 我们选取
c1=∣Bn(1)∣n2D2(2π)2, c2=∣Bn(1)∣n2d2(2π)2即可.