在 Matlab 中查找持续时间值的时间线答案

【问题标题】：Find timeline for duration values in Matlab在 Matlab 中查找持续时间值的时间线
【发布时间】：2016-03-17 14:52:34
【问题描述】：

我有以下时间序列：

b = [2 5 110 113 55 115 80 90 120 35 123];

b 中的每个数字都是一个数据点。我从b 计算了持续时间值。持续时间由b 内大于或等于 100 的所有数字表示，并连续排列（所有其他数字都被丢弃）。一个小于 100 的数字的最大间隙是允许的。持续时间的代码如下所示：

 N = 2;     % maximum allowed gap     
 duration = cellfun(@numel, regexp(char((b>=100)+'0'), [repmat('0',1,N) '+'],    'split'));

为b 提供以下持续时间值：

duration = [4 3];

我想为duration 中的每个值找到b 中的位置（时间线）。接下来，我想用零替换位于duration 之外的其他位置。结果将如下所示：

result = [0 0 3 4 5 6 0 0 9 10 11];

如果有人可以提供帮助，那就太好了。

【问题讨论】：

我从 b 计算了持续时间值发布你的代码，所以它也许可以扩展来做你想要的其他部分

标签： arrays matlab duration

【解决方案1】：

原始问题的答案：最多有一个值低于 100 的模式

这是一种使用正则表达式来检测所需模式的方法。我假设一个值 =100 之间（而不是之后）。所以模式是：一个或多个值 >=100，其中一个可能的值

b = [2 5 110 113 55 115 80 90 120 35 123]; %// data
B = char((b>=100)+'0'); %// convert to string of '0' and '1'
[s, e] = regexp(B, '1+(.1+|)', 'start', 'end'); %// find pattern
y = 1:numel(B);
c = any(bsxfun(@ge, y, s(:)) & bsxfun(@le, y, e(:))); %// filter by locations of pattern
y = y.*c; %// result

这给了

y =
     0     0     3     4     5     6     0     0     9    10    11

对已编辑问题的回答：在 100 以下的一行中最多有 n 个值的模式

正则表达式需要修改，它必须作为n的函数动态构建：

b = [2 5 110 113 55 115 80 90 120 35 123]; %// data
n = 2;
B = char((b>=100)+'0'); %// convert to string of '0' and '1'
r = sprintf('1+(.{1,%i}1+)*', n); %// build the regular expression from n
[s, e] = regexp(B, r, 'start', 'end'); %// find pattern
y = 1:numel(B);
c = any(bsxfun(@ge, y, s(:)) & bsxfun(@le, y, e(:))); %// filter by locations of pattern
y = y.*c; %// result

【讨论】：

我编辑了它。实际上是b>=100。如果b<100中有多个连续数字，则在计算时长时将忽略它们。
@Bowecho 已更正
如何在值 >= 100 之间从一个值
@Bowecho 这更棘手。这取决于您是允许所有这些值一次运行还是分散
我不知道你说的单次运行或分散是什么意思。

【解决方案2】：

这是另一个解决方案，不使用regexp。它自然地推广到任意间隙大小和阈值。不确定是否有更好的方法来填补空白。 cmets中的解释：

% maximum step size and threshold
N = 2;
threshold = 100;
% data
b = [2 5 110 113 55 115 80 90 120 35 123];

% find valid data
B = b >= threshold;
B_ind = find(B);
% find lengths of gaps
step_size = diff(B_ind);
% find acceptable steps (and ignore step size 1)
permissible_steps = 1 < step_size & step_size <= N;
% find beginning and end of runs
good_begin = B_ind([permissible_steps, false]);
good_end = good_begin + step_size(permissible_steps);
% fill gaps in B
for ii = 1:numel(good_begin)
    B(good_begin(ii):good_end(ii)) = true;
end
% find durations of runs in B. This finds points where we switch from 0 to
% 1 and vice versa. Due to padding the first match is always a start of a
% run, the last one always an end. There will be an even number of matches,
% so we can reshape and diff and thus fidn the durations
durations = diff(reshape(find(diff([false, B, false])), 2, []));

% get positions of 'good' data
outpos = zeros(size(b));
outpos(B) = find(B);

【讨论】：