Drum PDEs, Double Fourier Expansions, and Synthesis

fourier analysisfourier seriesmusic-theorypartial differential equationswave equation

In studying the $2$D Wave Equation, the application often encountered is the displacement of a drum. The main solution to the PDE is a double Fourier summation, either a double Fourier Sine Series, a double Fourier Cosine Series, or a mixture of the two.

However, synthesizing a drum sound seems to typically involve passing the output of a single Fourier summation through various envelopes and filters. For example, this tutorial uses a square wave as a foundation, which is the type of mathematical object that would be built using a single Fourier summation (although the tutorial chooses to use a phasor, it could have equally well produced the square wave with a sinesum that allows direct input of numeric Fourier Sine Series coefficients).

Why the disparity? Since it is possible to use a single Fourier summation (in conjunction with secondary envelopes and filters, which can probably be boiled down to modifying the Fourier coefficients when all is said and done) to create a drum sound, is it also possible to reduce the double Fourier summation solution to the $2$D Wave Equation to a single Fourier summation? What would be the analogous interpretation of this dimensional disparity in the case of the Heat Equation?

Best Answer

Context

The wave equation takes the following form $$\frac{\partial^2 f}{\partial t^2} = c^2\Delta f$$ where $c>0$.

This partial differential equation is usually for solutions on a given domain $\Omega$ and with given boundary and initial conditions.

If your drum membrane at rest as being a flat surface, then you can model it as $\Omega\in\mathbb R^2$. In that case, a solution $f(x,y,t)$ of the PDE represents the vertical displacement of point $(x,y)$ of the membrane at time $t$.

Why the solution of the 2D wave equation involves a double summation

Remember that $f(x,y,t)$ is a function of two coordinates $x$ and $y$, because the membrane itself is two dimensional, and one of time $t$. So there is no reason to expect a single Fourier summation as the function itself seemingly has 3 degrees of freedom. If anything you should expect 3 sums. However, because the time behavior is related to the space behavior by the wave equation, you can show that 2 sums are enough. But the bottom line is, it's more than just one sum.

To convince yourself, imaging the simplest case of the membrane being a flat rectangle at rest. The solutions to the wave equation can be written as $f(x,y,t)=u(x)v(y)w(t)$ as the equation is separable. And using Fourier series, the solution you get is the product of two Fourier series, which is a double sum over the Fourier modes in the $x$ and $y$ dimensions. The time oscillations are related to those same modes and the velocity $c$ (that explains why the 3rd sum is not needed).

Why the synthesis of sound via a drum or a speaker involves a single summation

That's because you're talking about the excitation function, not the solution of the wave equation. And you excite the drum in only one location (e.g. at its center $(0,0)$). Thus the excitation is a function of only one variable (time). Consequently, if you represent it via Fourier series or the Fourier transform, you'll only need one sum.

More details here, feel free to skip: To produce a sound, you need to excite the drum or speaker membrane. That means you need to provide it with energy. For a drum, this is done by hitting it with a drumstick. Now, think about what a simple quick hit of the drum looks like mathematically. It's short burst of energy passed on to the drum at a given specific location (where the drumstick hits the drum, typically the center). In this case, we usually model this by saying that the initial value of the displacement $f(x,y,0)$ is the Dirac delta. But this is for a very short burst of energy at $t=0$.

Next, imagine you want to excite the drum by hitting it at its center once every second. Then you can model it as a train to Dirac deltas.

Finally, instead of a simple pattern of hitting the drum once per second, you want to to continuously excite the drum, at least do it for some period of time. But the way you want to do it is by passing various amounts of energy, through its center, at different times. Then, this is like replacing the Dirac delta above, or the Dirac delta train, with a function $\varphi(t)$. This is exactly what happens for a speaker (= drum + a mechanism to continuously excite the membrane). Hopefully, by now, you see that the square wave function they mention in the tutorial is, very likely (I haven't watched the Youtube video, too long), the excitation signal. And by definition, this is a one-variable function (of time $t$), and can be represented by one sum only.

Conclusion

The various sums you are considering refer to 2 different objects:

  • The waves on the drum membrane = function of 2 space variables and one of time, and needs double sums. You usually cannot reduce it to a single sum.
  • The excitation function = function of 1 variable (time) and needs only one sum.

Consequently, there is no reason to expect that the two are equal.

Related Question