clear;clc;
%digital resolution is 100000 
resolution=100000;

i=1:resolution;
j=(2*pi)/resolution;

%The first function
amplitude=7;
frequency=3;

y1=amplitude*cos(frequency*i*j);

%The second function
amplitude=10;
frequency=11;

y2=amplitude*sin(frequency*i*j);

y=y1+y2;

plot(y)
FFTy=fft(y)';

i=1:21;
FFTy(i);

plot(y)
FFTy=fft(y)';
figure(1)
dy1=diff(y)*resolution;
plot(dy1)
ylabel('velocity')
figure(2)
dy2=diff(dy1)*resolution;
plot(dy2)
ylabel('acceleration')