更改

跳到导航 跳到搜索
删除147字节 、 2025年11月17日 (一) 21:15
第16行: 第16行:     
==时间离散==
 
==时间离散==
−
离散时间场景适用于状态演化由离散映射描述的系统,如采样数据系统、迭代过程、数字控制中的离散时间模型等,是实际工程中数据驱动建模的常用场景。 设离散时间自治动力系统的状态演化满足<ref name="koopman_1931" /><ref name="koop_v_1932" /><ref name="brunton_2022" />:
+
离散时间场景适用于状态演化由离散映射描述的系统,如采样数据系统、迭代过程、数字控制中的离散时间模型等,是实际工程中数据驱动建模的常用场景。 设离散时间自治动力系统的状态演化满足:
    
<math> x_{k+1} = F(x_k), \quad x_0 \in X </math>
 
<math> x_{k+1} = F(x_k), \quad x_0 \in X </math>
   −
其中 <math> k = 0, 1, 2, \dots </math>是离散时间步,<math> F: X \to X </math>是离散流映射(可能非线性), <math> x_{k+1} </math>是第<math> k+1 </math> 步的状态。离散时间Koopman算子 <math> \mathcal{K}: \mathcal{G}(X) \to \mathcal{G}(X) </math>定义为:
+
其中 <math> k = 0, 1, 2, \dots </math>是离散时间步,<math> F: X \to X </math>是离散流映射(可能非线性), <math> x_{k+1} </math>是第<math> k+1 </math> 步的状态。离散时间Koopman算子 <math> \mathcal{K}: \mathcal{G}(X) \to \mathcal{G}(X) </math>定义为<ref name="koopman_1931" /><ref name="koop_v_1932" /><ref name="brunton_2022" />:
    
<math> \mathcal{K} g(x) = g\left( F(x) \right) </math>
 
<math> \mathcal{K} g(x) = g\left( F(x) \right) </math>
第31行: 第31行:  
连续时间场景适用于状态演化由[[常微分方程(ODE)]]或[[偏微分方程(PDE)]]描述的系统,如连续运动的机械系统、流体流动等。  
 
连续时间场景适用于状态演化由[[常微分方程(ODE)]]或[[偏微分方程(PDE)]]描述的系统,如连续运动的机械系统、流体流动等。  
   −
设连续时间自治动力系统的状态演化满足<ref name="koopman_1931" /><ref name="koop_v_1932" /><ref name="brunton_2022" />:
+
设连续时间自治动力系统的状态演化满足:
    
<math> \dot{x}(t) = f(x(t)), \quad x(0) = x_0 \in X </math>  
 
<math> \dot{x}(t) = f(x(t)), \quad x(0) = x_0 \in X </math>  
第37行: 第37行:  
其中<math> \dot{x} = \frac{dx}{dt}  </math>是状态的时间导数,<math> f: X \to \mathbb{R}^n  </math>是描述动力学的向量场(可能非线性)。   
 
其中<math> \dot{x} = \frac{dx}{dt}  </math>是状态的时间导数,<math> f: X \to \mathbb{R}^n  </math>是描述动力学的向量场(可能非线性)。   
   −
连续时间Koopman算子<math> \mathcal{K}^t: \mathcal{G}(X) \to \mathcal{G}(X) </math> 定义为:
+
连续时间Koopman算子<math> \mathcal{K}^t: \mathcal{G}(X) \to \mathcal{G}(X) </math> 定义为<ref name="koopman_1931" /><ref name="koop_v_1932" /><ref name="brunton_2022" />:
    
<math> \mathcal{K}^t g(x) = g\left( F^t(x) \right) </math>   
 
<math> \mathcal{K}^t g(x) = g\left( F^t(x) \right) </math>   
第63行: 第63行:  
<math> dx(t) = f(x(t)) dt + \sigma(x(t)) dW(t),\quad x(0) = x_0 \in X  </math>
 
<math> dx(t) = f(x(t)) dt + \sigma(x(t)) dW(t),\quad x(0) = x_0 \in X  </math>
   −
其中 <math> W(t) </math>是维纳过程(Wiener process,白噪声),<math> \varphi </math>是噪声强度矩阵,<math> \varphi </math> 和 <math> \varphi </math>分别是时间微分和噪声微分。
+
其中 <math> W(t) </math>是维纳过程(Wiener process,白噪声),<math> \sigma(x(t)) </math>是噪声强度矩阵。
    
随机 Koopman 算子 <math> \mathcal{K}^t_{\text{stoch}}:\mathcal{G}(X) \to \mathcal{G}(X) </math>定义为期望算子<ref name="brunton_2022" /><ref>N. C\v  rnjaric\'-Z\v ic, S. Mac\'es\v ic\', and I. Mezi\'c, Koopman operator spectrum for random dynamical systems, J. Nonlinear Sci., 30 (2020), pp. 2007--2056, <nowiki>https://doi.org/10.1007/s00332019-09582-z</nowiki></ref>:
 
随机 Koopman 算子 <math> \mathcal{K}^t_{\text{stoch}}:\mathcal{G}(X) \to \mathcal{G}(X) </math>定义为期望算子<ref name="brunton_2022" /><ref>N. C\v  rnjaric\'-Z\v ic, S. Mac\'es\v ic\', and I. Mezi\'c, Koopman operator spectrum for random dynamical systems, J. Nonlinear Sci., 30 (2020), pp. 2007--2056, <nowiki>https://doi.org/10.1007/s00332019-09582-z</nowiki></ref>:
第69行: 第69行:  
<math> \mathcal{K}^t_{\text{stoch}} g(x) = \mathbb{E}_x \left[ g\left( F^t(x, W) \right) \right] </math>
 
