I have a signal sampled at $\Delta t$: $f_i(t_i=i\Delta t)$ where $i = 0,\ldots,n-1$. I want to find the first and second derivative of the signal: $f'(t)$ and $f''(t)$.
My first thought was to estimate the derivatives by central differences:
\begin{align} f'(t_{i})&=\frac{f(t_{i+1})-f(t_{i-1})}{2\Delta t}\\ f''(t_{i})&=\frac{f(t_{i+1})-2f(t_{i})+f(t_{i-1})}{(\Delta t)^2} \end{align}
However the signal may have a lot of high frequency noise that may cause quick fluctuations in $f'$ and $f''$.
What would be the best way to find "smoothed" estimates of $f'$ and $f''$?