如何对含分组数据的矩阵按组进行纵向排序？

看起来你是想把矩阵按连续3行作为一个分组，在每个分组内部对每一列的元素做纵向升序排序对吧？之前用splitapply(@sort,A,2)没得到想要的结果，问题出在分组索引和函数调用的逻辑上，我来给你两种可行的解决方法：

方法一：使用splitapply（直观易读） #

首先要构造正确的分组索引——因为每3行一组，300行的话就是100个分组，每个分组对应3行。我们可以用repelem生成和矩阵行数对应的分组标签，再用splitapply指定按列排序的处理函数：

% 示例矩阵（你的测试用例）
A = [ 1 7 4; 4 9 3; 8 5 7; 2 9 1; 7 4 4; 8 1 3];

% 生成分组索引：每3行对应同一个分组ID
groupIdx = repelem(1:(size(A,1)/3), 3);

% 对每个分组内的子矩阵按列（纵向）排序
B = splitapply(@(x) sort(x, 1), A, groupIdx);

运行这段代码后，B就会和你给出的预期结果完全一致。之前的splitapply(@sort,A,2)之所以不行，是因为第二个参数2是按列分组（把每一列当成一个独立分组），完全不符合你按行分组的需求，而且没有明确指定sort的排序维度。

方法二：使用reshape+sort（高效适合大矩阵） #

如果你的矩阵很大（比如300行8列），用矩阵重塑的方法会更高效，避免函数调用的开销：

% 将原矩阵重塑为三维数组：维度1=每组行数（3），维度2=列数（8），维度3=分组数（100）
A_reshaped = reshape(A, 3, size(A,2), []);

% 对每个分组的行维度（纵向）进行排序
A_sorted = sort(A_reshaped, 1);

% 再重塑回原来的二维矩阵格式
B = reshape(A_sorted, size(A,1), size(A,2));

这种方法利用Matlab的矩阵向量化操作，速度会比splitapply快不少，尤其适合大规模数据的处理。

两种方法的适用场景 #

• 如果你以后的分组逻辑可能变化（比如不是连续固定行数分组），选splitapply更灵活；
• 如果分组逻辑固定是连续N行，选reshape+sort的效率更高。

内容的提问来源于stack exchange，提问作者user9003011