如何在OpenModelica中实现迭代式最大功率跟踪控制器及修复代码异常？

需求说明 #

要实现一个带电气OnePort的控制器，通过迭代调节输入电压追踪最大输入功率。输入特性为：电压低于某未知阈值时电流恒定，超过阈值后电流指数下降；控制器需通过递增电压、测量功率并与前次对比的迭代方式，逐步确定阈值并锁定最大功率点。

现有代码与问题 #

用户设置Vstart(start = 30)，尝试的代码如下：

algorithm
  Vpv := Vstart;
  if P_pv_current > (P_pv_prev + dPthreshold) then
    P_pv_prev := P_pv_current;
    if increase then
      Vstart := Vstart + stepSize;
    else
      Vstart := Vstart - stepSize;
    end if;
  elseif P_pv_current < (P_pv_prev-dPthreshold) then
    P_pv_prev := P_pv_current;
    increase := not increase;
    if increase then
      Vstart := Vstart + stepSize;
    else
      Vstart := Vstart - stepSize;
    end if;
  end if;
  Vpv := Vstart;
equation
  v1 = Vpv;
  i1 = Ipv;
  
  P_pv_current = Vpv * Ipv;

模拟时仅进入一次if分支，之后Vpv、Vstart停留在30.01，P_pv_current与P_pv_prev相等，无法继续迭代追踪。

问题根源 #

1. 状态变量未正确声明：Vstart、P_pv_prev、increase这类需要跨时间步保留状态的变量，必须声明为离散状态变量，否则OpenModelica会将其视为瞬时变量，每次迭代后无法保留值。
2. 死区逻辑导致停滞：当功率变化处于[P_pv_prev - dPthreshold, P_pv_prev + dPthreshold]区间时，代码没有任何操作，初始调整后功率变化刚好落入死区，直接停止迭代。
3. 初始状态未初始化：P_pv_prev无初始值，第一次运行时默认取0，与初始功率对比触发一次调整后，后续因死区无法继续。

修正方案 #

1. 完整修正代码 #

model MPPTController
  extends Modelica.Electrical.Analog.Interfaces.OnePort;
  // 离散状态变量：跨时间步保留值
  discrete Real Vstart(start = 30, fixed = true);
  discrete Real P_pv_prev(start = 0, fixed = true);
  Boolean increase(start = true, fixed = true);
  // 参数配置
  parameter Real stepSize = 0.01;
  parameter Real dPthreshold = 0.001;
  // 中间变量
  Real P_pv_current;
  Real Vpv;
  Real Ipv;
equation
  // 从OnePort获取输入电流
  Ipv = i;
algorithm
  Vpv := Vstart;
  P_pv_current := Vpv * Ipv;
  
  // 初始化处理：模拟开始时设置初始功率基准
  if time < 1e-6 then
    P_pv_prev := P_pv_current;
  else
    if P_pv_current > P_pv_prev + dPthreshold then
      P_pv_prev := P_pv_current;
      // 功率上升，保持当前方向继续调压
      Vstart := increase ? Vstart + stepSize : Vstart - stepSize;
    elseif P_pv_current < P_pv_prev - dPthreshold then
      P_pv_prev := P_pv_current;
      // 功率下降，反转调压方向并调整
      increase := not increase;
      Vstart := increase ? Vstart + stepSize : Vstart - stepSize;
    else
      // 功率变化在阈值内，小幅续调避免停滞
      Vstart := increase ? Vstart + stepSize*0.5 : Vstart - stepSize*0.5;
    end if;
  end if;
  // 输出调压结果到OnePort
  v = Vpv;
end MPPTController;

2. 关键修正点 #

• 将Vstart、P_pv_prev声明为discrete Real，确保状态在时间步间保留。
• 添加初始功率基准初始化，避免第一次对比的异常值。
• 补充死区处理逻辑：功率变化在阈值内时，以减半步长沿当前方向微调，防止停滞在非最优功率点。
• 正确关联OnePort的输入电流i，保证功率计算的准确性。

3. 模拟设置注意事项 #

• 选择离散时间求解器（如Euler离散版本），适配控制器的离散逻辑。
• 根据系统功率量级调整stepSize和dPthreshold：步长过慢收敛慢，过快易振荡；阈值需匹配功率波动范围，避免过于敏感或迟钝。

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