Самое простое решение - использовать функцию QUANTILE , как юк предложил .
Y = quantile(X,0.01);
Однако вам понадобится Панель инструментов статистики , чтобы использовать функцию QUANTILE . Решение, которое не зависит от наборов инструментов, можно найти, отметив, что QUANTILE вызывает функцию PRCTILE , которая сама вызывает встроенную функцию INTERP1Q для выполнения первичное вычисление. Для общего случая двумерной матрицы, которая не содержит значений NaN , вы можете вычислить квантили каждого столбца, используя следующий код:
P = 0.01; %# Your probability
S = sort(X); %# Sort the columns of your data X
N = size(X,1); %# The number of rows of X
Y = interp1q([0 (0.5:(N-0.5))./N 1]',S([1 1:N N],:),P); %'# Get the quantiles
Это должно дать вам те же результаты, что и вызов QUANTILE , без каких-либо наборов инструментов.