第89行: |
第89行: |
| = 本征微观态及其演化 = | | = 本征微观态及其演化 = |
| 对于一个由<math>N</math>个主体组成的复杂系统,我们可以从实验测量或计算机模拟中获得主体的状态。依次使用时间<math>t=1, 2, … ,M</math>的状态,我们可以得到主体<math>i=1,2, … , N</math>的状态序列<math>S_i(t)</math>。 | | 对于一个由<math>N</math>个主体组成的复杂系统,我们可以从实验测量或计算机模拟中获得主体的状态。依次使用时间<math>t=1, 2, … ,M</math>的状态,我们可以得到主体<math>i=1,2, … , N</math>的状态序列<math>S_i(t)</math>。 |
| + | |
| | | |
| 一个主体<math>i</math>的平均状态是 | | 一个主体<math>i</math>的平均状态是 |
第94行: |
第95行: |
| <math>\begin{eqnarray}\delta {S}_{i}(t)={S}_{i}(t)-\langle {S}_{i}\rangle\end{eqnarray}</math>。这里定义一个具有所有主体的波动的微观态,它由一个<math>N</math>维矢量表示: | | <math>\begin{eqnarray}\delta {S}_{i}(t)={S}_{i}(t)-\langle {S}_{i}\rangle\end{eqnarray}</math>。这里定义一个具有所有主体的波动的微观态,它由一个<math>N</math>维矢量表示: |
| <math>\begin{eqnarray}\delta {\boldsymbol{S}}(t)=\left[\begin{array}{c}\delta {S}_{1}(t)\\ \delta {S}_{2}(t)\\ \vdots \\ \delta {S}_{N}(t)\end{array}\right]\end{eqnarray}</math>。 | | <math>\begin{eqnarray}\delta {\boldsymbol{S}}(t)=\left[\begin{array}{c}\delta {S}_{1}(t)\\ \delta {S}_{2}(t)\\ \vdots \\ \delta {S}_{N}(t)\end{array}\right]\end{eqnarray}</math>。 |
| + | |
| | | |
| 有了<math>M</math>个微观态,我们可以组成一个复杂系统的统计系综。这个系综由一个<math>N×M</math>的矩阵<math>\boldsymbol{A}</math>来描述,其元素为<math>\begin{eqnarray}{A}_{{it}}=\displaystyle \frac{\delta {S}_{i}(t)}{\sqrt{{C}_{0}}}\end{eqnarray}</math>,其中 | | 有了<math>M</math>个微观态,我们可以组成一个复杂系统的统计系综。这个系综由一个<math>N×M</math>的矩阵<math>\boldsymbol{A}</math>来描述,其元素为<math>\begin{eqnarray}{A}_{{it}}=\displaystyle \frac{\delta {S}_{i}(t)}{\sqrt{{C}_{0}}}\end{eqnarray}</math>,其中 |
| <math> {C}_{0}={\sum }_{t=1}^{M}{\sum }_{i=1}^{N}\delta {S}_{i}^{2}(t)</math>,<math>\boldsymbol{A}</math>的列序与微观态的演变相一致。 | | <math> {C}_{0}={\sum }_{t=1}^{M}{\sum }_{i=1}^{N}\delta {S}_{i}^{2}(t)</math>,<math>\boldsymbol{A}</math>的列序与微观态的演变相一致。 |
| | | |
− | <math>t</math>和<math>t^{\prime}</math>的微观态之间的相关性由它们的矢量乘积定义:<math>\begin{eqnarray}{C}_{{tt}^{\prime} }=\delta {\boldsymbol{S}}{\left(t\right)}^{{\rm{T}}}\cdot \delta {\boldsymbol{S}}(t^{\prime} )=\displaystyle \sum _{i=1}^{N}\delta {S}_{i}(t)\delta {S}_{i}(t^{\prime} )\end{eqnarray}</math>。以<math>{C}_{tt}^{\prime} </math>作为其元素,我们可以得到一个<math>M×M</math>的微观态相关矩阵:<math>\begin{eqnarray}{\boldsymbol{C}}={C}_{0}{{\boldsymbol{A}}}^{{\rm{T}}}\cdot {\boldsymbol{A}}\end{eqnarray}</math>,其轨迹<math>Tr\boldsymbol{C}={\sum }_{t=1}^{M}{C}_{tt}={C}_{0}</math>。
| |
| | | |
− | 统计系综也可以被视为主体的动态微观态的系综,它们由<math>M</math>维向量描述:<math>\begin{eqnarray}\delta {{\boldsymbol{S}}}_{i}=\left[\delta {S}_{i}(1),\delta {S}_{i}(2),\cdots ,\delta {S}_{i}(M)\right]\end{eqnarray}</math>,其中<math>i=1,2, … , N</math>。
| + | <math>t</math>和<math>t^{\prime}</math>的微观态之间的相关性由它们的矢量乘积定义:<math>\begin{eqnarray}{C}_{{tt}^{\prime} }=\delta {\boldsymbol{S}}{\left(t\right)}^{{\rm{T}}}\cdot \delta {\boldsymbol{S}}(t^{\prime} )=\displaystyle \sum _{i=1}^{N}\delta {S}_{i}(t)\delta {S}_{i}(t^{\prime} )\end{eqnarray}</math>。以<math>{C}_{tt}^{\prime} </math>作为其元素,我们可以得到一个<math>M×M</math>的微观态相关矩阵:<math>\begin{eqnarray}{\boldsymbol{C}}={C}_{0}{{\boldsymbol{A}}}^{{\rm{T}}}\cdot {\boldsymbol{A}}\end{eqnarray}</math>,其轨迹<math>Tr\boldsymbol{C}={\sum }_{t=1}^{M}{C}_{tt}={C}_{0}</math>。相关矩阵<math>\boldsymbol{C}</math>有<math>M</math>个特征向量<math>V_J</math>,其中<math>J = 1, 2, ⋯ , M</math>,我们可以用它们组成一个<math>M×M</math>的单元矩阵 |
| + | <math>\boldsymbol{V}=\left[\boldsymbol{V}_1 \boldsymbol{V}_2 \ldots \boldsymbol{V}_M\right]</math> |
| + | |
| + | |
| + | 此外,我们在这里考虑动态微观状态<math>δS_i</math>和<math>δS_j</math>之间的相关性:<math>\begin{eqnarray}{{\boldsymbol{S}}}_{j}^{{{T}}}=\displaystyle \sum _{t=1}^{M}\delta {S}_{i}(t)\delta {S}_{j}(t)\end{eqnarray}</math>。以<math>K_{ij}</math>为元素,我们可以得到一个<math>N×N</math>的动态微观态的相关矩阵: |
| + | <math>\begin{eqnarray}{\boldsymbol{K}}={C}_{0}{\boldsymbol{A}}\cdot {\boldsymbol{A}}^{T}\end{eqnarray}</math>,其轨迹<math>Tr\boldsymbol{K}={\sum }_{t=1}^{N}{K}_{ii}={C}_{0}</math>。相关矩阵<math>\boldsymbol{K}</math>有<math>N</math>个特征向量<math>U_I</math>,其中<math>I = 1, 2, ⋯ , N</math>,我们可以用它们组成一个<math>N×N</math>的单元矩阵 |
| + | <math>\boldsymbol{U}=\left[\boldsymbol{U}_1 \boldsymbol{U}_2 \ldots \boldsymbol{U}_N\right]</math>。 |
| + | |
| + | |
| + | 根据奇异值分解,系综矩阵<math>\boldsymbol{A}</math>可以被分解为<math>\begin{eqnarray}{\boldsymbol{A}}={\boldsymbol{U}}\cdot {\boldsymbol{\Sigma }}\cdot {{\boldsymbol{V}}}^{{\rm{T}}},\end{eqnarray}</math> |
| + | |
| + | 对<math>\boldsymbol{A}</math>进行奇异值分解,统计系综可被分解为: |
| + | |
| + | |
| + | <math>A=\sum_{i=1}^r \sigma_i U_i \otimes V_i</math> |
| + | |
| + | |
| + | 这里的<math>r=min (N, M)</math>,其中,<math>U_i</math>为本征微观态,<math>V_i</math>为该本征微观态遵循的时间演化,⊗表示克罗内克积,其中<math>\sigma_i</math>表示<math>U_i</math>在系综中的概率,且满足归一化条件 |
| | | |
− | 此外,我们在这里考虑动态微观状态<math>δS_i</math>和<math>δS_j</math>之间的相关性:<math>\begin{eqnarray}{{\boldsymbol{S}}}_{j}^{{rm{T}}}=\displaystyle \sum _{t=1}^{M}\delta {S}_{i}(t)\delta {S}_{j}(t)\end{eqnarray}</math>。以<math>K_{ij}</math>为元素,我们可以得到一个<math>N×N<\math>的动态微观态的相关矩阵:
| |
− | <math>\begin{eqnarray}{\boldsymbol{K}}={C}_{0}{\boldsymbol{A}}\cdot {{\boldsymbol{A}}}^{{{rm{T}}}\end{eqnarray}</math>,其轨迹<math>Tr\boldsymbol{K}={\sum }_{t=1}^{N}{K}_{ii}={C}_{0}</math>。
| |
| | | |
| + | <math>\sum_{i=1}^r \sigma_i^2=1</math> |
| | | |
| | | |
| + | 用这样的方法,可以将原来相互关联的微观态转变为相互独立的本征微观态,就可以将最初的微观态用本征微观态进行线性组合,这个线性组合的大小与本征值<math>\sigma_i</math>(即权重因子)成正比,即<math>\sigma_i</math>越大,占比越多。 |
| | | |
| + | <math></math> |
| | | |
| | | |