Глава 1 математические модели сигналов 4 страница

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

Вероятностные характеристики случайного процесса. В соответствии с определением случайный процесс U(t) может быть описан системой N обычно зависимых случайных величин U1 = U(t1), ..., Ui= U(ti), ..., UN = U(tN), взятых в различные моменты времени t1...ti...tN. При неограниченном увеличении числа N такая система эквивалентна рассматриваемому случайному процессу U(t).

Исчерпывающей характеристикой указанной системы является N-мерная плотность вероятности pN(U1, ..., Ui, ..., UN; t1, ..., tN). Она позволяет вычислить вероятность РN реализации, значения которой в моменты времени t1,t2, ...,tN будут находиться соответственно в интервалах (u1, u1+Δu1), ..., (ui, ui+ глава 1 математические модели сигналов 4 страница - student2.ru ui), ..., (uN, uN + глава 1 математические модели сигналов 4 страница - student2.ru uN), где ui(1 глава 1 математические модели сигналов 4 страница - student2.ru ) —значение, принимаемое случайной величиной Ui, (рис. 1.12).

глава 1 математические модели сигналов 4 страница - student2.ru

Если Δui, выбраны достаточно малыми, то справедливо соотношение

глава 1 математические модели сигналов 4 страница - student2.ru

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

На практике в таком подробном описании нет необходимости. Обычно ограничиваются одно- или двумерной плотностью вероятности.

Одномерная плотность вероятности p1(U1; t1) случайного процесса U(t) характеризует распределение одной случайной величины U1, взятой в произвольный момент времени t1. В ней не находит отражения зависимость случайных величин в различные моменты времени.

Двумерная плотность вероятности p2 = p2(U1, U2; t1, t2) позволяет определить вероятность совместной реализации любых двух значений случайных величин U1 и U2 в произвольные моменты времени t1 и t2 и, следовательно, оценить динамику развития процесса. Одномерную плотность вероятности случайного процесса U(t) можно получить из двумерной плотности, воспользовавшись соотношением

глава 1 математические модели сигналов 4 страница - student2.ru

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

Математическим ожиданием случайного процесса U(t) называют неслучайную функцию времени mu(t1), которая при любом аргументе t1 равна среднему значению случайной величины U(t\) по всему множеству возможных реализаций:

глава 1 математические модели сигналов 4 страница - student2.ru

Степень разброса случайных значений процесса U(t1) от своего среднего значения mu(t1) для каждого t1 характеризуется дисперсией Du(t1):

глава 1 математические модели сигналов 4 страница - student2.ru

где глава 1 математические модели сигналов 4 страница - student2.ru (t1) = U(t1) — mu(t1) — центрированная случайная величина.

Дисперсия Du(t1) в каждый момент времени t1 равна квадрату среднеквадратического отклонения su(t1):

глава 1 математические модели сигналов 4 страница - student2.ru

Случайные процессы могут иметь одинаковые математические ожидания и дисперсии (рис. 1.13, а, б), однако резко различаться по быстроте изменений своих значений во времени.

глава 1 математические модели сигналов 4 страница - student2.ru Для оценки степени статистической зависимости мгновенных значений процесса U(t) в произвольные моменты времени t1 и t2 используется неслучайная функция аргументов Ru(t1,t2), называемая автокорреляционной или просто корреляционной функцией.

При конкретных аргументах t1 и t2 она равна корреляционному моменту значений процесса U(t1) и U(t2):

глава 1 математические модели сигналов 4 страница - student2.ru

Через двумерную плотность вероятности выражение (1.71) представляется в виде

глава 1 математические модели сигналов 4 страница - student2.ru

В силу симметричности этой формулы относительно аргументов справедливо равенство

глава 1 математические модели сигналов 4 страница - student2.ru

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

глава 1 математические модели сигналов 4 страница - student2.ru

Из сравнения (1.69) и (1.70) следует, что при произвольном t1 = t2 автокорреляционная функция вырождается в дисперсию:

глава 1 математические модели сигналов 4 страница - student2.ru

а нормированная функция автокорреляции равна единице:

глава 1 математические модели сигналов 4 страница - student2.ru

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

