Rate of Change - Derivative of experimental data

78 views (last 30 days)
Does anyone have any advice on determining the instantaneous rate of change for strain or force data for multiple cycles? If there is a wasy or function to help with taking the derivative of a curve in MatLab? I've tried using the diff() function but the results give me a value in the thouands when the peak instantaneous loading rate I received from LabView is only in the hundreds. I can attach my current code if you would like to see it, but unfortanutely I can't post the experimental data I have.
  1 Comment
Sam Chak
Sam Chak on 24 Jun 2022
Perhaps considering to post only a fraction of the data, say 10 points, or 10% of the data (if not too large).

Sign in to comment.

Answers (2)

Star Strider
Star Strider on 24 Jun 2022
I usually do something like this —
dydx = gradient(y) ./ gradient(x);
This assumes the independent variable ‘x’ and the dependent variable ‘y’ are both vectors of the same length. See the documentation on gradient for details.
  4 Comments
Perri Johnson
Perri Johnson on 24 Jun 2022
I have used a lowpass filter to help smooth the data. The curves are relatively smooth but when the gauge is being unloaded and the signal settles before beginning a new cycle, there seems to be some noise there. I may need to figure out how to use a moving average to help smooth the data at those points.
Star Strider
Star Strider on 24 Jun 2022
There are several options, among which is filloutliers that detects outliers and interpolates them. Like all signal processing tasks, you will need to experiment to get the desired result.

Sign in to comment.


John D'Errico
John D'Errico on 24 Jun 2022
Edited: John D'Errico on 24 Jun 2022
If your data is smooth enough, then the direct gradient calls that Star shows are ok. The problem is, if there is any noise in your data, then using the finite differences from gradient may be a bad idea, as they will amplify the noise. Then you will benefit form a Savitsky-Golay style of filter. The idea is to use a moving polynomial approximation at each point to estimate the slope.
I posted movingslope on the File Exchange, which does exactly that, allowing you to choose the length of the moving filter, as well as the local polynomial order.
However, I recall it assumes your data curve is equally spaced in x. If you lack equally spaced data in x, then if you have any noise, you can still do something, but it will be way less efficient.

Categories

Find more on Signal Generation and Preprocessing in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!