駐定作用的一體多面：從變分抵消、局部動力到極限拮抗表徵

副標題：簡諧振子上的解析推導、離散變分與第二變分基本驗證
英文標題：The Many Faces of Stationary Action: From Variational Cancellation and Local Dynamics to an Extremal-Antagonistic Representation

作者：Neo.K（許筌崴）
機構：EveMissLab／一言諾科技有限公司
版本：v1.0
日期：2026 年 7 月 14 日
文件性質：理論物理方法論／變分原理／數值驗證／極限拮抗原理橋接論文
前置文件：《極限拮抗原理：邊界逼近、動態抵消、歷史連接與持續差異》v2.0
驗證範圍：一維簡諧振子、固定端點、離散作用量、第二變分與共軛時間
命題強度：解析等價、模型內數值驗證與橋接猜想；不宣稱已證明所有作用量原理皆源於極限拮抗

摘要

「最小作用量原理」常被理解為自然在所有可能路徑中選取作用量最小者。然而，現代變分力學更一般的條件是：

\boxed{ \delta S=0 }

即物理路徑使作用量的一階變分消失。該路徑可以是局部最小、局部最大、鞍點或退化駐定點；究竟屬於何者，必須由第二變分或作用量 Hessian 判斷。

本文提出「駐定作用的一體多面」觀點：作用量駐定、Euler–Lagrange 局部運動方程，以及不同作用量分項之間的一階變分抵消，不應被視為彼此競爭的三種理論，而是同一動力結構在路徑空間、局部時間與分項生成視角下的不同表徵。

本文以一維簡諧振子：

L(q,\dot q) = \frac12m\dot q^2 - \frac12kq^2

為最小驗證案例。首先解析推導：

\delta S=0 \Longleftrightarrow m\ddot q+kq=0.

其次將作用量離散化，直接求解離散駐定方程，並與解析 Euler–Lagrange 路徑比較。在 $N=240$ 個時間區間下，短區間 $T=2$ 的最大路徑誤差約為：

4.19\times10^{-6},

長區間 $T=4$ 的最大路徑誤差約為：

6.29\times10^{-5}.

離散作用量梯度與局部 Euler–Lagrange 殘差均接近數值零：

\max|\nabla S| \approx 1.16\times10^{-13},

\max|m\ddot q+kq| \approx 1.39\times10^{-11}.

本文進一步將作用量分解為動能作用與位能作用：

S=S_K+S_V.

在駐定路徑上，數值驗證：

\nabla S_K \approx -\nabla S_V,

顯示「作用量駐定」確可被觀察為不同變分貢獻的一階抵消。

第二變分分析則顯示，在第一個共軛時間：

T_c=\frac{\pi}{\omega}

之前，作用量 Hessian 正定，經典路徑是局部最小；通過該時間後，Hessian 出現負特徵值，經典路徑仍滿足 $\delta S=0$ ，卻已成為鞍點。於 $T=2$ 時：

\lambda_{\min}\approx0.01223>0,

而於 $T=4$ 時：

\lambda_{\min}\approx-0.006386<0.

由此可確認：

\boxed{ \text{最小作用只是駐定作用的一個條件性子類。} }

本文最後提出「變分—局部—拮抗三重表徵命題」：對可分解作用量，物理路徑可以同時被描述為作用量的一階駐定、局部 Euler–Lagrange 動力，以及分項變分的相互抵消。然而，這只能建立駐定作用與拮抗式表徵之間的結構橋梁，不能單靠簡諧振子證明所有作用量皆源於《極限拮抗原理》中的高強度邊界逼近、歷史連接與餘量結構。

本文的核心結論是：

\boxed{ \text{「最小」「駐定」「抵消」並非必然互斥，} }

\boxed{ \text{它們可能是同一動力系統在不同觀察算子下呈現的不同面。} }

關鍵詞： 最小作用量、駐定作用、Hamilton 原理、Euler–Lagrange 方程、第二變分、共軛點、鞍點、簡諧振子、變分抵消、一體多面、極限拮抗

1. 問題提出：自然究竟是在「最小化」，還是在「駐定」

1.1 歷史名稱造成的直覺偏差

「最小作用量原理」容易導致三種直覺：

自然會比較所有路徑；
自然會選擇數值最小的作用量；
所有經典運動都可理解為某種節省、經濟或最低消耗。

這些直覺在部分問題中具有啟發性，但不構成一般數學定義。

對路徑 $q(t)$ ，Hamilton 作用量為：

