/usr/share/octave/packages/3.2/nan-2.4.4/nanfilter.m is in octave-nan 2.4.4-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 | function [Y,Z] = nanfilter(B,A,X,z);
% NANFILTER is able to filter data with missing values encoded as NaN.
%
% [Y,Z] = nanfilter(B,A,X [, Z]);
%
% If X contains no missing data, NANFILTER should behave like FILTER.
% NaN-values are handled gracefully.
%
% WARNING: missing values can introduce aliasing - causing unintended results.
% Moreover, the behavior of bandpass and highpass filters in case of missing values
% is not fully understood, and might contain some pitfalls.
%
% see also: FILTER, SUMSKIPNAN, NANFFT, NANCONV, NANFILTER1UC
% $Id$
% Copyright (C) 2005,2011 by Alois Schloegl <alois.schloegl@gmail.com>
% This function is part of the NaN-toolbox available at
% http://pub.ist.ac.at/~schloegl/matlab/NaN/ and
% http://octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/extra/NaN/inst/
% This program is free software; you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation; either version 3 of the License, or
% (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program; if not, write to the Free Software
% Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
warning('NANFILTER is experimental. For more details see HELP NANFILTER');
na = length(A);
nb = length(B);
if any(size(X)==1)
nc = 1;
else
nc = size(X,2);
end;
if nargin<4,
[t,Z.S] = filter(B,A,zeros(na+nb,nc));
[t,Z.N] = filter(B,A,zeros(na+nb,nc));
elseif isnumeric(z),
Z.S = z;
[t, Z.N] = filter(B, A, zeros(na+nb,nc));
elseif isstruct(z),
Z = z;
end;
NX = isnan(X);
X(NX) = 0;
[Y , Z.S] = filter(B, A, X, Z.S);
[NY, Z.N] = filter(B, A, ~NX, Z.N);
Y = (sum(B)/sum(A)) * Y./NY;
|