Уменьшить размеры набора точек Соболь
Создайте семизмерный набор точек Соболя и скремблируйте точки.
p = sobolset(7);
ps = scramble(p,'MatousekAffineOwen')ps =
Sobol point set in 7 dimensions (9007199254740992 points)
Properties:
Skip : 0
Leap : 0
ScrambleMethod : MatousekAffineOwen
PointOrder : standard
Разделить первые 7168 точек в ps на семь уровней по 1024 балла каждый. Уменьшите первые 1024 точки до одномерных, вторые 1024 точки до двумерных и так далее. Для каждого уровня вычислите дисперсию значений точек в каждом размере.
variance = NaN(7); for level = 1:7 pr = reduceDimensions(ps,level); pr.Skip = (level-1)*1024; pts = pr(1:1024,:); variance(level,1:level) = var(pts); end
Постройте график отклонений. Темно-синие полосы показывают дисперсию точек в первом измерении, темно-оранжевые полосы показывают дисперсию точек во втором измерении и так далее.
bar(variance) xlabel('Level') ylabel('Variances')
