/usr/share/octave/packages/control-2.6.2/@tf/display.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 90 91 92 93 94 95 96 97 98 99 100 | ## 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 -*-
## Display routine for TF objects.
## Author: Lukas Reichlin <lukas.reichlin@gmail.com>
## Created: September 2009
## Version: 0.4
function display (sys)
sysname = inputname (1);
[inname, outname, tsam] = __lti_data__ (sys.lti);
[inname, m] = __labels__ (inname, "u");
[outname, p] = __labels__ (outname, "y");
disp ("");
if (sys.inv && ! isct (sys))
[num, den] = filtdata (sys);
else
num = sys.num;
den = sys.den;
endif
for nu = 1 : m
disp (["Transfer function '", sysname, "' from input '", inname{nu}, "' to output ..."]);
disp ("");
for ny = 1 : p
__disp_frac__ (num{ny, nu}, den{ny, nu}, sys.tfvar, outname{ny});
endfor
endfor
display (sys.lti); # display sampling time
if (tsam == -2)
disp ("Static gain.");
elseif (tsam == 0)
disp ("Continuous-time model.");
else
disp ("Discrete-time model.");
endif
endfunction
function __disp_frac__ (num, den, tfvar, name)
MAX_LEN = 12; # max length of output name
tfp = isa (num, "tfpoly");
if (num == tfpoly (0))
str = [" ", name, ": 0"];
elseif ((tfp && den == 1) || (! tfp && isequal (den, 1)))
## elseif (den == tfpoly (1)) doesn't work because it
## would mistakingly accept non-tfpoly denominators like [0, 1]
str = [" ", name, ": "];
numstr = tfpoly2str (num, tfvar);
str = [str, numstr];
##elseif (length (den) == 1) # de-comment for non-development use
## str = [" ", name, ": "];
## num = num * (1/get (den));
## numstr = tfpoly2str (num, tfvar);
## str = [str, numstr];
else
numstr = tfpoly2str (num, tfvar);
denstr = tfpoly2str (den, tfvar);
fracstr = repmat ("-", 1, max (length (numstr), length (denstr)));
str = strjust (strvcat (numstr, fracstr, denstr), "center");
namestr = name(:, 1 : min (MAX_LEN, end));
namestr = [namestr, ": "];
namestr = strjust (strvcat (" ", namestr, " "), "left");
namestr = horzcat (repmat (" ", 3, 1), namestr);
str = [namestr, str];
endif
disp (str);
disp ("");
endfunction
|