how I create improved lorenz code
4 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
hi everyone
I want to create improved lorenz code , but its doesnt work well
could someone help me , please
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/799809/image.jpeg)
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 commentaires
Jan
le 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.
Réponse acceptée
Jan
le 14 Nov 2021
Modifié(e) : Jan
le 15 Nov 2021
The test contains cos(dt), but your code only cos(t).
3 commentaires
Jan
le 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?
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Startup and Shutdown dans Help Center et File Exchange
Produits
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!