sobolset

Квазислучайная точка Sobol установлена

Описание

sobolset объект набора квазислучайной точки, который производит точки из последовательности Sobol. Последовательность Sobol является основой 2 цифровых последовательности, которые заполняют пробел очень универсальным способом.

Создание

Синтаксис

p = sobolset(d)

p = sobolset(d,Name,Value)

Описание

p = sobolset(d) создает d- размерная точка установила p, который является sobolset объект с настройками свойства по умолчанию. Входной параметр d соответствует Dimensions свойство p.

пример

p = sobolset(d,Name,Value) свойства наборов p использование одного или нескольких аргументов пары "имя-значение". Заключите каждое имя свойства в кавычки. Например, sobolset(5,'Leap',2) создает пятимерный набор точки из первой точки, четвертой точки, седьмой точки, десятой точки, и так далее.

Возвращенный объект p инкапсулирует свойства Sobol квазислучайная последовательность. Набор точки конечен с длиной, определенной Skip и Leap свойства и пределами на размере точки устанавливают индексы (максимальное значение ²⁵³). Значения набора точки сгенерированы каждый раз, когда вы получаете доступ к p использование net или индексация круглой скобки. Значения не хранятся в p.

Свойства

развернуть все

`Dimensions` — Количество размерностей
положительный целочисленный скаляр в интервале [1,1111]

Это свойство доступно только для чтения.

Количество размерностей точек в наборе точки, заданном как положительный целочисленный скаляр в интервале [1,1111]. Например, каждая точка в точке установила p с p.Dimensions = 5 имеет пять значений.

Используйте d входной параметр, чтобы задать количество размерностей, когда вы создаете набор точки с помощью sobolset функция.

Используйте reduceDimensions возразите функции, чтобы сократить количество размерностей после того, как вы создадите набор точки.

`Leap` — Интервал между точками
0 (значение по умолчанию) | положительный целочисленный скаляр

Интервал между точками в последовательности, заданной как положительный целочисленный скаляр. Другими словами, Leap свойство набора точки задает число точек в последовательности, чтобы прыгнуть и не использовать для каждой взятой точки. Leap по умолчанию значение 0, который соответствует взятию каждой точки от последовательности.

Прыгание является методом, используемым, чтобы улучшить качество набора точки. Однако необходимо выбрать Leap значения с осторожностью. Многие Leap значения создают последовательности, которым не удается затронуть большой sub-hyper-rectangles единичного гиперкуба и, поэтому, не удаться быть универсальным квазислучайным набором точки. Для получения дополнительной информации см. [4].

Пример: p = sobolset(__,'Leap',50);

Пример: p.Leap = 100;

`PointOrder` — Укажите метод генерации
`'standard'` (значение по умолчанию) | `'graycode'`

Укажите метод генерации, заданный как 'standard' или 'graycode'. PointOrder свойство задает порядок, в котором производятся точки последовательности Sobol. Когда PointOrder установлен в 'standard', точки произвели, совпадают с исходной реализацией последовательности Sobol. Когда PointOrder установлен в 'graycode', последовательность сгенерирована реализацией, которая использует Код Грея индекса вместо самого индекса.

Можно использовать 'graycode' опция для более быстрой генерации последовательности, но программное обеспечение затем изменяет порядок сгенерированных точек. Для получения дополнительной информации о реализации Кода Грея см. [1].

Пример: p = sobolset(__,'PointOrder','graycode');

Пример: p.PointOrder = 'standard';

`ScrambleMethod` — Настройки то скремблирование управления
`0x0` структура (значение по умолчанию) | структура с `Type` и `Options` поля

Настройки, которые управляют скремблированием последовательности, заданной как структура с этими полями:

Type — Вектор символов, содержащий имя скремблирования
Options — Массив ячеек значений параметров для скремблирования

Используйте scramble возразите функции, чтобы установить скремблирования. Для списка допустимых типов скремблирования смотрите type входной параметр scramble. Ошибка происходит, если вы устанавливаете недопустимый тип скремблирования для данного набора точки.

ScrambleMethod свойство также принимает пустую матрицу как значение. Программное обеспечение затем очищает все скремблирование и устанавливает свойство содержать 0x0 структура.

`Skip` — Количество начальных точек в последовательности, чтобы не использовать
0 (значение по умолчанию) | положительный целочисленный скаляр

Количество начальных точек в последовательности, чтобы не использовать от набора точки, заданного как положительный целочисленный скаляр.

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

Пример: p = sobolset(__,'Skip',2e3);

Пример: p.Skip = 1e3;

`Type` — Тип последовательности
`'Sobol'` (значение по умолчанию)

Это свойство доступно только для чтения.

Тип последовательности, на котором квазислучайная точка установила p базируется, заданный как 'Sobol'.

