Уменьшайте размерности набора точки Sobol
Сгенерируйте семимерный набор точки Sobol и скремблируйте точки.
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
Разделите первые 7 168 точек в ps на семь уровней 1 024 точек каждый. Уменьшайте первые 1 024 точки, чтобы быть одномерными, вторые 1 024 точки, чтобы быть двумерными, и так далее. Для каждого уровня вычислите отклонение значений точки в каждой размерности.
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')
