如何计算磁滞回线（hysteresis loop）内的面积？现有代码求助

• 核心逻辑：磁滞回线的面积对应单位体积的磁滞损耗，直接对闭合的B-H曲线做数值积分即可，最常用的是梯形积分法，用numpy的trapz函数就能实现，最后取绝对值保证面积为正。

• 结合代码的具体步骤：

  1. 先确认你的H（磁场强度）和B（磁感应强度）数据是闭合回线——也就是数据的起点和终点要重合，要是不重合，就截取到闭合部分或者补全数据。
  2. 直接用以下代码片段计算：

     import numpy as np
     # 替换成你实际的H、B数据数组
     hysteresis_area = np.trapz(H, B)
     # 取绝对值得到正的面积值
     hysteresis_area = abs(hysteresis_area)
     print(f"磁滞回线面积：{hysteresis_area}")
     

     要是积分结果不符合预期，就把参数反过来写成np.trapz(B, H)，这取决于你数据的采集顺序。
  3. 要是数据是分上升沿、下降沿两段存储的，先把两段数据合并成完整的闭合数组再计算。

• 额外提醒：

  • 注意对应H和B的单位，算出的面积单位就是磁滞损耗的单位（比如J/m³）。
  • 如果曲线有明显噪声，先用平滑工具处理数据，比如scipy的savgol_filter，避免积分误差过大。

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