1 @c /bode.texi/1.1/Mon Feb 27 22:09:17 2006//
3 * Definições para bode::
6 @node Definições para bode, , bode, bode
7 @section Definições para bode
10 @deffn {Função} bode_gain (@var{H}, @var{range}, ...@var{plot_opts}...)
11 Função para desenhar gráficos de ganho para Bode.
13 Exemplos (1 a 7 provenientes de
15 http://www.swarthmore.edu/NatSci/echeeve1/Ref/Bode/BodeHow.html,
17 8 proveniente de Ron Crummett):
21 (%i2) H1 (s) := 100 * (1 + s) / ((s + 10) * (s + 100))$
23 (%i3) bode_gain (H1 (s), [w, 1/1000, 1000])$
25 (%i4) H2 (s) := 1 / (1 + s/omega0)$
27 (%i5) bode_gain (H2 (s), [w, 1/1000, 1000]), omega0 = 10$
29 (%i6) H3 (s) := 1 / (1 + s/omega0)^2$
31 (%i7) bode_gain (H3 (s), [w, 1/1000, 1000]), omega0 = 10$
33 (%i8) H4 (s) := 1 + s/omega0$
35 (%i9) bode_gain (H4 (s), [w, 1/1000, 1000]), omega0 = 10$
39 (%i11) bode_gain (H5 (s), [w, 1/1000, 1000])$
41 (%i12) H6 (s) := 1/((s/omega0)^2 + 2 * zeta * (s/omega0) + 1)$
43 (%i13) bode_gain (H6 (s), [w, 1/1000, 1000]),
44 omega0 = 10, zeta = 1/10$
46 (%i14) H7 (s) := (s/omega0)^2 + 2 * zeta * (s/omega0) + 1$
48 (%i15) bode_gain (H7 (s), [w, 1/1000, 1000]),
49 omega0 = 10, zeta = 1/10$
51 (%i16) H8 (s) := 0.5 / (0.0001 * s^3 + 0.002 * s^2 + 0.01 * s)$
53 (%i17) bode_gain (H8 (s), [w, 1/1000, 1000])$
56 Para usar essa função escreva primeiramente @code{load("bode")}. Veja também @code{bode_phase}
60 @deffn {Função} bode_phase (@var{H}, @var{range}, ...@var{plot_opts}...)
61 Função para desenhar gráficos de fase para Bode
63 Exemplos (1 a 7 provenientes de
65 http://www.swarthmore.edu/NatSci/echeeve1/Ref/Bode/BodeHow.html,
67 8 proveniente de Ron Crummett):
71 (%i2) H1 (s) := 100 * (1 + s) / ((s + 10) * (s + 100))$
73 (%i3) bode_phase (H1 (s), [w, 1/1000, 1000])$
75 (%i4) H2 (s) := 1 / (1 + s/omega0)$
77 (%i5) bode_phase (H2 (s), [w, 1/1000, 1000]), omega0 = 10$
79 (%i6) H3 (s) := 1 / (1 + s/omega0)^2$
81 (%i7) bode_phase (H3 (s), [w, 1/1000, 1000]), omega0 = 10$
83 (%i8) H4 (s) := 1 + s/omega0$
85 (%i9) bode_phase (H4 (s), [w, 1/1000, 1000]), omega0 = 10$
89 (%i11) bode_phase (H5 (s), [w, 1/1000, 1000])$
91 (%i12) H6 (s) := 1/((s/omega0)^2 + 2 * zeta * (s/omega0) + 1)$
93 (%i13) bode_phase (H6 (s), [w, 1/1000, 1000]),
94 omega0 = 10, zeta = 1/10$
96 (%i14) H7 (s) := (s/omega0)^2 + 2 * zeta * (s/omega0) + 1$
98 (%i15) bode_phase (H7 (s), [w, 1/1000, 1000]),
99 omega0 = 10, zeta = 1/10$
101 (%i16) H8 (s) := 0.5 / (0.0001 * s^3 + 0.002 * s^2 + 0.01 * s)$
103 (%i17) bode_phase (H8 (s), [w, 1/1000, 1000])$
105 (%i18) block ([bode_phase_unwrap : false],
106 bode_phase (H8 (s), [w, 1/1000, 1000]));
108 (%i19) block ([bode_phase_unwrap : true],
109 bode_phase (H8 (s), [w, 1/1000, 1000]));
112 Para usar essa função escreva primeiramente @code{load("bode")}. Veja também @code{bode_gain}