周期分岔与混沌现象ppt课件_第1页
周期分岔与混沌现象ppt课件_第2页
周期分岔与混沌现象ppt课件_第3页
周期分岔与混沌现象ppt课件_第4页
周期分岔与混沌现象ppt课件_第5页
已阅读5页,还剩10页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、1周期分岔与混沌现象2主要内容 又称logistic回归分析,主要在流行病学中应用较多,比较常用的情形是探索某疾病的危险因素,根据危险因素预测某疾病发生的概率,等等。例如,想探讨胃癌发生的危险因素,可以选择两组人群,一组是胃癌组,一组是非胃癌组,两组人群肯定有不同的体征和生活方式等。这里的因变量就是是否胃癌,即“是”或“否”,为两分类变量,自变量就可以包括很多了,例如年龄、性别、饮食习惯、幽门螺杆菌感染等。自变量既可以是连续的,也可以是分类的。通过logistic回归分析,就可以大致了解到底哪些因素是胃癌的危险因素。 Logistic模型简介:Logistic模型是一个描述自然界生物种群中昆虫

2、数目变化的数学模型,它由上一年的昆虫数目 及增殖系数 预测次年的昆虫数目 其数学表达式为)(21nnnxxx) 10 , 40(x3初值相同,增殖系数不同时的结果4改变初值的情况5费根保姆图 为了研究增殖系数对分叉现象的影响,可以取相同的初值对所有的增殖系数进行计算。2.62.833.23.43.63.8400.10.20.30.40.50.60.70.80.916费根保姆图作图思路 方法是给定x的初值,对不同的u值计算新的x,共循环250次,循环计算150次后开始画图。计算所得的x值是矩阵,行标对应循环次数,列标对应u值。使用矢量化编程以后,对所有的u值同时计算一次新的x值,得到矩阵x中的一

3、列元素。所得图形即为费根鲍曼图。7两个做费根保姆图的程序u=2.6:0.001:4;x=0.6;for j=1:150, x=u.*(x-x.2);endfor i=1:100 x=u.*(x-x.2); plot(u,x,r.) hold on;end 不保留旧的不保留旧的x值,而是直值,而是直接用它作图,能节省内接用它作图,能节省内存存u=2.6:0.001:4;X=ones(250,1401);X(1,:)=0.6*X(1,:);for j=1:250, X(j+1,:)=u.*(X(j,:)-X(j,:).2);Endplot(u,X(150:end,:),r.)保留所有的X值,每次计

4、算的X值生成矩阵的一行元素,最后作图,程序可读性强,但占用内存较大891011初值的影响混沌现象有个特点,就是初值的微小变化将引起结果的完全不同,它说明混沌现象的不可预测性。0102030405060700.10.20.30.40.50.60.70.80.9112费根鲍姆常数利用分岔点的u值,计算相邻分岔点的间距之比,所得的极限值叫做费根鲍姆常数,这是一个普适常数,所有系统通过倍周期分岔进入混沌时,都会遵循这个规律。669201661. 4lim11nnnnnF13李雅普诺夫指数对于稳定的周期n,有李雅普诺夫指数小于0,对于倍周期分岔点李雅普诺夫指数等于0,混沌状态李雅普诺夫指数大于0,所以李雅普诺夫指数由负变正表明运动向混沌转变。14x=0.7;for j=1:4 subplot(2,2,j)u=input(输入增殖系数系数u=)for i=1:10 x=u*(x-x2); a=i,i+1 b(1,1)=x x2=u*(x-x2) b(1,2)=x2plot(a,b,r:*) hold onendfor i=11:60 x=u*(x-x2); plot(i,x,r:*) hold onendend15for j=1:2 x=input(输入初值=)u=3.8for i=1:60; x=u*(x-x.2); a=i,i+1; b(1,1)=x; x2=u*(x-x.2); b(

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论