OpenModelica中使用sample()函数无法为数组赋值的问题求助

你的代码遇到的问题主要来自两个核心点：未初始化的离散变量和数组元素的状态未明确定义，下面逐一拆解并给出修复方案：

1. 警告产生的原因 #

(1) 离散变量缺失初始值 #

你声明的iTick和y_buf都是离散变量，但没有设置初始值。OpenModelica编译时无法确定这些变量的起始状态，因此会自动假设固定起始值，同时抛出“Assuming fixed start value”的警告——这可能和你预期的初始化逻辑不符。

(2) 数组元素的状态依赖冲突 #

当iTick不在1到N的范围内时，你没有定义y_buf元素的更新规则。对于未被赋值的数组元素，模型会默认认为它们的当前值等于上一时刻的值（即y_buf[x] = $PRE.y_buf[x]），同时反向依赖也成立，形成了一个结构奇异的线性系统，这就是你看到“structurally or numerically singular”警告的根源。

2. 修复后的代码示例 #

block DFT
  import Modelica.Constants.pi;
  parameter Integer N = 360 "Total number of samples";
  Integer iTick(start=0) "Counter initialized to 0";
  Real y_buf[N](start=fill(0, N)) "Array initialized with zeros";
algorithm
  when sample(0, 0.1) then
    iTick := iTick + 1;
    if iTick >= 1 and iTick <= N then
      y_buf[iTick] := iTick;
    end if;
    // 可选：防止计数器无限递增，超出数组索引范围
    if iTick >= N then
      iTick := N;
    end if;
  end when;
end DFT;

3. 关键修改说明 #

• 初始化计数器iTick：添加(start=0)，确保计数器从0开始，第一次采样后变为1，正好匹配Modelica数组从1开始的索引规则。
• 初始化数组y_buf：用(start=fill(0, N))给所有数组元素设置初始值0，让未被赋值的元素有明确的初始状态，消除符号分析时的奇异系统问题。
• 限制计数器范围：增加额外条件，当iTick达到N后不再递增，避免计数器无限变大，也防止后续采样触发无效的数组赋值操作。

4. 扩展：采样正弦函数的场景 #

如果要回到你最初的需求——采样正弦函数并填充数组，只需修改y_buf[iTick]的赋值逻辑即可，示例如下：

// 先添加频率参数
parameter Real f = 1.0 "Frequency of sine wave";
// 替换赋值语句
y_buf[iTick] := sin(2*pi*f*(iTick-1)*0.1);

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