Guitar string

A guitar string of length l l stretched along the x x -axis is plucked in the middle. The initial deflection of the string has the shape u ( x , t = 0 ) = u 0 ( x ) = { 2 A 0 x l x < l 2 2 A 0 ( l x ) l x l 2 , u(x, t = 0) = u_0(x) = \begin{cases} \frac{2 A_0 x}l & x < \frac l2 \\\\ \frac{2 A_0 (l - x)}l & x \geq \frac l2, \end{cases} where A 0 A_0 is in the middle (see diagram). The string is then released and swings freely for time t > 0 , t > 0, and the resulting string vibration can be written as a superposition of standing waves u ( x , t ) = n = 1 A n cos ( 2 π f n t ) sin ( n π l x ) , u(x, t) = \sum_{n = 1}^\infty A_n \cos(2 \pi f_n t) \sin\left( \frac{n\pi}{l} x \right), with eigenfrequencies f n = n ν f_n = n \cdot \nu and fundamental frequency ν \nu .

What is the relative amount of vibrational energy stored in the fundamental mode f = ν ? f = \nu? To do this, you'll have to determine the amplitudes A n A_n .

Hint:

  • The set of sine functions has the property 0 l sin ( n π l x ) sin ( m π l x ) d x = { l 2 n = m 0 n m , n , m N . \int_0^l \sin\left( \frac{n\pi}{l} x \right) \sin\left( \frac{m\pi}{l} x \right) dx = \begin{cases}\frac l2 & n = m \\ 0 & n \not= m, \quad n, m \in \mathbb{N}. \end{cases}
  • The kinetic energy of a small section of the string at position x x of length d x dx is given by 1 2 ( m l d x ) mass ( d u ( x , t ) d t ) 2 , \frac12 \overbrace{\left(\frac{m}{l}dx\right)}^\textrm{mass} \left(\frac{du(x,t)}{dt}\right)^2, where m m is the mass of the string. Assume that this holds along the entire string.


The answer is 0.81.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

2 solutions

Markus Michelmann
Dec 18, 2017

The amplitudes A n A_n (fourier coefficients) results A n = 2 0 l u 0 ( x ) sin ( n π l x ) d x = 4 A 0 l [ 0 l / 2 x sin ( n π l x ) d x + l / 2 l ( l x ) sin ( n π l x ) d x ] A 2 k = 0 , A 2 k + 1 = 8 A 0 l 0 l / 2 x sin ( ( 2 k + 1 ) π l x ) d x = 8 A 0 l ( l ( 2 k + 1 ) π ) 2 [ sin ( x ) x cos ( x ) ] 0 ( 2 k + 1 ) π / 2 = ± 8 A 0 l ( 2 k + 1 ) 2 π 2 \begin{aligned} A_n &= 2 \int_{0}^l u_0(x) \sin\left( \frac{n\pi}{l} x \right) dx \\ &= \frac{4 A_0}{l} \left[ \int_{0}^{l/2} x \sin\left( \frac{n\pi}{l} x \right) dx + \int_{l/2}^{l} (l - x) \sin\left( \frac{n\pi}{l} x \right) dx \right]\\ \Rightarrow \qquad A_{2k} &= 0\,, \\ A_{2k+1} &= \frac{8 A_0}{l} \int_{0}^{l/2} x \sin\left( \frac{(2k+1) \pi}{l} x \right) dx \\ &= \frac{8 A_0}{l} \left( \frac{l}{(2k + 1) \pi} \right)^2 \left[\sin(x) - x \cos(x) \right]_0^{(2k + 1)\pi/2} \\ &= \pm \dfrac{8 A_0 l}{(2k + 1)^2 \pi^2} \end{aligned} For the even multiples n = 2 k n = 2k the amplitude is zero, because due to the symmetry the two partial integrals cancel each other out. The kinetic energy of the vibration for a string of mass m m is E kin = 1 2 m 1 l 0 l ( u t ) 2 d x = π 2 m n = 1 ( f n A n ) 2 sin 2 ( 2 π f n t ) = m ( 8 A 0 l ν π ) 2 k = 0 sin 2 ( 2 π ( 2 k + 1 ) ν t ) ( 2 k + 1 ) 2 E kin = 1 2 m ( 8 A 0 l ν π ) 2 k = 0 1 ( 2 k + 1 ) 2 π 2 / 8 \begin{aligned} E_\text{kin} &= \frac{1}{2} m \cdot \frac{1}{l} \int_{0}^l \left(\frac{\partial u}{\partial t}\right)^2 dx \\ &= \pi^2 m \sum_{n = 1}^\infty (f_n A_n)^2 \sin^2(2 \pi f_n t) \\ &= m \left(\frac{8 A_0 l \nu}{\pi}\right)^2 \sum_{k = 0}^\infty \frac{\sin^2(2 \pi (2 k + 1) \nu t)}{(2 k + 1)^2} \\ \Rightarrow \quad \langle E_\text{kin} \rangle &= \frac{1}{2} m \left(\frac{8 A_0 l \nu}{\pi}\right)^2 \underbrace{\sum_{k = 0}^\infty \frac{1}{(2 k + 1)^2}}_{\pi^2/8} \end{aligned} Hereby, . \langle . \rangle is the time average. The relative energy contribution of the fundamental mode gives E kin , k = 0 E kin = 1 k = 1 1 ( 2 k + 1 ) 2 = 8 π 2 81 % \frac{\langle E_{\text{kin},k=0} \rangle}{\langle E_\text{kin} \rangle} = \frac{1}{\sum_{k = 1}^\infty \frac{1}{(2 k + 1)^2} } = \frac{8}{\pi^2} \approx 81\,\text{\%} Finally, a plot of the energy spectrum of the string vibration:

Max Yuen
May 3, 2019

Here's a quick and dirty solution:

We only need to know the energy carried in each "mode" for the triangle wave and use orthogonality to show that the total energy is just sum of the energies of the modes. (Alternatively, one can just use Parseval's Theorem, but I will show the details here.)

Let's make life easy and denote the time dependent normal modes of vibration as vectors Ψ n ( x ) e i ω n t = 2 l sin ( n π x l ) e i ω n t \Psi_n(x)e^{-i\omega_n t}=\sqrt{\frac{2}{l}}\sin{\left(\frac{n\pi x}{l}\right)}e^{-i\omega_n t} .

Since they are normal modes, they satisfy the orthogonal property 0 l Ψ n ( x ) Ψ m ( x ) d x = δ m n \int_0^l\Psi_n(x)\Psi_m(x)dx=\delta_{mn} , where δ m n \delta_{mn} is the Kronecker delta function.

From this it is easily seen that the kinetic energy for each mode with amplitude A n A_n is given by integrating the kinetic energy per unit mass over the entire length of the string with mass per unit length μ = m / l \mu=m/l :

Here, the integral we need to compute is the modulus of the complex amplitude:

K E n = 1 2 μ A n 2 0 l ( i ω n ) ( + i ω n ) Ψ n ( x ) 2 e i ω n t e + i ω n t d x = 1 2 μ A n 2 ω n 2 KE_n = \frac{1}{2}\mu A_n^2\int_0^l{(-i\omega_n)(+i\omega_n)\Psi_n(x)^2e^{-i\omega_n t}e^{+i\omega_n t}}dx=\frac{1}{2}\mu A_n^2\omega_n^2

Here's the neat trick, when we compute the energy of the superposition of modes that makes up the triangle wave, we can ignore the cross-terms since they integral to zero by the orthogonality of the modes. Only the self-terms survive. Also, note that ω n = n ω 1 \omega_n = n\omega_1 .

K E t o t a l = n = 1 1 2 μ n 2 A n 2 ω 1 2 KE_{total}=\sum_{n=1}{\frac{1}{2}\mu n^2A_n^2\omega_1^2}

Thus, once the common factors are cancelled out, the ratio of the fundamental energy to the total energy is just

η 1 = A 1 2 n = 1 n 2 A n 2 = ( n = 1 n 2 ( A n A 1 ) 2 ) 1 \eta_1=\frac{A_1^2}{\sum_{n=1}n^2A_n^2}=\left(\sum_{n=1}{n^2(\frac{A_n}{A_1})^2}\right)^{-1} .

Now all that is left is to obtain the A n A_n 's for the triangle wave.

It turns out that it is easier to obtain the Fourier series for a square wave with 50/50 duty cycle, and take the integral of the entire wave to obtain a triangle wave. The actual amplitudes don't matter, since we just need the ratio to the fundamental. Since all integrals of the sines cancel between 0 to l/2 if the sine wave is even, only the odd integrals survive. It isn't hard to convince yourself that the A n A_n 's are proportional to 1 n \frac{1}{n} where n n is odd.

Now, comes the neat trick, since the triangle wave is really just the integral of the square wave, we see that term by term integration of the wave will just get you a set of Fourier coefficients that goes like 1 n 2 \frac{1}{n^2} where n n is odd. (I'll leave this part for you to check.)

Thus, in computing the energy ratio, we have n 2 ( A n A 1 ) 2 = 1 / n 2 n^2(\frac{A_n}{A_1})^2=1/n^2 where n n is odd.

η 1 = ( n = 1 , o d d 1 n 2 ) 1 \eta_1 = \left(\sum_{n=1,odd}\frac{1}{n^2}\right)^{-1}

Now we compute the sum:

n = 1 1 n 2 = n = 1 , o d d 1 n 2 + n = 2 , e v e n 1 n 2 \sum_{n=1}\frac{1}{n^2}=\sum_{n=1,odd}\frac{1}{n^2}+\sum_{n=2,even}\frac{1}{n^2} since the series is absolutely convergent.

It's is not hard to recognize that the second sum is just 1/4 of the original sum. This means the odd sum is just 3/4 of the original sum.

The LHS sum here is a special series and is the Riemann Zeta function for z = 2 z=2 which sums to π 2 6 \frac{\pi^2}{6} .

Thus, our odd terms will add up to 3 4 π 2 6 = π 2 8 \frac{3}{4}\frac{\pi^2}{6}=\frac{\pi^2}{8} .

Finally, we plug this back to get the ratio we were looking for earlier:

η 1 = 8 π 2 0.81 \eta_1 = \frac{8}{\pi^2} \approx 0.81

Voila, 81% of the energy is in the fundamental.

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...