自适应中值滤波器的设计与实现

自适应中值滤波器的设计与实现
自适应中值滤波器的设计与实现

北京邮电大学世纪学院

毕业设计(论文)

题目自适应中值滤波器的设计与实现

学号 08010321

学生姓名王立阳

专业名称通信工程

所在系(院)通信与信息工程系

指导教师鞠磊

2012年 6 月 1 日

北京邮电大学世纪学院毕业设计(论文)任务书

姓名王立阳学号08010321 专业通信工程系(院)通信与信息工程系设计(论文)题目自适应中值滤波器的设计与实现

□工程设计;■工程技术研究;□软件工程(如CAI课题等);□专题研究;□艺术设计;□其题目分类

□自然科学基金与部、省、市级以上科研课题;□企、事业单位委托课题;□院级课题;■自拟题目来源

课题□其他

指导教师(指导教组

职称工作单位备注

组长及成员姓名)

鞠磊讲师北京电子科技学院指导教师

毕业设计(论文)的内容和要求:

[注意:选题尽量与实际应用需求相结合。要求写明本设计(论文)所涉及的分析方法或技术手段(如定性、定量分析的方法);要求有学生独立的见解,设计内容要详细写明具体步骤和技术指标]。

图像滤波是图像处理的关键步骤,常用于图像增强图像分割前的预处理。当前滤波方法非常多,各种算法在特定的情况下会表现出不同的效果。通常滤波器将图像不加区别的作为一个整体处理,而不考虑图像细节差异,因此在滤除噪声的同时也不可避免的模糊了图像细节。

本课题主要研究针对图像细节特点使用中值图像滤波的方法,具体内容包括:

1、图像滤波器原理与基本方法研究;

2、分析标准中值滤波方法存在的不足;

3、提出2到3种改进方法,设计自适应中值滤波;

4、通过MATLAB编程实现,实现对具体图像的滤波,并与传统中值滤波器进行效果比较分析。

应完成的工作和提交材料要求(课题完成后应提交成果的种类、数量、质量等方面的要求):

(1)开题报告:调研自适应滤波器的相关技术,针对课题要求制定研究内容与实施计划,撰写3000

字左右开题报告;

(2)论文的中文摘要:200-300字左右,包含关键词,并译成英文。英文摘要以250个左右实词为宜;(3)论文正文不少于15000字;

(4)翻译1500汉字以上的与课题有关的技术资料或专业文献;

主要参考文献(参考文献不少于4篇,参考文献目录按GB/T7714—2005的要求填写):

[1] 1S.Haykin著.郑宝玉等译.自适应滤波器原理,北京:第四版.电子工业出版社.2003.7

[2] 何振亚,自适应信号处理,北京:科学出版社,200.2

[3] 邹国良,自适应滤波理论及应用,石家庄:河北大学出版社,199.7

[4] 胡广书,数字信号处理—理论、算法与实现,北京:清华大学出版社,1997.

[5] 赵力,语音信号处理,北京:机械工业出版社,2003.

[6] 姚天任,数字语音处理,武汉,华中科技大学出版社,199.2

[7] 葛良,陶智,基于自适应滤波的语音增强算法,江苏:苏州大学学报.2002.8

[8] 韩利竹,王华,MATALB电子仿真与应用(第2版),北京:国防工业出版社,2003.

毕业设计(论文)进度计划(从正式启动时间开始,以周为单位填写):

第 1周-第 2周课题调研、查资料、撰写开题报告

第 3周根据查询的资料确定总体设计思路,完成开题报告并上交.

第 4周-第 7周毕业设计单元部分研究,并设计出整体框架

第 8周完成论文中期检查报告

第 9周-第15周资料整理,撰写毕业论文;上交毕业设计论文,指导教师审查评阅设计报告,毕业设计答辩资格审查。毕业设计答辩,学生修改毕业设计论文,准备答辩。

第16周进行毕设答辩。

指导教师签字:

日期:年月日

教学

单位

意见

审核人签字:

系(院)(盖章)

年月日学

审核专家签字:

年月日

备注1、由指导教师撰写,可根据长度加页,一式三份,教务处、系(院)各留存一份,发给学生一份,任务完成后附在论文内;

2、凡审核不通过的任务书,请重新申报。

北京邮电大学世纪学院

毕业设计(论文)诚信声明

本人声明所呈交的毕业设计(论文),题目《自适应中值滤波器的设计与实现》是本人在指导教师的指导下,独立进行研究工作所取得的成果,除了文中特别加以标注和致谢中所罗列的内容以外,毕业设计(论文)中不包含其他人已经发表或撰写过的研究成果,也不包含为获得北京邮电大学或其他教育机构的学位或证书而使用过的材料。

申请学位论文与资料若有不实之处,本人承担一切相关责任。

本人签名:日期:

毕业设计(论文)使用权的说明

本人完全了解北京邮电大学世纪学院有关保管、使用论文的规定,其中包括:①学校有权保管、并向有关部门送交学位论文的原件与复印件;②学校可以采用影印、缩印或其它复制手段复制并保存论文;③学校可允许论文被查阅或借阅;④学校可以学术交流为目的,复制赠送和交换学位论文;⑤学校可以公布学位论文的全部或部分内容。

本人签名:日期:

指导教师签名:日期:

题目自适应中值滤波器的设计与实现

摘要

图像是一种重要的信息源,通过图像处理可以帮助人们了解信息的内涵。然而在图像使用传输过程中,不可避免会受到噪声的干扰。中值滤波器是能够有效滤除脉冲噪声的非线性滤波器,但标准中值滤波去脉冲噪声的性能受滤波窗口尺寸的影响较大, 在抑制图像噪声和保护细节两方面存在矛盾。本文在研究中值滤波器的基础上,给出了一种能够根据噪声位置进行中值滤波器的自适应中值滤波方法,该方法可以有效地克服中值滤波器造成图像边缘模糊的缺点。最后通过matlab对所提出方法进行了仿真测试,结果验证了所提方法的有效性。

关键词图像处理中值滤波自适应 MATLAB

Title Adaptive median filter of design and implementation

Abstract

The image is a kind of important information source. And the image-processing can help people understand the connotation of the information. However, in the process of image transmission, the image is inevitably affected by the noise. Median Filter is the nonlinear impulse noise filter that is able to effectively filter out noise. But the performance of standard Median Filter largely affected by the size filtering window. So there are contradictions between the image noise suppression and image detail protection .And this paper, based on the study of the basis of the Median Filter, presented a noise-positioning and self-adaptive median filtering method which can effectively overcome shortcomings of fuzzy edges in the images. At last the simulation test, based on a method from MATLAB , demonstrates the validity of the method.

Keywords Image Median Fltering Adaptive MATLAB

目录

1.前言 (1)

1.1自适应滤波器 (1)

1.2自适应滤波器研究现状 (1)

1.3论文安排 (2)

2. 数字图像处理简介 (3)

2.1数字图像处理技术的优点 (3)

2.2制约数字图像处理技术发展的因素如下 (4)

2.3数字图像处理的主要功能 (4)

2.4图像处理方法 (4)

2.5图像的类型 (5)

2.6图像的常见格式 (6)

2.6.1图像的常见格式和类型 (6)

2.7应用领域 (8)

2.8.现代研究方向 (10)

3. 图像滤波技术 (11)

3.1中值滤波 (12)

3.2图像降噪 (13)

3.2.1图象降噪的方法简介 (13)

4. 自适应中值滤波方法 (15)

4.1常规的中值滤波器 (15)

4.1.1一维中值滤波的原理 (15)

4.1.2二维中值滤波的原理 (17)

4.1.3中值滤波的主要特性 (18)

4.2自适应中值滤波器 (19)

4.2.1自适应中值滤波器算法 (20)

4.2.2噪声检测 (20)

4.2.3噪声滤波 (22)

5. 自适应滤波算法的结果分析 (25)

6. 结论 (27)

致谢 (28)

参考文献 (29)

附录 (30)

1.前言

1.1自适应滤波器

滤波技术广泛地用于图象的预处理, 抑制图象噪声, 增强对比度, 以及强化图象的边沿特征. 运用较为广泛的线性滤波器如平均值滤波器, 能较好地抑制图象中的加性噪声. 但是, 线性滤波器会引起图象的钝化或模糊, 使得图象中物体边界产生位移. 特别是, 在图象受到乘性噪声或脉冲噪声的干扰, 如超声波及雷达成像中普遍存在的斑点噪声, 线性滤波器就不能取得预期的效果. 如何设计合适的滤波器来去除脉冲噪声并且保持图象的局部细节, 是一个重要且实用的课题. 中值滤波是由Tukey 首先用于一维信号处理, 后来被用到二维图象的平滑中.该滤波器对脉冲干扰及椒盐噪声的抑制效果好, 同时能够保持图象的边沿特征. 中值滤波器是对一个滑动窗口内的象素灰度值排序, 用其中间值代替窗口中心象素的灰度, 它是一种非线性滤波器. 例如滤波窗口由3×3个象素组成, 则其中5 个象素的灰度值会小于等于该滤波器的输出灰度值, 同时5 个象素的灰度值会大于等于滤波器的输出. 由此可见, 对于离散的脉冲信号, 当其出现的次数小于窗口尺寸的一半时, 将被抑制掉. 中值滤波法也能较好地保证图象的边沿特征, 而且易于实现. 因此它被广泛地应用于图象处理。

