arrays - 根据条件在数组/向量中添加单元格 - Matlab

Question

我有一个数据向量 448X2，其中第一列对应于距离(从 -50 到 224)，第二列包含与每个距离关联的强度。

无论方向(正/负)如何，我都需要添加距离上的所有强度。例如，我想添加 1 的间隔。因此，如果在 [-1 1] 距离范围内有 5 个值，我需要添加它们的所有强度。然后，从 [-2 -1] & [1 2] 距离等等，等等，...

我首先将数据从 0 缩短到 224

datasort = sortrows(abs(data));

它取走了负数。现在我的数据是按距离升序排列的，我不知道如何在距离间隔内添加强度。

这些是我的 datasort 向量的前 10 行:

0.119865417480469   0.0115106217563152
0.669918060302734   0.356697350740433
0.909893035888672   0.659623980522156
1.45951080322266    1.69302666187286
1.70020294189453    2.30387997627258
2.24539566040039    3.98834538459778
2.49082183837891    4.94714736938477
3.02403259277344    7.16749382019043
3.28177261352539    8.59203147888184
3.79549789428711    11.1872768402100

我想以这样的方式结束:

1   1.02780000000000
2   3.99690000000000
3   8.93550000000000
4   26.9468000000000

其中第一列现在包含 1-224 (1:224) 的离散距离值，第二列包含距离范围内所有强度的总和。

希望我说得足够清楚。谢谢

更新:

我想我可以做到:

datasort(:,1) = round(datasort(:,1));

这给了我(前 10 分)，

0   0.0115106217563152
1   0.356697350740433
1   0.659623980522156
1   1.69302666187286
2   2.30387997627258
2   3.98834538459778
2   4.94714736938477
3   7.16749382019043
3   8.59203147888184
4   11.1872768402100

我仍然不知道如何将相同距离的强度相加!!!

Answer 1

我会首先使用 histc 来合并数据的第一列，然后将 histc 返回的索引传递给 accumarray对第二列进行操作。请注意，您甚至不需要为此对数据进行排序，尽管您仍然需要调用 abs 来删除负数。

>> data = [0.119865417480469   0.0115106217563152
0.669918060302734   0.356697350740433
0.909893035888672   0.659623980522156
1.45951080322266    1.69302666187286
1.70020294189453    2.30387997627258
2.24539566040039    3.98834538459778
2.49082183837891    4.94714736938477
3.02403259277344    7.16749382019043
3.28177261352539    8.59203147888184
3.79549789428711    11.1872768402100]

data =

         0.119865417480469        0.0115106217563152
         0.669918060302734         0.356697350740433
         0.909893035888672         0.659623980522156
          1.45951080322266          1.69302666187286
          1.70020294189453          2.30387997627258
          2.24539566040039          3.98834538459778
          2.49082183837891          4.94714736938477
          3.02403259277344          7.16749382019043
          3.28177261352539          8.59203147888184
          3.79549789428711            11.18727684021

>> [~, idx] = histc(data(:,1), [0:224]) 

idx =

     1
     1
     1
     2
     2
     3
     3
     4
     4
     4

>> result = accumarray(idx, data(:,2))

result =

           1.0278319530189
          3.99690663814544
          8.93549275398255
          26.9468021392823

>> 

希望对您有所帮助!

编辑:

只是为了确认这在不对数据进行排序的情况下确实有效

>> % Unsort your data snippet to illustrate
>> data = data([randperm(length(data))],:)

data =

          3.79549789428711            11.18727684021
          1.45951080322266          1.69302666187286
         0.909893035888672         0.659623980522156
          2.24539566040039          3.98834538459778
          1.70020294189453          2.30387997627258
         0.119865417480469        0.0115106217563152
         0.669918060302734         0.356697350740433
          3.28177261352539          8.59203147888184
          3.02403259277344          7.16749382019043
          2.49082183837891          4.94714736938477

>> [~, idx] = histc(data(:,1), [0:224])

idx =

     4
     2
     1
     3
     2
     1
     1
     4
     4
     3

>> result = accumarray(idx, data(:,2))

result =

           1.0278319530189
          3.99690663814544
          8.93549275398255
          26.9468021392823

>>