MMC最近电平逼近调制(NLC)在Simulink中的实现与优化

1. 单相MMC模型与最近电平逼近调制概述

在电力电子仿真领域，模块化多电平换流器(MMC)因其模块化结构、高输出电压质量和低开关损耗等优势，已成为高压直流输电(HVDC)系统的核心器件。Simulink自带的MMC案例默认采用载波移相调制(CPS-SPWM)策略，但在实际工程中，最近电平逼近调制(NLC)因其更低的开关频率和更高的效率，在大容量场合更具优势。

最近电平逼近调制的核心思想是根据当前时刻的调制波瞬时值，选择最接近的可用输出电平。与CPS-SPWM相比，NLC具有以下特点：

• 开关频率显著降低，减少了开关损耗
• 输出电压的阶梯特性更为明显
• 对电容电压平衡控制的要求更高

我在改造Simulink自带模型时发现，原系统的子模块电容电压平衡逻辑与NLC完全不兼容。这是因为CPS-SPWM采用循环投切策略，而NLC需要根据电压偏差动态调整投切顺序。这迫使我对模型进行了深度重构，包括调制波生成、电平计算和电容平衡三个核心环节。

2. 调制波生成与电平数计算实现

2.1 原模型调制部分的改造

Simulink自带的MMC案例采用标准的正弦脉宽调制(SPWM)架构，包含：

• 载波生成模块
• 比较器模块
• 驱动信号分配逻辑

这些模块在NLC策略下完全失效，必须彻底移除。我新建了以下功能单元：

1. 调制波生成器：产生标准正弦参考信号
2. 电平计算函数：实时确定所需输出电平数
3. 投切决策模块：根据电平数生成子模块投切命令

2.2 电平数计算的MATLAB函数实现

电平计算是NLC的核心算法，其数学表达式为：
[ n = round\left(\frac{u_{ref}}{V_{step}}\right) ]
其中，( V_{step} = \frac{V_{cap_sum}}{N} )，N为桥臂子模块总数。

在Simulink中实现的MATLAB Function代码如下：

matlab复制function n = NLC_Calc(u_ref, V_cap_sum)
    % 输入参数：
    % u_ref - 当前时刻的调制波电压值
    % V_cap_sum - 桥臂所有子模块电容电压总和
    % 输出参数：
    % n - 需要投入的子模块数量（可正可负）
    
    N = 40; % 单桥臂子模块总数
    V_step = V_cap_sum / N; % 计算单模块电压步长
    n = round(u_ref / V_step); % 四舍五入取整
    n = max(min(n, N/2), -N/2); % 限制在±N/2范围内
end

注意：Simulink的MATLAB Function模块不支持persistent变量存储状态，因此无法在函数内部实现变化率限制。必须在函数外接Rate Limiter模块，将电平变化率限制在每周期最多±2个电平，避免调制波突变时产生不连续跳变。

2.3 实际调试中的波形优化

初始实现时观察到的波形问题：

1. 电平跳变瞬间出现电压尖峰
2. 高次谐波含量异常增加
3. 电容电压波动明显

通过以下措施显著改善波形质量：

• 在电平计算模块后添加Rate Limiter，限制电平变化率
• 调整仿真步长至20μs以下
• 采用ode23t变步长算法

优化前后THD对比：

3. 子模块电容电压平衡策略重构

3.1 原模型排序算法的局限性

Simulink自带模型采用简单的循环投切策略，其特点是：

• 固定顺序轮流投切子模块
• 不考虑单个电容的电压偏差
• 仅适用于载波移相调制

这种策略在NLC下会导致：

• 部分电容长期处于过压状态
• 电压波动幅度超过±10%
• 系统动态响应变差

3.2 改进的带权重排序算法

我设计的新型排序算法包含以下关键改进：

1. 电压偏差加权计算：
   [ w_i = |v_{c_i} - V_{target}| \times (1 + k \times direction) ]
   其中direction反映电流方向（充电为+1，放电为-1）

2. 动态排序机制：

   • 按权重降序排列
   • 优先投入电压偏差大的电容
   • 根据电流方向调整排序权重

