矩阵力学，外文名matrix mechanics，是量子力学其中一种的表述形式，提出者是海森堡、玻恩、约尔丹，时间1925年。
凡是矩阵力学，皆可建于以下的假定：
所有的物理量，均以厄米矩阵表之。一个物理系统的哈密顿函数 H {displaystyle mathbf {H} ,} 是广义坐标矩阵 Q {displaystyle mathbf {Q} ,} 及其共轭动量矩阵 P {displaystyle mathbf {P} ,} 的函数。
一个物理量 F {displaystyle mathbf {F} ,} 的观察值，是该矩阵的本征值 f n 1 n 2 {displaystyle f_{{n_{1}}{n_{2}}},} 。而能量 E n 1 n 2 {displaystyle E_{{n_{1}}{n_{2}}},} 是哈密顿函数 H {displaystyle mathbf {H} ,} 的本征值。
一个物理系统的广义坐标矩阵及其共轭动量矩阵满足以下的对易关系，亦称为 强量子条件 ： P Q − − --> Q P = ℏ ℏ --> i I {displaystyle mathbf {PQ} -mathbf {QP} ={hbar over i}mathbf {I} ,} I {displaystyle mathbf {I} ,} 为单位矩阵。
一个物理系统(如原子)的频率 ν ν --> n 1 n 2 {displaystyle u _{{n_{1}}{n_{2}}},} ，由频率条件定之：
h ν ν --> n 1 n 2 = E n 1 n 1 − − --> E n 2 n 2 {displaystyle hu _{{n_{1}}{n_{2}}}=E_{{n_{1}}{n_{1}}}-E_{{n_{2}}{n_{2}}},}
对易关系的思想来源
h ν ν --> n 1 n 2 = E n 1 n 1 − − --> E n 2 n 2 {displaystyle hu _{{n_{1}}{n_{2}}}=E_{{n_{1}}{n_{1}}}-E_{{n_{2}}{n_{2}}},} 这个条件是由玻尔的频率条件直接得来;但对易关系是如何引进的呢?如何得知新的力学形式是用矩阵去表达的呢? 其实海森堡的思想来源是先来自周期系统的解;周期系统的解全都可用傅里叶级数去展示：
q n ( t ) = ∑ ∑ --> n = 0 ∞ ∞ --> ( a n cos --> ( 2 π π --> n ν ν --> t ) + b n sin --> ( 2 π π --> n ν ν --> t ) ) = ∑ ∑ --> − − --> ∞ ∞ --> ∞ ∞ --> q n e 2 π π --> i n ν ν --> t {displaystyle q_{n}(t)=sum _{n=0}^{infty }(a_{n}cos(2pi nu t)+b_{n}sin(2pi nu t))=sum _{-infty }^{infty }q_{n}e^{2pi inu t}}
在此的 q n = 1 2 ( a n − − --> i b n ) {displaystyle q_{n}={ rac {1}{2}}(a_{n}-ib_{n}),} , q − − --> n = q n ∗ ∗ --> {displaystyle q_{-n}=q_{n}^{*},} 。 傅里叶级数有一个特点，就是对它进行运算，例如相加、相乘或微分，都不会产生 n ν ν --> , n = 1 , 2 , ⋯ ⋯ --> {displaystyle nu ,quad n=1,2,cdots ,} 以外的新频率系列。 但原子系统的频率是不能用傅里叶级数去表示，而是有一个叫里兹组合原则的经验关系：
ν ν --> n 1 n 2 = ν ν --> n 1 n 3 + ν ν --> n 3 n 2 {displaystyle u _{{n_{1}}{n_{2}}}=u _{{n_{1}}{n_{3}}}+u _{{n_{3}}{n_{2}}},}
如果频率能表示为经验项之差(如氢原子的里德伯公式)：
ν ν --> n 1 n 2 = T n 1 − − --> T n 2 {displaystyle u _{{n_{1}}{n_{2}}}=T_{n_{1}}-T_{n_{2}},}
里兹组合原则即可满足，而在这里原子系统形成一个“二维”的系统;对于频率的“二维”本性，海森堡用“二维”的广义坐标
q n 1 n 2 o e 2 π π --> i ν ν --> n 1 n 2 t {displaystyle q_{{n_{1}}{n_{2}}}^{o}e^{2pi iu _{{n_{1}}{n_{2}}}t},}
去取代傅里叶分量 q n e 2 π π --> i n ν ν --> t {displaystyle q_{n}e^{2pi inu t},} 。而为了模拟傅里叶级数，要求“二维”数集有以下关系：
q n 1 n 2 o = q n 1 n 2 o ∗ ∗ --> {displaystyle q_{{n_{1}}{n_{2}}}^{o}=q_{{n_{1}}{n_{2}}}^{o*},}
至于谱线 ν ν --> n 1 n 2 {displaystyle u _{{n_{1}}{n_{2}}},} 的幅度及偏振分别由 | q n 1 n 2 | 2 {displaystyle |q_{{n_{1}}{n_{2}}}|^{2},} 及 q n 1 n 2 {displaystyle q_{{n_{1}}{n_{2}}},} 复数的相位去表示。从里兹组合原则及对应原理，可以知道这类“二维”数集的乘法规则是：
( x x ) n 1 n 2 = ∑ ∑ --> j x n 1 j x j n 2 {displaystyle (xx)_{{n_{1}}{n_{2}}}=sum _{j}x_{{n_{1}}{j}}x_{{j}{n_{2}}},}
以使“二维”数集的运算，都不会产生 ν ν --> n 1 n 2 {displaystyle u _{{n_{1}}{n_{2}}},} 以外的新频率，如
( q q ) n 1 n 2 = ∑ ∑ --> k q n 1 k e 2 π π --> i ν ν --> n 1 k t q k n 2 e 2 π π --> i ν ν --> k n 2 t = ∑ ∑ --> k q n 1 k q k n 2 e 2 π π --> i ( ν ν --> n 1 k + ν ν --> k n 2 ) t = ( ∑ ∑ --> k q n 1 k q k n 2 ) e 2 π π --> i ν ν --> n 1 n 2 t {displaystyle (qq)_{{n_{1}}{n_{2}}}=sum _{k}q_{{n_{1}}{k}}e^{2pi iu _{{n_{1}}{k}}t}q_{{k}{n_{2}}}e^{2pi iu _{{k}{n_{2}}}t}=sum _{k}q_{{n_{1}}{k}}q_{{k}{n_{2}}}e^{2pi i(u _{{n_{1}}{k}}+u _{{k}{n_{2}}})t}=(sum _{k}q_{{n_{1}}{k}}q_{{k}{n_{2}}})e^{2pi iu _{{n_{1}}{n_{2}}}t},}
( q ˙ ˙ --> ) n 1 n 2 = 2 π π --> i ν ν --> n 1 n 2 q n 1 n 2 e 2 π π --> i ν ν --> n 1 n 2 t {displaystyle ({dot {q}})_{{n_{1}}{n_{2}}}=2pi iu _{{n_{1}}{n_{2}}}q_{{n_{1}}{n_{2}}}e^{2pi iu _{{n_{1}}{n_{2}}}t},}
海森堡只凭这些结果，就能得到谐振子的零点能是 1 2 h ν ν --> {displaystyle { rac {1}{2}}hu ,} ，但计算其间要多次运用对应原理，先引入玻尔-索末菲量子条件 J = ∮ p d q = n h {displaystyle J=oint p,dq=nh,} ，利用经典物理去估算量子物理的结果。
接着海森堡将他的结果转寄给玻恩，玻恩对于这些“二维”数集初时亦大感不解，后来他便意识到这些数集的运算与一个矩阵的运算是一模一样的，于是玻恩便与海森堡和约尔丹开展矩阵力学的建立。 首先，任何两个矩阵的乘法是不对易的：
A B − − --> B A ≠ ≠ --> 0 {displaystyle mathbf {AB} -mathbf {BA} eq mathbf {0} ,}
所以一个物理系统的广义坐标矩阵及其共轭动量满矩阵的乘积是不对易的：
P Q − − --> Q P ≠ ≠ --> 0 {displaystyle mathbf {PQ} -mathbf {QP} eq mathbf {0} ,}
那么这个乘积会等于什么呢?其实这个乘积等于什么可从玻尔-索末菲量子条件 J = ∮ p d q = n h {displaystyle J=oint p,dq=nh,} 加上对应原理预示出来。 对于任何周期系统，作用量有：
J = ∮ p d q = ∮ 0 1 ν ν --> p q ˙ ˙ --> d t {displaystyle J=oint p,dq=oint _{0}^{ rac {1}{u }}p{dot {q}},dt,}
如 p , q {displaystyle p,quad q,} 都使用傅里叶级数表示，就有：
J = ∫ ∫ --> 0 1 ν ν --> ∑ ∑ --> n 1 p n 1 e 2 π π --> i ν ν --> t ∑ ∑ --> n 2 2 π π --> i ν ν --> q n 2 e 2 π π --> i ν ν --> t d t = 2 π π --> i ν ν --> ∑ ∑ --> n , k ∫ ∫ --> 0 1 ν ν --> p n q n − − --> k ( k − − --> n ) e 2 π π --> i ν ν --> t d t = − − --> 2 π π --> i ∑ ∑ --> τ τ --> = − − --> ∞ ∞ --> ∞ ∞ --> τ τ --> p τ τ --> q − − --> τ τ --> {displaystyle J=int _{0}^{ rac {1}{u }}sum _{n_{1}}p_{n_{1}}e^{2pi iu t}sum _{n_{2}}2pi iu q_{n_{2}}e^{2pi iu t},dt=2pi iu sum _{n,k}int _{0}^{ rac {1}{u }}p_{n}q_{n-k}(k-n)e^{2pi iu t},dt=-2pi isum _{au =-infty }^{infty }au p_{au }q_{-au },}
所以 1 = ∂ ∂ --> J ∂ ∂ --> J = − − --> 2 π π --> i ∑ ∑ --> τ τ --> = − − --> ∞ ∞ --> ∞ ∞ --> τ τ --> ∂ ∂ --> ∂ ∂ --> J p τ τ --> q − − --> τ τ --> {displaystyle 1={ rac {partial J}{partial J}}=-2pi isum _{au =-infty }^{infty }au { rac {partial }{partial J}}p_{au }q_{-au },} 。
在玻尔-索末菲的理论中，作用量被量子化：
J = n h {displaystyle J=nh,}
况且 Δ Δ --> J = ( Δ Δ --> n ) h = τ τ --> h , τ τ --> ≡ ≡ --> Δ Δ --> n {displaystyle Delta J=(Delta n)h=au h,quad au quiv Delta n,} 。
由对应原理可知，经典理论的任何一个物理量 F {displaystyle F,} 的导数 ∂ ∂ --> F ∂ ∂ --> J {displaystyle { rac {partial F}{partial J}},} ，在量子理论中可用 Δ Δ --> F Δ Δ --> J = Δ Δ --> F τ τ --> h {displaystyle { rac {Delta F}{Delta J}}={ rac {Delta F}{au h}},} ，所以 ∂ ∂ --> ∂ ∂ --> J p τ τ --> q − − --> τ τ --> {displaystyle { rac {partial }{partial J}}p_{au }q_{-au },} 可用 1 τ τ --> h Δ Δ --> ( p τ τ --> q − − --> τ τ --> ) {displaystyle { rac {1}{au h}}Delta (p_{au }q_{-au }),} 替代，在新的理论中又可用 P , Q {displaystyle mathbf {P} ,mathbf {Q} ,} 表达式替代，即
1 τ τ --> h Δ Δ --> ( p τ τ --> q − − --> τ τ --> ) → → --> 1 τ τ --> h Δ Δ --> ( p n , n − − --> τ τ --> q n − − --> τ τ --> , n ) = − − --> 1 τ τ --> h ( p n , n − − --> τ τ --> q n − − --> τ τ --> , n − − --> q n , n + τ τ --> p n + τ τ --> , n ) {displaystyle { rac {1}{au h}}Delta (p_{au }q_{-au })ightarrow { rac {1}{au h}}Delta (p_{n,n-au }q_{n-au ,n})=-{ rac {1}{au h}}(p_{n,n-au }q_{n-au ,n}-q_{n,n+au }p_{n+au ,n}),}
将此代入上述的 1 = ∂ ∂ --> J ∂ ∂ --> J = − − --> 2 π π --> i ∑ ∑ --> τ τ --> = − − --> ∞ ∞ --> ∞ ∞ --> τ τ --> ∂ ∂ --> ∂ ∂ --> J p τ τ --> q − − --> τ τ --> {displaystyle 1={ rac {partial J}{partial J}}=-2pi isum _{au =-infty }^{infty }au { rac {partial }{partial J}}p_{au }q_{-au },} ，他们就得到关系式：
1 = 2 π π --> i h ∑ ∑ --> τ τ --> ( p n , n − − --> τ τ --> q n − − --> τ τ --> , n − − --> q n , n + τ τ --> p n + τ τ --> , n ) {displaystyle 1={ rac {2pi i}{h}}sum _{au }(p_{n,n-au }q_{n-au ,n}-q_{n,n+au }p_{n+au ,n}),}
这可用矩阵重新写成：
( p q − − --> q p ) n n = h 2 π π --> i = ℏ ℏ --> i {displaystyle (pq-qp)_{nn}={ rac {h}{2pi i}}={ rac {hbar }{i}},}
他们便作以下的假定：一个物理系统的广义坐标矩阵及其共轭动量满矩阵足以下的 对易关系 ：
I {displaystyle mathbf {I} ,} 为单位矩阵。
注意，千万不要以为对易关系能用玻尔-索末菲量子条件“推导”出来，更不要以为它可从经典物理推导出来，总之，对易关系是一个全新的假定，只有实验才能确认它的真实性。
海森堡运动方程及量子泊松括号
根据上文的对易关系，如果有一个矩阵函数(哈密顿函数) H = H ( Q , P ) {displaystyle mathbf {H} =mathbf {H} (mathbf {Q} ,mathbf {P} ),} ，我们有以下的关系：
