matlab求最大李雅普诺夫(Lyapunov)指数程序.docVIP

下载本文档

515
0
约2.01千字
约 3页
2022-03-31 发布于浙江
举报
版权申诉

matlab求最大李雅普诺夫(Lyapunov)指数程序.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

Z1=[];Z2=[];Z3=[];global a;global b;global c;b=3;c=28;for a=linspace(32,40,100); y=[1;1;1;1;0;0;0;1;0;0;0;1]; lp=0; for k=1:200 [T,Y] = ode45(Chen2, 1, y); y = Y(size(Y,1),:); y0 = [y(4) y(7) y(10); y(5) y(8) y(11); y(6) y(9) y(12)]; y0=GS(y0); mod(1)=norm(y0(:,1)); mod(2)=norm(y0(:,2)); mod(3)=norm(y0(:,3)); lp = lp+log(abs(mod)); y0(:,1)=y0(:,1)/mod(1); y0(:,2)=y0(:,2)/mod(2); y0(:,3)=y0(:,3)/mod(3); y(4:12) = y0; end lp=lp/200; Z1=[Z1 lp(1)]; Z2=[Z2 lp(2)]; Z3=[Z3 lp(3)]; enda=linspace(32,40,100);plot(a,Z1,-,a,Z2,-,a,Z3,-);title(Lyapunov exponents of Chen)xlabel(b=3,c=28,parameter a),ylabel(lyapunov exponents)grid on 以上是三个变量的Lyapunov指数谱，下面是最大的Lyapunov指数谱： Z=[];d0=1e-8;for a=linspace(32,40,80)lsum=0;x=1;y=1;z=1;x1=1;y1=1;z1=1+d0;for i=1:100 [T1,Y1]=ode45(Chen,1,[x;y;z;a;3;28]); [T2,Y2]=ode45(Chen,1,[x1;y1;z1;a;3;28]); 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); x1=x+(d0/d1)*(x1-x); y1=y+(d0/d1)*(y1-y); z1=z+(d0/d1)*(z1-z); if i50 lsum=lsum+log(d1/d0); endendZ=[Z lsum/(i-50)];enda=linspace(32,40,80);plot(a,Z,-);title(Chen 系统最大lyapunov指数)xlabel(parameter a),ylabel(lyapunov exponents) 二模拟 Lorenz 系统最大lyapunov指数谱 function ly=jose_ly(b,k)% the largest lyapunov exponent of josephson% k 迭代步数，b 参数% 方程如下:% θ+G*θ+sinθ=I+A*sin(ωt)+αsin(βωt)% 变化：% dx=y% dy=-G*y-sin(x)+I+A*sin(w*t)+a*sin(b*w*t)%% Example:% ly=jose_ly(0,800)%% Author:LDYU% Authors email: ustb03-07@yahoo %d0=1e-8;ly=0;lsum=0;x=[0;2;b];x1=[d0;2;b];for t=1:k [T1,Y1]=ode45(Josephon,[t-1,t],x); [T2,Y2]=ode45(Josephon,[t-1,t],x1); x=Y1(end,:); x1=Y2(end,:); d1=norm(x-x1); x1=x+(d0/d1)*(x1-x); lsum=lsum+log(d1/d0);endly=lsum/k;