[问答]

怎样去设计一种CIC抽取滤波器并对其进行MATLAB仿真呢

1967 滤波器 CIC MATLAB仿真

问答对人有帮助，内容完整，我也想知道答案 0 为什么要设计一种理想滤波器呢？ CIC抽取滤波器是由哪些部分组成的？怎样去设计一种CIC抽取滤波器并对其进行matlab 仿真呢？ 0
2021-11-19 07:29:44　　评论淘帖0 邀请回答您可以邀请以下用户，快速回答问题 × heks 该类别下有 50 个回答。邀请回答 hgimtk 该类别下有 44 个回答。邀请回答新星之火12138 该类别下有 42 个回答。邀请回答 wang21cj 该类别下有 41 个回答。邀请回答 chm5 该类别下有 40 个回答。邀请回答 hjfjsdgfjdsf 该类别下有 37 个回答。邀请回答 fdjslkjd 该类别下有 35 个回答。邀请回答 uvysdfydad 该类别下有 35 个回答。邀请回答 werywer 该类别下有 35 个回答。邀请回答 h1654155957.9520 该类别下有 35 个回答。邀请回答 jenny042 该类别下有 34 个回答。邀请回答 dfzvzs 该类别下有 34 个回答。邀请回答维生素B2 该类别下有 34 个回答。邀请回答凤毛麟角该类别下有 34 个回答。邀请回答刘洋该类别下有 33 个回答。邀请回答 meihuacg 该类别下有 33 个回答。邀请回答 hfgdzc 该类别下有 33 个回答。邀请回答 ggfvxv 该类别下有 33 个回答。邀请回答 kpj3026 该类别下有 32 个回答。邀请回答江左盟该类别下有 32 个回答。邀请回答举报 carey123 相关推荐 • 怎样去设计一种基于加速度计与气压计的三阶卡尔曼滤波器呢 1545 • 怎样使用MATLAB去设计一种FIR滤波器呢 918 • 基于DSP Builder的CIC梳状滤波器该怎么设计？ 1468 • 怎样去计算IIR滤波器的系数呢 1802 • 如何使用matlab的FDAtool工具去设计一种FIR滤波器呢 1590 • 基于FPGA的DDC滤波器该怎么设计？ 2459 • 怎样利用MATLAB去设计IIR滤波器？ 1161 • 怎样去设计一种三阶环路滤波器？求过程 1710 • 如何去实现一种模糊PID并对其进行MATLAB仿真呢 971 • 怎么用XC2V1000型FPGA实现FIR抽取滤波器的设计 1310 1个回答

