Mercurial > hg > octave-nkf
diff scripts/control/sys2tf.m @ 3213:ba1c7cdc6090
[project @ 1998-11-06 16:15:36 by jwe]
author | jwe |
---|---|
date | Fri, 06 Nov 1998 16:16:31 +0000 |
parents | |
children | dbcc24961c44 |
line wrap: on
line diff
new file mode 100644 --- /dev/null +++ b/scripts/control/sys2tf.m @@ -0,0 +1,58 @@ +# Copyright (C) 1996 A. Scottedward Hodel +# +# This file is part of Octave. +# +# Octave 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 2, or (at your option) any +# later version. +# +# Octave 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 Octave; see the file COPYING. If not, write to the Free +# Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. + +function [num,den,tsam,inname,outname] = sys2tf(Asys) +# function [num,den,tsam,inname,outname] = sys2tf(Asys) +# Conversion from a system data structure format to a transfer function. The +# transfer function part of ASYS is returned to the user in the form: +# +# num(s) +# G(s)=------- +# den(s) +# +# The user can also have the sampling time (TSAM), the name of the input +# (INNAME), and the output name (OUTNAME) + +# Written by R. Bruce Tenison (June 24, 1994) btenison@eng.auburn.edu +# modified to make sys2tf by A. S. Hodel Aug 1995 +# modified again for updated system format by John Ingram July 1996 +# $Revision: 1.2 $ + + if(nargin != 1) + usage("[num,den,tsam,inname,outname] = sys2tf(Asys)"); + endif + + if( !is_struct(Asys)) + error("Asys must be a system data structure (see ss2sys, tf2sys, zp2sys)"); + elseif (! is_siso(Asys) ) + [n, nz, m, p] = sysdimensions(Asys); + error(["system is not SISO (",num2str(m)," inputs, ... + ", num2str(p)," outputs"]); + endif + + Asys = sysupdate(Asys,"tf"); # just in case + + num = Asys.num; + den = Asys.den; + + tsam = Asys.tsam; + inname = Asys.inname; + outname = Asys.outname; + +endfunction +