Scilab 6.1.1中Xcos搭建SIRD模型时自定义函数报错求助

阿华AIGC实验室

2026-4-28

排查Scilab Xcos中SIRD模型自定义函数报错的解决方案

你好，咱们一步步来解决你在Windows 10下用Scilab 6.1.1 Xcos仿真SIRD模型遇到的问题：

首先先给你吃个定心丸：你改用积分器模块来处理SIRD的4个状态变量（S、I、R、D）是完全正确的选择。Xcos里的基础1/S模块默认是标量输入，确实不支持直接设置向量初始条件，用多个独立积分器（或者向量积分器模块）分别处理每个状态变量的积分，是仿真多变量常微分方程组的标准操作。

接下来咱们从几个核心方向排查自定义函数的报错问题：

1. 先确认SIRD模型的微分方程是否正确映射到自定义函数

SIRD模型的标准微分方程组是：

dS/dt = -β·S·I/N
dI/dt = β·S·I/N - γ·I - μ·I
dR/dt = γ·I
dD/dt = μ·I
其中N=S+I+R+D（总人口），β为感染率，γ为康复率，μ为病死率。

请检查你的自定义函数是否准确实现了上述逻辑，重点关注：

函数输入是否包含所有必要变量（S、I、R、D、β、γ、μ等）
输出是否是对应4个微分结果的4维列向量
总人口N是否在函数内正确计算（N = S + I + R + D），有没有遗漏某个状态变量

比如符合要求的Scilab自定义函数格式应该是这样的：

function dY = sird_ode(t, Y, beta, gamma, mu)
    S = Y(1);
    I = Y(2);
    R = Y(3);
    D = Y(4);
    N = S + I + R + D;
    dS = -beta * S * I / N;
    dI = beta * S * I / N - gamma * I - mu * I;
    dR = gamma * I;
    dD = mu * I;
    dY = [dS; dI; dR; dD];
endfunction

2. 检查积分器模块的初始条件配置

每个积分器需要单独设置对应状态变量的初始值，千万不要搞混：

对应S的积分器：初始值设为9999000（也就是1e7 - 1000）
对应I的积分器：初始值设为1000
对应R的积分器：初始值设为0
对应D的积分器：初始值设为0

设置方法：双击积分器模块，在Initial condition栏输入对应数值，推荐用9999000或者1e7-1000，避免语法错误。

3. 自定义函数模块的配置检查

如果用的是Xcos里的User-defined function模块，重点检查：

输入端口数量要和你连接的变量数一致（比如需要传入S、I、R、D、β、γ、μ，那就是6个输入端口）
输出端口数量要对应4个微分方程的结果（4个输出端口）
模块中指定的函数名要和你定义的函数名完全一致，别出现拼写错误

4. 常见报错的快速排查方向

如果还是报错，可以从这几个常见问题入手：

维度不匹配：检查自定义函数输出的向量维度是否和积分器输入一致，Scilab里状态变量默认用列向量，所以输出要写成[dS; dI; dR; dD]而不是行向量
参数传递错误：确认β、γ、μ这些参数是否通过Constant模块正确传递，参数值是否合理（比如μ是病死率，应该是0到1之间的正数，不能为负）
仿真配置问题：打开菜单栏Simulation->Setup，确认仿真时间范围合理（比如0到100天），求解器选择合适的类型（比如RK4就完全能满足传染病模型的仿真需求）