fitting decaying oscillations

From inputs t, y the best fit to the decaying oscillation is determined.


Updated Thu, 09 Feb 2012 17:17:05 +0000

View License

[A, T1, T2, phase, offset, Ssq] = fit_decay_oscillations( t, y, options )
fitting decaying oscillations of the form
y(t) = A exp( -t/T1 ) cos( t / T2 + phi) + offset
where A is the amplitude, T1 is the decay time, T2 is the period, along with
phase and offest. Ssq is the sum of residuals.

options are 'notify' to show progress and 'plot' to give a plot
you can also supply a guess or series of guesses via a row vector or
matrix, m, where size( m, 2 ) = 5;

if you just run program without arguments it will do an example

examples are
fit_decay_oscillations( t, y, 'notify', 'plot')
[A, T1, T2, phase, offset] = fit_decay_oscillations( t, y, [1 2 3 4 5])

Note: this script uses LMFnlsq from MATLAB central file exchange 2/5/12

Cite As

Stephan Koehler (2023). fitting decaying oscillations (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R12
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

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

Start Hunting!
Version Published Release Notes

another bug fix

speed improvement, and fixed calculation of Jacobian

used the Jacobian matrix to speed calculations.

bug fix