Аналогично устанавливается мера связи между двумя случайными процессами U(t) и V(t). Она называется функцией взаимной корреляции:

глава 1 математические модели сигналов 4 страница - student2.ru

§ 1.9. СТАЦИОНАРНЫЕ И ЭРГОДИЧЕСКИЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ

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

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

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

Случайный процесс называют стационарным в узком смысле, если выражения для плотностей вероятности не зависят от начала отсчета времени, т. е. справедливо соотношение

глава 1 математические модели сигналов 4 страница - student2.ru

где U глава 1 математические модели сигналов 4 страница - student2.ru —случайная величина, отражающая значение процесса в момент времени t = ti + τ (τ — произвольное число).

Иначе говоря, стационарность процесса предполагает его существование и статистическую однородность во всем диапазоне времени от - глава 1 математические модели сигналов 4 страница - student2.ru до + глава 1 математические модели сигналов 4 страница - student2.ru .

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

При решении многих технических задач идут на дальнейшее упрощение модели, рассматривая случайный процесс стационарным в широком смысле. Процесс U(t) принято называть стационарным в широком смысле, если выполняется условие постоянства математического ожидания и дисперсии, а корреляционная функция не зависит от начала отсчета времени и является функцией только одного аргумента τ = t2 - t1, т.е.

глава 1 математические модели сигналов 4 страница - student2.ru

Так как условие постоянства дисперсии является частным случаем требования к корреляционной функции при τ = 0:

глава 1 математические модели сигналов 4 страница - student2.ru

то выполнения соотношений (1.79) и (1.81) достаточно, чтобы рассматривать случайный процесс U(t) как стационарный.

Всякий стационарный случайный процесс является стационарным в широком смысле. В дальнейшем, если это не оговорено особо, стационарность будем рассматривать в широком смысле.

Случайные процессы, наблюдаемые в устойчиво работающих реальных системах, имеют конечное время корреляции. Поэтому для стационарных процессов, представляющих практический интерес, справедливо соотношение

глава 1 математические модели сигналов 4 страница - student2.ru

Если для случайного процесса равенства (1.79), (1.81) не выдерживаются, но на интересующем нас интервале времени изменением указанных параметров можно пренебречь, его называют квазистационарным.

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

Следовательно, для стационарных эргодических процессов справедливы соотношения

глава 1 математические модели сигналов 4 страница - student2.ru

где u(t) — конкретная реализация случайного процесса U(t).

Результаты исследования случайных процессов в их временном представлении, т. е. с использованием формул (1.83) и (1.85), лежат в основе корреляционной теории сигналов.

Для облегчения практического определения корреляционных функций в соответствии с (1.85) серийно выпускаются специальные вычислительные устройства — коррелометры (корреляторы).

§ 1.10. СПЕКТРАЛЬНОЕ ПРЕДСТАВЛЕНИЕ СЛУЧАЙНЫХ СИГНАЛОВ

В § 1.2 была показана эффективность представления детерминированных сигналов совокупностью элементарных базисных сигналов для облегчения анализа прохождения их через линейные системы. Аналогичный подход может быть использован и в случае сигналов, описываемых случайными процессами [21].

Рассмотрим случайный процесс U(t), имеющий математическое ожидание mu(t). Соответствующий центрированный случайный процесс глава 1 математические модели сигналов 4 страница - student2.ru (t) характеризуется в любой момент времени t1 центрированной случайной величиной глава 1 математические модели сигналов 4 страница - student2.ru (t1):

глава 1 математические модели сигналов 4 страница - student2.ru

Центрированный случайный процесс глава 1 математические модели сигналов 4 страница - student2.ru (t) можно, как и ранее [см. (1.1)], выразить в виде конечной или бесконечной суммы ортогональных составляющих, каждая из которых представляет собой неслучайную базисную функцию jk(t) с коэффициентом Ck, являющимся случайной величиной. В результате имеем разложение центрированного случайного процесса глава 1 математические модели сигналов 4 страница - student2.ru глава 1 математические модели сигналов 4 страница - student2.ru (t):

глава 1 математические модели сигналов 4 страница - student2.ru глава 1 математические модели сигналов 4 страница - student2.ru

