function out = CICinterp(r,x); 
Hm = mfilt.cicinterp(r);    %使用默認的CIC階數和差分時延，以及字長
x_real = real(x);
x_imag = imag(x);
out_real = filter(Hm,x_real); %實部和虛部分開濾波
out_imag = filter(Hm,x_imag);
outt = (double(out_real) + j*double(out_imag))/gain(Hm);
out=outt*sqrt(var(x)/var(outt));   %使得輸入和輸出的功率相同以便比較