"....the cycle is defined as when the timeseries y value passes the positive threshold, then the negative threshold and then the positive threshold again in the time domain....." - To achieve this, you can try the following code -
% Example signal
y = [0, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2];
positive_threshold = 1;
negative_threshold = -1;
state = 'below_positive';
cycle_count = 0;
for i = 1:length(y)
switch state
case 'below_positive'
if y(i) > positive_threshold
state = 'above_positive';
end
case 'above_positive'
if y(i) < negative_threshold
state = 'below_negative';
end
case 'below_negative'
if y(i) > positive_threshold
state = 'above_positive';
cycle_count = cycle_count + 1;
end
end
end
disp(['Number of cycles: ', num2str(cycle_count)]);
Example of a cycle in the code -
- Initial State: Below the positive threshold.
- First Positive Threshold Crossing: The signal goes from 0 to 1 to 2 (crosses the positive threshold).
- Negative Threshold Crossing: The signal goes from 2 to 1 to -1 to -2 (crosses the negative threshold).
- Second Positive Threshold Crossing: The signal goes from -2 to -1 to 1 to 2 (crosses the positive threshold again).
This sequence completes one cycle. Hope it helps!