Функции объекта

`net`	Сгенерируйте набор квазислучайной точки
`reduceDimensions`	Уменьшайте размерности набора точки Sobol
`scramble`	Скремблируйте набор квазислучайной точки

Можно также использовать следующие функции MATLAB^® с sobolset объект. Программное обеспечение обрабатывает объект набора точки как матрица многомерных точек.

`length`	Длина самого большого измерения массива
`size`	Размер массивов

Примеры

свернуть все

Создайте набор точки Sobol

Скрипт Open Live Script

Сгенерируйте 3D набор точки Sobol, пропустите первые 1 000 значений, и затем сохраните каждую 101-ю точку.

p = sobolset(3,'Skip',1e3,'Leap',1e2)

p = 
Sobol point set in 3 dimensions (89180190640991 points)

Properties:
              Skip : 1000
              Leap : 100
    ScrambleMethod : none
        PointOrder : standard

Примените случайное линейное скремблирование, объединенное со случайным цифровым сдвигом при помощи scramble.

p = scramble(p,'MatousekAffineOwen')

p = 
Sobol point set in 3 dimensions (89180190640991 points)

Properties:
              Skip : 1000
              Leap : 100
    ScrambleMethod : MatousekAffineOwen
        PointOrder : standard

Сгенерируйте первые четыре точки при помощи net.

X0 = net(p,4)

X0 = 4×3

    0.7601    0.5919    0.9529
    0.1795    0.0856    0.0491
    0.5488    0.0785    0.8483
    0.3882    0.8771    0.8755

Сгенерируйте каждую третью точку, до одиннадцатой точки, при помощи индексации круглой скобки.

X = p(1:3:11,:)

X = 4×3

    0.7601    0.5919    0.9529
    0.3882    0.8771    0.8755
    0.6905    0.4951    0.8464
    0.1955    0.5679    0.3192

Советы

Skip и Leap свойства полезны для параллельных приложений. Например, если у вас есть лицензия Parallel Computing Toolbox™, можно разделить последовательность точек через N различные рабочие при помощи функционального labindex. На каждом n th рабочий, набор Skip свойство набора точки к n – 1 и Leap свойство к N – 1. Следующий код показывает, как разделить последовательность через трех рабочих.
```
Nworkers = 3;
p = sobolset(10,'Leap',Nworkers-1);
spmd(Nworkers)
    p.Skip = labindex - 1;

    % Compute something using points 1,4,7...
    % or points 2,5,8... or points 3,6,9...
end
```

Алгоритмы

развернуть все

Генерация последовательности Sobol

Рассмотрите sobolset по умолчанию объект p это содержит d- размерные точки. Каждый p(i,:) точка в последовательности Sobol. j th координата i th точка, p(i,j), равно

