巴特沃斯带通滤波器设计

基于DQ轴谐波提取器的永磁同步电机谐波抑制 PMSM 1.通过谐波提取器，直接提取DQ轴的谐波分量进行抑制，对五七次谐波电流抑制效果效果很好。 2.为了放大效果，采用主动注入谐波电压的方法，增大了电机中的谐波分量。 3.调制算法采用SVPWM，电流环处搭建了解耦补偿模块，控制效果更好。 3.纯手工搭建，可以提供参考资料。

咱们今天聊点硬核的——永磁同步电机里的谐波怎么治。玩过电机控制的都知道，五七次谐波就像蚊子一样烦人，传统陷波器削波容易把系统搞崩。这里有个邪道玩法：直接对DQ轴下手抽血，哦不，抽谐波。

先看这个谐波提取器的核心代码：

def design_harmonic_filter(order=4, f_low=295, f_high=305, fs=10e3): nyq = 0.5 * fs low = f_low / nyq high = f_high / nyq b, a = butter(order, [low, high], btype='band') return b, a # 在dq坐标系下并行运行 d_axis_filter = design_harmonic_filter() q_axis_filter = design_harmonic_filter()

这代码看着平平无奇？玄机在截止频率设置上。295Hz到305Hz正好卡住50Hz基波的5、6倍频，用四阶巴特沃斯保证相位延迟可控。实测中发现，直接这么干谐波分量还是太弱，这时候就得主动搞事情——往电压里掺沙子。

谐波注入这块有个骚操作：

// 谐波电压注入量计算 void inject_harmonic(float theta, float* vd, float* vq) { float harmonic_amp = 0.15 * (*vd); // 基波电压的15% *vd += harmonic_amp * sin(6 * theta); // 5次谐波对应6倍电角度 *vq += harmonic_amp * cos(6 * theta); }

注意这里用6倍电角度而不是5倍，因为DQ变换本身会引入频率搬移。相当于在旋转坐标系下反向补偿，实测注入量超过20%时电机就开始哼歌了，别问我怎么知道的...

说到SVPWM，老司机都懂的死角问题得处理。这里在电流环里塞了个解耦模块：

% 解耦补偿模块Simulink实现 function [v_dcomp, v_qcomp] = decoupling_coupling(i_d, i_q, Ld, Lq, omega) persistent prev_id prev_iq; if isempty(prev_id) prev_id = 0; prev_iq = 0; end diddt = (i_d - prev_id) / Ts; diqdt = (i_q - prev_iq) / Ts; v_dcomp = omega*Lq*i_q - Ld*diddt; v_qcomp = -omega*Ld*i_d - Lq*diqdt; prev_id = i_d; prev_iq = i_q; end

这个微分项的离散化处理暗藏杀机——用前向差分虽然简单，但实际调试时发现用二阶中心差分稳定性更好。曾经有个项目因为这块没处理好，电机启动直接跳闸，现场排查到凌晨三点...

整套系统纯手工搭建的最大好处是参数可调空间大。比如谐波提取器的Q值可以根据电机转速实时调整，这在预研阶段帮我们省了三个月时间。不过要小心别让各模块的延时堆叠起来，当年用DSP28335时就因为中断优先级设置不当，导致控制周期出现诡异的10us抖动。

最后说个血泪教训：实验室调试时记得给示波器电压量程留足余量。有次谐波注入过猛，电压尖峰直接干爆了采集卡的光耦隔离，烧掉的板子现在还躺在我的抽屉里当纪念品呢...（完）