肌电信号处理分析

为带通滤波器，需设置半阶数、高低截止频率和采样频率。保存为function函数的m文件。rar中为一肌电信号。

%带通滤波
%x-输入信号
%fs-输入信号采样频率
%M-滤波器半阶数
%fl-下截止频率
%fh-上截止频率
%xl-返回的滤波结果其长度为原序列长度减少2M去掉了前M和后M个不正确的点
function xl=beltfilter（xfsMflfh）
N=length（x）-2*M;
[i j]=size（x）;
if i~=1
    x=x‘;
end
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

if i~=1
    xl=xl‘;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    1001726  2008-11-06 15:45  为带通滤波器，需设置半阶数、高低截止频率和采样频率。保存为function函数的m文件。rar中为一肌电信号\ab.txt

     文件        525  2008-12-02 13:25  为带通滤波器，需设置半阶数、高低截止频率和采样频率。保存为function函数的m文件。rar中为一肌电信号\beltfilter.m

     文件        218  2007-06-05 03:14  为带通滤波器，需设置半阶数、高低截止频率和采样频率。保存为function函数的m文件。rar中为一肌电信号\www.pudn.com.txt

     目录          0  2013-04-23 20:59  为带通滤波器，需设置半阶数、高低截止频率和采样频率。保存为function函数的m文件。rar中为一肌电信号

----------- ---------  ---------- -----  ----

              1002469                    4