I am trying to replicate this:
https://math.stackexchange.com/questions/2280196/inverse-fourier-transform-of-a-partial-fraction-decomposition
with sympy.
With Mathematica I get just the same behavior as the above theory suggest.
With sympy I it looks wrong. At t<0, I get the result close to the correct
one for t>0. (I was not able to include all the assumptions in sympy).
Here is the sympy code:
import sympy as sym
sym.init_printing()
ω = sym.symbols('omega', real=True, positive=True)
R, λ = sym.symbols('R, lambda', complex=True)
t = sym.symbols('t', real=True, positive=True)
α = R/(sym.I*ω-λ)+sym.conjugate(R)/(sym.I*ω-sym.conjugate(λ))
α
sym.inverse_fourier_transform(α, ω, -t)
and the Mathematica:
a = InverseFourierTransform[ R/(I omega - lambda) + Conjugate[R]/(I omega
- Conjugate[lambda]), omega, t,
FourierParameters -> {1, -1}]
Simplify[a, {Re[lambda] < 0, t > 0}]
Is the sympy result really wrong?
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at https://groups.google.com/group/sympy.
To view this discussion on the web visit
https://groups.google.com/d/msgid/sympy/9a4c2b02-c786-4aee-826f-0983b35f5132%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
