how I create improved lorenz code
6 views (last 30 days)
Show older comments
hi everyone
I want to create improved lorenz code , but its doesnt work well
could someone help me , please
this is my code
function dy= lorenz(t,y)
a = 8/3;
b = 25;
c = 10;
d = 5.3;
e = 17.5;
f = 10.5;
g = 13.3;
h = 14;
if cos(5.3*t)>=0
p=1;
else
p=-1;
end
dy=zeros(3,1);
dy(1)=a*y(1)+p*y(2)*y(3);
dy(2)=y(3)*b-y(3)*c*cos(d*t) - b*y(2)+y(2)+c*cos(d*t);
dy(3)= -1*p*y(1)*y(2)-(e*cos(d*t)+f)*y(2)+(g-h*y(3)*(cos(d*t)));
end
6 Comments
Jan
on 16 Nov 2021
If I follow the instructions from scratch, I get the same output as you. So why to you assume that the diagram of the top is the correct solution? It looks like the standard Lorenz attractor, but the formula is something different.
Accepted Answer
Jan
on 14 Nov 2021
Edited: Jan
on 15 Nov 2021
The test contains cos(dt), but your code only cos(t).
3 Comments
Jan
on 18 Nov 2021
I still do not see a convincing reason to assume, that there is any problem. I've written some code to solve the answer and get the same output as you (as far as I can see). Why do you think that your code is wrong?
More Answers (0)
See Also
Categories
Find more on Startup and Shutdown in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!