Valores propios de cuerda pulsada/frecuencias armónicas: múltiplos enteros (o no)

Estoy tratando de derivar un modelo de una cuerda pulsada de la segunda ley de Newton. Mi derivación da como resultado

ω norte = C norte , norte = 1 , 2 , 3 entero
Creo que debería ser
ω norte = C norte , norte = 1 , 2 , 3 entero

Empecé con cuentas cada una con masa metro uniformemente espaciados en una cuerda con tensión de cuerda T . Talón norte se desplaza hacia arriba una distancia y . Las fuerzas en el talón norte me parece ser:

F norte = metro norte y ¨ norte = T pecado ( θ norte 1 , norte ) T pecado ( θ norte , norte + 1 )
pecado ( θ norte 1 , norte ) ( y norte y norte 1 ) / d
pecado ( θ norte , norte + 1 ) ( y norte y norte + 1 ) / d / d
Haciendo sustituciones y reorganizando:
y ¨ norte + T d metro norte ( y norte 1 + 2 y norte y norte + 1 ) = 0
Ahora deja
y norte = A norte mi i ω t , y ˙ norte = i ω y norte , y ¨ norte = w 2 y norte
también definir
C T d metro
Ahora
ω 2 + C 2 ( y norte 1 + 2 y norte y norte + 1 ) = 0
Expandiendo en forma matricial:
[ A ] y norte = ω 2 C 2 y norte
Dónde
| A | [ 2 1 0 0 0 0 1 2 1 0 0 0 0 1 2 1 0 0 0 0 1 2 1 0 0 0 0 1 2 1 0 0 0 0 1 2 ]
Una calculadora en línea muestra los valores propios de | A | ser: 1 , 2 , 3 , 4

Así que mi resultado parecería ser

ω norte = C norte , norte = 1 , 2 , 3 entero
¿Qué ocurre?

Relacionado si no duplicado: Ecuación de valor propio para energía cinética y potencial . La matriz tridiagonal es un caso especial de las llamadas matrices de Toeplitz [ver ecuación (19)] y sus valores propios y vectores propios podrían determinarse analíticamente [ver ecuaciones (20), (21) respectivamente].
¿Estás atado con cuerdas en los bordes a, digamos, paredes? ¿O hay infinitas cadenas?
norte ω 0 es correcto.
Parece que su modelo (por el cual se simula una cuerda de guitarra pulsada por norte cuentas idénticas distribuidas uniformemente a lo largo de la cuerda) es incorrecta.
Ustedes tienen razón. Utilicé un programa Schur (que había validado previamente) para calcular el sistema propio hasta 15 y los valores propios y no había números enteros allí (tenga en cuenta que había raíces cuadradas de números enteros). De acuerdo, mi modelo no funciona, pero la pregunta en mi mente aún permanece: ¿por qué no?

Respuestas (3)

¿Por qué pensarías ω norte = norte ω 0 o norte ω 0 ?

El "número cuántico" norte aquí hay impulso k ya que este es obviamente un sistema translacionalmente invariante. Con algunas transformaciones de Fourier, puedes obtener la energía ω norte = 2 ω 0 pecado ( k ) , dónde k = π / ( 2 norte ) i y i = norte , , norte , norte para el número total de sitios de celosía. Pequeñas desviaciones en k si el límite está abierto.

Esta dispersión concuerda con los valores propios de su matriz tridiagonal .

Pensé que mi resultado correspondería a la música real donde (pensé) los armónicos de un instrumento de cuerda pulsada son múltiplos enteros del tono fundamental.

Las respuestas publicadas por expertos anteriores dejan las dos preguntas fundamentales sin respuesta:

  1. ¿Por qué los valores propios de mi armónico uniforme, tridiagonal, simétrico (Toeplitz) no corresponden a las características conocidas de una cuerda de guitarra real?
  2. ¿Por qué los valores propios de matrices progresivamente más grandes no convergen a la solución continua del análisis lagrangiano?

Se puede encontrar ayuda en "Transporte de pulso casi sin dispersión en cadenas de masa de resorte cuasi uniformes largas ..." R. Via, 2 de agosto de 2018. Consulte las ecuaciones (38)-(41).

Los valores propios de la matriz ( A ) definido en mi pregunta (como se señaló en las respuestas anteriores y confirmado en el documento de referencia) son:

λ norte = 2 2 porque ( k norte )
dónde:
k norte ( norte π norte + 1 )
Entonces:
ω norte 2 C 2 = 2 2 porque ( k norte )

