Proper regression for determining correlations between derivatives of functions Say we have a permanent-magnet DC motor that roughly obeys the system equation $\ddot{x}(t) = \alpha \dot{x}(t) + \beta u(t) + \gamma$, where $x(t)$ is the displacement of the rotor, and $u(t)$ the applied voltage, at time $t$.
Say we wish to determine the values of $\alpha$, $\beta$, and $\gamma$ experimentally.  If we can only directly measure $x$, not $\dot{x}$ or $\ddot{x}$, how should we go about estimating these parameters from a set of time-series measurements of $u$ and $x$?
One naive approach is to compute the derivatives through some central finite difference scheme, and then perform an OLS regression - but it is unobvious how the derivative calculation interacts with the regression.  Additionally, I have found in practice that this suffers from a significant amount of regression dilution if the test is allowed to run too long at steady-state (the derivatives vanish here, and so all that's left is the noise).
Is there any more "complete" method for analyzing systems like this that handles the differentiation as part of the construction of the regression model?  Is there a good theory of correlations between derivatives of time-series data?
Update: I figured I'd share some actual plots to compare some of the methods that have been suggested.  In the time-domain plots below, the green line represents the path of a simulated motor that has been stepped along with the characterization constants determined from the regression.  Note that the constants here are presented in slightly different form - the presented parameters match the "feedforward equation" form:
$u(t) = K_a \ddot{x}(t) + K_v \dot{x}(t) + K_s \text{sgn}(\dot{x}(t))$
Strategy 1: Central finite difference followed by direct regression of acceleration against other variables:

Strategy 2: Central finite difference followed by direction regression of next velocity from previous velocity (slightly different UI - this is from an old version of the tool):

Strategy 3: "Accumulated" formulation - regression of the next position against previous position and integrals of other terms:

Alas, for our test procedure this methodology seems more or less useless.
I have yet to successfully run a Fourier-domain analysis; I might give that a shot eventually, but it's substantially more of a chore to implement than the ones I've tried above.

The differential equation can be rewritten
$$x_{t+1}=(\alpha+1)x_t+\beta\int^t u(x)\,\mathrm{d}x+\gamma t+C.$$
This eliminates many of the problematic aspects of the direct formulation and computing numerical derivatives, but it still leaves one untouched: because all the $x_t$ may be subject to measurement error, the explanatory variable $x_t$ on the right hand side might be subject to enough measurement error to bias the estimates.
There is a quick fix: fit the model using the data $(x_t, U_t, t, x_{t+1})$ (where $U(t) = \int_t u(x)\,\mathrm{d}x$) and then refit the model after replacing $x_t$ by the predicted values at those times.
Here, for example, is a dataset of $x_t$ represented by the $50$ dots on the right.

On the right, the derivative $\dot{x}_t$ is shown in blue and the second derivative $\ddot{x}_t$ in red.
Such data were independently generated $500$ times (with unit measurement variance), fit, and re-fit.  To assess these fits, I plotted histograms of the parameter estimates.

In each histogram a vertical red line marks the true parameter value while the vertical dashed black line indicates the mean of the estimated values.  In this example there is appreciable bias in evidence: the red lines are out in the tails of the distributions.
After refitting, the situation has improved, but at substantial cost in uncertainty: the estimates are much more spread out.

If you are limited to code for ordinary least squares solutions, as indicated in comments, you will have to pick your poison: precise biased estimates or much less precise unbiased estimates.  Perhaps an initial study of simulated data will suggest how much bias and imprecision are present in your setting and help you choose.
