教材
《变分法及其应用——物理、力学、工程中的经典建模》欧斐君 高等教育出版社
有限维到无限维
向量中有有限个元素,它们可以进行加法、数乘、定义范数、定义内积、定义夹角。比如,对于向量aaa和bbb,其夹角余弦值为
cos<a,b>=(a,b)∣∣a∣∣⋅∣∣b∣∣ cos<a,b> = \frac{(a,b)}{||a|| \cdot ||b||} cos<a,b>=∣∣a∣∣⋅∣∣b∣∣(a,b)
那么,如果把一个向量中的元素数量扩展到无穷,向量就变成了函数,我们尝试仿照向量,对函数定义范数、内积、夹角。首先是内积,它变成了
(f,g)=∫f(x)g(x)dx (f,g) = \int f(x) g(x) dx (f,g)=∫f(x)g(x)dx
有了内积,很容易定义出范数。
∣∣f∣∣=(f,f)=∫f2(x)dx ||f|| = \sqrt{(f,f)} = \sqrt{\int f^2(x) dx} ∣∣f∣∣=(f,f)=∫f2(x)dx
以及夹角。
cos<f,g>=(f,g)∣∣f∣∣⋅∣∣g∣∣ cos<f,g> = \frac{(f,g)}{||f|| \cdot ||g||} cos<f,g>=∣∣f∣∣⋅∣∣g∣∣(f,g)
有了夹角,我们就可以定义函数的正交,即如果cos<f,g>=0cos<f,g>=0cos<f,g>=0,则函数fff与ggg正交。正交函数族的例子有三角函数(傅里叶分解的基函数),切比雪夫多项式等等。函数在一组两两正交函数ϕi(x)\phi_i(x)ϕi(x)上的投影,被称为广义傅里叶级数。即
f=∑−∞∞ciϕi f = \sum_{-\infty}^{\infty} c_i \phi_i f=−∞∑∞ciϕi
其中,
ci=(fi,ϕi) c_i = (f_i, \phi_i) ci=(fi,ϕi)
如果我们只取级数中的NNN项和fNf_NfN作为对fff的逼近,那么就有误差
e=∣∣fN−f∣∣ e = || f_N - f || e=∣∣fN−f∣∣
对于这NNN项的基函数,按(fi,ϕi)(f_i, \phi_i)(fi,ϕi)选取cic_ici能让误差最小,因此fNf_NfN是一个最佳逼近。
线性泛函的例子
就像向量函数以向量为自变量,泛函以函数(无限维的向量)为自变量。其映射JJJ为函数空间→R函数空间 \rightarrow R函数空间→R。
最速降线
问题:从(0,0)(0,0)(0,0)到(x1,y1)(x_1,y_1)(x1,y1)构建一个光滑斜曲面,使得小球在重力中下降速度最快。
假设曲面的函数为y=y(x)y=y(x)y=y(x),那么有两端的约束y(0)=0,y(x1)=y1y(0)=0,y(x_1)=y_1y(0)=0,y(x1)=y1;当小球下落到(x,y)(x,y)(x,y)时,按照能量守恒,它的速度为v=2gyv=\sqrt{2gy}v=2gy,取一个微元,斜面长度为1+y′2\sqrt{1+y'^2}1+y′2,因此通过这个微元的时间为
t=1+y′22gy t = \frac{\sqrt{1+y'^2}}{\sqrt{2gy}} t=2gy1+y′2
因此通过整个斜面的总时间为
T=∫0x1tdx=∫0x11+y′22gydx T = \int_0^{x_1} t dx = \int_0^{x_1} \frac{\sqrt{1+y'^2}}{\sqrt{2gy}} dx T=∫0x1tdx=∫0x12gy1+y′2dx
假设这个斜面是连续的,那么解一定存在于如下的函数空间(C1C^1C1表示一阶连续):
A={
y∣y(x)∈C1(x,y)y(0)=0,y(x1)=y1} A=\left\{ y | \begin{matrix} y(x) \in C^1(x,y) \\ y(0)=0, y(x_1)=y_1 \end{matrix} \right\} A={
y∣y(x)∈C1(x,y)y(0)=0,y(x1)=y1}
假设y∗y^*y∗是最后的解,那么我们要求的问题就变成了一个优化问题
y∗=argminy∈AT y^* = {\arg\min}_{y \in A} T y∗=argminy∈AT
一会儿再考虑如何求解这个方程。
极小曲面问题
问题:求一段通过(x1,y1)(x_1,y_1)(x1,y1)和(x2,y2)(x_2,y_2)(x2,y2)的曲线,绕着xxx轴旋转后侧面表面积最小。
首先写出面积公式。按xxx轴分成无数个微元,每个微元都是一个周长为2πy2\pi y2πy,高为1+y′2\sqrt{1+y'^2}1+y′2的长方形,因此总面积为
S=∫x1x22πy1+y′2dx S = \int_{x_1}^{x_2} 2 \pi y \sqrt{1+y'^2} dx S=∫x1