Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Сигнали електрозв’язку- навчальний посібник з вивч

.pdf
Скачиваний:
163
Добавлен:
10.02.2016
Размер:
1.2 Mб
Скачать

41

Порівнювати значення KX(τ) і KX(0) для визначення статистичного зв'язку між значеннями процесу зсунутими за часом на τ не зовсім зручно, тому що ці значення залежать від дисперсії процесу (властивості 1 і 4). Визначення статистичного зв’язку зручно проводити за нормованою кореляційною функцією

RX (τ) = K X (τ) K X (0) .

(3.24)

Із співвідношення (3.22) випливає, що -1 £ RX(t) £ 1. Чим ближче значення RX(t) до 1, тим сильніше корельовані значення процесу зсунуті за часом на t. Якщо RX(t) = 0, то значення процесу зсунуті за часом на t, не корельовані, незалежні.

Для грубого опису кореляційних зв’язків уводять поняття інтервал (час)

кореляції процесу tк. Умовно вважать, що значення процесу, які зсунуті на

t £ tк істотно корельовані між собою, а значення процесу, які зсунуті на t > tк, не корельовані або корельовані несуттєво. Інтервал кореляції визначають по-

різному подібно тому, як оцінюється ширина спектра. Так, можна прийняти, що

 

tк =

 

RX (t)

 

dt .

(3.25)

 

 

0

 

 

 

 

 

Тут геометрично tк дорівнює основі прямокутника з висотою RX(0) = 1, що має ту ж площу, що й площа між кривою RX(t)½ при t > 0 і віссю абсцис.

Час кореляції tк можна визначити й інакше – як тривалість функціїRX(t)½ при t > 0 на рівні, наприклад, 0,1.

Аналогічно визначенню КФ процесу вводиться взаємна кореляційна фун- кція для характеристики зв’язку між значеннями двох випадкових процесів X(t) і Y(t), що зсунуті за часом на t,

 

 

 

K XY (τ) =

 

=

xyp2 (x, y, τ)dxdy ,

 

X (t)Y (t + τ)

(3.26)

−∞ −∞

де p2 (x, y, τ) – спільна густина імовірності значень стаціонарних процесів X(t) і

Y(t), що зсунуті за часом на τ.

Стаціонарний випадковий процес називається ергодичним, якщо будь- які статистичні характеристики його, знайдені шляхом усереднення по анса- мблю, збігаються з характеристиками, знайденими усередненням у часі однієї реалізації.

Середнє значення за часом позначають хвилястою лінією над залежніс- тю, що усереднюється за часом. Так, середнє в часі значення процесу X(t) ви-

значається усередненням за часом реалізації x(t)

X (t) = lim

1

T 2

 

 

x(t)dt.

(3.27)

 

 

T →∞ T T

2

 

Середня потужність процесу

 

 

 

 

42

 

 

1

T 2

 

PX

= lim

x2 (t)dt .

(3.28)

 

 

T →∞ T

T 2

 

Дисперсія процесу

 

1

T 2

 

D{X (t)} = lim

[x(t) - X (t)]2 dt.

(3.29)

 

T →∞ T

T 2

 

КФ ергодичного процесу визначається

 

 

1

T 2

 

K X

(t) = lim

x(t)x(t + t)dt .

(3.30)

 

 

T →∞ T

T 2

 

Два процеси X(t) і Y(t) називаються спільно ергодичними, якщо усереднення функцій від них по ансамблю приводить до того ж результату, що й усе-

реднення за часом. Взаємна кореляційна функція спільно ергодичних процесів

може бути визначена

 

 

1

T 2

 

K XY

(t) = lim

x(t) y(t + t)dt .

(3.31)

 

 

T →∞ T

T 2

 

Для рішення багатьох практичних задач досить розглядати лише середні

значення, дисперсії й КФ. Теорії, які обходяться лише цими характеристиками,

називаються кореляційними теоріями. У рамках кореляційної теорії природ-

но вважати стаціонарними такі випадкові процеси, у яких середнє значення й ди- сперсія процесу не залежать від часу, а КФ залежить лише від τ. Такі випад-

кові процеси називаються стаціонарними в широкому сенсі.

3.4 Спектральна густина потужності стаціонарного випадкового процесу

