MATLAB报错‘sparse must be same length’及旧版本兼容疑问

问题背景 #

运行基于最小曲率网格化算法的旧MATLAB脚本时，uniquept.m中的代码行a = sparse(ind,1,f(:,jj));抛出错误：sparse must be same length。当前使用MATLAB 2022a，其中唯一X/Y位置的数量（对应ind的长度）少于Z值的数量（对应f(:,jj)的长度），但原脚本曾可正常运行。

旧版MATLAB sparse的行为说明 #

是的，旧版本MATLAB（R2019b及更早）的sparse函数允许值向量长度大于索引向量长度——此时函数会自动截断值向量，仅保留与索引向量长度匹配的前N个元素。但从R2020a开始，MATLAB对sparse的输入校验规则进行了收紧，要求行索引、列索引和值向量三者的长度必须完全一致，否则触发上述错误。

修复方案 #

根据原脚本的逻辑需求，可选择以下修复方式：

方式1：截断值向量匹配索引长度 #

如果原逻辑仅需要每个唯一位置对应的第一个Z值，直接截断f(:,jj)以匹配ind的长度：

a = sparse(ind, 1, f(:,jj)(1:length(ind)));

方式2：聚合重复位置的Z值 #

如果原逻辑需要对同一X/Y位置的多个Z值做聚合（比如取均值、最大值），先分组计算再构建稀疏矩阵：

% 获取唯一索引并计算对应Z值的均值（可替换为mean以外的聚合函数）
unique_indices = unique(ind);
aggregated_z = arrayfun(@(idx) mean(f(ind == idx, jj)), unique_indices);
a = sparse(unique_indices, 1, aggregated_z);

方式3：修正索引生成逻辑 #

检查ind的生成代码，确保其长度与f(:,jj)一致。比如原脚本可能使用了旧版unique的输出，新版MATLAB中unique的参数或返回值格式有变化，需调整为：

% 确保ind为每个原始点对应的唯一位置索引，长度与R一致
[~, ind, ~] = unique(R, 'rows', 'stable');

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