使用MATLAB从列向量中获取中值

Question

我有两个列向量。

第一列向量的长度为数千个数据点，我需要从前四十个项目中获取中值，然后从下一个四十个项目中获取中值，依此类推。

第二列向量包含组ID（从1到3）。

我的目标是最终得到一堆中位数计算，并按组对它们进行排序。 我不确定如何在MATLAB中进行此操作。

Answer 1

将向量reshape为40xN的矩阵，然后使用median取每一列的中位数。

Answer 2

这里有一些代码可以帮助您入门。

如果两个向量都在一个命名变量中，并且列数可被40整除，请执行以下操作：

% column 1 = data, column 2 = groupID
test = rand(400,2);

% compute medians of data
medians = median( reshape(test(:,1), 40,[]) );

% make each entry correspond to the correct groupID
medians = repmat(medians, 40,1);
medians = medians(:);

如果您的数据不能被40整除，请使用一个简单的循环：

N = 40;

test = rand(10*N+4,2);

n = 1;
medians = zeros( ceil(size(test,1)/N), 1 );
for ii = 1:numel(medians)
    if n+N-1 > size(test,1)
        medians(ii) = median(test(n:end,1));
    else
        medians(ii) = median(test(n:n+N-1,1));
    end
    n = n+N;
end

并根据需要像以前一样复制。

如果您将groupID放在单独的变量中，或者如何根据groupID对这些内容进行排序，则对此代码的调整非常简单。

Answer 3

建立小组相当容易：

groupIDvec = groupID(1:40:end);% A vector with group numbers

通过使用reshape ，可以按照@Oli所述找到每个组的中位数

medianmat = reshape(datavector,40,[]); 
medianvec = median(medianmat); 

现在，您只需要对它们进行排序：

[groupIdvec,idx] = sort(groupIDvec)

这是您的排序结果，其中groupIDvec指示每个值在哪个组中：

result = medianvec(idx);

我手头没有Matlab，因此可能包含错误，但是应该没问题。