上个月调一台单相逆变器样机空载波形看着挺干净一带整流性负载示波器上的电压波形直接“炸毛”FFT一算THD干到了11%以上。折腾了两天最后问题不在控制环而是LC滤波器参数从一开始就没算明白。这类事情在逆变器调试里太常见了开关频率乱定、电容凭手感选、电感只看饱和电流结果就是高频谐波压不下去低频段还跟控制环耦合出振荡。这篇东西把我自己整理的一套单相逆变器LC滤波器设计流程完整写出来从开关频率怎么选、调制方式怎么定到参数手算、MATLAB仿真验证、THD计算代码最后是样机阶段踩过的几个坑全部串起来讲一遍。适合正在做单相逆变器、在线式UPS、光伏离网电源或者只是想搞明白滤波器到底怎么设计的工程师参考。1. 设计前先定三件事调制方式、开关频率、直流母线1.1 调制方式直接决定LC要滤什么很多新手拿到题目就想套公式算LC我的习惯是先停下来问一个问题你到底用哪种调制方式因为LC滤波器是针对输出电压里的开关纹波做衰减的而不同调制方式产生的谐波频谱完全不同。单相全桥逆变器最常见的两种SPWM策略是双极性调制和单极性调制。双极性SPWM下输出电压在Vdc和-Vdc之间来回切换一个载波周期内输出脉冲的等效频率就是开关频率本身谐波主要集中在载波频率及其整数倍附近而且低次谐波分量相对偏大。单极性SPWM则是输出电压在0到Vdc正半周或0到-Vdc负半周之间切换相当于每个桥臂单独做高频斩波全桥输出端得到的等效脉冲频率是开关频率的2倍。换句话说同样的20kHz载波频率单极性SPWM输出电压里第一次谐波群的频率点已经跑到40kHz附近了。对LC滤波器来说要压制的谐波频率越高同样的衰减目标下需要的电感电容就越小体积和成本都更友好。更重要的是THD表现单极性SPWM的谐波能量更集中、边带更窄输出THD通常比双极性低不少。所以只要不是成本极度敏感或者非要规避半桥驱动复杂度单相全桥逆变器我基本都建议用单极性SPWM。用IGBT做桥臂的话单极性还能减小开关损耗一举两得。1.2 开关频率不是越高越好死区占比是个隐藏杀手开关频率怎么定看起来是个选择题其实是多个约束的折中。频率高滤波电感可以缩小电流环带宽也能做高但高频开关带来的损耗、EMI、驱动要求和死区问题都会冒出来。频率低LC体积变大输出阻抗也变大动态响应变差。那到底取多少合适我给出一个工程判断思路先按功率器件选频率范围再按THD目标微调。IGBT适合工作在8~20kHz区间再高开关损耗会明显恶化MOSFET可以跑到30~50kHz但超过30kHz后磁芯损耗和驱动损耗上升很快。中小功率单相逆变器几百瓦到3kW用MOSFET开关频率选16~20kHz最常用用IGBT的话10~16kHz是稳妥区间。我这篇后面统一按20kHz展开因为1kW级单相全桥用MOSFET做到20kHz非常典型。这里有个新手最容易忽略的点死区时间。死区是防止桥臂直通而强制插入的空白时间它会造成输出电压脉冲误差引入低次谐波而且死区时间占开关周期的比例越大THD恶化越明显。同样2us死区10kHz开关频率下占空比误差是2%20kHz下就变成4%输出THD明显变差。所以你会发现一个反直觉的现象开关频率从10kHz升到20kHz滤波效果变好了但如果不做死区补偿低次谐波反而增加THD可能不降反升。这就是为什么调THD不能只看开关频率。1.3 直流母线电压要留足调制余量直流母线电压Vdc的取值直接决定调制比。单极性SPWM下输出电压基波幅值等于M×Vdc其中M是调制比最大理论值是1。要让输出达到220Vrms峰值311V母线电压至少要311V以上。考虑线路压降、死区损失、负载突变时的动态压降工程上一般把Vdc取到360~400V。母线电压越高相同功率下电流越小但器件耐压和开关损耗也上去了所以也是一个折中。实测中我见过有人把Vdc定得偏低结果满载时调制比逼近1波形顶部开始削顶THD飙到8%以上。把母线提到380V同样负载THD立刻降到3%左右。所以设计LC参数之前先把母线电压和调制比范围确定下来后面算纹波电流时要用到这两个值。2. LC参数手算全流程先电容后电感再校验谐振点2.1 三条约束条件无功电流、纹波电流、谐振频率LC滤波器就两个元件设计维度其实就三条线电容不能太大否则逆变器要额外输出很大的容性无功电流电感不能太小否则开关纹波电流压不住LC谐振频率必须落在合理区间既远离基波频率又要明显低于开关频率。下面一个一个看。第一条电容约束。输出滤波电容只有基波电流流入这部分电流不会对外做功只是在线路里循环会增加逆变桥的电流应力。工程习惯是让电容基波电流不超过额定电流的10%~20%我一般取10%留点裕量。公式是I_cap 2 × π × f0 × C × Vrms其中f0是基波频率50HzVrms是输出电压有效值220V。第二条电感纹波约束。电感的主要任务是把开关纹波电流限制在一定范围内。纹波电流取得太小电感体积巨大取得太大会增加输出电容的发热和器件应力。工程上通常按额定电流峰值的15%~25%来设计纹波峰峰值。单极性SPWM全桥在调制比M0.5附近出现最大纹波工程简化公式为L ≥ Vdc / (8 × fs × ΔiL_max)这里的fs是三角载波频率ΔiL_max是允许的电感电流纹波峰峰值。这个公式按Buck变换器等效过来是单极性调制的常用估算式实际调试基本够用。第三条谐振频率约束。LC滤波器的转折频率fc 1 / (2 × π × √(L×C))它的取值有两个硬性边界要远高于基波频率工程上至少大于10倍即500Hz以上否则输出电压的基波也会被衰减同时要远低于开关频率工程上取开关频率的1/10~1/5否则对高频谐波的衰减能力不够。取fs/10左右是个稳妥起点既能给开关谐波40dB以上的衰减又能跟电流环带宽拉开距离。2.2 一个1kW/220V实例的完整手算过程拿一个典型规格走一遍额定功率1kW输出电压220Vrms/50Hz直流母线电压360V开关频率20kHz。这四个参数定下来剩下全是计算。先算额定电流I_rms P / Vrms 1000 / 220 4.55A峰值电流 I_pk 6.43A。第一步定电容。按电容基波无功电流不超过10%额定电流C_max 0.1 × I_rms / (2 × π × f0 × Vrms) ≈ 6.6uF这个值是上限不是目标值。取一个标准容值4.7uF比较稳妥对应的基波无功电流约0.325A只占额定电流的7.1%。如果你非要把C取到6.8uF也不会有大问题只是空载时电容自身会有一点温升满载时逆变桥电流稍微大一点。第二步定电感。纹波按峰值电流的20%算ΔiL_max 0.2 × 6.43 ≈ 1.29A代入纹波约束公式L ≥ 360 / (8 × 20000 × 1.29) ≈ 1.74mH取标称值1.8mH满足要求。第三步校验谐振频率fc 1 / (2 × π × √(1.8e-3 × 4.7e-6)) ≈ 1730Hz这个值大约是开关频率的0.087倍大于10倍基波频率落在合理区间。如果算出来fc明显高于2kHz说明电感偏小或者电容偏小需要调整如果fc低于800Hz则要考虑是否离基波太近。2.3 手算结果的校验表把上面的结果整理成一张表后续仿真和样机调试都以这张表为基准参数符号数值说明额定容量P1kW纯阻性满载输出电压Vrms220V/50Hz单相直流母线Vdc360V满载调制比约0.864开关频率fs20kHz单极性SPWM滤波电容C4.7uFCBB/薄膜电容滤波电感L1.8mH饱和电流≥8A谐振频率fc约1.73kHz0.087×fs这里多说一句电感选型。模算出来的1.8mH只是高频纹波约束下的电感量真正的磁芯选择还要看两个电流一是满载电流有效值4.55A二是叠加纹波后的峰值电流约为7.7A6.43A峰电流加一半纹波。电感饱和电流一定要留到8~10A以上否则峰值电流处电感量掉下来纹波会突然增大THD跟着失控。磁芯材料我习惯用铁硅铝或铁氧体开气隙前者直流叠加特性好后者高频损耗小均适用20kHz这个频率点。3. MATLAB仿真建模与LC参数脚本化设计3.1 Simulink模型怎么搭才有参考价值参数算完先别急着打样用MATLAB/Simulink把整个逆变器跑一遍。仿真不是为了好看而是为了在动手焊板子之前确认LC参数没有低级错误。模型搭建的思路是这样的直流源给360V用四个MOSFET搭全桥Simulink里直接用Universal Bridge选Single-phase Bridge功率器件用MOSFET/Diode组合后面接LC滤波器再带一个48.4Ω的纯阻性负载模拟满载。调制部分如果用Ships? 我通常用PWM Generator模块内部设置载波频率20kHz输出两路互补带死区的PWM给桥臂调制波用正弦信号源幅值设0.864频率50Hz。之所以不用闭环控制是因为验证LC参数时开环就够了闭环反而会把控制环的响应混进来干扰对滤波器本身的判断。关键设置有几个容易被忽略的点powergui选离散仿真采样步长建议1e-6s。20kHz的开关脉冲步长太大会丢失开关细节THD计算结果失真步长太小仿真速度又慢得让人想砸电脑。1e-6s在这个规模下速度和精度都比较平衡。电压电流测量结果用To Workspace模块导出Save format选Structure with time变量名比如vout和iout后面脚本处理直接用。仿真时间至少跑0.2s也就是10个基波周期。只跑两三个周期的话后面第4章算FFT的时候频谱分辨率不够频率泄漏会搞得你怀疑人生。Simulink的电气模型建议直接从带initial variables的Simscape电源开始源用DC Voltage Source就行不需要额外稳压环节。模型里我给电感串联了一个0.1Ω的等效直流电阻仿真更接近真实磁芯铜损BC? 但这个电阻会让截止频率处的峰值稍微平滑一点不影响参数验证。3.2 LC参数计算与传递函数校验脚本除了Simulink模型我习惯先跑一段纯MATLAB脚本把LC参数算出来同时画一下滤波器传递函数Bode图确认谐振点位置和衰减能力。脚本不长但能把前面手算的内容固化下来以后改规格只要改开头几个变量。%% lc_design.m 单相全桥SPWM逆变器 LC参数设计与校验 % 设计规格220Vrms / 50Hz / 1kW开关频率20kHz clear; clc; % 输入参数 Vrms 220; % 输出有效值 f0 50; % 基波频率 P 1000; % 额定功率 Vdc 360; % 直流母线电压 fs 20000; % 载波/开关频率 % 第一步额定电流 Irms P / Vrms; Ipk sqrt(2) * Irms; fprintf(额定电流有效值%.2f A\n, Irms); fprintf(额定电流峰值%.2f A\n, Ipk); % 第二步电容约束无功电流取10%额定电流 Icb_target 0.1 * Irms; C_max Icb_target / (2*pi*f0*Vrms); C 4.7e-6; % 取标准容值4.7uF Icb_actual 2*pi*f0*Vrms*C; fprintf(电容上限%.2f uF\n, C_max*1e6); fprintf(实际电容%.2f uF基波无功电流%.3f A\n, C*1e6, Icb_actual); % 第三步电感纹波约束纹波取峰值电流20% dI_max 0.2 * Ipk; L_min Vdc / (8*fs*dI_max); L 1.8e-3; % 取标称1.8mH fprintf(电感下限%.3f mH\n, L_min*1e3); fprintf(实际电感%.3f mH\n, L*1e3); % 第四步谐振频率校验 fc 1 / (2*pi*sqrt(L*C)); fprintf(LC谐振频率%.2f Hz\n, fc); fprintf(约为开关频率的 %.2f 倍\n, fc/fs); % 第五步传递函数与Bode图 Rload Vrms^2 / P; % 满载等效电阻 s tf(s); H 1 / (L*C*s^2 (L/Rload)*s 1); figure; bode(H, {2*pi*10, 2*pi*1e6}); grid on; title(LC滤波器传递函数满载);脚本跑完就能看到Bode图在约1.7kHz处有一个谐振峰过了谐振点幅频曲线按-40dB/dec衰减在20kHz处已经有40dB以上的衰减。这跟单极性SPWM的等效纹波频率约40kHz叠加以后开关谐波会被削掉一两个数量级理论上输出波形应该是干净的。注意这里的传递函数H是“逆变桥输出电压→负载端电压”的表达式没有考虑负载本身在Bode图上表现出来的阻抗变化。空载或者轻载时谐振峰Q值更高需要格外注意控制环的抗扰动能力这个留在第5章细说。4. THD到底怎么算从FFT原理到MATLAB实现4.1 FFT测THD的三个前提条件很多人在MATLAB里用一句fft()然后取几个幅值就开始报THD结果每次算出来都不一样原因通常不在算法而在数据准备。算THD之前必须想清楚三个前提。第一采样时间必须是整数个基波周期。FFT的隐含假设是信号周期延拓如果截取的数据不是整数倍周期频谱会泄漏50Hz的能量会漏到45Hz、55Hz这些频点上THD计算结果偏大且不稳定。所以仿真数据一定要取完整周期段比如5个50Hz周期就是0.1s。第二采样率要够高。要分析到50次谐波2500Hz按奈奎斯特定律采样率至少5kHz但开关纹波在40kHz附近如果要观察开关频段的残余谐波采样率至少要80kHz以上。Simulink步长1e-6s等效采样率1MHz这没问题。实测数据用示波器导出时也要注意采样率设置别用1kHz去采50Hz波形那就什么都算不出来了。第三直流分量要剔除。逆变器输出理论上没有直流成分但仿真或者实测数据中可能有微小的直流偏置这个偏置在FFT里会体现在0Hz处。如果不去掉它不影响各次谐波的计算但会影响基波幅值的标定所以稳妥起见先做去均值处理。4.2 一个通用的THD计算函数我把自己常用的THD计算函数贴出来直接在MATLAB里调用即可。它接受一个波形数组和采样率输出THD百分比和各次谐波幅值。function [thd_out, harmonic_amps] calc_thd(signal, Fs, f0) % calc_thd 计算电压或电流波形的THD % 输入 % signal - 波形数据列向量建议包含整数倍基波周期 % Fs - 采样率单位Hz % f0 - 基波频率单位Hz % 输出 % thd_out - 总谐波畸变率百分比 % harmonic_amps - 1~N_harm次谐波幅值基波在第一位 % 去直流 x signal(:) - mean(signal(:)); N length(x); % FFT Y fft(x); df Fs / N; % 频率分辨率 freq (0:N-1). * df; amp abs(Y) * 2 / N; % 单边谱幅值 amp(1) amp(1) / 2; % 直流分量修正 N_harm 50; harmonic_amps zeros(1, N_harm); for k 1:N_harm f_target k * f0; idx round(f_target / df) 1; win_low max(1, idx - 5); win_high min(N, idx 5); [harmonic_amps(k), ~] max(amp(win_low:win_high)); end fundamental harmonic_amps(1); thd_out sqrt(sum(harmonic_amps(2:end).^2)) / fundamental * 100; end这个函数里在目标频率附近做了±5Hz的搜索窗口目的是吸收仿真或实测中基波频率的微小偏移。如果是严格的50Hz整周期采样窗口里最大值就是精确的谐波幅值如果频率有点偏差也不会漏掉。实际用下来比死板地取频点要稳得多。调用方式load sim_data.mat; % 包含vout、Fs等变量 [thd_v, harm] calc_thd(vout, Fs, 50); fprintf(输出电压THD %.2f%%\n, thd_v); % 画谐波柱状图 figure; bar(1:20, harm(1:20)); xlabel(谐波次数); ylabel(幅值); title(输出电压谐波频谱);4.3 仿真结果里能看到什么用第3章的Simulink模型跑满载vout存入workspace调用calc_thd函数正常情况下THD应该在2%~3%以下。谐波柱状图里可以看到低频段几乎没有明显的3、5、7次谐波能量主要出现在开关频率边带附近也就是40kHz±若干倍基波的位置。这就是LC滤波器正常工作的信号——高频谐波被滤掉低频段干净。如果THD偏高先看频谱分布落在哪里。低频3、5次谐波大十有八九是死区或者调制策略的问题开关频率附近残留多说明LC转折频率太高或者电感量不够完全看不出规律的高频毛刺要考虑模型采样步长太大或者测量环节的问题。这种“频谱定位法”在第6章还会展开。5. 样机阶段容易踩的坑死区、谐振阻尼、测量误差5.1 死区让THD“低不下去”的机理仿真里可以不管死区但样机不行死区是真实世界对THD最直接的打击。死区时间以内桥臂上下管都不导通电感电流通过续流二极管续流输出电压被钳位到母线或地上实际输出电压脉冲宽度比调制指令窄了一段。这种误差在过零点附近尤其严重因为它造成一种类似“台阶”的畸变反映在频谱上就是3、5、7次等低次谐波增大。最直观的现象是开关频率已经提到20kHz死区才2usTHD仍然在5%以上降不下去。这时第一步先算死区占比——2us×20kHz4%每周期损失4%的占空比这个量级足以把THD从仿真里的2%拉到实测的6%。解决思路分三层。第一层是优化调制逻辑在可接受的短路风险下尽量缩短死区第二层是死区补偿根据输出电流极性对调制波做正负修正抵消脉冲误差第三层是控制层兜底用PR控制器在50Hz处提供很大的增益把死区引起的基波畸变压下去。做实验调试的时候这三层可以叠加使用但是每改一步都要重新测THD别一次全改上去否则不知道谁起了作用。5.2 LC谐振峰为什么空载稳定、满载震荡LC滤波器是一个典型的二阶系统在谐振频率1.7kHz附近有一个幅值尖峰。满载时负载电阻本身提供了一定阻尼这个尖峰被压掉不少空载或者轻载时负载电阻无穷大阻尼几乎为零谐振峰非常高。如果控制环带宽覆盖到了这个谐振频率附近或者负载突然切出LC回路就可能起振输出电压波形出现很明显的1kHz~2kHz衰减振荡。因此样机调试时我强烈建议每次做负载突变测试前先核对谐振峰。给LC并联一个RC阻尼支路是最直接的办法阻尼电阻R取特征阻抗√(L/C)这个例子算出来是约20Ω阻尼电容Cd取滤波电容的5~10倍比如10uF到47uF。RC支路串联后并接在输出端能在谐振频率处提供损耗压低Q值。代价是阻尼电阻上有一点持续的损耗满载效率会掉零点几个百分点但对绝大多数应用完全可以接受。5.3 示波器测THD不准先查这三处仿真和样机的THD对不上不一定是设计问题也很可能是测量方法有问题。我踩过几次坑总结下来主要是三处。一是探头地线问题。普通无源探头的地线夹子很长探头地线本身引入的电感在高频下会和探头输入电容谐振测量结果里出现莫名其妙的高频成分。测输出电压谐波最好用差分探头或者至少把地线夹剪短直接接在输出端子附近的地上。二是示波器FFT的窗函数和采样周期。示波器本身自带FFT功能但它的窗函数通常默认是Hanning窗、记录长度和采样率很难刚好覆盖整数倍工频周期测出来的THD往往偏大。建议把波形数据以CSV格式导出来用第4章的函数离线计算这样至少能控制变量。三是注意电压探头和电流探头的相位偏差。如果用户要算功率因数或者同时分析电压电流谐波探头的相位延迟会导致计算误差。只分析电压THD的话这个问题不大但架不住调试验收时有人会对不准。5.4 常见问题快速排查表现象可能原因优先检查方向低频3、5次谐波偏高死区效应死区补偿、PR控制开关频率边带残留多LC转折频率偏高增大电感或电容输出波形有低频振荡LC谐振峰未阻尼加RC阻尼支路满载THD变化大电感饱和换更大饱和电流磁芯测量重复性差数据不满足整周期采样按第4章重新取数空载波形毛刺谐振Q值过高轻载阻尼或控制环动态6. 调优思路从参数到控制的组合拳6.1 通过频谱判断该动哪里LC参数定型之后如果THD仍然超目标不要盲目加大电感或者电容。我的习惯是先看频谱结构判断谐波来自哪个环节然后只动对应环节。低频段3、5、7次等占主导这是调制和死区问题主导LC滤波器很难压住低频谐波因为转折频率在1.7kHz。优先做死区补偿或者把反馈从开环改成闭环引入PR控制。高频段开关频率附近占主导这是LC滤波器的有效工作区间。如果残留太多说明转折频率偏高可以适当加大电感或者电容让转折频率往低走一点。但要注意别把转折频率降得太低否则动态响应变差控制环带宽也受限制。全频段都高先怀疑测量再怀疑控制环参数是否调试到位最后才怀疑滤波元件本身。样机阶段这几种情况我都遇到过“先频谱、后元件”的思路能帮你节省大量时间。6.2 闭环控制下LC谐振频率和电流环带宽的关系如果系统要做双闭环电压外环电流内环LC谐振频率和电流环带宽必须拉开距离否则电流内环会把LC谐振峰“激励”起来系统跟它形成正反馈振荡。工程经验是电流环穿越频率压到LC谐振频率的1/5~1/3以下。比如这个例子fc≈1.7kHz电流环带宽设计在300~500Hz比较稳妥再往上就要仔细做相位裕度校核。电压外环如果用PR控制器把谐振频率设在50Hz处基波增益可以做得很大但注意PR控制器的带宽要足够窄免得把50Hz附近的噪声放大。另一个常用手段是重复控制它对周期性畸变的抑制非常有效适合死区谐波和电网谐波这类固定次数的畸变。但重复控制动态响应慢通常和PI/PR并联使用。6.3 效率与体积的取舍什么时候可以考虑LCLLC滤波器结构简单、控制直观但在功率进一步增大或者THD要求极为苛刻时电感体积会变得很大。单相逆变器同样可以用LCL结构在电网侧再加一个电感构成三阶滤波高频衰减能力比LC强很多。代价是谐振峰变成两个参数设计和阻尼方案更复杂控制环的设计也更麻烦。我的判断边界是输出功率超过3kW或者对THD和EMI有严格标准时才有必要考虑LCL1kW这个量级LC加合理的阻尼已经完全够用。最后说一个我自己调试时养成的习惯每次拿到一台逆变器样机先不上满载直接用示波器看空载波形和FFT确认LC谐振点位置和理论计算对得上然后上半载看死区导致的低次谐波是否明显满载以后再对比一次。这三步走下来基本能在半小时内判断出THD问题的根源在滤波器还是控制环避免把时间浪费在来回拆板子上。这套方法经过几个项目验证对单相逆变器、UPS后级和光伏离网电源都适用希望对正在调机的你也有帮助。