Материал: методичка по лабам Рцис

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Определения Δωэф поясняет рис. 6.3, а: площадь штрихового прямо-

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

На практике Δωэф обычно измеряют на уровне 0,5 от максимального

значения СПМ, что удобнее. При этом Δωэф0,5

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

ше, чем Δωэф , но близкой к ней.

 

 

 

 

 

 

 

Ширину КФ (интервал корреляции τк ) определяют как

 

1

 

τк =

 

R(τ)

 

d τ ,

(6.10)

 

 

 

 

 

R(0)

 

0

 

 

 

 

 

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

(рис. 6.4, б).

WW((f))

WWmaxax

ff00,55

Wmaxmax/22

 

0

 

 

 

 

 

 

 

 

 

ffэфэф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f0,5

 

W(

(f)

 

 

 

0,5

 

 

 

 

 

 

 

 

W

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wmaxax

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wmaxax/22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

эфф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

rr(ττ)

11

кк

0,50,5

 

0

 

 

f

0

0

τ

 

 

 

τ

б

Рис. 6.3

rr(τ)

 

 

кк

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,50,5

00

f 0 ττ

б

Рис. 6.4

66

В силу свойств преобразований Фурье (6.7), (6.8) время корреляции и

ширина

СПМ находятся в обратной зависимости: τк 2π / Δωэф, или

τк 1 /

fэф , т. е. широкополосному случайному процессу соответствует ма-

лый интервал корреляции и наоборот.

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

W (ω) = W (ω)

 

Kɺ

(ω)

 

2

,

(6.11)

 

 

вых

 

 

 

 

 

 

 

где Wвых(ω), W (ω) – СПМ на выходе и входе цепи соответственно; Kɺ (ω)

комплексный коэффициент передачи. Справедливость (6.11) вытекает непосредственно из определения функции СПМ. Далее, используя соотношения

(6.7) или (6.8), находят Rвых(τ) – КФ выходного процесса.

При временном подходе по известным КФ входного процесса и импульсной характеристике цепи h(t) определяют КФ процесса на выходе

 

Rвых(τ) = R(λ)Bh (τ − λ)dλ

(6.12)

−∞

как свертку R(τ) и Bh (τ) – КФ импульсной характеристики цепи h(t), причем

Bh (τ) = h(t)h(t − τ)dt .

−∞

Справедливость выражения (6.12) следует из справедливости принципа суперпозиции для линейной цепи: выходные значения случайного процесса

ξвых(t) можно найти методом интеграла наложения, зная входные значения

ξ(t) и импульсную характеристику h(t). Затем, подставив ξвых(t) в общую формулу (6.2), получим (6.12). Далее, используя соотношения (6.7) или (6.8), находим Wвых(ω) .

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

67

входной шум без большой погрешности заменяют математической моделью – белым шумом. Белый шум – это случайный процесс, спектральная плотность мощности которого W(ω) постоянна на всех частотах: W(ω) = W0,

ω (−∞,∞) . Значение W0 выбирают равным значению СПМ реального вход-

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

Легко показать, что для белого шума R(τ) = W0δ(τ) , где δ(τ) –

дельта-

функция Дирака, откуда следуют соотношения:

 

Wвых(ω) = W0

 

Kɺ(ω)

 

2 , Rвых(τ) = W0Bh (τ) ,

(6.13)

 

 

имеющие непосредственное отношение к настоящей лабораторной работе.

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

графа и анализатора спектра обладает определенной «памятью», послесвечением: зафиксированная за период развертки прибора осциллограмма отрезка процесса не сразу исчезает; на нее накладывается осциллограмма следующего отрезка и т. д., так что наблюдается, в сущности, сумма осциллограмм, соответствующих различным отрезкам наблюдаемого сигнала (для анализатора спектра – СПМ этих отрезков); иногда это свойство, аналогичное свойству инерционности подвижной части электромеханического прибора, используют для приближенного усреднения (интегрирования) поступающего на осциллограф сигнала. Например, в режиме автозапуска развертки и при подаче на осциллограф непрерывного случайного сигнала распределение яркости «шумовой дорожки» по Y-оси экрана будет примерно соответствовать функции плотности распределения случайного процесса (известное правило «трех сигм» для оценки дисперсии нормального случайного процесса). Аналогично объясняется наблюдение на экране спектроанализатора именно усредненной СПМ.

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

лографа исследуемый сигнал, установить ждущий режим развертки осциллографа и отрегулировать уровень запуска развертки так, чтобы запуск происходил только от достаточно больших выбросов процесса (рис. 6.5, а), то изображение на экране ξ′(t) будет с точностью до постоянного множителя повторять КФ процесса, ξ′(t) ~ R(t). При этом начальный участок осцилло-

граммы совпадает с формой КФ практически точно, а по мере увеличения t расхождение увеличивается.

68

ξ(t)

УровеньУровень запускаапуска

 

ξ

t

 

 

 

 

U

ξ′ξ′(tt)

 

 

 

 

00

 

t

 

 

 

t

 

 

а

б

 

 

 

Рис. 6.5

На рис. 6.5, а представлена реализация случайного процесса ξ(t). На рис. 6.5, б показана наблюдаемая осциллограмма ξ′(t) , сформированная

наложением нескольких осциллограмм, соответствующих выделенным на рис. 6.5, а отрезкам случайного процесса ξ(t) и являющаяся их средним или математическим ожиданием, поскольку для эргодических процессов КФ (6.2), полученная усреднением по времени, совпадает с корреляционной функцией, полученной усреднением по реализациям с помощью оператора математического ожидания R(τ) = M [ξ(t)ξ(t − τ)] .

При запуске развертки осциллографа в момент времени t0 уровнем U

значение ξ(t0 ) = U становится детерминированным и дальнейшее поведение процесса определяется условной плотностью вероятности

 

 

pξ (ξ,t0 + τ

 

ξ=U , t0 ) =

p2 (ξ =U ,t0;ξ,t0 + τ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

pξ (ξ =U ,t0 )

 

 

 

или, сразу положив t0 = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pξ (ξ, τ

 

ξ(0)=U ) =

p2 (ξ(0) =U ;ξ, τ)

=

p2 (U ,ξ, τ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(6.14)

 

 

 

 

pξ (ξ(0) =U )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pU (U )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полагая процесс ξ(t) гауссовским, запишем при ξ(0) = U его двухмерную

плотность вероятности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ξ,U , τ) =

 

 

 

1

 

 

 

 

 

 

 

ξ2 − 2Ur

(τ)ξ + U

2

 

p2

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

2[1 − r2 (τ)]

 

 

 

2

 

 

 

2

(τ)

 

 

 

 

 

2πσ

1 − r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где r(τ) –

коэффициент корреляции,

определенный выражением (6.4);

σ2

дисперсия процесса, а плотность вероятности в знаменателе выражения

(6.14)

69

 

(U ) =

 

1

 

 

U

2

 

p

 

 

exp −

 

.

 

 

 

 

 

U

 

 

2πσ

 

2

 

 

 

 

 

 

 

 

 

 

Тогда условная плотность вероятности (6.14) определится как

 

 

 

 

 

 

1

 

 

 

 

2

 

 

 

p (ξ, τ

ξ(0) =U

) =

 

 

 

 

 

 

exp −

[ξ − Ur(τ)]

 

.

(6.15)

 

 

 

 

 

 

 

 

 

ξ

 

 

2πσ

2

[1 − r

2

(τ)]

 

2[1 − r2 (τ)]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формула (6.15) описывает гауссовский случайный процесс с математи-

ческим ожиданием Ur(τ)

и дисперсией σ2[1 − r(τ)], где τ отсчитывается от

максимума.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, наблюдаемая на экране осциллограмма

 

с точно-

ξ (t)

стью до постоянного множителя будет повторять КФ процесса ξ(t), ξ′(t) ~

Ur (t ) = UR (t )R (0) . Начальный участок осциллограммы совпадает с формой

КФ практически точно, так как дисперсия σ2[1 − r(0)] = 0 ; далее, осцилло-

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

lim σ2[1 − r(t)] = σ2 . Практически этот предел достигается при значении t, со-

t→∞

ставляющем от одного до двух-трех интервалов корреляции τк процесса ξ(t).

6.2. Описание лабораторной установки

Установка (рис. 6.6) состоит из лабораторного макета и подключенных к нему осциллографа и анализатора спектра.

Лабораторный макет включает в себя генератор шума, набор фильтров

(RC, RLC и ФНЧ с почти прямоугольной АЧХ) и коррелятор, позволяющий измерять нормированные КФ сигналов. Коррелятор реализует операцию,

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

го время τ задержки линии можно изменять от 0 до 9,9 мкс с шагом 0,1 мкс.

Переключатель S1 макета позволяет выбрать вид исследуемой цепи; S2 –

выбрать параметры выбранной цепи (кроме ФНЧ); S3 – подключать осцилло-

граф и анализатор спектра ко входу (положение 1) или выходу (положение 2)

исследуемой цепи.

70