在图像处理、环境科学、生物医学、通信与电子工程等方面的研究中, 系统的混沌性是人们非常关注的研究方向之一. 上世纪以来,学者们研究了许多复杂的非线性动力学行为, 比如,2011年,谭枫,张瑞丰[1]在动力系统(X,f)中讨论了敏感性、 F -敏感性和F -敏感对的关系, 证明了系统存在敏感对不一定意味着系统是Li-Yorke混沌的, 并给出了传递系统中敏感对存在的充分条件. 2012年,尹建东和 周作领[2]证明了传递的非极小系统是Tankens-Ruelle混沌的. 又如,2006年, Oprocha[3]证明了一个量子谐振子的湮没算符是分布ε -混沌的 (其中ε=116), 并且其准测度不小于116. 2011年,吴新星和 朱培勇[4]证明了这个算符能够达到最大意义上的分布混沌,并得到准测度的确切值为1. 2013年,黎日松等[5]则考查了耦合系数为1的 一类耦合映象格子的分布(p,q) -混沌性和准测度. 更多工程应用系统中的混沌研究见文献[6, 7, 8, 9]以及其参考文献.
然而,大多数文献都是在自治系统(X,f)中研究混沌性, 而现实中很多系统不具备自治系统的优良性质,不同的扰动需要用不同的函数来描述,即 同一个系统存在一系列不同的映射对其进行作用, 这表明工程实际中许多系统都属于非自治系统. 非自治系统最初由 Kolyada和 Snoha在文献[10]中提出, 随后他们又在[11, 12]等文献中讨论了非自治系统的极小性和区间上非自治分段单调动力系统的拓扑熵. 非自治系统的概念也与一些非自治不等式相联系,此处不赘述, 有兴趣的读者可以参见文献[13, 14]以及它们的参考文献等.
令I=[0,1],I上的度量记为ρ. 考虑一个映射序列fn:I→I,n∈N,将这个映射序列记为f1,∞=(f1,f2,⋯). 这个序列定义了一个非自治离散系统(I,f1,∞), 在这个映射序列之下,点x∈I的轨道为Orb(x,f1,∞)=(fn1(x)),这里fn1=fn∘⋯∘f1,f01表示恒等映射. 如果对∀n∈N, 有fn=f,则(I,f1,∞)即为自治离散系统(I,f).
这篇文章主要讨论了fn,∞的分布混沌性与乘积系统f[m]n,∞ (m为正整数)的分布混沌性的关系(第3部分). 在下文中, 我们总是假设fn(n∈N)为满射, 这个条件是处理这类系统时的常见条件,比如文献[12, 15, 16]等. 应用本文的主要结果, 容易得到文献[3, 4, 5]中的系统的混沌性在迭代之下是保持的.
混沌的严格数学定义最初于1975年由李天岩和 Yorke在文献[17]中给出, 即Li-Yorke混沌.后来,为了从不同角度描述系统的复杂性和不可预测性, 各种不同的混沌定义相继提出,比如Deveney混沌[18],本质混沌[19], 稠混沌[20],稠δ -混沌[21],Li-Yorke敏感[22]等等. Li-Yorke混沌定义的一个重要推广是分布混沌. 分布混沌的定义由 Schweizer和 Smital在文献[23]中提出, 类似自治系统中这个分布混沌的定义, 我们在非自治情形下给出一个分布混沌的定义.
定义2.1 对∀n∈N,映射序列fn,∞=(fn,fn+1,⋯)称为分布混沌的,如果存在不可数子集S⊂I使得∀x,y∈S,x≠y,有
(i)~ ∀t>0,F∗xy(t,fn,∞)=lim supk→∞1k−n+1k∑i=nχ[0,t)(ρ(fi∘⋯∘fn(x),fi∘⋯∘fn(y)))=1,
(ii)~ ∃t>0,Fxy(t,fn,∞)=lim infk→∞1k−n+1k∑i=nχ[0,t)(ρ(fi∘⋯∘fn(x),fi∘⋯∘fn(y)))=0, \\ 其中t∈R,χ[0,t)是集合[0,t)的特征函数,即 当s∈[0,t)时,χ[0,t)(s)=1; 当s∉[0,t)时,χ[0,t)(s)=0. F∗xy(t,fn,∞)和Fxy(t,fn,∞)分别称为fn,∞的上、下分布函数, 不可数子集S称为fn,∞在I中的分布攀援集.
现在,我们回顾一下非自治情形下Li-Yorke混沌的定义(类似文献[24]所给定义).
定义2.2 对∀n∈N,映射序列fn,∞称为Li-Yorke混沌的, 如果存在不可数子集S⊂I使得∀x,y∈S,x≠y,有 lim supk→∞ρ(fk∘⋯∘fn(x),fk∘⋯∘fn(y))>0 且 lim infk→∞ρ(fk∘⋯∘fn(x),fk∘⋯∘fn(y))=0. 不可数子集S称为fn,∞在I中的Li-Yorke攀援集.
定理2.3 如果映射序列fn,∞是分布混沌的, 则fn,∞是Li-Yorke混沌的.
证 假设对任意不可数子集S⊂I,存在两个不同的点x,y∈S使得 lim infk→∞ρ(fk∘⋯∘fn(x),fk∘⋯∘fn(y))=α>0. 则 ∃N∈N,∀k≥N,ρ(fk∘⋯∘fn(x),fk∘⋯∘fn(y))≥α2. 因此, lim supk→∞1k−n+1k∑i=nχ[0,α2)(ρ(fi∘⋯∘fn(x),fi∘⋯∘fn(y)))=0. 这与分布混沌定义的条件(i)矛盾.
类似的,如果任意的不可数子集S⊂I,存在两个不同的点x,y∈S使得 lim supk→∞ρ(fk∘⋯∘fn(x),fk∘⋯∘fn(y))=0. 这将与分布混沌的条件(ii)矛盾.
故,fn,∞是Li-Yorke混沌的.
定理2.4 映射序列f1,∞是分布混沌的当且仅当对∀n∈N,n≥2,映射序列fn,∞是分布混沌的.
证 只需考虑n=2的情形.
(充分性)~~ 设S⊂I是f2,∞的一个分布攀援集,则,∀x,y∈S,x≠y,有
(i)~ ∀t>0,F∗xy(t,f2,∞)=lim supk→∞1k−1|{2≤i≤k:ρ(fi∘⋯∘f2(x),fi∘⋯∘f2(y))<t}|=1,
(ii)~ ∃t0>0,Fxy(t0,f2,∞)=lim infk→∞1k−1|{ρ(fi∘⋯∘f2(x),fi∘⋯∘f2(y))<t0}|=0,
其中,|A|表示集合A的基数.
因为f1是满射,则对∀x,y∈S,x≠y,存在x∗,y∗∈I,x∗≠y∗使得f1(x∗)=x,f1(y∗)=y.
取S中每个元素在f1之下的一个逆象构成集合T, 则T是一个不可数集,并且
(a)~ 对∀t>0,若d(f1(x∗),f1(y∗))<t时,令c=1; 若d(f1(x∗),f1(y∗))≥t时,令c=0,则有 F∗x∗y∗(t,f1,∞)=lim supk→∞1k|{1≤i≤k:ρ(fi∘⋯∘f1(x∗),fi∘⋯∘f1(y∗))<t}|=lim supk→∞1k(|{2≤i≤k:ρ(fi∘⋯∘f1(x∗),fi∘⋯∘f1(y∗))<t}|+c)=lim supk→∞1k|{2≤i≤k:ρ(fi∘⋯∘f2(x),fi∘⋯∘f2(y))<t}|=lim supk→∞1k−1|{2≤i≤k:ρ(fi∘⋯∘f2(x),fi∘⋯∘f2(y))<t}|=1.
(b)~ 对(ii)中的t0, Fx∗y∗(t0,f1,∞)=lim infk→∞1k|{1≤i≤k:ρ(fi∘⋯∘f1(x∗),fi∘⋯∘f1(y∗))<t0}|=lim infk→∞1k(|{2≤i≤k:ρ(fi∘⋯∘f2(x),fi∘⋯∘f2(y))<t0}|+c)=lim infk→∞1k−1|{2≤i≤k:ρ(fi∘⋯∘f2(x),fi∘⋯∘f2(y))<t0}|=0 从而,f1,∞是分布混沌的.
(必要性)~~ 类似充分性的讨论,设S⊂I是f1,∞的一个分布攀援集, 则f1(S)是f2,∞的一个分布攀援集.
由定理2.3和定理2.4,下面的推论成立.
推论2.5 映射序列f1,∞是Li-Yorke混沌的当且仅当对∀n∈N,n≥2,fn,∞是Li-Yorke混沌的.
在这一节里,我们讨论非自治乘积系统的分布混沌性.
∀m,n∈N,令 g1=fn+m−1∘⋯∘fn,g2=fn+2m−1∘⋯∘fn+m,⋯, gp=fn+pm−1∘⋯∘fn+(p−1)m,⋯. (I,g1,∞)称为(I,fn,∞)的一个乘积系统. 为了更方便地看出系统(I,g1,∞)与系统(I,fn,∞)的关系, 乘积系统g1,∞也记作f[m]n,∞. 并且记 fkn=fn+k−1∘⋯∘fn (∀n≥1).
定理3.1 若映射序列f1,∞是分布混沌的,则对∀m∈N,f[m]1,∞也是分布混沌的.
证明设S0是f1,∞的一个分布攀援集,对∀m∈N,下证S0也是f[m]1,∞的一个分布攀援集.
(i)~ 首先,∀x,y∈S0:x≠y,∀t>0,我们证明 F∗xy(t,f[m]1,∞)=lim supk→∞1kk∑i=1χ[0,t)(ρ(f[im]1(x),f[im]1(y)))=1.
假设存在x0,y0∈S0:x0≠y0以及t0>0使得 F∗x0y0(t0,f[m]1,∞)<1,则 lim supk→∞1kmkm∑i=1χ[0,t0)(ρ(fi1(x),fi1(y)))=lim supk→∞1kmm∑j=1k−1∑i=0χ[0,t0)(ρ(fmi+j1(x),fmi+j1(y)))≤1mm∑j=1(lim supk→∞1kk−1∑i=0χ[0,t0)(ρ(fmi+j1(x),fmi+j1(y))))≤1m(m−1+F∗x0y0(t0,f[m]1,∞))<1.
由于S0是f1,∞的分布攀援集, 则F∗x0y0(t02,f1,∞)=1, 由上极限的定义,存在严格递增的序列{ns}s∈N使得 F∗x0y0(t02,f1,∞)=lims→∞1nsns∑i=1χ[0,t02)(ρ(fi1(x0),fi1(y0)))=1. 令 M0={km:k∈Z+}∩{ns:s∈N}, M1={km+1:k∈Z+}∩{ns:s∈N}, ⋯⋯ Mm−1={km+m−1:k∈Z+}∩{ns:s∈N}, 则存在l∈{0,1,⋯,m−1}使得Ml为无限集. 记Ml={nls:s∈N},我们有 lims→∞1nlsnls∑i=1χ[0,t02)(ρ(fi1(x0),fi1(y0)))=1. 从而 1≥lim sups→∞1nls+m−lnls+m−l∑i=1χ[0,t0)(ρ(fi1(x0),fi1(y0)))≥lim sups→∞1nls+m−lnls∑i=1χ[0,t0)(ρ(fi1(x0),fi1(y0)))≥lim sups→∞nlsnls+m−l1nlsnls∑i=1χ[0,t02)(ρ(fi1(x0),fi1(y0)))=1. 由Ml的定义,nls是除以m余l的数,因此m|nls+m−l,我们有 lim supk→∞1kmkm∑i=1χ[0,t0)(ρ(fi1(x0),fi1(y0)))≥lim sups→∞1nls+m−lnls+m−l∑i=1χ[0,t0)(ρ(fi1(x0),fi1(y0)))=1. 这与lim supk→∞1kmkm∑i=1χ[0,t0)(ρ(fi1(x),fi1(y)))<1矛盾.
(ii)~ 下证存在δ>0使得对∀x,y∈S0:x≠y都有 F∗xy(δ,f[m]1,∞)=lim infk→∞1kk∑i=1χ[0,δ)(ρ(fim1(x),fim1(y)))=0.
事实上,因为S0是f1,∞的分布攀援集, 则存在δ1>0 使得Fxy(δ1,f1,∞)=0对任意x,y∈S0:x≠y都成立.
假设存在xδ,yδ∈S0, 使得对∀δ>0都有Fxδyδ(δ,f[m]1,∞)>0. 特别地, 对δ∗=δ12>0,存在 xδ∗,yδ∗∈S0使得 Fxδ∗yδ∗(δ∗,f[m]1,∞)>0.
因为Fxδ∗yδ∗(δ1,f1,∞)=0, 则存在严格递增的序列{ms}s∈N使得 lims→∞1msms∑i=1χ[0,δ1)(ρ(fi1(xδ∗),fi1(yδ∗)))=Fxδ∗yδ∗(δ1,f1,∞)=0.
类似(i),可以找到c∈{0,1,⋯,m−1}使得mcs={km+c:k∈Z+}∩{ms:s∈N}为无限集. 从而 0=lims→∞1mcsmcs∑i=1χ[0,δ1)(ρ(fi1(xδ∗),fi1(yδ∗)))=lims→∞1mcs+m−cmcs+m−c∑i=1χ[0,δ1)(ρ(fi1(xδ∗),fi1(yδ∗)))=lims→∞1mcs+m−cm∑j=11m(mcs+m−c)−1∑i=0χ[0,δ1)(ρ(fim+j1(xδ∗),fim+j1(yδ∗)))≥m∑j=1lim infs→∞1mcs+m−c1m(mcs+m−c)−1∑i=0χ[0,δ1)(ρ(fim+j1(xδ∗),fim+j1(yδ∗)))≥lim infs→∞1mmmcs+m−c1m(mcs+m−c)−1∑i=0χ[0,δ1)(ρ(f(i+1)m1(xδ∗),f(i+1)m1(yδ∗)))≥1mFxδ∗yδ∗(δ1,f[m]1,∞)≥1mFxδ∗yδ∗(δ∗,f[m]1,∞)>0. 矛盾.
综上(i),(ii)知f[m]1,∞是分布混沌的.
这个定理表明,f1,∞的分布混沌性在乘积运算下是保持的. 那么,如果∀m∈N (m≥2),f[m]1,∞是分布混沌的,是否有f1,∞是分布混沌的呢? 为了回答这个问题,我们先给出三个引理.
引理3.2 若映射序列(fn)∞n=1一致收敛于f,则∀m∈N,m≥2,序列(fmn)∞n=1一致收敛于fm.
证 当m=2时. 因为序列(fn)∞n=1一致收敛于映射f, 则对∀ε>0,存在N0∈N使得对∀x∈I和∀n≥N0,都有ρ(fn(x),f(x))<ε.
由三角不等式, ρ(f2n(x),f2(x))≤ρ(fn+1(fn(x)),f(fn(x)))+ρ(f(fn(x)),f(f(x))), 又因为f一致连续,容易证明(f2n)∞n=1 一致收敛于f2.
假设(fkn)∞n=1一致收敛于fk. 由于对∀x∈X,都有 ρ(fk+1n(x),fk+1(x))≤ρ(fn+k(fkn(x)),f(fkn(x)))+ρ(f(fkn(x)),f(fk(x))) 可以得到(fk+1n)∞n=1一致收敛于fk+1.
由归纳法,对∀m∈N,序列 (fmn)∞n=1一致收敛于fm.
注意到闭区间上映射序列一致收敛等价于等度连续(等度连续的定义见文献[12]), 下面的引理成立.
引理3.3 若映射序列(fn)∞n=1等度连续,则∀m∈N,m≥2,序列(fmn)∞n=1等度连续.
从而,可以得到
引理3.4 若映射序列(fn)∞n=1等度连续(或一致收敛于f), 则对∀ε>0和∀m∈N, 存在δ(ε)和N(m)∈N,使得,对∀x,y∈I:ρ(x,y)<δ(ε)和n≥N(m),都有 ρ(fmn(x),fmn(y))<ε.
证由引理3.3和下面的不等式,结论显然成立. ρ(fmn(x),fmn(y))≤ρ(fmn(x),fm(x))+ρ(fm(x),fm(y))+ρ(fm(y),fmn(y)). 证毕.
现在我们来回答上面提出的问题.
定理3.5若(fn)∞n=1等度连续(或一致收敛于f), 如果∀m∈N,m≥2,f[m]1,∞是分布混沌的,那么f1,∞也是分布混沌的.
证 设D0是f[m]1,∞的一个分布攀援集,对x,y∈D0:x≠y,我们证明
(i)~ 若对∀δ>0有F∗x,y(δ,f[m]1,∞)=1,则对∀t>0有F∗x,y(t,f1,∞)=1.
事实上,∀t>0,由引理3.4,存在δ:0<δ<t以及N∈N,使得对∀x,y∈I:ρ(x,y)<δ和∀k≥N,有ρ(fik(x),fik(y))<t对任意i:1≤i≤m都成立. 从而,对∀j≥N:ρ(fjm1(x),fjm1(y))<δ和∀i:1≤i≤m,有 ρ(fjm+i1(x),fjm+i1(y))<t. 则 k∑j=1χ[0,δ)(ρ(fjm1(x),fjm1(y)))−N≤1m(k+1)m∑j=1χ[0,t)(ρ(fj1(x),fj1(y))). 因此, F∗x,y(t,f1,∞)=lim supk→∞1kk∑j=1χ[0,t)(ρ(fj1(x),fj1(y)))=lim supk→∞1(k+1)m(k+1)m∑j=1χ[0,t)(ρ(fj1(x),fj1(y)))=lim supk→∞1km(k+1)m∑j=1χ[0,t)(ρ(fj1(x),fj1(y)))=lim supk→∞1k(k∑j=1χ[0,δ)(ρ(fjm1(x),fjm1(y)))−N)=1.
下面我们证明
(ii)~ 若存在s>0使得Fx,y(s,f[m]1,∞)=0, 则存在δ>0使得Fx,y(δ,f1,∞)=0.
事实上,对上述s>0,存在δ:0<δ<s以及 N1∈N,使得对∀x,y∈I:ρ(x,y)<δ和∀k≥N1, 不等式ρ(fik(x),fik(y))<s对任意i:1≤i≤m都成立. 从而,对∀j≥N1:ρ(fjm1(x),fjm1(y))≥s和∀i:1≤i≤m,我们有ρ(fjm−i1(x),fjm−i1(y))≥δ.
因此, k∑j=1χ[s,∞)(ρ(fjm1(x),fjm1(y)))−N1≤1mkm∑j=1χ[δ,∞)(ρ(fj1(x),fj1(y))). 又因为 lim supk→∞1k(k∑j=1χ[s,∞)(ρ(fjm1(x),fjm1(y))))=1−lim infk→∞1k(k∑j=1χ[0,s)(ρ(fjm1(x),fjm1(y))))=1−Fx,y(s,f[m]1,∞)=1, 则 lim supk→∞1k(k∑j=1χ[δ,∞)(ρ(fj1(x),fj1(y))))=1. 因此, Fx,y(δ,f1,∞)=1−lim supk→∞1k(k∑j=1χ[δ,∞)(ρ(fj1(x),fj1(y))))=0. 故f1,∞是分布混沌的.
由定理2.3,定理3.1,定理3.5,下面两个推论成立.
推论3.6 若映射序列f1,∞是Li-Yorke混沌的,则∀m∈N,f[m]1,∞也是Li-Yorke混沌的.
推论3.7若(fn)∞n=1等度连续(或一致收敛于f),对∀m∈N (m≥2),如果f[m]1,∞是Li-Yorke混沌的, 那么f1,∞也是Li-Yorke混沌的.
注1考虑(fn)∞n=1一致收敛于f的情形. 2010年,Canovas[25]在紧度量空间上证明了f1,∞ 的Li-Yorke混沌性不一定意味着f的Li-Yorke混沌性. 2012年,Dvorakova[16]证明了f1,∞的分布攀援集的全 Lebesgue测度不能保证极限映射f的分布混沌性. 由本文第2部分的定理2.3和文献[25]的结论,我们可以知道, f1,∞的分布混沌性不一定意味着f的分布混沌性,这与文献[16]的结论一致. 那么,是否f1,∞的混沌性与极限映射f的混沌性没有任何关系呢?2006年, Canovas在文献[26]中告诉我们, f1,∞的ω -极限集也是f的ω -极限集,此外,早在1996年, Kolyada[10]就在紧度量空间上证明了若(fn)∞n=1一致收敛于f, 则有hhop(f1,∞)≤hhop(f) (其中hhop(⋅)表示拓扑熵). 因此,一个自然的问题是, 当(fn)∞n=1一致收敛于f时,除上述混沌行为之外, f1,∞和f的其他混沌性质有什么关系?
注2 在自治系统中, 关于分布混沌还定义了分布(p,q) -混沌、一致分布混沌、第一型分布混沌(DC1)、 第二型分布混沌(DC2)、第三型分布混沌(DC3)等等, 这些定义都可以被推广到非自治系统中, 那么它们又有些什么性质呢?这些都是值得我们进一步研究的问题.