Modelica含截断正态噪声代数运算的仿真收敛问题求助

核心问题原因 #

• 事件迭代收敛失效：55秒时BioPowerOutput刚好接近阈值489100，Noise.y的随机波动会让BioPowerOutput < 489100的条件反复触发切换，仿真器的事件迭代过程无法找到稳定状态，最终报错。
• 变量定义错误：Nmax是固定常量，却被定义为Real（状态变量），会引入不必要的计算开销，也可能干扰收敛逻辑。
• 噪声项的指数放大：exp(Noise.y*time)会随时间指数增长，且噪声的采样突变会让BioPowerOutput在阈值附近产生剧烈波动，进一步破坏迭代收敛性。

修复步骤 #

1. 修正变量类型：将Nmax改为parameter，明确其常量属性：

   parameter Real Nmax=25;
   

2. 添加滞回（Hysteresis）逻辑：给阈值判断加入滞回区间，避免因微小噪声波动导致条件反复切换。比如设置当BioPowerOutput低于489100-Δ时进入上升分支，高于489100+Δ时进入饱和分支（Δ可根据需求调整，比如设为1000）。
3. 调整噪声项的使用方式：避免直接将噪声与时间相乘后做指数运算，改为对饱和值添加小幅度噪声扰动，而非指数放大，比如：

   BioPowerOutput = 489100*(1 + Noise.y);
   

   这样既保留噪声的随机性，又不会让输出无限制增长，同时降低对事件迭代的干扰。

修改后的完整代码 #

model BioPowerSetup
  Modelica.Blocks.Interfaces.RealOutput BioPowerOutput;
  parameter Real alpha=0.05;
  parameter Real mu1=7/25;
  parameter Real Nmax=25; // 修正为parameter类型
  parameter Real D=0.2;
  parameter Real threshold = 489100;
  parameter Real hysteresis = 1000; // 滞回区间宽度，可按需调整
  Modelica.Blocks.Noise.TruncatedNormalNoise Noise(
    y_min=-0.028,
    y_max=0.028,
    samplePeriod=1);
equation 
  if BioPowerOutput < threshold - hysteresis then
    BioPowerOutput = alpha*mu1*(Nmax*exp(mu1*time));
  elif BioPowerOutput > threshold + hysteresis then
    BioPowerOutput = threshold*(1 + Noise.y); // 改为比例噪声，避免指数放大
  else
    // 滞回区间内保持当前值，避免条件反复切换
    der(BioPowerOutput) = 0;
  end if;
end BioPowerSetup;

额外说明 #

• 滞回区间的宽度hysteresis可以根据仿真需求调整，只要能覆盖噪声的波动范围即可，确保条件不会频繁切换。
• 如果必须保留指数形式的噪声影响，可以将噪声项改为与时间无关的小幅度扰动，或者限制指数增长的速率，避免输出在阈值附近产生剧烈跳变。

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