newchao.m

来自「Chaotic maps. Matlab simulation.」· M 代码 · 共 49 行

% OPTICAL LOGISTIC MAP BUILT USING MACH ZEHNDER INTERFEROMETERS.SO BEST IS TO GO FOR THIS TYPE OF COMMUNICATION USING CHAOTIC SEQUECNES.THIS PROGRAM
% IS BASED ON K.UMENOS PAPER PRESENTED AT 5th EXPERIMENTAL CHAOS CONFERENCE REGARDING OPTICAL IMPLEMENTATION OF CHAOTIC MAPS USING MACH ZEHNDER
% INTERFEROMETERS WHICH ARE N Numbers.THIS HAS BEEN SUCCESSFULLY DEMONSTRATED FOR PRODUCIG NECESSARY CHIP SEQUENCE FOR OPTICAL CDMA 

clc;  close all; clear all;
theta = 45;  k = 3;  X(1) = sin(theta*pi*(k^0))^2;
spc(1) = X(1);

for ii = 2:1:6000  
    spc(ii) = spc(ii-1)*(3 - (4*spc(ii-1)))^2;
    X(ii)   = (4*X(ii-1)*(1-X(ii-1)));
end

for it = 1:1:length(spc)-2   
   GSe(it) =  spc(it);
   GFi(it) =  spc(it+1);
   GTh(it) =  spc(it+2);
   Se(it)  =  X(it);
   Fi(it)  =  X(it+1);
   Th(it) =   X(it+2);   
end

is = 1;
for ss = 1:1:length(X)        
    if mod(ss,2) == 0    
    MNR(ss)  = spc(is);
    is = is + 1;
    else        
    MNR(ss)  = X(is);  
    is = is + 1;
    end        
end

for ll = 1:1:length(MNR)-2
   MFi(ll)  =  MNR(ll);
   MSe(ll)  =  MNR(ll+1);
   MTh(ll) =   MNR(ll+2);   
end

figure(1);   plot(spc);               title('\bf Cubic MAP');          figure(2);   plot(GSe,GFi,'r.');       title('\bf F R Map');
figure(3);   plot(X);                 title('\bf Logistic MAP');       figure(4);   plot(Se,Fi,'r.');         title('\bf F R Map');
figure(5);   plot3(Se,Fi,Th,'k.');    title('\bf CUBIC');    grid on;  figure(6);   plot3(GSe,GFi,GTh,'k.');  title('\bf LOGISTIC');  grid on;

figure(7);   plot3(MFi,MSe,MTh,'r.'); title('\bf MIXED MAP');  xlabel('\bf X(n)');    ylabel('\bf X(n+1)');  zlabel('\bf X(n+2)');
figure(8);   plot(MFi,MSe,'r.');      title('\bf MIXED MAP');  xlabel('\bf X(n)');    ylabel('\bf X(n+1)');
figure(9);   plot(MSe,MTh,'r.');      title('\bf MIXED MAP');  xlabel('\bf X(n)');    ylabel('\bf X(n+1)');
figure(10);  plot(MFi,MTh,'r.');      title('\bf MIXED MAP');  xlabel('\bf X(n+1)');  ylabel('\bf X(n+2)');