MATLAB Answers

0

fitting of a sinusoid

Asked by Anik Faisal on 13 Jun 2019
Latest activity Answered by Jesús Zambrano on 4 Aug 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

  0 Comments

Sign in to comment.

Products

1 Answer

Answer by Jesús Zambrano on 4 Aug 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.

  0 Comments

Sign in to comment.