Спектральное представление стационарных случайных процессов

Реализации стационарного СП x(t,w) по определению однородны на интервале tÎ(-¥,¥), и, следовательно, имеют неограниченную энергию Спектральное представление стационарных случайных процессов - student2.ru . Ограниченной является средняя по времени мощность реализации x(t,w)

Спектральное представление стационарных случайных процессов - student2.ru .(3.5.1)

Для таких процессов можно определить спектр реализации gx(f,w)=F{x(t,w)}. Расчет корреляционной функции Спектральное представление стационарных случайных процессов - student2.ru показывает (см. пример 3.5.4), что случайный спектр gx(f,w) является d - коррелированной нестационарной функцией частоты с неограниченной дисперсией, т.е. принадлежит к классу обобщенных функций.

В практических задачах используются усеченные реализации стационарных СП, отличные от нуля на конечном интервале tÎ(-T,T) (или 0£t£T)

Спектральное представление стационарных случайных процессов - student2.ru (3.5.2)

Спектр усеченной реализации

Спектральное представление стационарных случайных процессов - student2.ru (3.5.3)

является обычной функцией и может быть использован, например, для определения реакции линейной системы на воздействие данной реализации xт(t,w). Однако спектр Спектральное представление стационарных случайных процессов - student2.ru не является характеристикой процесса X(t)={xт(t,w), wÎW} в целом, так как для различных реализаций он получается различным. Кроме того, он изменяется при изменении интервала T.

Наиболее целесообразным, оказалось, характеризовать спектральные свойства стационарного СП распределением его мощности (1) на оси частот f(-¥,¥). Определим эту характеристику, для чего выразим среднюю мощность Px(w) реализации x(t,w) стационарного, неэргодического, СП X(t) через спектр Спектральное представление стационарных случайных процессов - student2.ru усеченной реализации

Спектральное представление стационарных случайных процессов - student2.ru . (3.5.4)

Изменив в сходящемся интеграле (4) порядок интегрирования и перехода к пределу, введем обозначение

Спектральное представление стационарных случайных процессов - student2.ru .(3.5.5)

В результате получим

Спектральное представление стационарных случайных процессов - student2.ru . (5.5.6)

Функция Wx(f,w) имеет ясный физический смысл - это спектральная плотность мощности px(w) конкретной реализации x(t,w), т.е. средняя мощность реализации, приходящаяся на 1Гц полосы частот на данной частоте f. Часть мощности реализации (1), сосредоточенной в спектре в произвольной полосе fÎ(f1, f2 ) определяется интегралом

Спектральное представление стационарных случайных процессов - student2.ru .

Спектральная плотность мощности Wx(f,w) является характеристикой реализации и для различных реализаций она различна (рис.3.8 штриховые линии). Заметим, что для эргодических СП средняя мощность для всех реализаций одинаковая Px(w)=Px . Однако спектральная плотность мощности Wx(f,w) при этом для различных реализаций остается различной. Иначе говоря, для эргодических СП эргодическим свойством обладает не спектр Wx(f,w), а интеграл (6) от этой функции.

Спектральное представление стационарных случайных процессов - student2.ru

Рис.3.8

Для получения спектральной плотности мощности (спектра мощности) Wx(f) процесса X(t) в целом (а не отдельной реализации) необходимо произвести усреднение Wx(f,w) по ансамблю реализаций (рис.3.8 сплошная линия)

Спектральное представление стационарных случайных процессов - student2.ru .(3.5.7)

Ученые А.Я.Хинчин и Н.Винер независимо установили, что спектральная плотность мощности Wx(f) стационарного СП X(t) и его основная характеристика - корреляционная функция Rx(t) взаимно однозначно связаны между собой преобразованием Фурье (теорема Хинчина-Винера). Для доказательства этой теоремы преобразуем, принимая во внимание (3), выражение (7) для спектральной плотности Wx(f), полагая, что x(t,w) – реализации центрированного СП,

Спектральное представление стационарных случайных процессов - student2.ru . (3.5.8)

Производим замену переменных (якобиан преобразования равен 1)

Спектральное представление стационарных случайных процессов - student2.ru (3.5.9)

и соответствующее преобразование области интегрирования (t1,t2)Î(-T,T) (рис.3.9).

