How to calculate average in a large matrix?

Question

afrya 2013-12-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/108268-how-to-calculate-average-in-a-large-matrix

评论： Jos (10584) 2013-12-3

Dear all,

I'm working with a large amount of data in a project.For each row of my matrix , I calculate the average. Here is an example that help you to understand my problem

   y1 y2 y3 y4 y5                                      
x1 2  4  4  6  7      Av1=4,6          
x2 1  2  3  4  5      Av2=3
x3 1  2  3  4  5      Av3=3
x4 1  2  4  5  6      Av4=3,6
x5 2  5  6  8  9      Av5=6
x6 4  7  8  4  9      Av6=6,4
x7 2  3  4  8  7      Av7=4,8
x8 1  2  4  6  4      Av8=3,4

I want to make my programm user friendly, that means the user enter the number of average that he wants.For instance, if the user enter the number 2 , the programm calculates the average of each 2 rows

   y1 y2 y3 y4 y5                                      
x1 2  4  4  6  7      NewAv1=Av1+Av2=3,8      
x2 1  2  3  4  5   
x3 1  2  3  4  5      NewAv2=Av3+Av4=3,3
x4 1  2  4  5  6    
x5 2  5  6  8  9      NewAv3=Av5+Av6=6,2
x6 4  7  8  4  9      
x7 2  3  4  8  7      NewAv4=Av7+Av8=4,1
x8 1  2  4  6  4

Thank you in advance

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Andrei Bobrov 2013-12-3

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/108268-how-to-calculate-average-in-a-large-matrix#answer_117073

编辑：Andrei Bobrov 2013-12-3

在 MATLAB Online 中打开

xy = [2  4  4  6  7     
1  2  3  4  5   
1  2  3  4  5
1  2  4  5  6    
2  5  6  8  9
4  7  8  4  9      
2  3  4  8  7
1  2  4  6  4];
k = 2;
i0 = zeros(size(xy,1),1);
i0(1:k:end) = 1;
out = accumarray(repmat(cumsum(i0),size(xy,2),1),xy(:),[],@mean);

ADD

xy2 = xy(1:floor(size(xy,1)/k)*k,:);
s = size(xy2);
i0 = zeros(s(1),1);
i0(1:k:end) = 1;
out = accumarray(repmat(cumsum(i0),s(2),1),xy2(:),[],@mean);

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

Andrei Bobrov 2013-12-3

see ADD

Jos (10584) 2013-12-3

My suggestion seems more efficient (no cumsum, for instance) ...

请先登录，再进行评论。

Answer 2

Azzi Abdelmalek 2013-12-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/108268-how-to-calculate-average-in-a-large-matrix#answer_116955

在 MATLAB Online 中打开

%Example--------------
A=rand(10)
p=3
%-----------------------
[n,m]=size(A);
q=ceil(n/p)*p
B=nan(q,m)
B(1:n,:)=A
ii2=p:p:q
ii1=[1 ii2(1:end-1)+1]
out=arrayfun(@(x) nanmean(nanmean(B(ii1(x):ii2(x),:))),1:numel(ii2))'

4 个评论
显示 2更早的评论隐藏 2更早的评论

afrya 2013-12-2

the size of the matrix is: 3500*579

Azzi Abdelmalek 2013-12-2

编辑：Azzi Abdelmalek 2013-12-2

在 MATLAB Online 中打开

%Example--------------
A=rand(3500,379);
p=3;
%-----------------------
tic
[n,m]=size(A);
q=ceil(n/p)*p;
B=nan(q,m);
B(1:n,:)=A;
B=reshape(B',p*m,[]);
out=nanmean(B);
toc

result

Elapsed time is 0.043406 sec

请先登录，再进行评论。

Answer 3

Jos (10584) 2013-12-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/108268-how-to-calculate-average-in-a-large-matrix#answer_116964

在 MATLAB Online 中打开

Here is a trick, using clever indexing with accumarray:

% data
  A = rand(10,5) ;
  K = 3 ; % average of so many rows
% engine
  ix = repmat(floor((0:(size(A,1)-1))/K),size(A,2),1).' + 1 % indices
  AV = accumarray(ix(:), A(:),  [ix(end) 1], @mean) % averages
% check
  abs(AV(1) - mean(reshape(A(1:K,:),1,[]))) < eps