Marix manipulation

Question

Charles 2011-8-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/13138-marix-manipulation

How do I perform this operation? I have a data set with recurring values at a given co-ordinate. I need to sum values of the repeating coordinates. I initially used a for loop to solve the problem but it is too slow. Here is a simple example of the problem:

val = rand(1);

 a = [1 1 1 val; 
      1 2 1 val; 
      1 1 1 val];

Since the 1st and 3rd rows have the same co-ordinates, I need to add their values, delete the entries and replace with the new summed value. I know it should be possible using sub2ind and a cumsum operation, but I haven't figure exactly how to combine them yet. If I haven't written an answer at the time you are reading this, it means I still haven't got the combination right. Thanks.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Pierre 2011-8-5

>> _Since the 1st and 3rd rows have the same co-ordinates, I need to add their values, delete the entries and replace with the new summed value._

I wonder, if anybody can understand what you're trying to do there... I can't. :S

Please post your actual solution (with the for loop), this would allow us to see what the output's supposed to look like, would you?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Daniel Shub 2011-8-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/13138-marix-manipulation#answer_17953

在 MATLAB Online 中打开

My solution has a for loop.

[b, i, j] = unique(a(:, 1:(end-1)), 'rows');
c = b;
for ii = 1:length(i)
  c(ii, 4) = sum(a(j==ii, 4));
end

My guess is that it can be eliminated (and potentially made faster) with a bsxfun or arrayfun. In order to optimize it, we would need to know the sizes and how often things repeat.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Charles 2011-8-5

ACCUMARRAY might be the better option as it doesn't involve a loop.

Daniel Shub 2011-8-5

While I am happy to have my answer accepted, my guess is that if you un-accept my answer (not even sure you can do that), once the people in the US wake up, you will get a bsxfun/arrayfun/accumarray solution.

请先登录，再进行评论。