doc/info/ja/bode.texi

   1 @menu
   2 * Functions and Variables for bode::
   3 @end menu
   4
   5 @node Functions and Variables for bode,  , bode, bode
   6 @section Functions and Variables for bode
   7
   8
   9 @deffn {関数} bode_gain (@var{H}, @var{range}, ...@var{plot_opts}...)
  10 Bodeのゲイン線図を描く関数。
  11
  12 例 (1から7までは
  13 @example
  14 @url{http://www.swarthmore.edu/NatSci/echeeve1/Ref/Bode/BodeHow.html},
  15 @end example
  16 から、8はRon Crummettから):
  17 @example
  18 (%i1) load("bode")$
  19
  20 (%i2) H1 (s) := 100 * (1 + s) / ((s + 10) * (s + 100))$
  21
  22 (%i3) bode_gain (H1 (s), [w, 1/1000, 1000])$
  23
  24 (%i4) H2 (s) := 1 / (1 + s/omega0)$
  25
  26 (%i5) bode_gain (H2 (s), [w, 1/1000, 1000]), omega0 = 10$
  27
  28 (%i6) H3 (s) := 1 / (1 + s/omega0)^2$
  29
  30 (%i7) bode_gain (H3 (s), [w, 1/1000, 1000]), omega0 = 10$
  31
  32 (%i8) H4 (s) := 1 + s/omega0$
  33
  34 (%i9) bode_gain (H4 (s), [w, 1/1000, 1000]), omega0 = 10$
  35
  36 (%i10) H5 (s) := 1/s$
  37
  38 (%i11) bode_gain (H5 (s), [w, 1/1000, 1000])$
  39
  40 (%i12) H6 (s) := 1/((s/omega0)^2 + 2 * zeta * (s/omega0) + 1)$
  41
  42 (%i13) bode_gain (H6 (s), [w, 1/1000, 1000]),
  43                   omega0 = 10, zeta = 1/10$
  44
  45 (%i14) H7 (s) := (s/omega0)^2 + 2 * zeta * (s/omega0) + 1$
  46
  47 (%i15) bode_gain (H7 (s), [w, 1/1000, 1000]),
  48                   omega0 = 10, zeta = 1/10$
  49
  50 (%i16) H8 (s) := 0.5 / (0.0001 * s^3 + 0.002 * s^2 + 0.01 * s)$
  51
  52 (%i17) bode_gain (H8 (s), [w, 1/1000, 1000])$
  53 @end example
  54
  55 この関数を使うためには、最初に@code{load("bode")}を書いてください。
  56 @code{bode_phase}も参照してください。
  57
  58 @opencatbox
  59 @category{Plotting}
  60 @category{Package bode}
  61 @closecatbox
  62
  63 @end deffn
  64
  65
  66 @deffn {関数} bode_phase (@var{H}, @var{range}, ...@var{plot_opts}...)
  67 Bodeの位相線図を描く関数。
  68
  69 例 (1から7までは
  70 @example
  71 @url{http://www.swarthmore.edu/NatSci/echeeve1/Ref/Bode/BodeHow.html},
  72 @end example
  73 から、8は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 この関数を使うためには、最初に@code{load("bode")}を書いてください。
 119 @code{bode_gain}も参照してください。
 120
 121 @opencatbox
 122 @category{Plotting}
 123 @category{Package bode}
 124 @closecatbox
 125
 126 @end deffn
 127