Python 辛逆拉普拉斯变换isn';不行吗?
当我尝试反变换(s+0.2)/(s*(s+0.2)+1时,我对Symphy中的逆拉普拉斯变换有一个大问题 我正在使用的代码:Python 辛逆拉普拉斯变换isn';不行吗?,python,sympy,Python,Sympy,当我尝试反变换(s+0.2)/(s*(s+0.2)+1时,我对Symphy中的逆拉普拉斯变换有一个大问题 我正在使用的代码: from sympy import * s = symbols('s') t = symbols('t') f = (s + 0.2)/(s*(s + 0.2) + 1) inverse_laplace_transform(f,s,t) 该代码破坏了程序,我不明白为什么,因为对于其他函数,该方法不会崩溃。显然,mpmath有一个问题,它是一个用于任意精度算术的模块,或者
from sympy import *
s = symbols('s')
t = symbols('t')
f = (s + 0.2)/(s*(s + 0.2) + 1)
inverse_laplace_transform(f,s,t)
该代码破坏了程序,我不明白为什么,因为对于其他函数,该方法不会崩溃。显然,mpmath有一个问题,它是一个用于任意精度算术的模块,或者该模块与sympy的集成 您可以通过将0.2替换为rational来规避该错误。像这样的事情对我很有用:
from sympy import *
s = symbols('s')
t = symbols('t')
a = Rational(1, 5)
f = (s + a)/(s*(s + a) + 1)
inverse_laplace_transform(f,s,t)
这将返回:
sqrt(11)*(sin(3*sqrt(11)*t/10) + 3*sqrt(11)*cos(3*sqrt(11)*t/10))*exp(-t/10)*Heaviside(t)/33