3

我无法组合我的 Matlab“数据”变量的重复元素。我可以使用 unique 和 sort 轻松组合这些值。

[sorted,idx] = sort(data);
[~,ij] = unique(sorted,'first');
Indx = (sort(idx(ij)));

但是,通过这样做,我组合了所有重复的值。我真正想做的是只组合重复元素的组。例如采取这个:

data = [1;1;1;2;2;2;3;3;3;4;4;4;4;4;3;3;2;2;2;2;1;1;1;1;4;4;4;4;]

合并重复的元素组:

data = [1;2;3;4;3;2;1;4;]

我需要在仍然保持顺序的情况下组合重复元素组。返回索引也很有帮助,因为我需要根据组合索引对另一个变量中的数据进行平均。

例如:

data  = [1;1;1;2;2;2;3;3;3;4;4;4;4;4;3;3;2;2;2;2;1;1;1;1;4;4;4;4;]
data2 = [7;2;4;5;3;4;6;8;5;3;5;7;4;2;4;6;8;4;3;6;7;8;4;2;9;3;2;0;]

dataCombined = [1;     2;  3;    4;   3;  2;     1;     4;   ]
data2average = [4.33;  4;  6.33  4.2  5;  5.25;  5.25;  3.5; ]

任何人都可以提出建议吗?


解决方案:

谢谢大家的答案。MZimmerman6 的解决方案对我来说效果很好。我想展示我做了什么来平均“data2”数组中的值。

data = [1;1;1;2;2;2;3;3;3;4;4;4;4;4;3;3;2;2;2;2;1;1;1;1;4;4;4;4;];
data2 = [7;2;4;5;3;4;6;8;5;3;5;7;4;2;4;6;8;4;3;6;7;8;4;2;9;3;2;0;];
change = diff(data)~=0;
indices = [1,find(change)'+1];
compressed = data(indices)';


numberOfRepeatingGroups = size(indices);


for i=1:numberOfRepeatingGroups(1,2)


  if(i == 1)   

      dataToAverage = data2(indices(1,1):(indices(1,2)-1));

  elseif (i == numberOfRepeatingGroups(1,2))

       dataToAverage = data2(indices(1,i):end);

  else

       dataToAverage = data2(indices(1,i):(indices(1,(i+1))-1));

  end

       data2Averaged(1,i) = mean(dataToAverage(:));

end   


data2Averaged =

4.3333    4.0000    6.3333    4.2000    5.0000    5.2500    5.2500    3.5000
4

2 回答 2

4

您可以使用导数来查找数据数组中的波动,这表明分组发生了变化。导数不为 0 的任何地方都会发生变化,无论是正的还是负的。找到这些变化发生的地方,然后抓取相应的索引。像下面的东西。

data = [1;1;1;2;2;2;3;3;3;4;4;4;4;4;3;3;2;2;2;2;1;1;1;1;4;4;4;4;];
change = diff(data)~=0;
indices = [1,find(change)'+1];
compressed = data(indices)';

结果将是

compressed =
     1     2     3     4     3     2     1     4

当然,您也可以根据indices需要使用该变量。

注意 在第三行,我们添加索引 1,因为从技术上讲,数组的开头是一个变化,然后我们在 find 命令中添加 1,因为我们在导数上使用 find,所以返回的change数组将比原始数组短 1 .

于 2013-09-09T16:02:51.843 回答
2

我永远不会停止从File Exchange 推荐这个运行长度编码/解码实用程序:rude() .

% Run-length encode preserving order
[len,val] = rude(data);
len =
     3     3     3     5     2     4     4     4
val =
     1     2     3     4     3     2     1     4

现在,要计算平均值,首先用 重新标记每个子序列rude(),然后使用accumarray()

% Decode and re-label each subsequence
subs = rude(len,1:numel(len))';

% Take average on each re-labelled subsequence
accumarray(subs,data2,[],@mean)
ans =
    4.3333
    4.0000
    6.3333
    4.2000
    5.0000
    5.2500
    5.2500
    3.5000
于 2013-09-09T16:31:19.387 回答