S[q] = \int_{t_1}^{t_2} L(q,\dot q,t)\,dt.

物理路徑滿足的是：

\boxed{ \delta S=0. }

這是一階變分條件。它只表示：在固定端點的允許路徑變分下，物理路徑附近的作用量沒有一階改變。

它不直接表示：

S[q_{\mathrm{cl}}] < S[q] \qquad \forall q\neq q_{\mathrm{cl}}.

1.2 「駐定」的最低數學含義

取：

q_\varepsilon(t) = q(t)+\varepsilon\eta(t),

其中：

\eta(t_1)=\eta(t_2)=0.

作用量展開：

S[q_\varepsilon] = S[q] + \varepsilon\delta S[q;\eta] + \frac{\varepsilon^2}{2} \delta^2S[q;\eta] + O(\varepsilon^3).

駐定要求：

\delta S[q;\eta]=0

對所有允許的 $\eta$ 成立。

若：

\delta^2S>0

對所有非零變分成立，則是局部最小。

若：

\delta^2S<0,

則是局部最大。

若第二變分有正方向也有負方向，則為鞍點。

所以：

\boxed{ \text{駐定是一階條件，最小是二階分類。} }

2. 一體多面：同一動力結構的三個觀察空間

本文將同一物理系統記為：

\mathfrak D.

對它可施加不同觀察算子。

2.1 路徑空間觀察

\mathcal O_{\mathrm{path}}(\mathfrak D) = \delta S=0.

這個視角一次處理完整歷史 $q:[t_1,t_2]\to Q$ 。

2.2 局部時間觀察

\mathcal O_{\mathrm{local}}(\mathfrak D) = \frac{d}{dt} \left( \frac{\partial L}{\partial\dot q} \right) - \frac{\partial L}{\partial q} = 0.

這個視角將完整路徑改寫為每一時刻成立的局部微分方程。

2.3 分項變分觀察

若：

S=S_1+S_2+\cdots+S_n,

則駐定條件為：

\delta S_1+\delta S_2+\cdots+\delta S_n=0.

此時同一駐定事件可被觀察為不同作用量貢獻的一階抵消。

2.4 不是任何分解都具有相同本體意義

作用量可以被多種方式分解。例如：

L=T-V,

也可以加上一個全導數：

L' = L+\frac{dF}{dt},

而不改變固定端點下的 Euler–Lagrange 方程。

所以：

\boxed{ \text{數學分項抵消} \not\Rightarrow \text{每一分項都必然是獨立基本實體。} }

「一體多面」需要轉換關係與觀測等價，不能只靠語言相似。

3. 從駐定作用到 Euler–Lagrange 方程

考慮：

S[q] = \int_{t_1}^{t_2} L(q,\dot q,t)\,dt.

其一階變分：

\delta S = \int_{t_1}^{t_2} \left[ \frac{\partial L}{\partial q}\eta + \frac{\partial L}{\partial\dot q}\dot\eta \right]dt.

對第二項分部積分：

\int_{t_1}^{t_2} \frac{\partial L}{\partial\dot q} \dot\eta\,dt = \left[ \frac{\partial L}{\partial\dot q}\eta \right]_{t_1}^{t_2} - \int_{t_1}^{t_2} \frac{d}{dt} \left( \frac{\partial L}{\partial\dot q} \right)\eta\,dt.

固定端點使邊界項消失，得到：

\delta S = \int_{t_1}^{t_2} \left[ \frac{\partial L}{\partial q} - \frac{d}{dt} \left( \frac{\partial L}{\partial\dot q} \right) \right] \eta(t)\,dt.

由於 $\eta$ 任意，故：

\boxed{ \frac{d}{dt} \left( \frac{\partial L}{\partial\dot q} \right) - \frac{\partial L}{\partial q} = 0. }

這顯示：

\boxed{ \mathcal O_{\mathrm{path}}(\mathfrak D) \Longleftrightarrow \mathcal O_{\mathrm{local}}(\mathfrak D) }

在正則變分條件下是解析等價，而非比喻相似。

4. 簡諧振子案例

4.1 拉格朗日量

令：

L = \frac12m\dot q^2 - \frac12kq^2,

並定義：

\omega^2=\frac{k}{m}.

Euler–Lagrange 方程：

m\ddot q+kq=0.

固定：

q(0)=q_0, \qquad q(T)=q_T.

在：

\sin(\omega T)\neq0

時，解析經典路徑為：

