[SciPy-User] Incoherent results with signal.impulse
Sun Nov 22 13:31:00 CST 2009
> >/ Hello,
> />/ I have scipy 0.7.1, python 2.6, and when I do :
> />/ n = scipy.array()
> />/ d = scipy.array([0.01, 0.2, 1.0])
> />/ T, yout = scipy.signal.impulse((n,d))
> />/ it gives incoherent results for yout.
> />/ And that doesn't occurs with [1.0, 2.0, 1.0] denominator.
> />/ Is it a bug ?
> />/ I'm doing something wrong ?
> />/ Is anybody knows a solution ?
> scipy.signal.impulse assumes that the state matrix A is diagonalizable,
> so it does not give a correct result when A is defective. I would call
> that a bug. :)
> The attached file contains the function impulse_response() that uses a
> different method to compute the impulse response. If run as a script,
> the code at the bottom of the file plots impulse responses computed by
> impulse_response() and by scipy.signal.impulse() for your example, and
> for two other values of the leading coefficient of your denominator.
Thank you very much, it works fine now !
Actualy, if i'm right, the solution is to use lsim2 ?
More information about the SciPy-User