matlab设计CIC滤波器

4阶，32倍插值CIC

% CIC 1M->32M   

data=[1;2;3;4;5;6;7;8;9];

    % Comb
    C1=diff([0;data]);
    C2=diff([0;C1]);
    C3=diff([0;C2]);
    C4=diff([0;C3]);
    % Interpolation and Intergrator
    IQ_temp=cumsum([0;C4]);
    E1=zeros(length(IQ_temp)*32,1);
    for i=1:32
        E1(i:32:end)=IQ_temp;
    end
    E2=cumsum([0;E1]);
    E3=cumsum([0;E2]);
    E4=cumsum([0;E3]);
    data_cic=round(E4/2^15);