fractional discrete henon map

1 次查看（过去 30 天）
sabrina
sabrina 2025-7-7
评论： sabrina 2025-7-29
I’m currently working on the fractional Hénon map, but I’m facing difficulties in drawing the bifurcation diagrams correctly. I’m not sure how to implement the fractional order in the iteration process or which method is most suitable. I would really appreciate any guidance or examples you could provide.
  4 个评论
显示 2更早的评论
John D'Errico
John D'Errico 2025-7-29
@sabrina
Note that you have now asked this exact same question 3 times. I closed the most recent one, as it asks or says nothing more than you have already asked twice before.

请先登录，再进行评论。

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

采纳的回答

Ruchika Parag
Ruchika Parag 2025-7-9
Hi @sabrina, when working with a fractional-order Hénon map, the main challenge is incorporating the fractional dynamics into the iteration process. Unlike the standard map, a fractional map includes memory effects—each state depends not only on the immediate previous value, but also on the full history of the system.
To implement this, you'll need to modify your iteration using a fractional difference method. A common approach is the Grünwald–Letnikov approximation, which sums previous states using fractional binomial weights.
Here's a basic outline to get started:
1. Compute fractional weights: For fractional order q (e.g., 0.9), the Grünwald–Letnikov weights can be computed as:
q = 0.9;
N = 1000;
w = zeros(N,1);
for j = 0:N-1
w(j+1) = (-1)^j * nchoosek(q, j);
end
2. Iterate using history: During each iteration, calculate the weighted sum over past states:
x = zeros(N,1); y = zeros(N,1);
x(1) = x0; y(1) = y0;
for n = 1:N-1
sumx = w(1:n+1)' * x(n:-1:0+1);
sumy = w(1:n+1)' * y(n:-1:0+1);
x(n+1) = 1 - a * sumx.^2 + y(n); % adjust formula as needed
y(n+1) = b * x(n);
end
3. Build the bifurcation diagram: Vary your bifurcation parameter (e.g., a) and plot the final values after discarding transients:
for a = a_min : a_step : a_max
[x, y] = fractionalHenon(a, b, q, N); % assuming you wrap the above loop in a function
plot(a * ones(100,1), x(end-99:end), '.k'); hold on;
end
If you're looking for further references or ready-to-use methods, some of Garrappa’s work on fractional differential equations may help adapt logic for discrete maps as well.
Hope this helps!

更多回答（0 个）

类别

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

标签

产品


版本

R2025a

Community Treasure Hunt

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

Start Hunting!

Translated by