Tengo un problema con la derivada numérica de un vector que es x: Nx1 con respecto a otro vector t (tiempo) que tiene el mismo tamaño de x. Hago lo siguiente (se elige x como función seno como ejemplo):

t=t0:ts:tf;
x=sin(t);
xd=diff(x)/ts;

pero la respuesta xd es (N-1) x1 y descubrí que no calcula la derivada correspondiente al primer elemento de x.

¿Hay alguna otra manera de calcular esta derivada?