IIR滤波器的设计与实现(内含设计IIR滤波器的高效方法)

写在前面:初学者学习这部分内容,要直接上手写代码可能会感到比较困难,我这里推荐一种高效快速的设计IIR,FIR滤波器的方法------MATLAB工具箱:filterDesigner。打开的方法很简单,就是在命令行键入:filterDesigner,就可以打开了。打开以后,界面是这样的:

直接输入设计的参数,想要的类型就可以很快的设计出滤波器了。点击"文件",还可以生成MATLAB代码。

还有一个很值得介绍的是:MATLAB FTVTool(滤波器可视化工具)

示例:(可视化工具会打卡下面这个窗口)

b = firls(80,[0 0.11 0.19 1],[1 1 0 0],[1 100]); % Design an FIR filter
hd = dfilt.dffir(b); % Create the direct-form FIR filter.
set(hd,'Arithmetic','fixed');  % Quantize filter using 16-bit
% fvtool(b,hd);      % Compare the fixed-point filter with reference
h1 = copy(hd);       % Copy hd to h1
set(h1,'CoeffWordLength',12); % Use 12 bits for coefficients
fvtool(hd, h1);      % Compare 12-bit & 16-bit filters
  • (笔记正文)

  • 背景知识

    • 理想滤波器的特性,实际滤波器的妥协

      • 滤波:去除信号中不想要的成分或者增强所需成分

      • 理想与实际的妥协------容许有一定的误差:1.通带可能有波纹;2.阻带可能有非零增益;3.通带和阻带之间可能有过渡带。

    • FIR vs IIR

      • FIR:可实现任意幅度和相位响应、可实现线性相位但通常需要较大的阶数𝑀。

      • IIR:基于模拟滤波器设计而转换得到;通常只能实现频率选择性滤波但通常只需较小的阶数M

  • IIR数字滤波器的设计

    • 基于模拟滤波器的设计

      • 模拟滤波器的4种就经典类型:

        • 通带阻带都单调的巴特沃斯

        • 通带等波纹,阻带单调的切比雪夫I型(过渡带比巴特沃斯窄)

        • 阻带等波纹,通带单调的切比雪夫II型(过渡带比巴特沃斯窄)

        • 通带阻带都等波纹的椭圆滤波器(过渡带最窄)

      • IIR数字滤波器设计流程:

        • 预畸变:将关键的频率指标由数字域的𝜔转换为模拟域的Ω

        • 频率转换:若设计高通、带通或带阻,可由原型模拟低通

        • 滤波器𝐻LP(𝑠)通过频率转换得到

        • 低通设计:设计出原型模拟低通滤波器𝐻LP(𝑠)

        • 双线性变换:从模拟𝑠域映射到数字𝑧域

        • 验证:验证所得𝐻(𝑧)的性能

      • 模拟原型低通滤波器:以以巴特沃兹滤波器为例(省略推导过程)

        • 由4个参数作为滤波器指标(或者是通带衰减,阻带衰减rs,rp)

        • 子主题 2

      • 模拟𝑠域的频率变换:从低通到高通、带通、带阻

        • 低通------>高通

        • 低通------>带通

        • 低通------>带阻

      • 双线性变换:从模拟𝑠域到数字𝑧域

        • 目的:从模拟滤波器𝐻(𝑠)映射到数字滤波器𝐻(𝑧)

        • 使用双线性变换前,需对指定频率𝜔 进行预畸变至Ω = 2 tan( 𝜔/2)

    • IIR滤波器设计案例


相关推荐
数据与后端架构提升之路25 分钟前
从神经元到神经网络:深度学习的进化之旅
人工智能·神经网络·学习
一行134 分钟前
电脑蓝屏debug学习
学习·电脑
咔叽布吉40 分钟前
【论文阅读笔记】CamoFormer: Masked Separable Attention for Camouflaged Object Detection
论文阅读·笔记·目标检测
johnny23341 分钟前
《大模型应用开发极简入门》笔记
笔记·chatgpt
亦枫Leonlew44 分钟前
微积分复习笔记 Calculus Volume 1 - 4.7 Applied Optimization Problems
笔记·数学·微积分·1024程序员节
小肥象不是小飞象1 小时前
(六千字心得笔记)零基础C语言入门第八课——函数(上)
c语言·开发语言·笔记·1024程序员节
星LZX1 小时前
WireShark入门学习笔记
笔记·学习·wireshark
阑梦清川1 小时前
在鱼皮的模拟面试里面学习有感
学习·面试·职场和发展
qq_433099401 小时前
Isaac Gym学习笔记——概述
学习
努力变厉害的小超超3 小时前
ArkTS中的组件基础、状态管理、样式处理、class语法以及界面渲染
笔记·鸿蒙