ButterWorth滤波
来源:互联网 发布:无法化成阶梯型矩阵 编辑:程序博客网 时间:2024/05/20 23:06
function y=highp(x,f1,f3,rp,rs,Fs)
%高通滤波
%使用注意事项:通带或阻带的截止频率的选取范围是不能超过采样率的一半
%即,f1,f3的值都要小于 Fs/2
%x:需要带通滤波的序列
% f 1:通带截止频率
% f 2:阻带截止频率
%rp:边带区衰减DB数设置
%rs:截止区衰减DB数设置
%FS:序列x的采样频率
% rp=0.1;rs=30;%通带边衰减DB值和阻带边衰减DB值
% Fs=2000;%采样率
%频率归一化
wp=2*f1/Fs;
ws=2*f3/Fs;
%{
% 设计切比雪夫滤波器; [n,wn]=cheb1ord(wp/pi,ws/pi,rp,rs);
% [bz1,az1]=cheby1(n,rp,wp/pi,'high');
%}
%% 设计巴特沃斯滤波器
%//注意切比雪夫和巴特沃斯传入参数的不同rp
[n,wn]=buttord(wp,ws,rp,rs);
[bz1,az1]=butter(n,wn,'high');
%{
%查看设计滤波器的曲线
% [h,w]=freqz(bz1,az1,256,Fs); %求离散系统频响特性的函数freqz()
% h=20*log10(abs(h));
% figure;plot(w,h);title('所设计滤波器的通带曲线');grid on;
%}
%% 滤波函数
y=filter(bz1,az1,x);
end
%高通滤波
%使用注意事项:通带或阻带的截止频率的选取范围是不能超过采样率的一半
%即,f1,f3的值都要小于 Fs/2
%x:需要带通滤波的序列
% f 1:通带截止频率
% f 2:阻带截止频率
%rp:边带区衰减DB数设置
%rs:截止区衰减DB数设置
%FS:序列x的采样频率
% rp=0.1;rs=30;%通带边衰减DB值和阻带边衰减DB值
% Fs=2000;%采样率
%频率归一化
wp=2*f1/Fs;
ws=2*f3/Fs;
%{
% 设计切比雪夫滤波器; [n,wn]=cheb1ord(wp/pi,ws/pi,rp,rs);
% [bz1,az1]=cheby1(n,rp,wp/pi,'high');
%}
%% 设计巴特沃斯滤波器
%//注意切比雪夫和巴特沃斯传入参数的不同rp
[n,wn]=buttord(wp,ws,rp,rs);
[bz1,az1]=butter(n,wn,'high');
%{
%查看设计滤波器的曲线
% [h,w]=freqz(bz1,az1,256,Fs); %求离散系统频响特性的函数freqz()
% h=20*log10(abs(h));
% figure;plot(w,h);title('所设计滤波器的通带曲线');grid on;
%}
%% 滤波函数
y=filter(bz1,az1,x);
end
0 0
- ButterWorth滤波
- 高通滤波,低通滤波,butterworth
- 利用matlab编写的Butterworth低通滤波程序
- Butterworth filter
- Matlab实现Butterworth滤波器
- butterWorth低通滤波器
- butterworth chebyshev bezier滤波器区别
- 经典滤波器之一:Butterworth滤波器
- [DSP] Butterworth (巴特沃斯)数字滤波器设计参考
- 滤波
- 滤波
- 滤波
- 滤波
- 滤波
- 滤波
- 滤波
- 软件滤波,硬件滤波
- 均值滤波,中值滤波,最大最小值滤波
- 2 最大子列和问题2
- 2017.4.18
- [最小割] SRM 590 div1 FoxAndCity
- RCTF-200 welpwn 地址泄露实例
- can总线控制器与单片机接口电路解析
- ButterWorth滤波
- LeetCode376. Wiggle Subsequence题解
- 为什么Android的AsyncTask不适合执行长时间操作的任务
- (25)折半插入排序
- 图片、文件表单上传以及异步上传
- ICN路由和转发流程
- Linux网络接口和VLAN配置
- C++ primer 5 5.41节练习
- 使用KMP算法求next数组