1.2自适应滤波器研究现状

图像在形成、传输过程中, 常因外界噪声干扰而导致其质量退化, 为减小噪声的影响, 可采取各种滤波方法对图像进行去噪处理,中值滤波由于可对长拖尾概率分布的噪声起到良好的平滑效果且可对图像中的某些细节起到保护作用, 因而在图像降噪处理中得到了比较广泛的应用,但标准中值滤波去脉冲噪声的性能受滤波窗口尺寸的影响较大, 在抑制图像噪声和保护细节两方存在一定的矛盾[1][2]:窗口滤波小, 可较好地保护图像中某些细节, 但滤除噪声的能力会受到限制; 滤波窗口大, 可加强噪声抑制能力, 但对细节的保护能力会减弱, 有时会滤去图像中的一些细线、尖锐边角等重要细节, 从而破坏图像的几何结构这种矛盾在图像中噪声干扰较大时表现得尤为明显。

自适应中值滤波器的滤波方式和常规的中值滤波器一样,都使用一个矩形区域,不同的是在滤波过程中,自适应滤波器会根据一定的设定条件改变(即增的窗口S

xy

加)滤波窗的大小,同时当判断滤波窗中心的像素是噪声时,该值用中值代替,否则

不改变其当前像素值,这样用滤波器的输出来替代像素(x,y)处(即目前滤波窗中心的坐标)的值。自适应中值滤波器可以处理噪声概率更大的脉冲噪声,同时能够更好地保持图像细节,这是常规中值滤波器做不到的。

自适应中值滤波总体上可以分为三步:(1)对图像各区域进行噪声检测;(2)根据各区域受噪声污染的状况确定滤波窗口的尺寸;(3)对检测出的噪声点进行滤波。1.3 论文安排

本文首先介绍数字图像处理的基本理论,然后具体研究了中值滤波器的理论。最后提出了自适应中值滤波方法, 并通过matlab对所提出方法进行了验证。

第一章对自适应中值滤波器进行了简介;

第二章介绍了图像处理的基本概念;

第三章重点研究了图像滤波技术;

第四章在前期研究基础上提出了一种自适应中值滤波方法及其仿真实现。

2. 数字图像处理简介

图像是自然界景物的客观反映,是人类认识世界和人类本身的重要源泉[1]。据计,人类获取的信息,绝大多数来自视觉,因此,图像作为视觉信息获取和传播的要载体,随着科技的飞速发展,也越来越受到人们的重视。目前,图像处理技术已为工程学、计算机科学、信息科学、统计学、物理学、化学、生物学、医学甚至社科学等领域各学科之间学习和研究的对象。它己给人类带来了巨大的经济和社会效,不久的将来它不仅在理论上会有更深入的发展,在应用上也是科学研究、社会生产至人类社会生活中不可缺少的强有力工具。

2.1数字图像处理技术的优点

数字图像处理技术的优点:

1.再现性好。数字图像处理与模拟图像处理的根本不同在于,它不会因图像的存储、传输或复制等一系列变换操作而导致图像质量的退化。只要图像在数字化准确地表现了原稿,则数字图像处理过程始终能保持图像的真实再现。

2.处理精度高。按目前的技术,几乎可将一幅模拟图像数字化为任意大小的二维数组,这主要取决于图像数字化设备的能力。现代扫描仪可以把每个像素的灰度级量化为16位甚至更高,这意味着图像的数字化精度可以满足任一应用需求。从原上讲不论图像的精度有多高,处理总是能实现的,只要在处理时改变程序中的数组数就可以了。

3.适用面宽。图像可以来自多种信息源。从图像反映的客观实体尺度看,以小到电子显微镜图像,大到航空照片、遥感图像甚至天文望远镜图像。这些来自同信息源的图像只要被变换为数字编码形式后,均是用二维数组表示的灰度图像组而成,因而均可用计算机来处理。

4.灵活性高。由于图像的光学处理从原理上讲只能进行线性运算,极大地制了光学图像处理能实现的目标;而数字图像处理不仅能完成线性运算,而且能实非线性处理,即凡是可以用数学公式或逻辑关系来表达的一切运算均可用数字图像处理实现。

5.信息压缩的潜力大。数字图像中各个像素是不独立的,其相关性大。在图像画面上,经常有很多像素有相同或接近的灰度。就电视画面而言,同一行中相邻两个像素或相邻两行间的像素,其相关系数可达0.9以上,而相邻两帧之间的相关性

和帧内相关性相比,一般来说相邻两帧之间的相关性还要大些。因此,图像处理中信息压缩的潜力很大。

2.2 制约数字图像处理技术发展的因素如下

1.数字图像处理占用的频带较宽。与音频信息相比,占用的频带要大几个数量级。如电视图像的带宽约5.6MHz,而语音带宽仅为4kHz左右。所以在成像、传输、存储、处理、显示等各个环节的实现上,技术难度较大,成本高。

2.由于图像是三维景物的二维投影,一幅图像本身不具备复现三维景物的全部几何信息的能力,三维景物背后部分信息在二维图像画面上是反映不出来的。因此,要分析和理解三维景物必须作合适的假定或附加新的测量;在理解三维景物时需要

知识导引,这也是人工智能中正在致力解决的知识工程问题。

3.图像处理出来的图片最终是给人观察和评价的,因此它受人的因素影响较大,由于人的视觉系统很复杂,受环境条件、视觉性能、人的情绪、兴趣及知识状况影响很大,因此,要求图像处理系统与人的视觉系统要有良好的“匹配”,但是目前要做到这一点是有很大困难的,因为到目前为止人们对人的视觉系统研究还很不充分,必须在进一步研究视觉的生理和心理特征和图像的统计规律之后,才能很好地解决这个问题。

2.3数字图像处理的主要功能

数字图像处理的主要功能是将获得的低质量图像(反差小,模糊,有噪声等)用计算机处理成更适合人眼观察、实际应用或仪器检测的图像。

随着计算机技术的飞速发展、高性能计算机的普及和计算机软件水平的提高,特别是多媒体技术、通信技术、信息存储技术和以Internet为代表的计算机网络技术的加速发展和广泛普及,数字图像处理获得飞速的发展,形成了许多重要的图像处理分支。总的来说,图像处理的分支可以分为三部分,对应于不同的科学研究目的。第一个目的是对图像的最小物理存储的寻找,这个学科叫做图像压缩;第二个目的是把一个带有噪声、模糊或有其它扰动的图像复原;第三个目的是分析(analysis),我们说的图像分析(image analysis),它的目的不是去除噪声或压缩,而是找出主要部分的不变量并且可以从大量的图形数据库中辨识物体。

2.4图像处理方法

图像处理方法也很多,主要包括图像数据编码与压缩、图像增强与复原、图像传

输、平滑、边缘锐化、图像分割、特征抽取、图像去噪、图像识别与理解等内容。经过这些处理,输出图像的质量得到相当程度的改善,既提高了图像的视觉效果,又便于计算机对图像进行分析、处理和识别。

图像的表现形式很多,但都有一个共同的特点,即图像是二维或三维空间信息。为了本文研究的方便,这里我们主要来介绍一下二维图像。平面上每一点仅包括两个坐标值,因此,平面图像函数是连续的二维函数,如式(2.1)所示:

()()()(){}y x fblue y x fgreen y x fred y x f ,,,,,,= 式(2.1) 数字图像是连续图像f(x,y)的一种近似表示,通常用由采样点的值所组成的矩

阵来表示,如图2-1所示:

()()()()()()()()()????????????------1,1...1,10,1.........

...1,1...1,10,11,...1,00,0M N f N f N f M f f f M o f f f 图2-1数字图像的矩阵表示图

2.5图像的类型

图像类型是指数组数值与像素颜色之间定义的关系,它与图像格式的概念有所不 同,常见的有四种类型的图像。

2.5.1二进制图像

在一幅二进制图像中,每一个像素将取两个离散数值(0或1)中的一个,从本质上说,这两个数值分别代表状态开或关。因为图像中的每个像素仅需1位信息,常把二进制图像称为是1位图像。单色图像具有比较简单的格式.一般由黑色区域和白色区域组成。

2.5.2灰度图像

灰度图像具有如下特征:

1.灰度图像的存储文件带有图像颜色表,此颜色表共有256项,图像颜色表中一表项由红、绿、蓝颜色分量组成,并且红、绿、蓝颜色分量值都相等。

2.每个像素由8位组成,其值范围从O-255,表示256种不同的灰度级。

