cic_filter.m

来自「关于软件无线电中CIC滤波器的实现与仿真」· M 代码 · 共 41 行

%积分梳状滤波器的设计
clear;
clc;
%积分滤波器的响应

b1=1;
a1=[1 -1];
% freqz(b1,a1,'whole');

hold on;
%梳状滤波器的响应
D=9;
b2=[1 zeros(1,D-1) -1];
a2=1;
% freqz(b2,a2,'whole');

%积分梳状滤波器响应
b3=b2;
a3=a1;
% freqz(b3/D,a3,'whole')

%两只内插因子为9的CIC级联响应
b4=conv(b3,b3);
a4=conv(a3,a3);
% freqz(b4/D^2,a4,'whole')

%三只内插因子为9的CIC级联响应
b5=conv(b4,b3);
a5=conv(a4,a3);
% freqz(b5/D^3,a5,'whole');

%四只内插因子为9的CIC级联响应
b6=conv(b4,b4);
a6=conv(a4,a4);
freqz(b6/D^4,a6,'whole');

%五只内插因子为9的CIC级联响应
b7=conv(b6,b3);
a7=conv(a6,a3);
freqz(b7/D^5,a7,'whole');