Случайные величины Сk называются коэффициентами разложения. В общем случае они статистически зависимы, и эта связь задается матрицей коэффициентов корреляции глава 1 математические модели сигналов 4 страница - student2.ru . Математические ожидания коэффициентов разложения равны нулю. Неслучайные базисные функции принято называть координатными функциями.

Для конкретной реализации коэффициенты разложения являются действительными величинами и определяются по формуле (1.7).

Предположив, что

глава 1 математические модели сигналов 4 страница - student2.ru

детерминированную функцию mu(f) в (1.86) на интервале —T<t<.T также можно разложить по функциям φk(t), представив в виде

глава 1 математические модели сигналов 4 страница - student2.ru

Подставляя (1.87 а) и (1.876) в (1.86) для случайного процесса U(t) с отличным от нуля средним, получим

глава 1 математические модели сигналов 4 страница - student2.ru

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

детерминированных функций [mu(t), глава 1 математические модели сигналов 4 страница - student2.ru jk(t)], а коэффициенты разложения, являющиеся случайными величинами, остаются неизменными.

Чтобы определить требования к координатным функциям, рассмотрим корреляционную функцию процесса глава 1 математические модели сигналов 4 страница - student2.ru (t), заданную разложением

глава 1 математические модели сигналов 4 страница - student2.ru

Так как

глава 1 математические модели сигналов 4 страница - student2.ru

то

глава 1 математические модели сигналов 4 страница - student2.ru

Соотношение (1.88) становится значительно проще, если коэффициенты {Ck} некоррелированы (Rkl = 0 при k глава 1 математические модели сигналов 4 страница - student2.ru l, Rkl = 1 при k = l):

глава 1 математические модели сигналов 4 страница - student2.ru

В частности, при t1 = t2 = t получим дисперсию случайного процесса U(t):

глава 1 математические модели сигналов 4 страница - student2.ru

Поэтому целесообразно выбирать такие координатные функции, которые обеспечивают некоррелированность случайных величин {Сk}. Разложение (1.87), удовлетворяющее этому условию, называют каноническим разложением.

Доказано [21], что по известному каноническому разложению корреляционной функции случайного процесса можно записать каноническое разложение самого случайного процесса с теми же координатными функциями, причем дисперсии коэффициентов этого разложения будут равны дисперсиям коэффициентов разложения корреляционной функции.

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

В каноническом разложении (1.87) этот спектр является дискретным (линейчатым) и может содержать как конечное, так и бесконечное число членов (линий).

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

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

§ 1.11. ЧАСТОТНОЕ ПРЕДСТАВЛЕНИЕ СТАЦИОНАРНЫХ

СЛУЧАЙНЫХ СИГНАЛОВ

Дискретные спектры. Корреляционную функцию Ru(t) (рис. 1.14) стационарного случайного процесса, заданного на конечном интервале времени [-Т, Т], можно разложить в ряд Фурье (1.15), условно считая ее периодически продолжающейся с периодом 4T (при —T<.t1, t2<T, -2Τ<τ<2Τ):

глава 1 математические модели сигналов 4 страница - student2.ru

где

глава 1 математические модели сигналов 4 страница - student2.ru

глава 1 математические модели сигналов 4 страница - student2.ru

Учитывая, что Ru(t) является четной функцией, имеем

глава 1 математические модели сигналов 4 страница - student2.ru

Положив τ = t1 - t2, находим

глава 1 математические модели сигналов 4 страница - student2.ru

что согласно (1.89) представляет собой каноническое разложение корреляционной функции. По нему, как было указано ранее, получаем каноническое разложение случайного процесса:

глава 1 математические модели сигналов 4 страница - student2.ru

причем

глава 1 математические модели сигналов 4 страница - student2.ru

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

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

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

глава 1 математические модели сигналов 4 страница - student2.ru

где

глава 1 математические модели сигналов 4 страница - student2.ru

На спектральной диаграмме такого процесса каждой гармонике ставится в соответствие вертикальный отрезок, длина которого пропорциональна дисперсии ее амплитуды, а расположение на оси абсцисс отвечает частоте (рис. 1.15).

глава 1 математические модели сигналов 4 страница - student2.ru

