matlab cic插值与抽取,CIC抽取滤波器和插值滤波器

一、CIC抽取滤波器

r =2;                  %抽取因子

hm = mfilt.cicdecim(r);

fs = 44.1e3;          %原始信号采样频率44.1khz

n = 1:10239;          % 10240个采样点

x = sin(2*pi*1e3/fs*n);   %原始信号

y_fi = filter(hm,x);      % 得到抽取的5120个采样点

x=double(x);

y=double(y_fi);

y=y/max(abs(y));

stem(n(1:44)/fs,x(2:45)); hold on;

stem(n(1:22)/(fs/r),y(3:24),'r','filled');

xlabel('时间(sec)'); ylabel('信号量');

二、CIC插值滤波器

r =2;                  %抽取因子

hm = mfilt.cicinterp(r);

fs = 22.05e3;          %原始信号采样频率22.05khz

n = 0:5119;          % 5120个采样点

x = sin(2*pi*1e3/fs*n);   %原始信号

y_fi = filter(hm,x);      % 得到抽取的  个采样点

x=double(x);

y=double(y_fi);

y=y/max(abs(y));

stem(n(1:22)/fs,x(1:22),'filled'); hold on;

stem(n(1:44)/(fs*r),y(4:47),'r');

xlabel('时间(sec)'); ylabel('信号量');