У меня проблема с числовой производной вектора x: Nx1 по отношению к другому вектору t (время) того же размера, что и x. Я делаю следующее (x выбран в качестве функции синуса в качестве примера):

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

но ответ xd равен (N-1) x1, и я понял, что он не вычисляет производную, соответствующую первому элементу x.

Есть ли другой способ вычислить эту производную?