数字滤波器

数字梳状滤波器

梳状滤波对于画面质量是非常重要的一个技术，因此我们有必要对其进行详细刨析。

那么具体什么是梳状滤波器呢？这就要从源头（信号源）开始讲起了，一开始，接收视频的Video端子是Composite端子（比如RF射频接口和AV接口），它所能接收的信号叫Composite Video Signal，即混合视频信号（也称复合信号），什么意思呢？因为这个Composite（混合）信号包括了亮度（Luminance，用字母Y表示）和色度/彩度（Chrominace）两方面的信号，视频电路要做的工作就是Y/C进行分离处理，目前的梳状滤波器是在保证图像细节的情况下解决视频信号亮色互窜的唯一方法，其内部有许多按一定频率间隔相同排列的通带和阻带，只让某些特定频率范围的信号通过，因为其特性曲线象梳子一样，故人们称之为梳状滤波器（Comb Filtering）。

梳状滤波器一般由延时、加法器、减法器、带通滤波器组成。对于静止图像，梳状滤波在帧间进行，即三维梳状滤波。对活动图像，梳状滤波在帧内进行，即二维梳状滤波。高档数字电视机采用行延迟的梳状滤波器与带通滤波器级联，构成Y/C分离方案就可获得满意的图像质量。使用梳状滤波器能使图像质量明显提高。解决了色串亮及亮串色造成的干扰光点、干扰花纹；消除了色度正交分量U、V色差信号混迭造成的彩色边缘蠕动；消除了亮、色镶边，消除了高频信号的色彩错误和灰度值表示错误。有一段时期国内很多工厂（为了节省成本）使用模拟的方式实现梳状滤波器，实际上效果很不好，原因有两个，一是延迟器件的带宽很难保证，二是解决行相关性差问题的自适应电路很复杂。而在数字电路里，只要有足够的存储器，就可以保证足够的延迟时间与信号带宽，且复杂的自适应电路很容易集成在芯片中硬件固化。

梳状滤波器原理及发展历史：

梳状滤波器采用频谱间置技术，理论上可以保证亮度和色度的无失真分离。如果我们好好回顾一下梳状滤波器的发展历程，将对其有个清醒的认识。

第一阶段：采用频率分离法将Y/C信号分开。这种方法是利用色度信号以副载波方式传输这一特点（PAL制副载波为4.43MHz，NTSC制副载波为3.58MHz），用选频电路将Y/C 信号分开。

内部由LC带通滤波器和陷波器组成，将视频信号通过一个中心频率（fsc）为色度信号窄带（比如PAL制式4.43MHz频率副载波）带通滤波器，取出色度信号。再将亮度信号经过一个中心频率为色度信号副载波4.43MHZ的色度陷波器，吸收色度信号，从而得到亮度信号。这种方法简单易行，采用元器件少且成本低，所以在早期彩电中应用得比较广泛。

但是，频率分离方法存在着一些严重的问题：在亮度通道中，色度陷波器在吸收色度信号的同时也将该频率范围内亮度信号的部分频率分量抑制掉了，这叫亮度信号的高频分量丢失，从而影响了亮度信号的清晰度——大家想想看，一段音乐哪怕再动听高音丢了那还能听吗？同时，残余的色度信号也可能进入亮度通道而引起串色干扰，通常在屏幕出现彩色测试卡是最后两条频带染色现象。

在色度通道中，L.C色度带通滤波器品质不高，取出色度信号，抑制亮度信号的同时，也把该频率范围内的亮度信号选出来了，高频亮度信号经色度解调器被解调出来（属于多余信号），使得一些细格子或条状区域出现闪烁的彩色干扰。可见利用传统的频率分离方法根本不能将Y、C信号作出彻底分离，必然存在着：亮串色、色串亮”的干扰，使图像质量难以令人满意。事实上这种方法在彩电中的应用，会出现图所示的影响。

左图为普通电视画面，右边为内置梳状滤波电路电视画面

第二阶段：采用梳状滤波器和频谱分离法进行亮色分离。它是根据视频信号频谱交叉的原理及梳状滤波器的梳齿滤波频率传输特性,以频谱分离的方式分离出亮度和色度信号，这种新的分离方法使Y/C信号分离比较干净彻底，从而大幅提高图像清晰度。通常梳状滤波器是由两行延迟线、加法器、减法器等部分组成。

事实上在大屏幕彩电中，又分为NTSC制Y/C分离和PAL制Y/C分离。例如对于NTSC －M制式，我们假设相邻两行的视频信号保持相关性以及延迟线无损耗，Y信号频谱与C 信号频谱以fH/2间隔交替出现（fH表示行频），副载波频率fsc为227.5fH，如果设计一个梳状滤波器电路，使V信号延时一行，再分别与未延时的信号进行加减。延时前后Y 信号相位不变，而C信号相位相反。延时信号与直通信号在加法器中相加后得到Y信号，即（Y＋C）＋（Y－C）＝2Y，在减法器中相减则得到C信号，即（Y＋C）－（Y－C）＝2C。从梳状滤波器幅频特性曲线分析，Y频谱落在加法器特性曲线峰点及减法器特性曲线谷点，所以比较彻底地使亮度信号与色度信号相互分离开来。

PAL梳状滤波器Y/C分离，常称之为二行分离法（也叫二元分离法）。它使用一个NTSC 制信号或两个PAL制信号行存储器与带通滤波器结合使用，组成垂直、水平二元带通滤波器，基本结构原理如上图所示，只是将1H延时线改为2H延时线。（这是因为NTSC制亮度信号、色度信号采用fH/2间置，而PAL制则采用fH/4间置（即副载频fsc＝283.75fH ＋25Hz≈283.75fH），因此PAL梳状滤波器Y/C分离电路要用2H延时线。）由亮度信

号Y和色度信号C组成的复合全电视信号，一方面直接加到加法器与减法器的输入端，称为直通信号，同时经两行延迟时间（2TH）使色度信号反相后得到的延迟信号（Y－C）也加到加法器和减法器的另一输入端。在加法器中直通信号（Y＋C）与延时信号（Y－C）相加得到亮度信号2Y；在减法器中直通信号（Y＋C）与延时信号（Y－C）相减即得到色度信号C，达到亮色分离的目的。

梳状滤波器Y/C分离法的特点可归纳为：由于加法器输出特性可选出亮度信号的高频分量，不会造成高频分量的丢失，并可将视频带宽全部加以利用，从而使图像清晰度大大提高，同时在亮度通道中将色度信号抑制得比较彻底，不致于产生残留色度信号干扰。而减法器输出中较好地抑制掉亮度信号，以最大传输系数选出色度信号C，并用带通滤波器对残留亮度信号作进一步衰减，解决了亮度信号对色度通道的串扰问题，从而提高了图像质量。

左右图分别是标准电视信号测试有无梳状滤波的拍摄画面，看看亮色互窜画面是多么“缤纷多彩”。

上述分析结果是基于信号相关性的假设,可将色度信号与亮度信号较彻底分离而获得较为理想的图象质量。但实际的视频信号并不是这样理想的，即会出现非相关情况，如垂直方向有色度跳变，那么在此处直通信号与延迟信号中的Y、C分量不再相同，加法器与减法器便不能将C或Y分量完全对消，造成Y与C分离不彻底。

如果不愿意看长篇大论的朋友直接看下面的表格：

表四：N 制/PAL 制梳状滤波器的异同

NTSC 制式梳状滤波器 PAL 制式梳状滤波器

延时

1H 2H

直通信号

Y ＋C 延时后反相信号

