Skip to content

17.2.1 吸引子上的概率测度 ​

17.2.1.1 不变测度 ​

1. 定义、支撑在吸引子上的测度 ​

设 {φt}t∈Γ 是 (M,ρ) 上动力系统, B 为 M 上所有博雷尔 (Borel) 集合构成的 σ 代数 (参见第 905 页 12.9.1,2.),令 μ:B→[0,+∞] 为 B 上的测度,且对任意 t∈Γ,φt 关于 μ 是可测的. 如果对任意 A∈B,t>0 ,有 μ(φ−t(A))=μ(A) , 则称 μ 为关于 {φt}t∈Γ 的不变测度. 如果系统 {φt}t∈Γ 可逆,则上述不变测度定义可叙述为 μ(φt(A))=μ(A),(A∈B,t>0) . 对博雷尔集合 A⊂M 上,如果 μ(M∖A)=0 ,则称 μ 支撑在 A 上. 设 Λ 为系统 {φt}t∈Γ 的一个吸引子, μ 为关于 {φt}t∈Γ 的不变测度,如果对任意博雷尔集合 B 满足 Λ∩B=∅ ,有 μ(B)=0 ,则称 μ 支撑在 Λ 上.

对测度 μ:B→[0,+∞] 定义 μ 的支撑为所有 μ 支撑于其上的闭集合的交. - A: 对 M=[0,1] 上的伯努利转移映射:

(17.28a)xℓ+1=2xℓ(mod1).

系统 φ:[0,1]→[0,1] 定义为

