{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Plotting Particle Displacement In Sound Fields\n", "\n", "[back to overview page](index.ipynb)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The following is based on this Mathematica notebook: http://library.wolfram.com/infocenter/MathSource/780/" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "particles = 1000\n", "length = 2 # along x-axis\n", "width = 0.2 # along y-axis\n", "wavelength = 0.3 * length\n", "amplitude = 0.1 * wavelength" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "t = 0 # TODO: animation from 0 to 1" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "x = np.random.uniform(0, length, particles)\n", "y = np.random.uniform(0, width, particles)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "x += amplitude * np.sin((2 * np.pi / wavelength) * x - 2 * np.pi * t)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", "\n" ], "text/plain": [ "