книги из ГПНТБ / Коновалов Г.В. Импульсные случайные процессы в электросвязи
.pdf
|
Г/2 |
|
f |
j и Л О е - " - ' dt = |
|
|
-Г/2 |
|
где |
Г/2 |
|
&,1 о |
|
|
j, |* |
{t) dt, |
|
|
Г/2 |
|
|
Г/2 |
ип (t) cos к (£>0tdt, |
|
I* |
|
|
—Г/2 |
|
_21_ (вда — І^к) |
При |
К > О, |
1 К * + і Ьпк) |
|
(4.6) |
при |
к < О, |
Г/2
bnK = ~Y I* ы„ (t) sin ,<cco0 1 dt.
—Г/2
В случае неограниченной длительности импульсов, когда k-я реализация случайной последовательности импульсов имеет вид
СО
l W (t)= У tnk){ t - t {nk)),
п=я— 0 0
а относительно функций %{k)(t) известно только, что они ограниче
ны и стремятся к нулю при |-/|-ч-оо, целесообразнее .'представлять в виде ряда не отдельные импульсы &-й реализации процесса, а достаточно длинные отрезки каждой реализации, на интервале разложения (—Г/2, Г/2) совпадающие с соответствующей реализа
цией рассматриваемого импульсного случайного процесса |
(см. па |
раграф 2 .1 ). |
|
Тогда, считая, что импульсный случайный процесс эргодический |
|
и характеризуется условием |
|
||ß * (t)|d t< оо, |
(4.7) |
-00 |
|
усеченную k-ю реализацию этого процесса можно представить ря дом Фурье на интервале усечения (—Г/2, Г/2) [127], т. е.
№ { t ) = |
|
Ч |
E(fe) |
1 К(ÖJ t |
|
У |
£ Г е |
|
(4.8) |
||
°Т |
|
|
|
|
|
где |
Г/2 |
|
|
|
|
|
|
|
|
||
—Г/2 |
|
|
|
||
s |
|
|
Г/2 |
|
|
Л/ |
|
і К (00 t W |
|
||
|
|
Ф j' i S (< ) (''“ |
* |
||
n=г=—ЛN |
|
|
|||
|
|
-Г /2 |
|
||
ЛГ |
|
— i к (Oo |
|
|
|
N |
|
n .(*) |
|
||
~ |
|
|
|
(4.9) |
|
= l i e |
|
gx0 (rt Cöo) |
|||
n——NS |
|
|
|
|
150
и Т= 2 п / ш о — период функции l(Tk) (/), тождественно равной |
(t) |
||
на интервале (—Т/2, |
Т/2) |
гари .веек значениях к. |
|
Коэффициенты |
в |
разложении (4.8) являются случайными |
комплексными величинами, так как функции b , T(k ) ( t ) различны при разных к в соответствии со сменной реализацией рассматриваемо
го случайного процесса. |
а для коэффициентов |
раз- |
Очевидно, что Ш ) = 1іт |
||
Т-+00 |
|
|
ложения (4.8) справедливы соотношения (см. [127]): |
|
|
1 іш т 1 { і к * ) } = 0 , 1 і т т 1 { [ ^ * )]2} |
= 0 |
(4.10) |
Т-*со Т-+со
В соответствии с соотношениями (4.4), (4.5) и (4.8) результат нелинейного преобразования импульсного случайного процесса мо жет быть представлен в виде детерминированных функций време ни и некоторого числа случайных параметров, определяемых ко эффициентами разложения преобразуемого процесса в функцио нальный ряд. Для каждой конкретной реализации преобразуемо го импульсного случайного процесса набор коэффициентов разло жения представляет собой набор детерминированных величин.
Представление импульсного случайного процесса в виде ряда упрощает решения как общей задачи определения функции рас пределения вероятностей нелинейных безынерционных преобразо ваний, так и частных задач вычисления моментных функций и энер гетических характеристик таких преобразований.
Следует заметить, что определение спектрально-корреляционных характеристик нелинейных преобразований случайных процессов задача, как правило, более сложная, чем определение функций раопределения.
В заключение этого параграфа отметим, что в практике проек тирования и исследования систем обмена информацией возникают задачи анализа нелинейных преобразований не всего импульсного процесса в целом, а только отдельных параметров импульсов (ам плитуды, длительности, временного положения и т. п.). С учетом возможности разделения операторов таких преобразований на не линейные безынерционные и инерционные линейные задача ис следования статистических характеристик импульсных случайных процессов на выходе таких преобразователей еще более упро щается.
Отмеченные особенности исследования нелинейных преобразо ваний импульсных случайных процессов подробнее будут рассмот рены в последующих разделах этой и других глав.
4.2. ВЕРОЯТНОСТНЫЕ ХАРАКТЕРИСТИКИ ИМПУЛЬСНОГО СЛУЧАЙНОГО ПРОЦЕССА НА ВЫХОДЕ НЕЛИНЕЙНОГО ЭЛЕМЕНТА
Преобразование амплитуды и формы импульсов. Рассмотрим результаты воздействия импульсного случайного процесса, k-я реа лизация которого имеет вид
151
|<*>(0 = | ] |
lnk) и |
t - t p \ |
|
< h) J |
|||
П — — СО |
|
на безынерционный нелинейный элемент с характеристикой у — =\f(x). Вероятностные характеристики |(7) заданы функциями рас пределения вероятностей случайных параметров In, хп и tn, а от носительно функции f(x) положим, что она ограничена и непрерыв на со всеми своими производными для каждого значения х, т. е. при любом значении х = Хо представима рядом Тейлора:
/(*) = / Ы Iх ■<■») Іт (х)\*-*,+ -і- (X— *>) /|!| W l~«. +
или |
|
|
|
|
f(x) = Y^bpix — x0)p, |
|
|
(4.11) |
|
p= 0 |
|
|
|
|
где |
|
|
|
|
/(P)W|л:=л-0 |
_ 1 |
dp/W |
(4.12) |
|
P! |
(d*)p |
|||
Pi |
|
Значение аргумента f(x), при котором строится ряд (4.11), в электросвязи называемое рабочей точкой, выбирается, исходя из требуемого режима работы нелинейного элемента (выпрямление, детектирование, усиление, ограничение и т. л.) в соответствии со статистическими характеристиками преобразуемого процесса.
Коэффициенты разложения Ьр в ряд (4.11) при заданных ха рактеристике преобразования y=f(x) и положении рабочей точки постоянны. Очевидно, что разным характеристикам f(x) соответ ствуют разные наборы коэффициентов {Ьр}, р = 0, 1 , 2 , ... даже при одном и том же значении аргумента хо, и, наоборот, при за данной характеристике нелинейного элемента y=f(x) величины f(-^(x)/x0 и, соответственно Ьр, определяются выбором рабочей точ
ки х0 (точки разложения f(x) в ряід |
Тейлора). Пусть, например, |
характеристика нелинейного элемента |
описывается функцией |
X |
|
У — f(x) — А j е~ 2*/2 dz, z — (х — а)/о. |
(4.13) |
На рис. 4.5 изображен график функции (4.13) для нескольких зна чений параметра о. Как известно [53, 76], функция (4.13) непре рывна вместе со всеми своими производными на интервале —о о < < х < о о . Ряд Тейлора для этой функции имеет вид f(x) =
152
—оэ
характеристики
- a \ 2
S'bp(x—xo)p,где коэффициенты fr — j L .
£o |
P Pl dxp er У 2л |
x = x a
и ершима ют разные значения ів разных точках оси абсцисс OM. рис. 4.5):
Очевидно, что положение рабочей точки влияет также на ха рактер аппроксимации іf(x), а значит и на характер преобразова ний (рис. 4.5).
По отношению к преобразуемому сигналу \ ( t ) рабочая точка обычно выбирается так, чтобы отклонения §(£) от х0 имели как можно меньшую вероятность, т. е.
Р { [ | ( 0 — *о] = е) = min
Этому условию, как известно, удовлетворяет математическое ожидание случайной функции ті{І^(і)}. Практически, рабочая
4) Предположение о непрерывности y = f ( x ) и всех ее производных в любой точке оси X не онижает общности дальнейшего анализа, так как полученные при этом результаты могут быть распространены и на функции, имеющие конеч ное число разрывов [112]. Разрыв непрерывности функции f ( x ) или ее произ
водных при этом могут быть устранены путем введения обобщенных функ ций |[:127].
153
точка х0 совмещается с постоянной составляющей преобразуемого процесса:
|
Г/2 |
____ |
|
х0 = 1іш іг |
f |
l {k)(t)dt = l {k){t) |
(4.14) |
Г-оо Т |
J |
|
|
|
-Г /2 |
|
|
Если процесс эргодический, то ^ k)( t ) — rrii{l,(-h'>(t)}, а положение ра бочей точки определяемое соотношением (4.14) неизменно во вре мени. Для нестационарных или неэргодических процессов вели чина х0 также может выбираться равной математическому ожида нию нестационарного процесса:
x0(t) = m1{ l ik) (О) = ml%{t) |
(4.15) |
или, в соответствии с (4.14), равной временному среднему неэргодического процесса. В этом случае положение рабочей точки будет разным в различные моменты времени для нестационарного, но эргодического процесса или разным для различных реализаций неэргодического процесса, и х0, выбранном в соответствии с (4.14).
Соответственно изменениям х0 будут меняться и значения про изводных функций tf(x) и коэффициенты Ьр, см. ф-лу (4.12), также будут детерминированными или даже случайными функциями вре мени.
Если Хо не зависит от времени, ход процесса х0 на выходе не линейного элемента определяется множителями вида (х—хо), где Xo=mi{gk)(t)}. Без потери общности примем Хо = 0. Тогда k-ю реа лизацию процесса на выходе нелинейного элемента можно пред ставить в виде:
г |
ч о = |
£>(*)]' |
|
Р- 0 |
|
|
0 Lfi=Еoo |
(4.16) |
” V |
r(‘) |
Заметим, что если импульсы преобразуемого процесса не пе рекрываются, то
0 при п Ф /,
(4.17)
Поэтому ф-лу (4.16) можно переписать в виде
5(ft) (0 = 5 ]б р |
(4.18) |
р=0 / |
|
154
или, меняя очередность суммирования по индексам р и п, обозна чая
|
|
Р |
|
|
(4.19) |
запишем |
|
|
l ik)(t)= |
&k) |
(4.20) |
|
П= — со |
■VI Г |
В соотношениях |
(4.19), (4.20) учтено, что в нелинейных безы |
нерционных системах реакция в любой момент времени определя ется только значением воздействующего процесса в тот же момент времени t.
Таким образом, импульсный характер воздействия сохраняет ся и інѳперѳкірыіваюіщ'иеся иміпульсы на входе остаются также інѳперекрывающимися на выходе.
Импульсы процесса £(7) на выходе нелинейного элемента опи сываются взвешенной суммой соответствующих функций вида [u(t)]p, взятых с весовыми коэффициентами Ьр (4.19). Каждый п-к импульс к-к реализации процесса t,(t) можно представить как ре
зультат умножения |
детерминированной |
функции v(t), |
которая |
|
тождественно равна нулю вне интервала |
на величину І,\к) |
|||
сдвига по оси |
времени на величину |
и изменения |
масштаба |
|
времени в т ^ ! |
раз. |
|
|
|
Таким образом, |
|
|
|
|
£ (fe)( 0 = |
|
t — t(k) |
|
(4.21) |
l (nk ) V |
%n n C |
|
||
|
|
T(*) |
|
|
|
|
C |
|
|
Учитывая, что функция u(t), определяющая форму импульса на входе нелинейного элемента, в какой-то момент времени t*
xhk)) принимает максимальное значение, равное единице, из ф-лы (4.19) найдем, что
f(k) |
(4.22) |
ЪП |
|
р = 0 |
— 1> *1пк) — 0> х„к) = 1> |
С другой стороны, если принять, что |
|
то |
|
оо |
|
* ( 0 = Х М “ (0 ]Р. |
(4-23) |
р = 0 |
|
Вероятностные характеристики амплитуды импульсов процесса на выходе нелинейного элемента в этом случае легко определяют-
155
ся, если известны вероятностные характеристики амплитуды импульсов процесса на входе. Так, моменты случайной величины £„ на основании (4.22) равны:
ПІг - ' ; Ш. ([ І “ ]'} = /« і р Д Й * '] ' |
(4.24) |
р=0 |
|
Видно, что для определения момента г-го порядка случайной величины £„ необходима одномерная функция распределения ве роятностей (по меньшей мере моментов до рг-го порядка).
Если известна двумерная функция распределения вероятностей случайных амплитуд п-іго и /-го импульсов процесса £(£), то мо жно определить двумерные моменты случайных амплитуд п-го и /-но импульсов процесса £(f) на выходе нелинейного элемента:
mrq { Ük) l) |
Й‘7 |
у ь Д б Г У |
(4.25) |
LP=O |
J |
L i=o |
|
В частности, среднее значение, дисперсия и коэффициент кор |
|||
реляции амплитуд импульсов процесса t,(t) согласно |
(4.24) и |
||
(4.25) равны: |
|
|
|
оо |
|
|
(4.26) |
а-. = ті { Ük)} = У,Ьртр{ lhk)}, |
|
||
р=о |
|
|
|
{ Ң*1) - <4 _ £ £ |
Ь,.Ь, {т„ { 6!»} - |
|
|
р= 0 1=0 |
|
|
(4.27) |
|
|
|
(4.28)
Ѵ м 2{ ^>}М 2 {£<*>}
где
(4.29)
р= 0/=0
Таким образом, для определения энергетического спектра ус редненной корреляционной функции нелинейного преобразования импульсов случайного процесса необходимо не только среднее зна чение, дисперсия и коэффициент корреляции амплитуд импульсов на входе, но и моменты более высоких порядков.
Очевидно, что, зная моменты функции распределения вероят ностей случайной амплитуды, можно построить по известным пра вилам [53, 81] ряд, сходящийся к самой функции распределения или к характеристической функции.
Рассмотрим теперь способы определения формы импульсов на выходе нелинейного элемента. Как было показано выше, форма
156
импульсов на выходе определяется как взвешенная сумма слагае мых вида {u(t)]P. Для нахождения функций [u(t)]p можно предста вить u(t) в виде экспоненциального ряда (4.4). Это позволяет представить {u(t)f в виде
Іи (О]2 = « 2 (0 = |
V £ |
АКі Ак, А****) 1 |
|
|
|
|
к,=0 к2=0 |
|
|
Тогда более высокие степени функции u(t) |
можно записать как |
|||
“•<*) =к ,£— 0 |
к £, — 0 |
к,=ОІ !• к А ‘ к А |
1 ) , ке(* , к ' і С " |
|
“ Р (0 = П |
X |
Л ^ е х р ^ О . |
(4-30) |
|
К=1 |
кх =0 |
|
|
Для исследований функций up(t), представленных рядом (4.30), можно использовать преобразование Лапласа [25]. В самом деле, обозначая преобразование Лапласа от f(t) как
F(s) = ] n f ) e - Stdt
о
в соответствии с (4.29) и (4 30), найдем
СО |
со |
|
|
|
(4.31) |
|
ЕЕ |
А* А* |
|
Л « Л * |
|
|
[»-(**.+ 8*,)] |
|
|
к,=0 /с,=0 |
|
Заметим, |
что представления (4.30), (4.31) и свойства преоб |
разования Лапласа [25] позволяют анализировать не только ре зультаты воздействия импульсных случайных процессов на безы нерционной, но и на многие виды инерционных нелинейных систем [104, 127].
В виде ряда Фурье может быть представлена любая функция f(t), удовлетворяющая на промежутке OsS^s^T условиям Дирих ле. Учитывая определение импульсного случайного процесса и ус ловие ограниченности характеристики нелинейного элемента f(x), функция v(t) аналогично (4.4) может быть представлена рядом
Фурье на том же, что и u(t) |
интервале Т, причем Т ^ х п, т. е. |
v(t) = у Акеіка°(\ |
(4.32) |
К— — СО |
|
157
где
Г/2
A . - T |
f » ( O e - " - ' . |
|
(4.33) |
|||
|
* |
J |
|
|
|
|
|
—7/2 |
|
|
|
|
|
|
Рассмотрим связь |
коэффициентов разложения в ряд |
Фурье |
|||
функции ц(Т) и и(Т). |
|
можно записать, что |
|
|||
Л |
На основании (4.33) и (4.23) |
|
||||
оо |
|
|
|
оо |
|
|
= |
* |
Г/2up(Oe-‘““»<Ä=5]6^Ä/,> |
(4.34) |
|||
|
ЛИ* |
J |
|
ши |
|
|
|
р = 0 |
|
- Г /2 |
|
р= 0 |
|
где |
Г/2 |
|
|
|
|
|
|
|
|
|
|
|
|
Л р |
= ^ г |
f |
ttp ( O e " “ |
“ , , Ä . |
|
(4.35) |
|
* |
и |
|
|
|
|
—Г/2
Для большей наглядности заменим непосредственное определе ние величин ЛкР, по ф-ле (4.34) возведением в р-ю степень раз ложения (4.5). Тогда
u2(t) |
V |
скеіка°* |
|
ZJ |
K |
|
/C=—CO |
|
u3(t) |
I |
i « (Oo* |
|
2 |
oo |
со |
|
І (*,+*„) (Ooi |
|
|
И |
Ъ |
Ск'Ck* |
|
|
3 |
/ C j = — oo K 2= — oo |
|
|
|
|
oo |
oo |
У |
cKl сКг сКзе |
' («і4«2+Кз) ®о< |
|
|
У |
У |
|
/С- — — 00 Ко= — 00 / С , = — 00
(4.36)
р |
со |
оо |
up(t) = X |
= у |
у ... у , KlcK2...S x |
Xexp I і £ /c r coo4
г = 1
Соотношения (4.36) позволяют сделать важный для практики вывод, а именно, в результате нелинейного преобразования в соста ве сигнала на выходе появляются новые гармонические составляю щие (так называемые комбинационные частоты). Пусть, например, в разложении (4.5) только со, си Съ С - с_ 2 отличны от нуля. Тогда, выбирая рабочую точку х0= с0 и считая, что Сі= с_і, с2= с_г, пред ставим (4.5) в виде
и (0 — с0 + 2 Сі cos CÖ0 t 2 c2 cos 2 co01,
а возведение {w(Y)—x0] в квадрат дает
и2(t) = |
4 cf cos2 (i)01 + 4 cf cos2 (j)01-f 8cx c2 cos a>0t cos 2 oo01 = |
|
= |
( cf + Сг) + |
4 Cj c2 cos coo ^ + 2 cf cos 2 co01 + |
-f 4 q c2 cos 3 co01 + |
2cfcos 4 co01. |
158
Таким образом, уже квадратичное преобразование импульса u(t) привело к появлению новых гармоник Зсо0 и 4(0<ь а также к уве личению уровня постоянной составляющей входного сигнала.
В более общем случае, при бесконечном числе слагаемых в (4.5), после нелинейного преобразования мощность перераспределяется, так что общая ширина спектра выходного сигнала увеличивается по
сравнению с шириной спектра входного сигнала. |
|
Для определения величин А кр |
заметим, что экспоненциальный |
р |
|
множитель ехр{і У Krmt}, определяющий частоту колебаний функ- |
|
г = 2 |
при изменении кг(г—\, 2 , . . ., |
ции uP(t) и множитель ехр{і/соѴ} |
р) и к изменяется дискретно. Поэтому при каждом р величина А кр
будет |
определяться |
как |
сумма |
произведений |
коэффициентов |
||||||||
Скг (Г= Ь 2, |
..., |
р), |
соответствующих |
таким |
наборам |
кг(г= 1 , 2 , |
|||||||
. . ., р), |
для которых |
выполняется |
|
р |
|
кг= к. |
Таким об- |
||||||
равенство V |
|
||||||||||||
разом, |
|
|
|
|
|
|
|
|
г=1 |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
о& |
|
|
|
|
|
со |
|
|
|
|
|
|
= |
У ЬрАйр- Ь 0~\- Ьі с0 + |
Ь2 у |
с*, С—К1А~Ь3у у |
у |
cKi cKt cKj+ ... |
||||||||
|
р =О |
|
|
|
|
« , = —«> |
|
К , к г |
к , |
|
|
||
|
|
|
|
|
|
|
|
|
|
к,+к2+к3=0' |
|
||
|
|
|
|
П |
ч |
+••• |
|
|
|
|
|
||
|
|
* 1 |
Kt |
Кр |
1 |
|
|
|
|
|
|
|
|
где |
р |
|
|
|
|
|
|
|
|
|
|
|
|
V кт = О |
|
|
|
|
|
|
|
|
|
|
|||
|
г= 1 |
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
= |
|
|
- Д |
Г К , |
|
|
|
(4.37) |
||
|
|
|
|
|
|
|
|
||||||
|
р—0 |
|
|
р—0 |
К і к2 |
|
Кр |
Г— 1 |
|
|
|
|
|
где |
Р |
|
|
|
|
|
|
|
|
|
|
|
|
у кг = к, |
|
|
|
|
|
|
|
|
|
|
|||
/■ =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ьр A-кр = £ /ѵ У |
1 J- |
• • |
I J П Ск<-. |
|
|
|
|
||||
|
р = 0 |
|
р= 0 к, к, |
|
кр г— 1 |
|
|
|
|
||||
|
Р |
/сл = — /с. |
|
|
|
|
|
|
|
|
|
||
где + у |
|
|
|
|
|
|
|
|
|
||||
|
г=1 |
|
|
|
|
|
|
|
|
|
|
|
|
Причем |
суммирование в |
(4.37) |
ведется по всем кг от — о о до |
||||||||||
|
, удовлетворяющим |
равенству |
р |
кг— к. |
|
|
;... |
||||||
- ) - о о |
у |
|
|
/•=1
159