Mercurial > octave-nkf
view scripts/control/series.m @ 3279:6dd06d525de6
[project @ 1999-10-12 16:52:40 by jwe]
author | jwe |
---|---|
date | Tue, 12 Oct 1999 16:53:06 +0000 |
parents | 041ea33fbbf4 |
children | f7e4a95916f2 |
line wrap: on
line source
# Copyright (C) 1996 Auburn University. All Rights Reserved. # # 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 [a,b,c,d] = series(a1,b1,c1,d1,a2,b2,c2,d2) # Forms the series connection of two systems. # # Superseded by sysmult. Do not use this routine! # used internally in zp2ss # # Type of input: Transfer functions # Command: [num,den]=series(num1,den1,num2,den2) # Forms the series representation of the two transfer functions. # # Type of input: State space systems # Command: [a,b,c,d]=series(a1,b1,c1,d1,a2,b2,c2,d2) # Forms the series representation of the two state space system arguments. # The series connected system will have the inputs of system 1 and the # outputs of system 2. # # Type of input: system data structure # Command: syst=series(syst1,syst2) # Forms the series representation of the two mu system arguments. # Written by David Clem August 15, 1994 # If two arguments input, take care of mu system case warning("series is superseded by sysmult; use sysmult instead.") muflag = 0; if(nargin == 2) temp=b1; [a1,b1,c1,d1]=sys2ss(a1); [a2,b2,c2,d2]=sys2ss(temp); muflag = 1; endif # If four arguments input, put two transfer functions in series if(nargin == 4) a = conv(a1,c1); % was conv1 b = conv(b1,d1); % was conv1 c = 0; d = 0; # Find series combination of 2 state space systems elseif((nargin == 8)||(muflag == 1)) # check matrix dimensions [n1,m1,p1] = abcddim(a1,b1,c1,d1); [n2,m2,p2] = abcddim(a2,b2,c2,d2); if((n1 == -1) || (n2 == -1)) error("Incorrect matrix dimensions"); endif # check to make sure the number of outputs of system1 equals the number # of inputs of system2 if(p1 ~= m2) error("System 1 output / System 2 input connection sizes do not match"); endif # put the two state space systems in series a = [a1, zeros(rows(a1),columns(a2));b2*c1, a2]; b = [b1;b2*d1]; c = [d2*c1, c2]; d = [d2*d1]; # take care of mu output if(muflag == 1) a=ss2sys(a,b,c,d); b=c=d=0; endif endif endfunction