Filter löschen
Filter löschen

Find value at which a line changes its slope

16 Ansichten (letzte 30 Tage)
Lina Koronfel
Lina Koronfel am 28 Mai 2023
Kommentiert: Lina Koronfel am 23 Jul. 2023
How can I find the x value of the point during which the line (significantly) changes direction? I tried inflection point but it doesn't work in this case. Basically, I want to find the exact value during which the line changes its slope (arrow). Thank you

Akzeptierte Antwort

Image Analyst
Image Analyst am 28 Mai 2023
See if my attached piecewise linear fit identifies it correctly for you.
If not, then attach your data after reading this: TUTORIAL: How to ask a question (on Answers) and get a fast answer
And as an alternative:
  1 Kommentar
Lina Koronfel
Lina Koronfel am 23 Jul. 2023
Thank you for the answer, I wrote a code that incorporated all suggestions from @John D'Errico @Torsten and @Image Analyst and it is working perfectly. Thank you for all your help!

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (2)

Torsten am 28 Mai 2023
Verschoben: Torsten am 28 Mai 2023
Plot gradient(y)./gradient(x) over x to get an approximation of the derivative of your curve and see if this curve can give you more information.

John D'Errico
John D'Errico am 28 Mai 2023
Bearbeitet: John D'Errico am 28 Mai 2023
There is NO exact point where that can be known to happen. Your curve is noisy.
At best, you can estimate where the curve starts to drift away from being a constant function. And of course, we lack any data from you, so it is difficult to give an example.
It seems your curve is essentially zero for a while, so why not just look for the point where y deviates from zero by some tolerance? That is the point you are essentially asking to find.
We can compute the approximate noise in your curve simply as:
noiseEst = std(diff(y))/sqrt(2);
This assumes the noise is additive and iid, and at least vaguely gaussian. It also assumes the signal underneath the noise in the curve is slowly varying.
Having done that, just find the point where
ind = find(abs(y) > 3*noiseEst,1,'first')
That gives you the first point where y deviates from zero by more than 3 sigma.


Community Treasure Hunt

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

Start Hunting!

Translated by