Dymola条件语句使用问题：房间人数模拟代码无法运行求助

云原生阿岚

2026-6-22

问题排查与修正方案

核心问题分析

变量类型错误：E用于表示离散状态（1-4），却声明为Real类型。浮点数的相等判断（如E == 1）易因精度误差失效，导致状态切换逻辑异常。
离散变量未标记：p和E均为事件触发更新的离散变量，但未添加discrete修饰符，求解器无法正确识别变量更新规则，引发求解错误。
随机数生成逻辑冗余：Modelica原生rand()函数直接返回[0,1)区间的实数，原代码中除以32767的操作多余（若原意图是将整数随机数归一化，当前写法也不符合逻辑）。

修正后的代码

model Trial
  discrete Real p(start=rand());
  discrete Integer E(start=1); // 改用Integer存储离散状态，避免浮点数相等判断误差

equation 
  when sample(1, 3600) then // 从1秒开始，每3600秒触发一次状态更新
    p = rand();
    
    if E == 1 then
      if p < 0.2 then
        E = 2;
      elseif p < 0.6 then
        E = 3;
      else
        E = 4;
      end if;
    elseif E == 2 then
      if p < 0.2 then
        E = 1;
      elseif p < 0.6 then
        E = 3;
      else
        E = 4;
      end if;
    elseif E == 3 then
      if p < 0.2 then
        E = 1;
      elseif p < 0.6 then
        E = 2;
      else
        E = 4;
      end if;
    elseif E == 4 then
      if p < 0.2 then
        E = 1;
      elseif p < 0.6 then
        E = 2;
      else
        E = 3;
      end if;
    end if;
  end when;

  annotation (
    Icon(coordinateSystem(preserveAspectRatio=false)),
    Diagram(coordinateSystem(preserveAspectRatio=false)),
    experiment(StopTime=86400, __Dymola_Algorithm="Dassl")
  );
end Trial;