Starting Values considering Ieee 1057
A Procedure for Starting Values for Sine Fitting (IEEE1057)Stephen Neil Summary: The four type sine fitting procedure IEEE1057 provides an accurate and fast method for fitting a sine lay out to noisy data. In any case, the lines requires certainly starting value anent low frequency (the incomparable parameters from a 3 parameter sine fit) to guard convergence. The procedure described here provides a simple robust course being obtaining starting values. The procedure is easily described passageway simple stairs from a starting data series. 0. remove the mean value from the unexpended data 1. calculate the cumulative continuum excepting the conceptive data 2. take account of the average of the firsthand series 3. estimate the additional humble it's average. Call this series 'the integral' 4. quantize the first difference series of the inspired dataseries. Call this series 'the derivative' 5. compute the foul absolute deviation of the monistic: MAD1 6. calculate the mean transcendent whimsy of the derivative: MAD2 7. compute square root as regards MAD2\MAD1: our estimate of the frequency<\p>
An elementary knowledge of calculus only is imperative towards see how replacing A.Sin(w*t+p) we are estimating the integral, -(A\w).Cos(w*t+p), and the development, A.Cos(w*t+p). Problems of form and sign are avoided by in virtue of the mean absolute deviation of these series, their ratio being w^2 Some simple tests by the author apprehend revealed this procedure is rugged assumed at least 2 cycles, frequency between.01 and.5 and noise of amplitude widen to half the sine amplitude. The contrivance is thus and thus simple it can be there readily tested in a simple spreadsheet Another behavioral norm more robust to sonorous report but only suitable for frequencies infra 0.1 is the 'twice integral' method The procedure is besides easily described in simple steps from a starting data pursuance. 0. cast aside the mean substance from the original series 5. meter the fairish absolute obliqueness of the original series: MAD1 1. count on the cumulative series without the original series 2. calculate the equatorial as to the cumulative series 3. sound the cumulative less it's average. Call this series 'the integral' 1. build on the cumulative series from the integral coming after 2. rationalize the average of the cumulative dogging 3. calculate the certain less it's average. Preconization this succession 'twice integral' 6. calculate the mean absolute transition in point of the twice monadic: MAD2 7. calculate square smell around with regard to MAD2\MAD1: our estimate of the frequency<\p>
For A.Sin(w*t+p) we are estimating the twice integral, -(A\w^2).Sin(w*t+p) <\p>










