%trajekt kirajzolása idõben (3D, 1D)

s=10; %paraméterek
r=28;
b=8/3;
Lorenz_egyenlet=@(t, xyz) [s*(xyz(2)-xyz(1)); r*xyz(1)-xyz(2)-xyz(1)*xyz(3); xyz(1)*xyz(2)-b*xyz(3)]; % Lorenz egyenletek

kezd1=[5 5 5]; % kezdõpont 1
kezd2=[4.99 5 5]; % kezdõpont 2

tmax=50; %eddig futtat
tlim=[0 tmax]; %határok ábrázoláshoz
xl=[-20,20];
yl=[-20 20];
zl=[-20,50];

figure;
%for i=0.001:0.1:tmax %lépésenkénti folyamatos kirajzolás
    %[t, x] = ode45(Lorenz_egyenlet,[0 i], kezd1); %mindig az adott idõig számolok
    
    [t, x] = ode45(Lorenz_egyenlet,[0 tmax], kezd1); %diffegy megold
 %   [t2, x2] = ode45(Lorenz_egyenlet,[0 tmax], kezd2); %2. pont ha kell
    
%x,y,z az idõ függvényében ábrázolva
%subplot(6,1,1);
subplot(6,2,1);
plot(t,x(:,1));%, t2,x2(:,1)); %2. pontot is kirajzol
xlabel('t');
ylabel('x');
xlim(tlim);
ylim(xl);

%subplot(6,1,2);
subplot(6,2,3);
plot(t,x(:,2));%, t2,x2(:,2));
xlabel('t');
ylabel('y');
xlim(tlim);
ylim(yl);

%subplot(6,1,3);
subplot(6,2,5);
plot(t,x(:,3));%, t2,x2(:,3));
xlabel('t');
ylabel('z');
xlim(tlim);
ylim(zl);

% 3D-s ábra
%subplot(6,1,[4,5,6]);
subplot(6,2,7:12);
plot3(x(:,1),x(:,2),x(:,3));%, x2(:,1),x2(:,2),x2(:,3)) %3D plotolás
xlabel('x');
ylabel('y');
zlabel('z');
xlim(xl);
ylim( yl);
zlim(zl);
%
%
subplot(6,2,[2,4,6]);
%Lorenz peak map megrajzolása
[mi, hol]=findpeaks(x(:,3)); %beépített függvény a signal processing toolboxból lokális maximumokat keres meg egy függvényen
%mi=x((find(diff(sign(diff(x(:,3))))==-2)+1),3); %"házilag" kiszámolva

peak_n1=mi(1:end-1); %ábrázoláshoz 2 adatsor ezeket ábrázolom egymás függvényében
peak_n2=mi(2:end);
%[mi2, hol]=findpeaks(x2(:,3));
%peak2_n1=mi2(1:end-1);
%peak2_n2=mi2(2:end);
plot(peak_n1, peak_n2, '.');%, peak2_n1, peak2_n2, '.'); %a lokális maximumok ábrája: az n+1.-et az n. függvényében
xlim([30 50]);
ylim([30 50]);
xlabel('Zn');
ylabel('Zn+1');
title('Lorenz peak map');
%}
%waitforbuttonpress; %gomb nyomásra rajzol
%pause(0.1); %késleltet rajzolást
%end %lépésenkénti folyamatos kirajzolás