Линейные преобразования случайных сигналов

При анализе прохождения стационарного СП через линейные электрические цепи (рис. 1) будем полагать, что режим цепи установившийся, т.е. после подачи на вход цепи сигнала все переходные процессы, связанные с включением, закончились. Тогда и выходной СП будет стационарным. Рассматриваемая задача будет состоять в том, чтобы по заданной корреляционной функции входного сигнала или его спектральной плотности мощности определить B (t) или G (w) выходного сигнала.

 

Рис. 1

 

Сначала рассмотрим решение этой задачи в частотной области. Входной СП задан своей спектральной плотностью мощности G_х (). Выходная спектральная плотность мощности G _y(w) определяется по формуле

G_y () = G_х () K ²(), (1)

 

где K ²() - квадрат модуля комплексной передаточной функции цепи. Возведение в квадрат модуля основано на том, что искомая характеристика является действительной функцией частоты и энергетической характеристикой выходного процесса.

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

B_x () = F ^-1 [ G_x ()]; F ^-1 [ K ²()] = B_h ()

- корреляционная функция импульсной характеристики исследуемой цепи:

B_h ()= h (t) h (t - ) dt.

 

Таким образом, корреляционная функция выходного СП есть

B_y () = B_x () B_h () = Bx(t) B_h (t - t) dt.

 

ПРИМЕР 1 прохождения стационарного случайного широкополосного сигнала через RC -цепь (фильтр нижних частот), представленную схемой на рис. 2.

Широкополосность понимается так, что энергетическая ширина спектра входного СП намного больше полосы пропускания цепи (рис. 3). При таком соотношении между формой K ²() и G_x () можно не рассматривать ход характеристики G_x () в области верхних частот.

 

Рис. 2

 

Учитывая, что в полосе частот, где K ²(w) существенно отличается от нуля, спектральная плотность мощности входного сигнала равномерна, можно без существенной погрешности входной сигнал аппроксимировать белым шумом, т.е. положить G_x () = G ₀ = const. Такое предположение существенно упрощает анализ. Тогда G_y () = G ₀ K ²()

Для заданной цепи

 

K ²() = 1/[1 + ( RC)²], тогда G_y () = G ₀/[1 + ( RC)²].

 

Рис. 3

 

Определим энергетическую ширину спектра выходного сигнала. Мощность выходного СП

 

P_y = s_y ² = (2p) ^-1 G_y () d  = G ₀/(2 RC), тогда

D _э = (G0)-1 Gy () d  = p/(2RC).

 

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

Спектральная плотность мощности по форме повторяет квадрат модуля комплексной передаточной функции цепи. Максимальное значение G_y () равно G ₀. Максимальное значение корреляционной функции выходного СП (его дисперсия) равна G ₀/(2 RC). Нетрудно определить площадь, ограниченную корреляционной функцией. Она равна значению спектральной плотности мощности при нулевой частоте, т.е. G ₀:

 

.

 

Рис. 4

 

Энергетической (шумовой) полосой пропускания электрической цепи называется полоса частот, численно совпадающая с энергетической шириной спектральной плотности мощности сигнала на выходе цепи при воздействии на вход цепи белого шума. В заданной случае D _э = p/(2 RC). Сравним ее с полосой пропускания _гр этой же цепи на уровне 0,707. Так как _гр = 1/(RC), то Dw_э = p/2 _гр, то есть D _э в p/2 раз больше _гр.

Определим корреляционную функцию сигнала на выходе RC -цепи при воздействии на ее вход белого шума.

Так как выходная спектральная плотность мощности уже определена, то можно вычислить искомую функцию обратным преобразованием Фурье. Но в рассматриваемом случае проще анализ выполнить во временной области, то есть B_y () = B_x () B_h (), но так как B_x () = W₀ d(), то B_y () = W ₀ B_h () (учитывая фильтрующее свойство дельта-функции).

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

h (t) = 1/(RC) exp[- t /(RC)], t ³ 0, то

B_h () = h(t)h (t - ) dt = 1/(2 RC)exp[-| |/(RC)], -¥ <  < ¥.

 

На рис. 5 представлены корреляционные функции (рис. 5а) и спектральные плотности мощности (рис. 5б) для двух значений постоянной времени заданной цепи (RC)₁ < (RC)₂. Дисперсия выходного СП _y ² = B_y (0) = = G ₀/(2 RC).

Площадь под кривой B_y () равна значению спектральной плотности мощности при  = 0, есть G ₀. Из сравнения графиков на рис. 5 следует, что с уменьшением полосы пропускания цепи начальное (максимальное) значение корреляционной функции B_y (0) уменьшается, что связано с уменьшением мощности выходного сигнала, и корреляционная функция изменяется медленнее с увеличением RC заданной цепи.

Нетрудно рассчитать интервал корреляции выходного СП

 

 = .

 

Откуда следует, что интервал корреляции выходного СП равен постоянной времени цепи.

 

Рис. 5

ПРИМЕР 2 анализа прохождения белого шума через колебательный контур (рис. 6). Чтобы придать этой задаче физический смысл, сводим задачу, как и предыдущую, к аппроксимации входного сигнала белым шумом.

 

Рис. 6

 

Чтобы использовать такой прием, входной сигнал должен иметь спектральную плотность мощности, неизменную в пределах практически значимых значений ординат АЧХ цепи. Тогда G_x () можно считать равной G ₀, а входной СП – белым шумом (рис. 7).

Передаточная функция такой цепи K (j )= K _р/[1+ j 2( - _р)/( _р Q _э)]; 0 <  < ¥, где K _р - коэффициент передачи цепи при резонансной частоте _р, то есть K _р = R _эр/(R _эр + R); Q _э = Q /(1 + R _эр/ R) - добротность шунтированного нагрузкой R колебательного контура, его постоянная времени _к = 2 Q _э/ _р = 2/(D )_0,7, то есть обратная половине полосы пропускания контура на уровне 0,707. Квадрат модуля передаточной функции K ²() = K _р²/[1 + (  - _р)² _к²]. Найдем дисперсию процесса на выходе цепи

 

_y ² = G ₀ K _р²/(p _к) (p/2 + arctg2 Q _э)» G ₀ K _р²/(2 _к) при Q _э >> 1.

 

Рис. 7

 

Оценим энергетическую полосу пропускания колебательного контура (рис. 6) D _э = (G ₀ K _р²)^-1 G_y () d  » p/ _к. Сравним с полосой пропускания по уровню 0,707 (-3 дБ). Так как _к = 2 Q _э/ _р, то D _э = p/2 (D )_0,7.

Вычислим корреляционную функцию выходного процесса (рис. 8):

B_y () = G ₀ K _р²/(2 _к) exp(-| |/ _к) cos _р ; -¥ <  < ¥.

 

Рис. 8

 

Если рассматривать анализ контуров с разными добротностями, то можно увидеть различия в реализациях выходных процессов: рис. 9 при добротности Q ₁ и рис. 10 при добротности Q ₂.

Рис. 9

 

 

Сравнительный анализ показывает, что увеличение добротности приводит к снижению полосы пропускания контура, а значит, к снижению средней скорости изменения огибающей во времени (можно сравнить с влиянием на огибающую АМК снижения частоты модулирующего колебания). белый шум сигнал линейный преобразование