Уменьшите размерности набора точек Соболь
Сгенерируйте семимерный набор точек Соболь и скремблируйте точки.
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')
