• / 2
  • 下载费用:10 金币  

混沌系统分岔图.doc

关 键 词:
混沌系统分岔图.doc
资源描述:
CHEN系统求解分岔图一:function dy=chen(t,y)global r;dy=zeros(3,1);dy(1)=35*(y(2)-y(1));dy(2)=-7*y(1)+28*y(2)-y(1)*y(3);dy(3)=y(1)*y(2)-r*y(3)end对 r的分岔图:function [Xmax] = getmax(y) a=length(y); j=1; for i=(a-1)/2:a b=(y(i,1)-y(i-2,1))/2; c=(y(i,1)+y(i-2,1))/2-y(i-1,1); if y(i-2,1)=y(i,1) elseif y(i-2,1)=y(i,1) Xmax(j)=y(i-1,1)-b^2/(4*c); j=j+1; end end求解:clear all global r t0=[0 100]; y0=[0.015,0.758,1.83]; for r=1:5.5:1000 [t,y]=ode45('FUN',t0,y0); [Xmax]=getmax(y(:,1)); plot(r,Xmax,'b','markersize',1) hold on clear Xmax end xlabel('r') ylabel('Xmax')求解分岔图二:clc;clearglobal r;Z=[];for r=linspace(0,0.01,100)options=odeset(’RelTol’,1e-6,’AbsTol’,[1e-4 1e-4 1e-5]);[T,Y]=ode45(’FUN’,[0,50],[0.015 0.758 1.83],options);n=length(Y);for k=round(n/2):nif abs(Y(k,1))1Z=[Z,r+abs(Y(k,2))*i];endendendfigureplot(Z,’.’,’markersize’,1)title(’chen映射分岔图’)xlabel(’r’),ylabel(’y’)
展开阅读全文
  微传网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:混沌系统分岔图.doc
链接地址:https://www.weizhuannet.com/p-9786368.html
微传网是一个办公文档、学习资料下载的在线文档分享平台!

微传网博客

网站资源均来自网络,如有侵权,请联系客服删除!

 网站客服QQ:80879498  会员QQ群:727456886

copyright@ 2018-2028 微传网络工作室版权所有

     经营许可证编号:冀ICP备18006529号-1 ,公安局备案号:13028102000124

收起
展开