变分法

想象一下，你要求一个函数的最小值，你会求导并令其为零．但如果问题是：在所有可能的函数（或路径）中，哪一个能使某个量（比如时间或能量）最小？ 这就是变分法要解决的问题．

变分法是处理 泛函(Functional) 的数学领域．泛函可以看作是「函数的函数」，它的输入是一个函数，输出一个数值．

• 函数: 输入一个数 $x$，输出一个数 $f(x)$．
• 泛函: 输入一个函数 $y(x)$，输出一个数 $J[y]$．

变分法的核心问题是：找到那个能使泛函 $J[y]$ 取极值（极大值或极小值）的特定函数 $y(x)$．

欧拉 - 拉格朗日方程

这是变分法中最核心的工具．

考虑一个由积分定义的泛函 $J[y]$：

$$J[y]={\int }_a^{b}L(x,y(x),y^{\prime }(x))dx$$

其中 $y(x)$ 是我们想找的函数，$y^{\prime }(x)$ 是它的导数，$L$ 是一个给定的函数，通常在物理中被称为 拉格朗日量．我们的目标是找到一个函数 $y(x)$，它满足边界条件 $y(a)=y_a$ 和 $y(b)=y_b$，并且能使泛函 $J[y]$ 取极值．

这个函数 $y(x)$ 必须满足下面的 欧拉 - 拉格朗日方程：

$$\frac{\partial L}{\partial y}-\frac{d}{dx}\left(\frac{\partial L}{\partial y^{\prime }}\right)=0$$

这个方程看起来可能有点吓人，但它只是普通函数求极值（令导数为零）在泛函上的推广．

例题 1：最速降线问题 (Brachistochrone Problem)

问题：在重力场中，一个珠子从点 A 沿着一根无摩擦的曲线轨道滑到点 B．我们应该如何设计轨道的形状，才能使珠子下滑的时间最短？

这是一个经典的变分问题．

1. 建立泛函：我们需要最小化的是总时间 $T$．时间是路径的泛函 $T[y]$．

   • 在轨道上任意一点，珠子的速度 $v$ 由能量守恒决定：$\frac{1}{2}m{v}^2=mgy$，所以 $v=\sqrt{2gy}$（假设从 $y=0$ 开始）．
   • 一小段弧长 $ds=\sqrt{d{x}^2+d{y}^2}=\sqrt{1+(y^{\prime }{)}^2}dx$．
   • 走过这段弧长的时间 $dt=\frac{ds}{v}=\frac{\sqrt{1+(y^{\prime }{)}^2}}{\sqrt{2gy}}dx$．

   • 总时间就是对 $dt$ 的积分：

     $$ T[y] = \int_A^B dt = \int_{x_A}^{x_B} \dfrac{\sqrt{1+(y')^2}}{\sqrt{2gy}} \,dx $$ 这就是我们要最小化的泛函．

2. 应用欧拉 - 拉格朗日方程：

   • 这里的拉格朗日量是 $L(y,y^{\prime })=\frac{\sqrt{1+(y^{\prime }{)}^2}}{\sqrt{y}}$（忽略常数）．
   • 将其代入欧拉 - 拉格朗日方程进行求解（过程比较复杂），最终会发现解出的曲线是一条 摆线 (Cycloid)．

这说明，连接两点最快的路径并不是直线，而是一段倒置的摆线．

物理中的应用：最小作用量原理

变分原理是理论物理的基石，它用一种极其深刻和优美的方式统一了物理学的各个分支．这个核心思想就是 最小作用量原理 (Principle of Least Action)．

该原理指出：一个物理系统从一个状态到另一个状态的真实演化路径，是使其「作用量」$S$ 取极值的路径．

作用量 $S$ 本身是一个泛函，定义为拉格朗日量 $L$ 对时间的积分：

$$S[q]={\int }_{t_1}^{t_2}L(q,\dot{q},t)dt$$

• $q$ 代表系统的广义坐标（可以是位置、角度等）．
• $\dot{q}$ 是广义速度．
• 拉格朗日量 $L=T-U$，即系统的总动能减去总势能．

从最小作用量到牛顿第二定律

让我们看看这个优美的原理是如何导出我们熟悉的牛顿定律的．

问题：对于一个在势能 $U(x)$ 中运动的一维粒子，其拉格朗日量为 $L=T-U=\frac{1}{2}m{\dot{x}}^2-U(x)$．请用欧拉 - 拉格朗日方程找出其运动方程．

解：

1. 确定变量和拉格朗日量:

   • 自变量是时间 $t$．
   • 函数是位置 $x(t)$．
   • 导数是速度 $\dot{x}(t)$．
   • $L(x,\dot{x})=\frac{1}{2}m{\dot{x}}^2-U(x)$．

2. 计算欧拉 - 拉格朗日方程的各项:

   • $\frac{\partial L}{\partial x}=\frac{\partial }{\partial x}(\frac{1}{2}m{\dot{x}}^2-U(x))=-\frac{dU(x)}{dx}$．我们知道力 $F=-dU/dx$．所以这一项就是力 $F$．
   • $\frac{\partial L}{\partial \dot{x}}=\frac{\partial }{\partial \dot{x}}(\frac{1}{2}m{\dot{x}}^2-U(x))=m\dot{x}$．这正是动量 $p$．
   • $\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{x}}\right)=\frac{d}{dt}(m\dot{x})=m\ddot{x}$．这是动量的变化率．

3. 组合成方程:

   $$\frac{\partial L}{\partial x}-\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{x}}\right)=0$$ $$F-m\ddot{x}=0⟹F=m\ddot{x}$$

   这正是牛顿第二定律 $F=ma$！

我们从一个抽象的、关于整个路径的积分（作用量）的极值原理出发，竟然推导出了描述每一瞬间的、局部的运动定律．这展示了变分法的强大威力．拉格朗日力学和哈密顿力学正是建立在这一原理之上，为描述复杂系统提供了比牛顿力学更普适和强大的框架．

本页面最近更新：2026/7/3 22:14:58，更新历史
发现错误？想一起完善？ 在 GitHub 上编辑此页！
本页面贡献者：Leafuke
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供，附加条款亦可能应用

评论