为了学习陷波器,几乎把全网的陷波器有关文章都看了一遍,整理下,加深印象,也给后来的朋友作参考,文章内容大部分直接复制过来,但是所有内容代码我都测试了一遍,其中有些错误之处,也作了修改。

陷波滤波器形式有很多种,一般简单的形式传递函数为:

这种形式滤波器陷波深度不可调

陷波滤波器差分方程推导(双线性变换法)

以100hz陷波频率为例,matlab代码如下:

f = 100;          % Hz
w = 2*pi*f;       % rad/s
zeta = 0.707;
num = [1 0  w*w];
den = [1 2*zeta*w w*w];
G = tf(num, den);
bode(G);
grid on;

从bode图中可以看出
a. 信号在100Hz处急速衰减
b. 调整阻尼比ζ,可以调节100Hz处的衰减特性

离散形式代码,很简单,我就不传了。

二参数陷波滤波器:

这种形式滤波器可以直接调节陷波带宽,缺点:陷波深度不可调

可参考这篇文章,推导很详细,代码可用:

https://blog.csdn.net/u013581448/article/details/116743786

三参数陷波滤波器:

这种滤波器陷波深度与带宽可调,传递函数:

wn:陷波频率;
k1:陷波因子1;
k2:陷波因子2;

采用双线性变换,其本质是一种数值积分法,采用梯形方法来近似计算积分。经过简要推导可以得到:

将式(2)带入式(1)可得,

将式(3)利用字符进行替换可得,

其中,

利用z变换的信号滞后特点可以得到以下形式:

最终可以得到离散形式如下:

通过式(7)的就可以得到输入信号经过陷波滤波器后的输出信号的完整离散形式。

陷波器的三个重要指标:

Wn:陷波中心频率

depth:陷波深度(中心频率处的衰减倍数)

Wb:陷波带宽(衰减到-3db时频率差值)

陷波器陷波系数K1,K2推导如下:

设陷波器的基本传递函数为:

①:将s域传递函数转化到G[jw](直接将s换为jw即可)

②:利用中心频率计算ζ1、ζ2之间的关系:(将w=wn带入方程求模)

由定义可知以下公式:

③:利用W1 W2的定义求解ζ1与陷波器指标的关系

因为

所以

即可得到公式:

由此我们得出两个重要公式:

k1=

k2=depth*k1

其中:B=2*pi*Wb

至此可求出离散化陷波滤波器6个系数,a1,a2,a3,b1,b2,b3,带入(7)式中可得差分方程。

至此只需要输入wn、depth、Wb通过关系式即可设计出一个陷波器。

 

设计了一个三参数滤波器,陷波频率100HZ,带宽40HZ,陷波幅度衰减100倍,matlab运行结果如下:

总结:

以上就是常见三种陷波滤波器离散化设计,前两种缺点已在文中说明,第三种是最值得使用的,matlab代码测试完美,有需要的请移步:

https://download.csdn.net/download/qczhao_10/87189346

博主花费了不少时间调试,感谢支持!

