Question

在以下稀疏矩阵中：

A=[1 1 1 3];
C = sparse(A',1:length(A),ones(length(A),1),4,4);
C =

   (1,1)        1
   (1,2)        1
   (1,3)        1
   (3,4)        1

>>full(C)

ans =

     1     1     1     0
     0     0     0     0
     0     0     0     1
     0     0     0     0

如何计算每行中非零元素的平均值？我无法在这些稀疏矩阵上使用 matlab 的内置均值函数。我找到了 this similar question 并且可以将其应用于我的问题

[row, ~, v] = find(C);
K>> rowmean = accumarray(row, v, [], @mean);
K>> rowmean 

rowmean =

     1
     0
     1

但是，我希望最后一行的值为零，而不是从答案中删除该行。

Answer 1

方法一

您可以将输出大小指定为accumarray的第三个输入：

[row, ~, v] = find(C);
rowmean = accumarray(row, v, [size(C,1), 1], @mean);

如果需要，您可以使用 accumarray 的第六个输入来获得 sparse 输出：

rowmean = accumarray(row, v, [size(C,1), 1], @mean, 0, true);

方法 2

你可以做一些像这样更简单的事情。结果是sparse：

rowmean = sum(C, 2) ./ sum(C~=0, 2); % mean of nonzeros in each row, manually
rowmean(isnan(rowmean)) = 0; % replace NaN (resulting from 0/0) by 0

稀疏矩阵每行非零元素的均值

1 个答案:

方法一

方法 2