efecto doppler de una fuente de sonido en un movimiento circular

Efecto Doppler para fuente de sonido en un movimiento circular- D.Tanoni y A. Krieger- UBA 2006 1

Efecto Doppler de una fuente de sonido en un movimiento circular: análisis con y sin reflexiones

Diego Tanoni y Alejandro Krieger e-mail: [email protected] y [email protected]

Laboratorio 5, Dpto. Física, FCEyN; UBA- 2006

En este trabajo se estudia el efecto Doppler de una fuente de sonido en un movimiento circular, para el caso libre de reflexiones y para el caso de reflexiones producidas en una pared cercana. Para el caso del movimiento libre, se estudió como varían distintos parámetros de la curva de frecuencia en función del tiempo variando la distancia mínima entre la fuente de sonido y el micrófono y la frecuencia angular de giro de la fuente. Para el caso del experimento en las cercanías de una pared, se trataron las reflexiones con el método de las imágenes. Se observa un buen acuerdo entre la teoría y el experimento y podemos decir que el método de las imágenes es adecuado para tratar las reflexiones de una onda sonora en una pared.

I. Introducción Cuando una fuente de sonido y un oyente están en movimiento relativo, la frecuencia del sonido oída por el oyente no es la misma que la frecuencia de la fuente f0. Este fenómeno, descrito por primera vez por el científico austriaco Christian Doppler, se llama efecto Doppler [1]. A continuación deduciremos una relación entre el cambio de frecuencia y las velocidades de la fuente y el oyente relativas al medio (usualmente aire). Consideremos el caso en que las velocidades de la fuente y el oyente están en la línea que los une y sea vf y voy las componentes de velocidad al lo largo de esta línea para la fuente y el oyente respectivamente, relativas al medio. Escogemos como dirección positiva la dirección que va desde la fuente al oyente. La rapidez del sonido relativa al medio, cs, siempre se considera positiva. Consideremos primero el caso de una fuente estacionaria y un oyente en movimiento (Figura 1.1). La fuente emite una onda sonora con frecuencia f0 . Las crestas de esta onda están separadas por una distancia λ igual a la longitud de onda λ = cs / f0 . Si un oyente se mueve con velocidad (-voy) hacia la fuente estacionaria, las crestas que se acercan al oyente en movimiento tienen una velocidad de propagación relativa al oyente de (cs – voy), así que la frecuencia con la que estas crestas llegan al oyente es

λoys

oy

vcf

−= (1.1)

OyenteFuente

λ

v oyente

v fuente

a)

b)

λ

Figura 1.1: La frecuencia que escucha el oyente es distinta si a) el oyente se acerca a la fuente o b) la fuente se acerca al oyente. La dirección que va de la fuente al oyente se toma como positiva.

Así, un oyente que se mueve hacia la fuente (voy < 0) escuchará una frecuencia más alta que un oyente estacionario y un oyente que se aleja (voy > 0) escuchará una frecuencia más baja.


Supongamos ahora que la fuente también se mueve con velocidad vf (Figura 1.1 b). La velocidad de la onda sonora relativa al medio (aire) sigue siendo cs . Está determinada por las propiedades del medio y no cambia por el movimiento de la fuente. Pero la longitud de onda ya no es igual a cs / f0.. El tiempo de emisión entre un máximo y el siguiente es el período de la onda T = 1/f0 . Durante este tiempo la onda viaja una distancia csT = cs / f0 y la fuente se mueve una distancia vf T = vf / f0 . Como la longitud de onda es la distancia entre dos máximos sucesivos, la longitud de onda para una fuente en movimiento depende del desplazamiento relativo de la fuente y la onda y es

000 fvc

fv

fc fsfs −

=−=λ (1.2)

Vemos que la longitud de onda no es la misma por delante y por detrás de una fuente móvil. Si la fuente se acerca al oyente, (vf>0) la longitud de onda es menor que en el caso de la fuente estacionaria, las crestas se juntan por delante de la fuente. Si la fuente se aleja del oyente (vf<0) la longitud de onda es mayor. Las ondas por delante y por detrás de la fuente se comprimen y estiran, respectivamente, por el movimiento de la fuente. Si sustituimos esta ecuación en (1.1), obtenemos una expresión que incluye todas las posibilidades de movimiento de la fuente y del oyente, relativas al medio, a lo largo de la línea que los une

0fvcvc

ffs

oysoy

−−

= (1.3)

