control.matlab.evalfr

control.matlab.evalfr(sys, x, squeeze=None)[source]

Evaluate transfer function of LTI system at complex frequency.

Returns the complex frequency response sys(x) where x is s for continuous-time systems and z for discrete-time systems, with m = sys.ninputs number of inputs and p = sys.noutputs number of outputs.

To evaluate at a frequency omega in radians per second, enter x = omega * 1j for continuous-time systems, or x = exp(1j * omega * dt) for discrete-time systems, or use freqresp(sys, omega).

Parameters:

sysStateSpace or TransferFunction: Linear system.
xcomplex scalar or 1D array_like: Complex frequency(s).
squeezebool, optional (default=True): If squeeze = True, remove single-dimensional entries from the shape of the output even if the system is not SISO. If squeeze = False, keep all indices (output, input and, if omega is array_like, frequency) even if the system is SISO. The default value can be set using config.defaults['control.squeeze_frequency_response'].

Returns:

frespcomplex ndarray: The frequency response of the system. If the system is SISO and squeeze is not True, the shape of the array matches the shape of omega. If the system is not SISO or squeeze is False, the first two dimensions of the array are indices for the output and input and the remaining dimensions match omega. If squeeze is True then single-dimensional axes are removed.