Are you sure you have a problem with your code? I increased N to 16000, and plotted your data. It sure looks like a Hénon attractor to me.
I am not getting unique random numbers from henon map??please help in this
2 views (last 30 days)
Show older comments
VINAY KUMAR SHARMA
on 19 Dec 2017
Commented: VINAY KUMAR SHARMA
on 20 Dec 2017
my code is here
if true
% code
N=16;
a=1.4;
b=0.3;
x=zeros(1,N);
y=zeros(1,N);
x(1)=0.63154*(10^(-11));
y(1)=0.18906*(10^(-11));
for i=1:N
x(i+1)=1+y(i)-a*((x(i))^2);
y(i+1)=b*x(i);
end
plot(x,y)
for j=0:(N/2)-1
q(2*j+1)=x(j+1);
q(2*j+2)=y(j+1);
end
d=floor(q*(10^15));
s=mod(d,N);
end
0 Comments
Accepted Answer
the cyclist
on 19 Dec 2017
Edited: the cyclist
on 19 Dec 2017
4 Comments
the cyclist
on 19 Dec 2017
Can you add some more detail? I don't understand why you expect non-repeating values in s, so it is not possible for me to help you.
Are you trying to produce pseudorandom numbers from the Hénon map? (Numbers if pseudorandom sequences can repeat. Just not in a predictable way, ideally.)
Do you have some reference? Are you trying replicate some results from a paper or a textbook?
It might be that your code is mathematically correct, but there are nuances with doing floating point math.
More Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!