混沌映射特性.zip_matlab 帐篷_tent映射_tent混沌_分岔 tent_混沌映射

%********************************************** %新复合序列值的分布 %********************************************** clear all;clc; x(1)=1; %a=2^7; %a=2^15; a=2^31;%可以选择不同长度a N=9999; x(1)=23; y(1)=2; x1(1)=36236; x2(1)=3690; Tent(1)=1; linear(1)=36529; linear1(1)=365; linear2(1)=3652; d=bitxor(uint8(82),uint8(160)); for i=1:N linear(i+1)=mod(linear(i)*16807,(2^31-1)); linear1(i+1)=mod(linear1(i)*16807,(2^31-1)); x(i)=linear1(i+1); if x(i)==0 x(i+1)=floor(4*x(i)^3/(a^2)-12*x(i)^2/a+9*x(i))+1; elseif x(i)==a x(i+1)=floor(4*x(i)^3/(a^2)-12*x(i)^2/a+9*x(i))+1; elseif x(i)==a/2 x(i+1)=floor(4*x(i)^3/(a^2)-12*x(i)^2/a+9*x(i))-1; else x(i+1)=floor(4*x(i)^3/(a^2)-12*x(i)^2/a+9*x(i)); end if Tent(i)<a Tent(i+1)=floor(2*Tent(i)+1); else Tent(i+1)=floor(2*(a*2-1-Tent(i))); end %y(i)=x(i+1); y(i+1)=floor(4*y(i)-2*y(i)^2/a-1); z(i)=bitxor(uint32(x(i+1)),uint32(y(i+1))); z(i)=bitxor(uint32(z(i)),uint32(linear(i+1))); end % m=mean(z); % for i=1:N % if z(i)<m % z(i)=-1; % else % z(i)=1; % end % end % z=uint16(z); % figure(1); % [a,b]=xcorr(z,'coeff'); % plot(b,a);title('自相关函数'); figure(2); plot(z(1,:),'k.');%用来画出序列的点的分布