Y －C （注明：亮度信号Y 不变，仅色度信号C 反相） 进入加法器

（Y ＋C ）＋（Y －C ）＝2Y ，解析出亮度信号 进入减法器

（Y ＋C ）－（Y －C ）＝2C ，解析出色度信号 优点 能够解决部分亮、色干扰问题，提高图像质量

缺点 要求前后两行信号要比较相似，即要求行相关性好。例如蓝色的天空、

黑沉的夜晚等静止画面图像信号，其亮度信号和色度信号可在这种Y/C

分离电路中彻底分离开来；但如果出现垂直方向有色度信号突变的不相

关情况（即前一行信号与后一行信号不同，尤其是动态画面），它会将

其当成相交信号来处理，于是在加法器中不能完全抵消色度信号，造成

Y/C 分离不彻底，仍出现色点干扰。 我们看到，PAL 制梳状滤波器比N 制梳状滤波器要略微复杂一些，不过不管怎么说，既然还有不足之处，那仍然需要做改进，所以梳状滤波器技术继续发展。

第三阶段：（模拟式）动态梳状滤波器

模拟动态梳状滤波器结构如图所示。它由两个延时线（PAL 为2H 延时线，NTSC 为1H 延时线）、三个带通滤波器、垂直相关性检测电路、加法器、减法器等组成。就PAL 制彩色信号而言，要对在动态梳状滤波器中直通信号（Y ＋C ）、延迟两行时间信号（Y －C ）、延迟四行时间信号（Y ＋C ）这三行信号进行垂直方向上的图像相关性检测，产生一个所需的彩色信号。

所谓相关性检测，实际上是检测场与场之间相关性的强弱，一般采用的是场差法或低通场差法，即对两场中各对应像素逐点相减并求和，以该值大小作为图像动态情况的描述。模拟动态梳状滤波器克服了普通梳状滤波器的缺点，改善了活动图像信号Y/C分离效果，从而进一步提高了图像质量。

区别普通梳状滤波器动态梳状滤波器相关性检测无有水平方向梳状滤波有有，效果更好

垂直方向梳状滤波无有

动态画面Y/C分离无有

成本一般设计复杂，成本较高

这是90年底初的技术，当时的大屏幕进口彩电，如29寸的东芝、松下中均采用了此种梳状滤波器。

第四阶段：动态数字式梳状滤波器。

上述的模拟动态梳状滤波器Y/C分离电路，虽然能有效解决活动图像信号的Y/C分离问题，但对单制式信号Y/C分离要用6个调节点，若要适应PAL、NTSC制Y/C信号分离，则需12个调节点，这样就存在调整繁杂的问题，如此多的调节点在生产中难以保证质量，况且批量生产会成为难题。于是，人们研制出更为先进的、精密的数字式动态梳状滤波器。

数字式动态梳状滤波器Y/C分离电路如下图所示。它主要由五块IC封装一起组成厚膜电路，即由A/D变换器CXD1176Q、延迟线CXK1202×2、数字式动态梳状滤波器CXD2011Q 及D/A变换器CXD1177Q组成，与时钟信号发生器配合完成Y/C分离。

新型三行数字化梳状滤波器结构图

它是利用三行彩色信号来完成垂直方向的相关检测，仅提取所需要的彩色信号，从而克服了模拟梳状滤波器的缺点，使图象的水平清晰度从350线一举提高到450线以上！

首先，全电视信号进来后，由模数转换器（ADC）转换成8bit的数字式信号，我们看上图中下半部分是典型的锁相环路，用于产生四倍于色副载波振荡取样频率，作为数字Y/C分离电路的时钟（这部分是一典型的时钟信号发生器），对PAL制为17.73MHZ，对NTSC 制为14.32MHZ。数字化后的视频信号送入动态梳状滤波器，在动态梳状滤波器中进行数字式动态梳状滤波Y/C分离，原理与前述大同小异，只是前、中、后三行视频信号经过色带通滤波处理后再进入逻辑运算电路。此电路中，每相邻的两行信号相减后都可取出色度信号C，而将此C信号再与中间一行视频信号（Y＋C）相加，则抵消了C信号而分离出Y信号。8bit数字亮度信号Y和数字色度信号C分别送入数模变换器（DAC），经其转换后输出模拟Y信号和C信号，分离效果极彻底且无需作任何调整，所以是效果较好，广泛应用的Y/C分离电路。

笔者手头有一份长虹经典2919PK的电路图，里面对梳状电路的原理讲解（信号流程）非常透彻，将来有机会再与各位同好共同探讨。

值得注意的是，在逻辑运算器中，仍然有一个垂直相关电路，用来比较判别三行信号的差别，当三行信号差别较大时，说明图像内容在垂直方向发生了变化，电路即进行运算，其运算过程是：先前两行信号进行运算取最小值，再对后两行信号进行运算并取最小值，然后把两个最小值进行比较取出大者；用同样方法取出此三行信号中相邻两行信号中的最大值，然后把两个最大值再作比较并取小者。最后把上述最小值中的大者和最大值中的小者平均后输出，此平均值应能代表相邻三行信号中的变化趋势，使输出接近于未来的状态。这就是数字式动态梳状滤波器的工作特点，多次比较和平均的结果，使数据非常精确，从而保证了分离效果，所以认为它是一种先进的Y/C信号分离电路。说到这儿，我们不仅想起显卡的AA全屏抗锯齿技术，不和梳状滤波有异曲同工之妙？技术的发展是能智能判断以求最少的代价得到最准确的结果。

小知识：为什么很多电视机宣传单里号称有梳状滤波器（N制），而没有PAL制？

很简单：因为单一制梳状滤波器成本低啊，PAL制梳状滤波虽然不复杂，但是却不能与N制滤波电路共用，因此势必增加成本。我们收看的电视节目一般是PAL制，所以没有PAL制的电视机就像MX440SE 用64bit DDR显存一样，属于低成本的产品，不推荐购买。还有我不理解的是为什么有的厂家吝惜得连这种对画质起至关重要作用的电路都要省掉，要知道一个普通的梳状滤波器仅仅价值3.5美元啊！

存储三个层面，比较并进行分割，东芝飞视采用的技术

目前专用的梳状滤波器已经越来越少了，因为独立数字梳状滤波器成本不菲，而且只能匹配模拟的彩色解码电路，我们知道一个好的电路，从头到脚不能存在瓶颈，否则效果会大打折扣——这就是单独做万金油般全能梳状滤波器不划算的地方，所以最新的数字电路解决方案是将数字化彩色解码芯片和数字梳状滤波器集成，至于做成单芯片还是多芯片具体要看厂商的设计功力了。常见的梳状滤波器有2D和3D之分，2D之中又有两行、三行、五行之分，它们具备不同的自适应效果，成本和价格差距很大。

最后我们总结成如下几点，以方便大家观看：

1、目前最先进的梳状滤波器是3D数字式梳状滤波器，它能够从空间（2D）、时间（第三维方向）将每组画面的亮度及色度信号精确地分离，有效消除影响信号中的杂波、斑点、色彩重叠现象，使画面更加清晰。

2、较好的梳状滤波器必须是双制式梳状滤波器——尤其是对国内消费者来说；

3、数字梳状滤波器有不同的自适应算法，因此各家效果也不同。但视频存贮能力（相当于显卡的显示缓存）越大，视频处理芯片越精密，算法越复杂其Y/C分离效果越好。

梳状滤波器等

级普通两行梳状滤

波器

动态模拟梳状滤

波器

三行动态数字梳

状滤波器（N制）

