matlab 陷波器,陷波滤波器—matlab实现
个人在处理工程信号的时候,会遇到电信号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实现相关推荐
- MATLAB:使用陷波器抑制LCL滤波器固有谐振
目录 0 背景 1 介绍 1.1 LCL滤波器的固有谐振 1.2 陷波器的功能 1.3 利用陷波器抑制LCL谐振 2 实例代码 PS: 0 背景 陷波器抑制谐振,这个方法不止适用于LCL滤波器,也适用 ...
- 【图像处理:频率域平滑与锐化】理想滤波器,巴特沃思滤波器,高斯滤波器
[频率域平滑.锐化滤波器]理想滤波器,巴特沃思滤波器,高斯滤波器 一.背景知识 二.理想滤波器原理及实现 1.理想低通滤波器 2.理想低通滤波器的实现: 3.理想高通滤波器: 三.巴特沃思滤波器原理及 ...
- matlab 陷波器,自己编的matlab陷波器的程序和大家分享
把以前学习Matlab编程的一些东东整理,觉得有些可能对大家有一点点帮助,拿出来和大家分享. 程序代码: %陷波器的设计 %陷波器的传输函数为 % B(1/z) (z-exp ...
- 陷波器的离散化及仿真验证
一.陷波器在连续域的传递函数 1.最基本的陷波器传函 (1) 其中,wo是所谓"中心频率",也就是你想要"陷掉 ...
- 啸叫抑制算法 ------陷波器
非常感谢好友 @Randolph .@泥沙的帮助,本文分享的简介很多都是来自于两位的讨论: 陷波器算法原理: 设计陷波器,降低啸叫频率点处增益,以破坏啸叫产生的增益条件,来达到啸叫抑制的效果. 陷波器 ...
- matlab ellip,IIR濾波器設計(調用MATLAB IIR函數來實現) | 學步園
% IIR濾波器設計 % 目的:設計一個採樣頻率為1000Hz.通帶截止頻率為50Hz.阻帶截止頻率為100Hz的低通濾波器,並要求通帶最大衰減為1dB,阻帶最小衰減為60dB. clc;clear; ...
- 50hz 60hz 级联 陷波器,卫星导航接收机陷波器级联抑制窄带干扰研究
限窄带干扰的有效抑制;8个II R 陷波器级联方法不仅具有较深陷波深度而且具有较为陡峭的带阻特 性,实现带限窄带干扰有效抑制的同时减小了卫星 导航信号的失真. 图7 干扰抑制前后信号的频谱 F i g ...
- iir陷波滤波器 matlab,IIR数字滤波器设计50Hz陷波器(MATLAB代码)
%% IIR陷波器设计 % 目的:设计一个陷波器阻带在50±1.5Hz以内,采样频率为400Hz的滤波器, % 并要求通带最大衰减为0.1dB,阻带最小衰减为60dB. clc; clear;clos ...
- matlab fdatool陷波器,基于MATLAB的简单50Hz数字陷波器.doc
基于MATLAB的简单50Hz数字陷波器 摘要:陷波器是一个二阶滤波器,它的幅度特性在ω=±ω0处为零,在其他频率上接近常数,是一个滤波单频干扰的滤波器.一般仪器都用50Hz的交流电源供电,因而信号中 ...
最新文章
- 如何让两个div在同一行显示?一个float搞定
- Android Studio编译NDK工程时,报ndk-build.cmd'' finished with non-zero exit value 2错误
- html5二次元插件,送你PS一键制作动漫二次元插件,80套背景+50种效果+20款中英字体...
- opencv学习笔记18:canny算子边缘检测原理及其函数使用
- linux 脚本 获取当前目录,Linux下获取脚本当前工作目录的一点感触
- rstudio 保存_Rstudio学习笔记
- 自定函数获取datagrid,datalist,rpeater控件中header,footer栏中控件
- Lnmp架构之PHP
- 启动zabbix server服务
- iOS 15 真机调试包 DeviceSupport
- NPOI导出一行多张图片
- android 自定义view实现仿QQ运动步数进度效果
- oracle move不能操作分区,关于Oracle分区表的move操作
- 棋牌游戏网站支付接口那些事儿
- 887. 鸡蛋掉落(困难)-动态规划
- linux操作系统 以下哪个命令,Linux认证考试基础试题及答案
- 看了他们的遭遇,我居然不想裸辞了
- 平方面积怎么测量?教你快速学会这个方法
- 比较好用的时间计算器有哪些?时间差计算器
- 微信小程序页面跳转刷新方式
热门文章
- 行为管理(锐捷云桌面2)
- 涂鸦 Wi-Fi SDK开发系列教程——4.烧录授权
- 浏览器兼容与图片整合
- 三轴运动控制系统西门子S7-200SMART型PLC和昆仑通泰触摸屏联机程序
- 8192EU无线网卡在ubuntu14.04下驱动编译错误解决
- 你知道怎样用Excel打印【条形码】吗?
- Catia V5高级图片造型视频教程 IMA模块宝马造型教程
- 1.仿QQ侧滑面板(对ViewGroup的自定义)
- 计算机考研面试翻译题,【北大计算机辅助翻译CAT考研】复试真题 复试参考书 复试面试真题...
- 从MTK的scat文件谈ROM和RAM的分配,管理和问题解决(拷贝)