Sinusoidal model

In statistics, signal processing, and time series analysis, a sinusoidal model is used to approximate a sequence y_i to a sine function: : y_i = c + \alpha\sin(\omega t_i + \varphi) + \varepsilon_i where c is constant defining a mean level, \alpha is an amplitude for the sine, \omega is the angular frequency, t_i is a time variable, \varphi is the phase-shift, and \varepsilon_i is the error sequence. This sinusoidal model can be fit using nonlinear least squares; to obtain a good fit, routines may require good starting values for the unknown parameters. Fitting a model with a single sinusoid is a special case of spectral density estimation and least-squares spectral analysis.

In statistics, signal processing, and time series analysis, a sinusoidal model is used to approximate a sequence y_i to a sine function: y_i = c + \alpha\sin(\omega t_i + \varphi) + \varepsilon_i where c is constant defining a mean level, \alpha is an amplitude for the sine, \omega is the angular frequency, t_i is a time variable, \varphi is the phase-shift, and \varepsilon_i is the error sequence. This sinusoidal model can be fit using nonlinear least squares; to obtain a good fit, routines may require good starting values for the unknown parameters. Fitting a model with a single sinusoid is a special case of spectral density estimation and least-squares spectral analysis. ==Good starting values== ===Good starting value for the mean=== 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_1t_i) + \alpha\sin(\omega t_i + \varphi) + \varepsilon_i or y_i = (b_0 + b_1t_i+b_2t_i^2) + \alpha\sin(\omega t_i + \varphi) + \varepsilon_i ===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. ===Good starting values for amplitude=== The root mean square of the detrended data can be scaled by the square root of two to obtain an estimate of the sinusoid 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(\omega t_i + \phi) + \varepsilon_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. ==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. ==Extensions== A different method consists in transforming the non-linear regression to a linear regression thanks to a convenient integral equation. Then, there is no need for initial guess and no need for iterative process : the fitting is directly obtained. ==See also== Pitch detection algorithm ==References== ==External links== Beam deflection case study

Text derived from https://en.wikipedia.org/w/index.php?oldid=1343510848 under CC BY-SA 4.0; modifications have been made.