Matlab中2D采样矩阵的1D无循环插值最优方案问询

问题背景 #

现有维度为120×N的温度矩阵T和高度矩阵H，需对N组独立序列（每列为一组），将温度插值到维度100×1的新高度向量Hnew。原循环实现因N约200万导致效率极低：

Tnew = zeros(length(Hnew),N);
for i=1:N
   Tnew(:,i) = interp1(H(:,i),T(:,i),Hnew);
end

无循环最优实现方案 #

方案1：直接利用interp1的向量化特性 #

Matlab的interp1原生支持多维输入的批量处理，当输入矩阵H和T为同维度（120×N），且Hnew为向量时，函数会自动对每一列独立执行插值操作，直接输出100×N的结果矩阵：

Tnew = interp1(H, T, Hnew);

该方法完全无需循环，是最简洁高效的实现方式，性能远优于显式循环。

方案2：使用griddedInterpolant批量求值 #

若需重复执行插值操作（多次更换Hnew），可先创建插值对象，再批量计算结果，进一步提升复用效率：

% 创建插值对象，默认采用线性插值（与interp1默认一致）
F = griddedInterpolant(H, T);
% 对Hnew执行插值，得到100×N的结果矩阵
Tnew = F(Hnew);

单次插值时，该方案与方案1性能接近，但多次插值时可避免重复初始化插值逻辑，效率更优。

注意事项 #

• 确保H的每一列元素为单调序列（interp1默认要求输入点单调，若存在非单调情况需添加'linear','extrap'或指定其他插值方法并处理边界）。
• 内存方面，120×200万的双精度矩阵约占192MB内存，需确保Matlab有足够内存空间分配。

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