巴特沃斯带阻滤波器仿真

,at]=lp2bs(b,a,f0,bw);%频率转换Р[Ht,ft]=freqs(bt,at,0:2:400);%计算带阻滤波器指定频率点的复数频率响应Рplot(ft,abs(Ht).^2);%%绘制平方幅度响应Рxlabel('f/hz');ylabel('Butterworth ︳H(jw) ︳^2');Рtitle('Butterworth 带阻滤波器(f0=200,bw=150)(付兴娜)');Рtext(68,0.9,'n=1') %作必要的标记Рtext(115,0.8,'n=2')Рtext(305,0.8,'n=4')Рtext(306,0.9,'n=6')Рgrid on;endР3.7 仿真波形Р 图 3.1 Butterworth型模拟原型低通滤波器Р图 3.2 Butterworth型模拟带阻滤波器Р3.8 结果分析Р 从图3.2可以看到,中心频率为,带宽为,满足设计要求。Р进一步分析我们可以得到:阶数越高,边沿越陡峭,越接近于理想情况(近似方波边沿),但是这就意味着设备的复杂程度加大,成本提高。从图3.1和图3.2可知,当和就已经很接近了,两者的收敛性相差无几,一阶系统最为简单,但效果极差,远远达不到设计效果。Р4 基于凯塞窗的FIR滤波器设计Р4.1 设计思路Р设欲设计的滤波器的理想频率响应为,单位脉冲响应为,与是一对傅式变换,因此有Р根据给定的求得的一般是无限长的且是非因果的。为了得到一个因果的有限长的滤波器h(n),最直接的方法是截断,或者说用一个窗口函数对进行加窗处理Рh(n)成为实际设计FIR滤波器的单位脉冲响应,其频率响应为为:Р其中N为窗口的长度。窗口函数的形状和窗口长度N决定了窗函数法设计出的FIR滤波器的性能。Р4.2 设计要求及方案Р设计一凯塞窗的低通FIR滤波器,要求如下:Р通带截止频率;Р阻带截止频率; Р阻带最小衰减;