При вивченні детермінованих сигналів для спектрального аналізу використовується перетворення Фур'є. Можна спробувати знайти перетворення Фур'є від k-ої реалізації процесу xk(t) і знайти її спектральну густину як межу

 

 

T / 2

 

lim Sk

( jw) = lim

xk (t)ejωt dt .

(3.32)

T →∞

T →∞

 

 

 

 

T / 2

 

Однак спроба може виявитися безуспішною, оскільки при T ® ¥ немає гарантії, що k-а реалізація задовольняє умові абсолютної інтегруємості. Якби межа інтеграла (3.32) для k-ї реалізації й існувала, то вона була б спектральною характеристикою лише k-ї реалізації, а не процесу в цілому.

Можна показати, що

 

1

 

 

 

 

lim

| Sk ( jw) |2

=

K X (t)ejωτ dt ,

(3.33)

 

T →∞ T

−∞

 

де пряма лінія означає усереднення по ансамблі реалізацій. Співвідношення (3.33) є не що інше, як перетворення Фур'є від КФ процесу. Оскільки КФ характеризує процес у цілому, те й ліва частина в (3.33) є спектральною характеристикою всього процесу. Її позначають як GX(w).

Фізичний зміст функції GX(w) легко з'ясувати, якщо врахувати, що ½Sk(jw)½2 відповідно до рівності Парсеваля (2.45) характеризує розподіл енергії

43

процесу по частоті. У результаті ділення цієї функції на T одержимо розподіл потужності процесу по частоті.

Співвідношення (3.33) можна переписати у вигляді прямого й зворотного перетворень Фур'є

G (ω) =

K (τ)ejωτ dτ,

 

 

 

X

 

X

 

 

 

 

−∞

 

 

 

(3.34)

 

1

 

K X (τ) =

GX (ω)e jωτ dω.

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

На підставі (3.34) можна записати

 

 

1

 

K X

(0) =

GX (ω)dω = GX ( f )df .

(3.35)

 

 

 

−∞

−∞

 

Нам вже відома властивість КФ KX(0) = PX, враховуючи її, стверджуємо,

що функція GX(ω) дійсно характеризує розподіл потужності процесу по час-

тоті на інтервалі (−∞, ∞), а значення функції GX(ω) або GX(f) дорівнює потужності процесу в смугах в 1 Гц навколо частот +f і −f. У зв’язку зі сказаним функцію GX(ω) називають спектральною густиною потужності процесу (СГП). Та-

ким чином, спектральна густина потужності й кореляційна функція стаціо- нарного випадкового процесу зв'язані перетвореннями Фурє. Це твердження відоме як теорема Хінчина-Вінера. Розмірність СГП В2/Гц, вона збігається з розмірністю енергії й, напевно тому, іноді функцію GX(ω) називають енергетичним спектром процесу.

Оскільки функції KX(τ) і GX(ω) парні, то замість пари перетворень (3.34) можна записати пару косинус-перетворень Фур’є, які, як правило, більш прості в обчисленнях, ніж співвідношення (3.34)

 

 

 

 

 

GX (ω) = 2

K X (τ) cos ωτdτ,

 

 

 

0

 

 

(3.36)

 

1

 

 

K (τ) =

G (ω) cos ωτdω.

 

 

 

 

 

 

 

 

X

π 0

X

 

 

 

 

 

 

Знаючи функцію GX(ω), можна визначити ширину спектра процесу тими самими способами, що і ширину спектра детермінованого сигналу, наприклад, як протяжність області додатних частот, поза якою значення функції не перевищують значення 0,1max{GX(ω)}. Якщо спектр примикає до нуля, то ширину спектра визначають величиною Fmax (рис. 3.4, а), а якщо спектр смуговий, то ширину спектра визначають величиною F (рис. 3.4, б).

Оскільки функції GX(ω) і KX(τ) зв’язані перетворенням Фур’є, то відповідно до властивості перетворенням Фур’є зміни масштабу ω і τ, чим менший інтервал кореляції процесу, тим ширший спектр процесу й навпаки. Інакше кажу-

чи, інтервал кореляції й ширина спектра процесів обернено пропорційні величи- ни, а їх добуток є величиною порядку 0,5:

τкFmax = 0,5.

44

 

Gx ( f )

 

 

 

 

 

 

 

 

 

 

 

Gx ( f )

 

 

max Gx ( f )

 

 

 

 

 

