本文借鉴 小墨学FPGA 的文章 《零基础学FPGA( 三十二) 写在京城,多级FIR半带滤波器的FPGA实现 》
地址是: http://bbs.elecfans.com/forum.php?mod=viewthread&tid=509022
假设我们的合成信号是由100,300和600hz的正弦波相加而成的信号,并将其进行8bit量化,初始采样频率为6400hz,现在我需要将采样信号降为800hz,并要求获得100hz有用信号,要求滤波器的通阻带容限为0.1,设计这样一个FIR半带滤波器。
clear all
clc
close all
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 通带截止频率近似3KHz,阻带截止频率近似5KHz
fp = 100; % 要求最后一级滤波器的通带截止频率
fc = 300; % 要求最后一级滤波器的阻带截止频率
fs = 6400; %原始采样频率
fs_out = 800; % 要求降频后的采样频率
dev=0.1; %要求的滤波器通阻带容限
D=fs/fs_out; %求总抽取倍数
N=log2(D) ; %求滤波器级数
devi=dev/N; % 求每级滤波器通阻带容限
% 设计前2级半带滤波器
hin1 = firhalfband('minorder',fp/(fs/2),dev); %designs a lowpass minimum-order filter,
% with passband edge fp. The peak