Path: Common/Control
% Computes a numerical transfer function from the impulse response. This function uses an fft but effectively it is integrating exp(-st)*f(t) dt from zero to infinity. If the system is not damped, the results will be poor. For fastest results the number of samples should be a power of two. For a demo type: NumericalTransferFunction -------------------------------------------------------------------------- Form: [l, w] = NumericalTransferFunction( f, dT ) -------------------------------------------------------------------------- ------ Inputs ------ f (1,n) Impulse response dT (1,1) Time step ------- Outputs ------- l (3,n) Complex output w (1,n) Frequencies --------------------------------------------------------------------------