混沌理论及应用

1、混沌理论及应用第1页，共84页，2022年，5月20日，2点13分，星期五混沌的概念：混沌（chaos）又称浑沌，人们通常用它来描述混乱、杂乱无章、乱七八糟的状态，在这个意义上它与无序的概念是相同的。 一、混沌的基本概念及特征第2页，共84页，2022年，5月20日，2点13分，星期五1确定性 在混沌系统中，描述系统演化的动力学方程的确定性，是指方程(常微分方程、差分方程、时滞微分方程)是非随机的，不含任何随机项。系统的未来(或过去)状态只与初始条件及确定的演化规则有关，即系统的演化完全是由内因决定的，与外在因素无关。这是至关重要的一条限制，所以我们现在讲的混沌也叫“确定性混沌”。正因为确定性

2、的系统出现了复杂行为，也叫内随机性，人们才兴奋起来，才一往倾心地钻研混沌。当然，从长远的观点来看，人们肯定会研究带有随机项的更复杂系统的非周期运动。然而，目前由于公众对混沌还有相当的误解，所以我们严格区分是否为确定性至关重要，还不能笼统地从现象的层次把一大堆似是而非的东西都称为混沌。总之，混沌概念的狭义化总比泛化好些。现在我们考虑的混沌主要是一种时间演化行为，不直接涉及空间分布变化，所以暂不考虑偏微分方程。第3页，共84页，2022年，5月20日，2点13分，星期五例：Lorenz系统Logistic 映射第4页，共84页，2022年，5月20日，2点13分，星期五2非线性 产生混沌的系统一定

3、含有非线性因素，有了非线性未必产生混沌，但没有非线性是肯定产生不了混沌的。也就是说，非线性是产生混沌的必要条件。从功能上看，非线性是通过线性来定义的，设G1和G2是任意两个(向量)函数，a和b是任意两个常数，若算子乙满足如下叠加原理: L(aGl+bG2) =aL(G1)+ bL(G2)，则称L是线性算子，否则L是非线性算子。包含非线性算子的系统称为非线性系统。应当注意的是线性与非线性也不是绝对分明的。对于某些复杂现象，在一定条件下，既可以把它视为非线性现象也可以把它视为线性现象，这与人们看问题的角度和所关心的变量的时空尺度不同有关。现在看来，非线性是普遍存在的，多数问题不能通过线性的办法或线

4、性化的办法来解决，因而直接面对非线性是不可避免的。第5页，共84页，2022年，5月20日，2点13分，星期五3对初始条件的敏感依赖性 1963年，洛伦兹发表了关于混沌理论的开创性研究，并提出了形象的“蝴蝶效应”。被冷落了12年之后，1975年数学家吕埃尔和塔肯斯建议了一种湍流发生机制，认为向湍流的转变是由少数自由度决定的，经过两三次突变，运动就到了维数不高的“奇怪吸引子”上。这里所谓“吸引子”是指运动轨迹经过长时间之后所采取的终极形态：它可能是稳定的平衡点，或周期性的轨道；但也可能是继续不断变化、没有明显规则或次序的许多回转曲线，这时它就称为“奇怪吸引子”。奇怪吸引子上的运动轨道，对轨道初始

5、位置的细小变化极其敏感，但吸引子的大轮廓却是相当稳定的。 第6页，共84页，2022年，5月20日，2点13分，星期五真实球虚拟球第7页，共84页，2022年，5月20日，2点13分，星期五今天，“蝴蝶效应”几乎成了混沌现象的代名词。1961年美国气象学家洛伦兹利用他的一台老爷计算机，根据他导出的描述气象演变的非线性动力学方程进行长期气象预报的模拟数值计算，探讨准确进行长期天气预报的可能性。有一次，洛伦兹为了检验上一次的计算结果，决定再算一遍。但他不是从上一次计算时的最初输入的数据开始验算，而是以一个中间结果作为验算的输入数据。他发现，经过一段重复过程后，计算开始偏离上次的结果，甚至大相径庭。

6、就好比一个计算结果预报几个月后的某天是晴空万里，另一个计算结果则告诉你这一天将电闪雷鸣！第8页，共84页，2022年，5月20日，2点13分，星期五后来洛伦兹发现两次计算的差别只是第二次输入中间数据时将原来的0.506127省略为0.506。洛伦兹意识到，因为他的方程是非线性的，非线性方程不同于线性方程，线性方程对初值的依赖不敏感，而非线性方程对初值的依赖极其敏感。正是初始条件的微小误差导致了计算结果的巨大偏离。由此洛伦兹断言：准确地作出长期天气预报是不可能的。对此，洛伦兹作了个形象的比喻：一只蝴蝶在巴西扇动一下翅膀会在美国的得克萨斯州引起一场龙卷风，这就是蝴蝶效应。第9页，共84页，2022

