目录
一.低通滤波器介绍
1 .概念
2 .原理
二、低通滤波器的设计
1 .技术指标
2 .设计步骤
三. MATLAB码
参考文献
另一方面,低通滤波器的概要1 .概念数字低通滤波器使低于截止频率的信号通过,高于截止频率的信号阻止通过。 通过这种操作方式,低通滤波器可以去除信号中的高频成分,留下低频成分,并可以去除高频噪声。
2 .原理数字滤波器可以由被称为滤波器系数函数的常数系数的线性差分方程表示。 通过建立理想的数字滤波器模型,重用数字滤波器的设计方法,可以得到与该模型相近的实际滤波器。
二、低通滤波器的设计1 .技术指标通带截止频率Wp、阻带起始频率Ws、通带纹波系数RP(db )、阻带衰减系数RP(db ) )。
2 .设计步骤(1)确定技术指标) Wp、Ws、Rp、Rs;
(2)计算滤波器最低阶n和3dB截止频率Wn;
(3)计算系数函数h ) z )的分子、分母多项式系数;
(4)计算系统函数h ) z )的幅度响应。
三、MATLAB码%%主函数-----------信号去噪clc; 清除全部关闭; %生成数据Fs=4000; %信号的采样频率t=(0:1/fs:(1-1/fs ) ); x1=3*sin(200*pi*TPI/4 ); x2=2*cos(1000*pi*t ); x=x1 x2; [m,n]=size(x ); % m表示信号点数Rp=4; %通带纹波系数Rs=20; %阻带衰减系数fp=250; %通带截止频率fs=400; %阻带开始频率[b,a]=buttord_filter(FP、fs、Rp、Rs、m、fs ); y=filter(b,a,x ); % filter函数实现IIR滤波器的直接滤波; subplot (2,1,1 ); plot(x ); subplot(2、1、2 ); plot(y; function [b,a]=buttord_filter(WP,Ws,Rp,Rs,n,Fs ) %设计微笑的锥形低通滤波器(滤波器设计参数(归一化角频率的计算) WP=WP/) %通带截止频率ws=ws/(fs/2; %阻带起始频率%滤波器的最低阶n和3dB截止频率Wn [n,wn]=buttord(WP,Ws,Rp,Rs ); %计算关系函数h(z )的分子,分母多项式系数[b,a]=Butter(n,Wn ); %计算关系函数h(z )的振幅响应: freqz(b ) b、a、计算点数、采样率); [H,f]=freqz(b,a,n,Fs ); %计算滤波器的相位PHA=angle(h ) *180/pi; %%低通滤波器振幅特性figure; subplot (2,1,1 ); plot(f,20*log10 ) ABS (h ); xlabel ('频率) Hz ) ); ylabel ('振幅) dB ) ); axis tight; 网格打开; subplot(2、1、2 ); plot(f,pha ); xlabel ('频率) Hz ) ); 相位; axis tight; 网格打开; 低通滤波器(title ); end参考文献[1]疏忽抽屉. MATLAB信号处理[M] .北京:清华大学出版社,2017:214233 .
[2]基于MATLAB的心电信号预处理