44 lines
1.2 KiB
Matlab
44 lines
1.2 KiB
Matlab
function y = smartmean(x,dim)
|
|
%SMARTMEAN Average or mean value ignoring NaN.
|
|
%
|
|
% Same as MEAN except that it returns the mean of the finite elements
|
|
% instead of propagating NaN and Inf
|
|
% Returns NaN if there are no finite elements
|
|
%
|
|
% For vectors, MEAN(X) is the mean value of the elements in X. For
|
|
% matrices, MEAN(X) is a row vector containing the mean value of
|
|
% each column. For N-D arrays, MEAN(X) is the mean value of the
|
|
% elements along the first non-singleton dimension of X.
|
|
%
|
|
% MEAN(X,DIM) takes the mean along the dimension DIM of X.
|
|
%
|
|
% Example: If X = [0 1 2
|
|
% 3 4 5]
|
|
%
|
|
% then mean(X,1) is [1.5 2.5 3.5] and mean(X,2) is [1
|
|
% 4]
|
|
%
|
|
% See also MEDIAN, STD, MIN, MAX, COV.
|
|
|
|
% Copyright (c) 1984-98 by The MathWorks, Inc.
|
|
% $Revision: 5.13 $ $Date: 1997/11/21 23:23:55 $
|
|
|
|
if nargin==1,
|
|
% Determine which dimension SUM will use
|
|
dim = min(find(size(x)~=1));
|
|
if isempty(dim), dim = 1; end
|
|
k=isfinite(x);
|
|
x(~k)=0;
|
|
warning off
|
|
y=sum(x)./sum(k,dim);
|
|
y(all(~isfinite(x)))=NaN;
|
|
warning on
|
|
else
|
|
k=isfinite(x);
|
|
x(~k)=0;
|
|
warning off
|
|
y=sum(x,dim)./sum(k,dim);
|
|
y(all(~isfinite(x), dim))=NaN;
|
|
warning on
|
|
end
|