ω norte = C 2 2 porque ( k norte ) = 2 C pecado ( k norte 2 ) = 2 C pecado ( norte π 2 ( norte + 1 ) )
Mi primera pregunta restante se responde con una oración clave en el documento de referencia:

Para k 1 (es decir norte norte ) ellos [los ω norte ] son ​​casi lineales en k e igualmente espaciados en norte .

Para norte norte

pecado ( norte π 2 ( norte + 1 ) ) ( norte π 2 ( norte + 1 ) )

Así que en este caso particular:

ω norte ( C π norte + 1 ) norte

Pero como se indica en mi pregunta:

C T d metro
Manteniendo la longitud de la cuerda como una constante L y la masa total de la cuerda como constante METRO :
metro = METRO norte , d = METRO norte + 1
para muy grande norte , es decir, como la cadena de masas se convierte en una cuerda continua:
C T L METRO norte
Ahora podemos definir:
ω 0 T L METRO norte norte π = π T L METRO
Finalmente:
ω norte ω 0 norte , norte = 1 , 2 , 3 entero

que es la respuesta clásica y resuelve mi última pregunta pendiente.

Thomas: Veo que editaste mi respuesta. Gracias, supongo). Como una oportunidad de aprendizaje, me gustaría que explicara por qué cambió lo que hizo.

Los resultados de OP son incorrectos. Más precisamente si

(01) Ξ = [ 2 1 0 0 0 1 2 1 0 0 0 1 2 0 0 0 0 0 2 1 0 0 0 1 2 ]
es un norte × norte matriz de Toeplitz simétrica, sus valores propios no son los del OP
(02) λ ρ = ρ , ρ = 1 , 2 , 3 , norte (equivocado)
pero
(03) ξ ρ = 4 pecado 2 [ ρ π 2 ( norte + 1 ) ] = 2 ( 1 porque [ ρ π ( norte + 1 ) ] ) , ρ = 1 , 2 , , norte

= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =

Contraejemplo 1: Considere el 2 × 2 matriz de Toeplitz simétrica

(04) Ξ = [ 2 1 a b 1 2 a b ]
De acuerdo con los resultados de OP sus valores propios son
(05) λ 1 = 1 , λ 2 = 2 (equivocado)
en lugar de
(06) ξ 1 = 4 pecado 2 ( π 6 ) = 1 ξ 2 = 4 pecado 2 ( π 3 ) = 3
de acuerdo con la ecuación (03) .

= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =

Contraejemplo 2: Considere el 3 × 3 matriz de Toeplitz simétrica

(07) Ξ = [ 2 1 0 a b 1 2 1 a b 0 1 2 a b ]
De acuerdo con los resultados de OP sus valores propios son
(08) λ 1 = 1 , λ 2 = 2 , λ 3 = 3 (equivocado)
en lugar de
(09) ξ 1 = 4 pecado 2 ( π 8 ) = 2 [ 1 porque ( π 4 ) ] = 2 2 ξ 2 = 4 pecado 2 ( π 4 ) = 2 [ 1 porque ( π 2 ) ] = 2 ξ 3 = 4 pecado 2 ( 3 π 8 ) = 2 [ 1 porque ( 3 π 4 ) ] = 2 + 2
de acuerdo con la ecuación (03) .

Tenga en cuenta que

(10) λ 1 λ 2 λ 3 = 1 2 3 = 6 4 = det ( Ξ )
mientras
(11) ξ 1 ξ 2 ξ 3 = ( 2 2 ) 2 ( 2 + 2 ) = 4 = det ( Ξ )

Otro cheque es

(12) ρ = 1 ρ = norte λ ρ = ρ = 1 ρ = norte ρ = norte ( norte + 1 ) 2 2 norte = Rastro ( Ξ )

Aprecio tus cálculos, pero mi pregunta de física es, ¿por qué los cálculos no arrojan algo que corresponda a lo que ves en un analizador de frecuencia para una cuerda de guitarra pulsada?
@Nick Boshaft: Entonces, lo que llamas cuerda pulsada no es simulado por norte cuentas idénticas distribuidas uniformemente a lo largo de la cuerda. Por favor, aclare.
Frobenius: ¿Ves alguna falla en mi análisis de “n cuentas idénticas distribuidas uniformemente a lo largo de una cuerda tensada sin peso? De lo contrario, parece tener razón en que el modelo no modela correctamente una cuerda pulsada, ni con sus valores propios ni con los míos, incluso si usé más cuentas más pequeñas más juntas que se aproximan al continuo. ¿Por qué no?