Мне кажется, вы просто хотите усреднить все значения в матрице P.
Чтобы сделать буквально то, что в вашей формуле, вы можете использовать
mu = 1/(size(P,1)*size(P,2)) *sum(sum(P));
Для любой матрицы P size (P, dim) возвращает свой размер по указанному размеру, т. Е. Ваш M или N для dim=1 или 2. Для матриц sum вернет вектор суммы значений каждого столбца матрицы, sum, примененный к вектору, вернет сумму всех его элементов.
Однако, то же самое может быть достигнуто легче:
mu = mean(P(:));
, где P(:) - P, рассматриваемый как один столбец. mean(P) снова вычислит среднее значение для каждого столбца P (таким образом, mean(mean(P)) - это еще один способ получить среднее значение для всех элементов P).
Редактировать : Если M и N не имеют размер P вдоль заданного измерения, т.е. если вы хотите рассмотреть только первые M строк и первые N столбцов, используйте P(1:M,1:N) для ссылки на соответствующую подматрицу.