Спектральное представление стационарных случайных процессов - student2.ru

Рис. 3.9

На рисунке зафиксированы характерные точки областей 1, 2, 3, 4. Интегрирование ведется сначала по заштрихованным полоскам (по t) - внутренний интеграл, затем по t - внешний интеграл

Спектральное представление стационарных случайных процессов - student2.ru (3.5.10)

После интегрирования по t (10) можно представить одним интегралом

Спектральное представление стационарных случайных процессов - student2.ru .(3.5.11)

Функция Y(t,T)=[1-|t|/(2T)], |t|£2T, входящая в (11), показана на рис.3.10. При T®¥ Y(t,T)®1. В результате получаем утверждение теоремы. Корреляционная функция Rx(t) - стационарного СП и его спектральная плотность (мощности) Wx(f) представляют пару преобразований Фурье

Спектральное представление стационарных случайных процессов - student2.ru . (3.5.12)

Спектральное представление стационарных случайных процессов - student2.ru . (3.5.13)

Спектральное представление стационарных случайных процессов - student2.ru

Рис. 3.10

В соответствии с определением (7) спектральная плотность стационарного СП Wx(f) является положительной четной функцией. При использовании экспоненциального базиса спектр Wx(f) определяется на всей оси частот fÎ(-f, f) и называется двусторонним спектром. Физический спектр, непосредственно измеряемый, определяется в базисе гармонических функций только на положительной полуоси частот fÎ(0, f) и называется односторонним спектром Спектральное представление стационарных случайных процессов - student2.ru

Спектральное представление стационарных случайных процессов - student2.ru . (3.5.14)

Для нецентрированных стационарных СП в (8) и далее корреляционная функция Спектральное представление стационарных случайных процессов - student2.ru , что приводит к выражению для спектральной плотности нецентрированного процесса Спектральное представление стационарных случайных процессов - student2.ru

Спектральное представление стационарных случайных процессов - student2.ru , (3.5.15)

которая отличается от (12) наличием спектральной линии на нулевой частоте с мощностью Спектральное представление стационарных случайных процессов - student2.ru .

Таким образом, в рамках корреляционной теории в качестве основной характеристики стационарного СП X(t) вместо корреляционной функции Rx(t) может использоваться спектральная плотность Wx(f). Если процесс, кроме того, гауссовский, то функция Wx(f) так же как и Rx(t) полностью определяет центрированный стационарный СП. Для нецентрированных процессов требуется еще задание mx.

Пример 3.5.1. Определить корреляционную функцию стационарного СП с равномерным ограниченным спектром мощности

Спектральное представление стационарных случайных процессов - student2.ru .(3.5.16)

N0 - односторонняя (N0/2 – двусторонняя) спектральная плотность мощности, являющаяся основной характеристикой процесса с равномерным спектром.

Полная мощность (дисперсия) процесса

Спектральное представление стационарных случайных процессов - student2.ru .(3.5.17)

Корреляционная функция

Спектральное представление стационарных случайных процессов - student2.ru (3.5.18)

совпадает (с точностью до коэффициента N0F) с нулевой функцией отсчетов g0(t). Корреляционная функция (18) обращается в нуль в точках tk=k/(2F), k=±1,±2,... Это значит, что сечения процесса X(t) и X[t+k/(2F)] со спектральной плотностью (16), отстоящие друг от друга на интервал, кратный величине 1/(2F), не коррелированны (для гауссовских СП - не зависимы). Время корреляции процесса величина порядка 1/(2F).

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

Спектральное представление стационарных случайных процессов - student2.ru .(3.5.19)

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

Корреляционная функция белого шума

Спектральное представление стационарных случайных процессов - student2.ru (3.5.20)

имеет вид d - функции. Это значит, что два сечения процесса X(t1) и X(t2) в любые несовпадающие моменты отсчета t1¹t2 не коррелированны (для гауссовских процессов не зависимы). Поэтому белый шум называют также процессом с некоррелированными значениями или некоррелированным СП. Соответственно белый шум можно рассматривать как идеализацию процесса с очень малым временем корреляции. В статистической теории РТС широко используется модель в виде гауссовского СП с характеристиками (19), (20), называемая гауссовским белым шумом (ГБШ).

