			beep off
			%% para
			Omega = 2*1.118;
			nb_steps = 100;
			nb_steps_excitation = 50;
			dt = 2pi/(Omeganb_steps_excitation);
			t0 = 0;
			t1 = nb_steps*dt;
			q0 = 0;
			eps = 10^(-8);
			q0_2 = q0 + eps;
			qDot0 = 10^(-6);
			fileName = 'input/pendulum.csv';

			% a = @(t,q)pendulum(t,q);
			% a = @(t,q,qDot)pendulumForced(t,q,qDot,0.01,12,0.5);
			a = @(t,q,qDot)aiguille(t,q,qDot,0.08,3,Omega);
			% [t,q]=eulerProg(nb_steps,t0,t1,q0,qDot0,a);
			[t,q,qDot]=RK2(nb_steps,t0,t1,q0,qDot0,a);
			% [t2,q2,qDot2]=RK2(nb_steps,t0,t1,q0_2,qDot0,a);
			[t2,q2,qDot2]=SV(nb_steps,t0,t1,q0_2,qDot0,a);
			generatedData = [t, q];
			% csvwrite(fileName,generatedData)
			%% plot
			figure(1)
			plot(t2,q2)
			% plot(t,q,t2,q2)
			% figure(2)
			% plot(t,abs(q-q2))
			figure(3)
			plot(q,qDot)

R2096/generate/gen.m
c755270c71b7master

gen.m
View Options