7、年，5月20日，2点13分，星期五逻辑斯蒂映射的形式为第10页，共84页，2022年，5月20日，2点13分，星期五Example : f(xn+1)=4xn(1-xn)brown: x0=0.6 green: x0=0.6001第11页，共84页，2022年，5月20日，2点13分，星期五Example : f(xn+1)=4xn(1-xn)brown: x0=0.37 green: x0=0.3701第12页，共84页，2022年，5月20日，2点13分，星期五Example : f(xn+1)=4xn(1-xn)i) 系统的变化看似毫无规则，但实际上是有迹可寻的。ii)系统的演化对初始条

8、件的选取非常敏感，初始条件极微小的分别（就例如0.6和0.6001仅仅相差六千分之一）， 在一段时间的演化后可带来南辕北辙的结果。第13页，共84页，2022年，5月20日，2点13分，星期五典型连续混沌系统Chen系统第14页，共84页，2022年，5月20日，2点13分，星期五典型连续混沌系统Lorenz系统第15页，共84页，2022年，5月20日，2点13分，星期五典型连续混沌系统Rssler系统第16页，共84页，2022年，5月20日，2点13分，星期五典型连续混沌系统Chua系统第17页，共84页，2022年，5月20日，2点13分，星期五典型离散混沌映射第18页，共84页，20

9、22年，5月20日，2点13分，星期五典型离散混沌映射第19页，共84页，2022年，5月20日，2点13分，星期五4非周期性 在数学和物理学中，周期性的定义是很明确的。对于函数f(x)，若能找到一个最小正数t满足关系f(x+t)=f(x)，则称f(x)是周期函数，t为其周期；否则f(x)就是非周期的, 非周期性意味着构成奇怪吸引子的积分曲线从不重复原曲线而封闭。这样，向着奇怪吸引子演化的系统，从来不以同样的状态重新经过。非周期性说明，混沌运动的每一瞬间都是“不可预见的创新”的发生器。应当注意的是“非周期性”这个概念比“混沌要广、要大的多。比如，准周期是非周期的，但不是混沌；遍历运动是非周期的

10、，但单纯遍历还不是混沌。混沌运动要求有“混合”的性质，即“对初始条件的敏感依赖性”。但这并不能因此说混沌运动就是杂乱而无用的，相反，混沌不是无序和紊乱。一提到有序，人们往往会想到周期排列或对称形状。第20页，共84页，2022年，5月20日，2点13分，星期五但是，混沌更像是没有周期性的次序。在理想模型中，它可能包含着无穷的内在层次，层次之间存在着“自相似性”或“不尽相似”。在观察手段的分辨率不高时，只能看到某一个层次的结构；提高分辨率之后，在原来不能识别之处又会出现更小尺度上的结构。第21页，共84页，2022年，5月20日，2点13分，星期五 分叉(bifurcation)是有序演化理论的

11、基本概念，这是混沌出现的先兆。在动态系统演化过程中的某些关节点上，系统的定态行为(稳定行为)可能发生定性的突然改变，即原来的稳定定态变为不稳定定态，同时出现新的定态，这种现象就是分叉。发生分叉现象的关节点叫做分叉点，在分叉点系统演化发生质的变化。动态系统演化中的分叉现象充分说明了量变引起质变的规律。分叉又是一种阈值行为，只要系统的非线性作用强到一定程度，就可能出现分叉。所以，凡是产生混沌的系统，总可以观察到分叉序列。5分叉第22页，共84页，2022年，5月20日，2点13分，星期五第23页，共84页，2022年，5月20日，2点13分，星期五第24页，共84页，2022年，5月20日，2点1

12、3分，星期五以参数a为横坐标、以x的稳定定态(stable steady states)为纵坐标作图， 得到1、图2等。从图中可以看出开始是周期加倍分岔(也称周期倍化分岔或周期倍分岔)，然后是混沌，混沌区中又有周期窗口。窗口放大后又可见到同样结构的一套东西。此 所谓无穷自相似结构。第25页，共84页，2022年，5月20日，2点13分，星期五分形性是指奇怪吸引子的结构具有自相似性和不可微性。它不是传统欧几里得几何中描述的直线、平面等整形几何形状所具有的可微性，而是分维的“分形”物，具有结构自相似性和不可微性(不连续性)。目前所发现的奇怪吸引子，如马蹄铁吸引子、洛伦兹吸引子、埃农(Michel

13、Henon)吸引子、若斯勒(Otto ROssler)吸引子等都具有分形性。所以分形并非纯数学抽象的产物，而是对普遍存在的复杂几何形态的科学概括。自然界中分形体无处不在，如起伏蜿蜒的山脉、凹凸不平的地面、曲曲折折的海岸线等等。它与混沌的内随机性、对初始条件的敏感依赖性有本质联系。所以我们说： “混沌本质上是非线性动力系统在一定控制参数范围内产生的对初始条件具有极度敏感依赖性的回复性的非周期性行为状态”。6.分形第26页，共84页，2022年，5月20日，2点13分，星期五分形(fractal)混沌世界的秩序结构：由不断的图形迭代而成利用简单的规则让系统复杂；从复杂不可解的系统中找到简单美妙的秩

