matlab低通滤波器库函数代码_FIR滤波器窗函数设计法详细步骤以及Matlab代码

clear

close all

global Fs

Fs = 360;

load '118m.mat'%mit数据库第118条数据

signal = val(1,100000:111600)/200;

%% 采用FIR I型设计20Hz以下的低通滤波器

fp=14; fs=18;

detap = 0.01;detas = 0.01;

[M,beta] = selectFirFilterN(fp,fs,detap,detas);

N = M+1;

w = kaiser(N,beta);

hd = FIRItypeIdealpulse(fp,fs,N,'low');

h = hd.*w';

% 设计的滤波器

omega = linspace(0,pi,512);

mag = freqz(h,[1],omega);

figure

plot(omega/(2*pi)*Fs,20*log10(abs(mag)));

title('FIR低通(14hz以下)滤波器频率相应');

xlabel('频率');

ylabel('增益(dB)');

%% 采用FIR I型设计8Hz以上的高通滤波器

fp2 = 8; fs2=4;

detap2 = 0.01; detas2 = 0.01;

[M2,beta2] = selectFirFilterN(fp2,fs2,detap2,detas2);

N2 = M2