can someone explain this code?

4 views (last 30 days)
Bhargavi
Bhargavi on 28 Jun 2012
x1= (Subject5_leftfall_1(:,44)-Subject5_leftfall_1(10,44));
x1_Y=find(x1 >= mean(x1(10:300))+3*std(x1(10:300)));
plot(x1(x1_Y(1):length(x1)),'g')
S5_Peak_left_1= max(x1(x1_Y(1):x1_Y(1)+70));
S5_index_peak_left_1=find(x1(x1_Y(1):x1_Y(1)+70)==S5_Peak_left_1)
I don't know what x1_Y is trying to do. Also, I don't know what the last two lines of the code are trying to do.

Accepted Answer

Geoff
Geoff on 28 Jun 2012
x1_Y is an array of indices. You're doing a boolean test on x1(10:300) to determine if a point is an outlier. That is, it is larger than (or equal to) 3 standard deviations over the mean. Actually, this is only testing outlying peaks - it's not testing standard deviations below the mean. Anyway, that's what find does. It takes a boolean expression and returns the index of every true value. So, if I had an array:
x = [ 0 1 1 0 0 0 1];
Then find(x) would return:
ans = [2 3 7]
Now, the last two lines are looking for the largest x1 value between the first detected outlier and 70 positions ahead. It then asks for the indices of any value equal to that in the same range. However, this is a little bogus. If there are two identical maximum peaks there, then find will return an array. I think this code intends to get a single index for the largest value. This is more correct:
[S5_Peak_left_1, S5_index_peak_left_1] = max( x1( x1_Y(1):x1_Y(1)+70 ) );
This code might still be wrong, depending on the reason for searching that 71-element range. Remember the index x1_Y is indexing into a slice of x1: (10:300), so if x1_Y(1) is 10, that means the actual position of that outlier in x1 is 20.
Hope that helps a little.
  1 Comment
Bhargavi
Bhargavi on 29 Jun 2012
Thanks a bunch Geoff. Also, there are no data on the first 9 lines on my data sheet. It just some frequency,measurement rate etc. This really helped. I had done it in the past, but since I hadn't commented, I sort of lost track of what was going on there.

Sign in to comment.

More Answers (0)

Categories

Find more on Get Started with MATLAB in Help Center and File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!