常见三种陷波滤波器(Notch Filter)的离散化设计相关推荐

  1. 温故知新(五)——三参数陷波滤波器离散化推导及MATLAB实现

    目录 1 三参数陷波滤波器 2 双线性变换 3 三参数计算 4 写在最后 1 三参数陷波滤波器 三参数型滤波器主要有以下一个参数: wn:滤波器带宽,单位rad/s=> * 2pi可转化为Hz: ...

  2. Simulink 窄带陷波滤波器(Notch filter)仿真到代码生成

    目录 窄带陷波滤波器(Notch filter)原理 设计一个Notch filter 生成代码 离散传递函数生成代码配置 差分方程生成代码 方程推导 构建差分方程模型 总结 窄带陷波滤波器(Notc ...

  3. 陷波滤波器(Notch Filter)的离散化设计

    陷波滤波器(Notch Filter)的离散化设计 符号说明 wbww_{bw}wbw​ - 陷波宽度,单位:rad/srad/srad/s wcw_cwc​ - 陷波中央频率,单位:rad/srad ...

  4. ArduCopter——ArduPilot——Harmonic Notch Filter(谐波陷波滤波器)

    版权声明:本文为博主原创博文,未经允许不得转载,若要转载,请说明出处并给出博文链接 之前,我们已经比较详细地学习和分析了Notch Filter陷波滤波器-->这里,点我. 那我们现在趁热打铁, ...

  5. 陷波滤波器(Notch filter)

    陷波滤波器: 指的是一种可以在某一个频率点迅速衰减输入信号,以达到阻碍此频率信号通过的滤波效果的滤波器.陷波滤波器属于带阻滤波器的一种,只是它的阻带非常狭窄,起阶数必须是二阶(含二阶)以上. 当带阻滤 ...

  6. 基于gnuradio的自适应陷波滤波器OOT模块(notch filter)

    基于gnuradio的自适应陷波滤波器OOT模块 Introduction 模型工具gr_modtool 创建OOT Module 创建OOT block 修改impl.h 修改impl.cc 修改a ...

  7. 自适应陷波滤波器(Notch滤波器)

    自适应滤波器以其自身调节参数的能力以及对信号和噪声的先验知识要求少的优势得到了广泛的应用. 自适应陷波滤波器相比于普通的数字滤波器具有较窄的阻带和通带,带外衰减较快.具有滤波输出(Y输出)和陷波输出( ...

  8. 【OpenCV 例程200篇】90. 频率域陷波滤波器

    [OpenCV 例程200篇]90. 频率域陷波滤波器 欢迎关注 『OpenCV 例程200篇』 系列,持续更新中 欢迎关注 『Python小白的OpenCV学习课』 系列,持续更新中 5.2 陷波滤 ...

  9. matlab提取电压基波分量,基于复系数陷波滤波器锁相环的电网正序分量检测

    摘  要: 为满足在电网电压不平衡情况下的系统控制需求,需要快速而准确地检测出基波正负序分量的幅值和相位.采用无限脉冲响应(IIR)复杂系数陷波滤波器结合锁相环来提取基波分量中的正序分量,利用MATL ...

最新文章

  1. 真赞!IDEA中这么玩MyBatis,让编码速度飞起!
  2. Python使用matplotlib函数subplot可视化多个不同颜色的折线图、自定义数据点的形状、自定义折线图的颜色
  3. PMcaff课堂:10年经验的产品大咖眼中的社交产品是这样的
  4. linux eclipse java_从Linux终端编译运行Eclipse Java项目
  5. 让Dapper支持读写分离
  6. 推荐免费的Windows Mobile截屏软件
  7. linux下vscode配置c++,使其和Visual Studio2019一样调试
  8. 大学excel题库含答案_Excel练习题及答案
  9. Centos7安装Caffe教程
  10. 上海房屋租赁合同登记备案证明办理流程
  11. 存储单元,字,MAR,存储字长,存储字,存储容量
  12. [画风突变] 藤原千花专属相册,C语言制作书记专属电子相册!
  13. qq修改实名认证已达上限_王者荣耀实名认证系统,怎么二次修改
  14. 黑苹果驱动hd4000
  15. qs计算机专业排名2017,2017qs世界大学学科专业排名top10完整版
  16. 生产制造企业用的ERP系统——流程管理
  17. 程序员如何提高自己的逻辑思维和技术水平
  18. html取服务器时间,客户端获取服务器时间记时
  19. 中国申请英国留学人数创新高,东伦敦成最受家长青睐“以房养学”地区
  20. 鸿蒙之境法有三乘,《神都夜行录》法有三乘副本组队通关攻略 游戏小伙伴都进来瞧瞧吧...

热门文章

  1. 百度网盘转存返回的错误码
  2. java防止闪屏代码_防止屏幕闪烁
  3. mysql简单应用题答案_MySQL数据库设计与应用题库免费期末考试2020答案
  4. Linux shell 脚本结构化命令 if-then
  5. 双非研究生选择延毕一年再战秋招,收割多家大厂SP岗
  6. QQ语音聊天记录也能轻松保存(转)
  7. selenium 隐藏窗口_查找Facebook用户的隐藏好友和社群
  8. if elseif执行顺序
  9. uniapp 文字无缝从右到左滚动
  10. Mysql中获取行号