求最大利亚普洛夫指数的matlab程序.docVIP

下载本文档

1
0
约7.2千字
约 6页
2018-05-17 发布于河南
举报
版权申诉

求最大利亚普洛夫指数的matlab程序.doc

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

求最大利亚普洛夫指数的matlab程序

程序【1】 function dx=Lorenz(t,x);dx(1,1)=10*(x(2)-x(1));dx(2,1)=x(1)*(30-x(3))-x(2);dx(3,1)=x(1)*x(2)-8/3*x(3);dx(4,1)=0;dx(5,1)=0;dx(6,1)=0;function lambda_1=lyapunov_wolf1(data,N,m,tau,P)%??该函数用来计算时间序列的最大Lyapunov 指数--Wolf 方法%??m: 嵌入维数%??tau:时间延迟%??data:时间序列%??N:时间序列长度%??P:时间序列的平均周期,选择演化相点距当前点的位置差，即若当前相点为I，则演化相点只能在|I－J|P的相点中搜寻%??lambda_1:返回最大lyapunov指数值%**************************************************************************%??ode计算整数阶系统的时间序列%******************************************************************delt_t1 = 0.001;t1 = 0:delt_t1:60;[tt1,y1]=ode45(@lorenz,t1,[-1,0,1]);xx1 = y1(:,1);x1 = spline(tt1, xx1, t1);data= x1(20000:10:60000);%采样N=length(data);m=3;tau=11;%*****************************************************% FFT计算平均周期%**********************************************************x=data;xPower=abs(fft(x)).^2;NN=length(xPower);xPower(1)=[];%去除直流分量NN=floor(NN/2);xPower=xPower(1:NN);freq=(1:NN)/NN*0.5;[mP,index]=max(xPower);P=index;%*************************************************************min_point=1??; %要求最少搜索到的点数MAX_CISHU=5 ;??%最大增加搜索范围次数%FLYINGHAWK%? ?求最大、最小和平均相点距离? ? max_d = 0;? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???%最大相点距离? ? min_d = 1.0e+100;? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? %最小相点距离? ? avg_dd = 0;? ? Y=reconstitution(data,N,m,tau);? ?? ?? ?? ?? ?? ???%相空间重构? ? M=N-(m-1)*tau;? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? %重构相空间中相点的个数? ? for i = 1 : (M-1)? ?? ???for j = i+1 : M? ?? ?? ?? ?d = 0;? ?? ?? ?? ?for k = 1 : m? ?? ?? ?? ?? ? d = d + (Y(k,i)-Y(k,j))*(Y(k,i)-Y(k,j));? ?? ?? ?? ?end? ?? ?? ?? ?d = sqrt(d);? ?? ?? ?? ?if max_d d? ?? ?? ?? ?? ?max_d = d;? ?? ?? ?? ?end? ?? ?? ?? ?if min_d d? ?? ?? ?? ?? ?min_d = d;? ?? ?? ?? ?end? ?? ?? ?? ?avg_dd = avg_dd + d;? ?? ???end? ? end? ? avg_d = 2*avg_dd/(M*(M-1));? ?? ?? ?? ?? ? %平均相点距离? ? ? ? dlt_eps = (avg_d - min_d) * 0.02 ;? ?? ?? ?%若在min_eps～max_eps中找不到演化相点时，对max_eps的放宽幅度? ? min_eps = min_d + dlt_eps / 2 ;? ?? ?? ?? ?%演化相点与当前相点距离的最小限?