how I create improved lorenz code

1 次查看（过去 30 天）
Ghofran Khaled
Ghofran Khaled 2021-11-14
评论： Jan 2021-11-18
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 个评论
显示 4更早的评论
Jan
Jan 2021-11-16
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.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Jan
Jan 2021-11-14
编辑：Jan 2021-11-15
The test contains cos(dt), but your code only cos(t).
  3 个评论
显示 1更早的评论
Ghofran Khaled
Ghofran Khaled 2021-11-15
@Jan Please can you help me this took me a week and I don't know where I went wrong I'm out of time
Jan
Jan 2021-11-18
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?

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Startup and Shutdown 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by