Пример 3.5.3. Определить корреляционную функцию стационарного СП с гауссовским (колокольным) спектром мощности

Спектральное представление стационарных случайных процессов - student2.ru , (3.5.21)

где 2F - эффективная ширина спектра, определяемая как основание прямоугольника П(f/2F) равновеликого (с равной площадью) экспоненте exp[-pf2/(2F)2]. Уровень отсчета полосы 2F равен exp(-p/4)=0,46.

Корреляционная функция (повторяется расчет, приведенный в примере 2.2.3)

Спектральное представление стационарных случайных процессов - student2.ru (3.5.22)

получается также гауссовской. Эффективная длительность корреляционной функции (удвоенное время корреляции tk) равна 1/(2F).

Пример 3.5.4. Определить основные характеристики случайной функции gx(f,w) - спектра реализации стационарного СП x(t,w). Математическое ожидание Спектральное представление стационарных случайных процессов - student2.ru gx(f,w) Спектральное представление стационарных случайных процессов - student2.ru =0, а корреляционная функция

Спектральное представление стационарных случайных процессов - student2.ru

после подстановки t1 =t, t2 =t-t, получается равной

Спектральное представление стационарных случайных процессов - student2.ru (3.5.23)

Следовательно, спектр gx(f,w) представляет собой нестационарную случайную функцию (корреляционная функция зависит от двух переменных f1 и f2), которая бесконечно быстро флюктуирует и имеет бесконечно большую дисперсию.

Пример 3.5.5. Определить корреляционную функцию Rx(t) стационарного узкополосного СП x(t) с равномерным спектром Wx(f) (рис.3.11а)

Спектральное представление стационарных случайных процессов - student2.ru . (3.5.24)

Спектральное представление стационарных случайных процессов - student2.ru

Принимая во внимание, что (sina–sinb)=2sin[1/2(a–b)]cos[1/2(a+b)], получаем (рис.3.11,б)

Спектральное представление стационарных случайных процессов - student2.ru

Рис. 3.11

Математические модели помех

Наблюдение, поступающее на вход приемного устройства РТС, представляет собой сумму полезного сигнала s(t,l,a), модель которого уже была изучена, и аддитивной (маскирующей) помехи n(t)

Спектральное представление стационарных случайных процессов - student2.ru .

Ближайшей задачей теперь является обоснование математической модели аддитивных помех n(t). В аддитивную помеху входят, во-первых, внутренние шумы аппаратуры и помехи природного происхождения: радиоизлучение Земли и атмосферы, радиоизлучение внеземного происхождения (Солнца, планет, космического пространства). Эти помехи являются гауссовскими, однородными на интервале наблюдения tÎ(0,T) и имеют сравнительно с сигналом широкий спектр.

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

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

Второй моделью является общая гауссовская помеха, нестационарная или стационарная со спектром неравномерным в пределах полосы частот, занятой сигналом. Эта модель описывает большинство помех второй группы. Иногда изучается оптимизация приема на фоне негауссовских помех. Аналитический аппарат при этом значительно усложняется. В связи с этим отметим, что в американской литературе [4] утверждается: "если система хорошо спроектирована для работы на фоне гауссовских помех, то она будет хорошо работать также на фоне негауссовских помех".

Рассмотрение начнем с модели основной гауссовской помехи. В дальнейшем, следуя практике, принятой в технической литературе, будем обозначать СП и его отдельные реализации одними и теми же символами: u(t), n(t), x(t) и др.

Рис. 3.12

На рис.3.12 заштрихована область частот fÎ(f0-Fc, f0+Fc), занятая на положительной полуоси fÎ(0,¥) спектром принимаемого сигнала s(t,l,a). На рисунке также показан примерный вид спектра мощности помехи Wп(f). Важно, что помеха обычно имеет более широкий по сравнению с сигналом спектр и в пределах сигнальной полосы fÎ(f0-Fc, f0+Fc) имеет примерно постоянную спектральную интенсивность. Мешающее действие оказывают только те частотные составляющие помехи, которые попадают в полосу частот, занятую сигналом. Наоборот, частотные составляющие помехи, расположенные вне полосы частот fÎ(f0-Fc, f0+Fc), практически приему сигнала не препятствуют. Поэтому, выбирая математическую модель помехи, предназначенную для анализа процесса приема, не требуется точное воспроизведение этой моделью частотных составляющих, которые не попадают в полосу частот, занятую сигналом. Удобно аппроксимировать помеху n(t) гауссовским стационарным СП с равномерным ограниченным спектром (рис.3.12)

