%szinkronizálódó trajekt kirajzolása idõben és káosz a rendszerben (3D, 1D)

% szinkronizációhoz, összekapcsolt rendszerek
kezd1=[-0.6 0 0.3]; % 1. Chua kör kezdõpontja
kezd2=[0.3 0.2 -0.7]; % 2. Chua kör kezdõpontja

%kezdeti értéktõl való érzékeny függés vizsgálatához másik rendszer
kezd3=[-0.61 0 0.3]; % 1. Chua kör kezdõpontja kicsit máshol
kezd4=[0.31 0.2 -0.7]; % 2. Chua kör kezdõpontja kicsit máshol

alfa=11.6;
beta=16.3;
Rg1=1.24;
Rg3=0.04;
gcs=0.7%0.7;%0.5;%1;%0;%
sync_chua_sys=@(t,x) [-alfa*x(1)+alfa*x(2)-alfa*((-Rg1)*x(1)+Rg3*x(1)^3);...
                  x(1)-x(2)+x(3)-gcs*(x(2)-x(5));...
                 -beta*x(2);...
                 -alfa*x(4)+alfa*x(5)-alfa*((-Rg1)*x(4)+Rg3*x(4)^3);...
                  x(4)-x(5)+x(6)+gcs*(x(2)-x(5));...
                 -beta*x(5) ];
             
n=200;
xl=[-1 n];
yl=[-6 6];
figure;
[t, x] = ode45(sync_chua_sys,[0 n], [kezd1 kezd2]);
[t2, x2] = ode45(sync_chua_sys,[0 n], [kezd3 kezd4]); %másik, független rendszer

%for i=0.001:1:n
%[t, x] = ode45(sync_chua_sys,[0 i], [kezd1 kezd2]);
%[t2, x2] = ode45(sync_chua_sys,[0 i], [kezd3 kezd4]); %másik, független rendszer 

subplot(6,1,[1,2,3]);
plot3(x(:,1),x(:,2),x(:,3),'r', x(:,4),x(:,5),x(:,6),'b'); %csatolás
%plot3(x(:,1),x(:,2),x(:,3),'b', x2(:,1),x2(:,2),x2(:,3),'c'); %kezdeti értéktõl való érzékeny függés   
%plot3(x(:,1),x(:,2),x(:,3),'r', x(:,4),x(:,5),x(:,6),'m', x2(:,1),x2(:,2),x2(:,3),'b', x2(:,4),x2(:,5),x2(:,6),'c'); %összes 1-ben
xlabel('x1');
ylabel('y1');
zlabel('z1');
%legend('Chua1 kp1', 'Chua2 kp1', 'Clua1 kp2', 'Chua1 kp2')

subplot(6,1,4);
plot(t,x(:,1),'r', t,x(:,4),'b'); %csatolás
%plot(t,x(:,1),'b', t2,x2(:,1),'c'); %kezdeti értéktõl való érzékeny függés 
%plot(t, x(:,1), 'r', t, x(:,4), 'm', t2, x2(:,1), 'b', t2, x2(:,4), 'c'); %összes 1-ben
xlabel('t');
ylabel('x1');
xlim(xl);

subplot(6,1,5);
plot(t,x(:,2),'r', t,x(:,5),'b'); %csatolás
%plot(t,x(:,2),'b', t2,x2(:,2),'c'); %kezdeti értéktõl való érzékeny függés 
%plot(t, x(:,2), 'r', t, x(:,5), 'm', t2, x2(:,2), 'b', t2, x2(:,5), 'c'); %összes 1-ben
xlabel('t');
ylabel('y1');
xlim(xl);

subplot(6,1,6);
plot(t,x(:,3),'r', t,x(:,6),'b'); %csatolás
%plot(t,x(:,3),'b', t2,x2(:,3),'c'); %kezdeti értéktõl való érzékeny függés 
%plot(t, x(:,3), 'r', t, x(:,6), 'm', t2, x2(:,3), 'b', t2, x2(:,6), 'c'); %összes 1-ben
xlabel('t');
ylabel('z1');
xlim(xl);

%waitforbuttonpress;
pause(0.1);
%end