Чтобы получить описание стационарного случайного процесса в точном смысле, т. е. справедливое для любого момента времени на бесконечном интервале - глава 1 математические модели сигналов 4 страница - student2.ru <t< глава 1 математические модели сигналов 4 страница - student2.ru , необходимо перейти к интегральному каноническому разложению.

Непрерывные спектры. Интегральное каноническое разложение для корреляционной функции получим из формулы (1.91) путем предельного перехода при Т глава 1 математические модели сигналов 4 страница - student2.ru . Увеличение интервала времени, на котором наблюдается случайный процесс, сопровождается уменьшением значений дисперсий, что следует из (1.92), а также сокращением расстояния между спектральными линиями, поскольку

глава 1 математические модели сигналов 4 страница - student2.ru

При достаточно большом, но конечном Т можно записать выражение для средней плотности распределения дисперсии по частоте:

S глава 1 математические модели сигналов 4 страница - student2.ru (wk) = Dk/(Δω) = 2DkT/ глава 1 математические модели сигналов 4 страница - student2.ru (k = 0,±l, ±2, ..) (1.98)

где S глава 1 математические модели сигналов 4 страница - student2.ru (wk) - средняя плотность дисперсии на участке, прилегающем к частоте ωk.

Теперь можно преобразовать формулы (1.94) и (1.98) к виду

глава 1 математические модели сигналов 4 страница - student2.ru

Переходя к пределу при Т глава 1 математические модели сигналов 4 страница - student2.ru , получаем

глава 1 математические модели сигналов 4 страница - student2.ru

где

глава 1 математические модели сигналов 4 страница - student2.ru

Так как величина S глава 1 математические модели сигналов 4 страница - student2.ruk)Δω являлась не только дисперсией Dk коэффициента разложения корреляционной функции Ru(t), но и дисперсией D[Ck] коэффициента разложения случайного процесса U(t), то величина Suu(w)dw, полученная в результате предельного перехода при Т глава 1 математические модели сигналов 4 страница - student2.ru , представляет собой дисперсию, приходящуюся на спектральные составляющие стационарного случайного процесса, занимающие бесконечно малый интервал частот (ω, ω + dw). Функцию Suu(w), характеризующую распределение дисперсии случайного процесса по частотам, называют спектральной плотностью стационарного случайного процесса U(t).

Выражение для интегрального канонического разложения корреляционной функции Ru(t) найдем, положив в формуле (1.101) τ = t1 - t2:

глава 1 математические модели сигналов 4 страница - student2.ru

Обозначив G глава 1 математические модели сигналов 4 страница - student2.ru (w) = Сk/( глава 1 математические модели сигналов 4 страница - student2.ru w) и повторив процедуру предельного перехода при T глава 1 математические модели сигналов 4 страница - student2.ru для соотношения (1.95), получим каноническое разложение стационарной случайной функции U(t):

глава 1 математические модели сигналов 4 страница - student2.ru

где дисперсией случайной функции G(w)dw является функция Suu(w)dw.

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

Основные свойства спектральной плотности. Отметим, что в формулах (1.101) и (1.102) Suu(w) определена как для положительных, так и для отрицательных частот. Перейдем к одностороннему спектру, ограничиваясь только положительными частотами. Воспользовавшись формулой Эйлера, представим соотношение (1.102) состоящим из двух слагаемых:

глава 1 математические модели сигналов 4 страница - student2.ru

В силу четности функции Ru(t) второе слагаемое равно нулю, а первое можно преобразовать к виду

глава 1 математические модели сигналов 4 страница - student2.ru

Из (1.105) следует, что Suu(w) является действительной и четной функцией, т. е.

глава 1 математические модели сигналов 4 страница - student2.ru

Это позволяет ограничиться положительными частотами и в (1.101):

глава 1 математические модели сигналов 4 страница - student2.ru

Соотношения (1 101) и (1.102), а также (1.105) и (1.107) являются парами интегрального преобразования Фурье, причем (1.105) и (1.107) для случая четной функции. Поэтому корреляционная функция и спектральная плотность подчинены закономерности: чем протяженнее кривая Suu(ω), тем уже корреляционная функция Ru(t) (тем меньше время корреляции), и наоборот.

