Matlab非线性方程与混沌

1、2021/4/21 一只蝴蝶在巴西扇动一只蝴蝶在巴西扇动翅膀，有可能在美国的德克翅膀，有可能在美国的德克萨斯引起一场龙卷风吗？萨斯引起一场龙卷风吗？2021/4/22Logistic方程与混沌方程与混沌n在生物学中，有一个刻画生物种群个体总量增长情在生物学中，有一个刻画生物种群个体总量增长情况的著名的方程况的著名的方程Logistic方程：方程：为为比比例例系系数数rxrxxnnn)1(1 n其中其中xn为某生物群体的第为某生物群体的第n代的个体总数与该群代的个体总数与该群体所能达到的最大保有量时的个体数之比。体所能达到的最大保有量时的个体数之比。n选定初值和比例系数选定初值和比例系数r的值后

2、，由方程就能生成一的值后，由方程就能生成一个数列：个数列：nxxx,212021/4/23 考察迭代格式考察迭代格式(Logistic方程方程 )1(1nnnxrxx 初值初值1 . 00 x1. 当参数当参数r取值分别为取值分别为1.2，3. 6时，时，考察其迭代序列的收敛情况考察其迭代序列的收敛情况2021/4/24clc;clf;x=0.1; y= ;r=1.2; % %改变取值得到相应的图形改变取值得到相应的图形hold onaxis(0 100 0 1)for k=1:100 x=r*x*(1-x);y=y,x; plot(k,x,k.,markersize,10) fprintf(

3、x(%d)=%.10fn,k,x)endt=1:100;plot(t,y,k-);grid程序程序010203040506070809010000.10.20.30.40.50.60.70.80.912021/4/25将参数将参数r取为取为3.6时迭代序列的收敛情况：时迭代序列的收敛情况：010203040506070809010000.10.20.30.40.50.60.70.80.912021/4/26clear;clf;hold onaxis(0,4,0,1);gridfor r=0:0.3:3.9 x=0.1; for k=2:150 x(k)=r*x(k-1)*(1-x(k-1);

4、end pause(0.05) for k=101:150 plot(r,x(k),k.); end text(r-0.1,max(x(101:150)+0.05,itr=,num2str(r)end程序程序2. 将参数将参数r取取0，0.3，0.6，0.9，1.2，3.9的迭代序列的迭代序列收敛情况放置到同一坐标系中观察其变化收敛情况放置到同一坐标系中观察其变化2021/4/2700.511.522.533.5400.10.20.30.40.50.60.70.80.91r=0r=0.3r=0.6 r=0.9r=1.2r=1.5r=1.8r=2.1r=2.4r=2.7r=3r=3.3r=3.6

5、r=3.9参 数 r迭 代 序 列 x2021/4/28clear;clf; axis(2.7,4,0,1);gridhold onfor r=2.7:0.005:3.9 x=0.1; for i=2:150 x(i)=r*x(i-1)*(1-x(i-1); end pause(0.1) fprintf(r=%.3fn,r) for i=101:150 plot(r,x(i),k.); endend程序程序3.现在对现在对r r取值在取值在2.72.7到到3.93.9之间进行加密迭代并作图，之间进行加密迭代并作图， 取步长为取步长为0.0050.005时时 2021/4/29参数参数r的微小变

6、化引起结果巨大的变化的微小变化引起结果巨大的变化2021/4/210请同学们再次加密请同学们再次加密r取值进行实验，回答取值进行实验，回答下面问题下面问题n（）是否由支分叉为支，并依次类（）是否由支分叉为支，并依次类推呢？推呢？n（）这些分叉点处（）这些分叉点处r的取值，是否有规律？的取值，是否有规律？混沌现象混沌现象2021/4/211n混沌混沌（译自英文（译自英文Chaos）的原意是指无序和混乱）的原意是指无序和混乱的状态。这些表面上看起来无规律、不可预测的的状态。这些表面上看起来无规律、不可预测的现象，实际上有它自己的规律。现象，实际上有它自己的规律。n 混沌学的任务：就是寻求混沌现象的

7、规律，加以混沌学的任务：就是寻求混沌现象的规律，加以处理和应用。处理和应用。n 60年代混沌学的研究热悄然兴起，渗透到物理学、年代混沌学的研究热悄然兴起，渗透到物理学、化学、生物学、生态学、力学、气象学、经济学、化学、生物学、生态学、力学、气象学、经济学、社会学等诸多领域，成为一门新兴学科。社会学等诸多领域，成为一门新兴学科。 什么是混沌呢？什么是混沌呢？ 2021/4/212n科学家给混沌下的定义是：混沌是指发生在确定性科学家给混沌下的定义是：混沌是指发生在确定性系统中的貌似随机的不规则运动，一个确定性理论系统中的貌似随机的不规则运动，一个确定性理论描述的系统，其行为却表现为不确定性一不可重

8、复、描述的系统，其行为却表现为不确定性一不可重复、不可预测，这就是混沌现象。不可预测，这就是混沌现象。n进一步研究表明，混沌是非线性动力系统的固有特进一步研究表明，混沌是非线性动力系统的固有特性，是非线性系统普遍存在的现象。牛顿确定性理性，是非线性系统普遍存在的现象。牛顿确定性理论能够充美处理的多为线性系统，而线性系统大多论能够充美处理的多为线性系统，而线性系统大多是由非线性系统简化来的。因此，在现实生活和实是由非线性系统简化来的。因此，在现实生活和实际工程技术问题中，混沌是无处不在的！际工程技术问题中，混沌是无处不在的！ 2021/4/213混沌的特征混沌的特征 1.差之毫厘，失之千里、牵一

9、发而动全身。差之毫厘，失之千里、牵一发而动全身。n一个小小初始条件的差异可以严重影响系统一个小小初始条件的差异可以严重影响系统长期的大变化。长期的大变化。 n2.对初始条件的敏感性。对初始条件的敏感性。n对原本西方的科学基本理念来说，如果你对原本西方的科学基本理念来说，如果你正在计算台面上的一颗撞球，就不用去理会正在计算台面上的一颗撞球，就不用去理会室外一片树叶的掉落。很轻微的影响可以被室外一片树叶的掉落。很轻微的影响可以被忽略，事物进行总会殊途同归，任意的小干忽略，事物进行总会殊途同归，任意的小干扰，并不致于膨胀到任意大的后果。扰，并不致于膨胀到任意大的后果。2021/4/214n 1960

10、年，美国麻省理工学院教授洛伦兹研究年，美国麻省理工学院教授洛伦兹研究“长期长期天气预报天气预报”问题时，在计算机上用一组简化模型模拟问题时，在计算机上用一组简化模型模拟天气的演变。他原本的意图是利用计算机的高速运算天气的演变。他原本的意图是利用计算机的高速运算来提高技期天气预报的准确性。但是，事与愿违，多来提高技期天气预报的准确性。但是，事与愿违，多次计算表明，初始条件的极微小差异，均会导致计算次计算表明，初始条件的极微小差异，均会导致计算结果的很大不同。结果的很大不同。 n 由于气候变化是十分复杂的，所以在预测天气时，由于气候变化是十分复杂的，所以在预测天气时，输入的初始条件不可能包含所有的

11、影响因素（通常的输入的初始条件不可能包含所有的影响因素（通常的简化方法是忽略次要因素，保留主要因素），而那些简化方法是忽略次要因素，保留主要因素），而那些被忽略的次要因素却可能对预报结果产生重大影响，被忽略的次要因素却可能对预报结果产生重大影响，导致错误的结论。由此，洛伦兹认定，尽管拥有高速导致错误的结论。由此，洛伦兹认定，尽管拥有高速计算机和精确的测量数据（温度、风速、气压等），计算机和精确的测量数据（温度、风速、气压等），也难以获得准确的长期天气预报。也难以获得准确的长期天气预报。 2021/4/215蝴蝶效应蝴蝶效应 n 19791979年年1212月，洛伦兹在华盛顿的美国科学促进会月，

12、洛伦兹在华盛顿的美国科学促进会的一次讲演中提出：一只蝴蝶在巴西扇动翅膀，的一次讲演中提出：一只蝴蝶在巴西扇动翅膀，有可能会在美国的德克萨斯引起一场龙卷风。他有可能会在美国的德克萨斯引起一场龙卷风。他的演讲和结论给人们留下了极其深刻的印象。从的演讲和结论给人们留下了极其深刻的印象。从此以后，所谓此以后，所谓“蝴蝶效应蝴蝶效应”之说就不胫而走，名之说就不胫而走，名声远扬了。声远扬了。 从科学的角度来看，从科学的角度来看，“蝴蝶效应蝴蝶效应”反映了混沌反映了混沌运动的一个重要特征：系统的长期行为对初始条运动的一个重要特征：系统的长期行为对初始条件的敏感依赖性。件的敏感依赖性。 2021/4/216n

13、 经典动力学的传统观点认为：系统的长期行为经典动力学的传统观点认为：系统的长期行为对初始条件是不敏感的，即初始条件的微小变化对初始条件是不敏感的，即初始条件的微小变化对未来状态所造成的差别也是很微小的。可混沌对未来状态所造成的差别也是很微小的。可混沌理论向传统观点提出了挑战。混沌理论认为在混理论向传统观点提出了挑战。混沌理论认为在混沌系统中，初始条件的十分微小的变化经过不断沌系统中，初始条件的十分微小的变化经过不断放大，对其未来状态会造成极其巨大的差别。放大，对其未来状态会造成极其巨大的差别。 2021/4/217一则西方寓言：一则西方寓言：丢失一个钉子，坏了一只蹄铁；丢失一个钉子，坏了一只蹄铁； 坏了一只蹄铁，折了一匹战马；坏了一只蹄铁，折了一匹战马；折了一匹战马，伤了一位骑士；折了一匹战马，伤了一位骑士；输了一场战斗，亡了一个帝国。输了一场战斗，亡了一个帝国。 伤了一位骑士，输了一场战斗；伤了一位骑士，输了一场战斗； 马蹄铁上一个钉子是否会丢失，本是初始条件的马蹄铁上一个钉子是否会丢失，本是初始条件的十分微小