Sine fitting

Determine parameters of a noisy sine function

sineFit is a function to detect the parameters of a noisy sine curve, even less than one period long.
It requires only x and y values and no additional parameters as input.
It is tested with R2016a and R2020a.
The mean calculation time is on my PC 13 ms with a maximum of 2400 ms.

Syntax:
[SineParams]=sineFit(x,y,optional)
optional: plot graphics if ommited. Do not plot if 0
Input:
x and y values, y = offs + amp * sin(2pi * f * x + phi) + noise
Output:
SineParams(1):offset (offs)
SineParams(2): amplitude (amp)
SineParams(3): frequency (f)
SineParams(4): phaseshift (phi)
SineParams(5): MSE , if negative then SineParams are from FFT

Method:
This is a brief and not exact description of the program flow.
• Estimate the offset by the mean of all y values.
• Build the FFT with heavy zero padding.
• Take the frequency, amplitude and phase of the largest FFT peak.
If the frequency is at the Nyquist limit or the period is less than one, add extra frequencies for evaluation.
• Take those values as initial values for the regressions.
• Take the resulting MSE as rating.
• Exclude results above Nyquist frequency.
• Depending on the number of samples and the MSE, set a limit for an accepted amplitude in relation to the FFT amplitude.
• If the amplitude from regression is higher than the accepted amplitude, take the FFT parameters.

A demonstration GUI is included.
For more information read the "ReadMe.pdf".