用Octave/MATLAB仿真分析CRM PFC：开关频率随功率、电压变化的完整代码与避坑指南

用Octave/MATLAB仿真分析CRM PFC：开关频率随功率、电压变化的完整代码与避坑指南

在电源设计领域，临界导通模式(CRM)下的功率因数校正(PFC)电路因其高效率、低谐波失真等优势，成为中低功率应用的热门选择。然而，CRM PFC的开关频率会随输入电压和输出功率动态变化，这一特性给电路设计和元件选型带来挑战。本文将手把手带你用Octave/MATLAB搭建完整的仿真模型，通过可视化分析揭示开关频率的变化规律，并提供可直接复用的代码模板和工程实践中的避坑技巧。

1. CRM PFC工作原理与关键公式解析

CRM PFC的核心在于电感电流在每个开关周期都从零开始上升，在达到输入电流包络后立即关断。这种工作模式自然实现了电流对电压的跟踪，但同时也带来了开关频率不固定的特性。

1.1 基本时序关系

CRM模式下每个开关周期包含两个阶段：

• 导通阶段(Ton)：开关管导通，电感电流线性上升

  ton = L * 2 * Iin_peak / Vin_peak; % 导通时间计算

• 关断阶段(Toff)：开关管关断，电感电流通过二极管续流

  toff = L * 2 * Iin_peak * sin(w*t) ./ (Vo - Vin_peak*sin(w*t));

总开关周期T为两者之和，开关频率fs=1/T。从公式可以看出，Toff与输入电压瞬时值密切相关，这直接导致了开关频率的变化。

1.2 关键参数影响分析

提示：实际设计中需要在这些参数间取得平衡，既要保证足够的频率范围，又要避免极端工况下的频率过高或过低。

2. 完整仿真模型搭建

下面我们构建一个完整的Octave/MATLAB仿真脚本，分析开关频率随时间和工况的变化。

2.1 基础参数设置

% 系统参数配置 Vo = 400; % 母线电压(V) L = 108e-6; % 电感(H) f_line = 50; % 电网频率(Hz) w = 2*pi*f_line; eta = 0.98; % 预估效率 % 仿真时间设置 t_step = 1e-6; % 时间步长(s) t_sim = 0:t_step:0.02; % 仿真2个工频周期

2.2 开关频率随时间变化分析

% 固定功率下的瞬时特性分析 P = 2000; % 额定功率(W) Vin_rms = 220; % 输入电压有效值(V) Iin_rms = P / (eta * Vin_rms); Iin_peak = Iin_rms * sqrt(2); Vin_peak = Vin_rms * sqrt(2); ton = L * 2 * Iin_peak / Vin_peak; toff = L * 2 * Iin_peak * sin(w*t_sim) ./ (Vo - Vin_peak*sin(w*t_sim)); T = ton + toff; fs = 1 ./ T; figure; yyaxis left; plot(t_sim, fs/1e3); ylabel('开关频率(kHz)'); yyaxis right; plot(t_sim, Vin_peak*sin(w*t_sim)); ylabel('输入电压(V)'); xlabel('时间(s)'); title('开关频率随输入电压瞬时值变化'); grid on;

运行这段代码，你将看到典型的"马鞍形"频率曲线——在输入电压过零处频率最高，在电压峰值处频率最低。

3. 功率与电压对开关频率的影响

3.1 功率变化的影响分析

固定输入电压为220V，观察功率从100W到3000W变化时的频率特性：

P_range = 100:100:3000; % 功率扫描范围 t_analyze = 0.005; % 选择电压峰值时刻分析 Iin_rms = P_range ./ (eta * Vin_rms); Iin_peak = Iin_rms * sqrt(2); ton = L * 2 * Iin_peak / Vin_peak; toff = L * 2 * Iin_peak * sin(w*t_analyze) / (Vo - Vin_peak*sin(w*t_analyze)); fs_power = 1 ./ (ton + toff); figure; plot(P_range, fs_power/1e3, 'LineWidth', 1.5); xlabel('输出功率(W)'); ylabel('开关频率(kHz)'); title('开关频率随输出功率变化(Vin=220V)'); grid on;

