# Are there transfer function, state space, and bode plots?

I am looking for something like tf(), ss(), and bode() in control package of Octave, but I don't see them in sage. Are there any equivalents?

edit retag close merge delete

las bases teóricas-matematicas son la equivalencia que buscas

( 2014-08-20 14:29:42 -0600 )edit

I'm sorry, but I can read nothing other than English and Chinese...

( 2014-08-20 16:58:33 -0600 )edit

Sort by » oldest newest most voted

There seems to be something like that in numpy/scipy, I will check that out.

UPDATE: Now I use Sage/Sympy to derive the symbolic function, and convert the numerator and denominator to coefficients of polynomials, and create LTI systems in Scipy to use functions like LTI.bode(), signal.convolve(), etc. Also, I really can't read most of your answers.

more

pues no se si son repetitivas, se define G(s) una ves y luego a partir de esa funcion se obtienen las graficas. Una solución es crear una funcion sage y pasarle a G(s) como argumento, esta funcion sage retorna las graficas de bode y nyquist, o grafica polar, tambien puedes obteber la de Nichols.

more

Luis Ratia has made some Sage notebook worksheets available: Bode plots, and just recently posted about it on the sage-support and sage-notebook Google groups.

more

Thanks, but I know how to do that. I am just not sure if there is an automatic way to help me do the s=I*omega and 20*log(abs(H),10) etc.

( 2014-08-20 13:19:38 -0600 )edit

En Sage, la graficas son exactas y las unidades son en decibelios en el eje de salida y radianes en la grafica de magnitud, en la fase son grados en el eje y y radianes en el eje x, en las dos graficas el eje x es semilogaritmico. Yo comparo estas con los libros de control automatico y los resultados son totalmente compatibles.

Edit: translation. In Sage, the graphics are exact and the units are decibels on the output axis and radians for the magnitude; in the phase, degrees in the $y$-axis and radians in the $x$-axis; in both graphics the $x$-axis is semilogarithmic. I compare these with books on automatic control and the results are fully compatible.

more

( 2014-08-20 10:27:26 -0600 )edit

I know how to plot if I type in the function manually. I am just trying to find a way so that I don't need to type lots of 20*log(abs(H),10) or H.subs(s=I*omega) things, they are all repetitive.

( 2014-08-20 13:18:12 -0600 )edit