matlab 符号值转为实数,Matlab：在实现二进制到实数转换的公式时,无法获得唯一的有理数Part1...

存在非线性动态系统x_n = f(x_n,eta),其函数形式为x [n 1] = 2 * x [n] mod 1.这是一个混沌动力系统,称为锯齿图或伯努利图.我在实现Eq(4)和Eq(5)给出的逆映射的两个表示时面临困难.以下是该问题的简要说明.

其中序列(s [n k])_ k = 1到N-1是状态x [n]的符号描述.该描述源于下面描述的单位间隔的划分.

设分区数M = 2,符号空间= {0,1},分配符号的规则是

s[n+1] = 1 if x[n] >= 0.5, otherwise s[n+1] = 0

本文作者：

Linear, Random Representations of Chaos

对于方程(5),我在反向之后没有获得相同的时间序列,在进行二进制到实数转换之后几个值不同.有人可以让我正确的程序吗？

我试图为Eqs(4)和(5)实现Bijective图,但它不起作用.

公式(5)的代码 – 我将二进制二进制化. x包含实数; s是每个实数的0/1二进制等价物; y是将s转换为real之后的答案. s1是x的1 / -1二进制等价物; b是转换为真实后的答案.在1 / -1的情况下,当我从符号表示转换为实数时,我将-1切换为0,然后在公式(5)中应用公式.从答案中可以看出,在进行转换后y和b与x不同.当原始的实数都是无符号的理性时,我也得到了负的实数！我怎样才能正确实施以使它们相同？

N =10;

x(1) = 0.1;

for i =1 : N

x(i+1) = mod(x(i)*2, 1);

end

y = x;

s = (y>=0.5); %generate 0/1 logicals

for n = 1: N

y(n) = 0.5*s(n+1) + 0.5*y(n+1);

end

b=x;

s1 = 2*(b>=0.5)-1; %Generate +1/-1

for k =1: N

if s1(k)== -1

s1(k) = 0;

end

b(k) = 0.5*s1(k+1) + 0.5*b(k+1);

end

设,x =

0.100000000000000 0.200000000000000 0.400000000000000 0.800000000000000 0.600000000000000 0.200000000000000 0.400000000000000 0.800000000000001 0.600000000000001 0.200000000000003 0.400000000000006

y =

0.100000000000000 0.200000000000000 0.900000000000000 0.800000000000000 0.100000000000000 0.200000000000000 0.900000000000000 0.800000000000001 0.100000000000001 0.200000000000003 0.400000000000006

b =

-0.400000000000000 0.700000000000000 0.900000000000000 -0.200000000000000 -0.400000000000000 0.700000000000000 0.900000000000000 -0.199999999999999 -0.399999999999999 -0.299999999999997 0.400000000000006