| 第381行: |
第381行: |
| | | | |
| | ===算法=== | | ===算法=== |
| − | 一、ResDMD计算无谱污染的特征对算法:从切片数据中提取Koopman算子的准确特征对(特征值+特征函数),核心解决传统EDMD的谱污染问题。
| + | 一、ResDMD计算无谱污染的特征对算法<ref name=Colbrook_2024>Colbrook, M. J. & Townsend, A. Rigorous data‐driven computation of spectral properties of koopman operators for dynamical systems. Commun. Pure Appl. Math. 77, 221–283 (2024).</ref>:从切片数据中提取Koopman算子的准确特征对(特征值+特征函数),核心解决传统EDMD的谱污染问题。 |
| | | | |
| | 输入:<math>\{x^{(j)}\}_{j=1}^M</math>(初始状态)、<math>\{y^{(j)}=F(x^{(j)})\}_{j=1}^M</math>(下一时刻状态)、<math>\{w_j\}_{j=1}^M</math>(数值积分权重)、<math>\{\psi_j\}_{j=1}^{N_K}</math>(观测函数)、<math>\epsilon>0</math>(残差阈值) | | 输入:<math>\{x^{(j)}\}_{j=1}^M</math>(初始状态)、<math>\{y^{(j)}=F(x^{(j)})\}_{j=1}^M</math>(下一时刻状态)、<math>\{w_j\}_{j=1}^M</math>(数值积分权重)、<math>\{\psi_j\}_{j=1}^{N_K}</math>(观测函数)、<math>\epsilon>0</math>(残差阈值) |
| 第395行: |
第395行: |
| | | | |
| | | | |
| − | 二、ResDMD计算ε-伪谱估计算法:扩展上个算法的功能,计算Koopman算子的<math>\epsilon</math>-伪谱,弥补特征值无法反映的非正规算子暂态行为和谱分布信息。 | + | 二、ResDMD计算ε-伪谱估计算法<ref name=Colbrook_2024>Colbrook, M. J. & Townsend, A. Rigorous data‐driven computation of spectral properties of koopman operators for dynamical systems. Commun. Pure Appl. Math. 77, 221–283 (2024).</ref>:扩展上个算法的功能,计算Koopman算子的<math>\epsilon</math>-伪谱,弥补特征值无法反映的非正规算子暂态行为和谱分布信息。 |
| | | | |
| | 输入:<math>\{x^{(j)}\}_{j=1}^M</math>(初始状态)、<math>\{y^{(j)}=F(x^{(j)})\}_{j=1}^M</math>(下一时刻状态)、<math>\{w_j\}_{j=1}^M</math>(数值积分权重)、<math>\{\psi_j\}_{j=1}^{N_K}</math>(观测函数)、<math>\epsilon>0</math>(残差阈值)、<math>\{z_1,...,z_k\} \subset \mathbb{C}</math>(待评估的谱区域网格点) | | 输入:<math>\{x^{(j)}\}_{j=1}^M</math>(初始状态)、<math>\{y^{(j)}=F(x^{(j)})\}_{j=1}^M</math>(下一时刻状态)、<math>\{w_j\}_{j=1}^M</math>(数值积分权重)、<math>\{\psi_j\}_{j=1}^{N_K}</math>(观测函数)、<math>\epsilon>0</math>(残差阈值)、<math>\{z_1,...,z_k\} \subset \mathbb{C}</math>(待评估的谱区域网格点) |
| 第410行: |
第410行: |
| | 高斯迭代映射是一类具有复杂非线性动力学行为的实值动力系统,因分岔图呈现类似 “鼠标” 的特征而成为验证非线性算法的经典案例。该映射的非线性源于指数项<math>\exp(-\alpha x^2)</math>,当参数<math>\alpha</math>和<math>\beta</math>取特定值时,系统会表现出从周期运动到复杂分岔的过渡,其 Koopman 算子的谱结构包含离散特征值与潜在的连续谱成分,对算法的谱分辨能力提出挑战。 | | 高斯迭代映射是一类具有复杂非线性动力学行为的实值动力系统,因分岔图呈现类似 “鼠标” 的特征而成为验证非线性算法的经典案例。该映射的非线性源于指数项<math>\exp(-\alpha x^2)</math>,当参数<math>\alpha</math>和<math>\beta</math>取特定值时,系统会表现出从周期运动到复杂分岔的过渡,其 Koopman 算子的谱结构包含离散特征值与潜在的连续谱成分,对算法的谱分辨能力提出挑战。 |
| | | | |
| − | 实验采用的高斯迭代映射定义为:<math>F(x) = \exp(-\alpha x^2) + \beta</math>
| + | 实验<ref name=Colbrook_2024>Colbrook, M. J. & Townsend, A. Rigorous data‐driven computation of spectral properties of koopman operators for dynamical systems. Commun. Pure Appl. Math. 77, 221–283 (2024).</ref>采用的高斯迭代映射定义为:<math>F(x) = \exp(-\alpha x^2) + \beta</math> |
| | | | |
| | 式中,<math>\alpha=2</math>;<math>\beta=-1-\exp(-\alpha)</math>。映射的状态空间限制为<math>\Omega=[-1,0]</math>,采用标准勒贝格测度(Lebesgue measure)<math>\omega</math>描述状态分布。 | | 式中,<math>\alpha=2</math>;<math>\beta=-1-\exp(-\alpha)</math>。映射的状态空间限制为<math>\Omega=[-1,0]</math>,采用标准勒贝格测度(Lebesgue measure)<math>\omega</math>描述状态分布。 |