Спектральное представление стационарных случайных процессов - student2.ru (3.6.1)

Верхняя частота спектра F может быть выбрана произвольно, но такой, что F³f0+Fc. Ниже будет показано, что при этом конечные формулы и выводы от величины F не зависят. N0/2 - двухсторонняя спектральная плотность, т.е. мощность идеализированной помехи, приходящаяся на каждый 1Гц полосы fÎ(-F,F), равная спектральной плотности мощности реальной помехи на несущей частоте сигнала Wп(f0). СП такого вида рассматривался ранее (пример 3.5.1)

Модель (1) является примером идеализации, хорошо согласующейся с реальным явлением и вместе с тем существенно облегчающей математический анализ. Для дальнейшего исследования закономерностей приема сигналов на фоне помех n(t) нужно ввести меру вероятности различных реализаций СП n(t) на интервале наблюдения tÎ(0,T). Перейдем к определению этой меры.

Если спектр мощности Wп(f) стационарного в широком смысле СП n(t) ограничен полосой |f|£F, в частности, имеет вид (1), то реализации этого процесса представимы рядом Котельникова

Спектральное представление стационарных случайных процессов - student2.ru .(3.6.1)

Особенность состоит в том, что оговоренный выше вид спектра (1) относится не к спектру мощности Wп(f,w) реализации, представленной рядом (2), а к математическому ожиданию этого спектра Спектральное представление стационарных случайных процессов - student2.ru Wп(f,w) Спектральное представление стационарных случайных процессов - student2.ru =Wп(f) - к спектру мощности процесса. Вследствие этой особенности равенство (2) является не абсолютным, как для детерминированных процессов, а соблюдается в вероятностном (в среднеквадратическом) смысле: математическое ожидание квадрата ошибки - разности между левой и правой частями (2), равно нулю.

Для представления реализации СП n(t) на большом, но конечном интервале наблюдения tÎ(0,T) с достаточной для технических приложений точностью, как и для детерминированных процессов, можно пользоваться усеченным рядом (2), в котором сохраняются только те отсчеты, которые попадают в интервал наблюдения tÎ(0,T). Таким образом, реализация {n(t),tÎ(0,T)} однозначно определяется совокупностью своих отсчетов n0=n(0), n1=n(T0), ..., nm=n(mT0) т.е. вектором Спектральное представление стационарных случайных процессов - student2.ru , где T0=1/2F, m=2FT. Появление реализации {n(t), tÎ(0,T)} тождественно появлению совокупности отсчетов Спектральное представление стационарных случайных процессов - student2.ru . Поэтому ПВ реализации {n(t),tÎ(0,T)}, обозначаемая далее P0[n(t)], совпадает с ПВ p0(n) вектора n, т.е. с (m+1)-мерной ПВ процесса n(t)

Спектральное представление стационарных случайных процессов - student2.ru .(3.6.3)

Аналитическое выражение для ПВ (3) легко получить, если учесть, что гауссовская помеха n(t) согласно выбранной модели имеет корреляционную функцию (пример 3.5.1)

Спектральное представление стационарных случайных процессов - student2.ru .(3.6.4)

При этом все фигурирующие в (3) отсчеты помехи n0,n1,...,nm, сдвинутые друг относительно друга на интервалы, кратные T0=1/2F, статистически независимы. Поэтому ПВ (m+1)-го порядка равна произведению ПВ 1-го порядка p1(×)

Спектральное представление стационарных случайных процессов - student2.ru .

Так как ПВ центированной гауссовской случайной величины nk с дисперсией Спектральное представление стационарных случайных процессов - student2.ru , равна

Спектральное представление стационарных случайных процессов - student2.ru ,

то

Спектральное представление стационарных случайных процессов - student2.ru

Принимая во внимание, что энергия Эn реализации {n(t), tÎ(0,T)}

Спектральное представление стационарных случайных процессов - student2.ru ,

