logistic.m

通过混沌理论系统产生混沌序列

N=100;                            %生成序列的长度
X=double(zeros(1,N));       %定义初始序列
u=4;                                %混沌系统初始条件
x0=0.9;
for(i=1:N)                         %产生混沌序列
    X(1,i)=u*x0*(1-x0);
    x0=X(1,i);
end
figure(1);
i=1:1:10;                           %绘图
plot(i,X(1,i));
%产生二值序列方法一
num0=0;
num1=0;
for (i=1:N)
    if X(1,i)<=0.5
        b(1,i)=0
        num0=num0+1;
    else 
       b(1,i)=1
       num1=num1+1;

   end
end
figure(2)
i=30:1:60
plot(i,b(1,i));
title('method 1');
num0;
num1;
%产生二值序列方法2
mean=0;
for i=1:N
    mean=mean+X(1,i);
    mean=mean/N;
    if X(1,i)<=mean;
        b2(1,i)=0
    else 
        b2(1,i)=1;
    end
end
figure(3);
i=30:1:60
plot(i,b2(1,i));
title('method 2');
%产生序列方法3
p=50;
for i=1:p
    if  X(1,i)<=0.5
        b3(1,i)=0
       else 
       b3(1,i)=1
    end
end
 for i=p+1:N
     if X(1,i)<=X(1,i-p)
         b3(1,i)=0;
     else
         b3(1,i)=1;
     end
 end
figure(4);
i=30:1:60
plot(i,b3(1,i));
title('method 3');