Segmenting sinusoidal time-series data

Ben
2016 2 2
1 个回答
回答已采纳
更新时间：2016 2 4
4 次查看（30 天）

0 个投票

I have some data that was collected by a mobile phone's gyroscope sensors. I would like to be able to chunk out one "repetition" of the motion. An example of a sample run is shown here (apparently the MathWorks image upload didn't work right): http://i.imgur.com/GuEgK2y.png
I would like to be able to determine how long a period is, and where they start (doesn't matter if the starting point is max or min amplitude or a center amplitude). I tried using xcorr, which resulted in a plot with some waves centered around 50, whereas manual inspection of the data shows a period of about 2000 between the first two peaks. What should I be looking into to figure out how to do this?

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

 采纳的回答

Star Strider
Star Strider 2016-2-2

0 个投票

Without your actual data, I can only guess at an approach. I would begin by detrending it (I used polyfit and polyval, there are several options), then use circshift to identify the zero-crossings, and from them the periods and frequencies:
t = 0:1.3E+4;
s = 60 + 5*exp(-5*t/1.3E+4) + sin(2*pi*t/3000); % Create Data
figure(1)
plot(t, s); % Original Data
p = polyfit(t, s, 3);
sd = s - polyval(p, t); % Detrend
xci = find(sd.*circshift(sd, [0 -1]) <= 0);
xci = xci(1:end-1); % Find Inices Of Zero-Crossings
figure(2)
plot(t, sd)
hold on
plot(t(xci), sd(xci), '+r')
hold off

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Ben
Ben 2016-2-4
Hmm. When I plot sd vs t (using my data), I see that the data has been adjusted as expected and has zero-crossings, but xci is an empty array.
Could you explain a little more what you're doing with find and circshift? I don't understand why you're shifting the array.
Star Strider
Star Strider 2016-2-4
If ‘xci’ is empty, that means your data are a column rather than a row vector. (This is the problem of not having your original data to work with.) The first approach is to switch the argument vector in the circshift call:
xci = find(sd.*circshift(sd, [-1 0]) <= 0);
Try it now.
The circshift call shifts the vector (in this instance, see the documentation for details) one position to the left (or up) (-1), then multiplies it by the unshifted vector. This creates negative values only where zero-crossings occur. The find call then locates the indices of them. It is possible to extend this to include a loop with an interp1 call at each iteration (for coarsely-sampled data) to get more accurate zero-crossing estimates. I kept it simple here.
Ben
Ben 2016-2-4
Ok, switching which dimension was shifted worked, thank you.

请先登录，再进行评论。

更多回答（0 个）

类别

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

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by