Файл: Уриг, Р. Статистические методы в физике ядерных реакторов.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 16.10.2024
Просмотров: 158
Скачиваний: 0
Построение обратной связи, необходимое для генерирования последовательности максимальной длины, было рассмотрено для широкого диапазона последовательностей. В некоторых случаях мо гут быть использованы несколько различных обратных связей, в других случаях может оказаться необходимым применять несколь ко полусумматоров для генерирования последовательности. В табл. 9.3 приведены обратные связи, которые требуют наличия только одного полусумматора для генерирования последователь ностей с количеством разрядов р до 33.
Одна из схем, не включенная в табл. 9.3, но весьма полезная в экспериментах с ядерными системами, состоит из 8-разрядного регистра сдвига, дающего последовательность с максимальной дли-
Сддиг
Рис. 9.10. Восьмиразрядный регистр сдвига.
ной 255Д временных интервалов. Она не представлена в табл. 9.3, поскольку не может быть выполнена с одной обратной связью. Одна из возможных схем для генерации такой последовательности представлена на рис. 9.10. Матрица для этой конфигурации имеет вид:
0 1 |
1 |
0 |
0 |
0 |
0 |
0 |
г |
|
|||||||
1 0 |
0 |
0 |
0 |
0 |
0 |
0 |
|
0 |
1 0 |
0 |
0 |
0 |
0 |
0 |
|
0 |
0 |
1 0 |
0 |
0 |
0 |
0 |
|
0 |
0 |
0 |
1 0 |
0 |
0 |
(9.23) |
|
0 |
|||||||
0 |
0 |
0 |
0 |
1 0 |
0 |
1 |
|
0 |
0 |
0 |
0 |
0 |
1 0 |
0 |
|
1о |
0 |
0 '0 |
0 |
0 |
1 0 |
Схемы с 12-, 13-, 14-, 16- и 19-разрядными регистрами сдвига возможны, однако они требуют более одного сумматора по модулю 2. 17- и 20-разрядные генераторы, использующие только один полу сумматор, обычно столь же полезны, как и 16- и 19-разрядные регистры сдвига, которые работают с несколькими полусуммато рами. 8-разрядный регистр сдвига упомянут здесь потому, что он очень удобен при проведении экспериментов на ядерных системах,
243
Т а б л и ц а |
9.3 |
|
|
|
|
|
Обратные связи для генерации m-последовательности с р |
|
|
|
|||
до 33 при использовании только одного полусумматора |
|
|
|
|||
|
|
Суммируемые разряды |
||||
2 |
3 |
1 И 2 |
2 |
и 3 |
||
3 |
7 |
1 И 3, |
||||
4 |
15 |
1 и 4, |
3 |
и 4 |
||
5 |
31 |
2 |
и 5, |
3 и 5 |
||
6 |
63 |
1 и 6, |
5 и 6 |
|||
7 |
127 |
1 и 7, |
6 и 7 |
|||
9 |
511 |
4 |
и 9, |
5 и 9 |
||
10 |
Г 023 |
3 и ю, 7 и 10 |
||||
11 |
2 047 |
2 |
и 11, |
9 и 1 |
||
15 |
32 767 |
1 и 15, |
|
14 И 15 |
||
17 |
131 071 |
3 |
и 17, |
|
14 и 17 |
|
18 |
262 143 |
7 |
и 18, |
|
11 |
и 18 |
20 |
1 043 575 |
3 |
и 20, |
|
17 и 20 |
|
21 |
2 097 151 |
2 |
и 21, |
|
19 и 21 |
|
22 |
4 194 303 |
1 и 22, |
21 |
и 22 |
||
23 |
8 388 607 |
5 |
и 23, |
|
18 и 23 |
|
25 |
33 554 431 |
3 и 25, |
22 |
и 25 |
||
28 |
268 435 455 |
3 и 28, |
25 и 28 |
|||
31 |
2 157 473 647 |
3 и 31, |
28 |
и 31 |
||
33 |
8 629 894 591 |
13 и 33, |
20 |
и 33 |
в которых информация в виде импульсов от детекторов накапли вается в многоканальном анализаторе, имеющем 256 каналов, или в которых проводятся 2 или 4 выборки на каждый сдвиг соответст венно в 512или 1024-канальных временных анализаторах.
Автокорреляционная функция m-последовательности. Одно из наиболее важных свойств максимальной последовательности регист ра сдвига состоит в том, что ее автокорреляционная функция до статочно близка к дельта-функции, а боковые полосы постоянны. Это постоянство позволяет учесть корреляцию боковых полос, если это оправдано условиями эксперимента. Автокорреляционная функция получена скорее путем рассуждений, чем математическим путем. Рассмотрим р-разрядный регистр сдвига с одной обратной связью, который генерирует максимальную последовательность (m-последовательность), имеющую период ZA, где Z определено уравнением (9.19)'. Поскольку каждый разряд регистра сдвига просто пересылает входной сигнал от предыдущего разряда, выход ные сигналы любых двух разрядов одинаковы, за исключением того что один выходной сигнал сдвинут на М , где k — число разрядов между двумя выходными функциями. Воспользовавшись свойством сдвига и суммирования m-последовательности, получим, что выход ной сигнал полусумматора, который подается на вход регистра
.244
сдвига и, следовательно, на выходы всех разрядов, является также /п-последовательностью.
Если положить, что в псевдослучайной переменной два уровня /«-последовательности равны +1 и —1, среднее значение перемен ной будет равно± 1/Z, причем знак определяется тем, какой уровень имеет больший временной интервал. Например, последовательность со 127 сдвигами с 63 положительными интервалами и 64 отрицатель ными интервалами имеет среднее значение— 1/127. Если среднее значение принять равным нулю, значения двух уровней в этом при мере должны быть равны +128/127 и —126/127.
Для иллюстрации автокорреляционной функции рассмотрим /«-последовательность с уровнями +1 и —1. Автокорреляционная функция имеет вид:
. |
z |
|
<Рхх№) = т |
2 xt xi+h, |
(9.24) |
где xt — значение /«-последовательности на г'-м интервале. Так как
xt и xi+h имеют |
значения +1 или —1, корреляционный процесс, |
|||||
согласно |
уравнению |
(9.24), |
превращается |
|
||
в определение |
числа |
интервалов, при ко |
Т а б л и ц а 9.4 |
|||
торых Xi |
и xi+h |
имеют одинаковый знак |
Таблица соответствия |
|||
(положительный |
или |
отрицательный), и |
||||
числа интервалов, в |
которых знаки раз |
для умножения |
||||
+ 1 и—1 |
||||||
личны. |
Табл. |
9.4 |
представляет собой |
|
||
таблицу |
соответствия |
для |
перемножения |
|
||
+ 1 или —1. Полное совпадение с таблицей |
|
|||||
соответствия для |
сложения |
по модулю 2, |
|
|||
приведенной в табл. 9.2, оказывается |
|
|||||
очень полезным при оценке автокорреля |
|
|||||
ционной |
функции. |
|
|
|
||
Процедура |
оценки автокорреляционной |
|
функции с помощью уравнения (9.24), когда xt и xi+h имеют значения +1 или —1, ока
зывается такой же, как и суммирование по^модулю 2 двух последо вательностей, которые сдвинуты на АД по отношению друг к другу. Поскольку результатом сдвига и суммирования является другая последовательность, у которой число интервалов одного уровня пре вышает число интервалов другого уровня на единицу, то число ин тервалов, в которых xt и xi+k имеют одинаковый знак, на единицу меньше, чем число интервалов, имеющих различные знаки. Тогда уравнение (9.24) принимает вид:
срхж(ЙД) = |
---- (А=+0, Z, |
2Z, 3Z, ...). |
(9.25) |
Для А = О, Z, 2Z, ... |
фхх = 1, так |
как две последовательности хг |
и xi+h идентичны. Если'придать А нецелочисленные значения в пре делах < Д из А = О, Z, 2Z, .... то фжя. (АД) оказывается линейной
V |
245 |
комбинацией двух значений 1 и —1/Z. Следовательно, автокорреля ция m-последовательности описывается выражением
' ф * * ( т ) = 1 — ^ - М - ' ( 0 < |т |< Д ) ;
Ф * * ( т ) = - у , ( 0 < | t | < ( Z - 1 ) A ) . ( 9 . 2 6 )
Поскольку псевдослучайная последовательность является перио дической с периодом ZA, автокорреляционная функция также по вторяется с периодом ZA. Автокорреляционная функция псевдо случайной двоичной переменной показана на рис. 9.11*.
Рис. 9.11. Автокорреляционная функция псевдослучайной двоич ной переменной.
Основное требование в отношении длины последовательности регистра сдвига состоит в том, чтобы период Z был достаточно боль шим и последовательные пики автокорреляционной функции не препятствовали свертыванию автокорреляционной функции с им пульсной переходной функцией. Это означает, что период должен быть больше, чем время, необходимое для затухания импульсной переходной функции до нуля. Применение 20или 25-разрядных регистров сдвига для получения очень больших периодов приводит к тому, что боковые полосы —1/Z становятся крайне малыми. Однако их присутствие необходимо при обработке информации для целого числа циклов, если используется специфическое свойство псевдослучайных двоичных последовательностей — свойство по-' стоянства боковых полос. Все это требует обработки огромного количества информации. Однако по разным причинам многие ис следователи предпочитают обрабатывать не целое число циклов, а не большую часть длинного цикла для многих циклов. В некоторых
* Если амплитуды двух уровней псевдослучайной двоичной последо вательности равны + А и —А, то правую часть уравнений (9.26) следует ум ножить на А2.
246
случаях количество получаемой информации при этом оказывается настолько ограниченным, что потеря информации, связанная с рассмотрением целого числа циклов, была бы недопустимой. В дру гих случаях используются очень длинные последовательности, по тому что величина боковых полос мала и ими можно пренебречь, не вводя какой-либо значительной ошибки. Вообще можно сказать, что отказ от целого числа циклов приводит .к отклонению боковых полос от идеала, показанного на рис. 9.11, но это отклонение уменьшается при увеличении числа циклов. Если при обработке информации рассматривается только часть одного цикла, ошибка гораздо больше, чем при использовании значительного, но не цело го количества циклов.
Спектральная плотность мощности m-последовательности*. Из теоремы Винера—Хинчина известно, что спектральная плотность мощности функции представляет собой преобразование Фурье ав
токорреляционной |
функции, т. е. |
|
|
СО |
|
0 |
^ (0 3 )= ^ ф*х(т) exp (— j' c o t ) dx. |
(9.27) |
|
— о о |
|
Поскольку фя:с (т) — периодическая величина с периодом Р, воз можно комплексное представление ряда Фурье для автокорреля ционной функции:
Ф „ М = |
|
2 Апe_Jt0n-т» |
(9.28) |
||
где |
|
П= -ОО |
|
||
Рj |
2 |
|
|||
|
(9.29) |
||||
An = - J |
|
]* |
Ф*х (т) е_ Jc°™т dx |
||
и |
—Р/2 |
|
|||
/шх = 2nnf1 -= 2зт/Р. |
(9.30) |
||||
соп - |
|||||
Тогда |
|
|
|
|
|
7 |
|
°° |
(9.31) |
||
Ф « (•>)=■ $ |
|
2 |
An elanz e~ m dt. ■ |
||
- о о Я = |
- ° ° |
|
Порядок интегрирования и суммирования можно изменить, так как ряд Фурье для автокорреляционной функции сходится равномерно:
|
о о |
с о |
|
|
|
2 Лп \ е|Ип те~,штс/т = |
|
||
|
П ~ — ОО |
— оо |
|
|
00 |
о о |
|
оо |
|
= 2 |
Ал S е - М “ - ^ А = |
2 Л п 8 ( с о — О . |
( 9 . 3 2 ) |
|
П = ----ОО |
----ОО |
п = |
— о о |
|
* Этот метод в основном аналогичен методу, использованному Хэмптоном
112].
247
Подставив уравнение (9.29) в (9.32), получим |
|
||||||||
ф * ,н = 2 |
|
_1_ |
Р / 2 |
Ф.та ( т ) |
|
|
|
|
|
|
Р |
J |
е х р |
( — |
/с о т ) с /т |
б ( с о — соп ) = |
|||
п = — оо |
|
- Р / 2 |
|
|
|
|
|
|
|
Р / 2 |
|
|
|
|
|
|
|
||
— |
j |
«Px*(T) e - J e ,rfT |
2 |
|
б (со — |
соп ) = |
|||
-р / |
|
|
|
|
П ~ |
— ОО |
|
|
|
|
|
|
|
|
|
|
|
||
= |
- ^ ф ;л» |
П= |
2 |
6 (w- |
“ n), |
(9.33) |
|||
где |
|
* |
|
— оо |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Р / 2 |
|
|
|
|
|
|
Фд-Л- М |
= |
5 |
Ф .та ( т ) e - j <0T с /т |
(9.34) |
||||
|
|
|
|
—Р / 2 |
|
|
|
|
|
представляет собой преобразование Фурье автокорреляционной функции, определенное на интервале в один период. Подстановка уравнения (9.26) в уравнение (9.34) в интервале от —Р/2 до + Р /2 дает:
Ф ;л(Сй): P(Z + 1) |
Г |
sin (coP/2Z) |
Р_ f sin (coP/2) I |
(9.35) |
|
Z2 |
L |
C0P/2Z |
' - f [ - |
coP/2 J |
|
Подставив это уравнение в уравнение (9.33), получим:
ф |
ж _ J _ |
у |
| |
Р (Z + 1) |
j~sin((Pn P/2Z) |
|
|
||
|
|
|
|
|
|
con P/2Z |
|
|
|
|
_Р |
sin (®т |
|
|
|
|
|
||
|
Z |
|
^ Ш б ( с о - с о „ ) = |
|
|
||||
|
соп Р/2 |
JJ |
|
|
|
|
|||
|
f Z+ l |
[sin (nn/Z)]2 |
1 |
[sin ял |
о , |
. |
|
||
л |
2ОО |
1_ nn/Zс |
J “ |
Т |
L — |
J / |
(9-36) |
||
и = --W |
|
|
|
|
|
|
|
|
|
Задающая частота / с равна |
|
|
|
|
|
|
|||
|
/ с = |
Zfn = |
Zd)nl2n - |
nu>xZI2n,, |
(9.37) |
||||
а период определяется |
как |
|
|
|
|
|
|
|
|
|
Р -- ZA = |
Z//c = |
Vfi — 2п/щ, |
(9.38) |
|||||
следовательно, равенство (9.36) можно представить |
в следующем |
||||||||
виде: |
|
|
|
|
|
|
|
|
|
|
|
Ф ~ ( /) = ~ 8 ( Я + |
|
|
|
||||
|
у |
Z -)^ Tsin_(njt/Z)‘J2 g Гг— п/Л ' |
(9 |
39) |
|||||
|
^ |
Z2 |
L |
nn/Z |
|
|
Z J |
v |
7 |
Я — — ОО
248