cicdecimation.m

来自「OFDM系统仿真」· M 代码 · 共 53 行

%-------cic decimation filter 
%-------time   :2007-12-06
%-------name   : Wanyiping
%----R ratio change factor
%----M delay order of comb filter
%----N cascade stages num.

function    data_out=cicdecimation(data,R,M,N);

data_len=max(size(data));

s(1,1)=data(1);
for i=2:data_len
    s(i,1)=data(i)+s(i-1,1);
end

if N>1
    for j=2:N
        s(1,j)=s(1,j-1);
        for i=2:data_len
            s(i,j)=s(i,j-1)+s(i-1,j);
        end
    end
end

for i=1:data_len
    data_intp(i)=s(i,N);
end

down_data=downsample(data_intp,R,2);

down_data_len=max(size(down_data));

for i=1:M
    ds(i,1)=down_data(i);
end
for i=M+1:down_data_len
    ds(i,1)=down_data(i)-down_data(i-M);
end

if N>1
    for j=2:N
        for k=1:M
            ds(k,j)=ds(k,j-1);
        end
        for i=M+1:down_data_len
            ds(i,j)=ds(i,j-1)-ds(i-M,j-1);
        end
    end
end

data_out = ds(:,N);