Установите случайный seed для воспроизводимости результатов.

rng('default')

Создайте матрицу экспоненциальных случайных чисел с 5 строками и 4 столбцами.

X = exprnd(1,5,4)

X = 5×4

    0.2049    2.3275    1.8476    1.9527
    0.0989    1.2783    0.0298    0.8633
    2.0637    0.6035    0.0438    0.0880
    0.0906    0.0434    0.7228    0.2329
    0.4583    0.0357    0.2228    0.0414

Вычислите средние гармонические и средние арифметические столбцов X.

harmonic = harmmean(X)

harmonic = 1×4

    0.1743    0.0928    0.0797    0.1205

arithmetic = mean(X)

arithmetic = 1×4

    0.5833    0.8577    0.5734    0.6357

Среднее арифметическое больше, чем среднее гармоническое для всех столбцов X.

Найдите среднее гармоническое всех значений в массиве.

Создайте 3 5 2 массивами X.

X = reshape(1:30,[3 5 2])

X = 
X(:,:,1) =

     1     4     7    10    13
     2     5     8    11    14
     3     6     9    12    15


X(:,:,2) =

    16    19    22    25    28
    17    20    23    26    29
    18    21    24    27    30

Найдите среднее гармоническое элементов X.

m = harmmean(X,'all')

m = 7.5094

Найдите среднее гармоническое по различным операционным измерениям и векторам размерностей для многомерного массива.

Создайте 3 5 2 массивами X.

X = reshape(1:30,[3 5 2])

X = 
X(:,:,1) =

     1     4     7    10    13
     2     5     8    11    14
     3     6     9    12    15


X(:,:,2) =

    16    19    22    25    28
    17    20    23    26    29
    18    21    24    27    30

Найдите среднее гармоническое X по измерению по умолчанию.

hmean1 = harmmean(X)

hmean1 = 
hmean1(:,:,1) =

    1.6364    4.8649    7.9162   10.9392   13.9523


hmean1(:,:,2) =

   16.9607   19.9666   22.9710   25.9743   28.9770

По умолчанию harmmean действует по первому измерению X, размер которого не равняется 1. В этом случае эта размерность является первой размерностью X. Поэтому hmean1 является 1 5 2 массивами.

Найдите среднее гармоническое X вдоль второго измерения.

hmean2 = harmmean(X,2)

hmean2 = 
hmean2(:,:,1) =

    3.1852
    5.0641
    6.5693


hmean2(:,:,2) =

   21.1595
   22.1979
   23.2329

hmean2 является 3 1 2 массивами.

Найдите среднее гармоническое X по третьему измерению.

hmean3 = harmmean(X,3)

hmean3 = 3×5

    1.8824    6.6087   10.6207   14.2857   17.7561
    3.5789    8.0000   11.8710   15.4595   18.8837
    5.1429    9.3333   13.0909   16.6154   20.0000

hmean3 является массивом 3 на 5.

Найдите среднее гармоническое каждой страницы X путем определения первых и вторых измерений с помощью входного параметра vecdim.

mpage = harmmean(X,[1 2])

mpage = 
mpage(:,:,1) =

    4.5205


mpage(:,:,2) =

   22.1645

Например, mpage(1,1,2) является средним гармоническим элементов в X(:,:,2).

Найдите среднее гармоническое элементов в каждом срезе X(i,:,:) путем определения вторых и третьих размерностей.

mrow = harmmean(X,[2 3])

mrow = 3×1

    5.5369
    8.2469
   10.2425

Например, mrow(3) является средним гармоническим элементов в X(3,:,:) и эквивалентен определению harmmean(X(3,:,:),'all').

Создайте вектор и вычислите его harmmean, исключая значения NaN.

x = 1:10;
x(3) = nan; % Replace the third element of x with a NaN value
n = harmmean(x,'omitnan')

n = 3.4674

Если вы не задаете 'omitnan', то harmmean(x) возвращает NaN.