${\begin{matrix} 0, & i = 1 \\ γ_{i} (1) v_{j} (1) \oplus γ_{i} (2) v_{j} (2) \oplus ..., & i > 1. \end{matrix}$

$γ_{i} (n)$ значения являются 0s или 1 с, таким образом что
$i - 1 = \sum_{n = 1} γ_{i} (n) 2^{n - 1} .$
Другими словами, $γ_{i} (n)$ значения являются двоичными цифрами целочисленного i – 1.
Значения _vj(n) называются числами направления. Они исключительно заданы для каждого координатного j. Для получения дополнительной информации об этих значениях смотрите Генерацию Чисел Направления.
⊕ оператор является поразрядным исключительным - или оператор. Для двух чисел, выраженных в двоичном файле, ⊕ оператор сравнивает цифры в каждом положении. Для данной позиции символа ⊕ оператор возвращает 1, если цифры в том положении отличаются, и возвращает 0, если цифры в том положении являются тем же самым.
- Например, $19 \oplus 24 = {(10011)}_{2} \oplus {(11000)}_{2} = {(01011)}_{2} = 11.$
- Точно так же $\frac{}{12} \oplus \frac{}{34} = {(0.1)}_{2} \oplus {(0.11)}_{2} = {(0.01)}_{2} = \frac{}{14} .$

Для получения дополнительной информации см. [3].

Генерация чисел направления

Набор чисел направления _vj(n) зависит от координатного j. Задайте числа направления в терминах значений _mj(n):

$v_{j} (n) : = \frac{m_{j} (n)}{2^{n}} .$

Для каждого j можно сгенерировать числа направления путем выбора следующего:

Примитивный полином в $ℤ_{2}$ из определенной степени _sj
$x^{s_{j}} + a_{j} (1) x^{s_{j} - 1} + a_{j} (2) x^{s_{j} - 2} + ... + a_{j} (s_{j} - 1) x + 1.$
Каждый коэффициент в полиноме или 0 или 1.
Числа направления начальной буквы _sj. Для каждого начального номера направления соответствующее значение _mj(n) должно быть или 1 или нечетное число меньше, чем 2ⁿ.

Остающиеся числа направления определяются следующим рекуррентным соотношением, которое использует коэффициенты примитивного полинома, предыдущих чисел направления и ⊕, поразрядно исключительного - или оператор.

$m_{j} (n) : = 2 a_{j} (1) m_{j} (n - 1) \oplus^{22} a_{j} (2) m_{j} (n - 2) \oplus ... \oplus 2^{s_{j} - 1} a_{j} (s_{j} - 1) m_{j} (n - s_{j} + 1) \oplus 2^{s_{j}} m_{j} (n - s_{j}) \oplus m_{j} (n - s_{j}) .$

sobolset использует те же примитивные полиномы и начальные числа направления, описанные в [3]. Эти параметры обеспечиваются для первых 1 111 размерностей.

Ссылки

[1] Bratley, P. и Б. Л. Фокс. “Алгоритм 659 Реализаций Квазислучайный Генератор Последовательности Собола”. Транзакции ACM на Mathematical Software. Издание 14, № 1, 1988, стр 88–100.

[2] Гонконг, H. S. и Ф. Дж. Хикернелл. “Алгоритм 823: Реализация Скремблированных Цифровых Последовательностей”. Транзакции ACM на Mathematical Software. Издание 29, № 2, 2003, стр 95–109.

[3] Джо, S. и Ф. И. Куо. “Отметьте относительно Алгоритма 659: Реализация Квазислучайного Генератора Последовательности Собола”. Транзакции ACM на Mathematical Software. Издание 29, № 1, 2003, стр 49–57.

[4] Kocis, L. и В. Дж. Уайтн. “Вычислительные Расследования Последовательностей Низкого Несоответствия”. Транзакции ACM на Mathematical Software. Издание 23, № 2, 1997, стр 266–294.

[5] Matousek, J. “На L2-несоответствии для Привязанных Полей”. Журнал Сложности. Издание 14, № 4, 1998, стр 527–556.

Документация

sobolset

Описание

Создание

Синтаксис

Описание

Свойства

`Dimensions` — Количество размерностей
положительный целочисленный скаляр в интервале [1,1111]

`Leap` — Интервал между точками
0 (значение по умолчанию) | положительный целочисленный скаляр

`PointOrder` — Укажите метод генерации
`'standard'` (значение по умолчанию) | `'graycode'`

`ScrambleMethod` — Настройки то скремблирование управления
`0x0` структура (значение по умолчанию) | структура с `Type` и `Options` поля

`Skip` — Количество начальных точек в последовательности, чтобы не использовать
0 (значение по умолчанию) | положительный целочисленный скаляр

`Type` — Тип последовательности
`'Sobol'` (значение по умолчанию)

Функции объекта

Примеры

Создайте набор точки Sobol

Советы

Алгоритмы

Генерация последовательности Sobol

Генерация чисел направления

Ссылки

Смотрите также

Темы

Введенный в R2008a

Документация Statistics and Machine Learning Toolbox

Поддержка

Документация

sobolset

Описание

Создание

Синтаксис

Описание

Свойства

Dimensions — Количество размерностей положительный целочисленный скаляр в интервале [1,1111]

Leap — Интервал между точками0 (значение по умолчанию) | положительный целочисленный скаляр

PointOrder — Укажите метод генерации 'standard' (значение по умолчанию) | 'graycode'

ScrambleMethod — Настройки то скремблирование управления 0x0 структура (значение по умолчанию) | структура с Type и Options поля

Skip — Количество начальных точек в последовательности, чтобы не использовать0 (значение по умолчанию) | положительный целочисленный скаляр

Type — Тип последовательности 'Sobol' (значение по умолчанию)

Функции объекта

Примеры

Создайте набор точки Sobol

Советы

Алгоритмы

Генерация последовательности Sobol

Генерация чисел направления

Ссылки

Смотрите также

Темы

Введенный в R2008a

Документация Statistics and Machine Learning Toolbox

Поддержка

`Dimensions` — Количество размерностей
положительный целочисленный скаляр в интервале [1,1111]

`Leap` — Интервал между точками
0 (значение по умолчанию) | положительный целочисленный скаляр

`PointOrder` — Укажите метод генерации
`'standard'` (значение по умолчанию) | `'graycode'`

`ScrambleMethod` — Настройки то скремблирование управления
`0x0` структура (значение по умолчанию) | структура с `Type` и `Options` поля

`Skip` — Количество начальных точек в последовательности, чтобы не использовать
0 (значение по умолчанию) | положительный целочисленный скаляр

`Type` — Тип последовательности
`'Sobol'` (значение по умолчанию)