三阶多名义分母有理函数的Laplace逆变换
三阶多名义分母有理函数的Laplace逆变换
提问于 2013-12-08 10:04:20
我试图用渐近法对形式为d/(A*s*2+ B*s + C)的广义有理函数做Laplace变换。
from sympy import *
from sympy import inverse_laplace_transform as ilt
from sympy.abc import s, t
var('A:D')
eq = D/(A*s**2 + B*s + C)
solution = ilt(eq, s, t)从数学分析中知道答案是：
-(A (e^(t (-1/2 sqrt(C^2-4 B)-C/2))-e^(t (1/2 sqrt(C^2-4 B)-C/2))))/sqrt(C^2-4 B)但是，交感不会产生一个解决方案，代码将卡在无限大的CPU进程中，而不会产生特定的增益。但是像这样把eq放在这里：
eq = B/((s - A)**2 + B**2)交感将产生如下形式的方程：
-I*(I*exp(t*im(B))*sin(t*(re(B) - im(A))) - exp(t*im(B))*cos(t*(re(B) - im(A)))
+ I*exp(-t*im(B))*sin(t*(re(B) + im(A))) + exp(-t*im(B))*cos(t*(re(B) + im(A))))
*exp(t*re(A))*Heaviside(t)/2这不是我所需要的。对于如何对这样一个人类可读的答案作出同情让步，有什么建议吗？
回答 1
Stack Overflow用户
回答已采纳
发布于 2013-12-10 06:42:33
看看追踪，这看起来就像一个合法的bug。你应该在SymPy问题跟踪器上报告。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/20452128
复制相关文章
相似问题
腾讯云开发者
