frey_value.m

frey map混沌随机序列的MATLAB产生源程序-frey map chaotic random sequences generated source MATLAB

%Elapsed time is 104367.048920 seconds.
close all;
clear all;
clc
tic
% N est la pr閏ision
N=16;
nnn=100000;
ku=0;

%eu(1)=22000;
eu(1)=20000;

eu(2)=400;


eun(1)=floor(eu(1)*(2^N))/2^N;
eu01(1)=eun(1)/(2^N);

es(1)=nonsigne_signe(eun(1),N);
esn(1)=es(1)/(2^(N-1));

yb(1)=floor(eu01(1).*(2^N));

eun(2)=floor(eu(2)*(2^N))/(2^N);
eu01(2)=eun(2)/(2^N);

es(2)=nonsigne_signe(eun(2),N);
esn(2)=es(2)/(2^(N-1));

yb(2)=floor(eu01(2).*(2^N));

for i = 3 :nnn
    if ( eun(i-2) < 2^(N-1) )
        su(i)=0;
    else
        su(i)=1;
    end
    
  eun(i)=ku+eun(i-1)+2*eun(i-2)+su(i);
  eun(i)=mod(eun(i),2^N);
  
eun(i)=floor(eun(i)*(2^N))/(2^N);
eu01(i)=eun(i)/(2^N);

es(i)=nonsigne_signe(eun(i),N);
esn(i)=es(i)/(2^(N-1));
end

yb=floor(eu01.*(2^N));
yb(nnn)
eun(nnn)

% Lecture et Ecriture d'un fichier binair
fid=fopen('frey3125000.bin','wb');
fwrite(fid,yb,'uint32');
fclose(fid);

fid=fopen('frey3125000.bin','rb');
yr=fread(fid,'uint32');
yr=yr';
fclose(fid);

% Lecture et Ecriture d'un fichier binair
y=dec2bin(uint32(yb));
[M1,N1]=size(y); 
y=y';
s1=M1*N1;
li1=reshape(y,1,s1);
fid=fopen('frey243200.txt','w');
fwrite(fid,li1);
fclose(fid);


yr=yr/(2^N);
save 'frey3125000.mat';
toc

% figure;
% plot(eu01);title('frey');
% 
% figure;
% plot(yr); title('frey ? partir de signal binaire');
% 
% xs=yr(1:(nnn-1));
% x1s=yr(2:nnn);
% 
% figure;
% stairs(xs,x1s);title('attracteur ? partir de fichier binaire');
% 
% xs1=eu01(1:(nnn-1));
% x1s1=eu01(2:nnn);
% 
% figure;
% stairs(xs1,x1s1);title('attracteur ');
% 
% figure;
% z2=corr_a_s(esn);
% title('autorcorrelation');
% 
% figure;
% z2=intercorr_a_s(esn);
% title('intercorrelation');
% 
% figure;
% yf2=fft(esn,2048);
% yf2=abs(yf2);
% plot(yf2);
% title('spectre');