14、序。第27页，共84页，2022年，5月20日，2点13分，星期五分形(fractal)混沌世界的秩序古典欧式几何：重视实际可测的量值 例如：长度、深度、厚度 分形：无法单纯用整数维度来描述第28页，共84页，2022年，5月20日，2点13分，星期五分形(fractal)混沌世界的秩序七十年代的数学家畢諾特曼德布洛特(Benoit Mandelbrot)提出一个问题：毛线团的维度是多少？Answer：看你的观点而异 第29页，共84页，2022年，5月20日，2点13分，星期五分形(fractal)混沌世界的秩序远距离來看，线团凝聚成点，维度为零；再近一点，看出来毛线团点据球形的空间，维度扩

15、展成三；再走近一些，看出毛线团是由一根根毛线所构成，他的维度为一，And then ?数据结果视观测者与其对象而改变。这种概念也正是这个世纪物理学的中心思想。毛线的维度?第30页，共84页，2022年，5月20日，2点13分，星期五分形几何的基本思想第31页，共84页，2022年，5月20日，2点13分，星期五欧几里得几何学的研究对象是具有特征长度的几何物体：一维空间：线段，有长度，没有宽度；二维空间：平行四边形，有周长、面积；三维空间：球，表面积、体积；自然界中很多的物体具有特征长度，诸如：人有高度、山有海拔高度等。第32页，共84页，2022年，5月20日，2点13分，星期五有一类问题却比

16、较特别，Mandelbrot就提出了这样一个问题：英国的海岸线有多长？第33页，共84页，2022年，5月20日，2点13分，星期五英国的海岸线地图第34页，共84页，2022年，5月20日，2点13分，星期五当你用一把固定长度的直尺(没有刻度)来测量时，对海岸线上两点间的小于尺子尺寸的曲线，只能用直线来近似。因此，测得的长度是不精确的。如果你用更小的尺子来刻画这些细小之处，就会发现，这些细小之处同样也是无数的曲线近似而成的。随着你不停地缩短你的尺子，你发现的细小曲线就越多，你测得的曲线长度也就越大。如果尺子小到无限，测得的长度也是无限。 第35页，共84页，2022年，5月20日，2点13分

17、，星期五得到的结论是：海岸线的长度是多少：决定与尺子的长短。海岸线的长度是无限的！而显然海岸线的面积为零；而我们确实看到了海岸线的存在，而且海岸线应该是有界的。海岸线什么有界？（长度、面积、体积显然无界）。第36页，共84页，2022年，5月20日，2点13分，星期五Koch 曲线第37页，共84页，2022年，5月20日，2点13分，星期五天空中的云朵植物的叶子第38页，共84页，2022年，5月20日，2点13分，星期五自然界中的分形山星 云第39页，共84页，2022年，5月20日，2点13分，星期五星 云第40页，共84页，2022年，5月20日，2点13分，星期五二、混沌在通信中的应

18、用第41页，共84页，2022年，5月20日，2点13分，星期五混沌同步混沌系统的同步是指一个系统的混沌轨道收敛于另一个混沌系统的轨道，它们之间步调一致。第42页，共84页，2022年，5月20日，2点13分，星期五第43页，共84页，2022年，5月20日，2点13分，星期五第44页，共84页，2022年，5月20日，2点13分，星期五第45页，共84页，2022年，5月20日，2点13分，星期五第46页，共84页，2022年，5月20日，2点13分，星期五第47页，共84页，2022年，5月20日，2点13分，星期五第48页，共84页，2022年，5月20日，2点13分，星期五第49页，共

19、84页，2022年，5月20日，2点13分，星期五第50页，共84页，2022年，5月20日，2点13分，星期五Shannon“好的混合变换通常是两个简单的非可交换运算的乘积。比如Hopf已经证明，如做馅皮的生面团可以通过下面的一系列操作进行混合：面团首先被揉搓成一个扁面皮，然后将它折叠，再搓揉，再折叠，如此往复。一个混合变换中的函数应该是复杂的，它的所有变量都应敏感，对任何一个变量来说，一个很小的变化都应引起输出的显著不同。”第51页，共84页，2022年，5月20日，2点13分，星期五混沌与密码学的关系随机密钥流产生器混沌与流密码学流密码的核心产生不可预测的混沌序列 混沌第52页，共84页