окончательно получаем

Спектральное представление стационарных случайных процессов - student2.ru , (3.6.5)

где C - независящий от реализации n(t) коэффициент

Спектральное представление стационарных случайных процессов - student2.ru .(3.6.6)

Выражение (5) представляет собой пример функциональной связи, ставящей в соответствие функции {n(t), tÎ(0,T)} число P0[n(t)]=Cexp(-Эn/N0). Такая функциональная связь называется функционалом. Поэтому в литературе ПВ P0[n(t)] называют функционалом ПВ помехи.

Если изменить выбранную в значительной мере произвольно верхнюю частоту F спектра помехи, то в функционале ПВ (5) изменится только величина коэффициента C. Однако в практических задачах представляет интерес не абсолютное значение ПВ P0[n(t)], а то, во сколько раз одна реализация помехи, положим n1(t), более или менее вероятна, чем другая - n2(t), т.е. отношение

Спектральное представление стационарных случайных процессов - student2.ru , (3.6.7)

которое от коэффициента C не зависит.

Часто в качестве модели помехи n(t) используется гауссовский стационарный процесс с равномерным неограниченным спектром (пример 3.5.2) - белый гауссовский шум (БГШ)

Спектральное представление стационарных случайных процессов - student2.ru .(3.6.8)

По соображениям, приведенным выше, мешающее действие этой помехи такое же, как и СП, выбранного нами в качестве предыдущей модели. Отличие состоит лишь в том, что для белого шума верхняя частота F®¥ . Выражение для функционала ПВ остается тем же (5). Отметим, при F®¥ согласно (6) коэффициент C®0, что в соответствии с пояснениями, приведшими к (7), не препятствует решению практических задач.

При использовании в качестве модели общей гауссовской помехи функционал ПВ помехи {n(t), tÎ(0,T)} определяется корреляционной функцией Rn(t1,t2) в общем случае нестационарного центрированного гауссовского СП. Запишем сначала для этой помехи многомерную ПВ любого конечного порядка m. Будем считать, что m сечений n(Dt)=n1,..,n(mDt)=nm равномерно распределены (Dt=T/m) на интервале tÎ(0,T). Согласно(3.3.5)

Спектральное представление стационарных случайных процессов - student2.ru , (3.6.9)

где C - коэффициент, не зависящий от реализации помехи n(t) и значений ее отсчетов nk и n, а [ Спектральное представление стационарных случайных процессов - student2.ru ]kℓ - элементы матрицы Спектральное представление стационарных случайных процессов - student2.ru , обратной к корреляционной R, с элементами Rkℓ=Rn(kDt,ℓDt) так что

Спектральное представление стационарных случайных процессов - student2.ru .(3.6.10)

Для получения функционала ПВ P0[n(t)] представим показатель экспоненты (9) в виде интегральной суммы и перейдем к пределу Dt®0 (m=T/Dt®¥). При этом kDt®tk , ℓDt®t , mDt®T и вводится обозначение

Спектральное представление стационарных случайных процессов - student2.ru .(3.6.11)

Тогда показатель экспоненты в выражении (9) можно преобразовать в интегральную сумму

Спектральное представление стационарных случайных процессов - student2.ru .

Функцию W(tk,t) называют обратной корреляционной функцией. Определяют ее из интегрального уравнения

Спектральное представление стационарных случайных процессов - student2.ru , (3.6.12)

получаемого преобразованием (10) в интегральную сумму

Спектральное представление стационарных случайных процессов - student2.ru

и перехода к пределу Dt®0 (m®¥ ), dkℓ/Dt®d(tk-t).

В результате функционал ПВ общей гауссовской помехи принимает вид (переменные интегрирования tk и t заменены на t1 и t2)

Спектральное представление стационарных случайных процессов - student2.ru .(3.6.13)

В частном случае, когда моделью помехи является белый шум Rn(t1,t2)=N0/2d(t1-t2), решением интегрального уравнения для обратной корреляционной функции

Спектральное представление стационарных случайных процессов - student2.ru

является

Спектральное представление стационарных случайных процессов - student2.ru (3.6.14)

и функционал ПВ помехи (13) преобразуется, как и следовало ожидать, в (5).

