Polynoms are not matching Fourier Fits for AD CMi (DSCT) - although its a good data set ???

Affiliation
Bundesdeutsche Arbeitsgemeinschaft fur Veranderliche Sterne e.V.(Germany) (BAV)
Fri, 03/04/2022 - 12:56

Hi,

I have a nice lightcurve of a DSCT named AD CMi. For determining the maximum i tried:

VSTAR Polynoms 9,12,15 degree

VStar Fourier Fit   3 grade

Bob Nelsons Minima27 Fourier Fit

Gnuplot Fourier Fit grade

Normally I take a mean of some of the nicest Polynoms.

 

Picture: V-C CMI AD gnuplot fourierfit.PNG  From the left, there are five arrows for the max.  (The solid blue fit is from the GNUPLOT Fourier Fit)

1st Arrow (from the left) is the VSTAR Polynom 15 degree

2nd Arrow (from the left) is the VSTAR Polynom 9 degree

3rd Arrow (from the left) is the Fourier Fit Bob Nelsons Minima 2.7

4th Arrow (from the left) is the Fourier VSTAR 3 grade

5th Arrow (from the left) is the Fourier FIT GNUPLOT 3 grade

 

QUESTION:

  • Why are the Fourierfits differing that much from the Polynoms?

 

  • How is the proper way to calculate a correct maximum?

 

see also:

https://www.aavso.org/LCGv2/index.htm?DateFormat=Julian&RequestedBands=&view=api.delim&ident=ad cmi&fromjd=2458913&tojd=2459643.047&delimiter=@@@

 

Here are my results:

JD: 2459624,37619, Mag: 9,184601 [Maximum] 9 Polynom

RMS: 0,010486

AIC: -1239,940579

BIC: -1213,595296

 

 

 

JD: 2459624,37359, Mag: 9,186395 [Maximum] 12 Polynom

RMS: 0,008823

AIC: -1281,591163

BIC: -1246,464119

 

 

JD: 2459624,37318, Mag: 9,188794 [Maximum] 15 Polynom

RMS: 0,008658

AIC: -1280,784898

BIC: -1236,876092

 

FOURIER FITS

0,38306 fourier MINIMA 2.7

0,3842 Fourier 3 grad VSTAR 2.21.3

0.389263 Fourier 3 grad GNUPOT

 

The DIFFERENCE form the different results is: 0,389263 - 0,37318 = 0.016083 HJD   are 0.385992 h = 23.15952 min

 

Affiliation
American Association of Variable Star Observers (AAVSO)
Fourier Fits for AD CMi

Bernhard:

I suspect that the value produced using Minima 2.7 is an outlier for some reason.   Your data produced the following results using Peranso 3.0.3.4: (Poly 6) 0.3744871 ±0.0014859 Brightness 9.188.   MAVKA (https://www.ta3.sk/caosp/Eedition/FullTexts/vol50no2/pp557-559.pdf) features a large number of fitting options which I've found to be very useful for HADS variables and totally eclipsing overcontact binaries.  In this case the best fit (0.0079) was obtained using the "parabolic spline approximation" in which the time of extremum was 0.37462 ±0.00123 at mag 9.19 ±0.003. 

The Fourier fits will depend on the starting and ending positions that have been selected for curve fitting.   I normally start and end the curve fit at half the amplitude but more importantly just be consistent with the software and fitting algorithm you are using for each variable star.

 

Kevin

Affiliation
Bundesdeutsche Arbeitsgemeinschaft fur Veranderliche Sterne e.V.(Germany) (BAV)
Thank you Kevin,For…

Thank you Kevin,

For looking into my data. Hmm.. Maybe Peranso is a try worth it...

I've downloaded MAVKA    https://uavso.org.ua/mavka/

It is sounding promising!

 

 InputFile=
max.dat
  OutFile=
max.da$
Parabolic spline approximation
linear coef 5
nonlinear coef 2
sigm 0.0121829292395493
time of extremum 2459624.3775955      -+0.000450692323739213
mag of extremum 9.18192711187491
-+0.00259315648316741
 

 

BernhardImage removed.

Bernhard,

I can comment…

Bernhard,

I can comment only on the VStar data. It seems that you may have selected the entire light curve on which to perform the polynomial fits. This is not necessary, since you are interested only in determining the time of maximum (TOM). I suggest that you experiment with selecting (displaying) a more limited parts of the light curve, for example about half of the ascending and descending limbs and the peak, then in the View menu select Filter from Plot. Fit the polynomials to this filtered dataset.

I then fit polynomials starting with 4th order, and increasing (5th, 6th, etc) until it is clear the AIC reaches a minimum value (e.g, if AIC values were -281, -283, -284, -282, then -284 is the minimum value). VStar gives an RMS error. This is the root mean square of the residuals, which is not an error estimation of the actual TOM. In order to obtain the latter, you need to use Peranso or some other software.

I use Peranso. If you note the range of HJD values you used for your polynomial fit in VStar, note the polynomial order, then use those values in Peranso, you will find a TOM (it should be identical to that obtained in VStar) and an error estimation of the time.

Just a note on AIC. I use this and not BIC in VStar because limited experiments I've done indicate that, after fitting polynomials to several different parts of a light curve around the peak or trough, noting the AIC and BIC, then noting the TOM, the range of times is smaller using AIC than it is for BIC. In other words, in that experiment, AIC was less sensitive to the actual part of the light curve selected for polynomial fitting than was BIC.

Roy