五行动态数字梳

状滤波器（P制）

N制智能梳状数

字滤波器

3D数字梳状滤

波器

清晰度350线400线450线480线500～520线直观看图辨优劣：

好，理论归理论，实际归实际，那么梳状滤波器对哪些画面更敏感呢？我们来看看下面的电视信号测试图：

竖直线条疏密相间，还有各种同心圆最能看出画面清晰度和亮色串扰问题，不过这种电视标准信号发生碟片可真不好弄。笔者手头只有一张低清晰度版本的，只能说明部分问题。

左右分别为无Y/C分离和3D梳状滤波画面，是不是大不一样？

DPTV-DX的实际效果图，梳状滤波还是很干净的。

这张图是Trident官方提供的资料，梳状滤波效果很明显。

其实呢大商场里比比皆是的实际演示碟很容易说明问题，关键是很多人不了解应该在哪些场景下进行判断：

我们注意上图色块交界处合，垂直线条处的细节是最容易出现问题的地方。

大面积色块交界处，所有有细小竖直或者横向条纹的画面细节、画面上的文本处，都是梳状滤波器大显身手的地方。梳状滤波器将亮色分离得更彻底，使画面更清晰、纯净和细腻完美。

看这种纵横交错的高楼大厦最能看出电视画面的好坏，比较差的（比如普通隔行无三行梳状滤波）电视机在这种地方原形毕露，模糊且闪烁得厉害。

再举一个例子当更能说明问题，看小包中间的竖直条纹是否清晰，好的梳状滤波器能有效去除亮度、色度串扰及色彩斑点。还有一个地方很说明问题，就是有大片郁郁葱葱的绿色植物时的电视画面……我不会骗你的，在这个地方是考验电视机真功夫的地方。

上面是举例说明梳状滤波器对哪些画面有影响，下面是相同画面的对比。

例一：

通过篮子的篾条可以很明显地看出，普通梳状滤波器的水平清晰度不错，但是在垂直方向上不是很好。

数字梳状滤波器更好的进行亮色分离，所以图像的垂直清晰度得到了提高。

而3D数字梳状滤波器更将斜线方向的信号干净地分离，这一点在篮子的提手上表现得很明显。

例二：

有没有梳状滤波器效果大不相同，所以说梳状滤波器不是可有可无的问题，而是没有梳状滤波器就一定是垃圾产品的问题。

各种梳状滤波器的区别：

采用技术观察效果

简单模拟梳状滤波器有限的亮度信号和色度信号分辨率，当观察高频率部份（如多波群信号，AVS frame 50816）时可以看到彩色串扰。即亮色分离不干净，高频成分（画面的细节）模糊。

模拟梳状滤波器（Glass或CCD）上述问题得到改善，但是在竖直彩色块的水平分界线和垂直分界线上出现不应该有的点。（测试信号：SMPTE color

bars，AVS Frame 50815）

数字2D自适应梳状滤波器上述在水平分界线上不应该出现的点被消除。垂直分界线上的点仍然存在，在斜线上仍然可以看到串色（测试信号：AVS Frame 50828）。

数字2D自适应梳状滤波器三菱改善型上述在垂直分界线上不应该出现的点被消除。（测试信号：AVS Frame 50815）

数字3D运动自适应梳状滤波器基本上没有亮度串扰和色度串扰。实际电视画面比较：

无梳状滤波器，亮色窜扰，高频信号丢失

有梳状滤波器，画面整洁干净，决无上述弊端

无梳状滤波器，白色字符上出现不应有的串色，字符边缘粗糙。

有梳状滤波器，白色字符上白色干净，字符边缘清晰。

想来大家看了这么多对比画面当能对梳状滤波（最简单的理解是它能像一把梳子一样把画面不良信号梳干净）有了清晰的概念

FIR数字滤波器设计与使用

实验报告 课程名称：数字信号处理指导老师：刘英成绩：_________________实验名称： FIR数字滤波器设计与使用同组学生姓名：__________ 一、实验目的和要求 设计和应用FIR低通滤波器。掌握FIR数字滤波器的窗函数设计法，了解设计参数（窗型、窗长）的影响。 二、实验内容和步骤 编写MATLAB程序，完成以下工作。 2-1 设计两个FIR低通滤波器，截止频率 C =0.5。 （1）用矩形窗，窗长N=41。得出第一个滤波器的单位抽样响应序列h 1(n)。记下h 1 (n) 的各个抽样值，显示h 1 (n)的图形（用stem(.)）。求出该滤波器的频率响应（的N 个抽样）H 1(k)，显示|H 1 (k)|的图形（用plot(.)）。 （2）用汉明窗，窗长N=41。得出第二个滤波器的单位抽样响应序列h 2(n)。记下h 2 (n) 的各个抽样值，显示h 2(n)的图形。求出滤波器的频率响应H 2 (k)，显示|H 2 (k)|的 图形。 （3）由图形，比较h 1(n)与h 2 (n)的差异，|H 1 (k)|与|H 2 (k)|的差异。 2-2 产生长度为200点、均值为零的随机信号序列x(n)（用rand(1,200)0.5）。显示x(n)。 求出并显示其幅度谱|X(k)|，观察特征。 2-3 滤波 （1）将x(n)作为输入，经过第一个滤波器后的输出序列记为y 1(n)，其幅度谱记为|Y 1 (k)|。 显示|X(k)|与|Y 1 (k)|，讨论滤波前后信号的频谱特征。 （2）将x(n)作为输入，经过第二个滤波器后的输出序列记为y 2(n)，其幅度谱记为|Y 2 (k)|。 比较|Y 1(k)|与|Y 2 (k)|的图形，讨论不同的窗函数设计出的滤波器的滤波效果。 2-4 设计第三个FIR低通滤波器，截止频率 C =0.5。用矩形窗，窗长N=127。用它对x(n)进行滤波。显示输出信号y

常用的8种数字滤波算法

常用的8种数字滤波算法 摘要：分析了采用数字滤波消除随机干扰的优点，详细论述了微机控制系统中常用的8种数字滤波算法，并讨论了各种数字滤波算法的适用范围。 关键词：数字滤波；控制系统；随机干扰；数字滤波算法 1 引言 在微机控制系统的模拟输入信号中，一般均含有各种噪声和干扰，他们来自被测信号源本身、传感器、外界干扰等。为了进行准确测量和控制，必须消除被测信号中的噪声和干扰。噪声有2大类：一类为周期性的，其典型代表为50 Hz 的工频干扰，对于这类信号，采用积分时间等于20 ms整倍数的双积分A/D转换器，可有效地消除其影响；另一类为非周期的不规则随机信号，对于随机干扰，可以用数字滤波方法予以削弱或滤除。所谓数字滤波，就是通过一定的计算或判断程序减少干扰信号在有用信号中的比重，因此他实际上是一个程序滤波。 数字滤波器克服了模拟滤波器的许多不足，他与模拟滤波器相比有以下优点： (1)数字滤波器是用软件实现的，不需要增加硬设备，因而可靠性高、稳定性好，不存在阻抗匹配问题。 (2)模拟滤波器通常是各通道专用，而数字滤波器则可多通道共享，从而降低了成本。 (3)数字滤波器可以对频率很低(如0.01 Hz)的信号进行滤波，而模拟滤波器由于受电容容量的限制，频率不可能太低。 (4)数字滤波器可以根据信号的不同，采用不同的滤波方法或滤波参数，具有灵活、方便、功能强的特点。 2 常用数字滤波算法 数字滤波器是将一组输入数字序列进行一定的运算而转换成另一组输出数字序列的装置。设数字滤波器的输入为X(n)，输出为Y(n)，则输入序列和输出序列之间的关系可用差分方程式表示为： 其中：输入信号X(n)可以是模拟信号经采样和A/D变换后得到的数字序列，也

