|
Раздел "Математика\Statistics Toolbox"
Список функций Statistics Toolbox
\ \ Функции планирования эксперимента
|
Функция генерации выборки чисел латинского гиперкуба |
Синтаксис
X = lhsdesign(n,p)
X = lhsdesign(...,'smooth','off')
X = lhsdesign(...,'criterion','c')
X = lhsdesign(...,'iterations',k)
Описание
X = lhsdesign(n,p) служит для генерации выборки латинского гиперкуба X, содержащую n значений для каждого значения переменной p (генерируется матрица размерностью n-p). Колонка матрицы содержит n случайных чисел полученных на одном из интервалов (0,1/n), (1/n,2/n), ..., (1-1/n,1) с последующей их случайной перестановкой.
Значение параметра 'smooth' определяет:
- X = lhsdesign(...,'smooth','off') - случайную генерацию чисел из интервалов (0,1/n), (1/n,2/n), ..., (1-1/n,1);
- X = lhsdesign(...,'smooth','on') - генерацию значений выборки как центральных точек указанных интервалов: 0.5/n, 1.5/n, ..., 1-0.5/n.
Значение параметра 'smooth' по умолчанию - 'on'.
Параметр 'criterion' задает итерационную процедуру генерации выборки латинского гиперкуба X в соответствии со значением 'c'. Возможные значения 'criterion' приведены в следующей таблице
'none' |
Отключение итерационной процедуры |
'maximin' |
Максимизация минимального расстояния между значениями выборки |
'correlation' |
Минимизация значения коэффициента корреляции |
Параметр 'iterations' позволяет задать число итераций k при генерации выборки X по итерационной процедуре. По умолчанию число итераций K=5.
План по латинскому гиперкубу используется в случаях, когда необходимо получить случайную многомерную выборку равномерно распределенную по каждому измерению.
Примеры использования функции генерации псевдослучайных чисел
Генерация матрицы латинского гиперкуба с размерностью 3-5.
>> n=3
n =
3
>> p=5
p =
5
>> X = lhsdesign(n,p)
X =
0.7987 0.3282 0.3561 0.3846 0.1093
0.5759 0.7511 0.8447 0.8249 0.7206
0.2671 0.5183 0.1938 0.2658 0.6601
Использование параметра 'smooth' для генерации матрицы латинского гиперкуба с размерностью 3-5.
>> X = lhsdesign(n,p,'smooth','off')
X =
0.5000 0.5000 0.8333 0.8333 0.1667
0.1667 0.8333 0.1667 0.1667 0.8333
0.8333 0.1667 0.5000 0.5000 0.5000
>> X = lhsdesign(n,p,'smooth','on')
X =
0.7205 0.0991 0.1018 0.6811 0.6090
0.4773 0.8178 0.7929 0.1591 0.6734
0.2099 0.5184 0.4017 0.3733 0.2429
Использование критерия генерации выборки латинского гиперкуба
>> X = lhsdesign(n,p,'criterion','none')
X =
0.1615 0.5914 0.1567 0.5401 0.8463
0.8887 0.1401 0.4532 0.7389 0.4774
0.5224 0.7465 0.9303 0.1064 0.0686
>> X = lhsdesign(n,p,'criterion', 'maximin')
X =
0.5600 0.5293 0.1867 0.9291 0.6221
0.6800 0.0851 0.3555 0.3869 0.2643
0.0911 0.9107 0.7722 0.1237 0.7976
>> X = lhsdesign(n,p,'criterion', 'correlation')
X =
0.9945 0.9808 0.4275 0.8485 0.9488
0.2700 0.5441 0.1024 0.5194 0.4415
0.4710 0.1228 0.9720 0.2156 0.1003
Генерация латинского гиперкуба с заданным числом итераций
X = lhsdesign(n,p,'criterion', 'correlation','iterations',10)
X =
0.1667 0.1667 0.8333 0.1667 0.5000
0.8333 0.8333 0.1667 0.8333 0.8333
0.5000 0.5000 0.5000 0.5000 0.1667
\ \
|