主題
Search

哈密頓對映

考慮一個一維哈密頓對映 的形式

H(p,q)=1/2p^2+V(q),
(1)

它滿足哈密頓方程

q^.=(partialH)/(partialp)
(2)
p^.=-(partialH)/(partialq).
(3)

現在，寫出

q^._i=((q_(i+1)-q_i))/(Deltat),
(4)

其中

q_i=q(t)
(5)
q_(i+1)=q(t+Deltat).
(6)

那麼運動方程變為

q_(i+1)=q_i+p_iDeltat
(7)
p_(i+1)=p_i-Deltat((partialV)/(partialq_i))_(q=q_i).
(8)

注意方程 (7) 和 (8) 不是保面積的，因為

(partial(q_(i+1),p_(i+1)))/(partial(q_i,p_i))=|1 -Deltat(partial^2V)/(partialq_i^2); Deltat 1|
(9)
=1+(Deltat)^2(partial^2V)/(partialq_i^2)
(10)
!=1.
(11)

然而，如果我們用以下公式代替 (9) 和 (10)，

q_(i+1)=q_i+p_iDeltat
(12)
p_(i+1)=p_i-Deltat((partialV)/(partialq_i))_(q=q_(i+1))
(13)
(partial(q_(i+1),p_(i+1)))/(partial(q_i,p_i))=|1 -Deltatpartial/(partialq_i)((partialV)/(partialq))_(q=q_(i+1)); Deltat 1|
(14)
=1+(Deltat)^2(partial^2V)/(partialq_i^2)=1,
(15)

它是保面積的。

另請參閱

保面積對映

使用 探索

請引用為

Weisstein, Eric W. “哈密頓對映。” 來自 Web 資源。 https://mathworld.tw/HamiltonianMap.html

主題分類