数字滤波器总结

1数字滤波器的应用领域 在信号处理过程中，所处理的信号往往混有噪音，从接收到的信号中消除或减弱噪音是信号传输和处理中十分重要的问题。根据有用信号和噪音的不同特性，提取有用信号的过程称为滤波，实现滤波功能的系统称为滤波器。在近代电信设备和各类控制系统中，数字滤波器应用极为广泛，这里只列举部分应用最成功的领域。 (1) 语音处理 语音处理是最早应用数字滤波器的领域之一，也是最早推动数字信号处理理论发展的领域之一。该领域主要包括5个方面的内容：第一，语音信号分析。即对语音信号的波形特征、统计特性、模型参数等进行分析计算；第二，语音合成。即利用专用数字硬件或在通用计算机上运行软件来产生语音；第三，语音识别。即用专用硬件或计算机识别人讲的话，或者识别说话的人；第四，语音增强。即从噪音或干扰中提取被掩盖的语音信号。第五，语音编码。主要用于语音数据压缩，目前已经建立了一系列语音编码的国际标准，大量用于通信和音频处理。近年来，这5个方面都取得了不少研究成果，并且，在市场上已出现了一些相关的软件和硬件产品，例如，盲人阅读机、哑人语音合成器、口授打印机、语音应答机，各种会说话的仪器和玩具，以及通信和视听产品大量使用的音频压缩编码技术。 (2) 图像处理 数字滤波技术以成功地应用于静止图像和活动图像的恢复和增强、数据压缩、去噪音和干扰、图像识别以及层析X射线摄影，还成功地应用于雷达、声纳、超声波和红外信号的可见图像成像。 (3) 通信 在现代通信技术领域内，几乎没有一个分支不受到数字滤波技术的影响。信源编码、信道编码、调制、多路复用、数据压缩以及自适应信道均衡等，都广泛地采用数字滤波器，特别是在数字通信、网络通信、图像通信、多媒体通信等应用中，离开了数字滤波器,几乎是寸步难行。其中，被认为是通信技术未来发展方向的软件无线电技术，更是以数字滤波技术为基础。 (4) 电视 数字电视取代模拟电视已是必然趋势。高清晰度电视的普及指日可待，与之配套的视频光盘技术已形成具有巨大市场的产业；可视电话和会议电视产品不断更新换代。视频压缩和音频压缩技术所取得的成就和标准化工作，促成了电视领域产业的蓬勃发展，而数字滤波器及其相关技术是视频压缩和音频压缩技术的重要基础。 (5) 雷达 雷达信号占有的频带非常宽，数据传输速率也非常高，因而压缩数据量和降低数据传输速率是雷达信号数字处理面临的首要问题。高速数字器件的出现促进了雷达信号处理技术的进步。在现代雷达系统中，数字信号处理部分是不可缺少的，因为从信号的产生、滤波、加工到目标参数的估计和目标成像显示都离不开数字滤波技术。雷达信号的数字滤波器是当今十分活跃的研究领域之一。 (6) 声纳

IIR数字滤波器设计及软件实现

实验四：IIR 数字滤波器设计及软件实现 一、实验原理与方法 1、设计IIR 数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法，其基本设计过程是： （1）将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； （2）设计过渡模拟滤波器； （3）将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。 本实验的数字滤波器的MATLAB 实现是指调用MATLAB 信号处理工具箱函数filter 对给定的输入信号x(n)进行滤波，得到滤波后的输出信号y(n ）。 二、实验内容 1、调用信号产生函数mstg 产生由三路抑制载波调幅信号相加构成的复合信号st ，该函数还会自动绘图显示st 的时域波形和幅频特性曲线，如图4.1所示。由图可见，三路信号时域混叠无法在时域分离。但频域是分离的，所以可以通过滤波的方法在频域分离，这就是本实验的目的。 图4.1 三路调幅信号st （即s （t ））的时域波形和幅频特性曲线 2、要求将st 中三路调幅信号分离，通过观察st 的幅频特性曲线，分别确定可以分离st 中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率。要求滤波器的通带最大衰减为0.1dB,阻带最小衰减为60dB 。实验结果如图4.2，程序见附录4.2。 提示：抑制载波单频调幅信号的数学表示式为 0001()cos(2)cos(2)[cos(2())cos(2())]2 c c c s t f t f t f f t f f t ππππ==-++ 其中，cos(2)c f t π称为载波，fc 为载波频率，0cos(2)f t π称为单频调制信号，f0为调制正弦波信号频率，且满足0c f f >。由上式可见，所谓抑制载波单频调幅信号，就是2个正弦信号相乘，它有2个频率成分：和频0c f f +和差频0c f f -，

数字滤波器的MATLAB设计与实现.

数字滤波器的MATLAB设计与实现 数字滤波器的MATLAB设计与实现 类别：电子综合 引言 随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。 1 数字滤波器的设计1.1 数字滤波器设计的基本步骤数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器。IIR滤波器的特征是，具有无限持续时间冲激响应。种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲激响应只能延续一定时间，在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法和Chebyshev逼近法等等。随着MATLAB软件尤其是MATLAB的信号处理工作箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。数字滤波器设计的基本步骤如下：（1）确定指标在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给出幅度和相位响应。幅度指标主要以两种方式给出。第一种是绝对指标。它提供对幅度响应函数的要求，一般应用于FIR滤波器的设计。第二种指标是相对指标。它以分贝值的形式给出要求。在工程实际中，这种指标最受欢迎。对于相位响应指标形式，通常希望系统在通频带中人有线性相位。运用线性相位响应指标进行滤波器设计具有如下优点：①只包含实数算法，不涉及复数运算；②不存在延迟失真，只有固定数量的延迟；③长度为N 的滤波器（阶数为N-1），计算量为N/2数量级。因此，本文中滤波器的设计就以线性相位FIR滤波器的设计为例。（2）逼近确定了技术指标后，就可以建立一个目标的数字滤波器模型。通常采用理想的数字滤波器模型。之后，利用数字滤波器的设计方法，设计出一个实际滤波器模型来逼近给定的目标。（3）性能分析和计算机仿真上两步的结果是得到以差分或系统函数或冲激响应描述的滤波器。根据这个描述就可以分析其频率特性和相位特性，以验证设计结果是否满足指标要求；或者利用计算机仿真实现设计的滤波器，再分析滤波结果来判断。 1.2 滤波器的MATLAB设计（1）MATLAB MATLAB是一套用于科学计算的可视化高性能语言与软件环境。它集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个界面友好的用户环境。它的信号处理工具箱包含了各种经典的和现代的数字信号处理技术，是一个非常优秀的算法研究与辅助设计的工具。在设计数字滤波器时，通常采用MATLAB来进行辅助设计和仿真。（2）FIR滤波器的MATLAB设计下面以设计线性相位FIR滤波器为例介绍具体的设计方法。线性相位FIR滤波器通常采用窗函数法设计。窗函数法设

数字滤波器设计与分析

吉林大学仪器科学与电气工程学院本科生实习报告 实习题目：信号分析和处理 实习时间：2012.09 专业：电气工程及其自动化 所在班级：65100615 学生姓名：王双伟 指导教师：朱凯光田宝凤林婷婷

信号实习报告 一.实验目的 加深对信号系统与信号处理理论的理解，学会信号处理的基本知识和方法，并在基本技能方面得到系统训练；熟悉MA TLAB编程环境，掌握MA TLAB编程基本技能，以及程序调试仿真方法，能够采用MATLAB语言和工具进行信号处理；掌握现代信号分析与处理技术，包括信号频谱分析和数字滤波器（FIR、IIR）设计，学会信号处理系统设计与系统功能检测的基本方法；将理论知识与实际应用结合，提高学生解决实际问题的动手能力，为信号系统与信号处理知识的应用、后续专业学习以及今后从事相关科学研究和实际工作打下坚实基础。二.实验工具 计算机，matlab软件 三.实验内容 设计FIR数字带通滤波器，对于给定函数s=sin(2πx100t)+sin(2πx200t)+sin(2πx400t)，设计带通滤波器滤除100和400赫兹的频率，并画出滤波前后的时频图及滤波器的增益图。 f1=100;f2=200;f3=400; fs=2000; m=(0.3*f1)/(fs/2); M=round(8/m); N=M-1; fc=[0.15,0.3]; b=fir1(N,fc); figure(1) [h,f]=freqz(b,1,1000); plot(f*fs/(2*pi),20*log10(abs(h))) xlabel('频率/赫兹'); ylabel('增益/分贝'); title('滤波器的增益响应'); figure(2) subplot(211) t=0:1/fs:0.5; s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t); plot(t,s); xlabel('时间/秒'); ylabel('幅度'); title('信号滤波前时域图');