从结果曲线可以看出，随着功率增加，开关频率几乎呈双曲线下降。这对散热设计提出了挑战——高功率时频率降低可能导致磁性元件体积增大，而低功率时频率升高又可能增加开关损耗。

3.2 输入电压变化的影响分析

固定输出功率为2000W，观察输入电压从110V到290V变化时的频率特性：

Vin_range = 110:5:290; % 输入电压扫描范围 P_fixed = 2000; % 固定功率 Iin_rms = P_fixed ./ (eta * Vin_range); Iin_peak = Iin_rms * sqrt(2); Vin_peak = Vin_range * sqrt(2); ton = L * 2 * Iin_peak ./ Vin_peak; toff = L * 2 * Iin_peak * sin(w*t_analyze) ./ (Vo - Vin_peak*sin(w*t_analyze)); fs_vin = 1 ./ (ton + toff); figure; plot(Vin_range, fs_vin/1e3, 'LineWidth', 1.5); xlabel('输入电压有效值(V)'); ylabel('开关频率(kHz)'); title('开关频率随输入电压变化(P=2000W)'); grid on;

有趣的是，开关频率随输入电压呈现先升高后降低的趋势，存在一个使频率最大的最优输入电压点。这一现象对宽电压输入设计尤为重要。

4. 工程实践中的关键问题与解决方案

4.1 频率变化范围控制策略

CRM PFC设计中，开关频率可能变化5-10倍，这带来诸多挑战：

1. EMI设计难题：宽频带噪声难以滤波
   • 解决方案：设置频率钳制电路或采用变频展频技术
2. MOSFET选型困境：高频时导通损耗占优，低频时开关损耗占优
   • 建议：使用品质因数(FOM)更优的新型器件如GaN HEMT

4.2 仿真中的常见错误

以下表格总结了初学者常犯的错误及修正方法：

4.3 参数优化实用技巧

% 电感量优化示例 L_test = [50, 100, 150, 200] * 1e-6; colors = ['r', 'g', 'b', 'k']; figure; hold on; for i = 1:length(L_test) ton = L_test(i) * 2 * Iin_peak / Vin_peak; toff = L_test(i) * 2 * Iin_peak * sin(w*t_sim) ./ (Vo - Vin_peak*sin(w*t_sim)); fs = 1 ./ (ton + toff); plot(t_sim, fs/1e3, colors(i), 'LineWidth', 1.5); end legend('50μH', '100μH', '150μH', '200μH'); xlabel('时间(s)'); ylabel('开关频率(kHz)'); title('不同电感量对开关频率的影响'); grid on;

通过这段代码可以直观比较不同电感量对频率范围的影响。实际工程中，建议：

• 先确定最高可接受频率(由MOSFET性能决定)
• 再确定最低允许频率(避免电感饱和)
• 最后折中选择电感值

5. 高级话题：频率特性改善方法

对于要求严格的应用，可以考虑以下进阶技术：

1. 混合模式控制：在轻载时切换到DCM模式，限制最高频率
2. 变电感设计：使用饱和电感，在高频时自动降低电感量
3. 数字控制实现：通过MCU动态调整控制参数

% 混合模式控制示例 fs_max = 150e3; % 设置最大频率限制 fs_limited = min(fs, fs_max); figure; plot(t_sim, fs/1e3, 'b', t_sim, fs_limited/1e3, 'r--'); legend('原始频率', '限频后'); xlabel('时间(s)'); ylabel('开关频率(kHz)'); title('频率限制技术效果对比'); grid on;

这种简单的频率钳制虽然会引入少量失真，但能有效降低开关损耗，特别适合对THD要求不高的应用。