Sinusoidal model

From Wikipedia, the free encyclopedia

This article does not cite any references or sources. (February 2008)
Please help improve this article by adding citations to reliable sources. Unverifiable material may be challenged and removed.

This article is orphaned as few or no other articles link to it.
Please help introduce links in articles on related topics. (February 2008)

In statistics, signal processing, and time series analysis, a sinusoidal model to approximate a sequence Y_i is:

Y i = C + αsin(ω T i + φ) + E i

where C is constant defining a mean level, α is an amplitude for the sine wave, ω is the frequency, T_i is a time variable, φ is the phase, and E_i is the error sequence in approximating the sequence Y_i by the model. This sinusoidal model can be fit using nonlinear least squares; to obtain a good fit, nonlinear least squares routines may require good starting values for the constant, the amplitude, and the frequency.

Fitting a model with a single sinusoid is a special case of least-squares spectral analysis.

1 Good starting value for C
2 Good starting value for frequency
3 Good starting values for amplitude
4 Model validation
5 External links
6 References

[edit] Good starting value for C

A good starting value for C can be obtained by calculating the mean of the data. If the data show a trend, i.e., the assumption of constant location is violated, one can replace C with a linear or quadratic least squares fit. That is, the model becomes

Y i = (B 0 + B 1 T i) + αsin(2πω T i + φ) + E i

$Y_i = (B_0 + B_1T_i+B_2T_i^2) + \alpha\sin(2\pi\omega T_i + \phi) + E_i$

[edit] Good starting value for frequency

The starting value for the frequency can be obtained from the dominant frequency in a periodogram. A complex demodulation phase plot can be used to refine this initial estimate for the frequency.^{[citation needed]}

[edit] Good starting values for amplitude

A complex demodulation amplitude plot can be used to find a good starting value for the amplitude. In addition, this plot can indicate whether or not the amplitude is constant over the entire range of the data or if it varies. If the plot is essentially flat, i.e., zero slope, then it is reasonable to assume a constant amplitude in the non-linear model. However, if the slope varies over the range of the plot, one may need to adjust the model to be:

Y i = C + (B 0 + B 1 T i)sin(2πω T i + φ) + E i

That is, one may replace α with a function of time. A linear fit is specified in the model above, but this can be replaced with a more elaborate function if needed.

[edit] Model validation

As with any statistical model, the fit should be subjected to graphical and quantitative techniques of model validation. For example, a run sequence plot to check for significant shifts in location, scale, start-up effects, and outliers. A lag plot can be used to verify the residuals are independent. The outliers also appear in the lag plot, and a histogram and normal probability plot to check for skewness or other non-normality in the residuals.