答案对人有帮助，有参考价值 0 　　CIC抽取滤波器MATLAB仿真和FPGA实现（1）　　设计理想滤波器目标：1、滤波器在有效频段内纹波满足设计要求。2、抽取或内插处理后在有效频段内不产生混叠。3、滤波器实现简单，需要资源较少。这个文章主要讲CIC理论基础，下个文章讲FPGA实现。　　级联积分梳状滤波器又称CIC。这是多速率信号处理中一种结构简单的滤波器，只需要加法器和寄存器即可实现，可以灵活的设置抽取因子和插值因子，并且CIC是一种基于零点和极点相消的FIR滤波器。　　CIC是由积分器和梳状滤波器组成。积分器可表示为（太多公式了只能截图发出来）　　　　根据上面的分析可以知道积分器对直流有着无限的的增益，在应用时这是一个不稳定的系统。　　梳状滤波器可表示为　　　　单级CIC在不同延时长度时冲击响应的频谱特性MATLAB代码如下：　　M=2; %滤波器长度　　b=ones（1，M）; 　　delta=［1，zeros（1，1023）］; 　　s=filter（b，1，delta）; %求取滤波器冲激响应　　Spec=20log10（abs（fft（s）））;%求取滤波器频谱特性　　Spec2=Spec-max（Spec）; 　　f=0:length（Spec）-1; 　　f=2f/（length（Spec）-1）; %对频率轴进行归一化处理　　M=5; %滤波器长度　　b=ones（1，M）; 　　delta=［1，zeros（1，1023）］; 　　s=filter（b，1，delta）; %求取滤波器冲激响应　　Spec=20log10（abs（fft（s）））;%求取滤波器频谱特性　　Spec5=Spec-max（Spec）; 　　M=7; %滤波器长度　　b=ones（1，M）; 　　delta=［1，zeros（1，1023）］; 　　s=filter（b，1，delta）; %求取滤波器冲激响应　　Spec=20log10（abs（fft（s）））;%求取滤波器频谱特性　　Spec7=Spec-max（Spec）; 　　M=8; %滤波器长度　　b=ones（1，M）; 　　delta=［1，zeros（1，1023）］; 　　s=filter（b，1，delta）; %求取滤波器冲激响应　　Spec=20log10（abs（fft（s）））;%求取滤波器频谱特性　　Spec8=Spec-max（Spec）; 　　%plot（f，Spec2，f，Spec5，f，Spec7，f，Spec8）;axis（［0 1 -50 0］）;grid; 　　%xlabel（‘归一化频率’）;ylabel（‘幅度（dB）’）; 　　%legend（‘Ｍ=２’，‘Ｍ=5’，‘Ｍ=7’，‘Ｍ=8’）; 　　%grid; 　　plot（f，Spec2，’-’，f，Spec5，’。’，f，Spec8，’–’）;axis（［0 1 -50 0］）; 　　xlabel（‘归一化频率’）;ylabel（‘幅度（dB）’）; 　　legend（‘Ｍ=2’，‘Ｍ=5’，‘Ｍ=8’）; 　　grid; 　　幅频响应如下图：　　　　　　理论计算与上图仿真结果相对应。　　然而在实际应用中第一级抑制比太低，可以通过级联CIC的方式得到想要的抑制比。在RM一定的情况下由于级联后传递函数是乘机的关系，因此抑制比是累加的关系，如四级CIC级联后第一旁瓣抑制比在53.84dB。但是增加滤波器级联数目虽然可以解决旁瓣电平衰减减小问题，但是会带来通带衰减较大的问题，并且在实现时会消耗更多的资源。下面是5级CIC级联不同RM的MATLAB代码和频谱特性图　　M=2; %滤波器长度　　b=ones（1，M）; 　　delta=［1，zeros（1，1023）］; 　　s1=filter（b，1，delta）; %求取滤波器冲激响应　　s2=filter（b，1，s1）; %求取滤波器冲激响应　　s3=filter（b，1，s2）; %求取滤波器冲激响应　　s4=filter（b，1，s3）; %求取滤波器冲激响应　　s=filter（b，1，s4）; %求取滤波器冲激响应　　Spec=20log10（abs（fft（s）））;%求取滤波器频谱特性　　Spec2=Spec-max（Spec）; 　　f=0:length（Spec）-1; 　　f=2f/（length（Spec）-1）; %对频率轴进行归一化处理　　M=5; %滤波器长度　　b=ones（1，M）; 　　delta=［1，zeros（1，1023）］; 　　s1=filter（b，1，delta）; %求取滤波器冲激响应　　s2=filter（b，1，s1）; %求取滤波器冲激响应　　s3=filter（b，1，s2）; %求取滤波器冲激响应　　s4=filter（b，1，s3）; %求取滤波器冲激响应　　s=filter（b，1，s4）; %求取滤波器冲激响应　　Spec=20log10（abs（fft（s）））;%求取滤波器频谱特性　　Spec5=Spec-max（Spec）; 　　M=7; %滤波器长度　　b=ones（1，M）; 　　delta=［1，zeros（1，1023）］; 　　s1=filter（b，1，delta）; %求取滤波器冲激响应　　s2=filter（b，1，s1）; %求取滤波器冲激响应　　s3=filter（b，1，s2）; %求取滤波器冲激响应　　s4=filter（b，1，s3）; %求取滤波器冲激响应　　s=filter（b，1，s4）; %求取滤波器冲激响应　　Spec=20log10（abs（fft（s）））;%求取滤波器频谱特性　　Spec7=Spec-max（Spec）; 　　M=8; %滤波器长度　　b=ones（1，M）; 　　delta=［1，zeros（1，1023）］; 　　s1=filter（b，1，delta）; %求取滤波器冲激响应　　s2=filter（b，1，s1）; %求取滤波器冲激响应　　s3=filter（b，1，s2）; %求取滤波器冲激响应　　s4=filter（b，1，s3）; %求取滤波器冲激响应　　s=filter（b，1，s4）; %求取滤波器冲激响应　　Spec=20log10（abs（fft（s）））;%求取滤波器频谱特性　　Spec8=Spec-max（Spec）; 　　%plot（f，Spec2，f，Spec5，f，Spec7，f，Spec8）;axis（［0 1 -50 0］）;grid; 　　%xlabel（‘归一化频率’）;ylabel（‘幅度（dB）’）; 　　%legend（‘Ｍ=２’，‘Ｍ=5’，‘Ｍ=7’，‘Ｍ=8’）; 　　%grid; 　　plot（f，Spec2，’-’，f，Spec5，’。’，f，Spec8，’–’）;axis（［0 　　1 -200 0］）; 　　xlabel（‘归一化频率’）;ylabel（‘幅度（dB）’）; 　　legend（‘Ｍ=2’，‘Ｍ=5’，‘Ｍ=8’）; 　　grid; 　　　　从上图可知第一旁瓣衰减67.3dB，但多级CIC主瓣内电平值下降也是较快，这就说明了在相同的通带内滤波器的通带衰减也明显增加。针对通带内衰减问题，可以通过在自定义的通带内来通过比较多级级联的衰减来说明问题，这里设置归一化通带为0.2，级联数目为D，进行MATLAB仿真，代码如下：　　M=5; %滤波器长度　　b=ones（1，M）; 　　delta=［1，zeros（1，1023）］; 　　s1=filter（b，1，delta）; %求取滤波器冲激响应　　s2=filter（b，1，s1）; %求取滤波器冲激响应　　s3=filter（b，1，s2）; %求取滤波器冲激响应　　s4=filter（b，1，s3）; %求取滤波器冲激响应　　s5=filter（b，1，s4）; %求取滤波器冲激响应　　Spec1=20log10（abs（fft（s1）））;%求取滤波器频谱特性　　Spec1=Spec1-max（Spec1）; 　　Spec2=20log10（abs（fft（s2）））;%求取滤波器频谱特性　　Spec2=Spec2-max（Spec2）; 　　Spec3=20log10（abs（fft（s3）））;%求取滤波器频谱特性　　Spec3=Spec3-max（Spec3）; 　　Spec4=20log10（abs（fft（s4）））;%求取滤波器频谱特性　　Spec4=Spec4-max（Spec4）; 　　Spec5=20log10（abs（fft（s5）））;%求取滤波器频谱特性　　Spec5=Spec5-max（Spec5）; 　　f=0:length（Spec1）-1; 　　f=2f/（length（Spec1）-1）; %对频率轴进行归一化处理　　plot（f，Spec1，’-’，f，Spec2，’。’，f，Spec3，’–’，f，Spec5，’-。’）;axis（［0 0.2 -25 0］）; 　　xlabel（‘归一化频率’）;ylabel（‘幅度（dB）’）; 　　legend（‘D=1’，‘D=2’，‘D=3’，‘D=5’）; 　　grid; 　　　　从图中可以看出随着级联数码的增加通带内衰减越大。　　因此在实际应用中要根据输入信号和输出信号的特性来进行使用和级联。需要注意的是CIC抽取并不会造成信号幅度的损失，但经过FFT后频谱幅度与所取的点数成正比，因此抽取后信号的幅度将变为原始信号频谱幅度发1/M（M是抽取系数）。抽取滤波器其实就是低通滤波器，因此低通滤波器的特性CIC抽取滤波器依然具有。

2021-11-19 11:27:22 评论举报张一珠