数字滤波器的DSP实现

摘要 当前我们正处于数字化时代，数字信号处理技术受到了人们的广泛关注，其理论及算法随着计算机技术和微电子技术的发展得到了飞速的发展，被广泛应用于语音图象处理、数字通信、谱分析、模式识别、自动控制等领域。数字滤波器是数字信号处理中最重要的组成部分之一，几乎出现在所有的数字信号处理系统中。数字滤波器是指完成信号滤波处理的功能，用有限精度算法实现的离散时间线性非时变系统，其输入是一组(由模拟信号取样和量化的)数字量，其输出是经过变换的另一组数字量。相对于模拟滤波器，数字滤波器没有漂移，能够处理低频信号，频率响应特性可做成非常接近于理想的特性，且精度可以达到很高，容易集成等，这些优势决定了数字滤波器的应用越来越广泛。同时DSP(数字信号处理器)的出现和FPGA的迅速发展也促进了数字滤波器的发展，并为数字滤波器的硬件实现提供了更多的选择。 本论文的主要研究了数字滤波器的基本理论及其算法。基于TI公司的数字信号处理器TMS320VC5509设计了一款稳定度高，低功耗的数字滤波器系统，并完成了软硬调试工作。主要工作如下： (1)研究了数字滤波器的基本理论，以及数字滤波器的实现方法。通过学习识字滤波器 的结构、数字滤波器的设计理论，掌握了各种数字滤波器的原理和特性。为实现数字滤波器奠定了理论基础。 (2)研究分析了如何利用MATLAB仿真软件来设计出符合各种要求的数字滤波器。并采用 了相关的函数设计了几款常用的数字滤波器，并得到了滤波器的相关系数，为利用DSP实现数字滤波做好了一些前期的工作。 (3)根据TI公司5000系列数字信号处理器的基本结构和特征，充分利用其片上资源t结 合MATLAB软件的仿真，用软件实现高性能稳定的数字滤波器。 关键字：数字滤波器，DSP，IIR(无限长单位脉冲响应)，FIR(有限长单位脉冲响应)

FIR数字滤波器设计及软件实现

实验五：FIR数字滤波器设计及软件实现 一、实验目的： （1）掌握用窗函数法设计FIR数字滤波器的原理和方法。 （2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。 （3）掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MATLAB函数设计与实现FIR滤波器。 二、实验容及步骤： （1）认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理； （2）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图1所示； 图1 具有加性噪声的信号x(t)及其频谱如图 （3）请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。 （4）根据滤波器指标选择合适的窗函数，计算窗函数的长度N，调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 （4）重复（3），滤波器指标不变，但改用等波纹最佳逼近法，调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。 友情提示： ○1MATLAB函数fir1和fftfilt的功能及其调用格式请查阅本课本；

○ 2采样频率Fs=1000Hz ，采样周期T=1/Fs ； ○ 3根据图10.6.1(b)和实验要求，可选择滤波器指标参数：通带截止频率fp=120Hz ，阻带截至频率fs=150Hz ，换算成数字频率，通带截止频率 p 20.24p f ωπ=T =π，通带最大衰为0.1dB ，阻带截至频率s 20.3s f ωπ=T =π，阻带最小衰为60dB 。] ○ 4实验程序框图如图2所示。 图2 实验程序框图 三、实验程序： 1、信号产生函数xtg 程序清单： %xt=xtg(N) 产生一个长度为N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz %载波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz. function xt=xtg N=1000;Fs=1000;T=1/Fs;Tp=N*T; t=0:T:(N-1)*T; fc=Fs/10;f0=fc/10; %载波频率fc=Fs/10，单频调制信号频率为f0=Fc/10;

数字滤波器的设计及实现