(17.28b)φ(x)={2x,0≤x≤12,2x−1,12<x≤1.

从定义可看出勒贝格测度是其不变测度. 若将 x∈[0,1) 写作二进制形式 x= ∑n=1∞⁡an⋅2−n (其中 an=0 或1),将其等同于 x=.a1a2a3 ,则其在算子 2x(mod1) 下的像为. a1′a2′a3′⋯ ,其中 ai′=ai+1 ,即所有数字 ak 都向左移动一位,同时首位数字去掉.

◼B: 定义映射 Ψ:[0,1]→[0,1] 为

(17.29)Ψ(y)={2y,0≤y<12,2(1−y),12≤y≤1.

该映射称为帐篷映射,勒贝格测度为该映射的一个不变测度. 通过同胚 h:[0,1)→ [0,1),h(x)=2πarcsin⁡x 可将系统 (17.5) 转变为系统 (17.29). 因此,当 α= 4 时, 系统 (17.5) 有一个绝对连续不变测度. 事实上, 对系统 (17.29) 的密度函数 ρ1(y)≡1 ,可找到系统 (17.5) 在 α=4 时对应密度函数 ρ(x) ,二者满足关系 ρ1(y)=ρ(h−1(y))|(h−1)(y)| ,即 ρ(x)=1πx(1−x).

◼C: 设 x0 为可逆离散系统 {φi} 的一个稳定 T 周期点,则 μ=1T∑i=0T−1⁡δφi(x0) 为 {φi} 的概率测度,其中 δx0 为支撑在 x0 点的狄拉克(Dirac)测度. (参见第 905 页 12.9.1,2.).

2. 自然测度 ​

设 Λ⊂M 为系统 {φt}t∈Γ 的一个吸引子,吸引域为 W . 对任意博雷尔集合 A⊂W 及任意 x0∈W ,定义

(17.30)μ(A;x0):=limT→∞⁡t(T,A,x0)T,

其中 t(T,A,x0) 表示轨道段 {φt(x0)}t=0T 落入集合 A 的所有时刻构成的集合. 若对 λ -几乎处处 x0∈W ,有 μ(A;x0)=α ,则定义 μ(A)=μ(A;x0) . 由于对几乎所有初始点 x0∈W 的轨道随着 t 趋向于 +∞ 都趋向于 Λ ,因此 μ 是支撑在集合 Λ 上的概率测度.

17.2.1.2 遍历论基础 ​

1. 遍历系统 ​

称 (M,ρ) 上具有不变测度 μ 的动力系统 {φt}t∈Γ 是遍历的 (或测度 μ 是遍历的),如果对所有满足 φ−t(A)=A(∀t>0) 的博雷尔集合 A ,有 μ(A)=0 或 μ(M∖A)=0 . 对离散系统 {φt}(17.3) ,其中 φ:M→M 是同胚, M 是紧致度量空间, 总存在不变的遍历测度.

◼A: 对圆周 S1 上的旋转映射

(17.31)xt+1=xt+Φ(mod2π),t=0,1,⋯,

φ:[0,2π)→[0,2π),φ(x)=x+Φ(mod2π) . 勒贝格测度是关于 φ 的不变测度. 如果 Φ2π 为无理数,则系统 (17.3) 是遍历的; 如果 Φ2π 为无理数,则系统 (17.3) 不是遍历的. _____ B :稳定平衡点或者稳定周期轨作为吸引子的动力系统关于自然测度是遍历的.

伯克霍夫 (Birkhoff) 遍历定理 设动力系统 {φt}t∈Γ 关于不变的概率测度 μ 是遍历的,则对任意可积函数 h∈L1(M,B,μ) ,对 μ -几乎处处 x0∈M,x0 点沿正半轨 {φtx0}t=0∞ 的时间平均 h¯(x0) 等于空间平均 ∫Mhdμ ,其中对连续系统, h¯(x0)= limT→+∞⁡1T∫0Th(φt(x0))dt ,对离散系统, h¯(x0)=limn→+∞⁡1n∑i=0n−1⁡h(φi(x0))dt .

2. 物理测度或 SRB 测度 ​

遍历定理叙述只有在测度 μ 的支撑集充分大时才有意义. 设 φ:M→M 是一个连续映射, μ:B→R 是一个不变测度. 如果对任何连续函数 h:M→R ,由所有

满足条件

(17.32a)limn→+∞⁡1n∑i=0n−1⁡h(φi(x0))dt=∫Mhdμ

的点 x0 构成的集合具有正勒贝格测度,则称 μ 是一个 SRB 测度(该命名源于赛奈 (Sinai),鲍恩 (Bowen) 和吕埃勒 (Ruelle),见 [17.6]). 如果对几乎所有 x∈M ,测

度序列

(17.32b)μn:=1n∑i=0n−1⁡δφi(x)

弱收敛于 μ ,即 ∫Mhdμn→∫Mhdμ,n→+∞ ,其中 δx 为支撑在 x 点的狄拉克测度,那么 μ 是一个 SRB 测度.

◼ 对一些重要吸引子, 如埃农吸引子, SRB 测度的存在性已被证明.

3. 混合系统 ​

称 (M,ρ) 上具有不变测度 μ 的动力系统 {φt}t∈Γ 是混合的,如果对任意博雷尔集合 A,B⊂M ,有 limt→+∞⁡μ(A∩φ−t(B))=μ(A)μ(B) . 对混合系统而言,由所有满足 t=0 时刻在集合 A ,经过充分长时间 φt 作用后落入集合 B 的点所构成的集合的测度只与乘积项 μ(A)μ(B) 有关.

混合系统也是遍历的: 设 {φt} 是一个混合系统,若博雷尔集合 A 满足 φ−t(A)= A(t>0) ,则有 μ(A)2=limt→+∞⁡μ(A∩φ−t(A))=μ(A) ,从而 μ(A)=0 或 1

系统 (17.1) 的流 {φt} 称为混合的,如果对任意 g,h∈L2(M,B,μ) 有

(17.33)limt→+∞⁡∫M[g(φt(x))−g¯][h(x)−h¯]dμ=0

成立,其中 g¯ 和 h¯ 表示空间平均,可由时间平均替换.

◼ 映射 (17.28a) 是混合的. 旋转映射 (17.31) 关于概率测度 λ2π 不是混合的.

4. 自相关函数 ​

设 M 上关于不变测度 μ 的动力系统 {φt}t∈Γ 是遍历的. 设 h:M→R 为某连续函数, {φt(x)}t≥0 为某半轨,将空间平均 h¯ 分别在连续和离散情形下被替换成时间平均,即连续情形替换为 limT→∞⁡1T∫0Th(φt(x))dt ,离散情形替换为 limn→∞⁡1n∑i=0n−1⁡h(φi(x)) , 则函数 h 沿半轨 {φt(x)}t≥0 到时间 τ≥0 的相关函数在流的情形定义为

(17.34a)Ch(τ)=limT→∞⁡1T∫0Th(φt+τ(x))h(φt(x))dt−h¯2,

在离散情形定义为

(17.34b)Ch(τ)=limn→∞⁡1n∑i=0n−1⁡h(φi+τ(x))h(φi(x))−h¯2.

自相关函数也可对负向时间迭代定义,这时 Ch(⋅) 看作 R 或 Z 上的偶函数.

周期或拟周期轨相应地导致 Ch 的周期或拟周期行为. 对任意测试函数 h,Ch(τ) 随 τ 的增加而快速减小的现象称为混沌. 如果 Ch(τ) 随 τ 的增加以指数速率衰减, 则说明系统混合.

5. 功率谱 ​

称 Ch(τ) 的傅里叶变换为功率谱 (参见第 1028 页 15.3.1.2,5.),记为 Ph(ω) . 在连续时间情形下,在假设 ∫−∞+∞|Ch(τ)|dτ<∞ 下,

(17.35a)Ph(ω)=∫−∞+∞Ch(τ)e−iωτdτ=20inftyCh(τ)cos⁡(ωτ)dτ

在离散时间情形下,若 ∑k=−∞+∞⁡|Ch(k)|<+∞ 成立,则

(17.35b)Ph(ω)=Ch(0)+2∑k=1+∞⁡Ch(k)cos⁡ωk.

如果 Ch(⋅) 的绝对可积或可加性假设不成立,则在大多数重要情形下, Ph 看作一个分布. 对应于动力系统周期行为的能量谱可被刻画为等距脉冲. 对拟周期行为, 能量谱也存在脉冲, 它们是拟周期行为基本脉冲的整系数线性组合. 在宽带谱中出现奇异峰值可被认为混沌行为标志.

◼A: 设 φ 是系统 (17.1) 的 T 周期轨道,试验函数 h 满足时间平均 h(φ(t)) 等于 0 . 若 h(φ(t)) 的傅里叶表示为

h(φ(t))=∑k=−∞+∞⁡αkeikω0t,ω0=2πT,

则有

Ch(τ)=∑k=−∞+∞⁡|αk|2cos⁡(kω0τ),Ph(ω)=2π∑k=−∞+∞⁡|αk|2δ(ω−kω0).

◼B: 设 φ 是系统 (17.1) 的一个拟周期轨道,试验函数 h 满足沿系统 φ 的时间平均等于 0 . 设 h(φ(t)) 可表示为 (双傅里叶级数)

h(φ(t))=∑k1=−∞+∞⁡∑k2=−∞+∞⁡αk1k2ei(k1ω1+k2ω2)t,

则

Ch(τ)=∑k1=−∞+∞⁡∑k2=−∞+∞⁡|αk1k2|2cos⁡(k1ω1+k2ω2)τ,Ph(ω)=2π∑k1=−∞+∞⁡∑k2=−∞+∞⁡|αk1k2|2δ(ω−k1ω1−k2ω2).

version 1.24.0