Формула (13) статистически исчерпывающе определяет модель общей гауссовской помехи, а также модели многих сигналов. Однако пользование этой моделью затруднено тем, что обратная корреляционная функция W(t1,t2) может быть определена аналитически в замкнутом виде из уравнения (12) для ограниченного вида функций корреляции Rn(t1,t2). В дальнейшем будем считать, что помеха n(t) имеет равномерный ограниченный или неограниченный спектр мощности и пользоваться моделью основной гауссовской помехи, т.е. функционалом ПВ (5). Только в некоторых специальных случаях будем пользоваться общей гауссовской моделью (13).

Функция правдоподобия

Запишем сначала функционал ПВ P{u(t), tÎ(0,T)}=P[u(t)], называемый для краткости ПВ процесса

Спектральное представление стационарных случайных процессов - student2.ru ,

представляющего собой сумму помехи n(t) с ПВ P0[n(t)], например, (3.6.5) и детерминированного процесса s(t). Процесс u(t) отличается от гауссовской помехи n(t) только тем, что имеет ненулевое математическое ожидание Спектральное представление стационарных случайных процессов - student2.ru u(t) Спектральное представление стационарных случайных процессов - student2.ru =s(t). Соответственно центрированное значение процесса u(t) совпадает с помехой Спектральное представление стационарных случайных процессов - student2.ru . Поэтому, если реализация процесса u(t) приняла некоторый конкретный вид Спектральное представление стационарных случайных процессов - student2.ru , то это значит, что реализация помехи приняла вид Спектральное представление стационарных случайных процессов - student2.ru . Следовательно, ПВ Спектральное представление стационарных случайных процессов - student2.ru того, что процесс u(t) примет конкретное значение Спектральное представление стационарных случайных процессов - student2.ru равна ПВ Спектральное представление стационарных случайных процессов - student2.ru того, что помеха n(t) примет значение Спектральное представление стационарных случайных процессов - student2.ru , т.е. Спектральное представление стационарных случайных процессов - student2.ru . Последнее равенство верно для любой реализации процесса u(t). Поэтому знак ^ можно опустить и окончательно записать

Спектральное представление стационарных случайных процессов - student2.ru . (3.7.2)

Таким образом, ПВ P[u(t)] гауссовского СП (1) с ненулевым математическим ожиданием Спектральное представление стационарных случайных процессов - student2.ru u(t) Спектральное представление стационарных случайных процессов - student2.ru =s(t) выражается через ПВ помехи P0[u(t)-s(t)].

Перейдем теперь к реальной модели наблюдения u(t). Будем считать, что наблюдение равно сумме

Спектральное представление стационарных случайных процессов - student2.ru ,(3.7.3)

основной гауссовской помехи n(t) и сигнала s(t,l1,..,l), являющегося квазидетерминированным СП, зависящим от совокупности случайных параметров l1,...,l. В данном случае не существенно, что собой представляют параметры l1,...,l , в частности, какие из них информационные, а какие - несущественные. Важно, что l1,...,l это вся совокупность случайных параметров сигнала (неизвестных на приемной стороне), и задание этой совокупности l1,...,l полностью определяет реализацию сигнала. При этом сигнал s(t,l1,...,l) становится детерминированной функцией времени s(t) и мы возвращаемся к модели (1).

Пример 3.7.1. Запишем простейший принимаемый радиолокационный сигнал

Спектральное представление стационарных случайных процессов - student2.ru .(3.7.4)

Параметры излучаемого сигнала tи и f0 заранее известны. Неизвестными на приемной стороне и воспринимаемые как случайные величины l1, l2, l3 являются параметры t, e, j . Задание конкретных значений этих параметров преобразует (3) в детерминированную функцию времени s(t).

Таким образом, при условии, что заданы конкретные значения всех случайных параметров сигнала l1,...,l функция s(t,l1,...,l) становится детерминированной функцией времени, равной математическому ожиданию наблюдения u(t). Поэтому условная ПВ P[u(t)/l1,...,l] реализации процесса u(t) при условии, что все параметры сигнала фиксированы и равны данным значениям l1,...,l представляет собой ПВ процесса u(t) вида (1) с математическим ожиданием Спектральное представление стационарных случайных процессов - student2.ru u(t) Спектральное представление стационарных случайных процессов - student2.ru =s(t, l1,...,l). В результате получаем

