matlab 生成滤波器系数

本文借鉴 小墨学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