声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1441|回复: 9

[滤波] 为什么不同滤波器,我滤波后的结果是一样的

[复制链接]
发表于 2010-3-17 22:19 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我有一个信号先是5Hz正弦持续1s,而后持续1s是15Hz正弦,然后又是1s时间的5Hz正弦。
我现在用8阶巴特沃斯滤波器滤波,第一次是带通[14 16],第二次是带通[15 17],但是最后出来的结果竟然是一样的,代码如下:

%先产生信号
fs=500;     %采样频率fs=500Hz.
t1=0:1/fs:1-1/fs;   
f1=sin(2*pi*5*t1);  
t2=1:1/fs:2-1/fs;  
f2=sin(2*pi*15*t2);  
t3=2:1/fs:3-1/fs;   
f3=sin(2*pi*5*t3);  
Data1=[f1 f2 f3];

%buterworth带通14-16Hz
[b,a]=butter(8,16/250,'low');
filter_Data11=filter(b,a,Data1);
[b,a]=butter(8,14/250,'high');
filter_Data2=filter(b,a,filter_Data11);

%buterworth带通15-17Hz
[b,a]=butter(8,17/250,'low');
filter_Data1=filter(b,a,Data1);
[b,a]=butter(8,15/250,'high');
filter_Data1=filter(b,a,filter_Data1);

subplot(3,1,1);plot(Data1);subplot(3,1,2);plot(filter_Data1,'r');subplot(3,1,3);plot(filter_Data2);
结果是,两个滤波结果基本一样,感觉好奇怪啊,请高手指点!

原始数据及滤波结果

原始数据及滤波结果
回复
分享到:

使用道具 举报

 楼主| 发表于 2010-3-17 22:22 | 显示全部楼层
我感觉用[15 17]的带通滤波,那信号应该被滤的没剩下什么了,为什么还有这么多?
发表于 2010-3-18 10:53 | 显示全部楼层
[16  18]再试试看
 楼主| 发表于 2010-3-25 20:45 | 显示全部楼层
16 18有差别,但是为什么14 16与15 17就没有什么差别,请指点
发表于 2010-3-25 21:42 | 显示全部楼层

回复 地板 tjucruiser 的帖子

14 16与15 17两种情况15Hz都在通频带内。即使16 18 对于15也不是滤干净的。
发表于 2010-3-26 18:30 | 显示全部楼层
第二次是带通[15 17],15Hz仍在通带之内,但按Butterwoth滤波器将有3dB的衰减,所以笫2张图的幅值比笫3张图小一些。
发表于 2010-3-27 12:13 | 显示全部楼层
楼主的方式是先16Hz低通, 再14Hz高通; 与直接14-16Hz带通有何不同?

fs=500; t1=0:1/fs:1-1/fs; f1=sin(2*pi*5*t1); t2=1:1/fs:2-1/fs;f2=sin(2*pi*15*t2);
t3=2:1/fs:3-1/fs; f3=sin(2*pi*5*t3); Data1=[f1 f2 f3]; fL=14; fH=16;

[b,a]=butter(8,fH/250,'low'); filter_Data11=filter(b,a,Data1);
[b,a]=butter(8,fL/250,'high'); filter_Data2=filter(b,a,filter_Data11);

[b,a]=butter(8,[fL/250,fH/250]); filter_Data1=filter(b,a,Data1);

subplot(3,1,1);plot(Data1);subplot(3,1,2);plot(filter_Data1,'r');subplot(3,1,3);plot(filter_Data2);

[ 本帖最后由 ChaChing 于 2010-3-27 16:28 编辑 ]
 楼主| 发表于 2010-3-27 15:37 | 显示全部楼层
楼上的,按照你的说法写的程序,结果都被滤波成直线了。请指点。
fs=500;     %采样频率fs=500Hz.
t1=0:1/fs:1-1/fs;   
f1=sin(2*pi*5*t1);  
t2=1:1/fs:2-1/fs;  
f2=sin(2*pi*15*t2);  
t3=2:1/fs:3-1/fs;   
f3=sin(2*pi*5*t3);  
Data1=[f1 f2 f3];

%buterworth带通14-16Hz
[b,a]=butter(8,[14/250,16/250]); filter_Data1=filter(b,a,Data1);
[b,a]=butter(8,[15/250,17/250]); filter_Data2=filter(b,a,Data1);

subplot(3,1,1);plot(Data1);
subplot(3,1,2);plot(filter_Data1);
subplot(3,1,3);plot(filter_Data2);

另外向各位请教,滤波器的特性steepness40 dB/decade是什么意思,还有别的什么滤波器特性吗
发表于 2010-3-28 22:02 | 显示全部楼层

回复 8楼 tjucruiser 的帖子

椭圆滤波器不是直接有带通吗?
发表于 2010-3-29 00:27 | 显示全部楼层

回复 9楼 ydiok 的帖子

我7F不就是请教"带通"与"先低通再高通"两者差异?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 04:43 , Processed in 0.066329 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表