Cómo encontrar la traducción horizontal de un gráfico senoidal, dada una tabla de datos

Suponiendo que los datos realmente siguen un patrón sinusoidal, entonces sí, es posible interpolarlos.

Para hacer esto, debemos conocer las características de los gráficos sinusoidales. Analicemos esto:

if [math] T = a \ sin {(b (tc))} + d [/ math], entonces

a es la amplitud de la función.

La amplitud es el rango de valores que puede tener un gráfico sinusoidal. Un gráfico sinusoidal puede tener valores de + A a -A. Entonces, el rango de los datos es el doble de la Amplitud.

[matemáticas] A = \ frac {Rango} {2} [/ matemáticas]

La temperatura varía de 14 (el mínimo) a 27 (el máximo). El rango será 27-14 = 13.

Así:

[matemáticas] a = A = \ frac {Rango} {2} = \ frac {13} {2} = 6.5 [/ matemáticas]

d es el desplazamiento de temperatura

Podemos establecer d en un valor arbitrario para que coincida con los datos de temperaturas. Por ejemplo, sabemos que el valor mínimo es 14. Esto ocurre cuando el gráfico sinusoidal está en el mínimo, entonces el valor de la onda sinusoidal está en -A. Esto significa:

[matemáticas] T = a \ sin {(b (tc))} + d [/ matemáticas]

[matemáticas] 14 = -6.5 + d [/ matemáticas]

[matemáticas] d = 20.5 [/ matemáticas]

b es la frecuencia de fase de la función

Una onda sinusoidal tiene una frecuencia o período que describe cuándo se repite el valor.

Para ilustrar mejor esto, una onda sinusoidal regular se repite después de que el ángulo alcanza [matemática] 2 \ pi [/ matemática]. A esto lo llamamos 1 período. Una frecuencia de fase (o período) siempre tiene la forma de [matemática] 2 \ pi f [/ matemática] o [matemática] \ frac {2 \ pi} {P} [/ matemática] donde f es la frecuencia y P es el período. Podemos interesarnos en detalles matemáticos si está interesado, pero por ahora vamos a explicarlo brevemente.

Ahora es más fácil para nosotros encontrar el período de los datos (más intuitivo). Podemos decir que un período es un tiempo que la función tarda desde el valor máximo para alcanzar nuevamente un valor máximo. En los datos, solo vimos el valor máximo una vez. Entonces, tenemos que improvisar.

