next previous
Up: Method of running

3. Derivatives of the smoothing function

According to definition (1), the smoothing function tex2html_wrap_inline1960 coincides with the function tex2html_wrap_inline2358 at points t=t0. However, this in not the case for the derivatives, i.e.
eqnarray547
Obviously, for the s-th derivative of a general parameter X,
equation565

Equation (31) allows to estimate accuracy tex2html_wrap_inline2366 of the s-th derivative tex2html_wrap_inline2370:
eqnarray569
Particularly, if t=t0 and polynomial basic functions (6),
equation582
and thus
equation587

Much more complicated is the determination of the derivatives by the argument t0. Most important are first and second derivatives, especially at the extrema. Let's determine the vectors h at a moment tex2html_wrap_inline2378. For small tex2html_wrap_inline2380, one may expand vectors into series restricting maximum order to 2:
equation592
And similarly for tex2html_wrap_inline2382 tex2html_wrap_inline2280 hk and tex2html_wrap_inline2388 Hereafter the last index s=0,1,2 corresponds to a coefficient at tex2html_wrap_inline2392 For the coefficient C0:
eqnarray602
The coefficients of the power series (35) for parameter X may be written as tex2html_wrap_inline2398, as all differences (t-t0) become tex2html_wrap_inline2402. However, for tex2html_wrap_inline2026 it is more suitable to use the expressions
eqnarray645
For further study of the polynomial fits, we will measure times in units of tex2html_wrap_inline1958, practically using dimensionless units tex2html_wrap_inline2408. Introducing the sums tex2html_wrap_inline2410 one may easily obtain
eqnarray662
for "unweighted" parabolic fits, and, for the weights (5) the matrices tex2html_wrap_inline2412 tex2html_wrap_inline2414 tex2html_wrap_inline2416 are equal to
eqnarray679
Other matrices tex2html_wrap_inline2418, tex2html_wrap_inline2262, tex2html_wrap_inline2258 may be consequently determined by using above mentioned expressions.

To determine extremum of the smoothing function, one has to solve an equation
equation693
After determining the root t0 for the given signal values xk, one has to obtain an accuracy estimate of it. Assuming small variations tex2html_wrap_inline2428 of the moment of the extremum caused by small tex2html_wrap_inline2430, one may write a linearized equation
equation699
or
equation707
where the first tex2html_wrap_inline2432 and second tex2html_wrap_inline2434 derivatives of the smoothing function tex2html_wrap_inline1950 are evaluated at argument t0. Assuming again tex2html_wrap_inline2440 tex2html_wrap_inline2442, one may obtain
equation719

As an illustration of the derived above expressions, we show in Fig. 2 (click here) the dependence of h[C*,k] on k for 19-point "wp" approximation.

  figure726
Figure 2: Projective 19-point vectors h[C*,k] for polynomial "wp" fits

Statistical properties of the test functions used for the period determination by using the moments of "characteristic events" are studied earlier (Andronov 1987, 1991).


next previous
Up: Method of running

Copyright by the European Southern Observatory (ESO)
This email address is being protected from spambots. You need JavaScript enabled to view it.