多元函数的驻点、Taylor 展开、Hessian 与二次型判别构成变分法的入门。它们训练的是同一种方法:选取允许的扰动,比较目标量在扰动前后的变化。已有计算基础可参考多元函数极值、梯度、海塞矩阵与正定二次型。
有限维的扰动展开
设 f 在内点 z∗∈Rn 的邻域内属于 C2。当 h→0 时,
f(z∗+h)=f(z∗)+∇f(z∗)Th+21hTHf(z∗)h+o(∥h∥2).
内部局部极小的必要条件为
∇f(z∗)=0,hTHf(z∗)h≥0对所有 h.
若驻点处 Hf(z∗) 正定,则得到严格局部极小;若半正定,二阶信息一般不足。三元及更高维问题需要判断整个二次型,不能只凭 Hessian 行列式的正负分类。
工程设计中的材料最少、表面积最小等问题,通常先把尺寸与约束写成函数关系。若能用约束消去一个变量,应同时记录剩余变量的定义域;内部驻点计算完成后,仍需检查边界、退化情况与全局比较。
把向量方向换成函数方向
对泛函
J[y]=∫abL(x,y,y′)dx,
选择容许方向 η,令 yε=y+εη。若 ε↦J[y+εη] 在零点邻域内二阶连续可微,就有沿固定方向的展开
J[y+εη]=J[y]+εδJ[y;η]+2ε2δ2J[y](η,η)+o(ε2).
| 有限维极值 | 变分法 |
|---|
| 未知向量 z | 未知函数 y |
| 向量扰动 h | 容许函数方向 η |
| Df(z)[h]=∇f(z)Th | 第一变分 δJ[y;η] |
| Hessian 二次型 hTHf(z)h | 第二变分 δ2J[y](η,η) |
| 一阶项对所有方向为零 | 第一变分对所有容许方向为零 |
固定端点要求 η(a)=η(b)=0。积分约束还会限制扰动方向,具体见等周约束变分问题。
二阶判别需要补上的条件
有限维中,正定矩阵在单位球面上具有严格正的最小值,所以其二次项能统一压过 Taylor 余项。无限维中的“每个非零方向都为正”不自动给出这样的统一下界;使用第二变分的充分条件时,必须同时明确范数、统一正下界及余项估计。
从这里进入第一变分与Euler-Lagrange方程,有限维极值的计算方法便成为变分法的学习起点。