这个问题我之前处理实验数据时也碰到过，Octave/Matlab 默认的 interp1 外插逻辑确实没法直接实现「两端各自保持端点值恒定」的需求——默认的 'extrap' 参数是跟着你选的插值方法趋势外推（比如线性插值就会延长直线），不是我们要的平端外插。不过有两种简单的思路可以搞定：

方法一：手动拆分查询点处理 #

这种方式逻辑最清晰，把查询点分成「左侧外插、中间插值、右侧外插」三部分分别处理，适合需要自定义插值方法（比如三次样条）的场景：

% 示例原始数据
x = [1, 2, 3, 4, 5];
y = [2, 4, 6, 8, 10];
xq = [-1, 0, 2.5, 6, 7]; % 包含需要外插的点

% 获取原始数据的x范围边界
x_min = min(x);
x_max = max(x);

% 初始化结果数组
yq = zeros(size(xq));

% 1. 处理左侧外插：xq小于最小x时，取第一个y值
left_mask = xq < x_min;
yq(left_mask) = y(1);

% 2. 处理中间插值：在x范围内的点，用你需要的插值方法计算
mid_mask = (xq >= x_min) & (xq <= x_max);
yq(mid_mask) = interp1(x, y, xq(mid_mask), 'linear'); % 这里可以换成'cubic'/'spline'等

% 3. 处理右侧外插：xq大于最大x时，取最后一个y值
right_mask = xq > x_max;
yq(right_mask) = y(end);

方法二：先插值再修正外插部分 #

如果觉得拆分太麻烦，也可以先让 interp1 完成中间插值和默认外插，再把外插的部分替换成端点值，代码更简洁：

% 先执行普通插值+默认外插
yq = interp1(x, y, xq, 'linear', 'extrap');

% 修正左侧外插区域
yq(xq < min(x)) = y(1);
% 修正右侧外插区域
yq(xq > max(x)) = y(end);

额外小技巧 #

如果你的插值需求刚好是「最近邻插值+两端恒定外插」，其实直接用 interp1(x, y, xq, 'nearest', 'extrap') 就可以了——因为最近邻外插的逻辑就是取距离最近的端点值，刚好符合我们要的两端恒定效果。

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