l=2;  %半帶濾波器的插值因子
n = 16;     %半帶濾波器階數
fp = 0.0001;     %半帶濾波器的通帶歸一化頻率
r=4;     %CIC濾波器的插值因子
interp1=r*l;   %總的插值因子
hb1 = halfinterp(n,fp,l,x);        %調用子函數halfinterp進行半帶濾波器的插值
x1 = CICinterp(r,hb1);            %調用子函數CICinterp進行CIC濾波器的插值
stem(real(x(1:len)));hold on;
stem(real(x3(1:interp1:interp1*len)),'r.');