带通滤波器

yangzj

前面发的滤波器，有人使用有些考虑不全的地方。现重新发一个改进后的。

1. %带通滤波
   
2. %x-输入信号
   
3. %fs-输入信号采样频率
   
4. %M-滤波器半阶数
   
5. %fl-下截止频率
   
6. %fh-上截止频率
   
7. %xl-返回的滤波结果,其长度为原序列长度减少2M,去掉了前M和后M个不正确的点
   
8. function xl=beltfilter(x,fs,M,fl,fh)
   
9. 
   
10. N=length(x)-2*M;
    
11. [i j]=size(x);
    
12. if i~=1
    
13.     x=x';
    
14. end
    
15. k=1:M;
    
16. w=0.5+0.5*cos(pi*k/M);
    
17. 
    
18. wl=2*pi*fl/fs;
    
19. wh=2*pi*fh/fs;
    
20. h(1)=(wh-wl)/pi;
    
21. h(2:M+1)=(sin(wh*k)-sin(wl*k))./(pi*k).*w;
    
22. 
    
23. for k=1:N
    
24.     kk=k-1+M;
    
25.     xl(k)=x(kk+1)*h(1)+sum(h(2:M+1).*(x(kk+2:kk+M+1)+x(kk:-1:kk-M+1)));
    
26. end
    
27. 
    
28. if i~=1
    
29.     xl=xl';
    
30. end

复制代码

wangwu1976

今天试 过了您的这个滤波器,对高频滤波效果很好,但对低频几乎没什么作用,不知道怎么回事,还请您指教

科技在线

本来滤波的范围就是有限的

小黑

这个程序怎么用？

不亦乐乎

原帖由 科技在线 于 2008-12-2 22:05 发表

                               
登录/注册后可看大图

本来滤波的范围就是有限的

没看懂，人家不是带通滤波吗？

Jnny_CN

今天用了一下，感觉很好，谢谢

请问哪位高手能够帮忙把楼主的程序改写为C语言呢？
我弄了半天不得要领，谢谢

[ 本帖最后由 ChaChing 于 2009-3-20 14:16 编辑 ]

Jnny_CN

有个问题：您这个滤波器在通带内的衰减系数不为1，按理说在通带内信号应该不衰减才对阿

参数如下
fs=20000;
fl=900;
fh=1100;
M=50;

希望楼主能帮忙看下

Jnny_CN

在上个采样频率40K的2k正弦波的数据对比

带通后衰减得非常厉害

fs=40000;
fl=1900;
fh=2100;
M=50;

[ 本帖最后由 Jnny_CN 于 2009-3-20 13:41 编辑 ]

purplejack

没有用matlab中自带的滤波器函数，厉害

太阳vib

弱弱的问一下，带通滤波和高通滤波有什么区别么？？？