Как вычислить среднее значение для конкретного столбца матрицы, за исключением заданного значения IN MATLAB

Question

У меня есть матрица 1000x7, в столбцах которой указаны значения, равные -99.

Я хочу, чтобы вычислить среднее значение каждого столбца в отдельности, потому что я создаю таблицу с помощью fprintf, но я хочу среднее быть вычислен без учета -99 значений.

Например, принимая матрицу A, я попытался

mean(A(A(:,1) ~= -99)) 

и она работает для вычисления первого столбца средней.

Однако, когда я пытаюсь

mean(A(A(:,2) ~= -99)) 

для второго столбца, результат точно равен результата первого столбца.

Что в этом плохого? Большое спасибо за ваши советы!

Answer 1

Вы принимаете значения индекса только в первом столбце.

Попробуйте

mean(A(A(:,2) ~= -99),2)% для 2-й колонке

Answer 2

Что вам нужно, для скажем, первый столбец:

mean(A(find(A(:,1) ~= -99),1)) 

Например:

>> A 

A = 

    1  2 
    1  3 
    -99  4 
    1  5 

Затем

>> find(A(:,1)~=-99) 

ans = 

    1 
    2 
    4 

И

>> mean(A(find(A(:,1)~=-99),1)) 

ans = 

    1 

Answer 3

Вы используете вектор, A(:,1)~=-99, к индексу в массиве таким образом, это только доступа к значениям в первом столбце. Вот пример того, как вы можете сделать это вычисление для каждого столбца в массиве:

A = [1 2 3; -99, 4, 3; 10 4 8]; 

for col = 1:3 
    index = A(:,col) ~= -99; 
    mean(A(index,col)) 
end 

Answer 4

Короче говоря: Вы должны указать столбец при чтении значения:

mean(A(A(:,2) ~= -99, 2)) 

Более длинная версия: A(:,1) ~= -99 возвращает логический массив 1000 x 1, также A(:,2) ~= -99. При индексировании с помощью логических массивов Matlab эффективно выравнивает оба массива и извлекает значения, для которых логический массив равен true.

Если логический массив B имеет M элементов и массив, который индексируется в, A, имеет N элементов, где M <= N Matlab будет работать только на первых M элементах B т.е.C = A(B) будет эквивалентно *

C = A(1:M); 
C = C(B); 

В вашем примере вы индексирование в 1000 x 7 массив с 1000 x 1 логического массива, который делает все вызовы только рассмотреть первые 1000 элементов, то есть первый столбец. Поэтому вам нужно указать, в каком столбце (столбцах) вы хотите применить логический индекс.

* Если B является 1xM массивом результирующего массив будет также изменен в 1xM массив, для всех остальных случаев Mx1 массива возвращается.