你当然可以在实数空间中计算它,但是你有可能遇到精度问题(取决于起点).如果您对研究轨道感兴趣,您可能更愿意使用合理的分数表示法.有更有效的方法可以做到这一点,但下面的代码说明了一种计算从该映射派生的系列的方法.您将在Link 2的第2页上看到period-n定义.您应该能够从此代码中看到如何在实数空间中轻松地工作作为替代(在这种情况下,matlab函数鼠将恢复从你的真实数字的理性近似).
[编辑]现在二进制序列明确!
% start at some point on period-n orbit
period = 6;
num = 3;
den = 2^period-1;
% compute for this many steps of the sequence
num_steps = 20;
% for each step
for n = 1:num_steps
% * 2
num = num * 2;
% mod 1
if num >= den
num = num - den;
end
% simplify rational fraction
g = gcd(num, den);
if g > 1
num = num / g;
den = den / g;
end
% recover 8-bit binary representation
bits = 8;
q = 2^bits;
x = num / den * q;
b = dec2bin(x, bits);
% dis