octave_packages/control-2.3.52/Boeing707.m

   1 ## Copyright (C) 1997, 2000, 2004, 2005, 2006, 2007 Kai P. Mueller
   2 ##
   3 ##
   4 ## This program is free software; you can redistribute it and/or modify it
   5 ## under the terms of the GNU General Public License as published by
   6 ## the Free Software Foundation; either version 3 of the License, or (at
   7 ## your option) any later version.
   8 ##
   9 ## This program is distributed in the hope that it will be useful, but
  10 ## WITHOUT ANY WARRANTY; without even the implied warranty of
  11 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 ## General Public License for more details.
  13 ##
  14 ## You should have received a copy of the GNU General Public License
  15 ## along with this program; see the file COPYING.  If not, see
  16 ## <http://www.gnu.org/licenses/>.
  17
  18 ## -*- texinfo -*-
  19 ## @deftypefn {Function File} {@var{sys} =} Boeing707 ()
  20 ## Creates a linearized state-space model of a Boeing 707-321 aircraft
  21 ## at @var{v}=80 m/s
  22 ## @iftex
  23 ## @tex
  24 ## ($M = 0.26$, $G_{a0} = -3^{\\circ}$, ${\\alpha}_0 = 4^{\\circ}$, ${\\kappa}= 50^{\\circ}$).
  25 ## @end tex
  26 ## @end iftex
  27 ## @ifnottex
  28 ## (@var{M} = 0.26, @var{Ga0} = -3 deg, @var{alpha0} = 4 deg, @var{kappa} = 50 deg).
  29 ## @end ifnottex
  30 ##
  31 ## System inputs: (1) thrust and (2) elevator angle.
  32 ##
  33 ## System outputs:  (1) airspeed and (2) pitch angle.
  34 ##
  35 ## @strong{Reference}: R. Brockhaus: @cite{Flugregelung} (Flight
  36 ## Control), Springer, 1994.
  37 ## @end deftypefn
  38
  39 ## Author: Kai P. Mueller <mueller@ifr.ing.tu-bs.de>
  40 ## Created: September 28, 1997
  41
  42 function outsys = Boeing707 ()
  43
  44   if (nargin != 0)
  45     print_usage ();
  46   endif
  47
  48   a = [-0.46E-01,            0.10681415316, 0.0,   -0.17121680433;
  49        -0.1675901504661613, -0.515,         1.0,    0.6420630320636088E-02;
  50         0.1543104215347786, -0.547945,     -0.906, -0.1521689385990753E-02;
  51         0.0,                 0.0,           1.0,    0.0];
  52
  53   b = [0.1602300107479095,      0.2111848453E-02;
  54        0.8196877780963616E-02, -0.3025E-01;
  55        0.9173594317692437E-01, -0.75283075;
  56        0.0,                     0.0];
  57
  58   c = [1.0, 0.0, 0.0, 0.0;
  59        0.0, 0.0, 0.0, 1.0];
  60
  61   d = zeros (2, 2);
  62
  63   inam = {"thrust"; "rudder"};
  64   onam = {"speed"; "pitch"};
  65   ## snam = {"x1"; "x2"; "x3"; "x4"};
  66
  67   outsys = ss (a, b, c, d, "inname", inam, "outname", onam);
  68
  69 endfunction