[求助]如何设计一带通(15KHZ-25KHZ)的滤波器

xiaming · 发表于 2006-5-18 18:40

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

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

x

如何设计一带通(15KHZ-25KHZ)的数字滤波器,来消除信号噪声,望各位大侠帮忙,在线等

yangzj · 发表于 2006-5-18 18:54

%带通滤波 %x-输入信号 %fs-输入信号采样频率 %N-返回点数 %M-滤波器半阶数 %fl-下截止频率 %fh-上截止频率 function xl=beltfilter(x,fs,N,M,fl,fh)
k=1:M; w=0.5+0.5*cos(pi*k/M);
wl=2*pi*fl/fs; wh=2*pi*fh/fs; h(1)=(wh-wl)/pi; h(2:M+1)=(sin(wh*k)-sin(wl*k))./(pi*k).*w;
for k=1:N kk=k-1+M; xl(k)=x(kk+1)*h(1)+sum(h(2:M+1).*(x(kk+2:kk+M+1)+x(kk:-1:kk-M+1))); end

xiaming · 发表于 2006-5-18 19:25

好象不行啊

yangzj · 发表于 2006-5-18 19:39

不会吧,我验证过的哦

xiaming · 发表于 2006-5-18 19:50

M怎么取合适呢,对数字信号有用吗

yangzj · 发表于 2006-5-18 19:55

越大滤波器性能越佳,但计算量增大

xiaming · 发表于 2006-5-18 20:04

数字信号能用吗

yangzj · 发表于 2006-5-18 20:09

我这本来就是数字滤波器呀. 还有要注意的是返回的序列是去掉了原序列的前M点和后M点,因此,如果要求初始相位的话要进行相应的补偿.

xiaming · 发表于 2006-5-18 20:33

??? Error using ==> unknown Matrix dimensions must agree.
Error in ==> beltfilter at 20 xl(k)=x(kk+1)*h(1)+sum(h(2:M+1).*(x(kk+2:kk+M+1)+x(kk:-1:kk-M+1))); 怎么报错啊

yangzj · 发表于 2006-5-18 20:50

你把参数x转置一下

lz1736 · 发表于 2006-5-19 11:58

出现错误： ??? Error using ==> unknown Matrix dimensions must agree.
Error in ==> kong11 at 16419 xl(k)=x(kk+1)*h(1)+sum(h(2:M+1).*(x(kk+2:kk+M+1)+x(kk:-1:kk-M+1))); x转置后还是出错， ??? Index exceeds matrix dimensions.
Error in ==> kong11 at 16419 xl(k)=(x(kk+1)*h(1))'+sum(h(2:M+1).*((x(kk+2:kk+M+1))'+(x(kk:-1:kk-M+1))'));

yangzj · 发表于 2006-5-19 12:01

你的序列长度不够长，最少为N+2*M

账号		自动登录	找回密码
密码			我要加入