
книги из ГПНТБ / Марченко Б.Г. Метод стохастических интегральных представлений и его приложения в радиотехнике
.pdfдал ему исчерпывающее объяснение, которое с несущественными доработками принято и теперь, В Шоттки [ 100]. Подобно дробовому эффекту могут быть построены модели других шумов, описывающих ся линейными случайными процессами. Некоторые из них будут рассмотрены ниже.
§ 2. Тепловые и полупроводниковые шумы. Реверберационные процессы и шумы торошения льдов
Остановимся кратко на рассмотрении некоторых наиболее часто встречающихся в приложениях шумов с целью установления при менимости для их описания математической модели линейных случайных процессов. К сожалению, не во всех случаях удается най ти прямое обоснование и полностью проследить связь для всех пара метров линейного процесса с исходными физическими предположе ниями так просто и полно, как это было сделано для дробового эф фекта-
Наряду с дробовым эффектом большое распространение в ра диотехнике и физике получила модель тепловых шумов, которая описывает флуктуации тока и напряжения в электрических цепях, обусловленные хаотическим движением электронов в проводнике. Хотя изучением этого вида шумов физики занимаются начиная с конца прошлого века, основополагающие результаты, удобные для применения в теории связи, были получены только в 1927—1928 го дах Джонсоном [91] и Найквистом [99]. Эти результаты, начиная с 1928 г., неоднократно подтверждались в экспериментальных рабо тах многих авторов.
Для случая теплового шума существует множество различных моделей, приводящих к результатам, полученным Найквистом. Наи более часто встречаются две из них:
а) модель, которая исходит из общих принципов статистической механики, и
б) модель, использующая теорему о равномерном распределе нии энергий применительно к электрической или магнитной энер гии цепи.
Так как целью настоящей главы является установление воз можности представления рассматриваемых здесь процессов с помо щью стохастического интеграла вида (1), то нет смысла подробно останавливаться на рассмотрении указанных выше моделей теплово го шума. Исходя из теоремы Найквиста 199] и близости теплового шума к броуновскому движению, что было установлено многими авторами [35] с использованием модели б), сразу приведем некото рые доводы в пользу возможности представления теплового шума в виде (1).
Согласно теореме Найквиста для электрических цепей [11, 73,
120
99], спектральная плотность шумовой э. д. с. в малом частотном
интервале df |
определяется |
выражением * |
|
|
|
|||
|
|
|
S(J)df = AkTRp (I/ 1) df, |
|
(409) |
|||
где f — частота, гц; S (f) |
спектральная плотность теплового шу- |
|||||||
ма; k — постоянная Больцмана ^1,38-10 23 |
|
j ; Т — абсолютная |
||||||
температура, |
°К; |
R — со |
|
|
|
|
|
|
противление, |
ом; |
р (f) — |
|
|
|
|
|
|
фактор Планка: |
|
|
|
|
|
|
||
= Л _ |
Р(/) = |
—1 |
|
|
|
|
|
|
exp |
hf_ |
|
|
|
|
|
||
kT |
kT |
— 1 |
|
|
|
|
|
|
где h = |
6,62 |
10 |
(410) |
|
|
|
|
|
34 дж X |
|
|
|
|
|
|||
X сек — постоянная План |
|
|
|
|
|
|||
ка . Обоснование выражения |
|
|
|
|
|
|||
(409) содержится в [11]. |
|
|
р (/), определенной со |
|||||
На рис. 6 приводится график функции |
||||||||
гласно (410) для значений /, удовлетворяющих условию 0 < |
< 4 . |
|||||||
При |
нормальной комнатной температуре |
в диапазоне |
частот, |
|||||
используемых в радиотехнике [73], имеем |
А|П |
1, поэтому в ука- |
||||||
|
kT |
|||||||
занном диапазоне р |
и выражение (409) в этом случае упро- |
|||||||
щается следующим образом: |
|
|
|
|
|
|||
|
|
|
S(f)df ^AkTRdf. |
|
|
(411) |
||
В работе [35] при помощи метода моментов в предположении |
||||||||
линейной внешней нагрузки и с использованием модели б) |
доказы |
вается, что тепловой шум имеет нормальное распределение, не за висящее от времени.
Теорема Найквиста и тот факт, что тепловой шум является ста ционарным нормальным процессом, позволяют сразу установить его принадлежность к классу линейных процессов и написать яв ное выражение для интегрального представления теплового шума.
Рассмотрим следующие три случая.
Случай 1. Если, исходя из приближенной формулы (411), пред положить, что тепловой шум содержит в полосе частот (—оо, оо) постоянную спектральную плотность, т. е. распространить эту формулу на бесконечную область частот, то получим, что тепловой шум является «непрерывным гауссовым процессом с независимы ми значениями» или, более точно, представляет собой обобщенную производную [15, 29] от винеровского процесса:
Z(t) = o ^ P - , |
(412) |
* В работе [ 11 ] в (409) вместо р (// |) приведено р (/), что не позволяет распро странить эту формулу на область / < 0.
121
где о2 = AkTR и w (t) — винеровский процесс. Так как g (t) яв ляется обобщенным случайным процессом, то выписать его плот ность распределения в явном виде в классе обычных функций не представляется возможным, поэтому под гауссовым процессом
подразумевается тот факт, что процесс |
[w (t -f- Д^) — w(f)] |
для |
|||||
At > 0, |
который с |
уменьшением At |
приближается к обобщенной |
||||
производной процесса ow (t), |
имеет |
плотность |
распределения |
||||
|
Раі {х) = |
|
ехр |
|
|
Д *> 0. |
|
Корреляционная |
функция |
этого |
процесса |
определяется |
так: |
||
М ! |
\w (t -f At) — w (*)] [w(t + |
s + A ) — w (t + s)J [ = |
|
||||
|
|
AP (A* — Is I) ПРИ |
И < А /, |
|
(413) |
||
|
|
О |
при \s\^> At. |
|
Следовательно, его спектральная плотность существует и опре деляется следующим образом:
A t |
|
|
,, , |
|
1 — cos 2nfAt |
||
—/гл/s аа . . . |
. |
|
|||||
« = ! |
AP (A t - \s \) d s |
= |
---- |
- |
° |
||
— A t |
|
|
|
|
|
|
|
При малых At эта |
спектральная |
плотность |
|
близка к а 2, т. е. |
при малых At функция 5д, (/) приближается к 5 (/), определенной согласно (411) при /£ (—оо, оо). Таким образом, хотя процесс £ (t), определенный согласно (412), и не существует в обычном смысле,
соответствующая ему |
спектральная плотность, полученная |
из |
по |
||||
следнего выражения |
при |
уменьшении |
At, ведет |
себя так, |
как |
||
будто бы |
I (t) существует |
и является |
стационарным гауссовым |
||||
процессом |
с постоянной спектральной плотностью о2. |
(412), |
|||||
В радиотехнике процесс g (t), определенный |
согласно |
обычно именуют «белым шумом». Для этого процесса Р (|£ {t)\~>A} = = 1 при любом положительном А . Физически осуществить процесс g (/) невозможно хотя бы потому, что он обладает бесконечной мощ ностью .
Таким образом, предположение о постоянстве спектральной плотности S if) на всем интервале частот/ £ (—оо, оо) для реального
физического теплового шума оказывается неприемлемым. |
еди |
|
Случай 2. Предположим, что в (409) функция р |/| равна |
||
нице в интервале(— /0, /0), f0> 0, равна нулю при |/ |l > f0, а |
зна |
|
чения в точках 1/1 = /0 |
могут быть любыми, но конечными. |
|
Не останавливаясь |
на вопросе единственности линейного |
ин- |
122
тегрального представления, покажем, что тепловой шум может быть описан в рассматриваемом случае выражением
V (f) — о \ |
sin 2я/о (t— т) |
dw (т), |
(414) |
я (t— т) |
где а 2 = 4kTR. Для этого достаточно рассмотреть спектральную плотность, соответствующую процессу (414) и сравнить ее с исход ной. Корреляционная функция процесса (414) согласно (298) опре деляется так:
со |
sin 2nf0t |
sin 2я/0 (t+ s) |
|
sin 2n,f„s (415) |
• (S) = ° 2 J |
dt |
|||
nt |
я {t-j- s) |
|
|
Выражение для спектральной плотности процесса (414) с учетом [19], 3.741, (2), получаем в следующем виде:
|
а2 |
при |
|/ |< / о , |
S (/) = а2 j sin^ f o e~i2nfs ds |
-т р |
при |
|/1 = /о, |
|
О |
при |
| / | > / 0. |
Таким образом, центрированный гауссов процесс (414) имеет спектральную плотность, постоянную в полосе (—/0, /„) и равную нулю вне этой полосы, а следовательно, может быть принят в ка честве интегрального представления теплового шума. Однако такое представление все же остается очень идеализированным, так как его ядро является физически нереализуемым.
Итак, тепловой шум в рассматриваемом случае описывается слу чайным процессом, допускающим представление (1). Заметим, что процесс I (t), определенный согласно (412), также является в неко тором смысле линейным, так как формально может быть выражен в виде предела процесса (414) следующим образом:
СО |
|
5(0 = а $ Пгп sin* f V - .$ -dw(T). |
(416) |
Процесс (416) можно рассматривать как элемент замыкания класса линейных процессов с использованием предельного перехода в их ядрах. Предел в (416) в обычном смысле не существует и под разумевается в обобщенном смысле, а его ядро не удовлетворяет условиям, наложенным на ядро (1), так как является дельта-функ цией (не существует интеграл от квадрата).
Случай 3. В наиболее общем случае, когда спектральная плот ность теплового шума определяется выражением (409), он может быть описан линейным процессом (1) по винеровскому:
CQ |
|
о0(0 = ц ] ф0(г — x)dw(x), |
(417) |
123
где функция ф0 (і) определяется из уравнения (65), которое в рас сматриваемом случае принимает вид
|
|
и |
С |
p—tZKfsI |
f I |
|
|
|
Фо(*)Фо(* + |
«)Я = |
j |
(h \f \ \ |
|
~ dt- |
(418) |
||
|
|
|
|
exp \~kT~) ~ |
1 |
|
||
Если дополнительно предположить, что функция ф0 (t) |
является |
|||||||
физически реализуемой, то уравнение (418) |
превращается |
в урав |
||||||
нение типа Винера — Хопфа [62]: |
|
( 2nkT |
\ |
|
||||
оо |
|
со |
|
|
||||
|
|
|
Xcos ( |
т |
XSI |
(419) |
||
J Фо (0 Фо (t + |
S) dt = |
j |
-------- г г - ----- dx, |
|||||
О |
|
о |
|
|
|
|
|
|
которое с учетом [17], 3.951, (5), можно записать еще и так: |
||||||||
оо |
кТ / |
h2 |
|
л2 |
, |
2 2n2kTs |
||
Г ... .. , . ,, |
|
|||||||
J *Ро (0 Фо (t + S) dt — |
( |
8n2fe2T2s2 |
2 cosech |
h |
|
|||
о |
|
|
|
|
|
|
|
|
Так как согласно (419) и с учетом [19], |
3.411, |
(1) |
|
|||||
|
О |
|
|
|
|
|
|
(421) |
|
|
|
|
|
|
|
|
то ядро в (417) интегрируемо в квадрате, т. е. и в самом общем слу чае тепловой шум является линейным случайным процессом, допу скающим представление вида (1).
Выражение (421) с учетом представления (417) позволяет опре делить дисперсию теплового шума (аналогично [73])
h(0) = Mvl]{t) = » 1,8928 • 1 0 (дж ■сек ■ом), (422)
а выражение (298) с учетом (417) и (420) определяет его корреляци онную функцию
h(s) = |
fiR |
2n2k2T2R |
cosech2 |
2n2kTs |
(423) |
2JI2S2 |
h |
h |
На рис. 7 приводится график нормированной корреляционной функции теплового шума, которая с учетом (422) и (423) определя ется следующим выражением:
ft(s) |
3ft.2 |
2n2kTs |
|
|
h |
|
|
||
А(0) |
4n*T2k2s2 |
2n4Ts |
|
|
7,2907 • IQ"24 |
0,4114 • 10127’s |
(424) |
||
|
T2s2 |
sh (0,4114 • 1012rs) |
||
|
|
|||
При этом диапазон |
изменения |
s выбран из |
условия |
0 < s < |
20 • ІО-12 |
|
|
интервале |
корреля |
•< -----j,-----. Как видно из графика, на этом |
ционная функция теплового шума убывает достаточно быстро-
124
Рассмотренные выше три случая позволяют сделать заключе
ние о пригодности математической модели линейных |
процессов |
|
для достаточно полного |
описания реальных тепловых шумов. |
|
Остановимся кратко |
на простейшей линейной модели |
полупро |
водниковых шумов. При анализе этих шумов обычно проводят раз личие между тепловым шумом, дробовым шумом и фликкер-шумом [11 ] Тепловой шум в полупроводниках обусловлен случайным дви жением носителей тока Он аналогичен тепловому шуму, возника ющему в проводниках, который, как было установлено выше, мож но описать интегральным представлением (1). Дробо вой шум в полупроводниках наиболее характерен, его спектр в области нижних частот равен постоянной ве личине независимо от ча стоты. По этому признаку он подобен дробовому шуму в электронных лампах, кото рый рассматривался ранее и, как было выяснено, мо жет быть описан интеграль
ным представлением (1). Фликкер-шум характеризуется своеобраз
ной |
спектральной |
плотностью |
[11, |
12, |
59] |
|
|
|||
|
|
|
|
5 (со) = с2со~а, со 6(0, |
оо), |
|
(425) |
|||
где |
а = 1 |
+ |
е, 8 — некоторая |
постоянная |
величина, 0 |
<; s <С 1, |
||||
со = |
2яД |
с — некоторая постоянная. |
По |
спектральной |
|
характе |
||||
ристике этот |
шум |
подобен фликкер-шуму |
в электронных |
лампах. |
Фликкер-шум в полупроводниках значительно превышает уровень теплового и дробового шумов. При его анализе с помощью линей ной * теории обычно делают предположение о дискретности электрон ного заряда подобно тому, как это делалось в § 1 настоящей главы, т. е. предполагают, что группы зарядов (в § 1 рассматривались не группы, а отдельные электроны) появляются в случайные моменты времени и удовлетворяют предположениям а) и б) из § 1 настоящей главы. Следовательно, моменты их появления распределены по за кону (121) Так как эти заряды носят случайный характер и имеют одинаковые функции распределения, то, следуя рассуждениям, приведенным в § I, придем к представлению фликкер-шума в виде
(403)Однако теперь функция aNe. (t) будет представлять собой
импульс, связанный с зарядом группы, состоящей из N синхрон
но |
возникающих носителей, каждый из которых имеет заряд q = |
|
со |
= |
j е(х) dx, aN — Nq, а не заряд, равный заряду электрона, как это |
|
о |
* Существует и нелинейная теория, развитая Моррисоном. Ее описание и ссылки на литературу даются в работе [11].
125
было в дробовом эффекте (е (т) — форма отдельного импульса). Следовательно, полупроводниковый фликкер-шум можно предста вить в виде
оо |
|
£(0 = J e (f--T ) cfatj (т). |
(426) |
о |
|
Теперь остается только определить е (т) таким |
образом, если это |
возможно, чтобы спектральная плотность процесса £ (t) описыва лась бы согласно (425).
Корреляционная функция процесса (426) с учетом (405) имеет следующий вид:
оо
|
h (s) — |
J е (т) е (т + s) dxf |
(427) |
|
со |
|
о |
|
|
|
|
|
|
|
где а% = J хЧЗ~(х), 3 (х) — функция распределения зарядов |
груп- |
|||
о |
а е (т) удовлетворяет |
условию физической |
реали |
|
пы носителей, |
||||
зуемости . |
|
|
|
|
Согласно теореме Винера — Хинчина [76] имеем |
|
|||
|
оо |
оо |
оо |
|
S (со) = |
j h (s)e~“°sds — Xa% j |
J e (r) e(r -f s) elasdsdx. |
(428) |
|
— co |
— oo 0 |
|
||
Функция |
|
00 |
|
|
|
|
|
|
|
|
Sx(too) = |
J e (s)е~ш ds |
(429) |
|
|
|
о |
|
|
аналитическая в правой полуплоскости. С учетом (429) выражение (428) принимает вид
S (со) = ка% Sx (іа) Sj (— tco), |
(430) |
т. e. последнее выражение представляет собой факторизацию спект ральной плотности 5 (со) на прямой. і
Из (429)
|
ОО |
■'* |
e(0 = - s r |
S 5х (ісо)еш <ѣ. |
(431) |
Следовательно, |
—00 |
|
00 |
|
|
|
|
|
5і(0) = |
J е(т) dx = q, |
(432) |
|
о |
|
поэтому должно выполняться условие |
|
|
2 |
5(0) |
(433) |
ан |
|
Если исходить из спектральной плотности, определенной вы ражением (425), то согласно (433) средняя интенсивность импуль
126
сов носителей равна бесконечности. Кроме того, |
функция |
(425) |
|
не интегрируема на интервале [0, оо). |
Все это |
не согласуется |
|
с физикой явления. Значит, спектральную |
плотность нельзя |
опре |
делить согласно (425). Заметим, что (425) получено в результате обработки экспериментальных наблюдений, а поэтому является од ним из возможных аппроксимирующих выражений для спектраль ной плотности и не исключает других подобных выражений для ее
аппроксимации. Исходя из |
того, что область определения |
(425) |
|
не включает в себя точку |
to = 0 и ее некоторую б-окрестность, |
||
совершаем в (425) |
замену переменной согласно соотношению |
со2 = |
|
= сох + б2. гДе 0 |
< б <С 1. |
В результате получаем еще одну ап |
проксимацию спектральной плотности фликкер-шума в полупро водниках:
а |
|
S (со) = с2 (б2 + со2)- ~ , (0 б [0, оо), |
(434) |
которая при малых, соответственно выбранных значениях б с на перед заданной степенью точности на любом конечном интервале из области со £ (0, оо) может аппроксимировать значения спектраль ной плотности, убывающие по гиперболическому закону, при этом постоянные с и п определены как в (425). Величину б при известных с и а можно определить экспериментально, путем измерения зна чения 5 (0) и вычислений по следующей формуле:
9 |
__1_ |
|
б = с а S |
а (0), |
(435) |
полученной из (434). В случае спектра, определенного выражением (434), имеем
S ( < o ) = t a £ ------
|
|
|
|
(б + |
|
«а) 2 |
(б — «а) 2 |
|
|
|
|||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
с, = |
|
|
— cqS |
г (0) |
|
|
|
||
|
|
|
|
V'-Ы{ |
|
|
|
|
|
|
|
|
|
При |
■ |
0 |
с |
учетом [6], |
|
3.2, (3), |
из последнего выражения |
||||||
находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
e(t) = |
- |
се~ ы |
1(0 |
|
|
gc exp { -[5 (0 ) с~2] |
а |
, |
(436) |
||||
|
|
1—- |
|
|
|
|
|
1—- |
|||||
|
|
|
|
|
|
|
|
|
|
|
|||
|
Ѵ Ч г ( у ) |
|
|
|
|
(0)Г |
т) |
|
|
|
|||
где 1 (/) — единичная функция |
|
Хевисайда: |
|
|
|
|
|||||||
|
|
|
|
1 (t) = |
[0 |
при |
^ с |
0, |
|
|
|
(437) |
|
|
|
|
|
|
, |
при |
_ |
п |
|
|
|
||
|
|
|
|
|
(1 |
I > |
0, |
|
|
|
|
||
которая |
терпит |
разрыв при t — 0 и в этой |
точке |
содержит |
пре |
||||||||
дел справа. |
|
|
|
|
|
|
|
|
|
|
|
|
127
Нетрудно видеть, |
что функция е (t), определенная согласно |
(436), отличается от |
функции Лагерра нулевого порядка только |
постоянным множителем, поэтому для нее выполняется условие (432).
Следовательно, |
в рассматриваемом случае форма импульса |
е (t) определяется |
спектральной плотностью S (со) с точностью до |
постоянного множителя независимо от вида cf (x).B свою очередь спектральная плотность не накладывает никаких существенных ограничений на распределение адт-суммарного заряда группы, за исключением условия (433), благодаря которому нам пришлось пе ресмотреть аппроксимацию спектральной плотности и придти к ее новому виду (434). Заметим, что вопрос, какая из аппроксимаций (425) или (434) лучше, подлежит экспериментальной проверке. Таким образом, исходя из физической линейной модели шума в по лупроводниках, приходим к заключению, что он также может быть описан интегральным представлением (1).
Рассмотрим кратко линейную модель шумов, обусловленных рассеянием звука на неоднородностях морской среды, именуе мых реверберационными шумами. Теоретические результаты пер вой главы позволяют описать при помощи характеристических функционалов поле реверберационной помехи. Рассмотрим вначале модель и свойства реверберационной помехи в некоторой точке шумового поля. Следуя работе [63], будем различать объемную ре верберацию, реверберацию от слоя и граничную реверберацию, а основное внимание сосредоточим на рассмотрении модели объем ной реверберации, так как остальные виды не имеют принципиаль ных отличий. На практике часто приходится иметь дело с линейной комбинацией сразу двух или трех видов реверберации [63].
Будем считать, что рассеиватели в морской среде расположе ны дискретно. Эффективные сечения рассеивателей практически оказываются намного меньше расстояний, занимаемых излучае мыми сигналами в пространстве. Если предположить, что рас сеиватели в среде расположены статистически независимо, для рассматриваемой достаточно большой области среды средняя плот ность рассеивателей постоянна, то количество элементарных рассеян ных сигналов, приходящих в точку приема за время Ат, т. е. попав ших в интервал (т, т -f- Дт), подчиняется распределению Пуассона (121). Это можно обосновать аналогично обоснованию выражения (121), проведенному для случая дробового эффекта.
Обозначим через ak случайные амплитуды, а через tk — слу чайные моменты возникновения рассеянных сигналов. Если из лучаемый сигнал описывается функцией s (t) и каждый элементар ный рассеянный сигнал s (tk, t), пришедший в момент времени tk, зависит от состояния рассеивателя в момент прохождения через него
фронта излучаемой волны, то процесс в точке приема для |
фикси |
рованного момента времени t может быть записан в виде |
|
£ ( 0 = Ü a*s(**,g. |
(438) |
— ОС |
|
128
При помощи определенного в § 2 гл. I процесса щ (т), прира щения которого связаны с ak соотношением (402), реверберационный шум можно записать в виде
t |
со |
|
|
Ш = Ш = i s (т, 0 djti (т) = |
j |
s(x, ОЛтДт), |
(439) |
— со |
— со |
|
|
где предполагается, что s (т, f) — 0 при т > |
t. |
|
Таким образом, модель реверберационного шума близка к модели дробового шума, только реверберационный процесс в об щем случае носит явно нестационарный характер, что и сказалось на виде ядра в (439)
Как уже упоминалось в § 1 настоящей главы, при к ->■ оо про цесс, полученный в результате центрирования и нормирования про цесса (т), стремится к винеровскому процессу, поэтому и процесс (439) после центрирования и нормирования также стремит
ся к гауссовому в общем |
нестационарному процессу, |
который в |
этом случае можно записать так: |
|
|
|
СО |
|
V (t) = |
^ s (т, t) dw (т). |
(440) |
|
— со |
|
Обоснование этого выражения можно построить аналогично обос нованию соотношения (408).
Ядро |
процесса (439) может |
быть найдено из уравнения |
|
|
ОО |
|
|
|
[ s (т, 4 )s(r, |
t%)dx — h(tv 4), |
(441) |
|
— СО |
|
|
где h (4, |
4) — корреляционная |
функция процесса (439), |
которая |
может быть определена из теоретических предпосылок или экспе риментально.
Если в (406) и (407) функции е (4 — т) и е (і — т) везде заме нить соответственно на s (т, tk) и s (т, і), то получим многомерную характеристическую функцию и характеристический функционал реверберационного шума для некоторой точки шумового поля.
Для описания «-компонентного вектора реверберационной по
мехи |
(^ (t), |
t,2 |
(t), |
..., £„(4 ), воспользовавшись теоремой 6, можно |
|||
сразу |
привести |
выражение характеристического функционала, |
|||||
который с учетом (127) можно записать в виде |
|||||||
|
|
|
|
|
(gif Sä, |
, |
gn) = |
|
|
|
ОО |
ОО |
п |
|
|
|
— ехр |
}Я |
j |
J |
[ехр(гл: 2 |
S * ( T , |
f)dgk {t)) — 1 ]d3r (x)dx). (442) |
|
|
|
— оо — oo |
f t = l |
|
|
В заключение этого параграфа заметим, что дискретная линейная модель шумов торошения льда может быть построена путем рассуждений, аналогичных приведенным выше при обосновании
9 3 — 6 6 2 |
129 |