Matlab调用函数实现CIC滤波器

matlab里设计cic滤波器的函数有以下两种：

1. fdesign.decimator

       例如：设定好采样频率Fs, 信号带宽Fp, 阻带衰减As, 差分时延m及降采样比D就可以得到cic滤波器的传输函数

d1 =fdesign.decimator(D,'CIC',m,Fpass,As,Fs);

Hcic =design(d1);

 2.mfilt.cicdecim (fixed-point CIC decimator, mfilt是matlab里专门用来设计多速率信号处理滤波器的一套函数)

hm =mfilt.cicdecim(decimation_factor,differential_delay,NumberofSections);

      decimation_factor为降采样比，differential_delay同上为差分时延，NumberofSections为cic滤波器的节数，与第一个函数相比，这个函数没有规定采样滤波，通带宽度、阻带衰减等

在信号处理中，信号发射时，信号通过载波，调制，以电磁波的形式发射出来

在接收端，射频信号通过天线接收，超外差式等方式进行处理，变频到中频

然后通过数字技术对中频信号进行处理，此时的中频信号采样率很高，多达几十兆赫

而要解调出来的信号很小，此时就要经过滤波来得到我们需要的信号

如果直接设计滤波器，阶数会非常庞大，软件，硬件设备都承受不了

因此，一般常用的技术是进行下变频处理，其主要技术就是通过抽取滤波进行下变频

常用的滤波器就是CIC抽取, 希望可以帮助一些坛友解决基本的问题。

下面以一个实例，对一个采样频率为45.5MHz的信号进行14倍的抽取滤波，同时进行

补偿滤波器的设计，并给出CIC滤波器、补偿滤波器和级联后的频谱图

1. %  Design a minimum-order CIC compensator that compensates...
2. %  for the droop in the passband for the CIC decimator.
3. Fs = 45.5e6;     % Input sampling frequency
4. Fpass = 0.5e6;   % Frequency band of interest
5. D = 14;          % Decimation factor of CIC
6. d1 = fdesign.decimator(D,'CIC',1,Fpass,65,Fs); %design a cic filter
7. Hcic = design(d1);
8. Hd(1) = cascade(dfilt.scalar(1/gain(Hcic)),Hcic);
9. d2 = fdesign.ciccomp(Hcic.DifferentialDelay, ...
10.     Hcic.NumberOfSections,Fpass,1.625e6,.005,66,Fs/D); % design a cic compensator filter
11. Hd(2) = design(d2);
12. fcfwrite([Hcic Hd(2)],'CICdesciption','dec'); % 其中，生成的.fcf文件描述滤波器的结构
13. hvt=fvtool(Hd(1),Hd(2),cascade(Hd(1),Hd(2)),'Fs',[Fs Fs/D Fs], ...   % plot whole response
14.           'ShowReference', 'off');
15. legend(hvt, 'CIC','CIC compensator', 'Whole response','Location', 'Northeast');

复制代码