2.5.3伪彩色图像

伪彩色图像与灰度图像相似,其存储文件中也带有图像颜色表,伪彩色图像具有

如下特点:

1.图像颜色表中的红、绿、蓝颜色分量值不全相等。

2.整幅图像仅有256种颜色,要表示256种不同的颜色,像素必须由8位组成,每个像素的像素值不是由每个基色分量的数值直接决定的,我们把具有256色的图像称为8位彩色图像。

2.5.4位真彩色图像

具有全彩色照片表达能力的图像为24位彩色图像,24位真彩色图像存储文件中不带有图像颜色表,其具有如下特征:

1.图像中每一像素由RGB三个分量组成,每个分量各占8位,每个像素需24位。

2.红、绿、蓝三个分量的取值范围是O-255。

2.6图像的常见格式

图像文件按其数字图像格式的不同一般具有不同的扩展名。最常见的是:BMP、GIF、TIFF、JPEG等[1]。

2.6.1图像的常见格式和类型

图像文件按其数字图像格式的不同一般具有不同的扩展名。最常见的是:BMP、GIF、TIFF、JPEG等[1][36]。下面简单介绍一下。

2.6.1.1 BMP(Bitmap)文件

BMP文件是Microsoft Windows所定义的图像文件格式,最早应用在微Microsoft Windows窗口系统中。BMP图像文件具有以下特点:(1)只存放一副图像;(2)只能存储单色、16色、256色和真彩色四种图像数据;(3)图像数据有压缩两种处理方式;(4)调色板的数据存储结构较为特殊,存储格式不是固定与文件头的某些具体参数(如像素位、压缩算法等)密切相关的。其中,Windows 设有RLE4和RLE8两种压缩方式,RLE4只能处理16色图像数据,RLE8则只能压缩256色图像数据。

BMP图像文件的文件结构由文件头、位图信息头、颜色信息和图像数据成其中,表头长度固定为54个字节,只有真彩色BMP图像文件内没有调色其余不超过256种颜色的图像文件都必须有调色板信息。

2.6.1.2 GIF文件

GIF图像文件格式是CompuServe公司最先在网络中用于在线传送图像数据。GIF

图像文件经常用于网页的动画等特技制作。具有以下特点:(1)文件具有多元化结构,能够储存多张图像,多图像的定序或覆盖,交错屏幕绘图以及文本覆盖;(2)调色板数据有通用调色板和局部调色板之分;(3)采用了LZW压缩法;(4)图像数据一个字节存储一点;(5)文件内的各种图像数据区和补充区多数没有固定的数据长度和存储位置,为了方便程序寻找数据区,就以数据区的第一个字节作为标识符,以使程序能够判断读到哪种数据区;(6)图像数据有顺序排列和交叉排列两种方式;(7)图像最多只能存储256色图像。

GIF图像文件结构一般由表头、通用调色板、图像数据区以及四个补充区共七个数据单元组成。其中,表头和图像数据区是文件不可缺少的单元,通用调色板和其余的四个补充区是可选择内容。

2.6.1.3 TIF文件

TIF图像文件格式是现有图像文件格式中最复杂的一种,其主要特点如下:(1)应用指针功能,实现多幅图像存储;(2)文件内数据区没有固定的排列顺序,但规定表头必须在文件前端,标识信息区和图像数据区在文件中可以任意存放;(3)可制定私人用的标识信息;(4)能够接受除了一般图像处理RGB模式之外的多种不同的图像模式;(5)可存储多份调色板数据,其调色板的数据类型和排列顺序较为特殊;(6)能够提供多种不同的压缩数据的方法;(7)图像数据可分割成几个部分进行分别存档。

TIF图像文件主要由表头、标识信息区和图像数据区三部分组成。其中,文件内固定只有一个位于文件前端的表头,表头有一个标志参数指出标识信息区在文件中的存储地址,标识信息区有多组标识信息用于存储图像数据区的地址。每组标识信息长度固定为12个字节,前8个字节分别表示标识信息的代号(2个字节)、数据类型(2个字节)、数据量(4个字节),最后4个字节用于存储数据值或标志信息。

2.6.1.4 JPEG格式

JPEG是对静止灰度或彩色图像的一种国际压缩标准,其正式的名称是“连续色调静态图像的数学压缩和编码”,己在数字照相机上得到广泛使用,当选用有损压缩方式时其可节省相当大的空间。JPGE标准只是定义了一个规范的编码数据流,并没有规定图像数据文件的格式。Cube Microsystems公司定义了一种JPGE文件交换格式JFIF。JFIF图像是一种或者使用灰度表示,或者使用Y,Cb,Cr分量彩色表示的JPEG图像。它包含一个与JPEG兼容的头。一个JFIF文件通常包含单个图像,图像

可以是灰度的,也可以是彩色的。

图像类型是指数组数值与像素颜色之间定义的关系,它与图像格式的概念有所不同,常见的有四种类型的图像。

2.7应用领域

图像是人类获取和交换信息的主要来源,因此,图像处理的应用领域必然涉及到人类生活和工作的方方面面。随着人类活动范围的不断扩大,图像处理的应用领域也将随之不断扩大。

1.航天和航空方面

航天和航空技术方面的应用数字图像处理技术在航天和航空技术方面的应用,除了JPL对月球、火星照片的处理之外,另一方面的应用是在飞机遥感和卫星遥感技术中。许多国家每天派出很多侦察飞机对地球上有兴趣的地区进行大量的空中摄影。对由此得来的照片进行处理分析,以前需要雇用几千人,而现在改用配备有高级计算机的图像处理系统来判读分析,既节省人力,又加快了速度,还可以从照片中提取人工所不能发现的大量有用情报。从60年代末以来,美国及一些国际组织发射了资源遥感卫星(如LANDSAT系列)和天空实验室(如SKYLAB),由于成像条件受飞行器位置、姿态、环境条件等影响,图像质量总不是很高。因此,以如此昂贵的代价进行简单直观的判读来获取图像是不合算的,而必须采用数字图像处理技术。如LANDSAT系列陆地卫星,采用多波段扫描器(MSS),在900km高空对地球每一个地区以18天为一周期进行扫描成像,其图像分辨率大致相当于地面上十几米或100米左右(如1983年发射的LANDSAT-4,分辨率为30m)。这些图像在空中先处理(数字化,编码)成数字信号存入磁带中,在卫星经过地面站上空时,再高速传送下来,然后由处理中心分析判读。这些图像无论是在成像、存储、传输过程中,还是在判读分析中,都必须采用很多数字图像处理方法。现在世界各国都在利用陆地卫星所获取的图像进行资源调查(如森林调查、海洋泥沙和渔业调查、水资源调查等),灾害检测(如病虫害检测、水火检测、环境污染检测等),资源勘察(如石油勘查、矿产量探测、大型工程地理位置勘探分析等),农业规划(如土壤营养、水份和农作物生长、产量的估算等),城市规划(如地质结构、水源及环境分析等)。中国也陆续开展了以上诸方面的一些实际应用,并获得了良好的效果。在气象预报和对太空其它星球研究方面,数字图像处理技术也发挥了相当大的作用。

2.生物医学工程方面

数字图像处理在生物医学工程方面的应用十分广泛,而且很有成效。除了上面介绍的CT技术之外,还有一类是对医用显微图像的处理分析,如红细胞、白细胞分类,染色体分析,癌细胞识别等。此外,在X光肺部图像增晰、超声波图像处理、心电图分析、立体定向放射治疗等医学诊断方面都广泛地应用图像处理技术。

3.通信工程方面

当前通信的主要发展方向是声音、文字、图像和数据结合的多媒体通信。具体地讲是将电话、电视和计算机以三网合一的方式在数字通信网上传输。其中以图像通信最为复杂和困难,因图像的数据量十分巨大,如传送彩色电视信号的速率达100Mbit/s以上。要将这样高速率的数据实时传送出去,必须采用编码技术来压缩信息的比特量。在一定意义上讲,编码压缩是这些技术成败的关键。除了已应用较广泛的熵编码、DPCM编码、变换编码外,目前国内外正在大力开发研究新的编码方法,如分行编码、自适应网络编码、小波变换图像压缩编码等。

4.工业和工程方面

在工业和工程领域中图像处理技术有着广泛的应用,如自动装配线中检测零件的质量、并对零件进行分类,印刷电路板疵病检查,弹性力学照片的应力分析,流体力学图片的阻力和升力分析,邮政信件的自动分拣,在一些有毒、放射性环境内识别工件及物体的形状和排列状态,先进的设计和制造技术中采用工业视觉等等。其中值得一提的是研制具备视觉、听觉和触觉功能的智能机器人,将会给工农业生产带来新的激励,目前已在工业生产中的喷漆、焊接、装配中得到有效的利用。

5.军事公安方面