Si inspecciona el gráfico sinusoidal, la función tarda la mitad del período en pasar del valor mínimo al valor máximo. Tenemos el valor mínimo y máximo, así que empecemos por eso. Si tomamos Jan como t = 1 y Feb como t = 2, entonces el valor mínimo ocurre en t = 2 y el valor máximo ocurre en t = 7. La diferencia en t (o [math] \ delta t [/ math] es 5. Esta es la mitad del período, por lo que el período en sí es 10.

Así:

[matemáticas] b = \ frac {2 \ pi} {P} = \ frac {2 \ pi} {10} = \ frac {\ pi} {5} [/ matemáticas]

c es el desfase de la función

Similar a d, si d es el desplazamiento de la temperatura, entonces c es el desplazamiento de la fase. O en este caso, el desplazamiento de la hora. Elegimos c para que el gráfico seno se alinee bien con los datos por tiempo.

¿Recuerdas nuestro valor mínimo? El valor mínimo ocurre en el gráfico seno cuando la fase es [matemática] \ frac {3 \ pi} {2} [/ matemática]. Es equivalente al decir que ocurre cuando el tiempo está en [math] \ frac {3} {4} [/ math] de su período. Podemos elegir cualquier definición. Por ejemplo, si vamos usando el ángulo de fase:

[matemáticas] \ frac {3 \ pi} {2} = \ frac {\ pi} {5} (tc) [/ matemáticas]

[matemáticas] c = t-7.5 [/ matemáticas]

el valor mínimo ocurre en febrero, entonces t = 2. Sustituya esto, obtenemos:

[matemáticas] c = -5.5 [/ matemáticas]

Si vamos usando periodo:

[matemáticas] \ frac {3} {4} P = tc [/ matemáticas]

[matemáticas] c = t- \ frac {3} {4} P = 2- \ frac {3} {4} 10 [/ matemáticas]

[matemáticas] c = -5.5 [/ matemáticas]

Para más información, c está realmente en el módulo P. Por lo tanto, es perfectamente válido tener:

[matemáticas] c = 10-5.5 = 4.5 [/ matemáticas]

o incluso:

[matemáticas] c = 10 \ cdot10-5.5 = 94.5 [/ matemáticas]

Esto se debe a que produce el mismo ángulo de fase para esta onda sinusoidal particular.

usemos c = 4.5, porque esperamos que c sea un número positivo

Resumen

Su función está en la forma:

[matemáticas] T = 6.5 \ sin {(\ frac {\ pi} {5} (t-4.5))} + 20.5 [/ matemáticas]

Solo tenga en cuenta que, se ajustará muy bien para el valor mínimo y máximo con los datos, pero puede ser un poco diferente para otros valores, ya que es solo una interpolación. Sería interesante calcular la varianza, pero no es el propósito de las preguntas.

Espero eso ayude.

Nota [EDITAR]
Como puede ver en la respuesta del usuario de Quora a ¿Cómo puede encontrar la traducción horizontal de un gráfico senoidal, dada una tabla de datos? Se mencionó que el período es 12, por lo tanto, produce [matemáticas] b = \ frac {\ pi} {6} [/ matemáticas], que es diferente de la respuesta a la que llegué.

Para repetirme de nuevo, esta respuesta aquí se ajustará perfectamente para el valor mínimo y máximo con los datos, pero podría estar desactivada para otros valores, ya que es solo una interpolación. Por otro lado, estaba perfectamente bien asumir que el período es 12. Los períodos que elijas dependerán de los supuestos que hagas.

Llegué a P = 10 con la suposición de que se necesitaron 5 t para que la función pasara del mínimo al máximo más cercano. Voy con este enfoque, porque solo quiero encontrar la P solo a partir de los datos y sin otras suposiciones. Para ilustrar esto, imagine si los datos realmente tienen una frecuencia mucho más alta (período más corto) pero aún se repiten después de 12 meses. Por ejemplo, si P = 6, obtendremos 2 valores máximos máximos y / o 2 valores mínimos máximos. Si en este caso elige P = 12, producirá un mal modelo a partir de suposiciones razonables. Es cierto que la función se repetirá después de 12 t, pero en algunos lugares estará muy apagada.

Entonces, en mi opinión, no quiero cambiar mi P a 12. Los lectores pueden tener más información sobre por qué llegamos a resultados diferentes. Solo recuerde probar el modelo con el conjunto de datos dado: p. Para estas preguntas en particular, creo que también está bien comparar el modelo haciendo algunas pruebas de regresión, eso sería interesante. Debido a que el conjunto de datos no es muy grande, es bastante posible calcular la suma residual de cuadrados sin usar tecnología adicional

Voy a tomar sin tecnología para significar papel y lápiz.

El primer paso es trazar los valores en una cuadrícula. Debería terminar con algo que se parece a un gráfico senoidal que se desplaza hacia arriba y a la derecha.

En la ecuación a representa la amplitud. Por lo tanto, la diferencia entre el pico y el valle será 2a.

Mi cálculo en la cabeza es que a = 6.5

El valor d el valor que se agrega para mover la curva hacia arriba. Lo encontraría promediando los números de pico y valle.

Mi cálculo en la cabeza es que d = 20.5

Ahora distribuya b al término en el lado para que tenga seno (bt + bc) Ahora suponiendo que la frecuencia es 1 ciclo / año entonces, dado que t está en meses y comienza con 1, quiere 12b = 2pi o b = pi / 6 . ¿Por qué 2 pi. Porque la función seno se repite cada 2 pi.

Ahora para encontrar bc, que es la cantidad de desplazamiento de la curva senoidal hacia la derecha. La forma en que lo haría es dibujar la línea T = d en el gráfico. Entonces soltaría un eje perpendicular al tiempo (x) esta vez es bc si tengo mis signos correctos. Sabes b para que puedas calcular c.

Mi cálculo en la cabeza es que bc = 5. por lo tanto c = 30 / pi

La verificación final es comenzar a sustituir los números en la ecuación para asegurarse de que no haya un error en algún lugar aquí. En particular, puedo haber cometido un error +/- en este cálculo en alguna parte. Voy a decir que usar una calculadora para verificar los números no es usar tecnología, ya que podría buscar los valores en una tabla o incluso usar la expansión de la serie Taylor para calcular los valores del seno de los distintos meses. Sería tedioso pero factible.

Mis números son una aproximación / conjetura.

Comienzan molestamente el tiempo en 1 en lugar de 0, por lo que [math] t = 1 [/ math] significa enero, [math] t = 2 [/ math] feb, etc. También están siendo imprecisos sobre si [math] t = 1 [/ math] significa el inicio o mediados de enero; Asumamos el medio.

De todos modos, este es un gráfico de temperatura anual, por lo que presumiblemente el período de la onda sinusoidal es de 1 año o 12 meses. El argumento seno tiene que atravesar [matemáticas] 2 \ pi [/ matemáticas] durante este período, por lo que vemos [matemáticas] b = \ frac {2 \ pi} {12} = \ frac {\ pi} {6} [/ matemáticas].

La función seno tiene su máximo en [matemáticas] 90 ^ \ circ = \ frac {\ pi} {2} [/ matemáticas] y el mínimo en [matemáticas] -90 ^ \ circ = – \ frac {\ pi} {2} [/matemáticas]. La tabla tiene el mínimo en febrero = 2 y el máximo en julio = 7, 5 o 7 meses de diferencia en lugar del esperado 6. No se supone que usemos tecnología aquí, por lo que podríamos decir min en enero = 1, max en Julio = 7 y listo. Intentemos ser un poco más precisos observando el promedio de cada par de meses adyacentes, dando un empate por minuto en enero / febrero ([matemáticas] t = 1.5 [/ matemáticas] según mi suposición al comienzo) y febrero / marzo ( [matemáticas] t = 2.5 [/ matemáticas]) y un máximo en julio / agosto ([matemáticas] t = 7.5 [/ matemáticas]). Así que vayamos por min en [math] t = 1.5 [/ math] y max en t = [math] 7.5 [/ math] y no nos preocupemos más.

Sabemos que [matemática] b (tc) = \ frac {\ pi} {2} [/ matemática] al máximo, entonces [matemática] c = t – \ frac {\ pi} {2b} = 7.5 – \ frac { \ pi} {2 \ pi / 6} = 4.5 [/ matemáticas]. Marcando [math] t = 1.5, b (tc) = \ frac {\ pi} {6} (1.5-4.5) = – \ frac {\ pi} {2} [/ math], el punto mínimo como se esperaba [math ] \ \ marca de verificación [/ math].

También podríamos terminar el problema. El rango de los datos va de 14 a 27 y el rango de pecado va de -1 a 1, entonces [matemáticas] (- 1) a + d = 14 [/ matemáticas] y [matemáticas] (1) a + d = 27. [/ Matemática] Sumando, [matemática] 2d = 41 [/ matemática] o [matemática] d = 20.5 [/ matemática], entonces [matemática] a = 6.5 [/ matemática].

Poniendolo todo junto:

[matemática] T = 6.5 \ sin (\ frac {\ pi} {6} (t-4.5)) + 20.5 [/ matemática]

Recreando la tabla de la fórmula, obtenemos

Ene [matemáticas] t = 1 \ \ T = 14.2 [/ matemáticas]

Feb [matemáticas] t = 2 \ \ T = 14.2 [/ matemáticas]

Mar [matemáticas] t = 3 \ \ T = 15.9 [/ matemáticas]

Abr [matemáticas] t = 4 \ \ T = 18.8 [/ matemáticas]

Mayo [matemáticas] t = 5 \ \ T = 22.2 [/ matemáticas]

Jun [matemáticas] t = 6 \ \ T = 25.1 [/ matemáticas]

Jul [matemáticas] t = 7 \ \ T = 26.8 [/ matemáticas]

Ago [matemáticas] t = 8 \ \ T = 26.8 [/ matemáticas]

Sep [matemáticas] t = 9 \ \ T = 25.1 [/ matemáticas]

Oct [matemáticas] t = 10 \ \ T = 22.2 [/ matemáticas]

Nov [matemáticas] t = 11 \ \ T = 18.8 [/ matemáticas]

Dic [matemáticas] t = 12 \ \ T = 15.9 [/ matemáticas]

Parece sesgado un poco alto, pero bastante bueno dado que no tiene tecnología, así que lo dejaré.

¡Transformada de Fourier!
Calcular
[matemáticas] f_1 = \ sum_ {k = 1} ^ {12} (t_k \ veces e ^ {i 2 \ pi k / 12}) [/ matemáticas]
que es un número complejo que representa la frecuencia principal de su patrón.
Calcule la fase de este número complejo y multiplíquelo con [matemáticas] 12 / (2 \ pi) [/ matemáticas] para obtener el mes en el que este patrón alcanza su punto máximo. El seno comienza tres meses antes.

Código Matlab

t = 1:12;
y = [15 14 15 18 21 25 27 26 24 20 18 16];
Fvect = exp (1i * 2 * pi * (1:12) / 12);
F0 = suma (y. * Fvect);
m = mod (ángulo (F0) / (2 * pi) * 12-1,12) +1;
figura;
barra (t, y); Espere;
a = abs (F0) / 6;
b = 2 * pi / 12;
c = (m-3);
d = media (y);
plot ([t], a * sin (b * (t – c)) + d, ‘r’, ‘LineWidth’, 3);

m = 7.45

Por lo tanto, alcanza su punto máximo justo antes de finales de julio.

El seno comienza tres meses antes, por lo tanto c = 4.45

(y a = 6.2028, b = 0.5236, c = 4.4504, d = 19.9167)

No quiero hacer la pregunta, sino dar algunas pistas e intuición:

Si [matemáticas] T = a sin (b (t – c)) + d [/ matemáticas]

Qué son:

a) los mayores cambios de

b) el período de los cambios en

c) la “fase” de los cambios

d) el promedio

valores de T