Si la fuente se acerca al oyente, vf > 0 y la frecuencia foy es mayor que f0. Si el oyente se aleja de la fuente, voy >0 y la frecuencia que percibe el oyente es menor que f0. Si tanto la fuente como el oyente están en reposo o tienen la misma velocidad relativa respecto del medio foy = f0 . Para el caso de la fuente en reposo, vf = 0 y recuperamos el caso de la ecuación (1.1). Si el oyente esta en reposo respecto del medio, voy = 0, y podemos escribir

00

1

1 f

cv

fvc

cf

s

ffs

soy

−=

−= (1.4)

Si la velocidad de la fuente vf es mucho menor que la velocidad del sonido cs, podemos aproximar la ecuación (1.4) y escribir

01 fcv

fs

foy

+≈ (1.5)


II. Efecto Doppler para una fuente en movimiento circular a) Caso sin reflexiones Si tenemos una fuente de sonido en el extremo de una barra que gira con velocidad angular uniforme ω , la velocidad de la fuente que debemos considerar para analizar el efecto Doppler es la componente de vf paralela a la dirección fuente-oyente que designaremos v// . De la figura 2.1 pude verse que esta componente es

)(// βω +−= tsenvv f (2.1)

donde vf = ω r es la magnitud de la velocidad tangencial. El signo menos corresponde al hecho de que la fuente se esta alejando del observador.

βωt

D

lr

v//

vf = ω r βωt

Figura 2.1: Diagrama de una fuente de sonido en un movimiento circular. La velocidad que debemos considerar para el efecto Doppler es la componente v// de la velocidad de la fuente vf

Como se muestra en el Apéndice A, la ecuación (2.1) puede rescribirse como

lDtsenvv f )(// ω−= (2.2)

donde d es la distancia del micrófono al centro de giro y l es la distancia desde el micrófono a la fuente l 2 = D 2 + r 2 – 2 rD cos(ωt). Si llamamos α a la relacion α = r / D podemos escribir la ecuación (2.3) como

)cos(211)(// t

tvsenvωαα

ω−+

−= (2.4)

De la ecuación (2.1) se ve que a medida que la fuente se mueve en la trayectoria circular, v// alcanza su valor máximo y mínimo cuando (ω t + β ) = ± π / 2. En ese momento r y l son perpendiculares y entonces

αω ==Drt)cos( (2.5)

Esto significa que v// alcanza sus valores máximo y mínimo en los puntos donde la dirección fuente – micrófono es tangente a la circunferencia (Figura 2.2)

βωt

v// min

v// max

Figura 2.2: La velocidad v// alcanza sus valores máximos y mínimos en los puntos donde la dirección fuente-micrófono es tangente a la circunferencia.


Para encontrar cuál es el corrimiento de la frecuencia percibida por el oyente, podemos usar la aproximación dada por la ecuación (1.5), ya que en este caso la velocidad de giro vf es mucho menor que la velocidad del sonido. Tenemos entonces

0)cos(211)(1 f

ttsen

cv

fs

foy

−+−≈

ωααω (2.6)

La frecuencia máxima percibida por el oyente se produce en el momento en que v// es máxima. Esta frecuencia es f max = (1 + ω r / cs ) f0 y la frecuencia mínima es f min = (1 - ω r / cs ) f0 , por lo que la máxima diferencia de frecuencia es

0minmax 2 fc

rfffs

ω=∆=− (2.7)

También podemos calcular el tiempo δ que transcurre para que la velocidad pase de su valor máximo a su valor mínimo. Usando la ecuación (2.5) vemos que

tar =ω

α )cos( (2.8)

Este es el tiempo que transcurre desde t = 0 hasta el momento que la velocidad es mínima, por lo que el tiempo δ buscado es dos veces ese valor

δω

α =)cos(2ar (2.9)

Cuanto más lejos se encuentre el micrófono, más tiempo transcurrirá entre el máximo y el mínimo de frecuencia. Como caso extremo, si el oyente es encuentra en el infinito, este tiempo será de π/ω, la mitad del periodo T. Por otro lado si el micrófono se encuentra sobre la trayectoria de la fuente, la diferencia de tiempo tiende a cero y se produce un salto instantáneo de la frecuencia, ya que la fuente se acerca con velocidad máxima y un instante después se aleja con la velocidad máxima (Figura 2.3) Estas mismas conclusiones se obtienen analizando el problema desde el punto de vista de los puntos tangentes. De la Figura 2.3 se ve que si d → ∞ la dirección fuente-oyente es tangente en los bordes superior e inferior de la trayectoria, mientras que si d → r ,

d → ∞ d → r

v// min

v// max

v// min

v// max

Figura 2.3: El tiempo δ entre la velocidad máxima y mínima depende de la distancia del micrófono y es igual a la mitad del periodo si el oyente se encuentra en el infinito y tiende a cero si el oyente se encuentra sobre la trayectoria de la fuente.


hay un solo punto tangente, lo que significa que la frecuencia se vuelve máxima e inmediatamente mínima. b) Análisis considerando reflexiones Si el experimento se lleva a cabo cerca de alguna pared, aparecerán junto con la señal de la fuente las reflexiones que se producen en esa pared. Es posible tratar estas reflexiones con el método de las imágenes, considerando que las reflexiones son las señales de una fuente imaginaria ubicada al otro lado de la pared, girando en sentido contrario a la fuente real [2] (Figura 2.4) En ese caso se trata de una fuente virtual de la misma frecuencia f0 que la fuente real, que gira en sentido contrario al de la fuente real y cuyo centro de giro se encuentra más lejos del micrófono que el centro de la fuente real. Esta distancia D’ es D’ = 2 s + D, por lo que el valor de α’ para la fuente virtual es

Dsr+

=2

'α (2.10)

pared

Fuentereal

Fuenteimaginaria

ss

D

D´

Figura 2.4: Es posible tratar las reflexiones de una pared con el método de las imágenes. El diagrama muestra la fuente real y una imaginaria ubicada del otro lado de la pared, girando en sentido contrario y ubicada más lejos del micrófono.

Por esta razón esperamos una señal de la fuente imaginaria que sea de menor intensidad y opuesta a la señal de la fuente real (ya que su sentido de giro es opuesto). Observando la figura 2.4 y la ecuación (2.9) también esperamos para la señal reflejada un δ’ mayor, ya que la distancia al micrófono es mayor. Si el modelo de las imágenes es el adecuado para tratar las reflexiones, las mediciones de frecuencia en función del tiempo deberían seguir la relación

0)cos('2'11)(1 f

ttsen

cv

fs

foy

−++≈

ωααω (2.11)

Esta es la misma ecuación que la ecuación (2.6), sólo que el signo de la velocidad de la fuente ha cambiado y ahora usamos α’ en lugar de α.


III. El experimento El dispositivo experimental utilizado para esta experiencia consistió en una barra de (150 ± 1) cm de largo, en cuyo extremo se colocó una fuente de sonido (“chicharra”) alimentada con una batería de 9V (Figura 3.1) El centro de la barra se colocó sobre el eje de una motor, cuya frecuencia de giro podía modificarse manualmente. La señal de la fuente se capturó utilizando un micrófono comercial conectado a la placa de sonido de una PC y se adquirió usando el programa de audio SoundForge

Figura 3.1 : Dispositivo experimental utilizado para estudiar el efecto Doppler de una fuente de sonido en un movimiento circular.

Con este programa podía adquirirse la amplitud de la señal del micrófono en función del tiempo y ver el espectro de frecuencias de la señal. El programa realizaba además un sonograma de la señal, que es un diagrama de la frecuencia en función del tiempo. Realizamos varias mediciones modificando la distancia d (“distancia mínima”) para una frecuencia de giro constante. También hicimos mediciones para una distancia d fija modificando la frecuencia de giro. Finalmente acercamos todo el dispositivo a una pared e hicimos mediciones en esas condiciones. IV. Resultados y Análisis A continuación se muestra una parte del espectro de la señal de la fuente de sonido. Aquí se aprecia sólo uno de los picos, a la frecuencia de 5,3 KHz. El espectro completo presenta varios picos de intensidad en distintas frecuencias. Sin embargo, nuestro interés en esta práctica esta centrado sólo en el pico mencionado, que es el de mayor intensidad, por lo que para el resto del trabajo se considerará que la frecuencia de la señal fue de 5,3 KHz. Cómo se evitan los demás picos de intensidad se comenta en la siguiente sección.

Figura 4.1: Espectro de la señal sonora en el rango de frecuencias de interés.

Fuente de sonido

micrófono

bateria

motor de frecuenciavariable

dr = 75 cm

D


Sonogramas8 Como se ha visto, en el presente trabajo surgió la necesidad de determinar las frecuencias en función del tiempo, pero primero hay que entender que significa la frase frecuencia en función del tiempo. Como es sabido, en la naturaleza no existen ondas perfectamente monocromáticas (el problema ya surge teóricamente, debido a que una onda monocromática tendría energía infinita al no ser una función de cuadrado integrable), sino que se trata de los denominados trenes de onda con una determinada longitud espacial. Lo que se emplea en general para tratar estos temas es la transformada de Fourier, que es una función que depende de la frecuencia y esta directamente relacionada con la señal original (cada una se deduce de la otra). El problema que surge es que la transformada de Fourier implica el cálculo de una integral que va de tiempo menos infinito hasta tiempo más infinito y que contiene a la señal en todo este rango. Para elaborar el sonograma,8 el programa calcula la transformada de Fourier en varios intervalos de tiempos finitos. Estos intervalos deben ser lo suficientemente grandes como para ser considerados como infinitos desde el punto de vista de la señal (una apreciación es que cuanto mayor sea el intervalo de tiempo de integración mejor será la resolución en frecuencia), y lo suficientemente pequeños como para que la fuente de donde se emite la señal no cambie apreciablemente de lugar. En nuestro experimento, este compromiso puede ser satisfecho en ambos extremos si la frecuencia de giro de la fuente es mucho menor que la frecuencia de la señal sonora, hecho que se verifica ampliamente debido a que la fuente giró con frecuencias inferiores a los 20 Hz, mientras que la señal sonora estudiada fue de una frecuencia superior a los 5 KHz. Los sonogramas son un mapa que contiene en un eje al tiempo y en otro a la frecuencia. En dicho mapa una escala de colores indica la intensidad de cada frecuencia en la transformada de Fourier en cada intervalo de tiempo. La escala de colores con las intensidades correspondientes esta dada en la parte inferior del sonograma. Debido a que la frecuencia sonora es alta, los intervalos de tiempo que pueden tomarse para satisfacer ambas relaciones de compromiso son lo suficientemente pequeños como para que en el sonograma la escala de tiempo de observe como si se tratase de valores continuos. La siguiente imagen (Figura 4.2) muestra un sonograma en un intervalo de frecuencias más amplio del que se estudio (0 a 10 Khz). Este sonograma correspondió a una distancia mínima entre micrófono y señal de 3 cm, un radio de giro de la fuente de 75 cm y una frecuencia angular de giro de 13,5 Hz.

Figura 4.2: Sonograma mostrando un amplio rango de frecuencias, más allá de las de interés.


Como se aprecia en la figura anterior, hay varias bandas de alta intensidad. Cada una de ellas corresponde a un pico distinto en el espectro de la señal en reposo. De esta manera, al realizar las mediciones correspondientes a la frecuencia deseada, sólo se tomó en cuenta un rango de frecuencias alrededor de los 5,3 KHz, que es la frecuencia pico mostrada en la sección anterior. Los demás sonogramas que se muestren sólo estarán en dichos intervalos. Además, en las escalas de tiempo, frecuencia, e intensidad sólo se mostrará el sentido en que las magnitudes aumentan. Estos sonogramas se emplearon para obtener la frecuencia en función del tiempo, ya que los máximos de intensidad indican como varia con el tiempo la señal correspondiente al pico de 5,3 KHz en reposo. Resumiendo, cuando se hable de frecuencia en función del tiempo, lo que se esta indicando son los valores de frecuencia de máxima intensidad para cada tiempo; además, en todas las expresiones empleadas se tomará como frecuencia sonora a 5,3 KHz. La siguiente figura muestra un sonograma cuando el sistema está en reposo. Observar que, como es de esperar, el sonograma no varía con el transcurso del tiempo. La línea de mayor intensidad corresponde a las frecuencias cercanas a los 5,3 KHz.

Figura 4.3: Sonograma correspondiente a la fuente en reposo.

a) Análisis de los resultados al variar la distancia d Según el comentario hecho para la figura 2.3, al incrementar la distancia mínima d entre el micrófono y la fuente (de las figuras 2.1 y 3.1 se ve que D = r + d ), el intervalo de tiempo entre una frecuencia máxima y otra mínima aumenta debido a que es mayor el camino que recorre la fuente entre estos dos puntos (considerando fija la frecuencia de giro), resultando que cuando la distancia tiende a infinito, el intervalo de tiempo entre estos dos puntos tiende la mitad del periodo de giro de la fuente. Los siguientes sonogramas muestran este comportamiento, siendo en todos los caso la frecuencia angular de giro igual a 13,5 Hz y el radio de giro de la fuente de 75 cm.

Figura 4.4a Figura 4.4b Figura 4.4c Figuras 4.4a, 4.4b, 4.4c: Sonogrmas correspondientes a distancias mínima entre la fuente y el micrófono de 4 cm, 21 cm y 81 cm respectivamente.


Los gráficos, de izquierda a derecha, tienen distancia mínima d de 4 cm, 21 cm y 81 cm respectivamente. Se aprecia que a medida que aumenta la distancia d el comportamiento predicho anteriormente se verifica. En la última figura, el intervalo de tiempo entre una frecuencia máxima y otra mínima (lo llamaremos, de aquí en más, intervalo de tiempo entre frecuencia máxima – mínima δ ) es comparable al intervalo entre una mínima y otra máxima ( lo llamaremos, de aquí en más, intervalo de tiempo entre frecuencia mínima – máxima). La suma de ambos intervalos da un periodo de giro. A continuación se mostrará un gráfico de las mediciones del intervalo de tiempo entre frecuencia máxima – mínima δ en función de la distancia D, comparando los resultados obtenidos con los predichos teóricamente. Los datos experimentales se obtuvieron como sigue: se registraron cuatro frecuencias máximas con las frecuencias mínimas que les seguían inmediatamente, y se tomo las diferencias de tiempo entre ellas (el tiempo entre cada frecuencia máxima con el tiempo de la frecuencia mínima que le seguía inmediatamente), obteniéndose así cuatro intervalos de tiempo. El promedio entre estos cuatro valores fue el intervalo de tiempo que se colocó en el gráfico, siendo la desviación estándar el error. La curva teórica se realizó utilizando la expresión (2.9), siendo la frecuencia angular de giro de 13,5 Hz y el radio de giro de la fuente de 75 cm. Se muestra en el gráfico el valor de χ2, que mide la calidad del ajuste de los datos experimentales y las predicciones teóricas (ver Apéndice B, ecuación B.5). Para hacer las mediciones sobre el sonograma, en ésta y en las demás secciones donde fue necesario, se empleó el programa de computación xyExtract,7 el cual permite obtener coordenadas de puntos de una imagen luego de fijar referencias en cada uno de los ejes.

1 10 100

0.00

0.05

0.10

0.15

0.20

0.25 χ2 = 0.04485

Inte

rval

o de

Tie

mpo

(s)

Distancia (cm)

Resultados obtenidos Curva Teórica

Figura 4.5: Intervalo de tiempo δ entre frecuencia máxima – mínima en función de la distancia D. Comparación entre resultados experimentales y predicción teórica dada por la ecuación (2.9)

Se aprecia una buena concordancia entre valores teóricos y experimentales dentro del rango del error estadístico. b) Análisis de los resultados al variar la frecuencia de giro ω Los siguientes gráficos muestran por un lado la variación del intervalo de tiempo δ entre frecuencia máxima – mínima al variar la frecuencia angular de giro ω y por el otro la diferencia entre estas frecuencias, también al variar la frecuencia angular de giro ω para una distancia d constante. Los


valores del intervalo de tiempo con sus respectivos errores se obtuvieron de igual forma que lo que se lo hizo en la sección anterior. La diferencia de amplitud es la diferencia entre la frecuencia máxima y la siguiente mínima. Los resultados se calcularon de igual forma que en los intervalos de tiempos, o sea, se hizo el promedio entre cuatro de estas diferencias de frecuencia y se tomó como error a la desviación estándar. La curva de diferencia de frecuencia se realizó debido a que esta magnitud varía al variar la frecuencia de giro, y no así, al modificar la distancia mínima entre la fuente y el micrófono. Las curvas teóricas se realizaron utilizando las expresiones (2.7) y (2.9) .La distancia mínima d entre la fuente y el micrófono fue de 3 cm y el radio de giro de la fuente de 75 cm por lo que la distancia D fue de 78cm. Se muestra en los gráficos el valor de χ2 relacionando los resultados experimentales con las predicciones teóricas.

2 4 6 8 10 12 14 16 180.00

0.02

0.04

0.06

0.08

0.10

0.12

0.14

0.16

0.18

0.20χ2 = 0.12

Inte

rval

o de

Tie

mpo

(s)

ω (Hz)

Resultados obtenidos Curva teórica

Figura 4.6: Intervalo de tiempo δ entre frecuencia máxima – mínima en función de la frecuencia angular de giro. Comparación entre resultados experimentales y predicción teórica dada por la ecuación (2.9)

2 4 6 8 10 12 14 16 18

0.05

0.10

0.15

0.20

0.25

0.30

0.35

0.40

χ2 = 0.03

Dife

renc

ia d

e Fr

ecue

ncia

(KH

z)

ω (Hz)

Resultados obtenidos Curva teórica

Figura 4.7: Diferencia entre frecuencia máxima y mínima en función de la frecuencia angular de giro. Comparación entre resultados experimentales y predicción teórica dada por la ecuación (2.7) Nuevamente se observa una buena concordancia entre las curvas teóricas y los resultados experimentales, siempre dentro del error estadístico.


c) Curvas de frecuencia en función del tiempo Se realizaron las curvas de frecuencia en función del tiempo para distancias mínimas entre fuente y micrófono de 6 cm, 31 cm, 51 cm y 81 cm. Se ajustó cada una de ellas con su correspondiente curva teórica, dada por la ecuación (2.6). Sin embargo, los resultados teóricos y experimentales mostraron una pequeña discrepancia. Para encontrar la frecuencia de giro que mejor ajustaba los datos experimentales, se realizó la minimización de la función χ2 variando el parámetro ω (ver Apéndice B, ecuación B.2). Se hizo esto sólo cuando la distancia mínima fue de 6 cm, encontrándose como frecuencia angular de giro a 13,4 ± 0,4Hz. Esta misma frecuencia se empleó luego para ajustar las demás curvas. Los gráficos muestran además los correspondientes valores de χ2. Inicialmente, con objetivos ilustrativos, se mostrará el sonograma de donde se obtuvieron los puntos superpuesto con la curva teórica en el caso en el cual la distancia mínima fue de 6 cm.

Figura 4.8: Sonograma en el caso de una distancia mínima entre fuente y microfono de 6 cm. Superpuesto en linea negra se muestra la curva teórica.

A continuación se muestran los gráficos en las distancias indicadas de frecuencia en función del tiempo con los puntos obtenidos empleando el programa xyExtract superpuestos a las correspondientes curvas teóricas.

5.0 5.5 6.0 6.5 7.0 7.55.10

5.15

5.20

5.25

5.30

5.35

5.40

5.45

5.50

5.55

5.60Distancia mínima = 6 cm

χ2 = 0.011

Frec

uenc

ia (K

Hz)

Tiempo (s)

Resultados experimentales Curva teórica

6.5 7.0 7.5 8.0 8.5 9.05.10

5.15

5.20

5.25

5.30

5.35

5.40

5.45

5.50

5.55


χ2 = 0.013

Frec

uenc

ia (K

z)

Tiempo (s)


Figura 4.9a Figura 4.9b


7.0 7.5 8.0 8.5 9.0 9.55.10

5.15

5.20

5.25

5.30

5.35

5.40

5.45

5.50

5.55


χ2 = 0.012Fr

ecue

ncia

(KH

z)

Tiempo (s)


7.4 7.6 7.8 8.0 8.2 8.4 8.6 8.8 9.0 9.2 9.4 9.6

5.15

5.20

5.25

5.30

5.35

5.40

5.45

5.50

5.55

5.60

5.65

Distancia mínima = 81 cm

χ2 = 0.010

Frec

uenc

ia (K

Hz)

Distancia (cm)


Figura 4.9c Figura 4.9d

Figuras 4.9a, 4.9b, 4.9c, 4.9d. Frecuencia en función del tiempo para distancias mínimas entre fuente y micrófono de 6cm, 31 cm,51 cm y 81 cm respectivamente. Comparación entre resultados experimentales y predicción teórica dada por la ecuación (2.6) Se aprecia una muy buena concordancia entre los valores experimentales y las curvas teóricas. También se aprecia aquí que el intervalo de tiempo δ entre frecuencia máxima – mínima aumenta al aumentar la distancia, siendo muy parecido el intervalo entre frecuencia máxima – mínima con el intervalo entre frecuencia mínima – máxima cuando la distancia mínima fue de 81 cm. d) Análisis de las reflexiones utilizando el método de las imágenes En esta sección se utiliza el método de imágenes para interpretar las reflexiones que se producen en una pared. Para esto se colocó una pared perpendicular a la dirección formada por el micrófono y el centro de giro de la fuente, con este último ubicado entre el micrófono y la pared. Se dirá directa a la señal proveniente directamente de la fuente e imagen a la que llega como reflejo desde la pared. La distancia mínima entre fuente y micrófono fue, para la señal directa de 3 cm, y debido a que la distancia mínima entre fuente y pared fue de 5 cm, la distancia mínima entre señal imagen y micrófono resulto ser de 1513 cm. El radio de giro fue de 75 cm (esto vale también para la señal imagen), y para ajustar con las curvas teóricas se tomó como frecuencia angular de giro a la obtenida en la sección anterior, que fue de 13,465 Hz. La bondad del ajuste del modelo de las imágenes se cuantificó usando la función χ2 . A continuación se mostrará el sonograma en esta situación, y superpuesto a él se mostrarán las curvas teóricas que se obtienen del método de imágenes, en negro la que corresponde a la señal directa y en rojo la que corresponde a la señal imagen.

Figura 4.10. Sonograma en el estudio del método de imágenes. Superpuesto en línea negra y en línea roja se muestran las curvas teóricas de las señales directa e imagen respectivamente.


A continuación se muestra el gráfico de frecuencia en función del tiempo, tanto para la señal directa como para la señal imagen, con los puntos obtenidos del sonograma anterior empleando el programa xyExtract, superpuestos a las correspondientes curvas teóricas (2.6) y (2.11)

10.5 11.0 11.5 12.0 12.5 13.05.10

5.15

5.20

5.25

5.30

5.35

5.40

5.45

5.50

5.55 Imagen (Experimental) Imagen (Teórica)

Directa: χ2 = 0.023Imagen: χ2 = 0.016

Y A

xis

Title

X Axis Title

Directa (Experimental) Directa (Teórica)

Figura 4.11. Frecuencia en función del tiempo para la señal directa y la señal imagen. Comparación entre resultados experimentales y predicción teórica.

Se aprecia una muy buena concordancia entre los valores experimentales y las curvas teóricas en ambos casos, validando la hipótesis en cuanto al método de imágenes para tratar las reflexiones en una pared. También se aprecia aquí muy claramente como varia el intervalo de tiempo entre frecuencia máxima – mínima al aumentar la distancia mínima entre fuente y micrófono, siendo este muy pequeño en el caso directo y prácticamente igual al de intervalo de tiempo entre frecuencia mínima –máxima en el caso imagen. IV. Conclusiones En este trabajo estudiamos el efecto Doppler de una fuente de sonido en un movimiento circular, para el caso libre de reflexiones y para el caso de reflexiones producidas en una pared cercana. Para saber cómo varía la frecuencia en función del tiempo se emplearon mapas de frecuencia – tiempo llamados sonogramas y se hicieron mediciones sobre estos diagramas. Para el caso del movimiento libre, se estudió como varían distintos parámetros de la curva de frecuencia en función del tiempo variando por un lado la distancia mínima d entre fuente y micrófono y por otro lado la frecuencia angular de giro ω de la fuente. Las mediciones del intervalo de tiempo δ en función de la distancia D se ajustaron con la ecuación (2.9) y se obtuvo para la función χ2 un valor de χ2 = 4,5% . Las mediciones de δ y de ∆f en función de la frecuencia ω se ajustaron utilizando las ecuaciones (2.9) y (2.7) obteniéndose para la función χ2 los valores χ2 = 12 % y χ2 = 3% respectivamente. Para las distancias d = 6 cm, 31 cm, 51 cm y 81 cm se hicieron mediciones de la frecuencia en función del tiempo. La frecuencia de giro se encontró minimizando la función χ2 para uno de los casos y utilizando este resultado en el resto de las mediciones. Se ajustaron los resultados usando la ecuación (2.6) obteniendo valores de χ2 del orden del 1% en todos los casos.


Para el caso del experimento en las cercanías de una pared, se trataron las reflexiones con el método de las imágenes. Las mediciones realizadas sobre el sonograma se ajustaron con la ecuación (2.11) y se obtuvo una buena correspondencia entre las predicciones de este modelo y las resultados (χ2 = 1,6 %). Los buenos valores obtenidos para la función Chi-cuadrado, del orden del 1% utilizando la frecuencia corregida y de 1,6% para el caso de las reflexiones, reflejan el buen acuerdo entre la teoría y el experimento y podemos decir que el método de las imágenes es adecuado para tratar las reflexiones de una onda sonora en una pared. Apéndice A: La componente paralela de la velocidad de rotación v// Podemos escribir la ecuación (2.1) como

( ))cos()cos()(// tsentsenvv f ωββω +−= (A.1) De la figura 2.1 se ve que l sen(β) = r sen(ωt), de manera que (A.1) se escribe como

+−=lrttsenvv f )cos()cos()(// ωβω

+−=l

rtltsenvv f)cos()cos()(//

ωβω (A.2)

Como se ve en la figura 2.1, cos(β) l + cos(ωt) r = D, así que (A.2) queda

−=lDtsenvv f )(// ω (A.3)

Como l 2 = r 2 + D 2 – 2rD cos(ωt) y α = r / D tenemos l 2 =D 2 (α 2 + 1 2 – 2α cos(ωt)), con lo que podemos escribir finalmente

)cos(211)(// t

tvsenvωαα

ω−+

−= (A.4)

Apéndice B: La función Chi-cuadrado Supongamos que tomamos una serie de N mediciones de dos magnitudes cuya relación deseamos determinar. Como resultado de nuestras mediciones obtenemos una serie de N ternas de la forma (xi, yi,σi) donde σi es la incertidumbre asociada a la medición de yi y suponemos que la incertidumbre de xi es despreciable. Supongamos tambien que el modelo teórico que ajusta los datos viene dado por la función y(xi) = f(xi;a,b,c,...) donde a,b,c etc. son los n par parámetros del modelo [5].


En el caso mas general definimos la función Chi-cuadrado como

( )∑=

−=N

iiii xyyw

1

22 )(χ (B.1)

La cantidad (yi – y(xi))2 representa las desviaciones entre el valor medido yi y el predicho por el modelo y(xi) y wi es un factor de peso de cada triada de datos (xi, yi,σi). De esta manera χ2 es una medida de la desviación total de los valores medidos respecto de los predichos teóricamente, ponderados por un factor wi. Suele usarse para wi el valor wi = 1/σi

2 , aunque a veces es útil emplear otros factores de peso como wi = 1/yi

2 o wi = 1/yi. Si los datos no tienen asignado un error de medición σi, puede elegirse tambien wi = 1 La funcion Chi-cuadrado se puede utilizar para encontrar el mejor valor de un parámetro y para cuantificar la bondad de un ajuste teórico de resultados experimentales. a) Encontrar el mejor valor de un parámetro Podemos usar la funcion Chi- cuadrado para encontrar el mejor valor de un parámetro de nuestro modelo. Supongamos que nuestro modelo depende de un solo parámetro y(xi) = f (xi,a) y deseamos obtener el valor de a que mejor ajuste los resultados experimentales. El mejor valor de a es aquel que minimiza la función χ2 [5]. Usando una planilla de cálculo puede calcularse la función χ2 y modificar el valor de a hasta encontrar el valor a0 que minimice la función χ2. Para ese valor a0

)( 022

min aχχ = (B.2) Un método aproximado para calcular la incertidumbre del parámetro a es considerar cuál es el incremento que debe sufrir a para que la función χ2 crezca en una unidad [5]. Tal incremento se toma como el error del parámetro a. Es decir

)(1 022

min aa ∆+=+ χχ (B.3) b) Bondad del ajuste Tambien podemos usar la función Chi-cuadrado para tener una medida de la bondad o calidad del ajuste de los resultados dado por el modelo. Una medida de la bondad del ajuste esta dada por el valor de χ2

ν [6]

( )∑=

−=N

iiii xyyw

1

22 )(1ν

χν (B.4)

donde ν = N – npar es el número de grados de libertad. Si tomaos wi = 1 / σi podemos decir que si todos los datos experimentales yi tienen desviaciones respecto del modelo y(xi) que no superen el error σi, nuestro modelo es el adecuado para describir nuestras observaciones. En ese caso cada término de la sumatoria será del orden de la unidad y podemos decir que si χ2

ν es del orden de la unidad o menor el modelo propuesto para explicar los datos experimentales es el adecuado.


Si los datos experimentales no tienen asignado un error, podemos elegir wi = 1 / yi2 y entonces

( )

∑=

−=

N

i

ii

iy

xyy1

2

22 )(1

νχν (B.5)

En este caso cada término de la sumatoria es el cuadrado de las desviaciones relativas y χ2

ν representa un promedio de esas desviaciones. Multiplicando por 100, tenemos una cantidad χ2

ν que mide el porcentaje promedio de las desviaciones relativas entre el modelo y los datos experimentales. REFERENCIAS [1] Sears, Zemansky; Fisica Universitaria, Vol 1, Addison Wesley Longman (1998)[2] S. Spicklemire, M. Coffaro, The treatment of reflections in a Doppler measurement using the method of images; Am. J. Phys, 74 (2006)

[3] M. Saba, R. Rosa; The Doppler effect of a sound source moving in a circle; The Physics Teacher; 41, 2003

[4] A. Romero, M. E. Capoulat; Análisis del efecto Doppler para una fuente en movimiento circular; Laboratorio 5, Dpto Física, FCEyN, UBA

[5] S. Gil, E. Rodríguez; Métodos cuantitativos de análisis gráficos; Física re-creativa

[6] S. Gil, E. Rodríguez; Análisis avanzado, Bondad de ajuste; Física re-creativa

[7] xyExtract Graph Digitizer, http://www.gold-software.com/download5149.html

[8] Existen muchos programas que generan un sonograma a partir de un archivo .wav. Por ejemplo Audacity (disponible en www.superarchivos.com )¨ o programa SpectroGram versión 11.2 (by Richard Horne http://www.visualizationsoftware.com/gram.html). Ver también: What are spectrograms? http://cslu.cse.ogi.edu/tutordemos/SpectrogramReading/spectrogram.html

http://www.gold-software.com/download5149.html

http://www.visualizationsoftware.com/gram.html

http://cslu.cse.ogi.edu/tutordemos/SpectrogramReading/spectrogram.html

efecto doppler de una fuente de sonido en un movimiento circular

Documents