基于MATLAB实现按x坐标子集替换/添加沙丘高度最大值列

我来帮你搞定这个需求，正好之前处理过类似的分组取最值的场景，咱们直接上可运行的代码，分两种方案供你选择：

第一步：读取原始数据 #

首先咱们先把txt里的数据读进来，MATLAB的load命令就能轻松搞定：

% 读取原始数据集
data = load('dune_h_data.txt');
x = data(:, 1);  % 提取x坐标列
h = data(:, 2);  % 提取沙丘高度列

第二步：获取分组标识与最大值 #

你提到的unique函数的三个输出正好派上用场，这里建议加上'stable'参数，保证唯一x的顺序和原始数据中首次出现的顺序一致（如果不需要保持顺序可以去掉）：

% 获取唯一x坐标、索引映射subs（每个原始x对应的分组序号）
[unique_x, ~, subs] = unique(x, 'stable');

% 按分组计算每个x对应的高度最大值
max_h = accumarray(subs, h, [], @max);

这里accumarray是MATLAB处理分组统计的神器，它会根据subs的分组标签，把对应位置的h值聚合并取最大值，结果max_h的每一项正好对应unique_x的最大高度。

方案1：替换原高度列为对应最大值 #

如果想直接把原始的高度列替换成该x分组的最大值，只需要把max_h按subs映射回原始数据的每一行：

% 替换原高度列
data(:, 2) = max_h(subs);

% 保存处理后的数据
save('dune_h_data_max_replaced.txt', 'data', '-ascii');

方案2：新增第三列存储最大值 #

如果不想修改原始高度列，而是新增一列存最大值，就把映射后的最大值和原数据拼接：

% 新增第三列存储对应x的最大高度
data_with_max = [data, max_h(subs)];

% 保存带最大值列的数据
save('dune_h_data_with_max.txt', 'data_with_max', '-ascii');

验证小技巧 #

你可以随便挑几个x值验证结果是否正确，比如：

% 比如取第一个唯一x，检查它的最大值是否正确
test_x = unique_x(1);
test_max = max_h(1);
disp(['x=', num2str(test_x), '对应的最大高度：', num2str(test_max)]);
disp('原始数据中该x的所有高度：');
disp(h(x == test_x));

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