?? function lyapunov(x).txt
字號:
function Lyapunov(x)
if length(x)<400
helpdlg('數據量需要大于400,請重新選定時間','提示信息');pause(3),
end
dm=26;tao=5;set_error=0.001;
for i=0:(dm-1)
for j=1:(length(x)/1.2)
k=j+i*tao;
Y(i+1,j)=x(k);
end
end
l=1; i=1; p=1; % i為Y的一列,j為Y的一列
while i<length(Y)
for j=p:length(Y) %求兩點之間距離
if j>i
e_normal=max(abs(Y(:,i)-Y(:,j)));
end
end
flag=e_normal-set_error;
if flag>0
d(l)=e_normal;l=l+1;
%i=i+1;
for j=1:length(Y) %找離點最小距離的點
l_normal(j)=max(abs(Y(:,i)-Y(:,j)));
end
p=find((l_normal-min(l_normal))==0); %找到最近點Y(:,j)
else
j=j+1;i=i+1; %范數小于set_error,Y(:,i)和Y(:,j)繼續向前一步
end
end
LE1=0;
for i=1:l
f1=d(i);f2=d(1);
LE1=log10(f1/f2)+LE1;
end
LE1=1/(l*0.014)*LE1,
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -