Неверные значения std и ожидаемые значения

Question

У меня есть матрица X с 3 колонками. Для пористого вопроса X=randn(5,3). Я хочу нормализовать столбцы X S.T. каждая колонка будет иметь значение 0 и 1 std. Я использую следующий код:

X=(X-mean(X))./std(X); 

Я получаю зЬй из 1. Мое значение, однако, очень небольшое значение, близкое к 0, но не обязательно 0. Я пытался играть немного с числами, чтобы найти объяснение:

X=1:15; 
X=reshape(X,[5 3]); 
mean(X-mean(X)); 

Который дает мне 0 значение для каждого столбца.

X=1:15; 
X=reshape(X,[5 3]); 
mean((X-mean(X))./std(X)); 

Какой нет. Но 0/что-то еще 0. Что мне не хватает?

• Почему я не получаю 0 значений?
• Являются ли значения, которые я получаю достаточно для нормализации алгоритма предварительного кластеризации?

Answer 1

Вот версия, которая делает то, что, как я думаю, вы пытаетесь сделать ... вам нужно реплицировать матрицу, потому что X-mean (X) недействительно (если вы используете стандартную реализацию) - вы не можете вычесть 1x3 из 5x3.

r = 5; c = 3; 
X=randn(r,c); 
Xm=repmat(mean(X),r,1); 
Xstd = repmat(std(X),r,1); 
Xn = (X-Xm)./Xstd; 
mean(Xn) 
std(Xn) 

Для меня это печатает

ans = 

    1.0e-16 * 

-0.6661   0 0.4441 


ans = 

    1.0000 1.0000 1.0000 

который, кажется, как именно то, что вы ищете ... обратите внимание на 1e-16 мультипликатора на средних значениях ... это, по существу, 0, с некоторой ошибкой с плавающей запятой.