Спектральное представление стационарных случайных процессов - student2.ru . (3.7.5)

ПВ P[u(t)|l1,...,l], как это видно из (4), представляет собой функцию наблюдения u(t) и параметров l1,...,l. Если значение параметров l1,...,l зафиксировано (обусловлено), то P[u(t)|l1,...,l] - условная ПВ. Если, наоборот, зафиксировать наблюдение {u(t),tÎ(0,T)}, то P[u(t)|l1,...,l] является функцией параметров l1,...,l и называется функцией правдоподобия (параметров l1,...,l). Смысл названия в том, что функция P[u(t)|l1,...,l] характеризует относительную правдоподобность предположения, что наблюдение u(t) образовалось из сигнала с данной совокупностью параметров l1,...,l. Если P[u(t)|l1(1),...,l(1)]>P[u(t)|l1(2),...,l(2)], то более правдоподобно предположение, что наблюдение {u(t),tÎ(0,T)} образовалось из сигнала s(t,l1(1),...,l(1)), чем из сигнала s(t,l1(2),...,l(2)).

Совершенно аналогично, если гауссовская помеха n(t) имеет произвольную корреляционную функцию, отличную от 0,5N0d(t1-t2) или N0Fsinc2pF(t1-t2), то условная ПВ наблюдения u(t) (функция правдоподобия) P[u(t)|l1,...,l] определяется ПВ (3.6.13)

Спектральное представление стационарных случайных процессов - student2.ru (3.7.6)

Во многих практически важных случаях требуется функция правдоподобия не для всей совокупности параметров сигнала l1,...,l , а только для информационных. Положим, что в принимаемом сигнале s(t, l1,...,l) только l1 и l2 - информационные параметры, тогда как все остальные параметры l3,...,l- несущественные. Требуется определить функцию правдоподобия информационных параметров P[u(t)|l1,l2]. Особенность этой ситуации состоит в том, что имеется аналитическое выражение (4) или (5) для случая, когда зафиксирована вся совокупность параметров сигнала l1,...,l и сигнал превращается в детерминированную функцию s(t), являющуюся математическим ожиданием СП u(t). Однако необходимо определить условную ПВ наблюдения u(t) при условии, что зафиксированы параметры l1 и l2, тогда как остальные параметры l3,...,l могут иметь произвольные значения. Задача эта решается непосредственным применением правил теории вероятностей. Совместную ПВ несущественных параметров l3,...,l обозначим p(l3,...,l). Тогда условная (обусловлено значение l1 и l2) совместная ПВ наблюдения u(t) и параметров l3,...,l согласно формуле произведения вероятностей равна

Спектральное представление стационарных случайных процессов - student2.ru . (3.7.6)

Далее, пользуясь условием согласованности многомерной ПВ, т.е. усредняя (6) по параметрам l3,...,l находим интересующую нас ПВ

Спектральное представление стационарных случайных процессов - student2.ru

или

Спектральное представление стационарных случайных процессов - student2.ru . (3.7.7)

Следовательно, функция правдоподобия информационных параметров определяется по функции правдоподобия для всей совокупности случайных параметров сигнала путем статистического усреднения последней, т.е. усреднения с весом p(l3,...,l), по несущественным параметрам l3,...,l.

Приведем еще один пример определения функции правдоподобия для наблюдения

Спектральное представление стационарных случайных процессов - student2.ru ,

в котором информационным является дискретный параметр i=1,..m, а a - несущественный параметр (возможно векторный), принадлежащий непрерывному интервалу A. Иначе говоря, сигнал на входе может иметь вид одной из m заранее известных функций s1(t,a),...,sm(t,a) времени t и параметра a. Требуется записать функцию правдоподобия для сигнала si (для информационного параметра i). На основании изложенного

Спектральное представление стационарных случайных процессов - student2.ru

или для основной гауссовской помехи

Спектральное представление стационарных случайных процессов - student2.ru . (3.7.8)

В дальнейшем статистический синтез и анализ РТС основывается в значительной мере на использовании различных условных функционалов ПВ наблюдения - функций правдоподобия.

РАЗЛИЧЕНИЕ СИГНАЛОВ

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

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