眾所周知,分析力學是建立在「最小作用量原理」之上的。它有很多名字,有很多形式和理解。這裡就用我喜歡、而且自認為比較容易的一種理解: 作用量:S=int_{t_!}^{t_2}Ldt , 作用量變分為0:delta S=0

它表示有一個可以描述系統狀態的量,它在時間上面的積分就是系統狀態的演化量,這樣的演化量需要滿足最小化原理。

因為這個量被要求描述系統的狀態,所以它的一組參數就被要求能完整描述系統狀態,我們把這樣能完整且不重複地描述系統狀態的一組參數叫做廣義坐標。

當然,拉格朗日選擇的是這樣一套坐標: L(mathop{q}limits ^. , q,t) ,其中 mathop{q}limits ^. 代表速度, q 代表位置, t 代表時間。大家可以感受一下,只用這三個量確實可以完全確定一個系統的狀態。(覺得記起來困難的朋友可以自動進行符號替換: mathop{q}limits ^.=v , q=r )。

首先介紹一下等時變分,它是指 滿足delta t=0 的變分,變分記號 delta 可以不嚴謹地理解成跟 微分記號d 一樣的東西。下面由等時變分導出拉格朗日方程( Lagrange equation ):

delta S=int_{t_1}^{t_2}(frac{partial L}{partial mathop{q}limits ^.}delta mathop{q}limits ^. + frac{partial L}{partial q}delta q)dt=0

在式 delta mathop{q}limits ^.=delta frac{dq}{dt} 中,對 q 交換記號(能夠交換記號的充分條件是,該函數被兩個記號分別作用後都連續),有 delta mathop{q}limits ^.=frac{d(delta q)}{dt}

對於積分中的第一項可以分部積分變為:

frac{partial L}{partial mathop{q}limits ^.}delta mathop{q}limits ^.=frac{partial L}{partial mathop{q}limits ^.} frac{d(delta q)}{dt}=frac{partial L}{partialmathop{q}limits ^. }delta q-delta qfrac{d}{dt}(frac{partial L}{partial mathop{q}limits ^.}) ,因為廣義坐標之間相互獨立,顯然 frac{partial L}{partialmathop{q}limits ^. }delta q0

所以 frac{partial L}{partial mathop{q}limits ^.}delta mathop{q}limits ^.=-delta qfrac{d}{dt}(frac{partial L}{partial mathop{q}limits ^.})

所以原積分變為:

delta S=int_{t_1}^{t_2}(-delta qfrac{d}{dt}(frac{partial L}{partial mathop{q}limits ^.})+ frac{partial L}{partial q}delta q)dt=0

delta qdt 提出: delta S=int_{t_1}^{t_2}( frac{partial L}{partial q}-frac{d}{dt}(frac{partial L}{partial mathop{q}limits ^.}))delta qdt

由此: frac{partial L}{partial q}-frac{d}{dt}(frac{partial L}{partial mathop{q}limits ^.})=0

即: frac{d}{dt}(frac{partial L}{partial mathop{q}limits ^.})=frac{partial L}{partial q}拉格朗日方程)。

對比牛頓第二定律 frac{dp}{dt}=ma=F

很容易看出, frac{partial L}{partial mathop{q}limits ^.}=pfrac{partial L}{partial q}=F 。所以我們把 frac{partial L}{partial mathop{q}limits ^.} 叫做廣義動量,把 frac{partial L}{partial q} 叫做廣義力

按照經典物理中動量式 p=mv ,由 frac{partial L}{partial mathop{q}limits ^.}=p 積分可以得到 :L=int_{}^{}mvdv=frac{1}{2}mv^2+C(q)

其中 C(q) 表示與 q (或者說 r )有關,而與 mathop{q}limits ^. (或者說 v )無關的函數。

再根據 frac{partial L}{partial q}=Ffrac{d}{dt}(frac{partial L}{partial mathop{q}limits ^.})-frac{partial L}{partial q}=0 可以確定 C(q) 的符號取負,並且 C(q) 就是勢能 U 。也就是說外力滿足: F=-frac{partial U}{partial r} (比如 F=qE=-frac{partial (qU)}{partial r} )。

所以我們有: L=T-V ,其中 T 代表動能, V 代表勢能。系統的拉格朗日量( Lagrangian )就是總動能減總勢能。

現在定義系統總能量 E=T+V 。顯然, E=2T-(T-V)=mv^2-L

現在,我們重新選取一組廣義坐標p,q,t ),其中, p 是動量, q 是位置, t 是時間。於是有 E=frac{p^2}{2m}+U(q) ,我們將其定義為哈密頓量 H(Hamiltonian)

這在物理上叫正則變換。之所以叫這個名字,是因為變換後定義的哈密頓量 H ,關於廣義坐標 p,q 高度對稱。如下是哈密頓正則方程

-frac{partial H}{partial q}=mathop{p}limits ^. (力方程) frac{partial H}{partial p}=mathop{q}limits ^. (速度方程)

他們與 F=-frac{partial U}{partial r}=frac{dp}{dt},frac{partial H}{partial p}=frac{d}{dp}(frac{p^2}{2m}+U(r))=v 等價。

這裡再簡單提一下泊松括弧Poisson Bracket ):

有一力學量 f ,可以用廣義坐標( p,q,t )描述。

frac{df}{dt}=frac{partial f}{partial t}+frac{partial f}{partial q}frac{dq}{dt}+frac{partial f}{partial p}frac{dp}{dt}

frac{df}{dt}=frac{partial f}{partial t}+frac{partial f}{partial q}frac{partial H}{partial p}-frac{partial f}{partial p}frac{partial H}{partial q}

定義泊松括弧[f,H]=frac{partial f}{partial q}frac{partial H}{partial p}-frac{partial f}{partial p}frac{partial H}{partial q}

frac{df}{dt}=frac{partial f}{partial t}+[f,H] 。當力學量不顯含時間時(通俗地說就是跟時間無關),有 [f,H]=0

值得一提的是, egin{aligned} &[q,H]=mathop{q}limits ^.\ &[p,H]=mathop{p}limits ^.\ end{aligned}


推薦閱讀:
相關文章