Площадь, ограниченная непрерывной кривой Suu(w) на спектральной диаграмме, очевидно, должна равняться дисперсии Du случайного процесса U(t). Действительно, положив в формуле (1.107) τ = 0, получим

глава 1 математические модели сигналов 4 страница - student2.ru

Подразумевая под случайным процессом U(t) напряжение, Du можно рассматривать как среднюю мощность, выделяемую этим напряжением на резисторе с сопротивлением в 1 Ом:

глава 1 математические модели сигналов 4 страница - student2.ru

Следовательно, величина

глава 1 математические модели сигналов 4 страница - student2.ru

представляет собой долю средней мощности, выделяемой составляющими спектра, относящимися к интервалу частот (ω, ω + dw).

В связи с этим спектральную плотность Suu(w) называют еще спектральной плотностью мощности, а также энергетическим спектром стационарного случайного процесса, поскольку Suu(w) имеет размерность энергии.

Спектральная плотность мощности случайного процесса является средней характеристикой множества реализаций. Ее можно получить и путем усреднения спектральной мощности реализации Ρk(ω) (1.62) по множеству реализаций.

Рассмотрим с этой целью одну реализацию u глава 1 математические модели сигналов 4 страница - student2.ru (t) стационарного случайного процесса U(t) сначала на ограниченном интервале времени -T<t<.T. Для нее можно записать преобразование Фурье:

глава 1 математические модели сигналов 4 страница - student2.ru

В соответствии с (1.63) спектральная плотность мощности этой реализации

глава 1 математические модели сигналов 4 страница - student2.ru

Найдем среднее значение Ρ глава 1 математические модели сигналов 4 страница - student2.ru (ω) по множеству реализации k. Имеем

глава 1 математические модели сигналов 4 страница - student2.ru

или

глава 1 математические модели сигналов 4 страница - student2.ru

Так как мы предполагаем, что случайный процесс U(t) стационарный, то

глава 1 математические модели сигналов 4 страница - student2.ru

где t1 - t2 = τ.

При выполнении условия (1.114) для выражения (1.113) существует предел при T глава 1 математические модели сигналов 4 страница - student2.ru :

глава 1 математические модели сигналов 4 страница - student2.ru

что и требовалось показать.

Пример 1.7. У центрированного стационарного случайного процесса спектральная плотность постоянна. Рассмотреть особенности такого процесса.

Пусть спектральная плотность Suu(ω) ограничена определенной полосой частот (рис. 1.16, а):

глава 1 математические модели сигналов 4 страница - student2.ru

В соответствии с (1.107) найдем автокорреляционную функцию процесса U(t):

глава 1 математические модели сигналов 4 страница - student2.ru

Вид функции Ru(t) приведен на рис. 1.16,6. Значение ее при τ = 0 равно дисперсии, а следовательно, средней мощности рассматриваемого процесса:

глава 1 математические модели сигналов 4 страница - student2.ru

Будем теперь расширять полосу частот, занимаемую энергетическим спектром (рис. 1.17, а). Интервал времени, на котором наблюдается существенная корреляционная связь значений процесса, при этом уменьшается, а дисперсия Du возрастает.

При w0 глава 1 математические модели сигналов 4 страница - student2.ru дисперсия становится безграничной, а корреляционная функция принимает вид дельта-функции (рис. 1.17,6).

глава 1 математические модели сигналов 4 страница - student2.ru

Идеализированный случайный процесс, энергетический спектр которого безграничен и равномерен, известен как «белый шум». Такое название возникло по аналогии с белым светом, имеющим равномерный и неограниченный спектр интенсивности. Основная особенность процесса в том, что его значения в любые два сколь угодно близкие моменты времени некоррелированы. Создать белый шум принципиально невозможно, так как реальные источники сигналов всегда имеют ограниченную мощность. Тем не менее, понятие «белый шум» нашло широкое применение в информационной технике. Такая модель может быть принята, например, для сигналов (шумов), имеющих равномерный энергетический спектр в пределах полосы пропускания входного блока системы, в которой они рассматриваются.

Наши рекомендации