Взаимная спектральная плотность мощности¶

Взаимная спектральная плотность мощности — это функция, описывающая распределение мощности совместных колебаний двух случайных процессов по частотам. Она характеризует степень статистической связи между двумя сигналами на каждой частоте, включая информацию о фазовом сдвиге, и используется для анализа систем передачи сигналов, вибраций и в радиотехнике.

Математически взаимная спектральная плотность мощности (ВСПМ) определяется как преобразование Фурье от взаимной корреляционной функции \( R_{xy}(\tau) \) двух стационарных случайных процессов \( x(t) \) и \( y(t) \):

\[ S_{xy}(f) = \int_{-\infty}^{\infty} R_{xy}(\tau) e^{-j2\pi f \tau} d\tau \]

В отличие от обычной спектральной плотности мощности (СПМ), которая всегда является вещественной и неотрицательной, ВСПМ является комплексной величиной. Её действительная часть называется коспектром, а мнимая — квадратурным спектром. Модуль ВСПМ показывает степень когерентности сигналов на данной частоте, а аргумент (фаза) — средний фазовый сдвиг между компонентами сигналов.

¶Свойства

Основные свойства взаимной спектральной плотности мощности:

• Симметрия: \( S_{xy}(f) = S_{yx}^*(f) \), где звёздочка означает комплексное сопряжение.
• Неравенство Коши — Буняковского: \( |S_{xy}(f)|^2 \le S_{xx}(f) S_{yy}(f) \), где \( S_{xx}(f) \) и \( S_{yy}(f) \) — автоспектральные плотности мощности соответствующих процессов.
• Нормировка: на основе ВСПМ вычисляется функция когерентности \( \gamma^2_{xy}(f) = \frac{|S_{xy}(f)|^2}{S_{xx}(f) S_{yy}(f)} \), принимающая значения от 0 до 1.

¶Методы оценки

На практике ВСПМ оценивается по конечным реализациям сигналов. Основные подходы:

• Периодограммный метод: вычисление преобразования Фурье от оценок взаимной корреляционной функции с применением оконных функций для сглаживания.
• Метод Уэлча: усреднение периодограмм по перекрывающимся сегментам данных, что снижает дисперсию оценки.
• Параметрические методы: авторегрессионное моделирование (метод Берга, ковариационный метод), позволяющее получить более гладкие оценки при коротких выборках.

Выбор метода зависит от длины записи, требований к частотному разрешению и вычислительных ресурсов.

¶Применение

ВСПМ широко используется в различных областях науки и техники:

• Анализ вибраций и акустика: выявление источников шума и вибраций, определение путей распространения колебаний в механических конструкциях.
• Радиотехника и связь: оценка характеристик каналов передачи, анализ интермодуляционных искажений, синхронизация сигналов.
• Геофизика и сейсмология: изучение распространения сейсмических волн, определение скорости и направления движения волновых фронтов.
• Биомедицинская инженерия: анализ электроэнцефалограмм (ЭЭГ) и электрокардиограмм (ЭКГ) для исследования функциональных связей между областями мозга или отделами сердца.
• Машиностроение: диагностика состояния подшипников и зубчатых передач по взаимному спектру вибраций корпуса и вала.

¶Связь с передаточными функциями

В теории линейных систем ВСПМ играет ключевую роль при идентификации динамических характеристик. Если \( y(t) \) — выход линейной стационарной системы с импульсной характеристикой \( h(t) \), а \( x(t) \) — вход, то справедливо соотношение:

\[ S_{xy}(f) = H(f) S_{xx}(f) \]

где \( H(f) \) — частотная характеристика системы. Это позволяет оценивать амплитудные и фазовые искажения, вносимые системой, без проведения прямых измерений передаточной функции.

¶Ограничения

Применение ВСПМ требует выполнения условий стационарности и эргодичности исследуемых процессов. Для нестационарных сигналов используются модифицированные подходы, например, на основе вейвлет-преобразования или короткооконного преобразования Фурье. Кроме того, при малом числе усреднений оценка ВСПМ может быть смещённой, что требует аккуратного выбора параметров анализа.

¶Источники

• Бендат Дж., Пирсол А. Прикладной анализ случайных данных. — М.: Мир, 1989.
• Марпл-мл. С. Л. Цифровой спектральный анализ и его приложения. — М.: Мир, 1990.
• Отнес Р., Эноксон Л. Прикладной анализ временных рядов. — М.: Мир, 1982.