Matlab画Lorenz系统的最大李雅普诺夫指数图学习资料
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
grid on;
4、“体验化”消费
除了“漂亮女生”形成的价格，优惠等条件的威胁外，还有“碧芝”的物品的新颖性，创意的独特性等，我们必须充分预见到。结果图
与此同时，上海市工商行政管理局也对大学生创业采取了政策倾斜：凡高校毕业生从事个体经营的，自批准经营日起，１年内免交登记注册费、个体户管理费、集贸市场管理费、经济合同鉴证费、经济合同示范文本工本费等，但此项优惠不适用于建筑、娱乐和广告等行业。
（2）文化优势plot(b,Z,'-');
此次调查以女生为主，男生只占很少比例，调查发现58％的学生月生活费基本在400元左右，其具体分布如（图1-1）title('JD_{1}系统最大lyapunov指数')
加拿大ｂｅａｄｗｏｒｋｓ公司就是根据年轻女性要充分展现自己个性的需求，将世界各地的珠类饰品汇集于“碧芝自制饰品店”内，由消费者自选、自组、自制，这样就能在每个消费者亲手制作、充分发挥她们的艺术想像力的基础上，创作出作品，达到展现个性的效果。xlabel('parameter b'),ylabel('The largest Lyapunov exponents');
x1=x+(d0/d1)*(x1-x);
y1=y+(d0/d1)*(y1-y);
z1=z+(d0/d1)*(z1-z);
if i>50
lsum=lsum+log(d1/d0);
当然，在竞争日益激烈的现代社会中，创业是件相当困难的事。我们认为，在实行我们的创业计划之前，我们首先要了解竞争对手，吸取别人的经验教训，制订相应竞争的策略。我相信只要我们的小店有自己独到的风格，价格优惠，服务热情周到，就一定能取得大多女孩的信任和喜爱。end
300-400元1632%end
Z=[Z lsum/(i-50)];
（3）年龄优势end
据调查统计，有近94%的人喜欢亲戚朋友送给自己一件手工艺品。无论是送人，个人兴趣，装饰还是想学手艺，DIY手工制作都能满足你的需求。下表反映了同学们购买手工艺制品的目的。如图（1-4）b=linspace(0,500,500);
Z=[];
global a;
global b;
global c;
a=10;
c=8/3;
d0=1e-7;for b=lin来自pace(0,500,500)
lsum=0;
x=1;y=1;z=1;
x1=1;y1=1;z1=1+d0;
for i=1:100
[T1,Y1]=ode45('Loren',1,[x;y;z;16;b;4]);
global b;
global c;
x=X(1); y=X(2); z=X(3);
% Y的三个列向量为相互正交的单位向量
%输出向量的初始化
dX = zeros(6,1);
% Lorenz吸引子
dX(1)=a*(y-x);
dX(2)=x*(b-z)-y;
dX(3)=x*y-c*z;
end
计算最大L指数文件
Lorenz系统
文档分两个文件方程m文件和计算L指数m文件分开写，复制粘贴即可运行matlab2012a，改写方程文件和参数即可算自己的系统，其中最大L指数用的是经典的柏内庭（G.Benettin）计算方法，准确快速无误！附计算结果图！！
方程m文件：
function dX = Loren(t,X)
global a; %变量不放入参数表中
[T2,Y2]=ode45('Loren',1,[x1;y1;z1;16;b;4]);
n1=length(Y1);n2=length(Y2);
x=Y1(n1,1);y=Y1(n1,2);z=Y1(n1,3);
x1=Y2(n2,1);y1=Y2(n2,2);z1=Y2(n2,3);
d1=sqrt((x-x1)^2+(y-y1)^2+(z-z1)^2);
4、“体验化”消费
除了“漂亮女生”形成的价格，优惠等条件的威胁外，还有“碧芝”的物品的新颖性，创意的独特性等，我们必须充分预见到。结果图
与此同时，上海市工商行政管理局也对大学生创业采取了政策倾斜：凡高校毕业生从事个体经营的，自批准经营日起，１年内免交登记注册费、个体户管理费、集贸市场管理费、经济合同鉴证费、经济合同示范文本工本费等，但此项优惠不适用于建筑、娱乐和广告等行业。
（2）文化优势plot(b,Z,'-');
此次调查以女生为主，男生只占很少比例，调查发现58％的学生月生活费基本在400元左右，其具体分布如（图1-1）title('JD_{1}系统最大lyapunov指数')
加拿大ｂｅａｄｗｏｒｋｓ公司就是根据年轻女性要充分展现自己个性的需求，将世界各地的珠类饰品汇集于“碧芝自制饰品店”内，由消费者自选、自组、自制，这样就能在每个消费者亲手制作、充分发挥她们的艺术想像力的基础上，创作出作品，达到展现个性的效果。xlabel('parameter b'),ylabel('The largest Lyapunov exponents');
x1=x+(d0/d1)*(x1-x);
y1=y+(d0/d1)*(y1-y);
z1=z+(d0/d1)*(z1-z);
if i>50
lsum=lsum+log(d1/d0);
当然，在竞争日益激烈的现代社会中，创业是件相当困难的事。我们认为，在实行我们的创业计划之前，我们首先要了解竞争对手，吸取别人的经验教训，制订相应竞争的策略。我相信只要我们的小店有自己独到的风格，价格优惠，服务热情周到，就一定能取得大多女孩的信任和喜爱。end
300-400元1632%end
Z=[Z lsum/(i-50)];
（3）年龄优势end
据调查统计，有近94%的人喜欢亲戚朋友送给自己一件手工艺品。无论是送人，个人兴趣，装饰还是想学手艺，DIY手工制作都能满足你的需求。下表反映了同学们购买手工艺制品的目的。如图（1-4）b=linspace(0,500,500);
Z=[];
global a;
global b;
global c;
a=10;
c=8/3;
d0=1e-7;for b=lin来自pace(0,500,500)
lsum=0;
x=1;y=1;z=1;
x1=1;y1=1;z1=1+d0;
for i=1:100
[T1,Y1]=ode45('Loren',1,[x;y;z;16;b;4]);
global b;
global c;
x=X(1); y=X(2); z=X(3);
% Y的三个列向量为相互正交的单位向量
%输出向量的初始化
dX = zeros(6,1);
% Lorenz吸引子
dX(1)=a*(y-x);
dX(2)=x*(b-z)-y;
dX(3)=x*y-c*z;
end
计算最大L指数文件
Lorenz系统
文档分两个文件方程m文件和计算L指数m文件分开写，复制粘贴即可运行matlab2012a，改写方程文件和参数即可算自己的系统，其中最大L指数用的是经典的柏内庭（G.Benettin）计算方法，准确快速无误！附计算结果图！！
方程m文件：
function dX = Loren(t,X)
global a; %变量不放入参数表中
[T2,Y2]=ode45('Loren',1,[x1;y1;z1;16;b;4]);
n1=length(Y1);n2=length(Y2);
x=Y1(n1,1);y=Y1(n1,2);z=Y1(n1,3);
x1=Y2(n2,1);y1=Y2(n2,2);z1=Y2(n2,3);
d1=sqrt((x-x1)^2+(y-y1)^2+(z-z1)^2);