在军事方面图像处理和识别主要用于导弹的精确末制导,各种侦察照片的判读,具有图像传输、存储和显示的军事自动化指挥系统,飞机、坦克和军舰模拟训练系统等;公安业务图片的判读分析,指纹识别,人脸鉴别,不完整图片的复原,以及交通监控、事故分析等。目前已投入运行的高速公路不停车自动收费系统中的车辆和车牌的自动识别都是图像处理技术成功应用的例子。

6.文化艺术方面

目前这类应用有电视画面的数字编辑,动画的制作,电子图像游戏,纺织工艺品设计,服装设计与制作,发型设计,文物资料照片的复制和修复,运动员动作分析和评分等等,现在已逐渐形成一门新的艺术--计算机美术。

7.机器人视觉

机器视觉作为智能机器人的重要感觉器官,主要进行三维景物理解和识别,是目前处于研究之中的开放课题。机器视觉主要用于军事侦察、危险环境的自主机器人,邮政、医院和家庭服务的智能机器人,装配线工件识别、定位,太空机器人的自动操作等。

8.视频和多媒体系统

目前,电视制作系统广泛使用的图像处理、变换、合成,多媒体系统中静止图像和动态图像的采集、压缩、处理、存贮和传输等。

9.科学可视化

图像处理和图形学紧密结合,形成了科学研究各个领域新型的研究工具。

10.电子商务

在当前呼声甚高的电子商务中,图像处理技术也大有可为,如身份认证、产品防伪、水印技术等。

总之,图像处理技术应用领域相当广泛,已在国家安全、经济发展、日常生活中充当越来越重要的角色,对国计民生的作用不可低估。

研究方向

2.8.现代研究方向

自20世纪60年代第三代数字计算机问世以后,数字图像处理技术出现了空前的发展,在该领域中需要进一步研究的问题主要有如下五个方向:

1.在进一步提高精度的同时着重解决处理速度问题;

2.加强软件研究,开发新的处理方法,特别要注意移植和借鉴其他学科的技术和研究成果,创造新的处理方法;

3.加强边缘学科的研究工作,促进图像处理技术的发展;

4.加强理论研究,逐步形成处理科学自身的理论体系;

5.时刻注意图像处理领域的标准化问题。

3. 图像滤波技术

人们根据实际图像的特点、噪声的统计特征和频谱分析的规律,发展了各式各样的去噪方法。其中最直观的方法是根据噪声能量一般集中于高频,而图像频谱分布于一个有限区间的这一特点,采用低通去噪方法,例如wiener 线性滤波器。此外,对图像进行平滑处理也是常用的方法。平滑的目的有两个:改善图像质量和抽出对象特征。平滑可以在空间域进行,也可以在频率域进行。由于噪声源众多(如光栅扫描、底片颗粒、机械元件、信道传输等),噪声种类复杂(如加性噪声、乘性噪声、量化噪声、高斯噪声、均衡噪声和脉冲噪声等),所以相应的平滑方法也多种多样。其中空域平滑是数字图像处理领域比较活跃的一个分支,己经经历了几十年的研究探索,并形成了较为系统的算法。总的来说,其基本思想是用所选的邻域中的各像素灰度的平均值来代替中心像素的灰度值。但是,空域平滑算法有一个共同的不足,就是它不仅平滑噪声,而且使图像中的细节模糊化。例如,邻域平均法就是空间域平滑噪声技术,它是一种线性滤波技术,它的基本原理是对于给定的图中的每个像点(m ,n),取其邻域S ,设邻域S 含有M 个像素,取其平均值作为处理后所得图像点(m ,n)处的灰度.用一个像素邻域内各像素灰度平均值来代替该像素原来的灰度。邻域S 的形状和大小根据图像特点确定,一般取的形状是正方形、矩形及十字形,如S 为3x3邻域,点m(,n)位于S 中心,则如式(3.1)所示:

()()∑∑==++=313

1

,91,i j j n i m f n m f 式(3-1) 但是这种线性滤波一般通过取模板做离散卷积实现,这种方法在平滑脉冲噪声点的同时也导致图像更加模糊,效果不是很好。如果邻域取得越大,图像的模糊程度也愈加严重。

1971年,J.W.TuKey 针对离散数据平滑问题首先提出了中值滤波的概念[4][7],其后这种新的滤波思想很快就被引入到一些重要的数字信号处理领域,其中图像处理是最主要的一个方面。当标准中值滤波器用于解决图像恢复问题时,人们注意到它有两个独特的性能:(1)能较好地保护图像细节;(2)有很好的韧性,能较好地抑制远偏离高斯型的,甚至不完全独立于有用信号的各种噪声。这些性能是人们所熟悉的线性滤波器所不具备的。为了解释标准中值滤波器的这些性能,人们开始研究标准中值滤波的特性,并且取得了许多重要结论。理论分析和实际应用的结果都使人们

最新自适应滤波器的设计开题报告

长江大学 毕业设计开题报告 题目名称自适应滤波器的设计与应用学院电信学院 专业班级信工10702班 学生姓名李雪利 指导教师王圆妹老师 辅导教师王圆妹老师 开题报告日期 2010年3月19日

自适应滤波器的设计与应用 学生:李雪利,长江大学电子信息学院 指导教师:王圆妹,长江大学电子信息学院 一、题目来源 来源于其他 二、研究目的和意义 滤波技术在当今信息处理领域中有着极其重要的应用。滤波是从连续的或离散的输入数据中除去噪音和干扰以提取有用信息的过程,相应的装置就称为滤波器。滤波器实际上是一种选频系统,他对某些频率的信号予以很小的衰减,使该部分信号顺利通过。而对其他不需要的频率信号予以很大的衰减,尽可能阻止这些信号通过。滤波器研究的一个目的就是:如何设计和制造最佳的(或最优的)滤波器。 在数字信号处理中,数字滤波是语音和图像处理、模式识别、频谱分析等应用中的一个基本处理算法。在许多应用场合,由于无法预先知道信号和噪声的特性或者它们是随时间变化的,仅仅用 FIR 和 IIR两种具有固定滤波系数的滤波器无法实现最优滤波。在这种情况下,必须设计自适应滤波器,以跟踪信号和噪声的变化。 自适应滤波器是利用前一时刻已获得的滤波器参数,自动地调节、更新现时刻的滤波器参数,以适应信号和噪声未知的统计特性,从而实现最优滤波。当在未知统计特性的环境下处理观测信号时,利用自适应滤波器可以获得令人满意的效果,其性能远超过通用方法所设计的固定参数滤波器。

三、阅读的主要参考文献及资料名称 1、《数字信号处理》刘益成(第二版)西安电子科技出版社 2、《数字信号处理》张小虹(第二版)机械工业出版社 3、自适应信号处理[M].西安:西安电子科技大学出版社,2001. 4.邹理和,数字信号处理, 国防工业出版社,1985 5.丁玉美等, 数字信号处理,西安电子科技大学出版社,1999 6.程佩青, 数字信号处理,清华大学出版社,2001 7. The MathWorks Inc, Signal Processing Toolbox For Use with MATLAB, Sept. 2000 8. vinay K.Ingle, John G.Proakis,数字信号处理及MATLAB实现,陈怀琛等译,电子工业出版社,1998.9 9、《MATLAB编程参考手册》 10、中国期刊网的相关文献 11、赫金,自适应滤波器原理第四版,西安工业出版社,2010-5-1 四、国内外现状和发展趋势与主攻方向 自适应滤波器的理论与技术是50年代末和60年代初发展起来的。它是现代信号处理技术的重要组成部分,对复杂信号的处理具有独特的功能。自适应滤波器在数字滤波器中试属于随机数字信号处理的范畴。对于随机数字信号的滤波处理,通常有维纳滤波,卡尔曼滤波和自适应滤波,维纳滤波的权系数是固定的,适用于平稳随机信号;卡尔曼滤波器的权系数是可变的,适用于非平稳随机信号中。但是,只有在对信号和噪声的统计特性先验

绝对经典的低通滤波器设计报告

经典 无源低通滤波器的设计

团队:梦知队 团结奋进,求知创新,追求卓越,放飞梦想 队员: 日期:2010.12.10 目录 第一章一阶无源RC低通滤波电路的构建 (3) 1.1 理论分析 (3) 1.2 电路组成 (4) 1.3 一阶无源RC低通滤波电路性能测试 (5) 1.3.1 正弦信号源仿真与实测 (5) 1.3.2 三角信号源仿真与实测 (10) 1.3.3 方波信号源仿真与实测 (15) 第二章二阶无源LC低通滤波电路的构建 (21) 2.1理论分析 (21) 2.2 电路组成 (22) 2.3 二阶无源LC带通滤波电路性能测试 (23) 2.3.1 正弦信号源仿真与实测 (23) 2.3.2 三角信号源仿真与实测 (28)

