参考文献
[Jost] | Jürgen Jost & Xianqing Li-Jost. Calculus of Variations. Cambridge Studies in Advanced Mathematics 64. Cambridge University Press, 1998. |
[Struwe] | Michael Struwe. Variational Methods. A Series of Modern Surveys in Mathematics 38. Springer, 2008. |
[Zhang] | 张恭庆. 非线性分析方法. 现代数学基础. 高等教育出版社, 2020. |
[Liu] | 刘宪高. 变分法与偏微分方程. 科学研究生教学丛书. 科学出版社, 2016. |
变分法入门
设 Ω⊂Rn 是有界区域, f∈C(Ω). 以下方程有解吗? {△uu∣∂Ω=f,=0.
为此, 考虑空间
X={u∈C1(Ω):u∣∂Ω=0},以及泛函
F:X→R,
F(u)=21∫Ω∣∇u∣2+∫Ωfu. 假设 v∈X 满足 F(v)=u∈XinfF(u), 则∫Ω(∇v⋅∇φ+fφ)=0,∀φ∈Cc∞(Ω).
证明. 首先注意到 XinfF>−∞, 毕竟 Poincaré 不等式说∫Ω∣∇u∣2≥C∫Ωu2,那么F(u)≥21∫Ω∣∇u∣2−2C∫Ωu2−2C1∫Ωf2≥−2C1∫Ωf2.(∗)更重要的是, dtd∣∣t=0F(v+tφ)=0 (不难验证 F(v+⋅φ)∈C1(R)) , 原因是F(v+tφ)≥F(v),∀t∈R.由此, 0=dtd∣∣t=0F(v+tφ)=∫Ω(∇v⋅∇φ+fφ).□
那么, 现在要关心的就是这样的
v 是否存在了. 对于一个
X 中的
极小化序列 {vn}n≥1, 即
n→∞limF(vn)=XinfF,{vn}n≥1 是否 (或在什么意义下) 会收敛于一个
v 呢?
v∈X 吗? 进一步,
F(v)=XinfF 成立吗?
回忆一下 Banach–Alaoglu 定理, 是说 Banach 空间的闭单位球是弱 * 紧的, 尤其是对于自反 Banach 空间, 闭单位球是弱紧的. 此定理能不能说明 {vn}n≥1 有子列弱收敛于 v∈X 呢? 这就要解决以下两个问题:
(i) X 是自反 Banach 空间吗?
明显不是, 但 X 关于 H1 范数的完备化是 Hilbert 空间, 因而是自反 Banach 空间, 记为 H01(Ω), 其中, ∥u∥H1:=∫Ω∣∇u∣2+∫Ωu2.
于是现在起改用 H01 为 X, vn 亦换成 H01(Ω) 下 F 的极小化序列。
不难看出, 命题 1.1 对 X=H01(Ω) 照样成立.
(ii) {vn}n≥1 关于 H1 范数有界吗?
确实如此, 这是因为改进上面的 (∗), 得到F(u)≥21∫Ω∣∇u∣2−4C∫Ωu2−C1∫Ωf2≥41∫Ω∣∇u∣2−C1∫Ωf2. 取 u 为 {vn}n≥1, 数列 {F(vn)}n≥1 有极限从而有上界, 所以 ∥∇vn∥L2 有界, 由 Poincaré 不等式, ∥vn∥H1 同样有界.
故 {vn}n≥1 有子列弱收敛于 v∈H01(Ω), 尽管目前还不知道是否 v∈X. 显然XinfF≤F(v)=n→∞limF(vn),那么只要说明 F(v)≤XinfF, 就有 F(v)=XinfF 了.
设 X 是 Banach 空间, F:X→R 是泛函.
• | 若 xn→x 必有 F(x)≤n→∞limF(xn), 则称 F 是下半连续的. |
• | 若 xn⇀x 必有 F(x)≤n→∞limF(xn), 则称 F 是弱下半连续的. (⇀ 是弱收敛) |
即是要说明 F 是弱下半连续的, 这是简单的: 由 Poincaré 不等式, 只需说明 ∥∇vn∥L2→∥∇v∥L2, 又由 Cauchy 不等式, n→∞lim∫Ω∣∇vn∣2≥∫Ω∣∇v∣2+2∫Ω∇v⋅∇vn−2∫Ω∣∇v∣2=∫Ω∣∇v∣2.
令人惊喜 (不爽) 的是, 本问题有简单得多的解法:
注意 φ↦−∫Ωfφ 是 H01(Ω) 上的有界线性泛函, 由 Riesz 表示定理, 存在 v∈H01(Ω) 使得∫Ω∇v⋅∇φ=−∫Ωfφ.
这里看上去未用变分法, 不过 Riesz 表示定理的证明本身即是求解一个变分问题.
设 Ω⊂Rn 是有界区域, g∈C∞(∂Ω). 以下方程有解吗? {△uu∣∂Ω=0,=g.
这次, 选取空间
X={u∈H1(Ω):u∣∂Ω=g},以及泛函
F:X→R,
F(u)=21∫Ω∣∇u∣2.再回忆 Mazur 定理, 是说 Banach 空间的闭凸集是弱闭的. 故对于凸泛函 F:X→R, 即F(λx+(1−λy))≤λF(x)+(1−λ)F(y)对 x,y∈X,λ∈(0,1),下半连续等价于弱下半连续. 这是因为, {F≤t} 是凸集对任意 t∈R, 那么由 Mazur 定理, {F≤t} 闭等价于弱闭, 而 F 下半连续就是 {F≤t} 是闭集对任意 t∈R, F 下半连续就是 {F≤t} 是弱闭集对任意 t∈R.
Sobolev 空间
对 u∈Lloc1(Ω), 称 Dαu(φ)=(−1)∣α∣∫uDαφ 对 φ∈Cc∞(Ω) 为 u 的广义导数.
对 k∈N,p∈[1,∞], (k,p) 阶 Sobolev 空间 Wk,p(Ω) 为Wk,p(Ω)={u∈Lp(Ω):Dαu∈Lp(Ω),∀∣α∣≤k}.W0k,p(Ω) 为 Cc∞(Ω) 在 ∥⋅∥k,p 范数下的完备化.
• | Wk,p(Ω)∩C∞(Ω) 在 Wk,p(Ω) 中稠密. |
• | 若 Ω 具有 C1 边界, 则 C∞(Ω) 在 Wk,p(Ω) 中稠密. |
W1−p1,p(∂Ω)=W1,p(Ω)/W01,p(∂Ω), 即∥{u+W01,p(Ω)}∥W1−p1,p(Ω)=v∈W01,p(Ω)inf∥u+v∥.
∂Ω∈C1, 存在有界线性算子 T:W1,p(Ω)→Lp(∂Ω), 使得
• | Tu=u∣∂Ω, u∈W1,p(Ω)∩C(Ω). |
• | ∥Tu∥Lp(∂Ω)≤C∥u∥W1,p(Ω). |
• | 对 φ∈C1(Rn,Rn), 有∫Ωudivφ=−∫ΩDu⋅φ+∫∂Ω(φ⋅n)Tu. |
∂Ω∈C1, 1<p<∞, 存在有界线性映射 ext:W1−p1,p(∂Ω)→W1,p(Ω) 使得 T∘ext=id.
Banach 空间上的微分
本讲基本来自于 [Zhang]. 以下 L(X,Y) 是 X 到 Y 的连续线性算子全体.
设 X,Y 是个 Banach 空间, U 是 X 的开子集, x0∈U, f:U→Y, f 在 h 方向可导是指存在 lh∈Y 使∥f(x0+th)−f(x0)−tlh∥Y=o(t),t→0,记为 δf(x0,h)=lh, 即δf(x0,h)=dtd∣∣t=0f(x0+th).
记号同上, f 在 x0 可微是指存在 A(x0)∈L(X,Y) 使∥f(x)−f(x0)−A(x0)(x−x0)∥Y=o(∥x−x0∥X),x→x0,记为 Df(x0)=A(x0). 此时对任意 h∈X, δf(x0,h)=Df(x0)h.
记号同上, 若 Df:U→L(X,Y) 是连续的 (L(X,Y) 用算子范数拓扑) , 则称 f 是 C1 映射.
设 f:U→Y G-可导, 对于任意 x∈U, 存在 A(x)∈L(X,Y), 使得 δf(x0,h)=A(x0)h 对一切 h∈X, 且x∈U↦A(x)∈L(X,Y)是连续的, 则 f 在 x0 处 F-可导, 且 Df(x0)=A(x0).
• | f:U→R, x0∈U⊂X, f 在 x0 可微, Df(x0)∈X∗. |
• | 当 X 是 Hilbert 空间, 存在 y0∈X 使 Df(x0)h=⟨y0,h⟩. |
• | 当 X 是 Hilbert 空间, f:X→R,x↦∥x∥2, 由于f(x0+h)−f(x0)=∥x0+h∥2−∥x0∥2=2⟨x0,h⟩+∥h∥2=2⟨x0,h⟩+o(∥h∥),有 A(x0)=2x0, 则 x↦A(x) 连续. |
• | Ω 是有界集, 对 u∈H01(Ω), 令f(u)=21∫Ω∣∇u∣2−∫Ω∣u∣p.根据 Sobolev 嵌入 H01(Ω)↪L2∗(Ω), 在 p≤2∗=n−22n 时 f(u) 是有意义的. 考察后一项是不是 C1 的, 记A(u):=dtd∣∣t=0∫Ω∣u+th∣p=p∫Ω∣u∣p−2u∈(H01(Ω))∗:=H−1(Ω),当 ui→u 于 H01(Ω) 时, 由于 ∥h∥Lp≲∥h∥H01, 有 (最后一步用了定理 3.7) ∥A(ui)−A(u)∥H−1(Ω)=h∈H01(Ω)∥h∥≤1supp∣∣∫Ω(∣ui∣p−2ui−∣u∣p−2u)h∣∣≲∥∣ui∣p−2ui−∣u∣p−2u∥Lp−1p→0. |
设 Ω 是测度空间, 称 φ:Ω×RN→R,(x,ξ)↦φ(x,ξ) 是 Carathéodory 函数, 是说
• | 对几乎处处的 x∈Ω, φ(x,⋅) 连续, |
• | 对所有的 ξ∈RN, φ(⋅,ξ) 可测. |
设 p1,p2≥1, a>0,b∈Lp2(Ω). 设 φ 是 Carathéodory 函数, 满足∣φ(x,ξ)∣≤b(x)+a∣ξ∣p2p1,则f:Lp1(Ω,RN)u(x)→Lp2(Ω,RN),↦φ(x,u(x))是有界连续函数.
设 Ω⊂Rn, φ:Ω×R→R 是 Carathéodory 函数, 对几乎处处的 x, φ(x,ξ)∈Cξ1, 且偏导数 φξ:Ω×R→R 亦是 Carathéodory 函数, 而且∣φ(x,ξ)∣≤b1(x)+a1∣ξ∣n−22n,∣φξ(x,ξ)∣≤b2(x)+a2∣ξ∣n−2n+2,b1∈L1(Ω),b2∈Ln+22n(Ω),那么 H1(Ω) 上的泛函f(u)=∫Ωφ(x,u(x))是 C1 的, 且δf(u,v)=∫Ωφξ(x,u(x))⋅v(x).
设 X,Y,Z 是 Banach 空间, U 是 X×Y 的开子集, f∈C1(U,Z). 若对 (x0,y0)∈U 有f(x0,y0)=0,(fy(x0,y0))−1∈L(Z,Y),则存在 r0,r1>0 及唯一的 u∈C(Br0(x0),Br1(y0)) 使得 Br0(x0)×Br1(y0)⊂U, 且u(x0)=y0,f(x,u(x))=0 对 x∈Br0(x0),此时Du(x)=−(fy(x,u(x)))−1∘fx(x,u(x))对 x∈Br0(x0).
Lagrange 乘子法
设 X 是 Banach 空间, 泛函 f,g:X→R 是 C1 的, 若 M={x:g(x)=0} 是 C1 子流形, 则在 f 限制于 M 的临界点 x0, 存在 λ∈R 使得 Df(x0)=λDg(x0).
M 在 x0 处的切空间是 X 的子空间Tx0M={h∈X:存在 M 中的 C1 曲线 γ(t),γ(0)=x0,γ′(0)=h}.
若 Dg(x0)=0, 则 Tx0M=kerDg(x0).
证明. 根据定义, 对任何的 h∈Tx0M, 存在 M 中的 C1 曲线使得 g(γ(t))=0. 那么 0=dtd∣∣t=0g(γ(t))=Dg(x0)h.
由
Dg(x0)=0, 有
ImDg(x0)=R, 则
X=kerDg(x0)⊕Rw, 即对任意
x∈X,
x=h1+(Dg(x0)x)w.找曲线
γ(t)=x0+th+α(t)w,
α(t) 待定, 要求
α′(0)=0.
G(t,α)=g(x0+th+α(t)w)=0.由于
G(0,0)=0,DαG(0,0)=Dg(x0)w=1,得
(0,0) 有邻域
Br0(0)×Br1(0), 对任意
t∈Br0, 存在
α 使得
G(t,α(t))=0.
设 f,g1,⋯,gN:X→R 是 C1 泛函, 则存在 λ1,⋯,λN∈R, 使Df(x0)=i=1∑NλiDgi(x0).
证明. 只证明
N=1 的情况. 由于
0=dtd∣∣t=0f(γ(t))=Df(x0)h,有
kerDg(x0)⊆kerDf(x0), 可证存在
Df(x0)=λDg(x0).
X=H01(Ω), {−△uu∣∂Ω=λu,=0.考虑泛函 f,g:H01(Ω)→R, f(u)=21∫Ω∣∇u∣2,g(u)=∫Ωu2−1.设 u0 是 f 在 M={g=0} 上的最小值点, 则 Df(u0,v)=λDg(u0,v), 即∫∇u0⋅∇v=λ∫u0v.由于 f(∣u0∣)≤f(u0), ∣u0∣ 也是最小值点, 故 u0≥0, 又由次调和函数的强极值原理, u0>0.
X=W01,p(Ω), 1<p<∞, 定义 p-Laplace 算子△pu=∇(∣∇u∣p−2∇u){−△puu∣∂Ω=λ∣u∣p−2u,=0.这个情况下变分法与 p=2 差别不大, 而特征值非常复杂, 未必是可数个.
不过 p=1 时 W01,1(Ω) 甚至不是自反的, u∈W01,1(Ω)\{0}min∫Ω∣u∣∫Ω∣∇u∣都不存在, 为了有最小值要在更大的有界变差空间 BV(Ω) 中做.
{−△uu∣∂Ω=∣u∣p−2u,=0.2<p<2∗ 时, 考虑泛函f(u)=21∫Ω∣∇u∣p−p1∫Ω∣u∣p,δf(u,v)=∫Ω∇u⋅∇v−∫Ω∣u∣p−2uv.p>2∗ 时, 若对 x∈∂Ω, x⋅n(x)>0, 方程没有正解. 这是因为, 两边乘以 x⋅Du 再积分, 得∫Ω(△u)x⋅Du=(2n−1)∫Ω∣∇u∣2−21∫∂Ωxj∣∇u∣2nj+∫∂Ω∂iu∂juxinj,∫Ωup−1x⋅Du=−pn∫Ω∣u∣p.其中∫∂Ω∂iu∂juxinj=∫∂Ω∣∇u∣2x⋅n(x),于是(2n−1−pn)∫Ω∣∇u∣2=−21∫∂Ω∣∇u∣2x⋅n≤0,而 2n−1−pn>0, 矛盾.
Noehther 定理
记 v(t)=x˙(t), 定义动量 p=DvL(v,x), 这里 x,v,p∈Rn 对任何 t. 例如 L=21m∣v∣2−H(x), 就有 p=mv.
设变换 {p=DvL(v,x)x=x 可逆, 记为 {v=v(p,x)x=x, 并定义H=p⋅v−L(v,x).
{x˙(t)=∂p∂Hp˙(t)=−∂x∂H, 且 H(p(t),v(t))= 常数.
证明. 由 Euler–Lagrange 方程, ⎩⎨⎧HpiHxi=vi+p⋅vpi−DvL⋅vpi=vi=x˙i,=p⋅vxi−DvL⋅vxi−DxiL=−DxiL=−dtdDviL=−p˙i.由此,dtdH=Hp⋅p˙+Hx⋅x˙=0.□
现有区域变分ητ:Rn→Rn,η0(x)=x,∂τητ∣∣τ=0=X,与函数变分uτ:Rn→Rn,u0(x)=u(x),∂τuτ∣∣τ=0=h.记FΩ(u)=∫ΩL(Du,u,x)dx.
设 u 是 FΩ 的临界点, FΩ 是关于 ητ,uτ 不变的, 则 div(Dpv⋅h−X)=0 当且仅当 ∑i∂xi(hDpiL−LXi)=0.
证明. 见 Evans 8.6.2 Theorem 1 (P513).
u 是调和函数, ητ:x→x+τek,τ∈(−ε,ε), uτ=u(x+τek), 取FΩ(u)=∫Ω∣∇u∣2.有FΩ(uτ)=∫Ωτ∣∇u∣2.此时0=div(Dpv⋅h−X)=i∑∂xi(uxiDuxk−21∣∇u∣2δik),即对任意 Ω, ∫∂Ω(uxiuxk−∣∇u∣2δik)ni=0.
∂ttu=△u,(t,x)∈R×Rn,取F(u)=∫(0,T)×Ω21(ut2−∣∇u∣2.)详见 Evans 8.6.2 Example 3 (P515-6).
设 u:M→N 是流形间的映射, u 为调和映照是说 Euler–Lagrange 方程成立. 如u:B1⊂Rn→Sn⊂Rn+1,x↦(u1(x),⋯,un+1(x)),则 i=1∑n+1ui(x)2=1. 令A={u∈H1(B1,Rm+1),u∣∂B1=g,u(x)∈Sn a.e. x∈B1.}{△u=−u∣∇u∣2,u∣∂B1=g.用 Noether 定理, 由球面的旋转对称性, uτ(x)=Aτu(x),A∈SOn+1, 守恒律为div(uj∇ui−ui∇uj)=0,∀i,j=1,⋯,n+1.
集中紧原理
参考 [Liu] 3.3 节.
回忆有嵌入 H1(Rn)↪Lp∗(Rn),2≤p≤2∗=n−22n. 当 2<p<2∗ 时, 考虑 Rn 上的方程−△u+u=∣u∣p−2u,(∗)其中 u∈H1(Rn). 假如泛函F(u)=21∫∣∇u∣2+∣u∣2−p1∫∣u∣p有最小值 u0, 则 u0 满足方程 (∗); 然而它并无下界. 转而取泛函F(u)=21∫∣∇u∣2+∣u∣2在 M={u∈H1(Rn):∥u∥Lp=1} 上的最小值, 若存在, 由 Lagrange 乘子法有−△u+u=λ∣u∣p−2u,取 v=au, 其中 λa/ap−1=1, 则 v 是 (∗) 的解了. 然而用极小化序列找 u 遇到的困难是, H1(Rn)↪Lp(Rn),2<p<2∗ 不是紧的: 有两种方法使得 H1(Rn) 的有界集不是 Lp(Rn) 的紧集, 一种是平移至无穷远处, 另一种是 u(x) 变成 uλ(x)=λnu(λx),λ→∞ (称为集中紧) . 下面的定理把这严格化.
设 {μm}m≥1 是一列 Rn 上的概率测度, 则存在子列, 仍记作 {μm}m≥1, 满足以下性质:
紧性 | 存在 Rn 的点列 {xm}, 使得 ∀ε>0, 存在 R>0 使得 μm(BR(xm))≥1−ε. |
消失 | 对任何 R>0, n→∞limx∈Rnsupμm(BR(x))=0. |
二分 | ∀ε>0, 存在 R>0 与点列 {xm} 和的数列 {Rm}→∞, 使得m→∞lim∣μm(BR(xm))−λ∣+∣μm(Rn\BRm(xm))−(1−λ)∣≤ε. |
证明. 见 [
Liu] 定理 3.17 (P112–5) .
回到前面的问题. 在定理中令
μm=∣um∣pdx, 则这三个性质之一成立. 假设 “紧性” 成立, 则对于极小化序列
vm(x)=um(x+xm), 有
∫BR∣vm∣pdx≥1−ε.由于
H1(Rn) 中
vm⇀v, 用紧嵌入
H1(BR)↪Lp(BR), 得
Lp(BR) 中有子列
vm→v, 则
∥v∥Lp(BR)p≥1−ε, 令
ε→0 得
v∈M. 下面要排除 “消失” 和 “二分” 的情况.
设 {um} 是 H1 的有界序列, 存在 R>0 使 m→∞lim∥um∥L2(BR)=0, 则存在子列 um→0 于 Lp(Rn) 对任意 p∈(2,2∗).
假设是消失的情况, 现在引理说明
∥um∥Lp→0, 则
∥u∥Lp=0,
u∈/M, 矛盾. 这排除了 “消失”.
假设是二分的情况, 构造 BRm 上的紧支函数 φm, 在 BRm/2 上为 1. (...) 这排除了 “二分”.