How to detect variations in data?

2019 12 23
2 个回答
回答已采纳
更新时间：2019 12 26
18 次查看（30 天）

0 个投票

Hello everyone,
I would like to determine the location at which there is data variation (like the highlighted points in the attached image).
Kindly find the attached code + .mat file.
Each sub-plot represents a variable on y-axis and the time on x-axis. Usually, the variation appear in several variables (plots) like this one.
data .jpg
Could you suggest the best method to do so?
and How can I determine the starting and ending points of that variation?
I appreciate your help.

请先登录，再回答此问题。

 采纳的回答

Image Analyst
Image Analyst 2019-12-24

1 个投票

You might try movstd() to do a sliding window standard deviation. Then threshold it to find elements where the standard deviation is what you deem "high". See attached code where I did it for your "B" signal. Adapt it to do any other signals you want.
0000 Screenshot.png

更多回答（1 个）

Ridwan Alam
Ridwan Alam 2019-12-23

1 个投票

Depends on how complex of a solution you are looking for. There are algorithms for change point detection that you can try.
https://www.mathworks.com/help/signal/ref/findchangepts.html
Or, you can have a sliding window (your defined size) over the signal and look for specific features (e.g mean crossing rate, variance, etc. depending your signal type) to find the windows that contain a "change".
Hope this helps.

9 个评论
显示 7更早的评论 隐藏 7更早的评论

Mohamed Nedal
Mohamed Nedal 2019-12-23
Can you show me an example plz?
It would make it much easier to understand.
Ridwan Alam
Ridwan Alam 2019-12-23
Please check the link above, there are some examples. Try on your own data, and if there is any difficulty, please feel free to ask about that.
Mohamed Nedal
Mohamed Nedal 2019-12-24
编辑：Mohamed Nedal 2019-12-24
Yes, it gives what I need for now. Thanks!
But I need to plot several findchangepts plots as subplots in a single figure. I found it can't be done unfortunately.
Do you know how to do that?
Ridwan Alam
Ridwan Alam 2019-12-24
编辑：Ridwan Alam 2019-12-24
Why can’t it be done? Please show your code and I will try to debug.
Mohamed Nedal
Mohamed Nedal 2019-12-25
Here's what I was trying.
% define 3 column vectors
x1 = rand(30,1);
x2 = rand(30,1);
x3 = rand(30,1);
% first subplot
subplot(3,1,1)
findchangepts(x1,'MaxNumChanges',5);
% second subplot
subplot(3,1,2)
findchangepts(x2,'MaxNumChanges',5);
% third subplot
subplot(3,1,3)
findchangepts(x3,'MaxNumChanges',5);
Image Analyst
Image Analyst 2019-12-25
Did you also give my code a try?
Mohamed Nedal
Mohamed Nedal 2019-12-26
Yes, I tried it and it works. Thanks!
But it doesn't depend on findchangepts function, which is fine, I just would like to compare your code's result with that of findchangepts results.
Image Analyst
Image Analyst 2019-12-26
I think findchangepts() locates points where the data goes from smooth to something drastically different, or vice versa. So it basically finds corners or elbows in the data.
movstd() on the other hand gives a high signal wherever the data has high variability regardless if the signal looks like a corner or elbow there.
Mohamed Nedal
Mohamed Nedal 2019-12-26
Okay, thanks a lot.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Signal Operations 的更多信息

产品

版本

R2017b

标签

Community Treasure Hunt

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

Start Hunting!

Translated by