Plotting Particle Displacement In Sound Fields¶

In [173]:

%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np

In [161]:

particles = 1000
length = 2  # along x-axis
width = 0.2  # along y-axis
wavelength = 0.3 * length
amplitude = 0.1 * wavelength

In [162]:

t = 0  # TODO: animation from 0 to 1

In [163]:

x = np.random.uniform(0, length, particles)
y = np.random.uniform(0, width, particles)

In [164]:

x += amplitude * np.sin((2 * np.pi / wavelength) * x - 2 * np.pi * t)

In [172]:

fig, ax = plt.subplots(figsize=(12, 1))
ax.scatter(x, y)
ax.axis('equal')
ax.axis('tight');

In [ ]:

External Links¶

To the extent possible under law, Matthias Geier has waived all copyright and related or neighboring rights to this work.