iwl = 16;  %CIC濾波器輸入信號的字長
owl = 32;  %CIC濾波器輸出信號的字長
wlps = 32;  %CIC濾波器每級濾波器中每個字的比特數
decf = 4;  %CIC濾波器的抽取率
diffd = 1;  %CIC濾波器中的差分時延
numsecs = 3;  %CIC濾波器級數
l = 2;      %半帶濾波器抽取因子
n = 16;       %半帶濾波器階數
fp = 0.0001;     %半帶濾波器的通帶歸一化頻率
decm=decf*l;
len=length(x);
x1 = CICdec(decf,diffd,numsecs,x,iwl,owl,wlps);   %調用子函數CICdec進行CIC的抽取
x2 = halfdec(n,fp,l,x1);          %調用子函數halfdec進行半帶濾波器的抽取
x3 = filter(h,1,x2);       %經過FIR濾波器，h為濾波器系數
stem(real(x(1:decm:decm*len)));hold on;
stem(real(x3(1:len)),'r.');