用Python推导傅里叶变换和拉普拉斯变换-CSDN博客

本文链接：https://blog.csdn.net/m0_37816922/article/details/147781662

文章目录

- Fourier变换
- 拉普拉斯变换

Fourier变换

【sympy】中封装了一系列积分变换的符号推导函数，其中，傅里叶变换和傅里叶变换的函数分别为

【fourier_transform】 $F(k)=\int^\infty_{-\infty}f(x)e^{-i2\pi ixk}\text dx$
【inverse_fourier_transform】 $f(x)=\int^\infty_{-\infty}F(x)e^{i2\pi ixk}\text dk$

则令 $f(x)=e^{-x^2}$ ，则 $F(k)=\sqrt{\pi} e^{- \pi^{2} k^{2}}$ ，变换过程如下。

from sympy import fourier_transform, inverse_fourier_transform
from sympy import exp, print_latex
from sympy.abc import x, k, a

F = fourier_transform(exp(-x**2), x, k)
print_latex(F)
f = inverse_fourier_transform(F)
print_latex(f)

拉普拉斯变换

Fourier变换存在两个问题，首先其作用对象为周期函数，其次要求这个函数不能发散。相比之下，拉普拉斯变换在变换过程中加入了一个衰减因子，从而在控制系统中更有优势。

【sympy】中封装了拉普拉斯变换及其逆变换的函数：

【laplace_transform】 $F(s)=\int^\infty_{0^-}e^{-st}f(t)\text dt$
【inverse_laplace_transform】 $f(t)=\frac{1}{2\pi i}\int^{c+i\infty}_{c-i\infty}$ e^{st}F(s)\text ds

其中，拉普拉斯变换函数的返回值是个元组，除了变换结果之外，还有结果的收敛于。例如，冲激信号 $\delta(t)$ 的拉普拉斯变换为常数，其返回值为 $\left( 1, \ -\infty, \ \text{True}\right)$ 。

from sympy import DiracDelta, laplace_transform
from sympy import inverse_laplace_transform, laplace_transform
from sympy.abc import t, s, a
F = laplace_transform(DiracDelta(t), t, s, simplify=True)
print_latex(F)
f = inverse_laplace_transform(F[0], s, t, simplify=True)
print_latex(f)