How to create a function to get bifurcation plot
6 次查看(过去 30 天)
显示 更早的评论
Study period-doubling in the Lorenz model by examining the behavior for r≤ 160. Claculate the bifurcation diagram and extract the value of Feigenbaum's δ parameter. Hint:While this problem can be done using the Euler method, it is probably advisable, in order to conserve computer time, to use the Runge-Kutta algorithm. Your matlab function should produce one clearly labeled bifurcation line-plot for period doubling->chaos transition in z,as r changes from 25->160.
Hi guys, I'm a totally new one to learn how to use matlab and met this problem. My personal thoughts about this question is firstly create:function [x, y, z, t] = problem(x0y0z0, tmax, dt, r) And my inputs should be x0y0z0 a 1x3 vector containing the initial values for x0, y0 and z0
tmax: the max time to iterate
dt: the step size
r: the range of r values to examine (I want to try the range of 100-160)
and my outputs can be Outputs: x, y, z as functions of time t, all vectors of length ~tmax/dt
But I truly don't know how to achieve all the details about codes and I m asking help here. Hoping you guys can show your hands and thanks for that so much!!!!
0 个评论
回答(2 个)
Ankitha Kollegal Arjun
2017-4-17
The following MATLAB Answers post might be a useful starting point for your question:
0 个评论
Yulia Verbitskaya
2019-4-23
m=[0:0.001:2]
A=[0:0.02:3]
for i = 1:length(m)
p(i)=(k3+(k_3*A))./K4*m(i);
MPF_MPF(i)=k1-(k2+k_2*APCon(i)*APCon);
MPF_APCon(i)=(((k3+k_3*A)/k4*m1)*(1-APCon(i))*(K4+APCon(i))/(APCon(i))*(K3+1-APCon(i)));
[cross,index]=min(abs(MPF_MPF-MPF_APCon));
MPF_APCon(i)=MPF_MPF(index);
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Communications Toolbox 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!