q_{\mathrm{cl}}(t) = q_0\cos(\omega t) + \frac{ q_T-q_0\cos(\omega T) }{ \sin(\omega T) } \sin(\omega t).

本文數值設定：

m=1, \qquad k=1, \qquad \omega=1,

q_0=0, \qquad q_T=1.

因此：

q_{\mathrm{cl}}(t) = \frac{\sin t}{\sin T}.

5. 離散作用量與局部方程

5.1 離散作用量

將區間分成 $N$ 段：

t_n=n\Delta t, \qquad \Delta t=\frac{T}{N}.

採梯形位能離散：

\boxed{ S_d[q] = \sum_{n=0}^{N-1} \left[ \frac{m}{2} \left( \frac{q_{n+1}-q_n}{\Delta t} \right)^2 - \frac{k}{4} \left( q_n^2+q_{n+1}^2 \right) \right]\Delta t. }

對內部節點 $q_n$ 求偏導：

\frac{\partial S_d}{\partial q_n} = \frac{m}{\Delta t} \left( 2q_n-q_{n-1}-q_{n+1} \right) - k\Delta t\,q_n.

駐定要求：

\frac{\partial S_d}{\partial q_n}=0,

等價於：

\boxed{ m \frac{ q_{n+1}-2q_n+q_{n-1} }{ \Delta t^2 } + kq_n = 0. }

這正是：

m\ddot q+kq=0

的中心差分形式。

所以離散模型再次顯示：

\boxed{ \text{離散路徑駐定} = \text{離散局部運動方程}. }

6. 數值驗證一：駐定路徑與局部解一致

本文以：

N=240

進行計算。

6.1 短時間區間

取：

T=2<\pi.

離散駐定路徑與解析路徑的最大誤差為：

\boxed{ \max_n |q_n-q_{\mathrm{cl}}(t_n)| = 4.19\times10^{-6}. }

6.2 長時間區間

取：

T=4>\pi.

最大誤差為：

\boxed{ 6.29\times10^{-5}. }

誤差主要來自有限網格離散。

6.3 作用量梯度

兩個案例的離散作用量梯度最大值約為：

1.07\times10^{-13} \sim 1.16\times10^{-13}.

因此：

\boxed{ \nabla S_d\approx0. }

6.4 局部 Euler–Lagrange 殘差

離散局部方程殘差：

r_n = m \frac{ q_{n+1}-2q_n+q_{n-1} }{ \Delta t^2 } + kq_n.

最大殘差約為：

6.42\times10^{-12} \sim 1.39\times10^{-11}.

所以：

\boxed{ \text{同一條數值路徑同時滿足全域駐定與局部動力。} }

7. 數值驗證二：動能與位能變分抵消

將作用量拆為：

S=S_K+S_V,

其中：

S_K = \int_0^T \frac12m\dot q^2\,dt,

S_V = -\int_0^T \frac12kq^2\,dt.

離散梯度分別為：

(\nabla S_K)_n = \frac{m}{\Delta t} \left( 2q_n-q_{n-1}-q_{n+1} \right),

(\nabla S_V)_n = -k\Delta t\,q_n.

在駐定路徑上：

\boxed{ \nabla S_K+\nabla S_V=0. }

數值最大剩餘：

\boxed{ \max_n | (\nabla S_K)_n+(\nabla S_V)_n | \approx 1.16\times10^{-13}. }

因此同一件事可以被表述為：

路徑空間

\delta S=0.

分項空間

\delta S_K=-\delta S_V.

局部時間

m\ddot q=-kq.

這是本文對「一體多面」最直接的模型內驗證。

8. 第二變分：為什麼駐定不等於最小

8.1 第二變分

簡諧振子作用量的第二變分為：

\boxed{ \delta^2S[\eta] = \int_0^T \left[ m\dot\eta^2 - k\eta^2 \right]dt. }

令：

\eta_n(t) = \sin \left( \frac{n\pi t}{T} \right).

則：

\delta^2S[\eta_n] = \frac{T}{2} \left[ m \left( \frac{n\pi}{T} \right)^2 - k \right].

用：

k=m\omega^2

可寫成：

\delta^2S[\eta_n] = \frac{mT}{2} \left[ \left( \frac{n\pi}{T} \right)^2 - \omega^2 \right].

最低模態 $n=1$ 的符號在：

T=\frac{\pi}{\omega}

改變。

8.2 共軛時間以前

若：

T<\frac{\pi}{\omega},