2.3.3 方波信号源仿真与实测 (33) 第三章结论与误差分析 (39) 3.1 结论 (39) 3.2 误差分析 (40) 第一章一阶无源RC低通滤波电路的构建1.1理论分析 滤波器是频率选择电路,只允许输入信号中的某些频率成分通过,而阻止其他频率成分到达输出端。也就是所有的频率成分中,只是选中的部分经过滤波器到达输出端。 低通滤波器是允许输入信号中较低频率的分量通过而阻止较高频率的分量。 图1 RC低通滤波器基本原理图 当输入是直流时,输出电压等于输入电压,因为Xc无限大。当输入

频率增加时,Xc减小,也导致Vout逐渐减小,直到Xc=R。此时的频率为滤波器的特征频率fc。 解出,得: 在任何频率下,应用分压公式可得输出电压大小为: 因为在=时,Xc=R,特征频率下的输出电压用分压公式可以表述为: 这些计算说明当Xc=R时,输出为输入的70.7%。按照定义,此时的频率称为特征频率。 1.2电路组成

IIR数字滤波器设计原理

IIR 数字滤波器设计原理 利用双线性变换设计IIR 滤波器(只介绍巴特沃斯数字低通滤波器的设计),首先要设计出满足指标要求的模拟滤波器的传递函数)(s H a ,然后由)(s H a 通过双线性变换可得所要设计的IIR 滤波器的系统函数)(z H 。 如果给定的指标为数字滤波器的指标,则首先要转换成模拟滤波器的技术指标,这里主要是边界频率 s p w w 和的转换,对s p αα和指标不作变化。边界频率的转换关系为)21tan(2w T =Ω。接着,按照模拟低通滤波器的技术指标根据相应 设计公式求出滤波器的阶数N 和dB 3截止频率c Ω;根据阶数N 查巴特沃斯归一 化低通滤波器参数表,得到归一化传输函数 )(p H a ;最后,将c s p Ω=代入)(p H a 去归一,得到实际的模拟滤波器传输函数)(s H a 。之后,通过双线性变换法转换公式 11 112--+-=z z T s ,得到所要设计的IIR 滤波器的系统函数)(z H 。 步骤及内容 1) 用双线性变换法设计一个巴特沃斯IIR 低通数字滤波器。设计指标参数为: 在通带内频率低于π2.0时,最大衰减小于dB 1;在阻带内[]ππ,3.0频率区间上,最小衰减大于dB 15。 2) 以π02.0为采样间隔,绘制出数字滤波器在频率区间[]2/,0π上的幅频响应特 性曲线。 3) 程序及图形 程序及实验结果如下: %%%%%%%%%%%%%%%%%%

%iir_1.m %lskyp %%%%%%%%%%%%%%%%%% rp=1;rs=15; wp=.2*pi;ws=.3*pi; wap=tan(wp/2);was=tan(ws/2); [n,wn]=buttord(wap,was,rp,rs,'s'); [z,p,k]=buttap(n); [bp,ap]=zp2tf(z,p,k); [bs,as]=lp2lp(bp,ap,wap); [bz,az]=bilinear(bs,as,.5); [h,f]=freqz(bz,az,256,1); plot(f,abs(h)); title('双线性z 变换法获得数字低通滤波器,归一化频率轴'); xlabel('\omega/2\pi'); ylabel('低通滤波器的幅频相应');grid; figure; [h,f]=freqz(bz,az,256,100); ff=2*pi*f/100; absh=abs(h); plot(ff(1:128),absh(1:128)); title('双线性z 变换法获得数字低通滤波器,频率轴取[0,\pi/2]'); xlabel('\omega'); ylabel('低通滤波器的幅频相应');grid on; 运行结果: 00.050.10.150.20.25 0.30.350.40.450.500.1 0.2 0.3 0.40.50.60.70.8 0.9 1 双线性z 变换法获得数字低通滤波器,归一化频率轴 ω/2π低通滤波器的幅频相应

自适应滤波器的dsp实现

学号: 课程设计 学院 专业 年级 姓名 论文题目 指导教师职称 成绩 2013年 1 月 10 日

目录 摘要 (1) 关键词 (1) Abstract (1) Key words (1) 引言 (1) 1 自适应滤波器原理 (2) 2 自适应滤波器算法 (3) 3 自适应滤波算法的理论仿真与DSP实现 (5) 3.1 MATLAB仿真 (5) 3.2 DSP的理论基础 (7) 3.3 自适应滤波算法的DSP实现 (9) 4 结论 ............................................... 错误!未定义书签。致谢 ................................................. 错误!未定义书签。参考文献 ............................................. 错误!未定义书签。

自适应滤波器算法的DSP实现 学生姓名:学号: 学院:专业: 指导教师:职称: 摘要:本文从自适应滤波器的基本原理、算法及设计方法入手。本设计最终采用改进的LMS算法设计FIR结构自适应滤波器,并采用MATLAB进行仿真,最后用DSP 实现了自适应滤波器。 关键词:DSP(数字信号处理器);自适应滤波器;LMS算法;FIR结构滤波器 DSP implementation of the adaptive filter algorithm Abstract:In this article, starting from the basic principles of adaptive filter and algorithms and design methods. Eventually the design use improved the LMS algorithm for FIR adaptive filter,and use MATLAB simulation, adaptive filter using DSP. Key words:DSP;adaptive filter algorithm;LMS algorithm;FIR structure adaptive filter 引言 滤波是电子信息处理领域的一种最基本而又极其重要的技术。在有用信号的传输过程中,通常会受到噪声或干扰的污染。利用滤波技术可以从复杂的信号中提取所需要的信号,同时抑制噪声或干扰信号,以便更有效地利用原始信号。滤波器实际上是一种选频系统,它对某些频率的信号予以很小的衰减,让该部分信号顺利通过;而对其他不需要的频率信号则予以很大的衰减,尽可能阻止这些信号通过。在电子系统中滤波器是一种基本的单元电路,使用很多,技术也较为复杂,有时滤波器的优劣直接决定产品的性能,所以很多国家非常重视滤波器的理论研究和产品开发[1]。近年来,尤其数字滤波技术使用广泛,数字滤波理论的研究及其产品的开发一直受到很多国家的重视。从总的来说滤波可分为经典滤波和现代滤波。经典滤波要求已知信号和噪声的统计特性,如维纳滤波和卡尔曼滤波。现代滤波则不要求己知信号和噪声的统计特性,如自适应滤波。 自适应滤波器是统计信号处理的一个重要组成部分。在实际应用中,由于没有充足的信息来设计固定系数的数字滤波器,或者设计规则会在滤波器正常运行时改变,因此我们需要研究自适应滤波器。凡是需要处理未知统计环境下运算结果所产生的信

简单低通滤波器设计及matlab仿真

东北大学 研究生考试试卷 考试科目: 课程编号: 阅卷人: 考试日期: 姓名:xl 学号: 注意事项 1.考前研究生将上述项目填写清楚. 2.字迹要清楚,保持卷面清洁. 3.交卷时请将本试卷和题签一起上交. 4.课程考试后二周内授课教师完成评卷工作,公共课成绩单与试卷交研究生院培养办公室, 专业课成绩单与试卷交各学院,各学院把成绩单交研究生院培养办公室. 东北大学研究生院培养办公室

数字滤波器设计 技术指标: 通带最大衰减: =3dB , 通带边界频率: =100Hz 阻带最小衰减: =20dB 阻带边界频率: =200Hz 采样频率:Fs=200Hz 目标: 1、根据性能指标设计一个巴特沃斯低通模拟滤波器。 2、通过双线性变换将该模拟滤波器转变为数字滤波器。 原理: 一、模拟滤波器设计 每一个滤波器的频率范围将直接取决于应用目的,因此必然是千差万别。为了使设计规范化,需要将滤波器的频率参数作归一化处理。设所给的实际频 率为Ω(或f ),归一化后的频率为λ,对低通模拟滤波器令λ=p ΩΩ/,则1 =p λ, p s s ΩΩ=/λ。令归一化复数变量为p ,λj p =,则p p s j j p Ω=ΩΩ==//λ。所以巴 特沃思模拟低通滤波器的设计可按以下三个步骤来进行。 (1)将实际频率Ω规一化 (2)求Ωc 和N 11010/2-=P C α s p s N λααlg 1 10 110lg 10 /10/--= 这样Ωc 和N 可求。 p x fp s x s f

