[求助]论坛中的一段滤波算法,请求高手转化为C语言,能够脱离matlab运行
[求助]论坛中的一段滤波算法,请求高手转化为C语言,能够脱离matlab运行=====================================================================
%带通滤波
%x-输入信号
%fs-输入信号采样频率
%M-滤波器半阶数
%fl-下截止频率
%fh-上截止频率
%xl-返回的滤波结果,其长度为原序列长度减少2M,去掉了前M和后M个不正确的点
function xl=bp(x,fs,M,fl,fh)
N=length(x)-2*M;
=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
页:
[1]