genAiguille.m
61f82e71adb9master

genAiguille.m
View Options

			beep off
			%% para
			Omega = 1.118;
			nu = 0.008;
			B1 = 3;

			t0 = 0;
			nb_steps = 40000;
			% % %small oscillations
			% t1 = 50;
			% % %chaos
			nb_steps_Pmap = 50;
			dt = 2pi/(Omeganb_steps_Pmap);
			t1 = nb_steps*dt;
			q0 = 0;
			eps = 10^(-8);
			q0_2 = q0 + eps;
			qDot0 = 10^(-6);
			%% integrator
			a = @(t,q,qDot)aiguille(t,q,qDot,nu,B1,Omega);
			[t, q, qDot ] = SV(nb_steps,t0,t1,q0 ,qDot0,a);
			[t2,q2,qDot2] = SV(nb_steps,t0,t1,q0_2,qDot0,a);
			%% write in file
			% fileName = 'input/pendulum.csv';
			% generatedData = [t, q];
			% csvwrite(fileName,generatedData)
			%% plot
			figure(1)
			plot(t,q,'b',t2,q2,'r')
			figure(2)
			semilogy(t,abs(q-q2),'.')
			figure(3)
			qP = q - 2pifloor( (q+pi)/(2*pi) );
			qDotP = qDot - 2pifloor( (qDot+pi)/(2*pi) );
			qP = qP(3*nb_steps_Pmap+1:nb_steps_Pmap:end);
			qDotP = qDotP(3*nb_steps_Pmap+1:nb_steps_Pmap:end);
			t_P = t(3*nb_steps_Pmap+1:nb_steps_Pmap:end);
			check = sin(Omega*t_P);%should be 0
			plot(qP,qDotP,'.')