根据滤波器设计要求=3dB ,则C =1,这样巴特沃思滤波器的设计就只剩一个参数N ,这时 N p N j G 222 )/(11 11)(ΩΩ+= += λλ (3)确定)(s G 因为λj p =,根据上面公式有 N N N p j p p G p G 22)1(11 )/(11)()(-+= += - 由 0)1(12=-+N N p 解得 )221 2exp(πN N k j p k -+=,k =1,2, (2) 这样可得 1 )21 2cos(21 ) )((1 )(21+-+-= --= -+πN N k p p p p p p p G k N k k 求得)(p G 后,用p s Ω/代替变量p ,即得实际需要得)(s G 。 二、双线性变换法 双线性变换法是将s 平面压缩变换到某一中介1s 平面的一条横带里,再通过标准变换关系)*1exp(T s z =将此带变换到整个z 平面上去,这样就使s 平面与z 平面之间建立一一对应的单值关系,消除了多值变换性。 为了将s 平面的Ωj 轴压缩到1s 平面的1Ωj 轴上的pi -到pi 一段上,可以通过以下的正切变换来实现: )21 tan(21T T Ω= Ω 这样当1Ω由T pi -经0变化到T pi 时,Ω由∞-经过0变化到∞+,也映射到了整个Ωj 轴。将这个关系延拓到整个s 平面和1s 平面,则可以得到

自适应滤波器毕业设计论文

大学 数字信号处理课程要求论文 基于LMS的自适应滤波器设计及应用 学院名称: 专业班级: 学生姓名: 学号: 2013年6月

摘要自适应滤波在统计信号处理领域占有重要地位,自适应滤波算法直接决定着滤波器性能的优劣。目前针对它的研究是自适应信号处理领域中最为活跃的研究课题之一。收敛速度快、计算复杂性低、稳健的自适应滤波算法是研究人员不断努力追求的目标。 自适应滤波器是能够根据输入信号自动调整性能进行数字信号处理的数字滤波器。作为对比,非自适应滤波器有静态的滤波器系数,这些静态系数一起组成传递函数。研究自适应滤波器可以去除输出信号中噪声和无用信息,得到失真较小或者完全不失真的输出信号。本文介绍了自适应滤波器的理论基础,重点讲述了自适应滤波器的实现结构,然后重点介绍了一种自适应滤波算法最小均方误差(LMS)算法,并对LMS算法性能进行了详细的分析。最后本文对基于LMS算法自适应滤波器进行MATLAB仿真应用,实验表明:在自适应信号处理中,自适应滤波信号占有很重要的地位,自适应滤波器应用领域广泛;另外LMS算法有优也有缺点,LMS算法因其鲁棒性强特点而应用于自回归预测器。 关键词:自适应滤波器,LMS算法,Matlab,仿真

1.引言 滤波技术在当今信息处理领域中有着极其重要的应用。滤波是从连续的或离散的输入数据中除去噪音和干扰以提取有用信息的过程,相应的装置就称为滤波器。滤波器实际上是一种选频系统,他对某些频率的信号予以很小的衰减,使该部分信号顺利通过;而对其他不需要的频率信号予以很大的衰减,尽可能阻止这些信号通过。滤波器研究的一个目的就是:如何设计和制造最佳的(或最优的)滤波器。Wiener于20世纪40年代提出了最佳滤波器的概念,即假定线性滤波器的输入为有用信号和噪音之和,两者均为广义平稳过程且己知他们的二阶统计过程,则根据最小均方误差准则(滤波器的输出信号与期望信号之差的均方值最小)求出最佳线性滤波器的参数,称之为Wiener滤波器。同时还发现,在一定条件下,这些最佳滤波器与Wiener滤波器是等价的。然而,由于输入过程取决于外界的信号、干扰环境,这种环境的统计特性常常是未知的、变化的,因而不能满足上述两个要求,设计不出最佳滤波器。这就促使人们开始研究自适应滤波器。自适应滤波器由可编程滤波器(滤波部分)和自适应算法两部分组成。可编程滤波器是参数可变的滤波器,自适应算法对其参数进行控制以实现最佳工作。自适应滤波器的参数随着输入信号的变化而变化,因而是非线性和时变的。 2. 自适应滤波器的基础理论 所谓自适应滤波,就是利用前一时刻已获得的滤波器参数等结果,自动地调节现时刻的滤波器参数,以适应信号和噪声未知的或随时间变化的统计特性,从而实现最优滤波。所谓“最优”是以一定的准则来衡量的,最常用的两种准则是最小均方误差准则和最小二乘准则。最小均方误差准则是使误差的均方值最小,它包含了输入数据的统计特性,准则将在下面章节中讨论;最小二乘准则是使误差的平方和最小。 自适应滤波器由数字结构、自适应处理器和自适应算法三部分组成。数字结构是指自适应滤波器中各组成部分之间的联系。自适应处理器是前面介绍的数字滤波器(FIR或IIR),所不同的是,这里的数字滤波器是参数可变的。自适应算法则用来控制数字滤波器参数的变化。 自适应滤波器可以从不同的角度进行分类,按其自适应算法可以分为LMS自适应滤波

MATLAB课程设计(自适应中值滤波)

信息工程系课程设计报告 课程MATLAB课程设计 专业通信工程 班级 2级本科二班 学生姓名1 景学号114 学生姓名2 学号1414 学生姓名3 王学号6 学生姓名4 学号31 学生姓名4 学号02 二〇一四年十二月

目录 目录 (1) 摘要: (2) 关键词: (2) 1.算法描述 (2) 1.1 噪声点 (3) 1.2 窗口尺寸选择 (3) 1.3求滤波窗口内中值,并替换像素点。 (3) 2程序实现 (4) 2.1准备和描述 (4) 2.2扩大窗口、确定窗口 (5) 2.3 确定最大、最小值和中值 (6) 2.4中值替换像素点、输出图像 (7) 实验结果 (9) 参考文献 (9)

摘要:通过本次课程设计,主要训练和培养学生综合应用所学MATLAB课程的自适应中值的相关知识,独立学习自适应中值滤波的原理及处理方式。学会扩大窗口并找到其区域内的中值、最小值、以及最大值,然后用中值代替像素点。通过自主学习和查阅资料来了解程序的编写及改进,并用MATLAB进行仿真。 关键词:自适应中值滤波灰度值椒盐噪声像素点.

1.算法描述 1.1 噪声点 脉冲噪声是图像处理中常见的一类,中值滤波器对消除脉冲噪声非常有效。噪声脉冲可以是正的(盐点),也可以是负的(胡椒点),所以也称这种噪声为“椒盐噪声”。椒盐噪声一般总表现为图像局部区域的最大值或最小值,并且受污染像素的位置是随机分布的,正负噪声点出现的概率通常相等。图像噪声点往往对应于局部区域的极值。 1.2窗口尺寸选择 滤波窗口尺寸的选择影响滤波效果,大尺寸窗口滤波能力强,但细节保持能力较弱;小尺寸窗口能保持图像大量细节但其滤波性能较低。根据噪声密度的大小自适应地选择滤波窗口可以缓和滤波性能与细节保持之间的矛盾,同时也增加了算法的时间复杂度。从形状看来窗口方向要沿着边缘和细节的方向,不能穿过它们也不能把它们和周围相差很大的像素包含在同一窗口中否则边缘和细节会被周围像素模糊。 1.3求滤波窗口内中值,并替换像素点。 设f ij为点(i,j)的灰度,A i,j为当前工作窗口,f min、f max 和f med分别为A i,j中的灰度最小值、灰度最大值和灰度中值, A

有源低通滤波器设计报告要点

课程设计(论文)说明书 题目:有源低通滤波器 院(系):信息与通信学院 专业:通信工程 学生姓名: 学号: 指导教师: 职称: 2010年 12 月 19 日

摘要 低通滤波器是一个通过低频信号而衰减或抑制高频信号的部件。理想滤波器电路的频响在通带内应具有一定幅值和线性相移,而在阻带内其幅值应为零。有源滤波器是指由放大电路及RC网络构成的滤波器电路,它实际上是一种具有特定频率响应的放大器。滤波器的阶数越高,幅频特性衰减的速率越快,但RC网络节数越多,元件参数计算越繁琐,电路的调试越困难。根据指标,本次设计选用二阶有源低通滤波器。 关键词:低通滤波器;集成运放UA741;RC网络 Abstract Low-pass filter is a component which can only pass the low frequency signal and attenuation or inhibit the high frequency signal . Ideal frequency response of the filter circuit in the pass band should have a certain amplitude and linear phase shift, and amplitude of the resistance band to be zero. Active filter is composed of the RC network and the amplifier, it actually has a specific frequency response of the amplifier. Higher the order of the filter, the rate of amplitude-frequency characteristic decay faster, but more the number of RC network section, the more complicated calculation of device parameters, circuit debugging more difficult. According to indicators ,second-order active low-pass filter is used in this design . Key words:Low-pass filter;Integrated operational amplifier UA741;RC network,

自适应滤波器的原理与设计

实验二 自适应滤波信号 一、实验目的: 1.利用自适应LMS 算法实现FIR 最佳维纳滤波器。 2.观察影响自适应LMS算法收敛性,收敛速度以及失调量的各种因素,领会自适应信号处理方法的优缺点。 3.通过实现AR 模型参数的自适应估计,了解自适应信号处理方法的应用。 二、实验原理及方法 自适应滤波是一种自适应最小均方误差算法(LMS ),这种算法不像维纳滤波器需要事先知道输入和输出信号的自相关和互相关矩阵,它所得到的观察值 ,滤波器等价于自动“学习”所需要的相关函数,从而调整FIR 滤波器的权系数,并最终使之收敛于最佳值,即维纳解。 )(n y 下面是自适应FIR 维纳滤波器的LMS 算法公式: (2-1) )()()(0 ^ ^ m n y n h n x M m m -=∑= (2-2) ^ )()()(n x n x n e -=M m m n y n e n h n h m m ?=-?+=+,1) ()(2)()1(^ ^ μ (2-3) 其中FIR 滤波器共有M+1个权系数,表示FIR 滤波器第m 个权系数在第n 步的估计值。 ),0)((^ M m n h m ?=因此,给定初始值)M ,0(),0(?=m h m ,每得到一个样本,可以递归得到一组新的滤波器权系数,只要步长)(n y μ满足 max 1 0λμ< < (2-4) 其中max λ为矩阵R 的最大特征值,当∞→n 时,)M ,0(),0(?=m h m 收敛于维纳解。

现在我们首先考察只有一个权系数h 的滤波器,如图2.1所示。假如信号由下式确定: )(n y )()()(y n w n s n += (2-5) )()(n hx n s = (2-6) 其中h 为标量常数,与互不相关,我们希望利用和得到 )(n x )(n w )(n y )(n x )(n s 图1 利用公式(2-1),(2-2),(2-3),我们可以得到下面的自适应估计算法: (2-7) )()()(^ ^n x n h n s = (2-8) )())()()((2)()1(^ ^ ^ n x n x n h n y n h n h -+=+μ其框图如图所示。 图2 选择的初始值为,对式2-8取数学期望可得 ^)(n h ^ )0(h (2-9) ))0(()21(])([^ ^ h h R h n h E n --+=μ其中

word完整版自适应滤波器原理 带图带总结word版推荐文档

第二章自适应滤波器原理 2.1 基本原理 2.1.1 自适应滤波器的发展 在解决线性滤波问题的统计方法中,通常假设已知有用信号及其附加噪声的某些统计参数(例如,均值和自相关函数) ,而且需要设计含噪数据作为其输入的线性滤波器,使得根据某种统计准则噪声对滤波器的影响最小。实现该滤波器优化问题的一个有用方法是使误差信号(定义为期望响应与滤波器实际输出之差)的均方值最小化。对于平稳输入,通常采用所谓维纳滤波器( Wiener filter) 的解决方案。该滤波器在均方误差意义上使最优的。误差信号均方值相对于滤波器可调参数的曲线通常称为误差性能曲面。该曲面的极小点即为维纳解。 维纳滤波器不适合于应对信号和/或噪声非平稳问题。在这种情况下,必须假设最优滤波器为时变形式。对于这个更加困难的问题,十分成功的一个解决方案使采用卡尔曼滤波器 (Kalman filter )。该滤波器在各种工程应用中式一个强有力的系统。 维纳滤波器的设计要求所要处理的数据统计方面的先验知识。只有当输入数据的统计特性与滤波器设计所依赖的某一先验知识匹配时,该滤波器才是最优的。当这个信息完全未知时,就不可能设计维纳滤波器,或者该设计不再是最优的。而且维纳滤波器的参数是固定的。 在这种情况下,可采用的一个直接方法是“估计和插入过程”。该过程包含两个步骤,首先是“估计”有关信号的统计参数,然后将所得到的结果“插入( plug into)”非递归公式以计算滤波器参数。对于实时运算,该过程的缺点是要求特别精心制作,而且要求价格昂贵的硬件。为了消除这个限制,可采用自适应滤波器(adaptive filter)。采用这样一种系统,意味着滤波器是自设计的,即自适应滤波器依靠递归算法进行其计算,这样使它有可能在无法获得有关信号特征完整知识的环境下,玩完满地完成滤波运算。该算法将从某些预先确定的初始条件集出发,这些初始条件代表了人们所知道的上述环境的任何一种情况。我们还发现,在平稳环境下,该运算经一些成功迭代后收敛于某种统计意义上的最优维纳解。在非平稳环境下,该算法提供了一种跟踪能力,即跟踪输入数据统计特性随时间的变化,只要这种变化时足够缓慢的。 40年代,N.维纳用最小均方原则设计最佳线性滤波器,用来处理平稳随机

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 滤波器。

自适应滤波器的设计与实现毕业论文

自适应滤波器的设计与实现毕业论文 目录 第一章前言 (1) 1.1 自适应滤波器简介 (1) 1.2 选题背景及研究意义 (1) 1.3 国外研究发展现状 (2) 第二章自适应滤波器的基础理论 (4) 2.1 滤波器概述 (4) 2.1.1 滤波器简介 (4) 2.1.2 滤波器分类 (4) 2.1.3 数字滤波器概述 (4) 2.2 自适应滤波器基本理论 (7) 2.3 自适应滤波器的结构 (9) 第三章自适应滤波器递归最小二乘算法 (11) 3.1 递归最小二乘算法 (11) 3.1.1 递归最小二乘算法简介 (11) 3.1.2 正则方程 (11) 3.1.3 加权因子和正则化 (16) 3.1.4 递归计算 (18) 3.2递归最小二乘(RLS)算法的性能分析 (22) 第四章基于MATLAB自适应滤波器仿真 (23) 4.1 正弦波去噪实验 (23) 4.2 滤波器正则化参数的确定 (28) 4.2.1 高信噪比 (28) 4.2.2 低信噪比 (31) 4.2.3 结论 (33) 4.3 输入信号不同对滤波效果的影响 (33)

4.3.1 输入信号为周期信号 (33) 4.3.2 输入信号为非周期信号 (38) 第五章结论与展望 (44) 5.1 结论 (44) 5.2 对进一步研究的展望 (44) 参考文献 (45) 致谢 (46) 附录 (46) 声明 (58)

第一章前言 1.1自适应滤波器简介 自适应滤波器属于现代滤波的畴,它是40年代发展起来的自适应信号处理领域的一个重要应用,自适应信号处理主要是研究结构可变或可调整的系统,可以通过自身与外界的接触来改善自身对信号处理的性能,通常这类系统是时变的非线性系统,可以自动适应信号传输的环境和要求,无须详细的知道信号的结构和实际知识,无须精确设计处理系统本身。 自适应系统的非线性特性主要是由系统对不同的信号环境实现自身参数的调整来确定的。自适应系统的时变特性主要是由其自适应响应或自适应学习过程来确定的,当自适应过程结束和系统不再进行时,有一类自适应系统可成为线性系统,并称为线性自适应系统,因为这类系统便于设计且易于数学处理,所以实际应用广泛。本文研究的自适应滤波器就是这类滤波器。 自适应滤波器是相对固定滤波器而言的,固定滤波器属于经典滤波器,它滤波的频率是固定的,自适应滤波器的频率则是自动适应输入信号而变化的,所以其适用围更广。在没有任何信号和噪声的先验知识的条件下,自适应滤波器利用前一时刻已获得的滤波器参数来自动调节现时刻的滤波器参数,以适应信号和噪声未知或随机变化的统计特性,从而实现最优滤波。1.2选题背景及研究意义 伴随着移动通信事业的飞速发展,自适应滤波技术应用的围也日益扩大。早在20世纪40年代,就对平稳随机信号建立了维纳滤波理论。根据有用信号和干扰噪声的统计特性(自相关函数或功率谱),用线性最小均方误差估计准则设计的最佳滤波器,称为维纳滤波器。这种滤波器能最大程度地滤除干扰噪声,提取有用信号。但是,当输入信号的统计特性偏离设计条件,则它就不是最佳的了,这在实际应用中受到了限制。到60年代初,由于空间技术的发展,出现了卡尔曼滤波理论,即利用状态变量模型对非平稳、多输入多输出随机序列作最优估计。现在,卡尔曼滤波器已成功地应用到许多领域,它既可对平稳的和非平稳的随机信号作线性最佳滤波,也可作非线性滤波。实质上,维纳滤波器是卡尔曼滤波器的一个特例。 在设计卡尔曼滤波器时,必须知道产生输入过程的系统的状态方程和测量方程,即要求对信号和噪声的统计特性有先验知识,但在实际中,往往难以预知这些统计特性,因此实现不了真正的最佳滤波。 Widrow.B等于1967年提出的自适应滤波理论,可使自适应滤波系统的参数自动地调整而

低通滤波器的设计

低通滤波器的设计 模拟滤波器在各种预处理电路中几乎是必不可少的,已成为生物医学仪器中的基本单元电路。有源滤波器实质上是有源选频电路,它的功能是允许指定频段的信号通过,而将其余频段上的信号加以抑制或使其急剧衰减。各种生物信号的低噪声放大,都是首先严格限定在所包含的频谱范围之内。 最常用的全极点滤波器有巴特沃斯滤波器和切比雪夫滤波器。就靠近ω=0处的幅频特性而言,巴特沃斯滤波器比切比雪夫滤波器平直,即在频率的低端巴特沃斯滤波器幅频特性更接近理想情况。但在接近截止频率和在阻带内,巴特沃斯滤波器则较切比雪夫滤波器差得多。本设计中要保证低频信号不被衰减,而对高频要求不高,因此选择了巴特沃斯滤波器。巴特沃思滤波电路(又叫最平幅度滤波电路)是最简单也是最常用的滤波电路,这种滤波电路对幅频响应的要求是:在小于截止频率ωc。的范围内,具有最平幅度响应,而在ω>ωc。后,幅频响应迅速下降。 因为本设计中要保证低频信号不被衰减,而对高频要求不高,所以选择 二阶滤波器即可。本系统采用二阶Butterworth低通滤波器,截止频率f H=100HZ,其电路原理图如1: 图1 低通滤波器图 根据matlab软件算得该设计适合二阶低通滤波器,FSF=628选Z=10000,则

Z R R FSF Z ?=?=的归一值的归一值 C C 3.2脉象信号的的前置放大 由于人体信号的频率和幅度都比较低,很容易受到空间电磁波以及人体其它生理信号的干扰,因此在对其进行变换、分析、存储、记录之前,应该进行一些预处理,以保证测量结果的准确性。因此需要对信号进行放大,“放大”在信号预处理中是第一位的。根据所测参数和所用传感器的不同,放大电路也不同。用于测量生物电位的放大器称为生物电放大器,生物电放大器比一般放大器有更严格的要求。 在本研究中放在传感器后面的电路就是前置放大电路,由于从传感器取得的信号很微弱,且混杂了一些其他的干扰信号。因此前置放大电路的主要功能是,滤除一些共模干扰信号,同时进行一定的放大。该电路由4部分构成:并联型双运放仪器放大器,阻容耦合电路,由集成仪用放大器构成的后继放大器和共模信号取样电路。并联型双运放仪器放大器的优点是不需要精密的匹配电阻,理论上它的共模抑制比为无穷大,且与其外围电阻的匹配程度无关。集成仪用放大器将由并联型双运放仪器放大器输出的双端差动信号转变为单端输出信号,并采用阻容耦合电路隔离直流信号,可以使集成仪用放大器取得较高的差模增益,从而得到很高的共模抑制比。共模取样驱动电路由两个等值电阻和一只由运放构成的跟随器构成,能够使共模信号不经阻容耦合电路的分压直接加在集成放大器的输入端,避免了由于阻容耦合电路的不匹配而降低电路整体的共模抑制比。此电路中也采用了右腿驱动电路来抑制位移电流的影响。前置放大电路参数选择:此部分总的增益取为1000,其中并联型双运放仪器放大器的增益为5,集成仪用放大器的增益为200。具体设计电路如图2所示

自适应中值滤波器matlab实现

将下面代码直接贴入matlab中,并将读入图像修改成自己机子上的,就可以运行了。可以按照“%%”顺序分步来运行 %% function 自适应中值滤波器 %%%%%%%%%%%%%%% %实现两个功能: %1.对高密度的椒盐噪声有好的滤除效果; %2.滤波时减少对图像的模糊; %%%%%%%%%%%%%%% %%%%%%%%%%%%%%% %原理: %1.椒盐噪声概率越大,滤波器窗口需越大。故若滤波器窗口随噪声概率自适应变化,才能有好的滤除效果 %2.为减少对图像的模糊,需在得出原图像值并非椒盐噪声点时,保留原图像值不变; %3.椒盐噪声点的特点:该点的值为该点领域上的最大或最小;%%%%%%%%%%%%%%% %%%%%%%%%%%%%%% %步骤(得到图像中某点(x,y)(即窗口中心点)的值的步骤): %1.设定一个起始窗口,以及窗口的最大尺寸; %2.(此步用于确定窗口大小)对窗口内像素排序,判断中值是否是噪声点,若不是,继续第3步,若是,转到第5步; %3.判断中心点是否是噪声点,若不是,则输出该点的值(即图像中该点的原值不变);若是,则输出中值; %4.窗口尺寸增大,若新窗口尺寸小于设定好的最大值,重复第2步,若大于,则滤波器输出前一个窗口的中值; %%%%%%%%%%%%%%% %%%%%%%%%%%%%%% %参数说明:

%被噪声污染的图像(即退化图像也即待处理图像):Inoise %滤波器输出图像:Imf %起始窗口尺寸:nmin*nmin(只取奇数),窗口尺寸最大值:nmax*nmax %图像大小:Im*In %窗口内图像的最大值Smax,中值Smed,最小值Smin %%%%%%%%%%%%%%%%%%%% %% clear clf %% 读入图像I I=imread('e:/photo/cat.jpg'); %转化为灰度图Ig Ig=rgb2gray(I); %被密度为0.2的椒盐噪声污染的图像Inoise Inoise=imnoise(Ig,'salt & pepper',0.2); %或者是被方差为0.2的高斯噪声污染的图像Inoise %Inoise=imnoise(Ig,'gaussian',0.2); %显示原图的灰度图Ig和噪声图像Inoise subplot(2,2,1),imshow(Ig);xlabel('a.原始灰度图像'); subplot(2,2,2),imshow(Inoise);xlabel('b.被噪声污染的图像'); %% 定义参数 %获取图像尺寸:Im,In [Im,In]=size(Inoise); %起始窗口尺寸:nmin*nmin(窗口尺寸始终取奇数) nmin=3; %最大窗口尺寸:nmax*nmax nmax=9; %定义复原后的图像Imf Imf=Inoise; %为了处理到图像的边界点,需将图像扩充

自适应滤波器的设计(终极版)

目录 摘要…………………..………………………………………………………..….............I 第1章绪论....................................................................................................................错误!未定义书签。 1.1引言……………………………………………...…..…………...……………...错误!未定义书签。 1.2课题研究意义和目的 (1) 1.3国内外研究发展状况 (2) 1.4本文研究思路与主要工作 (4) 第2章自适应滤波器理论基础 (5) 2.1自适应滤波器简介 (5) 2.2自适应滤波器的原理 (5) 2.3自适应滤波算法 (7) 2.4TMS320VC5402的简介 (8) 第3章总体方案设计 (10) 3.1无限冲激响应(IIR)滤波器 (10) 3.2有限冲激响应(FIR)滤波器 (11) 3.3电路设计 (11) 4基于软件设计及仿真 (17) 4.3 DSP的理论基础 (17) 4.4自适应滤波算法的DSP实现 (18) 5总结 (21) 参考文献 (22) 致谢 (23) 附录自适应滤波源代码 (24)

第1章绪论 1.1引言 随着微电子技术和计算机技术的迅速发展,具备了实现自适应滤波器技术的各种软硬件条件,有关自适应滤波器的新算法、新理论和新的实施方法不断涌现,对自适应滤波的稳定性、收敛速度和跟踪特性的研究也不断深入,这一切使该技术越来越成熟,并且在系统辨识、通信均衡、回波抵消、谱线增强、噪声抑制、系统模拟语音信号处理、生物医学电子等方面都获得了广泛应用口。自适应滤波器实现的复杂性通常用它所需的乘法次数和阶数来衡量,而DSP强大的数据吞吐量和数据处理能力使得自适应滤波器的实现更容易。目前绝大多数的自适应滤波器应用是基于最新发展的DSP 来设计的. 滤波技术是信号处理中的一种基本方法和技术,尤其数字滤波技术使用广泛,数字滤波理论的研究及其产品的开发一直受到很多国家的重视。从总的来说滤波可分为经典滤波和现代滤波。经典滤波要求已知信号和噪声的统计特性,如维纳滤波和卡尔曼滤波。现代滤波则不要求己知信号和噪声的统计特性,如自适应滤波。自适应滤波的原理就是利用前一时刻己获得的滤波参数等结果,自动地调节现时刻的滤波参数,从而达到最优化滤波。自适应滤波具有很强的自学习、自跟踪能力,适用于平稳和非平稳随机信号的检测和估计。自适应滤波一般包括3个模块:滤波结构、性能判据和自适应算法。其中,自适应滤波算法一直是人们的研究热点,包括线性自适应算法和非线性自适应算法,非线性自适应算法具有更强的信号处理能力,但计算比较复杂,实际应用最多的仍然是线性自适应滤波算法。线性自适应滤波算法的种类很多,有LMS自适应滤波算法、R路自适应滤波算法、变换域自适应滤波算法、仿射投影算法、共扼梯度算法等。 1.2课题研究意义和目的 自适应滤波理论与技术是现代信号处理技术的重要组成部分,对复杂信号的处理具有独特的功能,对自适应滤波算法的研究是当今自适应信号处理中最为活跃的研究课题之一。自适应滤波器与普通滤波器不同,它的冲激响应或滤波参数是随外部环境的变化而变化的,经过一段自动调节的收敛时间达到最佳滤波的要求。自适应滤波器本身有一个重要的自适应算法,这个算法可以根据输入、输出及原参量信号按照一定准则修改滤波参量,以使它本身能有效的跟踪外部环境的变化。因此,自适应数字系

相关文档
最新文档