通带截止频率近似3KHz,阻带截止频率近似5KHz
fs = 16e3; % 采样频率
n = 0:5119; % 采样点数
x = sin(2* pi*1e3/fs* n); %输入原始信号
b =
firhalfband(18,0.00245,'dev'); % 获得设计的半带滤波器抽头系数
impz(b); %
半带滤波器的冲击响应
h = mfilt.firdecim(2,b)
figure
y_fi = filter(h,x);
x = double(x);
y = double(y_fi);
y = y/max(abs(y));
stem(n(1:44)/fs,x(1:44),'r');
xlabel('时间(sec)');
ylabel('输入信号');
hold on
figure
stem(n(1:2:44)/fs,y(5:26),'filled');
xlabel('信号值');
ylabel('时间(sec)');