subplot(212)

plot(f,180/pi*unwrap(angle(h)));

xlabel('频率Hz');ylabel('相位');

title('相位响应')

grid on

图4-1窗函数法设计的幅度响应和单位冲激响应特性图

4。3约束最小二乘法优化设计FIR数字滤波器

  采用约束最小二乘法（CLS）来设计FIR滤波器最大的优势就在于可以使用不明确的过渡带幅值响应，只需给定截止频率，通带或阻带的边缘频率。在一些特殊的情况下，这种不用确切指定过渡带的位置函数而只要间接定义过渡带功能是十分有效的[8]。文献综述

约束最小二乘法整体看来就是在给定所要求的滤波器幅频响应Hd(ejω)的情况下，探索一个与期望的幅频响应的相近响应函数H(^)d(ejω)，从而使在整个频域上，得到的方差∫[H(^)d(ejω)- Hd(ejω)]2dω达到最小。

从这里可以看出，这种优化设计主要在要使总的误差最小，但是却不能确保设计出来的滤波器在每个局部频段上误差最小。CLS要领在于使最大允许波纹的上下限值包括在幅值响应之中。

1）约束最小二乘法原理

一般如果优化的设计要求中包含约束条件，那么把约束方法与最小二乘法结合是个好办法。约束最小二乘法即将最小二乘结合拉格朗日乘子法[9]。