doc/info/bode.texi

   1 @menu
   2 * Functions and Variables for bode::
   3 @end menu
   4
   5 @node Functions and Variables for bode,  , Package bode, Package bode
   6 @section Functions and Variables for bode
   7
   8
   9 @anchor{bode_gain}
  10 @deffn {Function} bode_gain (@var{H}, @var{range}, ...@var{plot_opts}...)
  11 Function to draw Bode gain plots.
  12
  13 Examples (1 through 7 from
  14 @example
  15 @url{http://www.swarthmore.edu/NatSci/echeeve1/Ref/Bode/BodeHow.html},
  16 @end example
  17 8 from Ron Crummett):
  18 @example
  19 (%i1) load("bode")$
  20
  21 (%i2) H1 (s) := 100 * (1 + s) / ((s + 10) * (s + 100))$
  22
  23 (%i3) bode_gain (H1 (s), [w, 1/1000, 1000])$
  24
  25 (%i4) H2 (s) := 1 / (1 + s/omega0)$
  26
  27 (%i5) bode_gain (H2 (s), [w, 1/1000, 1000]), omega0 = 10$
  28
  29 (%i6) H3 (s) := 1 / (1 + s/omega0)^2$
  30
  31 (%i7) bode_gain (H3 (s), [w, 1/1000, 1000]), omega0 = 10$
  32
  33 (%i8) H4 (s) := 1 + s/omega0$
  34
  35 (%i9) bode_gain (H4 (s), [w, 1/1000, 1000]), omega0 = 10$
  36
  37 (%i10) H5 (s) := 1/s$
  38
  39 (%i11) bode_gain (H5 (s), [w, 1/1000, 1000])$
  40
  41 (%i12) H6 (s) := 1/((s/omega0)^2 + 2 * zeta * (s/omega0) + 1)$
  42
  43 (%i13) bode_gain (H6 (s), [w, 1/1000, 1000]),
  44                   omega0 = 10, zeta = 1/10$
  45
  46 (%i14) H7 (s) := (s/omega0)^2 + 2 * zeta * (s/omega0) + 1$
  47
  48 (%i15) bode_gain (H7 (s), [w, 1/1000, 1000]),
  49                   omega0 = 10, zeta = 1/10$
  50
  51 (%i16) H8 (s) := 0.5 / (0.0001 * s^3 + 0.002 * s^2 + 0.01 * s)$
  52
  53 (%i17) bode_gain (H8 (s), [w, 1/1000, 1000])$
  54 @end example
  55
  56 To use this function write first @code{load("bode")}. See also @mrefdot{bode_phase}
  57
  58 @opencatbox{Categories:}
  59 @category{Plotting}
  60 @category{Package bode}
  61 @closecatbox
  62
  63 @end deffn
  64
  65 @anchor{bode_phase}
  66 @deffn {Function} bode_phase (@var{H}, @var{range}, ...@var{plot_opts}...)
  67 Function to draw Bode phase plots.
  68
  69 Examples (1 through 7 from
  70 @example
  71 @url{http://www.swarthmore.edu/NatSci/echeeve1/Ref/Bode/BodeHow.html},
  72 @end example
  73 8 from Ron Crummett):
  74 @example
  75 (%i1) load("bode")$
  76
  77 (%i2) H1 (s) := 100 * (1 + s) / ((s + 10) * (s + 100))$
  78
  79 (%i3) bode_phase (H1 (s), [w, 1/1000, 1000])$
  80
  81 (%i4) H2 (s) := 1 / (1 + s/omega0)$
  82
  83 (%i5) bode_phase (H2 (s), [w, 1/1000, 1000]), omega0 = 10$
  84
  85 (%i6) H3 (s) := 1 / (1 + s/omega0)^2$
  86
  87 (%i7) bode_phase (H3 (s), [w, 1/1000, 1000]), omega0 = 10$
  88
  89 (%i8) H4 (s) := 1 + s/omega0$
  90
  91 (%i9) bode_phase (H4 (s), [w, 1/1000, 1000]), omega0 = 10$
  92
  93 (%i10) H5 (s) := 1/s$
  94
  95 (%i11) bode_phase (H5 (s), [w, 1/1000, 1000])$
  96
  97 (%i12) H6 (s) := 1/((s/omega0)^2 + 2 * zeta * (s/omega0) + 1)$
  98
  99 (%i13) bode_phase (H6 (s), [w, 1/1000, 1000]),
 100                    omega0 = 10, zeta = 1/10$
 101
 102 (%i14) H7 (s) := (s/omega0)^2 + 2 * zeta * (s/omega0) + 1$
 103
 104 (%i15) bode_phase (H7 (s), [w, 1/1000, 1000]),
 105                    omega0 = 10, zeta = 1/10$
 106
 107 (%i16) H8 (s) := 0.5 / (0.0001 * s^3 + 0.002 * s^2 + 0.01 * s)$
 108
 109 (%i17) bode_phase (H8 (s), [w, 1/1000, 1000])$
 110
 111 (%i18) block ([bode_phase_unwrap : false],
 112               bode_phase (H8 (s), [w, 1/1000, 1000]));
 113
 114 (%i19) block ([bode_phase_unwrap : true],
 115               bode_phase (H8 (s), [w, 1/1000, 1000]));
 116 @end example
 117
 118 To use this function write first @code{load("bode")}. See also @mrefdot{bode_gain}
 119
 120 @opencatbox{Categories:}
 121 @category{Plotting}
 122 @category{Package bode}
 123 @closecatbox
 124
 125 @end deffn
 126