Hamilton力学的辛算法简介

  • Hamilton力学的辛算法简介
  • 冯康

  • 我的熟人和我
  • 我的朋友和我
  • 只能说 嗨 蛮好的

外微分形式与辛几何

外微分形式

  • 1-形式
  • 2-形式
  • 闭2-形式(辛构造)

Euclid Space

  • 符合如下内积定义的线性空间V称为Euclid空间
    • 对称性 (a, b) = (b, a)
    • 线性
      • (a, kb) = k(a, b)
      • (a+c, b) = (a, b) + (c, b)
    • 非简并性
      • (a, a)>= 0 当且仅当 a=0 时 (a, a) = 0

Simplectic Space

defination 

  • 具有如下内积定义的线性空间W称为 辛空间
    • 反对称性 (a, b) = -(b, a)
    • 双线性
      • (a_1+a_2, b) = (a_1, b) + (a_2, b)
      • (a, b_1+b_2) = (a, b_1) + (a, b_2)
    • 非简并性
      • 若向量a 对于W中的任意向量b均匀(a, b)=0,则a=0

Simplectic Space 

  • 度量
    • 做功
    • 面积 体积
    • 流量
  • 辛内积
    • 单位辛矩阵 J_{2n}\equiv \begin{bmatrix} 0 &1_n \\ -1_n & 0 \end{bmatrix}
    • <a, b>\equiv (a,J_{2n}b)\equiv \sum_{i=1}^n (a_ib_{n+i}-a_{n+i}b_i)

单位辛矩阵的性质

\begin{matrix} J^2=-1\\ J^T=J^{-1}=-J\\ \end{matrix}

decision theorem of simplectic matrix 

Euclid space and Simplectic space 

  • 内积
  • 单位阵与单位辛矩阵
  • 正交与辛正交
  • 正交归一基于共轭正交归一基
  • 对称变换与Hamilton变换
    • 实对称矩阵的本征值为实数
    • 实对称矩阵不同本征值的本征向量正交
    • 实对称矩阵所有的本征向量构成一组正交归一基

Simplectic structure of Simplectic space 

q\equiv [q_1,q_2,...,q_f]^T\,\,\,p \equiv [p_1,p_2,...,p_f]^T

\dot{q_i}=\frac{\partial H}{\partial p_i}\,\,\, \dot{p_i}=\frac{\partial H}{\partial q_i}

z\equiv \begin{bmatrix} p\\ q \end{bmatrix}\,\,\, \dot{z}\equiv \begin{bmatrix} \dot{p}\\ \dot{q} \end{bmatrix}

J^{-1}\frac{\partial H}{\partial z}=\dot{z}

Simplectic structure of canonical transformation

basic concept

  1. J^{-1}\frac{\partial H}{\partial z}=\dot{z}
  2. y=Mz\,\,\, \dot{y}=M\dot{z}
  3. J^{-1}\frac{\partial H}{\partial z}=\dot{z}  J^{-1}\frac{\partial H}{\partial y}=\dot{y}
  4. simplectic transformation MJM^T=J

Properties of Simplectic transformation 

\left\{\begin{matrix} MJM^T=J\\ (det \,M)^2=1 \end{matrix}\right.

infinitesimal Simplectic matrix

define

  • B^TJ_{2n}+J_{2n}B=0,则该阵为无穷小辛阵

Dcision Theorem

  1. 若C为对称阵,当且仅当B = J_{2n}C 时,B为无穷小辛阵

Advantage of Simplectic algorithm 

  • 保证相空间体积不变
  • 保证力学守恒量守恒

你可能感兴趣的:(#,Hamilton,系统,数学建模,抽象代数,python,numpy,mathematica)