則所有 $n\geq1$ 都有：

\delta^2S[\eta_n]>0.

所以經典路徑是嚴格局部最小。

8.3 共軛時間

若：

T=\frac{\pi}{\omega},

則最低模態：

\delta^2S[\eta_1]=0.

Hessian 出現零模態，表示局部唯一性或最小性條件退化。

8.4 共軛時間以後

若：

T>\frac{\pi}{\omega},

最低模態具有：

\delta^2S[\eta_1]<0,

但高頻模態仍為正。

因此作用量既有下降方向，也有上升方向：

\boxed{ q_{\mathrm{cl}} \text{ 是鞍點。} }

9. 數值驗證三：最小值轉為鞍點

9.1 $T=2$

因：

2<\pi,

離散 Hessian 最小特徵值為：

\boxed{ \lambda_{\min} \approx 0.012228>0. }

負特徵值數：

0.

沿最低模態擾動時：

S[q_{\mathrm{cl}}+\varepsilon\eta] - S[q_{\mathrm{cl}}] >0

對小非零 $\varepsilon$ 成立。

所以此時可以稱為「最小作用」。

9.2 $T=4$

因：

4>\pi,

離散 Hessian 最小特徵值為：

\boxed{ \lambda_{\min} \approx -0.006386<0. }

負特徵值數：

1.

沿負曲率模態：

S[q_{\mathrm{cl}}+\varepsilon\eta] - S[q_{\mathrm{cl}}] <0.

但其他方向仍可上升，所以不是最大，而是鞍點。

9.3 名稱上的結論

因此：

\boxed{ \text{最小作用量原理} }

只在二階形式正定的區域準確。

更一般的名稱必須是：

\boxed{ \text{駐定作用原理}. }

10. 變分—局部—拮抗三重表徵命題

命題 10.1

對作用量：

S[q] = \sum_{a=1}^n S_a[q]

以及滿足固定端點條件的正則駐定路徑 $q_{\mathrm{cl}}$ ，若：

\delta S[q_{\mathrm{cl}}]=0,

則：

\boxed{ \sum_{a=1}^n \delta S_a[q_{\mathrm{cl}}] = 0. }

同時，若作用量局部且滿足標準正則條件，則：

\boxed{ \delta S=0 \Longleftrightarrow \mathcal E_L[q]=0, }

其中 $\mathcal E_L$ 是 Euler–Lagrange 算子。

所以可建立：

\boxed{ \begin{aligned} \mathcal R_{\mathrm{path}} &: \delta S=0, \\ \mathcal R_{\mathrm{parts}} &: \sum_a\delta S_a=0, \\ \mathcal R_{\mathrm{local}} &: \mathcal E_L[q]=0. \end{aligned} }

它們是同一作用量動力的三個表徵。

命題地位

\Lambda = \mathrm{CON}.

其解析部分是標準變分法結果；數值部分已由本案例驗證。

11. 與《極限拮抗原理》的橋接

11.1 可以直接成立的弱連接

《極限拮抗原理》指出，某些穩定結構可以由相反生成項的持續抵消維持。

在變分層，簡諧振子滿足：

\delta S_K+\delta S_V=0.

所以可以說：

\boxed{ \text{駐定作用具有可被拮抗式解讀的變分結構。} }

這是一種表徵相容性。

11.2 尚未證明的中強命題

不能由：

\delta S_K+\delta S_V=0

直接推出：

$S_K$ 與 $S_V$ 是兩個接近物理上界的通道；
系統處於極限邊界；
存在非零調節餘量；
存在歷史連接泛函；
所有作用量皆來自 ETN。

因此：

\boxed{ \text{變分抵消} \neq \text{完整極限拮抗條件}. }

11.3 變分陰影猜想

可以保留以下較強猜想：

猜想 VAC

對某些可由作用量描述的持續系統，作用量駐定可能是底層正反生成、歷史連接與約束結構在路徑空間中的有效投影。

形式為：

\boxed{ \left( \mathcal U_+, \mathcal U_-, \mathcal C, \mathcal Q \right) \longrightarrow S_{\mathrm{eff}} \longrightarrow \delta S_{\mathrm{eff}}=0. }

或：

S_{\mathrm{eff}} = S_+ - S_- + S_C + S_{\mathcal Q}.

駐定條件：

\delta S_+ - \delta S_- + \delta S_C + \delta S_{\mathcal Q} = 0.

命題狀態：

\boxed{ \mathrm{CONJ/HYP}. }

