简单的CIC滤波器matlab程序

%%%%%%%%%
%Hogenauer CIC
%编写时间：2018/6/15   15：14
%%%%%%%%%
f1=1000;       %信号频率为1KHz
f2=30000;      %信号频率为40KHz
f3=40000;
Fs=200*f1;     %采样频率为200KHz
M=5;           %CIC滤波器阶数
D=M;           %抽取因子
C=3;           %滤波器级数
N=10;          %量化位数

%产生信号
t=0:1/Fs:0.02;
c1=2*pi*f1*t;    %信号频率为1KHz
c2=2*pi*f2*t;     %信号频率为30KHz
c3=2*pi*f3*t;     %信号频率为40KHz
si=sin（c1）+sin（c2）+sin（c3）;%产生正弦波
 
%%%%%%%%%%%%%%%%%%%%%%%%%%
%3级积分器
a=[1