max Gx ( f )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

0,1

 

 

 

 

 

0,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

Fmax

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б

 

 

 

 

 

а

 

 

 

 

Рисунок 3.4 – Визначення ширини спектра процесу:

а– спектр примикає до нульової частоти; б – смуговий спектр

3.5Гауссів випадковий процес

Найбільш часто в теорії й техніці зв'язку зустрічається так званий гауссів

(або нормальний) випадковий процес. Випадковий стаціонарний процес X(t) на- зивається гауссовим, якщо його одновимірна й двовимірна густини ймовірно-

сті визначаються наступними виразами

 

 

 

 

 

 

 

 

 

p(x) =

 

 

 

1

 

 

(x a)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

,

 

 

 

 

 

 

(3.37)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2πσ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, τ) =

 

 

1

 

 

 

(x a)2 − 2R

X

(τ)(x a)(x

2

a) + (x

2

a)2

 

p

 

(x , x

 

 

 

 

 

 

exp

1

 

 

 

 

 

1

 

 

 

 

, (3.38)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

2

 

 

2πσ

2

2

 

 

 

 

 

 

 

 

 

 

 

2

 

2

(τ)]

 

 

 

 

 

 

 

 

 

 

 

1 − RX (τ)

 

 

 

 

 

 

 

 

 

 

 

