rizwana28282
Newbie level 1
- Joined
- Jan 5, 2012
- Messages
- 0
- Helped
- 0
- Reputation
- 0
- Reaction score
- 0
- Trophy points
- 1,280
- Location
- chennai, India
- Activity points
- 1,280
i m in need to plot the lyapunov exponent 'f vs lyapunov,' but i am unable to run the program due to some error in the program. can anyone pls help me to get the plot and any suggestions to make it perfect.
with kind regards
matlab codes:
clc
global f
for f=0.0:0.0001:0.5
f
initialx(1)=0.0100;initialx(2)=0.01;
lyap(1)=.01;
[t,x]=ode45('loreeneq',0:1:300,[initialx(1),initialx(2)]);
lyap1=lyap([t,x,1)]);
plot(f,lyap1,'r','MarkerSize',0.01);
hold on
drawnow
end
xlabel ('f');ylabel('x')
--------------------------------------------------------------------------
function dxdt=loreeneq(t,x)
global r
dxdt1=10*(x(2)-x(1));
dxdt2=(r*x(1))-x(2)-(x(1)*x(3));
dxdt3=(x(1)*x(2))-(2.67*x(3));
dxdt=[dxdt1;dxdt2;dxdt3];
with kind regards
matlab codes:
clc
global f
for f=0.0:0.0001:0.5
f
initialx(1)=0.0100;initialx(2)=0.01;
lyap(1)=.01;
[t,x]=ode45('loreeneq',0:1:300,[initialx(1),initialx(2)]);
lyap1=lyap([t,x,1)]);
plot(f,lyap1,'r','MarkerSize',0.01);
hold on
drawnow
end
xlabel ('f');ylabel('x')
--------------------------------------------------------------------------
function dxdt=loreeneq(t,x)
global r
dxdt1=10*(x(2)-x(1));
dxdt2=(r*x(1))-x(2)-(x(1)*x(3));
dxdt3=(x(1)*x(2))-(2.67*x(3));
dxdt=[dxdt1;dxdt2;dxdt3];