数字滤波器的设计及实现 【一】设计目的 1. 熟悉IIR 数字滤波器和FIR 数字滤波器的设计原理和方法； 2. 学会调用MATLAB 信号处理工具箱中的滤波器设计函数设计各种IIR 和FIR 数字滤波器，学会根据滤波要求确定滤波器指标参数； 3. 掌握用IIR 和FIR 数字滤波器的MA TLAB 实现方法，并能绘制滤波器的幅频特性、相频特性； 4. 通过观察滤波器的输入、输出信号的时域波形及其频谱，建立数字滤波的概念。 【二】设计原理 抑制载波单频调幅信号的数学表达式为 []))(2cos())(2cos(2 1)2cos()2cos()(000t f f t f f t f t f t s c c c ++-==ππππ （2.1） 其中，)2cos(t f c π称为载波，c f 为载波频率，)2cos(0t f π称为单频调制信号，0f 为调制正弦波信号频率，且满足0c f f >。由（2.1）式可见，所谓抑制载波单频调制信号，就是两个正弦信号相乘，它有2个频率成分：和频c f +0f ，差频c f -0f ，这两个频率成分关于载波频率c f 对称。所以，1路抑制载波单频调幅信号的频谱图是关于载波频率c f 对称的两根谱线。 复合信号st 产生函数mstg 清单： function st=mstg %产生信号序列st ，并显示st 的时域波形和频谱 %st=mstg 返回三路调幅信号相加形成的混合信号，长度N=800 N=800; %信号长度N 为800 Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz ，Tp 为采样时间 t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; %第1路调幅信号载波频率fc1=1000Hz fm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100Hz fc2=Fs/20; %第2路调幅信号载波频率fc2=500Hz fm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hz fc3=Fs/40; %第3路调幅信号载波频率fc3=250Hz fm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hz xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第1路调幅信号 xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第2路调幅信号 xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第3路调幅信号 st=xt1+xt2+xt3; %三路信号相加，得到复合信号

FIR数字滤波器设计与软件实现(精)讲解学习

实验二:FIR 数字滤波器设计与软件实现 一、实验指导 1.实验目的 (1掌握用窗函数法设计 FIR 数字滤波器的原理和方法。 (2掌握用等波纹最佳逼近法设计 FIR 数字滤波器的原理和方法。 (3掌握 FIR 滤波器的快速卷积实现原理。 (4学会调用 MA TLAB 函数设计与实现 FIR 滤波器。 2. 实验内容及步骤 (1认真复习第七章中用窗函数法和等波纹最佳逼近法设计 FIR 数字滤波器的原理; (2调用信号产生函数 xtg 产生具有加性噪声的信号 xt ,并自动显示 xt 及其频谱,如图 1所示;

图 1 具有加性噪声的信号 x(t及其频谱如图 (3请设计低通滤波器,从高频噪声中提取 xt 中的单频调幅信号,要求信号幅频失真小于 0.1dB ,将噪声频谱衰减 60dB 。先观察 xt 的频谱,确定滤波器指标参数。 (4根据滤波器指标选择合适的窗函数,计算窗函数的长度 N ,调用 MATLAB 函数 fir1设计一个 FIR 低通滤波器。并编写程序,调用 MATLAB 快速卷积函数 fftfilt 实现对 xt 的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 (5 重复 (3 , 滤波器指标不变, 但改用等波纹最佳逼近法, 调用MA TLAB 函数 remezord 和 remez 设计 FIR 数字滤波器。并比较两种设计方法设计的滤波器阶数。 提示:○ 1MA TLAB 函数 fir1的功能及其调用格式请查阅教材; ○ 2采样频率 Fs=1000Hz,采样周期 T=1/Fs;

○ 3根据图 1(b和实验要求,可选择滤波器指标参数:通带截止频率 fp=120Hz,阻带截 至频率 fs=150Hz, 换算成数字频率, 通带截止频率 p 20.24 p f ωπ =T=π, 通带最大衰为 0.1dB , 阻带截至频率 s 20.3 s f ωπ =T=π,阻带最小衰为 60dB 。 3、实验程序框图如图 2所示,供读者参考。 图 2 实验程序框图 4.信号产生函数 xtg 程序清单(见教材 二、滤波器参数及实验程序清单 1、滤波器参数选取 根据实验指导的提示③选择滤波器指标参数: 通带截止频率 fp=120Hz,阻带截至频率 fs=150Hz。代入采样频率 Fs=1000Hz,换算成 数字频率,通带截止频率 p 20.24 p f

FIR数字滤波器设计与实现

FIR 数字滤波器设计与实现 一.摘要：数字滤波器是一种具有频率选择性的离散线性系统，在信号数字处理中有着广泛的应 用。其中FIR 滤波器是一种常用的滤波器，它在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性，在语音分析、图像处理、雷达监测等对信号相位要求高的领域有着广泛的应用，能实现IIR 滤波器不能实现的许多功能。 二.关键词：FIR 窗函数系统函数MATLAB 三.内容提要： 数字滤波器的功能就是把输入序列通过一定的运算变换成输出序列，因此数字滤波器的结构系 统中就必须包括一定数量和性能的运算器件和运算单元，而运算器件和运算单元的配置必须由数字滤波器的结构特点和性能特点来决定，因此在进行FIR 数字滤波器的设计之前，有必要介绍和总结FIR 数字滤波器的基本结构和相关特性（包括频响曲线（幅度和相位），单位冲激响应等），在介绍完其基本结构和相关特性后，就进行FIR 数字滤波器的设计和实现。 （一）FIR 滤波器的基本结构 在讨论任何一种滤波器时，都要着重分析其系统函数，FIR 滤波器的系统函数为： n N n z n h z H ∑-==1 0)()(。从该系统函数可看出，FIR 滤波器有以下特点： 1)系统的单位冲激响应h(n)在有限个n 值处不为零； 2)系统函数H(z)在|z|>0处收敛，极点全部在z=0处(稳定系统)； 3)结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中(例如频率抽样结构)也包 含有反馈的递归部分。 1.FIR 滤波器实现的基本结构有： 1） 横截型（卷积型、直接型） a.一般FIR 滤波器的横截型(直接型、卷积型)结构： 若给定差分方程为： 。则可以直接由差分方程得出FIR 滤波器结构如 下图所示： 这就是FIR 滤波器的横截型结构，又称直接型或卷积型结构。 b .线性相位FIR 滤波器的横截型结构

FIR数字滤波器设计与软件实现

实验二：FIR数字滤波器设计与软件实现 一、实验指导 1．实验目的 （1）掌握用窗函数法设计FIR数字滤波器的原理和方法。 （2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。（3）掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MATLAB函数设计与实现FIR滤波器。 2．实验内容及步骤 （1）认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理； （2）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图1所示； 图1 具有加性噪声的信号x(t)及其频谱如图 （3）请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。

（4）根据滤波器指标选择合适的窗函数，计算窗函数的长度N，调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。（4）重复（3），滤波器指标不变，但改用等波纹最佳逼近法，调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。 提示：○1MATLAB函数fir1的功能及其调用格式请查阅教材； ○2采样频率Fs=1000Hz，采样周期T=1/Fs； ○3根据图1(b)和实验要求，可选择滤波器指标参数：通带截止频率fp=120Hz，阻带截至频率fs=150Hz，换算成数字频率，通带截止 频率 p 20.24 p f ωπ =T=π，通带最大衰为0.1dB，阻带截至频率 s 20.3 s f ωπ =T=π，阻带最小衰为60dB。 ○4实验程序框图如图2所示，供读者参考。

FIR数字滤波器设计与软件实现

一、实验目的 （1）掌握用窗函数法设计FIR数字滤波器的原理和方法。 （2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。 （3）掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MA TLAB函数设计与实现FIR滤波器。 二、实验内容及步骤 （1）认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理； （2）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图10.5.1所示； 图10.5.1 具有加性噪声的信号x(t)及其频谱如图 （3）请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。 （4）根据滤波器指标选择合适的窗函数，计算窗函数的长度N，调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 （4）重复（3），滤波器指标不变，但改用等波纹最佳逼近法，调用MATLAB函数remezord和remez 设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。 提示：○1MA TLAB函数fir1和fftfilt的功能及其调用格式请查阅本书第7章和第？章； ○2采样频率Fs=1000Hz，采样周期T=1/Fs； ○3根据图10.6.1(b)和实验要求，可选择滤波器指标参数：通带截止频率fp=120Hz，阻带截至频率 fs=150Hz，换算成数字频率，通带截止频率 p 20.24 p f ωπ =T=π，通带最大衰为0.1dB，阻带截至频率 s 20.3 s f ωπ =T=π，阻带最小衰为60dB。] ○4实验程序框图如图10.5.2所示。

经典数字滤波器

数字滤波器是现在电视中常用的电路元件之一。数字滤波器(digital filter)是由数字乘法器、加法器和延时单元组成的一种装置。其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。由于电子计算机技术和大规模集成电路的发展，数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。数字滤波器是一个离散时间系统（按预定的算法，将输入离散时间信号转换为所要求的输出离散时间信号的特定功能装置）。应用数字滤波器处理模拟信号时，首先须对输入模拟信号进行限带、抽样和模数转换。数字滤波器输入信号的抽样率应大于被处理信号带宽的两倍，其频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率即1／2抽样频率点呈镜像对称。为得到模拟信号，数字滤波器处理的输出数字信号须经数模转换、平滑。数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。数字滤波器在语言信号处理、图像信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用。数字滤波器有低通、高通、带通、带阻和全通等类型。它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。应用最广的是线性、时不变数字滤波器. 阶数越高，截止频率等参数越精确，但是电路结构也越复杂。简单说比如你的截止频率是100HZ，你只有2阶的话可能实际的截止平率是95-1000HZ，衰减比较慢，但如果是20阶的话，可能截止

频率就变成了95-105HZ，衰减很快。但是阶数上升，实际电路的结构就会非常的复杂，浪费资源。 先听我慢慢说啊，先说傅里叶变换，然后再说滤波器，就懂了。 周期信号可以用一系列的不同频率不同幅度的正弦信号表示出来，就是傅里叶级数。 而非周期信号亦可以，比如门信号，它的傅氏变换是抽样信号，意思就是，它可以用的一系列不同频率的正弦信号表示，比如有：频率为0.1Hz幅度为2的正弦，频率为0.2Hz幅度为1的正弦，频率为0.25幅度为a的正弦……这些无数个的所谓的“频率为某Hz幅度为某”的正弦波叠加之后，就成了门信号。 从门信号的频谱图可看出：用来表示门信号的一系列频率连续的无数个的正弦波幅度是不同的，甚至有些是0 。尤其频率越高的正弦波，它们的幅度普遍很小，因为这些频率成分是表示细节（门信号的棱角）的。另一方面，低频成分显示的是门信号的轮廓。

实验五IIR数字滤波器设计及软件实现

实验四：IIR数字滤波器设计及软件实现 一、实验内容及步骤 1、调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号 st，三路信号在时域混叠无法在时域分离，但频域是可分离的，所以可以通过滤波的方法在频域分离。 2、要求将st中三路调幅信号分离，通过观察st的幅频特性曲线，分别确定可 以分离st中三路抑制载波单频调幅信号的三个滤波器（低通、高通、带通）的通带截止频率和阻带截止频率。要求滤波器的通带最大衰减为0.1db，阻带最小衰减为60db. 3、编程序调用MATLAB滤波器设计函数ellipord和ellip分别设计三个椭圆滤 波器，并绘图显示其损耗函数曲线。 4、调用滤波器实现函数filter，用三个滤波器分别对信号产生函数mstg产生 的信号st进行滤波，分离出st中的三路不同载波频率的调幅信号yn1、yn2、yn3的，并绘图显示其时域波形，观察分离效果。 二、实验结果显示 原信号图形：

高通滤波器 输出波形

带通滤波器输出波形

低通滤波器输出波形

带阻滤波器输出波形

三、实验结论：由上面所绘图形可知，利用数字滤波器完全可以将时域混叠而频域未混叠的波形分开，达到滤波目的。 四、思考题 （1）请阅读信号产生函数mstg，确定三路调幅信号的载波频率和调制信号频率。答：第一路调幅信号的调制信号频率为100HZ，载波频率为1000HZ；第二路调幅信号的调制信号频率为50HZ，载波频率为500HZ；第三路调幅信号的调整信号频率为25HZ，载波频率为250HZ。 （2）信号产生函数mstg中采样点数N=1600，对st进行N点FFT可以得到6根理想谱线。如果取N=1800，可否得到6根理想谱线？为什么？N=2000呢？请改变函数mstg 中采样点数N的值，观察频谱图验证您的判断是否正确？ 答： 因为信号st是周期序列，谱分析时要求观察时间为整数倍周期。分析可知，st的每个频率成分都是25Hz的整数倍。采样频率Fs=10kHz=25×400Hz，即在25Hz 的正弦波的1个周期中采样400点。所以，当N为400的整数倍时一定为st的整数个周期。因此，采样点数N=800和N=2000时，对st进行N点FFT可以得到6根理想谱线。如果取N=1000，不是400的整数倍，不能得到6根理想谱线。

数字滤波器的一般概念

数字滤波器的一般概念 滤波器可广义地理解为一个信号选择系统。它让某些信号成分通过又阻止或衰减另一些成分。在更多地情况下，被窄义地理解为选频系统，如低通、高通、带通、带阻。频域与时域均衡器也是一种滤波器，通信系统的传输媒介如明线、电缆等从特性看也是滤波器。滤波器如系统一样可分为三类：模拟滤波器、采样滤波器和数字滤波器.模拟滤波器（AF）可以是由RLC构成的无源滤波器，也可以是加上运放的有源滤波器，它们是连续时间系统。采样滤波器（SF）由电阻、电容、电荷转移器件、运放等组成，属于离散时间系统，其幅度是连续的。开关电容滤波器、电荷耦合滤波器军属这类滤波器。数字滤波器（DF）由加法器、乘法器、存储延迟单元、时钟脉冲滤波器及逻辑单元等数字电路构成。它精度高，稳定性好，不存在阻抗匹配问题，可以时分复用，能够完成一些模拟滤波器完成不了的滤波任务。其缺点是需要抽样、量化、编码，以及手时钟频率所限，所能处理的信号最高频率还不够高。另外，由于有限字长效应会造成域设计值的频率偏差、量化和运算噪声及极限环振荡。 本章讨论的是数字滤波器。 5.1.1 数字滤波器的分类 下面从各种不同角度对数字滤波器分类： 1.按冲激响应h(n)的长度分类 分为有限冲激响应（FIR）DF和无限冲激响应（IIR）DF两种。冲 激响应本来是用于模拟系统，指系统对冲激函数δ（t）的响应。 发展到数字滤波器后，工程上仍沿用这个名称，与单位抽样响应和 单位脉冲响应的说法通用。 FFR DF的冲激响应h(n)为有限长序列，其差分方程为 y(n)= (5.1) 系统函数为 H(z)= (5.2) IIR DF 的冲激响应h(n)为无限长序列，其差分方程为

FIR数字滤波器的设计与实现

FIR 滤波器的设计 一.摘 要：数字滤波器是一种具有频率选择性的离散线性系统，在信号数字处理中有 着广泛的应用。其中FIR 滤波器是一种常用的滤波器，它在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性，在语音分析、图像处理、雷达监测等对信号相位要求高的领域有着广泛的应用，能实现IIR 滤波器不能实现的许多功能。 二.关键词：FIR 窗函数 系统函数 MATLAB 三.引言： 数字滤波器的功能就是把输入序列通过一定的运算变换成输出序列，因此数字滤波器的结构系统中就必须包括一定数量和性能的运算器件和运算单元，而运算器件和运算单元的配置必须由数字滤波器的结构特点和性能特点来决定，因此在进行FIR 数字滤波器的设计之前，有必要介绍和总结FIR 数字滤波器的基本结构和相关特性（包括频响曲线（幅度和相位），单位冲激响应等），在介绍完其基本结构和相关特性后，就进行FIR 数字滤波器的设计和实现。 （1）.FIR 滤波器的基本结构 在讨论任何一种滤波器时，都要着重分析其系统函数，FIR 滤波器的系统函数为： n N n z n h z H ∑-== 1 )()( 。从该系统函数可看出，FIR 滤波器有以下特点： 1)系统的单位冲激响应h(n)在有限个n 值处不为零； 2)系统函数H(z)在|z|>0处收敛，极点全部在z=0处(稳定系统)； 3)结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中(例如频率抽样结构)也包含有反馈的递归部分。 1.FIR 滤波器实现的基本结构有： 1） 横截型（卷积型、直接型） a.一般FIR 滤波器的横截型(直接型、卷积型)结构： 若给定差分方程为： 。 则可以直接由差分方程得出FIR 滤波器 结构如下图所示： 这就是FIR 滤波器的横截型结构，又称直接型或卷积型结构。 b .线性相位FIR 滤波器的横截型结构 若h(n)呈现对称特性，即此FIR 滤波器具有线性相位，则可以简化成横截型结构，下面分情况讨论： ①N 为奇数时线性相位FIR 滤波器实现结构如图所示：

FIR数字滤波器设计实验_完整版

班级： 姓名： 学号： FIR 数字滤波器设计实验报告 一、实验目的 1．掌握FIR 数字滤波器的设计方法； 2．熟悉MATLAB 信号处理工具箱的使用； 3．熟悉利用MATLAB 软件进行FIR 数字滤波器设计，以及对所设计的滤波器 进行分析； 4．了解FIR 滤波器可实现严格线性相位的条件和特点； 5．熟悉FIR 数字滤波器窗函数设计法的MATLAB 设计，并了解利用窗函数法 设计FIR 滤波器的优缺点； 6．熟悉FIR 数字滤波器频率采样设计法的MATLAB 设计，并了解利用频率采 样法设计FIR 滤波器的优缺点； 7．熟悉FIR 数字滤波器切比雪夫逼近设计法的MATLAB 设计，并了解利用切 比雪夫逼近法设计FIR 滤波器的优缺点。 二、实验设备及环境 1．硬件：PC 机一台； 2．软件：MATLAB （6.0版以上）软件环境。 三、实验内容及要求 1．实验内容：基于窗函数设计法、频率采样设计法和切比雪夫逼近设计法，利用MATLAB 软件设计满足各自设计要求的FIR 数字低通滤波器，并对采用不同设计法设计的低滤波器进行比较。 2．实验要求： （1）要求利用窗函数设计法和频率采样法分别设计FIR 数字低通滤波 器，滤波器参数要求均为：0.3c w π=。其中，窗函数设计法要求分别利用矩形窗、汉宁窗和布莱克曼窗来设计数字低通滤波器，且 21N ≥,同时要求给出滤波器的幅频特性和对数幅频特性； 频率

采样法要求分别利用采样点数21N =和63N =设计数字低通滤波器，同时要求给出滤波器采样前后的幅频特性，以及脉冲响应及对数幅频特性。 （2）要求利用窗函数设计法和切比雪夫逼近法分别设计FIR 数字低通 滤波器，滤波器参数要求均为： 0.2π, 0.25dB, 0.3π, 50dB p p s s ωαωα==== 其中，窗函数设计法要求利用汉明窗来设计数字低通滤波器，且 66N ≥，同时要求给出滤波器理想脉冲响应和实际脉冲响应，汉 名窗和对数幅频特性； 切比雪夫逼近法要求采用切比雪夫Ⅰ型，同时要求给出滤波器的脉冲响应、幅频特性和误差特性。 （3）将要求（1）和（2）中设计的具有相同参数要求，但采用不同设 计方法的滤波器进行比较，并以图的形式直观显示不同设计设计方法得到的数字低通滤波器的幅频特性的区别。 四、实验步骤 1．熟悉MATLAB 运行环境，命令窗口、工作变量窗口、命令历史记录窗口，FIR 常用基本函数； 2．熟悉MATLAB 文件格式，m 文件建立、编辑、调试； 3．根据要求（1）的内容，设计FIR 数字低通滤波器，建立M 文件，编写、调试、运行程序； 4．根据要求（2）的内容，设计FIR 数字低通滤波器，建立M 文件，编写、调试、运行程序； 5．将要求（1）和（2）中设计的具有相同参数要求，但采用不同设计方法的滤波器进行比较分析； 6．记录实验结果； 7．分析实验结果； 8．书写实验报告。 五、实验预习思考题 1．FIR 滤波器有几种常用设计方法？这些方法各有什么特点？

matlab中关于数字滤波器的函数介绍1

MATLAB下的数字信号处理实现示例 一信号、系统和系统响应 1、理想采样信号序列 （1）首先产生信号x(n),0<=n<=50 n=0:50;%定义序列的长度是50 A=444.128;%设置信号有关的参数 a=50*sqrt(2.0)*pi; T=0.001;%采样率 w0=50*sqrt(2.0)*pi; x=A*exp(-a*n*T).*sin(w0*n*T);%pi是MATLAB定义的π，信号乘可采用“.*”close all%清除已经绘制的x(n)图形 subplot(3,1,1);stem(x);%绘制x(n)的图形 title(‘理想采样信号序列’); （2）绘制信号x(n)的幅度谱和相位谱 k=-25:25; W=(pi/12.5)*k; X=x*(exp(-j*pi/12.5)).^(n’*k); magX=abs(X);%绘制x(n)的幅度谱 subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’); angX=angle(X);%绘制x(n)的相位谱 subplot(3,1,3);stem(angX);title(‘理想采样信号序列的相位谱’)（3）改变参数为：1,0734.2,4.0,10==Ω==TAα n=0:50;%定义序列的长度是50 A=1;%设置信号有关的参数 a=0.4; T=1;%采样率 w0=2.0734; x=A*exp(-a*n*T).*sin(w0*n*T);%pi是MATLAB定义的π，信号乘可采用“.*”close all%清除已经绘制的x(n)图形 subplot(3,1,1);stem(x);%绘制x(n)的图形 title(‘理想采样信号序列’); k=-25:25; W=(pi/12.5)*k; X=x*(exp(-j*pi/12.5)).^(n’*k); magX=abs(X);%绘制x(n)的幅度谱 subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’); angX=angle(X);%绘制x(n)的相位谱 subplot(3,1,3);stem(angX);title(‘理想采样信号序列的相位谱’) 2、单位脉冲序列 在MatLab中，这一函数可以用zeros函数实现： n=1:50;%定义序列的长度是50 x=zeros(1,50);%注意：MATLAB中数组下标从1开始

实验四(IIR数字滤波器设计及软件实现)讲解学习

实验四(I I R数字滤波器设计及软件实现)

10.4实验四IIR数字滤波器设计及软件实现 10.4.1 实验指导 1．实验目的 （1）熟悉用双线性变换法设计IIR数字滤波器的原理与方法； （2）学会调用MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具fdatool）设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。 （3）掌握IIR数字滤波器的MATLAB实现方法。 （4）通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。 2．实验原理 设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。基本设计过程是：①先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标；②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。第六章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫2和椭圆模拟和数字滤波器。本实验要求读者调用如上函数直接设计IIR数字滤波器。 本实验的数字滤波器的MATLAB实现是指调用MATLAB信号处理工具箱函数filter对给定的输入信号x(n)进行滤波，得到滤波后的输出信号y(n）。 3. 实验内容及步骤 （1）调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st，该函数还会自动绘图显示st的时域波形和幅频特性曲线，如图

10.4.1所示。由图可见，三路信号时域混叠无法在时域分离。但频域是分离的，所以可以通过滤波的方法在频域分离，这就是本实验的目的。 图10.4.1 三路调幅信号st 的时域波形和幅频特性曲线 （2）要求将st 中三路调幅信号分离，通过观察st 的幅频特性曲线，分别确定可以分离st 中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率。要求滤波器的通带最大衰减为0.1dB,阻带最小衰减为60dB 。 提示：抑制载波单频调幅信号的数学表示式为 0001()cos(2)cos(2)[cos(2())cos(2())]2 c c c s t f t f t f f t f f t ππππ==-++ 其中，cos(2)c f t π称为载波，f c 为载波频率，0cos(2)f t π称为单频调制信号，f 0为调制正弦波信号频率，且满足0c f f >。由上式可见，所谓抑制载波单频调幅信号，就是2个正弦信号相乘，它有2个频率成分：和频0c f f +和差频0c f f -，这2个频率成分关于载波频率f c 对称。所以，1路抑制载波单频调幅 信号的频谱图是关于载波频率f c 对称的2根谱线，其中没有载频成分，故取名 为抑制载波单频调幅信号。容易看出，图10.4.1中三路调幅信号的载波频率分