[Octave-cvsupdate] SF.net SVN: octave:[10805] trunk/octave-forge/extra/control-devel

[<par...@us...>]

Revision: 10805
          http://octave.svn.sourceforge.net/octave/?rev=10805&view=rev
Author:   paramaniac
Date:     2012-08-03 15:38:45 +0000 (Fri, 03 Aug 2012)
Log Message:
-----------
control-devel: add filter method

Modified Paths:
--------------
    trunk/octave-forge/extra/control-devel/INDEX

Added Paths:
-----------
    trunk/octave-forge/extra/control-devel/inst/@iddata/filter.m

Modified: trunk/octave-forge/extra/control-devel/INDEX
===================================================================
--- trunk/octave-forge/extra/control-devel/INDEX	2012-08-03 15:19:41 UTC (rev 10804)
+++ trunk/octave-forge/extra/control-devel/INDEX	2012-08-03 15:38:45 UTC (rev 10805)
@@ -6,6 +6,7 @@
   @iddata/detrend
   @iddata/diff
   @iddata/fft
+  @iddata/filter
   @iddata/get
   @iddata/ifft
   @iddata/merge

Added: trunk/octave-forge/extra/control-devel/inst/@iddata/filter.m
===================================================================
--- trunk/octave-forge/extra/control-devel/inst/@iddata/filter.m	                        (rev 0)
+++ trunk/octave-forge/extra/control-devel/inst/@iddata/filter.m	2012-08-03 15:38:45 UTC (rev 10805)
@@ -0,0 +1,59 @@
+## Copyright (C) 2012   Lukas F. Reichlin
+##
+## This file is part of LTI Syncope.
+##
+## LTI Syncope 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.
+##
+## LTI Syncope 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 LTI Syncope.  If not, see <http://www.gnu.org/licenses/>.
+
+## -*- texinfo -*-
+## @deftypefn {Function File} {@var{dat} =} filter (@var{dat}, @var{sys})
+## @deftypefnx {Function File} {@var{dat} =} filter (@var{dat}, @var{b}, @var{a})
+## Detrend outputs and inputs of dataset @var{dat} by
+## removing the best fit of a polynomial of order @var{ord}.
+## If @var{ord} is not specified, default value 0 is taken.
+## This corresponds to removing a constant.
+## @end deftypefn
+
+## Author: Lukas Reichlin <luk...@gm...>
+## Created: April 2012
+## Version: 0.1
+
+function dat = filter (dat, b, a = [], si = [])
+
+  if (nargin < 2 || nargin > 4 || ! isa (dat, "iddata"))
+    print_usage ();
+  endif
+
+  if (isa (b, "lti") && issiso (b))
+    si = a;
+    if (isct (b))
+      tsam = dat.tsam;
+      b = c2d (b, tsam{1});         # does this make sense?
+    endif
+    [b, a] = filtdata (b, "vector");
+  endif
+
+  dat.y = cellfun (@(y) filter (b, a, y, si, 1), dat.y, "uniformoutput", false);
+  dat.u = cellfun (@(u) filter (b, a, u, si, 1), dat.u, "uniformoutput", false);
+
+endfunction
+
+
+%!shared DATD, Z
+%! DAT = iddata ({[(1:10).', (1:2:20).'], [(10:-1:1).', (20:-2:1).']}, {[(41:50).', (46:55).'], [(61:70).', (-66:-1:-75).']});
+%! DATD = detrend (DAT, "linear");
+%! Z = zeros (10, 2);
+%!assert (DATD.y{1}, Z, 1e-10);
+%!assert (DATD.y{2}, Z, 1e-10);
+%!assert (DATD.u{1}, Z, 1e-10);
+%!assert (DATD.u{2}, Z, 1e-10);

This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.