本論文尚未給出從五輪極限拮抗動力到唯一作用量的完整逆變分證明。

12. 「一體多面」的正式判準

不能因兩個公式都含有減號，就稱為同一理論。

要把兩個描述視為同一動力的不同面，至少要求：

12.1 狀態對應

存在映射：

\Phi_{ij}: \mathcal R_i \to \mathcal R_j.

12.2 解集合對應

\Phi_{ij} \left( \operatorname{Sol}(\mathcal R_i) \right) = \operatorname{Sol}(\mathcal R_j).

12.3 可觀測一致

對允許觀測量 $\mathcal O$ ：

\boxed{ \mathcal O \left( \mathcal R_i \right) = \mathcal O \left( \mathcal R_j \right). }

12.4 邊界與約束一致

兩種表徵必須使用相容的：

初始條件；
邊界條件；
約束；
規範自由度；
時間區間。

12.5 失效域可說明

若轉換只在局部、保守或 Markov 條件下成立，必須標示失效域。

所以本文使用：

\boxed{ \text{表徵等價} }

而不輕率使用：

\boxed{ \text{本體完全同一}. }

13. 作用量分解的非唯一性

13.1 全導數自由

若：

L' = L+\frac{dF(q,t)}{dt},

則：

S' = S+F(q_T,T)-F(q_0,0).

固定端點時：

\delta S'=\delta S.

所以同一運動方程可以對應不同拉格朗日量。

13.2 正反分項不是唯一

即使：

L=T-V,

也可以重新分組：

L=(T+G)-(V+G).

因此「兩項抵消」本身不足以指定唯一底層實體。

13.3 對極限拮抗的要求

若要把：

S_+, S_-, S_C

賦予物理地位，必須額外給出：

對稱性；
可觀測映射；
尺度行為；
守恆律；
可識別參數；
不同分解間的不變量。

否則它只是方便的數學分解。

14. 非保守、耗散與歷史系統

標準：

S=\int L(q,\dot q,t)\,dt

最自然地描述封閉保守系統。

但現實中存在：

摩擦；
開放環境；
延遲；
記憶核；
熱流；
不可逆耗散。

這些系統可能需要：

Rayleigh 耗散函數；
倍增自由度；
非局部有效作用量；
Schwinger–Keldysh 或閉時路徑形式；
非保守作用項；
廣義變分原理。

若歷史連接為：

\mathcal C_t = \int_0^t K(t-s)q(s)\,ds,

則可考慮非局部作用量：

S_{\mathrm{nl}}[q] = \int dt\,L_0(q,\dot q) + \frac12 \int dt\,ds\, q(t)K(t,s)q(s).

其變分會產生時間非局部項。

因此，《極限拮抗原理》的歷史連接並不必然與駐定作用衝突；它可能要求更一般的作用量空間。

15. 可證偽與失效條件

F1：離散駐定與局部方程不一致

若在收斂網格下，離散駐定解不趨近 Euler–Lagrange 解，則離散方法或邊界設定有錯。

本案例未出現此問題。

F2：分項梯度不抵消

若：

\nabla S_K+\nabla S_V \neq0

於駐定解上顯著成立，則分解或數值解不正確。

F3：Hessian 不在共軛時間改變符號

若網格加密後最低 Hessian 特徵值不在：

T=\frac{\pi}{\omega}

附近穿越零，則第二變分推導或離散化有誤。

F4：駐定必然最小

本案例已構成反例： $T=4$ 時 $\delta S=0$ ，但存在負二階方向。

F5：變分抵消足以證明極限拮抗

若無法定義通道上界、歷史連接與餘量，則只能證明一般變分抵消，不能宣稱完整極限拮抗。

F6：不同表徵產生不同可觀測預測

若作用量形式與局部方程在相同條件下產生不同解，則不能稱為一體多面，而是不同理論。

16. 限制

本文只驗證一維簡諧振子；
只使用固定端點 Hamilton 變分；
未處理規範場、約束系統與奇異拉格朗日量；
未處理量子路徑積分中的測度與相位問題；
未處理一般非保守系統；
離散作用量採特定差分格式；
數值誤差未進行完整收斂階分析；
動能與位能分解並非唯一；
變分抵消不自動代表高強度物理拮抗；
尚未從《極限拮抗原理》五輪模型推導唯一作用量；
尚未證明 ETN 與 Hamilton 原理的必要關係；
「一體多面」目前是表徵方法論，不是本體同一性定理。

17. 結論