[1 − RX

 

 

 

 

де σ2

дисперсія процесу;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

середнє значення процесу;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

RХ(τ) –

значення нормованої кореляційної функції процесу.

 

 

 

 

 

 

Щоб визначити двовимірну густину ймовірності нормального випадково-

го стаціонарного процесу, досить знати лише його КФ. Таким чином, нормальні стаціонарні процеси можуть відрізнятися один від іншого видами КФ й, відповідно, СГП.

Одновимірна функція розподілу ймовірностей нормального процесу на-

ступна

 

 

 

 

 

x a

 

 

F (x) = 1 − Q

 

 

,

(3.39)

 

σ

 

 

 

 

 

 

 

 

 

 

 

1

 

 

t

2

 

 

де

Q(z) =

 

 

exp(−

 

)dt

(3.40)

 

 

 

 

 

 

 

 

z

2

 

 

– гауссова Q-функція або доповнення до функції розподілу ймовірностей. Графіки функцій (3.39) і (3.40) наведені на рис. 3.5.

 

 

 

 

 

 

 

 

 

 

 

45

1,0

 

 

 

 

 

 

 

 

 

 

 

F(x)

 

 

 

 

 

0,4

 

 

 

σ = 1

 

 

a = 0

a = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a = 1

 

0,5

 

 

 

 

 

 

a = 0

 

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

σ = 2

 

σ = 1

 

 

σ = 2

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

–4

–2

0

2

x

4

–4

–2

0

2

x

4

 

 

а

 

 

 

 

 

 

б

 

 

Рисунок 3.5 – Гауссів розподіл:

а – функція розподілу ймовірностей; б – густина ймовірності

Гауссів смуговий процес зручно представити через квадратурні складові

X (t) = A(t) cos(ω0 t + Φ(t)) = Ac (t) cos ω0 t + As (t) sin ω0 t = X c (t) + X s (t) , (3.41)

де A(t) і Φ(t) – обвідна і фаза процесу;

Xc(t) і Xs(t) – квадратурні складові процесу;

ω0 – деяка частота, що належить смузі частот процесу X(t).

Квадратурні складові – некорельовані процеси, що мають гауссів розподіл імовірностей, їхні дисперсії однакові й дорівнюють половині дисперсії процесу X(t).

Обвідна A(t) і фаза Φ(t) також є некорельованими процесами. Обвідна A(t)

має релеєвський розподіл імовірностей (рис. 3.6)

 

a

 

 

 

a2

 

 

 

 

 

 

 

 

exp −

 

 

 

 

,

a > 0,

 

σ

2

2

 

 

p(a) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a ≤ 0;

 

0,

 

 

 

 

 

 

 

(3.42)

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

a2

a > 0,

 

1− exp

 

 

 

 

,

 

 

2

 

F (a) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

 

 

 

 

 

 

 

a ≤ 0.

 

У виразах (3.42) σ2

дисперсія процесу X(t). Числові характеристики ре-

 

 

 

 

 

=

 

 

 

леєвського

процесу:

середнє

значення

A(t)

π / 2σ ,

дисперсія

D[ A(t)] = (2 − π / 2)σ2 , середня потужність PA = 2σ2.

Фаза Φ(t) має рівномірний розподіл імовірностей на інтервалі (0, 2π) (рис. 3.7)

 

 

 

 

 

46

 

1

,

0 £ j < 2p,

 

 

p(j) =

2p

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

 

j < 0, j ³ 2p;

(3.43)

0,

 

j < 0,

 

 

 

 

 

 

 

0 £ j < 2p,

 

p(j) = j/2p,

 

 

 

 

 

j ³ 2p.

 

 

1,

 

 

 

1

 

 

 

 

F(a)

σ = 1

σ = 2

 

 

0,5

 

 

 

 

0

 

 

 

 

0

2

4

a

6

 

 

а

 

 

0,6

σ = 1

 

 

 

p(a)

 

 

 

 

 

 

 

0,4

 

 

 

 

0,2

 

σ = 2

 

 

 

 

 

 

0

 

 

 

 

0

2

4

a

6

 

 

б

 

 

 

 

Рисунок 3.6 – Розподіл Релея:

 

 

а

функція розподілу ймовірностей; б

густина ймовірності

 

1

 

 

 

 

 

 

F(φ)

 

 

p(φ)

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

0,5

 

 

 

 

 

0

 

 

0

0

φ

0

φ

 

 

а

 

 

 

б

 

 

 

Рисунок 3.7 –

Рівномірний розподіл:

 

 

а– функція розподілу ймовірностей; б – густина ймовірності

3.6Білий шум

Білим шумом називається випадковий процес, СГП якого є постійною

величиною

G(w) =

N 0

,

- ¥ < w < ¥ або

G(ω) = N0 , 0 < ω < ∞ ,

(3.44)

 

2

 

 

 

 

де N0 – потужність процесу в смузі 1 Гц.

 

 

Визначенням

(3.44) відповідають

графічні залежності, наведені на

рис. 3.8.

КФ білого шуму визначається як зворотне перетворення Фур'є від (3.44)

K (τ) =

1

N0

 

jωτ

dω =

N0

δ(τ).

 

−∞

 

e

 

 

2

 

2

На рис. 3.9 наведено графік КФ білого шуму.

G(ω)

 

G(f)

 

N0/2

 

N0

 

 

 

 

 

 

 

 

 

0

ω

0

f

а

 

 

 

б

Рисунок 3.8 – Спектральна густина потужності білого шуму: а – двосторонній спектр; б – односторонній спектр

47

(3.45)

K(τ)

N0 δ(τ)

2

0 τ

Рисунок 3.9 – Кореляційна функція білого шуму

За виразом (3.45) та рис. 3.8 бачимо, що інтервал кореляції прямує до ну-

ля, тобто два значення білого шуму для яких τ ≠ 0 є незалежними. Очевидно, що білий шум має нескінченну потужність і фізично існувати не може, однак є зручною математичною моделлю, яка часто використовується в теорії і техніці зв'язку.

3.7 Перетворення випадкових процесів лінійними електричними ко-

лами

Нагадаємо, що за визначенням лінійним колом називається таке, для кого виконується принцип суперпозиції реакція кола на суму дій дорівнює сумі

реакцій на кожну дію окремо.

При дослідженні проходження випадкових процесів через лінійні кола вважається, що відомі статистичні характеристики вхідного випадкового про-

 

 

 

 

 

цесу X(t) та передатна функція

лінійного

кола

X(t)

 

 

 

 

Y(t)

 

 

 

 

Лінійне коло

 

 

H(jω). Необхідно знайти характеристики вихідно-

 

 

 

 

Рисунок 3.10

До вивчення

го процесу Y(t) (рис. 3.10).

лінійного

кола

проходження випадкових

СГП процесу на виході

звязана з СГП вхідного процесу через квадрат

процесів через лінійні кола

 

 

 

 

 

 

АЧХ кола

 

 

 

 

 

 

 

 

GY(ω) = GX(ω)H 2(ω).

 

(3.46)

Зокрема, якщо вхідний процес білий шум, то СГП вихідного процесу повторює

квадрат АЧХ лінійного кола.

КФ процесу на виході лінійного кола визначається як перетворення Фур'є від СГП процесу

KY(τ) =

1

G (ω) cos ωτ dω .

(3.47)

 

 

π 0

 

 

 

Y

 

Нехай для наочності X(t) – білий шум з однобічною СГП GX(f) = N0, 0 ≤ f < ∞, подається на вхід ідеального ФНЧ з АЧХ

 

 

 

 

 

 

 

 

 

48

 

 

H 0 , 0 ≤ f < Fзр ,

 

 

 

(3.48)

 

 

H(f) =

 

³ Fзр ,

 

 

 

 

 

0,

f

 

 

 

 

де Fзр – частота зрізу ФНЧ. Тоді СГП процесу Y(t):

 

 

 

 

 

 

 

N

0 H 02 , 0 £ f < Fзр ,

 

 

(3.49)

 

GY(f) = GX(f)×H2(f) =

0,

f ³ Fзр .

 

 

 

 

 

 

 

 

 

 

 

СГП процесу Y(t) показано на рис 3.11, а. Шум з СГП такого виду нази-

вають квазібілим, оскільки він має частотні складові однакової амплітуди,

але в обмеженому діапазоні частот.

 

 

 

 

 

 

 

КФ процесу Y(t):

 

 

 

 

 

 

 

 

 

 

KY(t) = N0 H 02

Fзр sin 2pFзрt .

 

 

 

(3.50)

 

 

 

 

2pFзрt

 

 

 

 

На рис 3,11, б показана нормована КФ RY(t) = KY(t)/KY(0). Інтервал коре-

ляції процесу Y(t)

tк = 1/(2Fзр).

 

 

 

 

 

 

 

GY(f)

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N0H02

 

RY(τ)

 

 

 

 

 

 

 

Fзр

f

0

 

 

τк

 

 

 

 

 

 

 

 

 

 

 

 

–0,3

 

 

 

 

 

 

 

 

− 1

 

 

 

 

 

 

 

− 2

 

0

1

t

2

a

 

 

Fзр

 

Fзр

 

Fзр

Fзр

 

 

 

б

 

 

 

 

 

 

 

 

 

 

Рисунок 3.11 –

Характеристики процесу Y(t) при фільтрації ФНЧ: а

СГП; б – КФ

Середня потужність процесу Y(t):

Fзр

 

PY = GY ( f ) df

= N0 H 02 df = N0 H 02 Fзр

(3.51)

0

0

 

В інженерних розрахунках часто виникає задача визначення потужності шуму на виході довільного лінійного кола (неідеального) за умови дії білого шуму з однобічною СГП N0 на його вході1. Для легкого вирішення цієї задачі одноразово визначається так звана шумова смуга кола, як інтеграл від квад-

рата нормованої АЧХ кола

H 2 ( f )

 

 

Fш =

 

 

df ,

(3.52)

H

2

0

max

 

 

 

 

 

де Hmax – максимальне значення АЧХ.

1 Насправді шум є квазібілим, але при цьому займає смугу частот значно ширше смуги пропускання лінійного кола і тому може трактуватися, як білий

Рисунок 3.11 – До вивчення проходження випадкових процесів через нелінійні кола
Нелінійне коло у = f(x)

 

49

За відомої шумової смуги PY обчислюється виразом

 

PY = N0×Fш× H max2 .

(3.53)

У ідеального ФНЧ шумова смуга фільтра Fш = Fзр.

Щодо розподілу ймовірностей процесу, що пройшов через лінійне коло,

відомо наступне. Якщо на вході лінійного кола процес гауссів, то вихідний про- цес також гауссів вид розподілу не змінюється, змінюються тільки його па- раметри. Якщо на вході кола процес не гауссів, то вид розподілу змінюється, і вихідний процес має розподіл імовірностей більш близький до гауссового, ніж

розподіл вхідного процесу.

Фільтрація є вузькосмуговою, якщо смуга пропускання кола значно ме-

нша ширини спектра вхідного процесу. При вузькосмуговій фільтрації має місце явище нормалізації процесу, що полягає в наступному – незалежно від виду

розподілу вхідного процесу, розподіл імовірностей процесу на виході кола є га-

уссовим. Доказ цього явища базується на центральній граничній теоремі теорії ймовірностей.

3.8 Перетворення випадкових процесів нелінійними електричними колами

Тут розглядаються лише безінерційні електричні кола, тобто такі, ре-

акція яких на дію миттєво припиняється після її закінчення. Коло називається нелінійним, якщо для нього не виконується принцип суперпозиції. Варто відзначити, що присутність одного нелінійного елементу робить усе електричне коло нелінійним.

При дослідженні проходження X(t) Y(t) випадкових процесів через нелінійні безінерційні кола вважається, що відомі статистичні характеристики вхідного процесу X(t) і залежність y = f(x) між миттєвими значеннями вхідного й вихідного процесів. Необхідно знайти

характеристики вихідного процесу Y(t) (рис. 3.11).

Найпоширенішою функцією f(x) для опису нелінійних перетворень є полі- ном степеня n

f(x) = a0 + a1x + a2 x2 + ... + an xn,

(3.54)

де a0, a1, a2,..., an – коефіцієнти полінома.

Коефіцієнти й степінь поліному визначаються в результаті апроксимації характеристики реального електричного кола або виходячи з деяких припущень. Крім поліноміальної залежності (3.54) використовуються й інші залежності.

Кожна із складових (3.54) вносить свій внесок у формування значень реакції нелінійного кола на вхідну дію. Так, a0 описує появу постійної складової при х = 0; a1x – лінійний доданок, що забезпечує пропорційне відображення

50

значень х в y; a2x2 – квадратичний доданок, a3x3 – кубічний доданок і т.д. забезпечують внески, пропорційні х2, х3 і т.д.

Найпростіша дія – гармонічне коливання x(t) = A1cos2pf1t. У цьому випадку сигнал на виході нелінійного кола буде наступним:

y(t) = a0 + a1A1cos2pf1t + a2A12cos22pf1t + ... + an A1ncosn2pf1t.

(3.55)

Якщо скористатися формулами кратних аргументів, то одержимо

 

y(t) = Y0 + Y1cos 2pf1t + Y2cos2p2f1t + ... + Yncos2pnf1t,

(3.56)

де Y0 – постійна складова реакції;

Y1, Y2, ..., Yn – амплітуди першої, другої, ..., n-ї гармонік реакції.

Таким чином, реакція на гармонічну дію містить постійну складову й га- рмоніки частоти дії це принципово відрізняє нелінійні кола від лінійних, у яких

нові складові не виникають.

У випадку бігармонічної дії

x(t) = A1cos2pf1t + A2cos2pf2t,

(3.57)

підхід до визначення реакції такий самий, як і використаний вище – вираз для x(t) підставляється в поліном (3.54). При зведенні суми (3.57) у квадрат, куб і т.д. з’являться степені косинусоїд частот f1 і f2, що після перетворень дає вираз виду (3.56) для коливань частот f1 і f2. Але з’являються ще й добутки косинусоїд і їхніх степенів. Добуток косинусоїд приводить до появи складових сумарних і різницевих частот.

У загальному випадку будуть мати місце складові комбінаційних час-

тот

fкомб =½pf1 ± qf2½,

(3.58)

де p, q – цілі числа 0, 1, 2, ..., але такі, що p + q £ n. Їхня сума N = p + q назива-

ється порядком комбінаційної частоти.

Так, якщо n = 3, то в спектрі реакції можуть бути складові частот f1, f2, 2f1, 2f2, ½f1±f2½, 3f1, 3f2, ½2f1 ± f2½, ½f1 ± 2f2½ і постійна складова. Амплітуди складо-

вих залежать від амплітуд А1 і A2 і коефіцієнтів полінома (3.54). Якщо амплі-

туди і фази складових дії на коло є випадковими, то, відповідно, випадковими будуть амплітуди і фази складових на комбінаційних частотах реакції.

Визначити СГП вихідного процесу GY(f) можна в такий спосіб: визначити спочатку КФ вихідного процесу KY(t), а потім виконати над нею перетворення Фур’є. Виходячи з визначення КФ

∞ ∞

 

KY (τ) = ∫ ∫ f (x1 ) f (x2 ) p2 (x1 , x2 , τ)dx1dx2 ,

(3.60)

−∞−∞

 

де f(x) – функція, що описує нелінійне коло;

р2(х1, х2, t) – двомірна густина імовірності вхідного процесу.

При проходженні випадкового процесу через нелінійне коло вид розподілу миттєвих значень істотно змінюється.