如何正确使用try-catch使MATLAB循环正常进入下一次迭代？

我懂try-catch的基本逻辑，但在批量处理7758个TLE文件的MATLAB循环里，碰到异常进入catch块后，下一次迭代计算的position和velocity结果会停留在首次出错的数值上。已知[position, velocity] = states(sat, time, "CoordinateFrame", "ecef");这行代码会因为部分TLE文件数据问题报错，但文件太多没法逐个排查，所以把循环核心逻辑放进了try块，希望异常后自动跳到下一次迭代，现在达不到这个效果，求解决办法。

原循环代码如下：

% Loop through each TLE file    
for i = 1:7758
    tleFile = tleFiles{i};
    tleFilePath = tleFile;
    [epochYearArray, month, day] = extractTLEData(tleFilePath);
    
    try
        % Create the satellite object using the TLE file
        sat = satellite(sc, tleFile);
        
        % Calculate the year for the specified date
        year = 2000 + epochYearArray;
        
        % Specify the time for which you want to get position and velocity
        time = datetime(year, month, day, 12, 0, 0);

        % Compute position and velocity in ECEF coordinate frame
        [position, velocity] = states(sat, time, "CoordinateFrame", "ecef");

        % Display the results
        disp("TLE File: " + tleFile);
        position_realn(i) = (norm(position') / 1000) - 6371;
        velocity_realn(i) = norm(velocity') / 1000;
        disp(time);
        disp("Position (ECEF) [km]: " + position_realn(i));
        disp("Velocity (ECEF) [km/s]: " + velocity_realn(i));
        disp("-------------------------------");

    catch exception
        fprintf('Inconsistent data in iteration %d, skipped.\n', i);
        disp("Error Message: " + exception.message);
    end
end

解决方法 #

问题根源是结果数组未初始化+局部变量残留：position_realn和velocity_realn未提前初始化，异常时没有给对应索引赋值，导致后续迭代可能读取到之前的残留值；同时局部变量sat、position等未重置，可能影响下一次计算。

1. 提前初始化结果数组（核心修复） #

循环开始前先初始化数组为NaN，异常时自动保留空标记，避免残留旧值：

% 提前初始化结果数组，预设NaN避免残留值
position_realn = NaN(1, 7758);
velocity_realn = NaN(1, 7758);

2. 重置局部变量（可选，增强稳定性） #

在每次迭代的try块开头，清空局部变量，彻底清除上一次迭代的残留：

try
    % 重置局部变量
    sat = [];
    position = [];
    velocity = [];
    
    % 后续原有逻辑...
end

3. 把extractTLEData纳入try块（潜在风险修复） #

原代码中extractTLEData在try块外，如果该函数处理异常TLE文件时也报错，会直接终止循环，建议把这部分也放进try块：

完整修复后代码 #

% 提前初始化结果数组
position_realn = NaN(1, 7758);
velocity_realn = NaN(1, 7758);

% Loop through each TLE file    
for i = 1:7758
    tleFile = tleFiles{i};
    tleFilePath = tleFile;
    
    try
        % 重置局部变量
        sat = [];
        position = [];
        velocity = [];
        
        [epochYearArray, month, day] = extractTLEData(tleFilePath);
        
        sat = satellite(sc, tleFile);
        year = 2000 + epochYearArray;
        time = datetime(year, month, day, 12, 0, 0);
        [position, velocity] = states(sat, time, "CoordinateFrame", "ecef");

        disp("TLE File: " + tleFile);
        position_realn(i) = (norm(position') / 1000) - 6371;
        velocity_realn(i) = norm(velocity') / 1000;
        disp(time);
        disp("Position (ECEF) [km]: " + position_realn(i));
        disp("Velocity (ECEF) [km/s]: " + velocity_realn(i));
        disp("-------------------------------");

    catch exception
        fprintf('Inconsistent data in iteration %d, skipped.\n', i);
        disp("Error Message: " + exception.message);
        % 显式标记异常位置结果为NaN（可选，和初始化效果一致）
        position_realn(i) = NaN;
        velocity_realn(i) = NaN;
    end
end

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