<< Chapter < Page Chapter >> Page >
σ ˜ k ( v ˜ k ( t ) - v ˜ k + 1 ( t ) ) + σ ˜ k - 1 ( v ˜ k ( t ) - v ˜ k - 1 ( t ) ) + m ˜ k v ˜ k ' ' = 0 , k = 1 , 2 , , n 2

By assuming the displacement of bead m k behaves as u k e λ t , the equations of motion for a beaded string give recurrence relations for the amplitudes of each individual bead:

u k = R 2 k - 2 ( λ 2 ) u 1 , u ˜ k = R ˜ 2 k - 2 ( λ 2 ) u ˜ 1

Boyko and Pivovarchik show that the characteristic polynomial of the system of equations governing the motion of the beads is given by

φ ( λ 2 ) = R 2 n 1 ( λ 2 ) R ˜ 2 n 2 - 1 ( λ 2 ) + R 2 n 1 - 1 ( λ 2 ) R ˜ 2 n 2 ( λ 2 ) .

The zeros of φ ( λ 2 ) , λ k 2 { k = 1 , . . . , n 1 + n 2 }, are the eigenvalues of the system of ordinary differential equations that describe the motion of the string and are the frequencies at which the whole string vibrates. The polynomials R 2 n 1 ( λ 2 ) and R ˜ 2 n 2 ( λ 2 ) have roots that are the eigenvalues of the fixed-fixed boundary value problem for the left and right strings formed by clamping the string at the point where the left and right parts meet. The polynomials R 2 n 1 - 1 ( λ 2 ) and R ˜ 2 n 2 - 1 ( λ 2 ) have roots that are the corresponding eigenvalues of the fixed-flat boundary value problem. Although for the forward problem it is usually simpler to solve the eigenvalue problem generated directly from the ODE's given by equations and and the assumption that the displacements are complex exponentials, it is important to consider these R polynomials for the inverse problem.

The three-spectral inverse problem

Denote the spectra of the unclamped string by λ k , and the spectra of the left and right parts by ν k , and ν k , r . L is the length of the whole string and L and L r are the lengths of the separate parts. From this information we immediately construct three polynomials:

p w ( λ 2 ) = L k = 1 n 1 + n 2 1 - λ 2 λ k 2 , p ( λ 2 ) = L k = 1 n 1 1 - λ 2 ν k , 2 , p r ( λ 2 ) = L r k = 1 n 2 1 - λ 2 ν k , r 2

Note that p w ( λ 2 ) is proportional to φ ( λ 2 ) , p ( λ 2 ) is proportional to R 2 n 1 ( λ 2 ) , and p r is proportional to R ˜ 2 n 2 ( λ 2 ) . It is known that the ratio of polynomials R 2 n 1 ( z ) / R 2 n 1 - 1 ( z ) has the continued fraction expansion:

R 2 n 1 ( λ 2 ) R 2 n 1 - 1 ( λ 2 ) = n + 1 - m n σ λ 2 + 1 n - 1 + 1 - m n - 1 σ λ 2 + . . . + 1 1 + 1 - m 1 σ λ 2 + 1 0

which reveals the masses of the beads and the lengths between them that we are looking for. Since p ( λ 2 ) is proportional to R 2 n 1 ( λ 2 ) , we search for a second polynomial q ( λ 2 ) such that p ( λ 2 ) / q ( λ 2 ) gives the sought after continued fraction expansion. The same reasoning applies to the right part of the string.

q ( λ 2 ) = k = 1 n 1 λ 2 p w ( ν k , 2 ) ν k , 2 p r ( ν k , 2 ) j = 1 , j k n 1 ( λ 2 - ν j , 2 ) ( ν k , 2 - ν j , 2 ) + k = 1 n 1 ν k , 2 - λ 2 ν k , 2
q r ( λ 2 ) = k = 1 n 2 λ 2 p w ( ν k , r 2 ) ν k , r 2 p ( ν k , r 2 ) j = 1 , j k n 2 ( λ 2 - ν j , r 2 ) ( ν k , r 2 - ν j , r 2 ) + k = 1 n 2 ν k , r 2 - λ 2 ν k , r 2

After constructing these polynomials, we find their roots, which give us the second set of spectra for each of the left and right string parts. We now have all of the necessary information (fixed-fixed and fixed-flat spectra) to consider the problem reduced to two two-spectral problems. We use the algorithm presented in to recover the lengths and masses from the continued fraction expansion using only the roots of the polynomials.

Experimental setup

In the laboratory we have a beaded monochord consisting of steel piano wire held taut between two 5C collet fixtures. Behind one of the collets is a force transducer, and behind the other is a tensioner. The collets are initially open while the string is tensioned. The tension is monitored during this process via the force transducer. The collets are then closed and the system is ready for collecting data. Two phototransistors are positioned at each end of the string in perpendicular axes. In this way we can monitor the transverse vibrations both parallel and perpendicular to the ground at both ends of the string. The clamp is placed somewhere in the middle of the string such that at least one bead is left on each side of it. Both the tension and the phototransistors are monitored via a data acquisition card in a PC, and the information is collected and processed by MATLAB scripts.

Beaded monochord setup used to record data
The clamping mechanism
A two axis photodetector

Experimental results

Experimental testing has given very positive results. Despite the difficulties in collecting real data, the inverse algorithm just described has produced reasonable bead masses and lengths of string between the beads.

The peaks of the FFT plots shown in Figures , , and , give the observed frequencies for each section of the string. These are the eigenvalues that are used to construct the polynomials used during the inversion procedure.

Two sets of lenghts and masses that we recovered are depicted visually in figures and .

FFT plot for the whole string
FFT plot for the left clamped string
FFT plot for the right clamped string

Results for a string with 3 beads

Measured
Recovered
m 1 m ˜ 1 m ˜ 2 1 2 ˜ 1 ˜ 2 ˜ 3
Measured 30.8 30.8 17.8 0.203 0.245 0.184 0.203 0.289
Recovered 29.5 29.7 16.1 0.207 0.240 0.192 0.211 0.273

Results for a string with 4 beads

Measured
Recovered
m 1 m 2 m ˜ 1 m ˜ 2 1 2 3 ˜ 1 ˜ 2 ˜ 3
Measured 30.8 17.8 30.8 30.8 0.127 0.229 0.092 0.362 0.152 0.162
Recovered 27.7 17.2 29.6 29.0 0.139 0.226 0.082 0.368 0.150 0.159

Acknowledgements

This Connexions module describes work conducted as part of Rice University's VIGRE program, supported by National Science Foundationgrant DMS–0739420.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, The art of the pfug. OpenStax CNX. Jun 05, 2013 Download for free at http://cnx.org/content/col10523/1.34
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'The art of the pfug' conversation and receive update notifications?

Ask