/usr/share/octave/packages/control-2.6.2/@lti/mconnect.m is in octave-control 2.6.2-1build1.
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 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 | ## Copyright (C) 2009-2014 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{sys} =} mconnect (@var{sys}, @var{m})
## @deftypefnx {Function File} {@var{sys} =} mconnect (@var{sys}, @var{m}, @var{inputs}, @var{outputs})
## Arbitrary interconnections between the inputs and outputs of an @acronym{LTI} model.
##
## @strong{Inputs}
## @table @var
## @item sys
## @acronym{LTI} system.
## @item m
## Connection matrix. Each row belongs to an input and each column represents an output.
## @item inputs
## Vector of indices of those inputs which are retained. If not specified, all inputs are kept.
## @item outputs
## Vector of indices of those outputs which are retained. If not specified, all outputs are kept.
## @end table
##
## @strong{Outputs}
## @table @var
## @item sys
## Interconnected system.
## @end table
##
## @strong{Example}
## @example
## @group
## Solve the system equations of
## y(t) = G e(t)
## e(t) = u(t) + M y(t)
## in order to build
## y(t) = H u(t)
## The matrix M for a (p-by-m) system G
## has m rows and p columns (m-by-p).
##
## Example for a 3x2 system:
## u1 = -1*y1 + 5*y2 + 0*y3
## u2 = pi*y1 + 0*y2 - 7*y3
##
## | -1 5 0 |
## M = | pi 0 7 |
## @end group
## @end example
## @end deftypefn
## Author: Lukas Reichlin <lukas.reichlin@gmail.com>
## Created: October 2009
## Version: 0.2
function sys = mconnect (sys, M, in_idx, out_idx = ":")
if (nargin < 2 || nargin > 4)
print_usage ();
endif
[p, m] = size (sys);
[mrows, mcols] = size (M);
if (p != mcols || m != mrows)
error ("mconnect: second argument must be a (%dx%d) matrix", m, p);
endif
if (! is_real_matrix (M))
error ("mconnect: second argument must be a matrix with real-valued coefficients");
endif
sys = __sys_connect__ (sys, M);
if (nargin > 2)
sys = __sys_prune__ (sys, out_idx, in_idx);
endif
endfunction
|