Sinusoidal model

Sinusoidal model

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

:Y_i = C + alphasin(omega T_i + phi) + E_i

where "C" is constant defining a mean level, α is an amplitude for the sine wave, ω is the frequency, "Ti" is a time variable, φ is the phase, and "Ei" is the error sequence in approximating the sequence "Yi" 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.

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_1T_i) + alphasin(2piomega T_i + phi) + E_i

or

:Y_i = (B_0 + B_1T_i+B_2T_i^2) + alphasin(2piomega T_i + phi) + E_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.Fact|date=February 2008

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(2piomega T_i + phi) + 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.

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.

External links

* [http://www.itl.nist.gov/div898/handbook/eda/section4/eda425.htm Beam deflection case study]

References


Wikimedia Foundation. 2010.

Игры ⚽ Поможем написать курсовую

Look at other dictionaries:

  • Dbx Model 700 Digital Audio Processor — The dbx Model 700 Digital Audio Processor was a professional audio ADC/DAC combination unit, which digitized a stereo analog audio input into a bitstream, which was then encoded and encapsulated in an analog composite video signal, for recording… …   Wikipedia

  • dbx Model 700 Digital Audio Processor — The dbx Model 700 Digital Audio Processor was a professional audio ADC/DAC combination unit, which digitized a stereo analog audio input into a bitstream, which was then encoded and encapsulated in an analog composite video signal, for recording… …   Wikipedia

  • Fermi–Ulam model — The Fermi ndash;Ulam model (FUM) is a dynamical system that was introduced by Polish mathematician Stanislaw Ulam in 1961.FUM is a variant of Enrico Fermi s primary work on acceleration of cosmic rays, namely Fermi acceleration. The system… …   Wikipedia

  • List of mathematics articles (S) — NOTOC S S duality S matrix S plane S transform S unit S.O.S. Mathematics SA subgroup Saccheri quadrilateral Sacks spiral Sacred geometry Saddle node bifurcation Saddle point Saddle surface Sadleirian Professor of Pure Mathematics Safe prime Safe… …   Wikipedia

  • Audio timescale-pitch modification — Time stretching is the process of changing the speed or duration of an audio signal without affecting its pitch.Pitch scaling or pitch shifting is the reverse: the process of changing the pitch without affecting the speed. There are also more… …   Wikipedia

  • List of statistics topics — Please add any Wikipedia articles related to statistics that are not already on this list.The Related changes link in the margin of this page (below search) leads to a list of the most recent changes to the articles listed below. To see the most… …   Wikipedia

  • Sine wave — Sinusoid redirects here. For the blood vessel, see Sinusoid (blood vessel). The graphs of the sine and cosine functions are sinusoids of different phases. The sine wave or sinusoid is a mathematical function that describes a smooth repetitive… …   Wikipedia

  • Thaw (weather) — Thaw is the period when the snow and ice melt, at the end of the winter, in cold climates.A special thaw, outside normal season, is January thaw. January thaw The January thaw is an observed but unexplained temperature rise in mid winter found in …   Wikipedia

  • Sine — For other uses, see Sine (disambiguation). Sine Basic features Parity odd Domain ( ∞,∞) Codomain [ 1,1] P …   Wikipedia

  • Map projection — A medieval depiction of the Ecumene (1482, Johannes Schnitzer, engraver), constructed after the coordinates in Ptolemy s Geography and using his second map projection A map projection is any method of representing the surface of a sphere or other …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”