4. 求解器选择与配置优化

在 Simulink 控制模型里，求解器不是“点一下运行”的附属选项，而是决定仿真速度、精度和稳定性的核心配置。同一个控制器模型，换一个求解器，可能从 2 秒跑完变成 20 秒，也可能从结果平滑变成剧烈振荡。本节围绕四个问题展开：定步长与变步长怎么选、刚性系统该用谁、容差和步长怎么调、性能与精度如何权衡。

一、定步长与变步长：先看系统的时间尺度是否均匀

定步长求解器以固定间隔 h 推进，每一步计算量可预测，适合实时仿真、代码生成和硬件在环测试。变步长求解器则根据误差估计动态调整步长：系统变化剧烈时自动缩小步长，平稳时放大步长，因此在宽时间尺度系统中效率更高。

判断依据很直接：如果模型包含快速开关、PWM、脉冲或刚性环节，定步长必须小到能捕捉最快动态，否则会漏掉事件；而变步长会在这些区间自动加密，在稳态区间稀疏采样。对于纯连续、时间尺度均匀的系统，定步长配合合适步长反而更稳定、更可复现。

二、刚性系统：显式求解器为什么会“卡死”

刚性系统的特征是：系统中同时存在变化极快和极慢的动态，且快模态衰减很快但并未消失。典型例子是电机电流环（时间常数毫秒级）与转速环（时间常数秒级）耦合，或化学反应中快反应与慢反应共存。

如果用显式求解器（如 ode45、ode1、ode5）处理刚性系统，为了保持数值稳定，步长会被最快模态“绑架”，被迫缩到极小，导致仿真时间急剧增加，甚至出现数值发散。此时应改用隐式或刚性专用求解器：ode15s、ode23s、ode23t、ode23tb。它们通过隐式公式和雅可比矩阵处理快模态，允许更大步长。

选择顺序可参考：一般刚性用 ode15s；中等刚性、需要较好精度用 ode23t；刚性较强且允许低精度用 ode23tb；如果模型可解析雅可比，ode15s 配合 Jacobian 设置为 on 或提供解析式，速度会明显提升。

三、容差与步长：精度和速度的两个旋钮

变步长求解器的核心参数是相对容差 RelTol 和绝对容差 AbsTol。误差估计要求每一步满足：|e_i| ≤ max(RelTol·|y_i|, AbsTol_i)。也就是说，RelTol 控制相对精度，AbsTol 控制接近零时的绝对精度。

默认 RelTol=1e-3、AbsTol=1e-6 对多数控制模型够用。若发现结果有高频抖动或能量不守恒，可先收紧到 RelTol=1e-4、AbsTol=1e-7；若仿真太慢，可放宽到 RelTol=1e-2，但要观察关键信号是否失真。注意：AbsTol 对接近零的信号影响很大，例如电流过零点，过大会导致零点附近误差累积。

定步长求解器的核心参数是步长 FixedStep。经验上，步长应小于系统最快时间常数的 1/10～1/20。例如 PWM 频率 10 kHz，周期 100 μs，步长取 1～5 μs 较合理。若模型含离散控制器，步长最好与采样周期成整数倍关系，避免采样时刻错位。

% 变步长配置示例（ode15s，刚性系统）
set_param('ctrl_model', 'Solver', 'ode15s');
set_param('ctrl_model', 'RelTol', '1e-4');
set_param('ctrl_model', 'AbsTol', '1e-7');
set_param('ctrl_model', 'MaxStep', '1e-3');   % 限制最大步长，防止漏事件
set_param('ctrl_model', 'Jacobian', 'on');    % 启用雅可比加速

% 定步长配置示例（离散控制器）
set_param('ctrl_model', 'Solver', 'FixedStepDiscrete');
set_param('ctrl_model', 'FixedStep', '1e-5'); % 10 μs

四、性能权衡实战：一次电机控制模型的调优过程

假设有一个直流电机双闭环模型：电流环采样 10 kHz，转速环采样 1 kHz，PWM 开关频率 10 kHz，仿真时长 2 s。初始配置为 ode45、RelTol=1e-3，运行后发现仿真耗时约 45 s，且电流波形在换向处有毛刺。

第一步，判断刚性：电流环时间常数约 1 ms，转速环约 100 ms，两者相差 100 倍，且电流环快速衰减，属于典型刚性。将求解器改为 ode15s，RelTol=1e-4，AbsTol=1e-7，MaxStep=1e-4。重新运行，耗时降到约 12 s，电流毛刺明显减小。

第二步，检查离散模块：模型中 PWM 和采样保持是离散模块，变步长求解器在离散时刻需要额外处理。将 MaxStep 设为 1e-4（对应 10 kHz 的 1/10），保证每个 PWM 周期至少采样 10 次。若仍偏慢，可尝试 ode23tb，在允许稍低精度时进一步提速。

第三步，对比验证：用 ode15s 和 ode23tb 分别运行，记录仿真时间、电流峰值误差、转速稳态误差。若 ode23tb 的电流峰值误差在 2% 以内且耗时减少 30%，则可接受。最终选择需结合实时性要求和精度指标。

总结本节要点：定步长适合实时和确定性场景，变步长适合宽时间尺度和刚性系统；刚性系统优先选 ode15s 等隐式求解器；容差和步长是精度与速度的调节旋钮，需结合最快动态和关键信号误差来定；最终选择应通过耗时与误差的对比实验确定，而不是凭感觉。