个人在处理工程信号的时候,会遇到电信号50Hz及其倍频干扰特别强的情况,如何直接滤掉这些特定信号,而不影响周边有用信号呢,下面,就给大家推荐一个能精确剔除掉指定频率的滤波器——陷波滤波器。

input=textread('input2(t) Mar 21,2013 12-53-14.txt');%原信号

fs=2560;%采样频率

Ts=1/fs;

x=input(15*fs:17*fs,2);%提取信号第二列要处理的那一段,即要处理的时域信号。

f0=100;%要滤掉的信号,单位Hz

NLen=length(x(:,1));%信号长度

n=0:NLen-1;%陷波器的设计

apha=-2*cos(2*pi*f0*Ts);

beta=0.96;

b=[1 apha 1];

a=[1 apha*beta beta^2];

figure(1);

freqz(b,a,NLen,fs);%陷波器幅频特性显示

y=dlsim(b,a,x);%陷波器滤波处理,y即为滤波后的时域信号%对信号进行频域变换。

xfft=fft(x,NLen);

xfft=xfft.*conj(xfft)/NLen;​

y1=fft(y,NLen);

y2=y1.*conj(y1)/NLen;

figure(2);%滤除前后的信号对比。

subplot(2,2,1);plot(n,x);

grid;

xlabel('Time (s)');

ylabel('Amplitude');

title('Input signal');

subplot(2,2,3);plot(n,y);

grid;

xlabel('Time (s)');

ylabel('Amplitude');

title('Filter output');

subplot(2,2,2);

plot(n*fs/NLen,xfft);

axis([0 fs/2 min(xfft) max(xfft)]);

grid;xlabel('Frequency (Hz)');

ylabel('Magnitude (dB)');

title('Input signal');

subplot(2,2,4);

plot(n*fs/NLen,y2);

axis([0 fs/2 min(y2) max(y2)]);

grid;

xlabel('Frequency (Hz)');

ylabel('Magnitude (dB)');

title('Filter output');

我们看到100Hz的信号被干净利落的剔除了

matlab 陷波器,陷波滤波器—matlab实现相关推荐

  1. MATLAB:使用陷波器抑制LCL滤波器固有谐振

    目录 0 背景 1 介绍 1.1 LCL滤波器的固有谐振 1.2 陷波器的功能 1.3 利用陷波器抑制LCL谐振 2 实例代码 PS: 0 背景 陷波器抑制谐振,这个方法不止适用于LCL滤波器,也适用 ...

  2. 【图像处理:频率域平滑与锐化】理想滤波器,巴特沃思滤波器,高斯滤波器

    [频率域平滑.锐化滤波器]理想滤波器,巴特沃思滤波器,高斯滤波器 一.背景知识 二.理想滤波器原理及实现 1.理想低通滤波器 2.理想低通滤波器的实现: 3.理想高通滤波器: 三.巴特沃思滤波器原理及 ...

  3. matlab 陷波器,自己编的matlab陷波器的程序和大家分享

    把以前学习Matlab编程的一些东东整理,觉得有些可能对大家有一点点帮助,拿出来和大家分享. 程序代码: %陷波器的设计 %陷波器的传输函数为 %        B(1/z)       (z-exp ...

  4. 陷波器的离散化及仿真验证

    一.陷波器在连续域的传递函数 1.最基本的陷波器传函                              (1) 其中,wo​是所谓"中心频率",也就是你想要"陷掉 ...

  5. 啸叫抑制算法 ------陷波器

    非常感谢好友 @Randolph .@泥沙的帮助,本文分享的简介很多都是来自于两位的讨论: 陷波器算法原理: 设计陷波器,降低啸叫频率点处增益,以破坏啸叫产生的增益条件,来达到啸叫抑制的效果. 陷波器 ...

  6. matlab ellip,IIR濾波器設計(調用MATLAB IIR函數來實現) | 學步園

    % IIR濾波器設計 % 目的:設計一個採樣頻率為1000Hz.通帶截止頻率為50Hz.阻帶截止頻率為100Hz的低通濾波器,並要求通帶最大衰減為1dB,阻帶最小衰減為60dB. clc;clear; ...

  7. 50hz 60hz 级联 陷波器,卫星导航接收机陷波器级联抑制窄带干扰研究

    限窄带干扰的有效抑制;8个II R 陷波器级联方法不仅具有较深陷波深度而且具有较为陡峭的带阻特 性,实现带限窄带干扰有效抑制的同时减小了卫星 导航信号的失真. 图7 干扰抑制前后信号的频谱 F i g ...

  8. iir陷波滤波器 matlab,IIR数字滤波器设计50Hz陷波器(MATLAB代码)

    %% IIR陷波器设计 % 目的:设计一个陷波器阻带在50±1.5Hz以内,采样频率为400Hz的滤波器, % 并要求通带最大衰减为0.1dB,阻带最小衰减为60dB. clc; clear;clos ...

  9. matlab fdatool陷波器,基于MATLAB的简单50Hz数字陷波器.doc

    基于MATLAB的简单50Hz数字陷波器 摘要:陷波器是一个二阶滤波器,它的幅度特性在ω=±ω0处为零,在其他频率上接近常数,是一个滤波单频干扰的滤波器.一般仪器都用50Hz的交流电源供电,因而信号中 ...

最新文章

  1. 如何让两个div在同一行显示?一个float搞定
  2. Android Studio编译NDK工程时,报ndk-build.cmd'' finished with non-zero exit value 2错误
  3. html5二次元插件,送你PS一键制作动漫二次元插件,80套背景+50种效果+20款中英字体...
  4. opencv学习笔记18:canny算子边缘检测原理及其函数使用
  5. linux 脚本 获取当前目录,Linux下获取脚本当前工作目录的一点感触
  6. rstudio 保存_Rstudio学习笔记
  7. 自定函数获取datagrid,datalist,rpeater控件中header,footer栏中控件
  8. Lnmp架构之PHP
  9. 启动zabbix server服务
  10. iOS 15 真机调试包 DeviceSupport
  11. NPOI导出一行多张图片
  12. android 自定义view实现仿QQ运动步数进度效果
  13. oracle move不能操作分区,关于Oracle分区表的move操作
  14. 棋牌游戏网站支付接口那些事儿
  15. 887. 鸡蛋掉落(困难)-动态规划
  16. linux操作系统 以下哪个命令,Linux认证考试基础试题及答案
  17. 看了他们的遭遇,我居然不想裸辞了
  18. 平方面积怎么测量?教你快速学会这个方法
  19. 比较好用的时间计算器有哪些?时间差计算器
  20. 微信小程序页面跳转刷新方式

热门文章

  1. 行为管理(锐捷云桌面2)
  2. 涂鸦 Wi-Fi SDK开发系列教程——4.烧录授权
  3. 浏览器兼容与图片整合
  4. 三轴运动控制系统西门子S7-200SMART型PLC和昆仑通泰触摸屏联机程序
  5. 8192EU无线网卡在ubuntu14.04下驱动编译错误解决
  6. 你知道怎样用Excel打印【条形码】吗?
  7. Catia V5高级图片造型视频教程 IMA模块宝马造型教程
  8. 1.仿QQ侧滑面板(对ViewGroup的自定义)
  9. 计算机考研面试翻译题,【北大计算机辅助翻译CAT考研】复试真题 复试参考书 复试面试真题...
  10. 从MTK的scat文件谈ROM和RAM的分配,管理和问题解决(拷贝)