octave_packages/control-2.3.52/@ss/__ctranspose__.m

   1 ## Copyright (C) 2012   Lukas F. Reichlin
   2 ##
   3 ## This file is part of LTI Syncope.
   4 ##
   5 ## LTI Syncope is free software: you can redistribute it and/or modify
   6 ## it under the terms of the GNU General Public License as published by
   7 ## the Free Software Foundation, either version 3 of the License, or
   8 ## (at your option) any later version.
   9 ##
  10 ## LTI Syncope is distributed in the hope that it will be useful,
  11 ## but WITHOUT ANY WARRANTY; without even the implied warranty of
  12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  13 ## GNU General Public License for more details.
  14 ##
  15 ## You should have received a copy of the GNU General Public License
  16 ## along with LTI Syncope.  If not, see <http://www.gnu.org/licenses/>.
  17
  18 ## -*- texinfo -*-
  19 ## Conjugate transpose of SS models.
  20
  21 ## Author: Lukas Reichlin <lukas.reichlin@gmail.com>
  22 ## Created: May 2012
  23 ## Version: 0.1
  24
  25 function sys = __ctranspose__ (sys, ct)
  26
  27   a = sys.a;
  28   b = sys.b;
  29   c = sys.c;
  30   d = sys.d;
  31   e = sys.e;
  32
  33   if (ct)       # continuous-time
  34     sys.a = -a.';
  35     sys.b = -c.';
  36     sys.c = b.';
  37     sys.d = d.';
  38     sys.e = e.';
  39     sys.stname = repmat ({""}, rows (a), 1);
  40   else          # discrete-time
  41     [n, m] = size (b);
  42     p = rows (c);
  43     if (isempty (e))
  44       e = eye (n);
  45     endif
  46     sys.a = blkdiag (e.', eye (p));
  47     sys.b = [zeros(n, p); -eye(p)];
  48     sys.c = [b.', zeros(m, p)];
  49     sys.d = d.';
  50     sys.e = [a.', c.'; zeros(p, n+p)];
  51     sys.stname = repmat ({""}, n+p, 1);
  52   endif
  53
  54 endfunction