FIR滤波器设计之最大误差最小准则

原理

考虑设计一个低通FIR滤波器，通带截止频率是${\omega }_p$，阻带截止频率是${\omega }_s$。即满足：

$$1-{\delta }_1\le H_r(\omega )\le 1+{\delta }_1,|\omega |\le {\omega }_p$$
$$-{\delta }_2\le H_r(\omega )\le {\delta }_2,{\omega }_s<|\omega |<\pi$$

接下来，$H_r(\omega )$经过变换可写成如下形式：
$$H_r(\omega )=Q(\omega )P(\omega )$$
其中：
$$Q(\omega )=\{\begin{array}{ll}1& \text{h(n)对称，h(n)长度是奇数}\\ cos\frac{\omega }{2}& \text{h(n)对称，h(n)长度是偶数}\\ sin\omega & \text{h(n)反对称，h(n)长度是奇数}\\ sin\frac{\omega }{2}& \text{h(n)反对称，h(n)长度是偶数}\end{array}$$
h(n)是FIR滤波器的单位脉冲响应。

$$P(\omega )=\sum _{k=0}^L\alpha (k)cos\omega k$$
L与FIR滤波器h(n)的长度和对称性有关。具体可见《数字信号处理》John G.Proakis

接下来：

定义：理想低通滤波器为$H_{dr}(\omega )$,在通带范围内幅度=1，阻带内幅度=0。

那么我们可以写出理想低通滤波器与我们设计的FIR滤波器的误差$E(\omega )$：
$$\begin{array}{rl}E(\omega )& =W(\omega )[H_{dr}(\omega )-H_r(\omega )]\\ & =W(\omega )[H_{dr}(\omega )-Q(\omega )P(\omega )]\\ & =W(\omega )Q(\omega )[\frac{H_{dr}(\omega )}{Q(\omega )}-P(\omega )]\end{array}$$
改下如下：
$$E(\omega )=\hat{W}(\omega )[{\hat{H}}_{dr}(\omega )-P(\omega )]$$
其中，$\hat{W}(\omega )$和 ${\hat{H}}_{dr}(\omega )$称为$W(\omega )$和$H_{dr}(\omega )$的修正，$W(\omega )$是加权因子。

将$P(\omega )={\sum }_{k=0}^L\alpha (k)cos\omega k$带入上式：
$$E(\omega )=$$