本文以一個最小但足以反駁常見誤解的模型，驗證了四件事。

第一：

\boxed{ \delta S=0 }

與：

\boxed{ m\ddot q+kq=0 }

是同一簡諧振子動力在路徑空間與局部時間中的等價表述。

第二，在：

S=S_K+S_V

的分解下，駐定路徑滿足：

\boxed{ \delta S_K=-\delta S_V. }

所以作用量駐定確可被觀察為不同變分貢獻的抵消。

第三，作用量是否為最小值取決於第二變分。於：

T<\frac{\pi}{\omega}

時，經典路徑是局部最小；於：

T>\frac{\pi}{\omega}

時，它仍然駐定，卻成為鞍點。

所以：

\boxed{ \text{最小作用量原理} \subset \text{駐定作用原理}. }

第四，駐定、局部動力與拮抗式抵消可以是同一體系的不同面，但這種等價只在給出明確轉換、相同解集合與可觀測一致時成立。

因此，對《極限拮抗原理》最穩健的橋接不是：

\text{極限拮抗取代駐定作用}.

而是：

\boxed{ \text{駐定作用可能是某些拮抗結構在路徑空間中的表徵。} }

目前已被驗證的是變分抵消與局部動力的等價。

尚待驗證的是：

\boxed{ \text{邊界逼近、歷史連接與極限餘量，} }

\boxed{ \text{能否共同生成一個具有物理差異性預測的有效作用量。} }

本文最終保留的「一體多面」原理是：

\boxed{ \text{同一動力結構在不同觀察空間中，} }

\boxed{ \text{可以被表述為駐定、局部運動、分項抵消或關係生成；} }

\boxed{ \text{表述不同不必意味世界不同，但等價必須被證明。} }

附錄 A：數值結果總表

案例	$T$	最大路徑誤差	最大作用量梯度	最大局部殘差	最小 Hessian 特徵值	負模態數
共軛時間前	$2$	$4.19\times10^{-6}$	$1.16\times10^{-13}$	$1.39\times10^{-11}$	$0.012228$	$0$
共軛時間後	$4$	$6.29\times10^{-5}$	$1.07\times10^{-13}$	$6.42\times10^{-12}$	$-0.006386$	$1$

附錄 B：核心命題帳本

命題	狀態
Hamilton 駐定條件導出 Euler–Lagrange 方程	`CON`
離散駐定路徑收斂至解析經典路徑	`CON`，本模型數值支持
動能與位能變分在駐定路徑抵消	`CON`
駐定不必是最小	`CON`
第一共軛時間後出現負二階方向	`CON`
駐定作用可被拮抗式表徵	`CON`，表徵層
所有駐定作用皆源於極限拮抗	`CONJ/HYP`
ETN 必然生成現有作用量	`UNK/HYP`

附錄 C：可重現性

資料包包含：

本文 Markdown；
完整 Python 驗證程式；
數值總表 CSV；
Hessian 時間掃描 CSV；
作用量曲率 CSV；
駐定解與 Euler–Lagrange 解比較圖；
動能與位能變分抵消圖；
最小特徵值跨越零點圖；
作用量最小／鞍點曲率比較圖；
SHA-256 檔案清單。

參考文獻

Lanczos, C. (1970). The Variational Principles of Mechanics. Dover Publications.
Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
Arnold, V. I. (1989). Mathematical Methods of Classical Mechanics (2nd ed.). Springer.
Gelfand, I. M., & Fomin, S. V. (2000). Calculus of Variations. Dover Publications.
Courant, R., & Hilbert, D. (1953). Methods of Mathematical Physics, Volume I. Wiley.
Feynman, R. P., & Hibbs, A. R. (1965). Quantum Mechanics and Path Integrals. McGraw-Hill.
Marsden, J. E., & West, M. (2001). Discrete mechanics and variational integrators. Acta Numerica, 10, 357–514. DOI: 10.1017/S096249290100006X.
Whiting, B. F. (2022). When is the principle of least action really a principle of least action? arXiv:2205.10318.
Galley, C. R. (2013). Classical mechanics of nonconservative systems. Physical Review Letters, 110, 174301. DOI: 10.1103/PhysRevLett.110.174301.
Neo.K（許筌崴）. (2026). 《極限拮抗原理：邊界逼近、動態抵消、歷史連接與持續差異》v2.0，未發表稿。

全文完

駐定作用的一體多面_從變分抵消局部動力到極限拮抗表徵_v1.0