¿Por qué el biespectro no se usa comúnmente en física experimental?

Los espectros de potencia, los espectros de coherencia y las funciones de transferencia lineal son herramientas omnipresentes de la física experimental. Sin embargo, nuestros instrumentos a menudo retienen pequeños efectos no lineales que pueden contaminar las mediciones. Parece que los espectros de orden superior, en particular el biespectro , serían herramientas ideales para investigar interacciones no lineales. No obstante, en realidad nunca los he visto ponerlos en práctica en la física experimental.

Por ejemplo, considere la coherencia (dominio de frecuencia):

C X y ( F ) = X ( F ) Y ( F ) X ( F ) X ( F ) Y ( F ) Y ( F )

La bicoherencia considera no dos sino tres señales, y busca correlaciones entre oscilaciones en frecuencias F 1 y F 2 combinando no linealmente para producir una señal en F 1 + F 2 :

C X y z ( F 1 , F 2 ) = X ( F 1 ) Y ( F 2 ) Z ( F 1 + F 2 ) X ( F 1 ) X ( F 1 ) Y ( F 2 ) Y ( F 2 ) Z ( F 1 + F 2 ) Z ( F 1 + F 2 )

...que parece algo útil.

¿Por qué el biespectro y la bicoherencia no se usan con más frecuencia en la física experimental?

Estoy pensando específicamente en el dominio del tiempo, sistemas de múltiples entradas / múltiples salidas donde uno busca acoplamientos no lineales entre varias señales.

Uno de los principales éxitos de Google sobre el tema es para la caja de herramientas Matlab Higher Order Spectral Analysis (HOSA) , que parece un buen recurso (aunque parece que ya no se mantiene y ahora sufre de bit-rot).

Respuestas (3)

Se usa mucho en cosmología. A menudo, en una aproximación decente, las cantidades que tratamos de medir en cosmología (por ejemplo, mapas de polarización y temperatura CMB, distribuciones de galaxias) son realizaciones de procesos aleatorios gaussianos en una aproximación decente, pero tienen (o se predice que tendrán) interesantes no- Características gaussianas en algún nivel bajo. La gente estima el biespectro de estas cosas todo el tiempo para controlar varios efectos no gaussianos.

Gracias por la respuesta. ¿Estás hablando, por ejemplo, de una función de correlación de tres puntos en la esfera? ¿Tiene alguna referencia a documentos que muestren estas cosas?
Sí, en el contexto del fondo de microondas, una función de 3 puntos en la esfera es exactamente lo que se calcula. Por lo general, se hace en el espacio armónico esférico, aunque no en el espacio real. Consulte, por ejemplo, prd.aps.org/abstract/PRD/v80/i4/e043510 . Creo que la gente habla sobre el biespectro de varias funciones (perturbaciones de densidad, por ejemplo) en el espacio tridimensional también. Una vez más, esto generalmente se hace en el espacio de Fourier, no en el espacio real.

Estos objetos se utilizan todo el tiempo. Las matemáticas se hacen en términos de campos cuánticos, lo que hasta cierto punto oculta lo que está pasando.

Por ejemplo, su "coherencia (dominio de frecuencia)" es un coeficiente de correlación , que está normalizado, mientras que los físicos suelen trabajar en términos de funciones de correlación, que normalmente no lo son, pero en gran medida equivalen a lo mismo. tus observables X ( F 1 ) , etc., se construyen como funciones de frecuencia, sin embargo, este es un objeto singular en la teoría cuántica de campos. En la teoría cuántica de campos, en cambio, construimos observables ϕ ( F 1 ) , etc., como funcionales de funciones de prueba F 1 ( X ) , etc. Una opción singular sería F 1 ( X ) = Exp ( i F 1 X ) , que hace ϕ ( F 1 ) esencialmente el mismo objeto que su X ( F 1 ) ; es singular , sin embargo, porque F 1 ( X ) no es integrable al cuadrado.

Otra opción de función de prueba singular es, por supuesto, d ( X y ) , que da el valor del campo en un punto, que podríamos escribir en sus términos como algo así X ( y ) . Para un campo cuántico, este también es un objeto bastante singular.

De hecho, cuando dices X ( F 1 ) , lo que realmente quieres decir es X ( F ) d F , en un pequeño rango de frecuencias, y en los detalles matemáticos y experimentales esto debe tenerse en cuenta. Hacer que todo sea preciso requiere que sepamos cuáles son los rangos de frecuencia de cada una de las mediciones, que un experimentador debe caracterizar o debe leer de las hojas de datos del fabricante. Con más detalle aún, tendremos que construir una función de peso, diciendo que las frecuencias cercanas a F 1 todavía están registrados por el dispositivo de medición, pero no tanto como cerca F 1 . Bien podemos tomar la función de peso, como primera aproximación, como gaussiana. Esto corresponde a tomar la función de prueba F 1 ( X ) ser tan gaussiano. El análisis de señales generalmente llama a la función de prueba una función de ventana . Puede ser difícil familiarizarse con las funciones de prueba o ventana, pero creo que vale la pena llegar allí.

En estos términos, su C X y z es una elección particular de función de 3 puntos (normalizada). La elección de F 1 + F 2 para la tercera frecuencia, por supuesto, no es necesaria, podemos considerar correlaciones de 3 puntos entre tres frecuencias cualesquiera F 1 , F 2 , F 3 (y sus alrededores). En la teoría cuántica de campos, representaríamos la función de correlación de 3 puntos, en el estado de vacío, como 0 | ϕ ( F 1 ) ϕ ( F 2 ) ϕ ( F 3 ) | 0 . Reemplace el vector de vacío por algún otro vector de estado, si lo desea.

En el caso particular en que los observables de campo cuántico son mutuamente conmutativos, se puede entender que generan medidas de probabilidad que corresponden a una descripción equivalente en términos de medidas de probabilidad sobre variables aleatorias clásicas y, por lo tanto, con bastante precisión a un análisis de señal estocástica. Cuando los observables de campo cuántico no conmutan, todo se vuelve mucho más complicado, pero se puede mantener un remanente del punto de vista del procesamiento de señales.

Hay una matemática de campos aleatorios que se usa en cosmología porque generalmente no es necesario preocuparse por la incompatibilidad de medición en ese contexto. Los matemáticos generalmente presentan el análisis de señales en términos de espacio de Hilbert a menos que estén escribiendo para una audiencia de ingeniería.

La perspectiva del análisis de señales en campos cuánticos no se discute mucho, AFAIK, pero lo encuentro extremadamente efectivo. Es en gran medida una perspectiva de campos en lugar de partículas en QFT, lo que lo pone en desacuerdo con la presentación tradicional de libros de texto de QFT en física. Estoy trabajando en la construcción de campos cuánticos que interactúan de manera motivada por este enfoque (aunque hay tantos enfoques por ahí, y mis cosas rara vez funcionan lo suficientemente bien). Parte de mi perspectiva sobre los campos aleatorios/cuánticos que no interactúan se puede encontrar en mis artículos publicados (aunque son defectuosos).

Mis conjeturas son:

  1. Se necesitan muchos más datos para producir un biespectro de bajo ruido que un espectro de potencia (?).
  2. Los resultados son confusos de interpretar (o al menos desconocidos). Parte de esto se debe a la gran cantidad de redundancia en la salida.
  3. La gran cantidad de datos en el biespectro puede ser enorme. ¡Un biespectro con la misma resolución de frecuencia que un espectro de potencia requiere el cuadrado de la cantidad de memoria!
  4. Cualquier buen reloj parecerá ser biespectralmente coherente. Los transitorios también crean una bicoherencia espuria.