next previous
Up: Method of running

5. Fits of the harmonic signals

If the input signal is sinusoidal:
equation976
then the corresponding smoothing curve is
eqnarray978
where tex2html_wrap_inline2546 and
eqnarray985
Asymmetric fits change not only the amplitude, but the phase as well. Such phase distortions occur at the temporal edges of the observations. For intermediate values of the argument, z1=-1 and z2=1. Because of the symmetry of the function tex2html_wrap_inline2552, in this case tex2html_wrap_inline2554 for all tex2html_wrap_inline2556.

  figure997
Figure 4: Dependence of the amplitude of the harmonic fit smoothing a model sine function of unit amplitude and a period P on tex2html_wrap_inline2560

For our test symmetric fits,
eqnarray1002
These functions are shown in Fig. 4 (click here). A function tex2html_wrap_inline2562 crosses zero at the value tex2html_wrap_inline2564, much larger as compared with tex2html_wrap_inline2566 corresponding to a zero of tex2html_wrap_inline2568. Some useful values of tex2html_wrap_inline2560 are 0.7047 and 0.8553, which correspond to tex2html_wrap_inline2572 and 1/2, respectively.

The power spectra S(x(t),f) of the smoothed signal are proportional to r2 (e.g. Box & Jenkins 1970), thus for harmonic signal with frequency f=1/P one may obtain tex2html_wrap_inline2582. However, for the deviations from the fit
eqnarray1023
At high frequencies tex2html_wrap_inline2584, thus fast variations are not affected by removal of a "slow trend" such as the slope tex2html_wrap_inline2586 of the power spectra tex2html_wrap_inline2588 detected in some cataclysmic variables (Andronov 1993). Applications corresponding to the integer values tex2html_wrap_inline2590 are described by Terebizh (1992). Influence of the trend removal on the shape of the autocorrelation function (ACF) was studied by Andronov (1994).

The functions tex2html_wrap_inline2592 cross zero for all fits with finite z1 and z2. The "ideal" "rectangular" shape H=1, if tex2html_wrap_inline2600 corresponds to an "ideal low frequency" signal tex2html_wrap_inline2602 usually described in radiotechnics (e.g. Baskakov 1983). Here tex2html_wrap_inline2604 tex2html_wrap_inline2606 tex2html_wrap_inline2608. For real functions tex2html_wrap_inline2610 limiting values of z are finite, causing for some frequencies the negative values of tex2html_wrap_inline2592. In this case, the "output" smoothing function is in anti-phase with the "input" signal, thus the difference tex2html_wrap_inline2616 will have even larger amplitude than the "true" one. The maximum "amplification" occurs for 4 tested methods at tex2html_wrap_inline2618 0.715, 1.112, 1.112, 1.489 with corresponding values of tex2html_wrap_inline2620 1.4816, 1.0839, 1.5421, 1.1446. For the "white noise", mathematical expectation of the amplitude does not depend on frequency, thus a power spectrum of the detrended observations may show an apparent peak due to negative values of tex2html_wrap_inline2622 For symmetric approximations, values of tex2html_wrap_inline2624 are much smaller for weighted fits, than for unweighted. Case "wp" is not significantly worse than "wm" according to this criterion, but has a strong advantage due to wider shape of tex2html_wrap_inline2592.

For noisy sinusoidal signal tex2html_wrap_inline2628 and a "continuous" approximation one may estimate a "signal/noise" ratio
equation1037
where tex2html_wrap_inline2630, r0 is the amplitude of a harmonic component and tex2html_wrap_inline2634 is an accuracy estimate of an individual values which may be set to tex2html_wrap_inline2192. For fixed n (number of observations per period P), this ratio increases with increasing tex2html_wrap_inline1958 proportionally to tex2html_wrap_inline2644 (as the number of the observations inside an interval increases). However, for large tex2html_wrap_inline1958, the amplitude of the smoothed function decreases tex2html_wrap_inline2648. Thus exists an "optimal" value of tex2html_wrap_inline2650 where "signal/noise" reaches its maximum. For continuous approximation, this occurs at tex2html_wrap_inline2652 0.1855, 0.2895, 0.4123, 0.5450 for 4 fits with corresponding tex2html_wrap_inline2654 0.7885, 0.7843, 0.8749, 0.8716 and a value tex2html_wrap_inline2656 0.3396, 0.3531, 0.37450, 0.38353. This factor decreases by tex2html_wrap_inline2658 times at tex2html_wrap_inline2618 0.0608, 0.0936, 0.1589, 0.2093, i.e. there is a long enough nearly "standstill", if "signal/noise" is plotted vs. tex2html_wrap_inline2380.

In other words, approximations by a parabola give better results than that by a constant, the "optimal" fit corresponding to "wp" allowing to use more wide intervals than other fits.

Estimates of extremum time errors are possible only for m>0, as
equation1049
Minimum of this function for fits "up" and "wp" occurs at tex2html_wrap_inline2618 0.5515 and 0.7440 with corresponding tex2html_wrap_inline2668 0.6725 and 0.65669 and values 0.15929 and 0.2928 tex2html_wrap_inline2670. Here an accuracy estimate for "wp" is worse than for "up". But "up" fit is a discontinuous fit, and sometimes its derivative may be infinite. The accuracy estimate is twice larger than an "optimal" value at tex2html_wrap_inline2618 0.2712, 0.3612.

For small noise tex2html_wrap_inline2008 which is comparable with tex2html_wrap_inline2676 and/or discrete signal one may determine numerically the parameter tex2html_wrap_inline2560 optimizing a fixed characteristic.


next previous
Up: Method of running

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr