Modelica中控制系统启停的Boolean量开关频率计算方法咨询

方法1：利用Modelica内置模块统计计算 #

• 用Modelica.Blocks.Logical.Edge模块捕获Boolean信号的上升沿（启动）和下降沿（停止），每次跳变会输出一个脉冲信号
• 搭配Modelica.Blocks.Math.Counter模块统计脉冲总数，由于每次完整的启停对应两次跳变（开→关或关→开），总开关次数为脉冲总数的1/2
• 仿真结束后，将总开关次数除以仿真总时长，即可得到平均开关频率

方法2：自定义算法实时计算 #

在模型中添加自定义逻辑，记录跳变时间与次数：

// 假设你的启停信号名为startStopSignal
Boolean startStopSignal;
Real lastSwitchTime(start=0);
Integer totalSwitchEvents(start=0);
Real averageFrequency;

algorithm
// 检测信号跳变（无论上升还是下降沿）
when pre(startStopSignal) != startStopSignal then
  lastSwitchTime := time;
  totalSwitchEvents := totalSwitchEvents + 1;
end when;

// 计算平均频率：总开关次数为事件数的1/2，除以总仿真时间
averageFrequency := totalSwitchEvents / (2 * time);

• 这段代码会实时统计跳变事件，最终得到的averageFrequency就是整个仿真周期内的平均开关频率
• 如果需要瞬时频率，可以在跳变时计算1/(time - lastSwitchTime)并存储

方法3：仿真后处理计算 #

如果不想修改模型，可以导出仿真的时间序列数据（如.mat或.csv格式），用Python/Matlab等工具处理：

• 读取时间数组t和Boolean信号数组s
• 遍历数组，找出所有满足s[i] != s[i-1]的索引，统计这类索引的总数（即跳变次数）
• 总开关次数为跳变次数的1/2，除以t[-1] - t[0]得到平均频率；若要瞬时频率，计算每两次相邻跳变的时间间隔的倒数

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