<math> \mathcal{K}^t_{\text{stoch}} g(x) = \mathbb{E}_x \left[ g\left( F^t(x, W) \right) \right] </math>
   −
其中,<math> \mathcal{K}^t_{\text{stoch}} </math>为随机 Koopman 算子,参数 t 为演化时间;<math> \varphi </math>为条件期望算子,下标 <math> \varphi </math> 表示初始状态为 <math> x </math>,期望针对噪声过程 <math> W(t) </math>计算;<math> F^t(x, W) </math>为随机流映射,表示初始状态<math> x </math>在噪声 <math> W(t) </math>作用下经时间 t 的演化结果(SDE 的解,依赖噪声实现);<math> g\left( F^t(x, W) \right) </math>为随机观测结果,随噪声 <math> W(t) </math>变化。
+
其中,<math> \mathcal{K}^t_{\text{stoch}} </math>为随机 Koopman 算子,参数 t 为演化时间;<math> \mathbb{E}_x \left[ \right] </math>为条件期望算子,下标表示初始状态为 <math> x </math>,期望针对噪声过程 <math> W(t) </math>计算;<math> F^t(x, W) </math>为随机流映射,表示初始状态<math> x </math>在噪声 <math> W(t) </math>作用下经时间 t 的演化结果(SDE 的解,依赖噪声实现);<math> g\left( F^t(x, W) \right) </math>为随机观测结果,随噪声 <math> W(t) </math>变化。
    
公式核心是通过期望消除噪声随机性,得到观测函数的平均演化 —— 随机 Koopman 算子作用后,输出为确定性函数(平均观测值)。
 
公式核心是通过期望消除噪声随机性,得到观测函数的平均演化 —— 随机 Koopman 算子作用后,输出为确定性函数(平均观测值)。
第77行: 第77行:  
<math> \mathcal{L}_{\text{stoch}} g(x) = \nabla g(x) \cdot f(x) + \frac{1}{2} \text{tr} \left( \sigma(x) \sigma(x)^T \nabla^2 g(x) \right) </math>
 
<math> \mathcal{L}_{\text{stoch}} g(x) = \nabla g(x) \cdot f(x) + \frac{1}{2} \text{tr} \left( \sigma(x) \sigma(x)^T \nabla^2 g(x) \right) </math>
   −
其中, <math> \nabla g \cdot f </math>为与确定性系统的生成元一致,对应无噪声时观测函数的瞬时变化率;<math> \varphi </math>为噪声贡献项,<math> \varphi </math>是矩阵迹(对角线元素之和),<math> \varphi </math>是 <math> g </math>的 Hessian 矩阵(二阶偏导数矩阵)。
+
其中, <math> \nabla g \cdot f </math>为与确定性系统的生成元一致,对应无噪声时观测函数的瞬时变化率;<math> σ(x) </math>为噪声贡献项,<math> Tr(⋅) </math>是矩阵迹(对角线元素之和),<math> \nabla^2 </math>是 <math> g </math>的 Hessian 矩阵(二阶偏导数矩阵)。
    
噪声项通过 Hessian 矩阵影响观测函数的演化,反映噪声对观测不确定性的贡献;当噪声强度<math> \sigma(x) = 0 </math>时,随机算子退化为确定性连续时间 Koopman 算子;随机 Koopman 算子的特征值可能包含衰减项(对应噪声导致的不确定性耗散),特征函数可刻画系统的 “平均线性化” 演化。
 
噪声项通过 Hessian 矩阵影响观测函数的演化,反映噪声对观测不确定性的贡献;当噪声强度<math> \sigma(x) = 0 </math>时,随机算子退化为确定性连续时间 Koopman 算子;随机 Koopman 算子的特征值可能包含衰减项(对应噪声导致的不确定性耗散),特征函数可刻画系统的 “平均线性化” 演化。
    
== 算子性质 ==
 
== 算子性质 ==
−
Koopman算子主要具有三种特性:线性性、半群性质、与Koopman特征函数的关系<ref name="brunton_2022" />。
+
Koopman算子主要具有两种特性:线性性、半群性质<ref name="brunton_2022" />。
    
1. 线性性:满足<math> \mathcal{K} (\alpha_1 g_1 + \alpha_2 g_2) = \alpha_1 \mathcal{K} g_1 + \alpha_2 \mathcal{K} g_2 </math> ;  
 
1. 线性性:满足<math> \mathcal{K} (\alpha_1 g_1 + \alpha_2 g_2) = \alpha_1 \mathcal{K} g_1 + \alpha_2 \mathcal{K} g_2 </math> ;  
   −
2. 半群性质:<math> \mathcal{K}^{m+n} = \mathcal{K}^m \circ \mathcal{K}^n </math>( m 步+ n 步演化等价于 m+n 步演化);
+
2. 半群性质:<math> \mathcal{K}^{m+n} = \mathcal{K}^m \circ \mathcal{K}^n </math>( m 步+ n 步演化等价于 m+n 步演化)。若<math> \varphi </math>是<math> \mathcal{K} </math> 的特征函数,满足<math> \mathcal{K} \varphi = \lambda \varphi </math>,则<math> \mathcal{K}^m \varphi = \lambda^m \varphi </math> (详细过程见下文Koopman模式分解)。  
−
 
  −
3. 特征函数关系:若<math> \varphi </math>是<math> \mathcal{K} </math> 的特征函数,满足<math> \mathcal{K} \varphi = \lambda \varphi </math>,则<math> \mathcal{K}^m \varphi = \lambda^m \varphi </math> (详细过程见下文Koopman模式分解)。
      
=Koopman模式分解=
 
=Koopman模式分解=
761

个编辑

导航菜单