matlab 级联cic,Matlab中CIC滤波器的应用
CIC滤波器基本原理
CIC(积分梳状级联)滤波器是工程上经常用的滤波器,因为CIC滤波器不需要乘法,CIC滤波器往往在级联抽取滤波器的第一级和级联插值滤波器的最后一级。这一节我们以CIC抽取滤波器为例来讲CIC滤波器的使用方法,CIC滤波器的数学模式为:
image-20201002173656346
可以将其分解为两部分:
image-20201002173729012
其中第一级是积分器,就是做累加;第二级是级联梳状滤波器。
image-20201002173755692
CIC滤波器有三个参数:
R:抽取倍数
M:差分延迟,在工程中一般取1或者2
N:滤波器级数,一级CIC的阻带有13.46dB的衰减,N级就是有13.46*N dB的衰减。
若R=8,M=1,N=5,那么CIC滤波器的频响如下所示:
image-20201002192819206
Matlab中CIC滤波器的使用
在Matlab中,有两个函数可以生成CIC滤波器。我们以CIC抽取滤波器为例,一个是fdesign.decimator,第二个是dsp.CICDecimator。下面我们分别来看这两个函数的使用方法。
生成滤波器模型
首先来看如何用fdesign.decimator生成一个CIC滤波器和其补偿滤波器,并查看CIC滤波器、补偿滤波器、级联滤波器的频响。下面代码中hcic即为构造的CIC滤波器,抽取倍数为4,差分延时为1;cic_comp为CIC补偿滤波器;其中Ap表示通带衰减,Astp表示阻带衰减。
Fs = 1; % sample rateR = 4; % decimator factorD = 1; % differential delayFp = 0.05; % pass bandFstp = 0.075; % stop bandAp = 0.1; % attenuation in pass bandAstp = 60; % attenuation in stop bandhcic = design(fdesign.decimator(R,'cic',D, Fp, Astp, Fs),'SystemObject',true);cic_comp = design(fdesign.ciccomp(hcic.DifferentialDelay, ...hcic.NumSections,Fp,Fstp,Ap,Astp,Fs/R), 'SystemObject',true);fvtool(hcic,cic_comp,...cascade(hcic,cic_comp),'ShowReference','off','Fs',[Fs Fs/R Fs])legend('CIC Decimator','CIC Compensator','Resulting Cascade Filter');
image-20210207151607757
dsp.CICDecimator的用法也非常类似,我们可以直接指定阶数;CIC补偿滤波器也直接使用dsp.CICCompensationDecimator函数来设计;
Fs = 1; % sample rateR = 4; % decimator factorD = 1; % differential delayN = 5; % number of stageFp = 0.05; % pass bandFstp = 0.075; % stop bandAp = 0.1; % attenuation in pass bandAstp = 60; % attenuation in stop bandCICDecim = dsp.CICDecimator(R, D, N);CICCompDecim = dsp.CICCompensationDecimator(CICDecim, ...'DecimationFactor',2,'PassbandFrequency',Fp, ...'StopbandFrequency',Fstp,'SampleRate',Fs/R);fvtool(CICDecim,CICCompDecim,...cascade(CICDecim,CICCompDecim),'ShowReference','off','Fs',[Fs Fs/R Fs])legend('CIC Decimator','CIC Compensator','Resulting Cascade Filter');
频响如下:
image-20210207161734382
滤波操作
如果要用上面生成的CIC模型来完成滤波操作,可以用dsp.SignalSource函数,该函数用于生成信号源,函数的第一个参数是原信号,第二个参数是每帧的点数。(至于为什么要分成不同的帧本人也不是很清楚,只是按照Mathworks提供的例程来完成)
下面的例子中,我们要完成的是对输入的1024点的信号进行抽4的CIC滤波,在dsp.SignalSource中我们指定一帧为64个点,因此总共有1024/64=16帧,在抽4之后,每帧的长度就是64/4=16;所以定义的y是一个16x16的矩阵,其中每一行为一帧的滤波结果,将16行拼接起来,就是输出的1024/4=256点的结果。
dsp.CICDecimator和fdesign.decimator的滤波的方法是一样的:
F0 = Fs/40;N = 1024;t = (0:N-1)'; % length: 1024x = fi(sin(2*pi*F0/Fs*t),true,16,15);SamplesPerFrame = 64;src = dsp.SignalSource(x,SamplesPerFrame);y = zeros(N/SamplesPerFrame,SamplesPerFrame/R);for ii = 1:length(x)/SamplesPerFramey(ii,:) = hcic(src());%y(ii,:) = CICDecim(src());endyy = y.';figure;plot(yy(:))
滤波后的时域结果如下所示:
image-20210207155034229
matlab 级联cic,Matlab中CIC滤波器的应用相关推荐
- cic matlab,Matlab中CIC滤波器的应用
CIC滤波器基本原理 CIC(积分梳状级联)滤波器是工程上经常用的滤波器,因为CIC滤波器不需要乘法,CIC滤波器往往在级联抽取滤波器的第一级和级联插值滤波器的最后一级.这一节我们以CIC抽取滤波器为 ...
- 用matlab画多普勒加宽线性函数,MTD雷达中多普勒滤波器组的设计与实现
合肥工业大学理学院电子科学与技术2006届毕业论文集 目 录 中文摘要1 英文摘要2 1 引言3 1.1 研究背景及意义3 1.2 国内外研究现状4 1.3 本设计的指导思想和主要工作4 2 动目标检 ...
- 【MATLAB教程案例1】通信系统中成形滤波器原理的MATLAB设计实现
FPGA教程目录 MATLAB教程目录 -------------------------------------------------------------------------------- ...
- MATLAB中FIR滤波器的时延溢出问题详解:线性相位对信号造成的时延溢出及其消除方法,以及fir1等函数的使用
1.问题由来: 前段时间在对用MATLAB处理试验数据时,需要对多路信号进行滤波后做同步,在这个过程中使用MATLAB中自带的 fir1 函数以及Filter Designer工具箱(FDATool) ...
- matlab中SG滤波器的说明
matlab中SG滤波器的说明 项目需要最近在进行光谱预处理方面研究,看了很多论文,发现其中涉及到的方法也就那么几种,比如SG filter.MSC.SNV等等.今天对SG filter 进行了学习, ...
- cic matlab,matlab生成verilog 语言的CIC interpolation filter
matlab生成verilog 语言的CIC interpolation filter 今天在测试filter偶然发现这个就试了下,还不错. clc clear cic_decim = mfilt.c ...
- fir抽取滤波器matlab,Xilinx FPGA之抽取FIR滤波器的实现
过去半年有几位客户请我帮助他们设计和实现数字下变频器所用的降采样(即"抽取")滤波器,这种滤波器在软件无线电与数据采集类应用中都很常见.这项工作即便对于经验丰富的设计师也不是一件小 ...
- 用MATLAB绘制国债NSS模型,Matlab在数字信号处理中的运用.ppt
<Matlab在数字信号处理中的运用.ppt>由会员分享,可在线阅读,更多相关<Matlab在数字信号处理中的运用.ppt(68页珍藏版)>请在装配图网上搜索. 1.第七讲 M ...
- 【 MATLAB 】使用 MATLAB 求由差分方程表示的滤波器的响应的两种方法
例题: 一个3阶低通滤波器由下面差分方程描述: y(n) = 0.0181 x(n) + 0.0543 x(n-1) + 0.0543 x(n-2) + 0.0181 x(n-3) + 1.76 y( ...
最新文章
- 用漫画了解 Linux 内核到底长啥样!
- poj 1265 Area(pick 定理)
- 电脑运行慢的原因汇总
- centos7.3修改mysql密码_Centos7.3下mysql5.7.18安装并修改初始密码的方法
- ceph部署出现错误及解决
- dp入门 专题记录 2017-7-26
- 【笔记】c++中opencv的使用
- pyqt5讲解11:自定义信号和槽
- 业务id转密文短链的一种实现思路
- 每日一笑 | 史上最贴心的骗子
- Centos7 使用Docker MySQL部署_01
- 数据3分钟丨Gartner宣布明年12大战略性技术趋势;PolarDB-X正式开源;OceanBase 3.2发布...
- 图片添加对话气泡app_微信公众号、对话框、看一看功能又更新!儿童版微信也要来了?!...
- 药品管理系统php,数据库 药品管理系统+源代码
- Pta——敲笨钟,制定位置输出字符串
- 计算机运用基础2020四川传媒学院考题,2020年计算机基础考试题EY[含答案](15页)-原创力文档...
- 一周疫情对行业影响观察 | 对香港IPO冲击大于内地,纺织服装新增企业同比下滑超五成...
- 2019 前端年终总结(干货满满)
- 梦幻西游手游:工坊进阶考试题目攻略—考古、乐艺篇
- UCK Network 为开发者提供全方位孵化,打造区块链爆款应用
热门文章
- linux alpine 中telnet转移至busybox-extras
- ifconfig命令实例
- python中@修饰符用法
- 鸟哥的Linux私房菜(基础学习,服务器架设)
- java浮点数化为整数_如何在JavaScript中将浮点数转换为整数?
- oracle日志版本不同,Oracle重做日志文件版本不一致问题处理
- TCP校验和的设计与实现
- thonny python ide_学习用 Thonny 写代码:一个面向初学者的Python IDE
- android java服务器文件传输_java – 使用FTPS将文件从android传输到服务器
- java之父求职_Java求职实战之继承和多态