死磕带通滤波器

大家好，又见面了，我是你们的朋友全栈君。

带通滤波器的作用

与陷波器类似，带通滤波器在数字电源控制领域有重要作用。比如在三相LCL逆变器的谐振抑制控制方面，通过带通滤波器可以提取谐振点附近的频谱做进一步的控制策略。在有源电力滤波器利用带通滤波器可以提取电网信号的基波频率从而做进一步的控制。

带通滤波器传递函数

带通滤波器的传递函数是：
$h(s)=\frac{Aw_oBs}{s^2+Bs+w_o^2}$
其中， $w_o$ 是带通的“中心频率”，也就是想要通过频率的中心点频率。 $B$ 是带通的频宽比，注意此处频宽比是一个相对于中心频率的比例，比如:
$w_o=50*2*pi$
$B = 0.2$
表达的意义是设定中心频率为50Hz，带通的带宽为50*0.2=10Hz。

带通滤波器的伯德图

设定“中心频率”为50Hz，频宽比为0.4。用matlab绘制伯德图，如下：
带通滤波器的幅频特性
可见，仅仅在50Hz附近有大于0的增益，其他频率点都被抑制了。于是就有了“带通”的效果。

离散化

上述都是在连续域中分析的，但是对于数字控制应用，它是无法落地实现的，所以我们需要对连续域模型进行离散化分析。

Z变换

利用Z变换可以离散化。也可以利用matlab对S函数进行Z变换，选定离散时间Ts=0.0002，则其Z变换如下：
$F(z)=\frac{0.0012557z-0.0012557}{z^2-1.996z+0.999920}$

差分方程

z变换后很自然能得到差分方程,只需要对分子分母除以 $z$ 的最高次幂：
$\frac{Y}{X}=\frac{0.0012557X_{k-1}-0.0012557X_{k-2}}{1-1.996X_{k-1}+0.999920X_{k-2}}$
有了差分方程，程序的实现可以落到实地。在Matlab的m文件中编写matlab function为例说明：

function Y = BandFilter(X)
%#codegen
%% 中间变量定义及初始化
Num0 = 0;
Num1 = 0.0012557;
Num2 = -0.0012557;

Den0 = 1;
Den1 = -1.996;
Den2 = 0.999920;

persistent Xk_1; % 上1次的输入
persistent Xk_2; % 上2次的输入

persistent Yk_1; % 上1次的输出
persistent Yk_2; % 上2次的输出

if isempty(Xk_1)   Xk_1 = 0;
end
if isempty(Xk_2)   Xk_2 = 0;
end
if isempty(Yk_1)   Yk_1 = 0;
end
if isempty(Yk_2)   Yk_2 = 0;
end
%% 执行计算
Temp = Num0*X + Num1*Xk_1 + Num2*Xk_2 - (Den1*Yk_1 + Den2*Yk_2);
Y = Temp/Den0;
Xk_2 = Xk_1;
Xk_1 = X;
Yk_2 = Yk_1;
Yk_1 = Y;