Numerical laplace in python

Laplace transform python

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.cm as cm

def s_plane_plot(sfunc, limits = [3,3,10], nsamp = 500):
    fig = plt.figure()
    ax = fig.gca(projection = '3d')

    sigma = np.linspace(-limits[0], limits[0], nsamp)
    omega = sigma.copy()

    sigma, omega = np.meshgrid(sigma, omega)

    s = sigma + 1j*omega 
    
    surf = ax.plot_surface(sigma, omega, np.abs(sfunc(s)), cmap = cm.flag)
    ax.set_zlim(0, limits[2])
    plt.xlabel('$\sigma$')
    plt.ylabel('$j\omega$')
    fig.tight_layout()

def X(s): return 1/((s + .2+.5j)*(s + .2-.5j))

s_plane_plot(X, limits = [1,1,4], nsamp =40)

Source: dsp.stackexchange.com

Laplace transform code in python

sympy.Heaviside(t-2)

Source: dynamics-and-control.readthedocs.io

Similar codes

0 vote277 views1 answer
Instead of: newlist = [] for i in range(1, 100): if i % 2 == 0: newlist.append(i**2)
Python
0 vote250 views1 answer
Django filter queryset by date
Python
3 votes260 views2 answers
Randomly shuffle array python
Python
0 vote248 views1 answer
Calculate volume of mask
Python
10 votes249 views3 answers
Python reduce
Python
1 vote311 views1 answer
Django queryset first element
Python Django
21 votes272 views2 answers
How to stop running code in python
Python
1 vote287 views2 answers
Python read parquet
Python
3 votes254 views6 answers
Squash commits git
Python
1 vote280 views2 answers
Program to print duplicates from a list of integers in python
Python