具体实现代码如下：

matlab复制function [gate] = CapacitorBalance(capVoltages, direction)
    % 输入参数：
    % capVoltages - 各子模块电容电压数组
    % direction - 电流方向（1充电，-1放电）
    % 输出参数：
    % gate - 子模块投切命令（1投入，0切除）
    
    targetVoltage = 5000; % 目标电压值(V)
    k = 0.1; % 方向权重系数
    
    % 计算带权重的电压偏差
    weights = abs(capVoltages - targetVoltage) .* (1 + k*direction);
    
    % 按权重降序排序
    [~, index] = sort(weights, 'descend');
    
    % 生成投切命令（前20个投入）
    gate = zeros(size(capVoltages));
    gate(index(1:20)) = 1; 
end

3.3 平衡算法参数优化

通过大量仿真实验确定的优化参数：

算法执行周期对系统性能的影响：

• 周期缩短→电压波动减小→CPU负载增加
• 需要根据仿真精度要求折中选择
• 实际工程中通常采用FPGA实现硬件加速

4. 系统级调试与性能优化

4.1 仿真步长选择策略

不同仿真算法的对比测试结果：

选择建议：

1. 初期调试使用固定步长50μs快速验证
2. 波形分析时切换至ode23t变步长
3. 需要与硬件对接时采用20μs固定步长

4.2 谐波特性分析与抑制

NLC调制固有的谐波分布特点：

• 主要谐波集中在nN±1次（N为电平数）
• 本模型中40电平对应的主要谐波为39/41次
• 谐波幅值与调制比密切相关

实测谐波频谱数据：

抑制措施：

1. 增加输出滤波器（L=5mH，C=50μF）
2. 优化调制波生成算法（添加预畸变）
3. 采用混合调制策略（NLC+PWM）

4.3 实际工程应用建议

基于仿真经验的工程实施要点：

1. 硬件设计：

   • 子模块电容容值不小于10mF
   • IGBT开关频率按1kHz设计余量
   • 电压采样精度需达到0.1%

2. 控制参数：

   • 排序算法周期≤50μs
   • 电平变化率限制在2/ms
   • 电容电压均衡系数k=0.05~0.2

3. 系统保护：

   • 设置电容电压偏差报警阈值(±10%)
   • 检测并抑制电平跳变过大情况
   • 过流保护响应时间<100μs

5. 模型扩展与进阶应用

5.1 三相MMC系统扩展

单相模型可扩展至三相系统的要点：

1. 相间电容电压协调控制
2. 环流抑制策略实现
3. 三相调制波同步生成

关键修改部分：

• 增加相间平衡控制器
• 实现120°相位差的三相调制波
• 共享排序算法资源

5.2 混合调制策略实现

结合NLC与PWM优势的混合方案：

1. 高调制比区域使用纯NLC
2. 低调制比区域切换至PWM
3. 过渡区域采用NLC+PWM

实现方法：

matlab复制function [n_final] = HybridModulation(u_ref, V_cap_sum)
    m = abs(u_ref)/(V_cap_sum/2); % 计算调制比
    
    if m > 0.8
        n_final = NLC_Calc(u_ref, V_cap_sum); % 纯NLC
    elseif m < 0.3
        n_final = PWM_Calc(u_ref); % 纯PWM
    else
        n_NLC = NLC_Calc(u_ref, V_cap_sum);
        n_PWM = PWM_Calc(u_ref);
        n_final = round(0.5*(n_NLC+n_PWM)); % 混合输出
    end
end

5.3 硬件在环测试准备

将Simulink模型部署至实时仿真器的要点：

1. 模型离散化设置：

   • 固定步长20μs
   • 使用RT模块替换连续模块
   • 禁用变步长算法

2. 代码生成优化：

   • 启用SIMD指令集
   • 设置单精度浮点运算
   • 优化内存访问模式

3. 实时性能监测：

   • 跟踪最坏情况执行时间(WCET)
   • 监控CPU负载率
   • 分析任务调度延迟