《混沌学》课件.pptVIP

* * 从Malthus模型 上海交通大学数学系 到浑沌 数学实验 巴西的一只蝴蝶扇动翅膀会引起 — 洛伦兹 数学的伟大使命在于从混沌中 发现秩序。 — 倍尔 在那个混沌的体制中，结构上的微 小差异几乎都会造成行为方式上的巨大 变化，可控制的行为似乎已被排除。 — 斯图尔特.考夫曼 明年在得克萨斯的大风暴吗？ ● 函数的迭代，不动点和有关的作图 ● 介绍浑沌，用数值迭代、蛛网迭代和密度 ● 浑沌的倍周期分叉、遍历性和某些 ● 计算机与科学研究(即使是数学) 分布等方法来研究浑沌 普适结构 实验将告诉我们什么？ 问题的提出 出现在各个领域的一种现象:数学、物理、 由此引起的复杂而有趣的现象 ? “侏罗纪公园”中的恐龙重现 ● 我们的讨论 ● 什么是浑沌? 生物、金融、经济、管理等等： ? 宇宙的起源 ? 龙卷风的产生、厄尔尼诺现象 ? 东南亚金融危机爆发 从某些简单的离散的数学模型开始,进而讨论 数学模型 ● Malthus 模型 口数成正比,从而 xn+1＝ xn ＋r xn xn+1 = a xn 其中 a=r+1. 设xn是某人类群体在第n个时间段(例如年)末 时的总数，若在单位时间段内人口相对增长率为r （出生率与死亡率之差）,那么人口增长数与原人 即是 记 g (x) = a x，则是函数迭代 xn= a xn－1= a2xn－2 =…= an x0 于是Malthus的结论：人口增长呈几何级数 约35年增加一倍，与1700－1961年世界人口 与近年统计结果有误差，由a 1,xn趋向无穷， xn＋1＝ g( xn ) 容易得到 统计结果一致 模型在人口长期预测方面必定是失效的. . ● Logistic模型 生存资源是重要的因素，修改模型为： xn+1 － xn= r xn－ b xn2 － b xn2为竞争(约束)项，r、b 称生命系数,则 xn+1= a xn- bxn2 ， （a=r+1） 这是一个如下非线性映射的迭代 f1(x)= ax - bx2 数据观察（利用Mathematica） In[1]:= a=1.029; b=1.48654*10^(-11); f1[x_]:=a*x-bx^2; For[n=1979; x[1979]=9.7542*10^8, n=2002, n++, x[n+1]=f1[x[n]]; Print[n+1, “ ”，x[n+1]/10^8]] 与统计数字接近 9.89564 10.037 10.1784 10.3195 10.4605 10.6012 10.7416 10.8815 11.0211 11.1601 11.2986 11.4365 11.5738 11.7103 1994 11.8460 1995 11.9809 1996 12.1150 1997 12.2482 1998 12.3803 1999 12.5115 2000 12.6417 2001 12.7701 2002 12.8986 2003 13.0245 了解混沌 ● Logistic映射 (Robert.May的研究） f(x)= a x(1- x)， x 在[0,1]内变化 xn+1= f(xn) 从[0,1]内点x0出发，由Logistic映射的迭代形成 xn= f n(x0), n = 0,1,2,… 序列{xn}称为x0的轨道 种群数的模型简化： 相应的迭代为 了一个序列，即 ● 数值迭代（ a 逐渐增加，迭代会有何结果） 1．倍周期分叉现象 ■ 当0a 1时，由于0xnaxn+1 ■ 当1a3时，任何(0,1)中初始值的轨道趋于 两个不动点x1*, x2* ,一个稳定(吸引),另一个 xn →0 物种逐渐灭亡 x*=1-1/a 其中x*是方程f(x)=x的解，为映射f 的不动点 （周期1点)例：a =1.5时 xn → 1/3. 不稳定，轨道{xn}趋向稳定点 这两个数满足 ■ 当3a1+61/2时， xn 绕着两个数 x3*，x4*振动, x2k-1 →0.799455 x2k →o.513045 ■ 当1+61/2a3.5440903506…时, 从任意的点x