1

estoy teniendo problemas para obtener la transformada inversa de laplace de la siguiente funcion de transferencia wtl:

                            -2⋅s         
         (1.0e-8⋅s + 0.02)⋅ℯ             
────────────────────────────────────────
  ⎛         2                           ⎞
s⋅⎝1.0e-11⋅s  + 2.0000003e-5⋅s + 0.00177⎠

Estoy usado sympy, y he intentando con la siguiente expresion:

wtl_t = sp.inverse_laplace_transform(wtl, s, t, noconds=True)

Tambien probe con expandir por fracciones simples, o el denominador factorizado, pero no hay forma..

El error, a decir verdad, es bastante extenso, y creo que tiene algo que ver con la expresion a antitranformar:

Traceback (most recent call last):
  File "/usr/local/lib/python3.9/dist-packages/sympy/polys/polyutils.py", line 211, in _parallel_dict_from_expr_if_gens
    monom[indices[base]] = exp
KeyError: RisingFactorial(_t + 1, 1.0)

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/usr/local/lib/python3.9/dist-packages/sympy/polys/partfrac.py", line 81, in apart
    (P, Q), opt = parallel_poly_from_expr((P, Q), x, **options)
  File "/usr/local/lib/python3.9/dist-packages/sympy/polys/polytools.py", line 4340, in parallel_poly_from_expr
    return _parallel_poly_from_expr(exprs, opt)
  File "/usr/local/lib/python3.9/dist-packages/sympy/polys/polytools.py", line 4393, in _parallel_poly_from_expr
    reps, opt = _parallel_dict_from_expr(exprs, opt)
  File "/usr/local/lib/python3.9/dist-packages/sympy/polys/polyutils.py", line 332, in _parallel_dict_from_expr
    reps, gens = _parallel_dict_from_expr_if_gens(exprs, opt)
  File "/usr/local/lib/python3.9/dist-packages/sympy/polys/polyutils.py", line 216, in _parallel_dict_from_expr_if_gens
    raise PolynomialError("%s contains an element of "
sympy.polys.polyerrors.PolynomialError: RisingFactorial(_t + 1, 1.0) contains an element of the set of generators.

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/mnt/hgfs/FACULTAD/3er año/SDC/S&C_withPython/SDC_ex/GTP2_ex3.py", line 62, in <module>
    wtl_t = sp.inverse_laplace_transform(wtl, s, t, noconds=True)
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 1341, in inverse_laplace_transform
    return InverseLaplaceTransform(F, s, t, plane).doit(**hints)
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 130, in doit
    return self._compute_transform(self.function,
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 1287, in _compute_transform
    return _inverse_laplace_transform(F, s, t, self.fundamental_plane, **hints)
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 218, in wrapper
    res = func(*args, **kwargs)
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 1213, in _inverse_laplace_transform
    f, cond = inverse_mellin_transform(F, s, exp(-t), (None, oo),
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 908, in inverse_mellin_transform
    return InverseMellinTransform(F, s, x, strip[0], strip[1]).doit(**hints)
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 130, in doit
    return self._compute_transform(self.function,
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 846, in _compute_transform
    return _inverse_mellin_transform(F, s, x, strip, **hints)
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 218, in wrapper
    res = func(*args, **kwargs)
  File "/usr/local/lib/python3.9/dist-packages/sympy/integrals/transforms.py", line 773, in _inverse_mellin_transform
    h = hyperexpand(G)
  File "/usr/local/lib/python3.9/dist-packages/sympy/simplify/hyperexpand.py", line 2513, in hyperexpand
    return f.replace(hyper, do_replace).replace(meijerg, do_meijer)
  File "/usr/local/lib/python3.9/dist-packages/sympy/core/basic.py", line 1497, in replace
    rv = walk(self, rec_replace)
  File "/usr/local/lib/python3.9/dist-packages/sympy/core/basic.py", line 1481, in walk
    rv = F(rv)
  File "/usr/local/lib/python3.9/dist-packages/sympy/core/basic.py", line 1490, in rec_replace
    v = _value(expr, result)
  File "/usr/local/lib/python3.9/dist-packages/sympy/core/basic.py", line 1419, in <lambda>
    _value = lambda expr, result: value(*expr.args)
  File "/usr/local/lib/python3.9/dist-packages/sympy/simplify/hyperexpand.py", line 2509, in do_meijer
    r = _meijergexpand(G_Function(ap[0], ap[1], bq[0], bq[1]), z,
  File "/usr/local/lib/python3.9/dist-packages/sympy/simplify/hyperexpand.py", line 2394, in _meijergexpand
    slater2, cond2 = do_slater(tr(func.bm), tr(func.an), tr(func.bq), tr(func.ap),
  File "/usr/local/lib/python3.9/dist-packages/sympy/simplify/hyperexpand.py", line 2320, in do_slater
    hyp = _hyperexpand(Hyper_Function(nap, nbq), harg, ops,
  File "/usr/local/lib/python3.9/dist-packages/sympy/simplify/hyperexpand.py", line 2061, in _hyperexpand
    formula = try_lerchphi(func)
  File "/usr/local/lib/python3.9/dist-packages/sympy/simplify/hyperexpand.py", line 1799, in try_lerchphi
    part = apart(numer/denom, t)
  File "/usr/local/lib/python3.9/dist-packages/sympy/utilities/decorator.py", line 34, in threaded_func
    return func(expr, *args, **kwargs)
  File "/usr/local/lib/python3.9/dist-packages/sympy/polys/partfrac.py", line 84, in apart
    raise PolynomialError(msg)
sympy.polys.polyerrors.PolynomialError: RisingFactorial(_t + 1, 1.0) contains an element of the set of generators.

AGradezco cualquier sugenrencia.

2
  • 1
    Versión de sympy? Podría ser que estás chocando con este bug. Ahí mismo hay workarounds para hacer el cálculo que intentas...
    – Alfabravo
    el 30 ago. 2021 a las 17:29
  • 1
    Tengo sympy en la version 1.8 y python en la version 3.9.6. A ese foro habia entrado, pero tampoco encontre solucion, ni con evalf(). Ademas, defini s = sp.symbols('s') t = sp.symbols('t', positive=True) pero no se si influye en algo. Me parece algo tonto, porque es una FT sencilla.
    – Gino
    el 30 ago. 2021 a las 21:47

0

Tu Respuesta

Al pulsar en “Publica tu respuesta”, muestras tu consentimiento a nuestros términos de servicio, política de privacidad y política de cookies

Examina otras preguntas con la etiqueta o formula tu propia pregunta.