20、，2022年，5月20日，2点13分，星期五混沌与密码学的关系混沌与分组密码学混淆和扩散 分组密码 混沌对密钥敏感 对明文敏感增加信源的熵反复压缩和拉伸的混沌变换混沌对初始条件和参数的敏感 混沌具有遍历性的性质 混沌具有拓扑传递性 第53页，共84页，2022年，5月20日，2点13分，星期五混沌与密码学的关系单向函数混沌与公钥密码学 公钥密码已知部分结构重构出全部高维混沌系统；未知部分参数同步两个超混沌系统或时空混沌系统； 混沌第54页，共84页，2022年，5月20日，2点13分，星期五3.混沌掩盖第55页，共84页，2022年，5月20日，2点13分，星期五混沌掩盖第56页，共84页，2

21、022年，5月20日，2点13分，星期五4.混沌开关第57页，共84页，2022年，5月20日，2点13分，星期五混沌开关第58页，共84页，2022年，5月20日，2点13分，星期五5.混沌调制第59页，共84页，2022年，5月20日，2点13分，星期五混沌调制第60页，共84页，2022年，5月20日，2点13分，星期五混沌在图像加密中的应用原图像加密后图像第61页，共84页，2022年，5月20日，2点13分，星期五混沌在图像加密中的应用正确解密图像错误解密图像第62页，共84页，2022年，5月20日，2点13分，星期五信息伪装：第63页，共84页，2022年，5月20日，2点13分

22、，星期五作业第64页，共84页，2022年，5月20日，2点13分，星期五一、函数迭代 给定一函数 以及初始点 ，定义数列 称为函数 的迭代序列。 满足 的点 称为 的不动点，记之为 。如果所有附近的点在迭代过程中都趋向于某一不动点，则该不动点称为吸引点。如果所有附近的第65页，共84页，2022年，5月20日，2点13分，星期五点都远离它，则它是排斥点。例如，0 与 1 是 的不动点。0 是吸引点，1是排斥点。如果则点集 形成一个 k 循环。 称为 k 周期点。k称为周期。第66页，共84页，2022年，5月20日，2点13分，星期五类似地，周期点也可以分吸引点与排斥点。如果点 最终归宿于某

23、个循环中，则称它为预周期点。如 1 是 的预周期点。迭代序列 的收敛与发散性质不仅与函数 有关， 而且与初值的选择有关。 例如，对于迭代第67页，共84页，2022年，5月20日，2点13分，星期五 当初值 时， 迭代序列收敛，否则发散。第68页，共84页，2022年，5月20日，2点13分，星期五二、二次函数的迭代对二次函数 做迭代: 迭代的几何直观图第69页，共84页，2022年，5月20日，2点13分，星期五练习 1 对几组不同的参数值 （如 )以及不同的初值 ，观察迭代是否收敛。练习 2 取参数 ，用不同的初值做迭代。你能找到一个吸引的不动点吗？一个排斥的不动点吗？哪些初值收敛到吸引的

24、不动点？哪些初值使序列发散？取不动的参数 回答同样的问题。第70页，共84页，2022年，5月20日，2点13分，星期五练习 3 找出一个参数 使它对应的迭代具有2周期点。这种性质依赖于初值吗？练习 4 对任意的整数 ，你能找到一个 值使得它对应的迭代具有 周期点吗？对哪些 值能给出 周期点？在每种情况下，结果是否依赖于初值？（对 和 的值进行验证）第71页，共84页，2022年，5月20日，2点13分，星期五练习 5 如果某个值能给出周期点，它是否一定是吸引的周期点？你能否找到排斥的周期点？练习 6 根据前面的练习，试着从理论上分析：如何求不动点？对哪些值对应吸引的不动点？哪些值对应排斥的不

25、动点？初值对结果有什么影响？对周期点做类似的分析。第72页，共84页，2022年，5月20日，2点13分，星期五不动点的计算 从 得到 及第73页，共84页，2022年，5月20日，2点13分，星期五 吸引的不动点与排斥的不动点 定理设 是 的不动点，如果在 附近有 ，则 是 的吸引的不动点；否则， 是 的排斥的不动点。由于 故当 0a1时，为吸引点，(a-1)/a为排斥点。当1a3, 为排斥点，(a-1)/ a为吸引点。第74页，共84页，2022年，5月20日，2点13分，星期五 2 周期点 得第75页，共84页，2022年，5月20日，2点13分，星期五三、Feigenbaum图将区间（0, 4 以某个步长 （如 ）离散化。对每个离散的 值做迭代。忽略前50个迭代值，而把点 显示在坐标平面上，最后形成的图形称为 Feigenbaum图。第76页，共84页，2022年，5月20日，2点13分，星期五练习 7 观察Feigenbaum图。（1）它的左部有一条曲线，这表示什么意义？（2） 从某一点 开始，这条曲线分成两支，这说明了迭代的什么