761
个编辑
更改
→理论
:<math>[\mathcal{K}g](\boldsymbol{x}) = \Psi(F(\boldsymbol{x})) \boldsymbol{g} = \Psi(\boldsymbol{x}) (K \boldsymbol{g}) + \underbrace{\sum_{j=1}^{N_K} \psi_j(F(\boldsymbol{x})) g_j - \Psi(\boldsymbol{x}) (K \boldsymbol{g})}_{r(g, \boldsymbol{x})}.</math>
:<math>[\mathcal{K}g](\boldsymbol{x}) = \Psi(F(\boldsymbol{x})) \boldsymbol{g} = \Psi(\boldsymbol{x}) (K \boldsymbol{g}) + \underbrace{\sum_{j=1}^{N_K} \psi_j(F(\boldsymbol{x})) g_j - \Psi(\boldsymbol{x}) (K \boldsymbol{g})}_{r(g, \boldsymbol{x})}.</math>
残差项衡量的是'''在子空间<math>V_{N_K}</math>上近似Koopman算子'''的误差大小。如果残差项为0，说明我们完美地用一个线性系统重构出了原本的非线性系统。通常情况下，<math>V_{N_K}</math>不是<math> \mathcal{K} </math>的不变子空间，因此不存在能使<math>\displaystyle{ r(g, \boldsymbol{x}) }</math>对所有<math>\displaystyle{ g \in V_{N_K} }</math>和<math>\boldsymbol{x} \in \mathcal{X}</math>均为零的<math>\displaystyle{ K }</math>。
'''优化目标：'''因此，我们的核心任务就变成了找到一个最优的矩阵<math>K</math>，让这个残差项的整体误差最小。为了在系统整个定义域上量化整体误差，我们要让残差项在定义域上求积分。对应的优化问题就是：
'''优化目标：'''因此，我们的核心任务就变成了找到一个最优的矩阵<math>K</math>，让这个残差项的整体误差最小。为了在系统整个定义域上量化整体误差，我们要让残差项在定义域上求积分。对应的优化问题就是：
:<math>\displaystyle{ \arg\min_{B \in \mathbb{C}^{N_K \times N_K}} \left\{ \int_{\mathcal{X}} \max_{\substack{g \in \mathbb{C}^{N_K}, \\ \|\boldsymbol{g}\|=1}} |r(g, \boldsymbol{x})|^2 d\omega(\boldsymbol{x}) = \int_{\mathcal{X}} \|\Psi(F(\boldsymbol{x})) - \Psi(\boldsymbol{x}) B\|_{\ell^2}^2 d\omega(\boldsymbol{x}) \right\} }</math>
:<math>\displaystyle{ \arg\min_{B \in \mathbb{C}^{N_K \times N_K}} \left\{ \int_{\mathcal{X}} \max_{\substack{g \in \mathbb{C}^{N_K}, \\ \|\boldsymbol{g}\|=1}} |r(g, \boldsymbol{x})|^2 d\omega(\boldsymbol{x}) = \int_{\mathcal{X}} \|\Psi(F(\boldsymbol{x})) - \Psi(\boldsymbol{x}) B\|_{\ell^2}^2 d\omega(\boldsymbol{x}) \right\} }</math>
其中<math>\|\mathbf{ · }\|_2</math>表示[[向量]]的标准欧几里得范数。不过实际计算中无法直接对上式的求积分（状态空间连续，无法遍历所有状态），因此通过 “数值积分规则”（以<math>\displaystyle{ \{x_j\}_{j=1}^M }</math>为[[节点]]，<math>\boldsymbol{y}^{(j)} = F(\boldsymbol{x}^{(j)}) </math>，<math>\displaystyle{ \{w_j\}_{j=1}^M }</math>为权重，<math>M </math>为[[初始条件]]数量）近似积分。通过这种方式，把连续的积分转化为离散的加权和，原来的优化问题就变成了更容易计算的“[[加权最小二乘问题]]”：
其中<math>\|\mathbf{ · }\|_2</math>表示[[向量]]的标准欧几里得范数。不过实际计算中无法直接对上式的求积分（状态空间连续，无法遍历所有状态），因此通过 '''数值积分规则'''（以<math>\displaystyle{ \{x_j\}_{j=1}^M }</math>为[[节点]]，<math>\boldsymbol{y}^{(j)} = F(\boldsymbol{x}^{(j)}) </math>，<math>\displaystyle{ \{w_j\}_{j=1}^M }</math>为权重，<math>M </math>为[[初始条件]]数量）近似积分。通过这种方式，把连续的积分转化为离散的加权和，原来的优化问题就变成了更容易计算的'''[[加权最小二乘问题]]'''：
:<math>\arg\min_{B\in\mathbb{C}^{N_K \times N_K}} \sum_{j=1}^{M} w_j \|\Psi(y^{(j)}) - \Psi(x^{(j)})B\|_{\ell_2}^2</math>
:<math>\arg\min_{B\in\mathbb{C}^{N_K \times N_K}} \sum_{j=1}^{M} w_j \|\Psi(y^{(j)}) - \Psi(x^{(j)})B\|_{\ell_2}^2</math>