变分法

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

变分法是处理泛函 (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$ 的积分：

1
2
3
4
5

    $$
    T[y] = \int_A^B dt = \int_{x_A}^{x_B} \frac{\sqrt{1+(y')^2}}{\sqrt{2gy}} \,dx
    $$

这就是我们要最小化的泛函。

1. 应用欧拉-拉格朗日方程：
   • 这里的拉格朗日量是 $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)$。

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

   • $\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}$ 。这是动量的变化率。

2. 组合成方程:

   $$\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$！

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

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

评论