nanmean究竟是怎么样的一个求平均方法

zane

我在论坛搜索不了相关的帖子，所以才在此发帖求助的。
一般数组中存在nan的时候，大多数情况下会用nanmean来求平均值，但是如果先把不是非数的数据求和并计算其个数，然后用求和的结果除以不是非数的个数之后的值，与前者会是不一样。下面是我编的小程序，其中b=3.0833是用nanmean求出来的结果，c=3是用后者方法求出来的结果。若看客知道个中缘由的，望告知，先行谢过。:)
a=[1 2;3 nan;4 5;nan nan];
b=nanmean(nanmean(a));
[k,l]=size(a);
sum=0;h=0;
for i=1:k
    for j=1:l
        if isnan(a(i,j))==0
            sum=sum+a(i,j);
            h=h+1;
        end
    end
end
c=sum/h;

二爷名声在外

将b=nanmean(nanmean(a));
改写成
b1=nanmean(a)
b=nanmean(b1)
你就会明白了！

zane

二爷名声在外 发表于 2014-5-25 10:09
将b=nanmean(nanmean(a));
改写成
b1=nanmean(a)

我根据nanmean的程序算了一下，发现算到这一步老是报错‘Index exceeds matrix dimensions.’。而且这程序后面部分看起来也觉得怪怪的，不明白真正用意，望解析。

function m = nanmean(x,dim)
nans = isnan(x);
x(nans) = 0;

if nargin == 1 % let sum deal with figuring out which dimension to use
    % Count up non-NaNs.
    n = sum(~nans);
    n(n==0) = NaN; % prevent divideByZero warnings
    % Sum up non-NaNs, and divide by the number of non-NaNs.
    m = sum(x) ./ n;
else
    % Count up non-NaNs.
    n = sum(~nans,dim);
    n(n==0) = NaN; % prevent divideByZero warnings
    % Sum up non-NaNs, and divide by the number of non-NaNs.
    m = sum(x,dim) ./ n;
end

二爷名声在外

简单的看了看，第一部分说的是只有一个参数的时候怎么计算；第二部分讲的是有两个参数的时候怎么计算……
区别就是第二段是由你限定维数了，你的错误也和维数有关，试着加上维数这个限定值试试……