时常仰望星空,做到脚踏实地!
matlab提供了用于设计贝塞尔滤波器的函数:besself。调用时,[b,a] = besself(n,Wn,'ftype')用于设计高通和带阻滤波器,默认[b,a] = besself(n,Wn)是设计贝塞尔低通滤波器。截止频率为Wn(单位为弧度/秒),切记转化为数字频率,不然设计容易出错。
心电图处理基线漂移时一般选用高通0.5HZ来滤除基线漂移。本例设计一个5阶贝塞尔高通滤波器。设计该滤波器的过程如下:
1.设计模拟高通贝塞尔滤波器:[b,a]=besself(5,2*pi*0.5,'high');
2.[num,den]=bilinear(b,a,500); %500为采样频率