fitting critically damped simple harmonic oscillator

Four-parameter fit to the equation A*(1 + t/T2).*exp(-t/T1) + offset
431 Downloads
Updated 9 Feb 2012

View License

[A, T1, T2, offset, Ssq] = fit_critical_damping( t, y, options )
fitting critically damped oscillations to the equation
y(t) = A (1 + t/T2) exp( -t/T1 ) + offset
where A is the amplitude, T1 is a time constant, T2 is the decay time, along with
the offest. Ssq is the sum of residuals (i.e. fitting error).

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 ) = 4 (i.e. there are four columns);

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

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

Note: this script uses LMFnlsq from MATLAB central file exchange
sak@wpi.edu 2/8/12

Cite As

Stephan Koehler (2024). fitting critically damped simple harmonic oscillator (https://www.mathworks.com/matlabcentral/fileexchange/35015-fitting-critically-damped-simple-harmonic-oscillator), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R13
Compatible with any release
Platform Compatibility
Windows macOS Linux
Categories
Find more on Polynomials in Help Center and MATLAB Answers

Community Treasure Hunt

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

Start Hunting!
Version Published Release Notes
1.0.0.0