Matlab interp1函数输出最后一行为NaN的问题求助

嘿，我之前处理发动机试验数据的时候也碰到过一模一样的问题！那种最后一行全是NaN的情况真的让人头大，不过搞清楚原因就好办了——核心就是曲柄转角的循环特性没被你用的interp1函数正确识别到。

为什么会出现最后一行NaN？ #

曲柄转角是循环的（0°~360°循环往复），但interp1默认是线性插值，只会在你原始数据的转角范围内插值。当你用行数最多的文件创建目标数组时，它的最后一个转角大概率接近360°，甚至可能刚好等于360°，而你的原始数据里可能没有360°之后的点（毕竟循环回到0°），这时候interp1找不到插值依据，就直接返回NaN了。

具体解决方案 #

1. 先处理数据的循环连续性 #

把每个文件的首尾数据衔接起来，相当于给循环“补个尾巴”，让插值函数能识别到360°之后其实就是0°的重复。比如把第一个点的转角加360°，对应的行数据复制到文件末尾：

% 假设你已经把所有文件数据存在all_files_data这个cell数组里，每个元素是[crank_angle, data_cols]
for idx = 1:length(all_files_data)
    ca = all_files_data{idx}(:,1);
    data = all_files_data{idx}(:,2:end);
    
    % 检查原始数据是否覆盖完整0~360°循环
    if abs(ca(end)-360) < 0.1 && abs(ca(1)) < 0.1
        % 扩展数据：把第一个点的转角+360°，数据复制到末尾
        ca_extended = [ca; ca(1)+360];
        data_extended = [data; data(1,:)];
        all_files_data{idx} = [ca_extended, data_extended];
    else
        warning('第%d个文件的曲柄转角没有覆盖完整0~360°循环，检查一下数据！', idx);
    end
end

2. 生成合理的目标转角序列 #

别直接用最大行数的原始转角，而是生成一个覆盖0~360°但避开刚好360°的序列（因为360°和0°是同一个位置）：

max_rows = max(size(file,1) for file in all_files_data);
% 生成从0°到接近360°的均匀序列，避免最后一个点刚好是360°
target_ca = linspace(0, 360-eps, max_rows);

3. 用循环插值或扩展数据插值 #

现在有两种插值方式可选：

方式一：用interp1的循环插值选项（Matlab R2017b及以上支持）

直接指定'circular'作为外插选项，它会自动把数据当成循环的：

% 初始化结果数组
result_array = NaN(max_rows, 20);
result_array(:,1) = target_ca;

for idx = 1:length(all_files_data)
    src_ca = all_files_data{idx}(:,1);
    src_data = all_files_data{idx}(:,2:end);
    % 对每一列数据插值
    for col = 2:20
        result_array(:,col) = interp1(src_ca, src_data(:,col-1), target_ca, 'linear', 'circular');
    end
end

方式二：用我们扩展后的插值数据

如果你的Matlab版本比较老，不支持'circular'，就用之前扩展的首尾数据来插值：

result_array = NaN(max_rows, 20);
result_array(:,1) = target_ca;

for idx = 1:length(all_files_data)
    src_ca = all_files_data{idx}(:,1);
    src_data = all_files_data{idx}(:,2:end);
    for col = 2:20
        result_array(:,col) = interp1(src_ca, src_data(:,col-1), target_ca, 'linear');
    end
end

这样处理后，最后一行就不会再出现NaN了——因为插值函数能找到360°附近的衔接数据，完美适配发动机循环的特性！

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