fitting of a sinusoid

1 visualización (últimos 30 días)
Anik Faisal
Anik Faisal el 13 de Jun. de 2019
Respondida: Jesús Zambrano el 4 de Ag. de 2019
Hi Everyone!
I have generated energy change (gammas.mat) of a simulation cell as a function of rigid displacement (ubp.mat) of upper half of the simulation cell with respect to to the lower half from atomistic calculations. I was trying to fit the energy vs displacement data to continous function of the displacement. I have tried using the curve fitting tool and the sum of sine model with 6 terms shows good agreement except at zero displacement where the fitting produces a negative value for the energy which is not physical since no change in energy before the displacement is initiated. I have attached a screenshot of the curve fitting tool's attempt.
I am looking for suggestions if I the data fitting can be done with a zero energy at zero displacement. Thanks!
Screen Shot 2019-06-12 at 10.31.09 PM.png

Respuestas (1)

Jesús Zambrano
Jesús Zambrano el 4 de Ag. de 2019
Hi Anik,
Try to use the Piecewise Cubic Hermite Interpolating Polynomial (PCHIP) function. The function pchip does not include oscillations between sample points. You can see the comparison between pchip and spline in the link below, with some examples:
Hope this asnwer can help you to solve your question.

Categorías

Más información sobre Get Started with Curve Fitting Toolbox en Help Center y File Exchange.

Etiquetas

Productos

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by