2.4. ОГИБАЮЩАЯ И ФАЗА СЛУЧАЙНОГО ПРОЦЕССА

Квазигармоническое представление возможно для любого сигнала, однако наиболее полезно оно при изучении относительно узкополосных сигналов, у которых ширина спектра значительно меньше средней частоты Для таких сигналов огибающая является низкочастотным процессом, т. е. изменяется медленно по сравнению с изменениями Это же верно и для мгновенной начальной фазы и для ее производной при условии, что частота выбрана вблизи средней круговой частоты спектра При наблюдении на осциллографе такой процесс имеет вид колебания, промодулированного по амплитуде и фазе.

Выражение (2.67) можно переписать также следующим образом:

где комплексная огибающая сигнала. В отличие от действительной огибающей комплексная огибающая не определяется однозначно сигналом а зависит также от выбора произвольной частоты Из очевидного равенства

следует, что комплексную огибающую можно рассматривать как результат сдвига спектра аналитического сигнала вниз на

Введем обозначения для действительной и мнимой частей комплексной огибающей:

С помощью этих функций можно получить еще одно полезное представление сигнала

Таким образом, сигнал разложен на две квадратурные составляющие:

Функции зависят от выбора частоты Для узкополосного сигнала эти функции, при надлежащем выборе являются, как и медленно меняющимися (низкочастотными). Их, однако, не следует называть огибающими квадратурных составляющих, поскольку, в отличие от они принимают не только положительные, по и отрицательные значения.

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

Рассмотрим вначале центрированный стационарный гауссовский процесс Поскольку преобразование Гильберта линейно, то, учитывая инвариантность нормального распределения при любых линейных операциях, можно утверждать, что также гауссовский процесс.

Па основании того, что отличается от только сдвигом фаз всех составляющих, вспоминая, что спектральная плотность мощности не зависит от фазовых соотношений,

приходим к выводу, что Из этого на основании теоремы Хинчина-Винера следует равенство корреляционных функций а отсюда и равенство дисперсий В совпадающие моменты времени между собой не коррелированы а поскольку они гауссовские, то и независимы. Комплексная огибающая полученная из посредством линейной операций (2.70а), обладает теми же свойствами: также являются центрированными гауссовскими процессами, независимыми в совпадающие моменты времени. Дисперсия аналитического сигнала

Такова же, как легко проверить, и дисперсия комплексной огибающей Нетрудно также убедиться в том, что

Совместную плотность вероятности в одном сечении можно записать на основании (2.12) так:

Здесь дисперсия, одинаковая для X и для

Процессы и , полученные из с помощью нелинейных операций, стационарны, но конечно не являются гауссовскими. Найдем совместную функцию распределения в некотором сечении т. е. функцию

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

Для нахождения функции распределения можно проинтегрировать

Рис. 2.10. К вычислению распределения вероятностей огибающей и фазы центрированного нормального процесса

(2.73) по области Переходя при этом к полярной системе координат по известным формулам

получим

Таким образом, двумерная функция распределения оказалась произведением двух одномерных функций:

Это значит, что огибающая и фаза в одном сечении независимы.

Плотность распределения огибающей найдем как производную функции распределения:

Это известное распределение Рэлея, график которого показан на рис. 2.11 вместе с графиком гауссовского распределения мгновенных значений Наиболее вероятное значение равно . Математическое ожидание Математическое ожидание квадрата огибающей откуда дисперсия

Рис. 2.11. Плотность вероятности мгновенного значения и огибающей нормального центрированного процесса

Рис. 2.12. Плотность вероятности фазы нормального центрированного процесса

Для плотности вероятности фазы имеем (рис. 2.12)

Следовательно, фаза стационарного гауссовского центрированного процесса распределена равномерно на интервале протяженностью

Мгновенная начальная фаза рассматриваемого процесса также равномерно распределена на интервале протяженностью (если ее определять по модулю Действительно, можно отыскать совместную функцию распределений исходя из плотности процессов Повторяя буквально те же рассуждения и выкладки, только с заменой на найдем, как и следовало ожидать, плотность (2.76) для и равномерную плотность вероятностей для

Воспользуемся этим методом для рассмотрения несколько более сложного случая суммы стационарного центрированного гауссовского процесса и гармонического сигнала:

где центрированный стационарный гауссовский процесс; — некоторые постоянные. Этот случай нередко встречается на практике, когда представляет регулярный сигнал, флуктуациовный шум. Заметим, что нестационарный процесс, так как его математическое ожидание и зависит от времени.

Представив в форме (2.72) и выбирая можем записать:

где

Отсюда видно, что являются стационарными гауссовскими процессами с одинаковыми дисперсиями взаимно некоррелированными в совпадающие моменты времени. При этом

Для нахождения функции распределения

рассмотрим совместную плотность распределения вероятностей в некотором сечении

Функция распределения (2.80) равна интегралу от (2.81) по площади, определяемой (2.80). Переходя к полярным координатам находим

Не будем вычислять этот интеграл. Заметим лишь, что в этом случае не выражается произведением двух функций распределения т. е. не являются независимыми. В дальнейшем нам потребуется только распределение вероятностей Для нахождения нужно в (2.82) положить так как по определению и вероятность того, что равна 1, поэтому при совпадает с Итак,

Из теории бесселевых функций известно, что

где модифицированная функция Бесселя нулевого порядка:

При эта функция равна единице и монотонно возрастает с х. Таким образом,

и плотность вероятностей для суммы гауссовского шума и гармонического сигнала

Распределение вероятностей, выраженное формулой (2.86), называется обобщенным распределением Рэлея, или распределением Рэлея-Райса. Очевидно, обычное распределение Рэлея (2.76) является частным случаем распределения (2.86) при

Заметим без дсжазательства, что и в более общем случае суммы гауссовского шума с синусоидальным сигналом

где любая детерминированная или случайная начальная фаза, огибающая подчиняется распределению (2.86). В частности, если равномерно распределена на интервале протяженностью то начальная фаза также равномерно

Рис. 2.13. Плотность вероятности огибающей суммы гармоиического сигнала и гауссовского шума

распределена на том же интервале и не зависит от Процесс о этом случае сташюнарен.

На рис. 2.13 показаны графини плотности (2.86) при разных значениях и. Легко видеть, что чем больше т. е. чем больше отношение амплитуды гармонического сигнала к средней квадратичной величине шума, тем больше и наиболее вероятное значение