巴特沃斯滤波器带通滤波阶数问题

神瑛侍者

想按照图片所示的要求进行滤波，遇到问题，请各位帮看看，指点指点，多谢啦~

到[b a]=butter(n,wn,ft);这一步提示Filter order too large.看来一下，阶数达到了9088，不知道怎么改怎么修改啊？Acce(:,1)这一列的数据见附件txt。 Acce.txt (193.97 KB, 下载次数: 1)

2013-11-30 10:09 上传

点击文件名下载附件

1. 
   
2. fs=512;
   
3. wp=[0.1,25];
   
4. ws=[0.09,25.01];
   
5. rp=2;
   
6. rs=30;
   
7. ft='bandpass';
   
8. [n wn]=buttord(wp/(fs/2),ws/(fs/2),rp,rs);
   
9. [b a]=butter(n,wn,ft);
   
10. [h w]=freqz(b,a,1024,fs);
    
11. l=length(Acce(:,1));
    
12. t=0:1/fs:(l-1)/fs;
    
13. y=filter(b,a,x);
    
14. subplot(2,1,1);
    
15. plot(w,abs(h));
    
16. xlabel('Frequency(Hz)')
    
17. ylabel('Amplitude')
    
18. legend('Frequency-Amplitude Curve')

复制代码

神瑛侍者

滤波后感觉时程曲线发生了时间上的平移，请问是什么原因啊，怎么解决啊？

江树空

正常，勿惊慌！你用的是IIR滤波器，IIR滤波器不经过额外处理，是不能保证其相位的。

神瑛侍者

江树空 发表于 2013-11-30 20:01

                               
登录/注册后可看大图

正常，勿惊慌！你用的是IIR滤波器，IIR滤波器不经过额外处理，是不能保证其相位的。

哦哦，谢谢，额外处理该怎么处理啊，求教~

lragy

零相位滤波！！！！！！！！！！！！！！！！！

神瑛侍者

lragy 发表于 2013-12-1 01:23

                               
登录/注册后可看大图

零相位滤波！！！！！！！！！！！！！！！！！

谢谢啊，正解啊~