how I create improved lorenz code
4 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
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 Commenti
Jan
il 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.
Risposta accettata
Jan
il 14 Nov 2021
Modificato: Jan
il 15 Nov 2021
The test contains cos(dt), but your code only cos(t).
3 Commenti
Jan
il 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?
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Startup and Shutdown in Help Center e File Exchange
Prodotti
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!