Matlab对于含有nan值矩阵求均值的误解

MATLAB中，对于矩阵，如果含有nan值（或者其他异常值）时，若对其中正常值如何进行求均值。如下

a=[1 2 3 0; 2 3 1 NaN; 1 2 NaN 3; NaN 1 2 3];

当我知道MATLAB中如nanmean函数时欣喜若狂，太好用了直接nanmean(nanmean(a))不就可以了吗？！！结果错了，大错特错，

如果用手动计算，那么我们很容易知道结果是：24/13=1.8462

而如果你用nanmean方法计算，你会发现结果是nanmean(nanmean(a))=1.8333，原因很简单，就是你两次nanmean方法时不对等性。正确的方法是先求出所有非nan值的和，然后除以所有非nan值的数量。正确方法一如下：

b=a;          %复制一份a

n=isnan(a);  %n为0,1矩阵，a中nan值位置为1，非nan值为0

b(n)=0;       %b矩阵中是nan值的负值为0，这样求和就不影响了，因为变为0了

c=sum(sum(b,2))./sum(sum(~n,2));%分子是求所有矩阵元素的和，原来的nan值为0，所以不影响求和结果。分母是求出所有非nan值的个数，这样就得到了矩阵均值c=1.8462

或者方法二以下代码也可以：

notNanValues=a(~isnan(a));

meanvalue=sum(notNanValues)./length(notNanValues); 或者直接 用   mean(notNanValues)

或者直接用 nanmean(a(:))