CIC滤波器可以先对有用信号进行滤波
DDC 由数控整荡器，数字混频器和低通滤波器组成，原理上是输入信号与本地振荡信号混频，然后由低通滤波器滤除高频分量；数字下变频的主要功能包括三个方面：第一是变频，数字混频器将数字中频信号和数控振荡器（Numerical Control Oscillator — NCO ）产生的正 交本振信号相乘
，生成 I/Q 两路混频信号，将感兴趣的信号下变频至零中频；第二是低通滤波，滤除带外信号，提取有用信号；第三是采样速率转换，降低采样速率，大抽取因子范围提供了可设计成宽带或窄带数字信道的能力；
CIC 滤波器可以先对有用信号进行滤波，再抽取；
CIC 抽取滤波器由N 级积分器，抽取器，N 级梳状滤波器三部分组成；N 级积分器工作在Fs 下，每级积分器都是一个反馈系数为1的单级点IIR 滤波器，其传递函数为： CIC 滤波器的梳妆部分工作在较低的频率Fs/ D.，由N 级梳状滤波器组成，每级微分延迟M 个样本；其单级梳状滤波器的传递函数为：
DM C z H --=1，
单级CIC 积分梳状滤波器的传递函数为： 这是单级CIC 的实现方式：
由上式可知，H （z ）有DMN 个零点(M 决定抽取滤波器频率响应中零点个数）和N 个极点，由积分器引人的N 个位于z ＝1处的极点被梳状滤波器的同样位于z=1处的N 个零点抵消； 其单级CIC 频率响应为：
()()()jw jw jw I C H e H e H e ==
sin(/2)sin(/2)wDM w =1()()22
wDM w
DM Sa Sa -??
其中x x x Sa /)sin()(=为抽样函数，且1)0(=Sa ，所以CIC 滤波器在0=ω处的幅
度最大值为DM ，即：DM e H j =)(0；在1...2,1,0,2-==
DM k k DM
w π
处为零；可知当抽取倍数确定后，M 决定CIC 滤波器的零点位置，影响着幅频特性。 N 级CIC 的实现框图：
图3、CIC 抽取滤波器结构
()()()jw jw jw I C H e H e H e ==
sin(/2)sin(/2)wDM w =1()()22
wDM w
DM Sa Sa -??
由其频率响应图可以看出其主瓣电平最大为DM ，旁瓣电平为
21 1.51
()
sin(3/2)/sin(3/2)sin(3/2)
j DM
A H e DM DM ωπ
ωπππ==?
==
，
主瓣与旁瓣的差值 （用dB 数表示）为：
可计算出主瓣与旁瓣的差值约为，意味着阻带衰减很差，单级级联时旁瓣电平很大，为降低旁瓣电平，增加阻带衰减采用级联的方式，N 级频率响应为：
)2()2()()2/sin()2/sin()(ωωωωωN N N N
j N Sa DM Sa DM DM e H -??=??
????=， 可得到N 级CIC 的旁瓣抑制
假如级联5级，此时可以达到的旁瓣衰减；
分析一下发现在N 级联CIC 多出了N DM 这个处理增益，因此分析一下尽量减少带内容差(通带波纹），即，在通带内，幅度应尽量平缓； 由一般数字滤波器的指标，
下面分析一下其频率响应所涉及的一些设计指标： w1 w2 设在w1处抽取的信号带宽很窄，为无混叠信号的带宽,能很好的对窄带信号进行滤波，去除掉高频信号噪声；且在122w DM
w -=
π
处衰减值足够大，则在其信号带宽内，信号在w1到w2所引起的，给CIC 滤波器带来的混叠就可以忽略，计算此时阻带衰减：
)
2/sin()
2/sin(lg 20)()(lg 2022012w DM w DM e H e H A jw j ==·
引入带宽比例因子b=B/（fs/DM ）, B 为抽取信号的带宽，D 为抽取因子，M 为延时因子；fs 为输入端采样率，则w1=b*2pi/DM ；带入可化简得：b A lg 201-≈； （假设b=;即fs=100MHz ，D=20，信号带宽为50khz,此时衰减为40dB);可见单级的CIC 滤波器的无混叠信号带宽内的阻带衰减能达到40dB;；并不怎么大，适用于较粗略的滤波，适合放在第一级抽取；如果采用级联的方式可以加大无混叠信号带宽；但是满足的通带不够窄；
在w1处幅度不能下降太多，通带内幅值容差不能太大，否则会引起高频失真；设该带内容
(0)20lg
()
(0)
20lg
()
a p a p a s a s H j H j H j H j αα=Ω=Ω通带最大衰减阻带最小衰减
差为s δ，则,)()(lg 201
0jw j s e H e H =δ将w1带入可简化得)
sin(lg 20b b
s ππδ≈，当N 级时，其带内容差也会增大；由上面分析可知，阻带衰减和带内容差，只与带宽比例因子b 有关，
D
f B
b s /=
,分析可知，在信号带宽一定的前提下，应尽可能采用小的抽取因子，或增大输入采样率；故一般把它放在抽取系统的第一级，所以在配置CIC 时，信号带宽B ，采样率Fs ，抽取因子D ，延时因子M,需要结合考虑； 总结：
CIC 滤波器本身的性能由M ，D ，N 决定
1、M 决定零点的位置，M=1,2时，比较下面两个图；
（M=1,N=1,D=1) (M=3,N=1,D=1)
发现M=1时通带主瓣衰减较小，而延时因子M=3，衰减很大，但相对带宽减小；M 一般不会取大于2的值；
2、D 决定抽取后信号的采样频率，和M 一起决定主瓣和旁瓣的宽度；
3、N
假设
样率为32MHz 的信号，即抽取倍数D=100=25*4；先将正弦信号加噪声，然后将其与10M 信号混至零频；然后用CIC 进行25倍抽取，HB 滤波器进行4倍抽取，32
1
25/3200400/===
M k D f B b s ； 设计带宽比例因子b=1/32；通带带内容差衰减可忽略；则CIC 滤波器可根据要求通带衰减最大不得超过3dB ，阻带衰减不得低于60dB 。由以上公式及上表查出，可级联6级，差分延迟为1，可完成25倍抽取；半带滤波器完成级联2级的2倍抽取，完成4倍抽取；然后经过FIF 滤波器进行整形；下面是经CIC 滤波器的仿真：
窗函数设计低通滤波器 电信课设
XXXX大学 课程设计报告 学生：xxx 学号：xxx 专业班级：电子信息工程 课程名称：数字信号处理课程设计 学年学期20XX——20XX 学年第X学期指导教师：xxx 2014年6月
课程设计成绩评定表
目录 1. 窗函数设计低通滤波器 1.1设计目的 (1) 1.2设计原理推导与计算 (1) 1.3设计容与要求 (2) 1.4设计源程序与运行结果 (3) 1.5思考题 (10) 2. 用哈明窗设计FIR带通数字滤波器 2.1设计要求 (14) 2.2设计原理和分析 (14) 2.3详细设计 (15) 2.4调试分析及运行结果 (15) 2.5心得体会 (17) 参考文献 (17)
1.窗函数设计低通滤波器 1.1设计目的 1. 熟悉设计线性相位数字滤波器的一般步骤。 2. 掌握用窗函数法设计FIR 数字滤波器的原理和方法。 3. 熟悉各种窗函数的作用以及各种窗函数对滤波器特性的影响。 4. 学会根据指标要求选择合适的窗函数。 1.2设计原理推导与计算 如果所希望的滤波器的理想的频率响应函数为() ωj d e H ，则其对应的单位脉冲响应为 ()() ωπ ωωπ π d e e H n h j j d d ?- = 21 （4.1） 窗函数设计法的基本原理是设计设计低通FIR 数字滤波器时，一般以理想低通滤波特性为逼近函数() ωj e H ，即 ()?????≤<≤=-π ωωωωωα ω c c j j d ，，e e H 0,其中21-=N α ()() ()[]() a n a n d e e d e e H n h c j j j j d d c c --= = = ??- -- πωωπ ωπ ωαωω ωαω π π ω sin 21 21 用有限长单位脉冲响应序列()n h 逼近()n h d 。由于()n h d 往往是无限长序列，而且是非因果的，所以用窗函数()n ω将()n h d 截断，并进行加权处理，得到： ()()()n n h n h d ω= （4.2） ()n h 就作为实际设计的FIR 数字滤波器的单位脉冲响应序列，其频率响应函 数() ωj e H 为 ()()n j N n j e n h e H ωω ∑-==1 （4.3） 式中，N 为所选窗函数()n ω的长度。 用窗函数法设计的滤波器性能取决于窗函数()n ω的类型及窗口长度N 的取
信号与系统试题附答案
信科0801《信号与系统》复习参考练习题一、单项选择题：
14、已知连续时间信号,) 2(100)2(50sin )(--=t t t f 则信号t t f 410cos ·)(所占有的频带宽度为（） A ．400rad ／s B 。200 rad ／s C 。100 rad ／s D 。50 rad ／s
f如下图（a）所示，其反转右移的信号f1(t) 是（） 15、已知信号)(t f如下图所示，其表达式是（） 16、已知信号)(1t A、ε（t）＋2ε(t－2)－ε(t－3) B、ε(t－1)＋ε(t－2)－2ε(t－3) C、ε(t)＋ε(t－2)－ε(t－3) D、ε(t－1)＋ε(t－2)－ε(t－3) 17、如图所示：f（t）为原始信号，f1(t)为变换信号，则f1(t)的表达式是（） A、f(－t+1) B、f(t+1) C、f(－2t+1) D、f(－t/2+1)
18、若系统的冲激响应为h(t),输入信号为f(t),系统的零状态响应是（ ） 19。信号)2(4sin 3)2(4cos 2)(++-=t t t f π π 与冲激函数)2(-t δ之积为（ ） A 、2 B 、2)2(-t δ C 、3)2(-t δ D 、5)2(-t δ ，则该系统是（）＞－系统的系统函数．已知2]Re[,6 51)(LTI 202s s s s s H +++= A 、因果不稳定系统 B 、非因果稳定系统 C 、因果稳定系统 D 、非因果不稳定系统 21、线性时不变系统的冲激响应曲线如图所示，该系统微分方程的特征根是（ ） A 、常数 B 、 实数 C 、复数 D 、实数+复数 22、线性时不变系统零状态响应曲线如图所示，则系统的输入应当是（ ） A 、阶跃信号 B 、正弦信号 C 、冲激信号 D 、斜升信号
数字信号处理-低通滤波器设计实验
实验报告 课程名称：数字信号处理 实验名称：低通滤波器设计实验 院（系）： 专业班级： 姓名： 学号： 指导教师： 一、实验目的： 掌握IIR数字低通滤波器的设计方法。 二、实验原理： 2.1设计巴特沃斯IIR滤波器 在MATLAB下，设计巴特沃斯IIR滤波器可使用butter 函数。 Butter函数可设计低通、高通、带通和带阻的数字和模拟IIR滤波器，其特性为使通带内的幅度响应最大限度地平坦，但同时损失截止频率处的下降斜度。在期望通带平滑的情况下，可使用butter函数。butter函数的用法为：
[b,a]=butter(n,Wn)其中n代表滤波器阶数，W n代表滤波器的截止频率，这两个参数可使用buttord函数来确定。buttord函数可在给定滤波器性能的情况下，求出巴特沃斯滤波器的最小阶数n，同时给出对应的截止频率Wn。buttord函数的用法为：[n,Wn]= buttord(Wp,Ws,Rp,Rs)其中Wp和Ws分别是通带和阻带的拐角频率（截止频率），其取值范围为0至1之间。当其值为1时代表采样频率的一半。Rp和Rs分别是通带和阻带区的波纹系数。 2.2契比雪夫I型IIR滤波器。 在MATLAB下可使用cheby1函数设计出契比雪夫I 型IIR滤波器。 cheby1函数可设计低通、高通、带通和带阻契比雪夫I 型滤IIR波器，其通带内为等波纹，阻带内为单调。契比雪夫I型的下降斜度比II型大，但其代价是通带内波纹较大。cheby1函数的用法为：[b,a]=cheby1(n,Rp,Wn,/ftype/)在使用cheby1函数设计IIR滤波器之前，可使用cheblord 函数求出滤波器阶数n和截止频率Wn。cheblord函数可在给定滤波器性能的情况下，选择契比雪夫I型滤波器的最小阶和截止频率Wn。cheblord函数的用法为： [n,Wn]=cheblord(Wp,Ws,Rp,Rs)其中Wp和Ws分别是通带和阻带的拐角频率（截止频率），其取值范围为0至1之间。当其值为1时代表采样频率的一半。Rp和Rs分别是通带和阻带区的波纹系数。 三、实验要求： 利用Matlab设计一个数字低通滤波器，指标要求如下：
数字信号处理滤波器
1.设计物理可实现的低通滤波器 设计思路：因为要设计FIR有限脉冲响应滤波器，通常的理想滤波器的单位脉冲响应h是无限长的，所以需要通过窗来截断它，从而变成可实现的低通滤波器。程序如下： clc;clear all; omga_d=pi/5; omga=0:pi/30:pi; for N=3:4:51; w1= window(@blackman,N); w2 = window(@hamming,N); w3= window(@kaiser,N,2.5); w4= window(@hann,N); w5 = window(@rectwin,N); M=floor(N/2); subplot(311);plot(-M:M,[w1,w2,w3,w4,w5]); axis([-M M 0 1]); legend('Blackman','Hamming','kaiser','hann','rectwin'); n=1:M; hd=sin(n*omga_d)./(n*omga_d)*omga_d/pi; hd=[fliplr(hd),1/omga_d,hd]; h_d1=hd.*w1';h_d2=hd.*w2';h_d3=hd.*w3';h_d4=hd.*w4';h_d5=hd.*w5'; m=1:M; H_d1=2*cos(omga'*m)*h_d1(M+2:N)'+h_d1(M+1); H_d2=2*cos(omga'*m)*h_d2(M+2:N)'+h_d2(M+1); H_d3=2*cos(omga'*m)*h_d3(M+2:N)'+h_d3(M+1); H_d4=2*cos(omga'*m)*h_d4(M+2:N)'+h_d4(M+1); H_d5=2*cos(omga'*m)*h_d5(M+2:N)'+h_d5(M+1); subplot(312);plot(omga,[H_d1,H_d2,H_d3,H_d4,H_d5]); legend('Blackman','Hamming','kaiser','hann','rectwin'); subplot(313);plot(abs([fft(h_d1);fft(h_d2);fft(h_d3);fft(h_d4);fft(h_ d5)])'); pause(); end 程序分析： 整个对称窗的长度为N，然而为了在MATLAB中看到窗函数在负值时的形状需将N变为它的一半，即为2M+1个长度。窗长设置为从3开始以4为间隔一直跳动51。则长度相同的不同窗函数在时域[-M,M]的形状如第一个图所示。 对窗函数进行傅里叶变换时，将零点跳过去先构造一个一半的理想滤波器的脉冲响应hd，再将零点位置求导得出的数赋值进去。将生成的hd左右颠倒形成了一个理想的滤波器的脉冲响应。将构造的理想滤波器的脉冲响应依次与之前定义的窗函数相乘，相乘出来的为列向量，用转置将其变成行向量，形成的h_d就是非理想的低通滤波器的脉冲响应序列。因为h_d为对称奇数长度序列，它的DTFT 可以是二倍的离散余弦变化，而零点的位置则直接带入求出，两者相加则是H_d。则第二个图表示的是五个矩阵向量在频域的变化，而第三个图表示的是五个非理想低通滤波器的傅里叶变换，图三FFT给出的结果永远是对称的，因为它显示
fir低通滤波器设计(完整版)
电子科技大学信息与软件工程学院学院标准实验报告 （实验）课程名称数字信号处理 电子科技大学教务处制表
电 子 科 技 大 学 实 验 报 告 学生姓名： 学 号： 指导教师： 实验地点： 实验时间：14-18 一、实验室名称：计算机学院机房 二、实验项目名称：fir 低通滤波器的设计 三、实验学时： 四、实验原理： 1. FIR 滤波器 FIR 滤波器是指在有限范围内系统的单位脉冲响应h[k]仅有非零值的滤波器。M 阶FIR 滤波器的系统函数H(z)为 ()[]M k k H z h k z -==∑ 其中H(z)是k z -的M 阶多项式，在有限的z 平面内H(z)有M 个零点，在z 平面原点z=0有M 个极点. FIR 滤波器的频率响应 ()j H e Ω 为 0 ()[]M j jk k H e h k e Ω -Ω ==∑ 它的另外一种表示方法为 () ()()j j j H e H e e φΩΩΩ=
其中 () j H e Ω和()φΩ分别为系统的幅度响应和相位响应。 若系统的相位响应()φΩ满足下面的条件 ()φαΩ=-Ω 即系统的群延迟是一个与Ω没有关系的常数α，称为系统H(z)具有严格线性相位。由于严格线性相位条件在数学层面上处理起来较为困难，因此在FIR 滤波器设计中一般使用广义线性相位。 如果一个离散系统的频率响应 ()j H e Ω 可以表示为 ()()()j j H e A e αβΩ-Ω+=Ω 其中α和β是与Ω无关联的常数，()A Ω是可正可负的实函数，则称系统是广义线性相位的。 如果M 阶FIR 滤波器的单位脉冲响应h[k]是实数，则可以证明系统是线性相位的充要条件为 [][]h k h M k =±- 当h[k]满足h[k]=h[M-k],称h[k]偶对称。当h[k]满足h[k]=-h[M-k],称h[k]奇对称。按阶数h[k]又可分为M 奇数和M 偶数，所以线性相位的FIR 滤波器可以有四种类型。 2. 窗函数法设计FIR 滤波器 窗函数设计法又称为傅里叶级数法。这种方法首先给出()j d H e Ω， ()j d H e Ω 表示要逼近的理想滤波器的频率响应，则由IDTFT 可得出滤波器的单位脉冲响应为 1 []()2j jk d d h k H e e d π π π ΩΩ-= Ω ? 由于是理想滤波器，故 []d h k 是无限长序列。但是我们所要设计的FIR 滤波 器，其h[k]是有限长的。为了能用FIR 滤波器近似理想滤波器，需将理想滤波器的无线长单位脉冲响应 []d h k 分别从左右进行截断。 当截断后的单位脉冲响应 []d h k 不是因果系统的时候，可将其右移从而获得因果的FIR 滤波器。
北京邮电大学数字信号处理习题库选择题附加答案重点
13.下列关于冲激响应不变法描述错误的是 ( C A.S 平面的每一个单极点 s=sk 变换到 Z 平面上 z= e skT 处的单极点 B.如果模拟滤波器是因果稳定的，则其数字滤波器也是因果稳定的 C.Ha(s和 H(z的部分分式的系数是相同的 D.S 平面极点与Z 平面极点都有 z= e s kT 的对应关系 14．下面关于 IIR 滤波器设计说法正确的是( C A. 双线性变换法的优点是数字频率和模拟频率成线性关系 B. 冲激响应不变法无频率混叠现象 C. 冲激响应不变法不适合设计高通滤波器 D. 双线性变换法只适合设计低通、带通滤波器 15.以下关于用双线性变换法设计 IIR 滤波器的论述中正确的是( B 。 A.数字频率与模拟频率之间呈线性关系 B.总是将稳定的模拟滤波器映射为一个稳定的数字滤波器 C.使用的变换是 s 平面到 z 平面的多值映射 D.不宜用来设计高通和带阻滤波器 16.以下对双线性变换的描述中不正确的是 ( D 。 A.双线性变换是一种非线性变换 B.双线性变换可以用来进行数字频率与模拟频率间的变换C.双线性变换把 s 平面的左半平面单值映射到 z 平面的单位圆内 D.以上说法都不对17．以下对双线性变换的描述中正确的是 ( B 。 A．双线性变换是一种线性变换B．双线性变换可以用来进行数字频率与模拟频率间的变换 C．双线性变换是一种分段线性变换 D．以上说法都不对 18．双线性变换法的最重要优点是：；主要缺点是 A 。 A. 无频率混叠现象；模拟域频率与数字域频率间为非线性关系 B. 无频率混叠现象；二次转换造成较大幅度失真 C. 无频率失真；模拟域频率与数字域频率间为非线性关系 D. 无频率失真；二次转换造成较大幅度失真 19．利用模拟滤波器设计法设计 IIR 数字滤波器的方法是先设计满足相应指标的模拟滤波器，再按 某种方法将模拟滤波器转换成数字滤波器。双线性变换法是一种二次变换方法，即它 C 。 A. 通过付氏变换和 Z 变换二次变换实现 B. 通过指标变换和频谱变换二次变换实现 C. 通过二次变换，使得变换后 S 平面与 Z 平面间为一种单值映射关系 D. 通过模拟频率变换和数字频率变换二次变换实现 20．下列对 IIR 滤波器特点的论述中错误的是( C 。 A．系统的单位冲激响应 h(n是无限长的 B.结构必是递归型的C.肯定是稳定的 D.系统函数 H(z在有限 z 平面（0<|z|<∞）上有极点 21．在数字信号处理中通常定义的数字频率ω是归一化频率，归一化因子为 C 。 A．采样周期B. 模拟采样角频率 C. 模拟采样频率 D. 任意频率 22．信号数字频谱与模拟频谱间的一个显著区别在于数字频谱具有 A 。 A．周期性 B. 更大的精确度 C. 更好的稳
常见的信号处理滤波方法
低通滤波：又叫一阶惯性滤波，或一阶低通滤波。是使用软件编程实现普通硬件RC 低通滤波器的功能。 适用范围：单个信号，有高频干扰信号。 一阶低通滤波的算法公式为： Y(n)X(n)(1)Y(n 1)αα=+-- 式中： α是滤波系数；X(n)是本次采样值；Y(n 1)-是上次滤波输出值；Y(n)是本次滤波输出值。 滤波效果1： 红色线是滤波前数据（matlab 中生成的正弦波加高斯白噪声信号） 黄色线是滤波后结果。 滤波效果2：
matlab中函数，相当于一阶滤波，蓝色是原始数据（GPS采集到的x（北）方向数据，单位m），红色是滤波结果。 一阶滤波算法的不足： 一阶滤波无法完美地兼顾灵敏度和平稳度。有时，我们只能寻找一个平衡，在可接受的灵敏度范围内取得尽可能好的平稳度。
互补滤波：适用于两种传感器进行融合的场合。必须是一种传感器高频特性好（动态响应好但有累积误差，比如陀螺仪。），另一传感器低频特性好（动态响应差但是没有累积误差，比如加速度计）。他们在频域上互补，所以进行互补滤波融合可以提高测量精度和系统动态性能。 应用：陀螺仪数据和加速度计数据的融合。 互补滤波的算法公式为： 1122Y(n)X (n)(X (n)Y(n 1))αα+=+-- 式中：1α和2α是滤波系数；1X (n)和2X (n)是本次采样值；Y(n 1)-是上次滤 波输出值；Y(n)是本次滤波输出值。 滤波效果 （测试数据）： 蓝色是陀螺仪 信号，红色是加 速度计信号，黄 色是滤波后的 角度。
. 互补滤波实际效果： .
卡尔曼滤波：卡尔曼滤波器是一个“optimal recursive data processing algorithm （最优化自回归数据处理算法）”。对于解决很大部分的问题，它是最优，效率最高甚至是最有用的。他的广泛应用已经超过30年，包括机器人导航，控制，传感器数据融合甚至在军事方面的雷达系统以及导弹追踪等等。近来更被应用于计算机图像处理，例如头脸识别，图像分割，图像边缘检测。 首先，用于测量的系统必须是线性的。 (k)(k 1)(k)(k)X AX BU w =-++ (k)(k)(k)Z HX v =+ (k)X 是系统k 时刻的状态，(k)U 是系统k 时刻的控制量。(k)Z 是系统k 时 刻的测量值。A 和B 为系统参数，(k)w 和(k)v 分别表示过程和测量的噪声，H 是测量系统参数。 在进行卡尔曼滤波时： 首先进行先验预测： (k 1|k)(k |k)(k)(k)X AX BU w +=++ 计算先验预测方差： '(k 1|k)(k |k)(k)P AP A Q +=+ 计算增益矩阵： (k 1)(k 1|k)'/((k 1|k)'(k 1))Kg P H HP H R +=++++ 后验估计值： (k 1|k 1)(k 1|k)(k 1)(Z(k 1)(k 1|k))X X Kg HX ++=++++-+ 后验预测方差： (k 1|k 1)(1(k 1))(k 1|k)P Kg H P ++=-++ 其中，(k)Q 是系统过程激励噪声协方差，(k)R 是测量噪声协方差。 举例说明： （下文中加粗的是专有名词，需要理解） 预测小车的位置和速度的例子（博客+自己理解）:
巴特沃斯数字低通滤波器要点说明
目录 1.题目........................................................ .................................. .2 2.要求........................................................ (2) 3.设计原理........................................................ . (2) 3.1 数字滤波器基本概念......................................................... (2) 3.2 数字滤波器工作原理......................................................... (2) 3.3 巴特沃斯滤波器设计原理 (2) 3.4脉冲响应不法......................................................... .. (4) 3.5实验所用MATLAB函数说
明 (5) 4.设计思路........................................................ .. (6) 5、实验内容........................................................ . (6) 5.1实验程序......................................................... . (6) 5.2实验结果分析......................................................... . (10) 6.心得体会........................................................ . (10) 7.参考文献........................................................ . (10) 一、题目：巴特沃斯数字低通滤波器 二、要求：利用脉冲响应不变法设计巴特沃斯数字低通滤波器，通带截止频率100HZ,采样频率1000HZ，通带最大衰减为0.5HZ，阻带最小衰减为10HZ，画出幅频、相频相应相应曲线。并假设一个信号x(t)=sin(2*pi*f1*t)+sin(2*pi*f2*t),其中f1=50HZ,f2=200HZ。用此信号验证
08级数字信号处理试卷A及参考答案1
2008 ~2009《数字信号处理》考试试卷（A ） 一、 填空题（本题满分30分，共含4道小题，每空2分） 1. 两个有限长序列x 1(n)，0≤n ≤33和x 2(n)，0≤n ≤36，做线性卷积后结果的长度 是 ，若对这两个序列做64点圆周卷积，则圆周卷积结果中n= 至 为线性卷积结果。 2. DFT 是利用nk N W 的 、 和 三个固有特性来实现FFT 快速运算的。 3. IIR 数字滤波器设计指标一般由 、 、 和 等四项组成。 4. FIR 数字滤波器有 和 两种设计方法，其结构 有 、 和 等多种结构。 二、 判断题（本题满分16分，共含8道小题，每小题2分，正确打√，错误打×） 1. 相同的Z 变换表达式一定对应相同的时间序列。（ ） 2. Chirp-Z 变换的频率采样点数M 可以不等于时域采样点数N 。（ ） 3. 按频率抽取基2 FFT 首先将序列x(n)分成奇数序列和偶数序列。（ ） 4. 冲激响应不变法不适于设计数字带阻滤波器。（ ） 5. 双线性变换法的模拟角频率Ω与数字角频率ω成线性关系。（ ） 6. 巴特沃思滤波器的幅度特性必在一个频带中（通带或阻带）具有等波纹特性。（ ） 7. 只有FIR 滤波器才能做到线性相位，对于IIR 滤波器做不到线性相位。（ ） 8. 在只要求相同的幅频特性时，用IIR 滤波器实现其阶数一定低于FIR 阶数。（ ） 三、 综合题（本题满分18分，每小问6分） 若x (n)= {3，2，1，2，1，2 }，0≤n≤5， 1) 求序列x(n)的6点DFT ，X (k)=？ 2) 若)()]([)(26k X W n g DFT k G k ==，试确定6点序列g(n)=?
CIC滤波器的原理与设计
CIC 的冲击响应{ 1,010,()n D h n ≤≤-= 其他 ，D 为CIC 滤波器的阶数（即抽取因子）， Z 变换后 1 1()1 D z H z z ---=-， 当积分梳状滤波器的阶数不等于抽取器的抽取倍数时，令N=DM(N 为滤波器的 阶数，D 为抽取倍数) 则积分梳状滤波器的传递函数为：)1(11 )(1 DM z z z H ----= M 是梳状滤波器中的延时因子，故称M 为差分延时因子； 其频率总响应为12()()()jw jw jw H e H e H e == sin(/2)sin(/2)wDM w =1()()22 wDM w DM Sa Sa -?? x x x Sa /)sin()(=为抽样函数，且1)0(=Sa ，所以CIC 滤波器在0=ω处的幅度值 为N ，即：DM e H j =)(0； 一般数字滤波器的指标： ()20lg ()()20lg () a p a p a s a s H j H j H j H j ααΩ=ΩΩ=Ω通带最大衰减阻带最小衰减
即： CIC 幅频特性响应曲线图 由其频率响应函数可以看出其主瓣电平最大为D ，旁瓣电平为 21.51 () sin(3/2)/sin(3/2)sin(3/2) j DM H e DM DM ωπωπππ=? == ， 旁瓣与主瓣的差值 （用dB 数表示）为： dB A DM s 46.132 3lg 20lg 201===π α 可计算出旁瓣与主瓣的差值约为13.46，意味着阻带衰减很差，单级级联时旁瓣 电平很大，为降低旁瓣电平，增加阻带衰减采用级联的方式，N 级频率响应为： )2()2()()2/sin()2/sin()(ωωωωωQ Q Q Q j Q Sa DM Sa DM DM e H -??=?? ????=， 可得到N 级CIC 的旁瓣抑制 dB Q Q A DM Q Q s )46.13(2 3lg 20)lg( 201?=?==π α 分析一下发现在Q 级联时多出了Q DM 这个处理增益，因此分析一下尽量减少带容差(通带衰减），即，在通带，幅度应尽量平缓；下面就它的幅平响应曲线来分 析： 00()20lg () ()20lg () p s j a p jw a j a s jw a H e H e H e H e αα==
低通滤波器
实验名称：FIR 低通滤波器的DSP 实现 一、实验目的 1、掌握用窗函数法设计FIR 数字滤波器的原理和方法。 2、熟悉线性相位FIR 数字滤波器的特性。 3、了解各种窗函数对数字滤波器的特性的影响。 二、实验设备 PC 兼容机一台（操作系统为Windows XP ），安装Code Composer Studio 2.2.1软件和MATLAB 6.5.1(含SIMULINK 工具包)软件。 三、实验内容 FIR 低通滤波器的DSP 实现 1、实验要求：用窗口法设计线性相位的FIR 低通滤波器，截止频率为1kHz ，采样速率为8kHz 。窗口大小N=11，分别加矩形窗和海明窗。检验结果并计算峰值的位置和过渡带宽度。 2、对设计要求的理解 （1）要设计的滤波器为理想低通滤波器，便于FIR 低通滤波器的实现。 （2）FIR 滤波器满足线性相位应具有以下两个条件：一是)1,,1,0()(-=N n n h 为 实数；二是h(n)满足以2 1 -= N n 的偶对称或奇对称，即)1()(n N h n h --±=。 故使用窗函数设计滤波器时，所加的窗都以原点对称，设计的理想滤波器的单位脉冲响应d(k)都以原点为偶对称或奇对称。同时保证了d(k)加窗平移后得到的h(n)以M=(N-1)/2对称，也就是保证了设计出的FIR 滤波器具有线性相位。 （3）低通滤波器的系数采用MATLAB 软件仿真工具产生，并把仿真产生的系数导出成头文件，运用到CCS 程序中。课题一中要求采用两种窗函数设计滤波器，两者仅仅是窗函数不同，相应的滤波器系数不同。运用MATLAB 产生两个窗函数对应的系数文件。同时对两种窗函数滤波效果进行对比。 （4）峰值和过滤带的宽度通过理论计算后，再和实际的信号波形进行对比。 3、窗函数法设计FIR 滤波器的思路 首先从窗口大小N 中计算出M ，其中2 1 -=N M ；其次是利用离散时间傅里叶 反变换，从)(ωd 中计算出滤波器系数d(k)；最后考虑到滤波器的因果性，把d(k)延迟M 个单位得到因果的滤波器系数h(n)，其中1,,1,0,)()(-=-=N n M n d n h 。 四、实验原理
(完整word版)数字信号处理题库(附答案)
数字信号处理复习题 一、选择题 1、某系统)(),()()(n g n x n g n y =有界，则该系统（ A ）。 A.因果稳定 B.非因果稳定 C.因果不稳定 D. 非因果不稳定 2、一个离散系统（ D ）。 A.若因果必稳定 B.若稳定必因果 C.因果与稳定有关 D.因果与稳定无关 3、某系统),()(n nx n y =则该系统（ A ）。 A.线性时变 B. 线性非时变 C. 非线性非时变 D. 非线性时变 4.因果稳定系统的系统函数)(z H 的收敛域是（ D ）。 A.9.0 计算机仿真技术实验指导书 河南科技大学电子信息工程学院 二〇〇八年二月 计算机仿真技术实验指导书 MATLAB是一种交互式的以矩阵为基本数据结构的系统。在生成矩阵对象时，不要求明确的维数说明。所谓交互式，是指MATLAB的草稿纸编程环境。 与C语言或FORTRON语言作科学数值计算的程序设计相比较，利用MATLAB可节省大量的编程时间。 本实验指导书主要讨论四个实验。 实验一信号与系统的时域分析以及信号合成与分解 1. 实验目的 (1) 连续时间信号的向量表示法和符号运算表示法，典型离散信号表示； (2) 连续信号和离散信号的时域运算与时域变换； (3) 连续系统和离散系统的卷积，以及冲激响应、阶跃响应、单位响应、零状态响应； (4) 周期信号的傅立叶级数分解与综合（以周期方波为例）； 2. 实验原理与方法 (1) 信号在MATLAB中的表示方法 MATLAB用两种方法来表示连续信号，一种是用向量的方法来表示信号，另一种则是符号运算的方法来表示信号。用适当的MATLAB语句表示出信号后，就可以利用MATLAB的绘图命令绘制出直观的信号时域波形。 向量表示法表示信号的方法是：MATLAB用一个向量表示连续信号的时间范围，另一个向量表示连续信号在该时间范围内的对应样值。如下列代码p=0.001; t=-pi:p:pi; f=1+cos(t); plot(t,f) title('f(t)=1+cos(t)') xlabel('t') axis([-pi,pi,-0.2,2.4]) 执行后即可绘制连续信号1+cos(t)的时域波形。 借助于符号运算以及符号绘图函数ezplot，也可以绘制连续信号时域波形。如下列代码 syms t f=sym('1+cos(t)') %定义符号表达式 ezplot(f,[-pi,pi]) %绘制符号表达式波形 set(gcf,'color','w') %设置当前图形背景颜色为白色 执行后即可绘制连续信号1+cos(t)的时域波形。 与连续信号的表示相似，在MATLAB中，离散信号也需要用两个向量来表示，其中一个向量表示离散信号的时间范围，另一个向量表示该离散信号在该时间范围内的对应样值。但与连续信号表示有所不同的是，表示离散信号时间范围向量的元素必须为整数。如下列代码 n=[-3,-2,-1,0,1,2,3]; x=[-3,2,-1,3,1,-2,1]; stem(n,x,'filled') set(gcf,'color','w') title('x(n)') xlabel('n') 执行后即可绘制离散信号x(n)={ -3,2,-1,3,1,-2,1}的时域波形。 ↑ n=0 (2) 连续信号和离散信号的时域运算与时域变换 对连续信号而言，其基本时域变换有反褶、平移、尺度变换、倒相。 利用MATLAB的符号运算功能以及符号绘图函数ezplot，可以直观的观察和分析连续信号的时域运算与时域变换。如下列代码 syms t; f=sym('(t+1)*(heaviside(t+1)-heaviside(t))'); f=f+sym('(heaviside(t)-heaviside(t-1))'); %定义信号符号表达式 ezplot(f,[-3,3]) %绘制信号波形 axis([-3,3,-1.2,1.2]) set(gcf,'color','w') 实验报告 课程名称：信号分析与处理指导老师：项基成绩：__________________ 实验名称：________滤波器_____实验类型：___研究型________同组学生姓名：__________ 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 一、 实验目的和要求（必填） 1、了解有源滤波器的种类、基本结构、工作原理及其特性。 2、学会测量有源滤波器的幅频特性。 二、 实验内容和原理（必填） 有源滤波器具有体积小、性能好、调整方便等优点，在信号处理方面得到了广泛的应用。 通常高阶的有源滤波器都可由一阶和二阶的滤波器串联而成，其中一阶滤波器只需一只电阻和一只电容构成一级RC 无源网络即可。本实验研究二阶RC 有源滤波器的有关问题。 1．二阶低通有源滤波器 二阶低通有源滤波器的实验电路如图2－7－1（a ）所示。图中将1C 接地端改接到输出端是为了改善 10=ωω附近的滤波器性能。因为在10 <ωω 且接近1的范围内，o u 和i u 相位差小于 90，1C 起正反馈作用， 因而有利于提高这段范围内的输出幅度，而在频带外即10 ??ωω 时，o u 和i u 基本相同，1C 起促进带外衰减 的作用。 当R R R ==21时，该滤波器电路的传递函数为 2 12 122 122002200121 )(C C R s RC s C C R s Q s K s H ++=++=ωωω 截止频率为 2 101C C R = ω 品质因数为 2 1 21C C Q = 通带增益为10=K 。 该电路的优点是改变电阻R 即可改变截止角频率而不影响品质因数Q ，因此，调整时应先调1C 或2C ，使Q 满足要求，然后通过调节电阻R 将0ω调准确。 实验四 基于窗函数的FIR DF 的设计 提示： 1. Matlab 中提供了很多常用的窗函数，其中一些窗函数的调用形式为： 矩形窗：w=boxcar(N) 三角形窗：w=bartlett(N) 汉宁窗：w=hanning(N) 哈明窗：w=hamming(N) 布莱克曼窗：w=blackman(N) 其中，输入参数N 表示窗口的长度，返回的变量w 是一个长度为N 的列向量，表示窗函数在这N 点的取值。 2. b=fir1(N,Wc,'ftype',Window) fir1函数用来设计FIR 滤波器。其中N 为滤波器的阶数；Wc 是截止频率，其取值在0～1之间，它是以π为基准频率的标称值，设计低通和高通滤波器时，Wc 是标量，设计带通和带阻滤波器时，Wc 是1×2的向量；设计低通和带通滤波器时，无需 'ftype'，当ftype=high 时，设计高通滤波器，当ftype=stop 时，设计带阻滤波器；Window 表示设计滤波器所采用的窗函数类型，Window 的长度为N+1，若Window 缺省，则fir1默认使用哈明窗；b 对应设计好的滤波器的系数h(n)，即单位冲激响应，h(n)的长度为N+1。 需注意)(n h 的长度与滤波器的阶数间的关系。FIR 滤波器的系统函数可表示为： ∑-=-=1 )()(N n n z n h z H )(n h 的长度为N ，而滤波器的阶数为1-N 阶。 3. 求数字滤波器的频率响应 h=freqz(b,a,w) 其中，b 和a 分别为系统函数)(z H 的分子多项式和分母多项式的系数。对于FIR 滤波器，此处的b 即为h(n)，a 可看作1。 实验题目： 1. 分别用矩形窗和哈明窗设计FIR 低通滤波器，设窗宽11=N ，截止频率rad c πω 2.0=，要求绘出两种窗函数设计的滤波器幅频曲线，并进行比较。 滤波器组框架理论及其在图信号处理中的应用 摘要：传统滤波器组框架理论通常用来处理低维规则结构数据,如时间信号、空 间信号和时空信号等。随着现代科技高速发展,高维非规则化数据信息大量涌现, 如社交网络、能源网络、交通运输网络、神经元网络等。如何对高维图结构数据 进行处理成为一个备受关注且亟待解决的问题。借助代数图论和谱图理论,图信号 处理成为近年来兴起的研究方向,用来处理高维加权图上的信号。众多学者从各自 角度出发,将传统滤波器组框架理论推广到图滤波器组框架中,取得了一系列成果。 关键词：滤波器组；框架理论；图信号；图滤波器 引言：滤波器组框架理论是应用数学、信号处理、图像处理和数字通信等领 域的重要问题之一,对滤波器组框架的分析和设计问题进行研究有着重要的科学意 义和应用前景。近年来,随着高维非规则化数据信息大量涌现,很多学者开始研究 图信号处理的滤波器组方法。因此对滤波器组框架理论及其在图信号处理中的应 用进行研究。 一、滤波器组框架理论 在各种框架中,实际应用最广泛的是由滤波器组实现的框架。有限维框架、离 散小波框架和离散Gabor框架都属于滤波器组框架。接下来介绍滤波器组基础知识、滤波器组框架理论及应用。 （一）滤波器组基础 滤波器组是一组有着共同输入或共同输出的带通滤波器。典型滤波器组的结 构如下图所示。其中左边部分为分析滤波器组,右边部分为综合滤波器组。分析滤 波器组有一个输入多个输出,其将输入信号分解成不同的子带信号,每个分析滤波 器Hi(z)有不同的频率特性,输入信号x(n)通过M个分析滤波器Hi(z)后,得到M个不 同的子带信号。信号在子带分解后,对每个通道Mi下采样，可降低信号的采样率。下采样后的子带信号可以被编码、处理或者传输。综合滤波器组具有多个输入一 个输出,其将处理后的子带信号通过带通滤波后再组合起来,重构原始信号。为保 证重构信号x?(n)与原信号x(n)具有相同的采样频率,在综合滤波器组前对各子带信 号Mi上采样(Upsampling)。也有论文将下采样称为抽取(Decimation),将上采样称 为内插(Interpolation),两者实际并无区别,本文统一称为下采样、上采样。 M通道滤波器组： 将每个通道的下采样因子Mi相同的滤波器组称为均匀滤波器组；将下采样因子不同的滤波器组称为非均匀滤波器组.将下采样因子和通道数相同的滤波器组称 为临界采样滤波器组；将下采样因子小于通道数的滤波器组称为过采样滤波器组。如果滤波器组由理想滤波器构成,没有混叠产生,则可以完全重构原始信号。由于 理想滤波器是不可实现的,为了消除混叠,需要选择合适的Hi(z)和Fi(z),使得 x?(n)=x(n?m),这样的滤波器组称为完全重构滤波器组。多采样率信号处理的核心 是信号采样率的转换和滤波器组。信号的上/下采样是多采样率信号处理的基本操作。多相(Polyphase)结构是滤波器组的一种基本表示方法。 （二）滤波器组框架 框架理论最先由Duffin等在研究非谐波Fourier序列时创立的,小波框架和Gabor框架是应用最广泛的两类框架。二十世纪八九十年代,与小波理论并行发展 的滤波器组分析和设计方法使得小波的物理实现成为现实,此后小波在信号处理、 数据压缩与编码等领域得到了飞速发展和巨大应用。目前滤波器组框架理论在采 实验六《理想滤波器、原型模拟滤波器和窗函数的特性》1、实验内容 1、计算下列理想数字滤波器的单位冲激响应，并画出其频率响应和单位冲激响应，观察单位冲激响应波形的对称特性 1）理想低通滤波器，截止频率0.3π，群延时10 2）理想高通滤波器，截止频率0.65 π，群延时20 3）理想带通滤波器，下、上截止频率0.35 π、0.7 π，群延时15 2、画出下列原型模拟滤波器的幅度响应特性、相位响应特性和衰减特性，频率范围0—12000Hz （调用freqs），观察它们在通带、阻带、过渡带宽度、相位特性等方面的特点。 1）巴特沃斯低通滤波器，截止频率5000π，阶数5，调用butter 2）切比雪夫I型低通滤波器，截止频率5000 π，阶数5，通带波纹0.5dB，调用cheby1 3）切比雪夫II型低通滤波器，截止频率5000 π，阶数5，阻带衰减50dB，调用cheby2 4）椭圆滤波器，截止频率5000 π，阶数5，通带波纹0.5dB，阻带衰减50dB，调用ellip 3、编写程序画出下列窗函数的时域图形和频域特性（幅度dB表示和相位），与矩形窗函数 相比，观察它们在阻带最小衰减、主瓣宽带等方面的特点。 1）矩形窗，长度40 2）三角窗，长度40 3）升余弦窗，长度40 4）Blackman，长度40 2、编程原理、思路和公式 1、首先写出理想低通、高通、带通滤波器的频率响应，画出其频谱图，然后根据计算 得到的各滤波器的脉冲响应，写出它们的Matlab表达形式，画出脉冲响应图形。 三者的程序类似，只是在具体的频率响应和脉冲响应的形式上有所差别。 低通单位脉冲响应： 1,|| () 0,|| ja j c LP c e H e ω ω ωω ωωπ - ?? ? ?? ≤ = <≤ 对应的单位脉冲响应为： 1 ()() 2 sin[()] 11 2() j j n LP LP ja j n c h n H e e d n a e e d n a πωω π πωω π ω π ω ω ππ - - - = - == - ? ? 高通单位脉冲响应为： sin[()] () () n a n a c h n HP n a ω π = - =- - 对应的单位脉冲响应为： sin[()]sin[()] () () H L n a n a h n BP n a ωω π --- = - 2、以butterworth低通滤波器为例，其余三种只是调用的函数不同而已，原理相同。首 先写出滤波器的阶数、截至频率，然后调用butter函数得到滤波器的系统函数，再调用freqs函数得到0-12000Hz范围内的频谱函数，最后画出幅度响应特性、相位响应特性和衰减特性的图形。 3、首先调用各种窗的Matlab函数生成各自的时域函数并画出时域图形，然后进行fft 变换得到频谱特性，再转化为幅度特性（db表示）和相位特性。各种窗的Matlab 函数是：矩形窗rectwin，三角窗triang，升余弦窗hanning，blackman窗blackman。数字信号处理和滤波器设计
滤波器信号分析与处理实验
信号与系统实验四 答案
滤波器组框架理论及其在图信号处理中的应用
理想滤波器、原型模拟滤波器和窗函数的特性matlab6