下载文档原格式
巴特沃斯有源低通滤波器的设计摘要随着社会科学技术的飞速发展，各种科技产品在人类社会中随处可见，极大的丰富了人们的日常生活。
物联设备、可穿戴设备以及虚拟仪器产品在各种应用和消费场合变得极为普遍。
就目前而言，在几乎所有的电子产品中，各种增益、带宽以及高性能的滤波器都发挥着至关重要的作用，例如可穿戴设备的语音信号输入系统中，运用高性能的低通滤波器进行语音信号的降噪、滤波、回声消除，来提高系统的音质和语音识别精准度等。
本论文通过对各种低通滤波器的通频带、增益和截止频率的分析，采用通频带最大扁平度技术（巴特沃斯技术）来设计实现四阶高性能低通滤波器，通过Multisum仿真软件，验证了设计的正确性。
在这基础上，本文还对如何提高该滤波器的响应速度进行了研究，提出了一种有效的提高响应速度的方案，并通过仿真软件得以验证。
这在低通滤波器的理论以及实际工程应用中，都具有非常重要的意义。
关键词：有源低通滤波器，巴特沃斯，运算放大器Design of Butterworth Active Low Pass FilterABSTRACTWith the rapid development of social science and technology, various technological products can be seen everywhere in human society, which greatly enriches people's daily lives. IoT devices, wearable devices, and virtual instrument products have become extremely common in various applications and consumer occasions. For now, in almost all electronic products, various gains, bandwidths, and high-performance filters play a vital role. For example, in the voice signal input system of wearable devices, the use of high-performance low-pass The filter performs noise reduction, filtering, and echo cancellation of the speech signal to improve the sound quality of the system and the accuracy of speech recognition.In this paper, through the analysis of the passband, gain and cutoff frequency of various low-pass filters, the maximum flatness of the passband technology (Butterworth technology) is used to design and implement a fourth-order high-performance low-pass filter, through Multisum simulation software To verify the correctness of the design. On this basis, this paper also studies how to improve the response speed of the filter, and puts forward an effective scheme to improve the response speed, which is verified by simulation software. This is of great significance in the theory of low-pass filters and in practical engineering applications.KEYWORDS：active low-pass filter，butterworth，amplifier1绪论1.1 引言在近现代的科技发展中，滤波器作为一种必不可少的组成成分，在仪器仪表、智能控制、计算机科学、通信技术、电子应用技术和现代信号处理等领域有着十分重要的作用。
用MATLAB 设计巴特沃斯低通滤波器1 巴特沃斯低通滤波器的特性一个理想低通滤波器的幅频特性如图3-80的阴影部分所示。
为了实现这个理想低通特性，需要在从0~ωC 的整个频带内增强增益，在ω>ωC 增益要降到0。
实际上，理想滤波器是不可能实现的。
图3-78是实际滤波器的幅频特性。
但是实际滤波器的特性愈接近理想特性愈好，巴特沃斯（Butterworth ）滤波器就是解决这个问题的方法之一。
巴特沃斯滤波器以巴特沃斯函数来近似滤波器的系统函数，巴特沃斯的低通模平方函数为：221|()|1,2,,1(/)NC H j N j j ωωω==+ (3-138)式中以C ω是滤波器的电压-3dB 点或半功率点。
不同阶次的巴特沃斯滤波器特性如图3-79(a)所示。
4阶巴特沃斯滤波器的极点分布如图3-79(b)所示。
巴特沃斯滤波器幅频响应有以下特点：• 最大平坦性：在0=ω附近一段范围内是非常平直的，它以原点的最大平坦性来逼近理想低通滤波器。
• 通带、阻带下降的单调性。
这种滤波器具有良好的相频特性。
• 3dB 的不变性：随着N 的增加，频带边缘下降越陡峭，越接近理想特性。
但不管N 是多少，幅频特性都通过-3dB 点。
• 极点配置在半径为ωC 的圆上，并且均匀分布。
左半平面上的N 个极点是)(s H 的极点，右半平面上的N 个极点是)(s H -的极点。
2 巴特沃斯低通滤波器的实现为使巴特沃斯滤波器实用，我们必须能够实现它。
一个较好的方法是将巴特沃斯滤波器函数化成若干二阶节级联，其中每一节实现一对共轭复极点。
通过将极点以共轭复数的形式配对，对所有的每一个二阶节都具有实系数。
1图3-78 低通滤波器的幅频特性图3-80所示运算放大器电路为实现一对共轭极点提供了很好的方法。
电路的系统函数为202202121121122121)(1)11(1)(ωωω++=+++=s Qs C C R R s C R C R s C C R R s H (3-139)式中，ω0是S 平面原点与极点之间的距离，Q 被称为电路的“品质因数”，它提供了对响应峰值尖锐程度的一种度量。
信息工程学院实验报告课程名称：数字图像处理Array实验项目名称：实验四图像增强实验时间：班级：姓名：学号：一、实验目的1.了解图像增强的目的及意义，加深对图像增强的感性认识，巩固所学理论知识。
2. 掌握图像空域增强算法的基本原理。
3. 掌握图像空域增强的实际应用及MATLAB实现。
4. 掌握频域滤波的概念及方法。
5. 熟练掌握频域空间的各类滤波器。
6.掌握怎样利用傅立叶变换进行频域滤波。
7. 掌握图像频域增强增强的实际应用及MATLAB实现。
二、实验步骤及结果分析1. 基于幂次变换的图像增强程序代码：clear all;close all;I{1}=double(imread('fig534b.tif'));I{1}=I{1}/255;figure,subplot(2,4,1);imshow(I{1},[]);hold onI{2}=double(imread('room.tif'));I{2}=I{2}/255;subplot(2,4,5);imshow(I{2},[]);hold onfor m=1:2Index=0;for lemta=[0.5 5]Index=Index+1;F{m}{Index}=I{m}.^lemta;subplot(2,4,(m-1)*4+Index+1),imshow(F{m}{Index},[])endend执行结果：图1 幂次变换增强结果实验结果分析：由实验结果可知，当r<1时，黑色区域被扩展，变的清晰；当r>1时，黑色区域被压缩，变的几乎不可见。
2.直方图规定化处理程序代码：clear allclcclose all%0.读图像I=double(imread('lena.tiff'));subplot(2,4,1);imshow(I,[]);title('原图')N=32;Hist_image=hist(I(:),N);Hist_image=Hist_image/sum(Hist_image);Hist_image_cumulation=cumsum(Hist_image);%累计直方图subplot(245);stem(0:N-1,Hist_image);title('原直方图');%1.设计目标直方图Index=0:N-1;%正态分布直方图Hist{1}=exp(-(Index-N/2).^2/N);Hist{1}=Hist{1}/sum(Hist{1});Hist_cumulation{1}=cumsum(Hist{1});subplot(242);stem([0:N-1],Hist{1});title('规定化直方图1');%倒三角形状直方图Hist{2}=abs(2*N-1-2*Index);Hist{2}=Hist{2}/sum(Hist{2});Hist_cumulation{2}=cumsum(Hist{2});subplot(246);stem(0:N-1,Hist{2});title('规定化直方图2');%2. 规定化处理Project{1}=zeros(N);Project{2}=zeros(N);Hist_result{1}=zeros(N);Hist_result{2}=zeros(N);for m=1:2Image=I;%SML处理(SML,Single Mapping Law单映射规则for k=1:NTemp=abs(Hist_image_cumulation(k)-Hist_cumulation{m});[Temp1,Project{m}(k)]=min(Temp);end%2.2 变换后直方图for k=1:NTemp=find(Project{m}==k);if isempty(Temp)Hist_result{m}(k)=0;elseHist_result{m}(k)=sum(Hist_image(T emp));endendsubplot(2,4,(m-1)*4+3);stem(0:N-1,Hist_result{m}); title(['变换后的直方图',num2str(m)]);%2.3结果图Step=256/N;for K=1:NIndex=find(I>=Step*(k-1)&I<Step*k) ;Image(Index)=Project{m}(k);endsubplot(2,4,(m-1)*4+4),imshow(Imag e,[]);title(['变换后的结果图',num2str(m)]);end执行结果：原图规定化直方图2变换后的直方图1变换后的结果图1变换后的直方图2变换后的结果图2图2 直方图规定化实验结果分析：由实验结果可知，采用直方图规定化技术后，原图的直方图逼近规定化的直方图，从而有相应的变换后的结果图1和变换后的结果图2。
实验报告六频域图像增强的方法实验报告六姓名：学号：班级：实验日期： 2016.5.13 实验成绩：实验题目：频域图像增强的方法一．实验目的（1）熟练掌握频域滤波增强的各类滤波器的原理及实现。
（2）分析不同用途的滤波器对频域滤波增强效果的影响，并分析不同的滤波器截止频率对频域滤波增强效果的影响。
二．实验原理变换最慢的频率分量与图像的平均灰度成正比，当远离变换的原点时，低频对应于图像中变换缓慢的灰度分量，当从原点离开得更远时，较高的频率开始对应图像中越来越快的灰度变换，频率域滤波是通过傅里叶变换在频域上对频谱进行修改后再回到空间域的一种方法，在频域中直流项决定了图像的平均灰度，衰减高频通过低频的低通滤波器会模糊一副图像，而衰减低频通过高频的高通滤波器则会增加尖锐的细节，但会导致图像对比度的降低。
三．实验内容及结果（1）选择图像fig620.jpg，对其进行傅里叶变换，在频率域中实现五种不同半径（截止频率）的butterworth低通滤波器的平滑作用。
显示原始图像和滤波图像。
图 1 不同半径巴特沃斯低通滤波图（2）选择图像fig620.jpg，对其进行傅里叶变换，在频率域中实现五种不同半径（截止频率）的butterworth高通滤波的锐化效果，显示原始图像和滤波图像。
图2 不同半径巴特沃斯高通滤波图四．结果分析(1)观察图1，可以发现巴特沃斯低通滤波器半径越小，图像越模糊，但图像的背景亮度大小和原图像别无二致，这是因为低通滤波器实现的是滤除高频分量，保留低频分量的功能，所以半径越小，通过的低频分量越少，所以越模糊，但不论半径多小，它的整体亮度不变，这是由于决定图像平均灰度的直流分量处于图像中点（经过fftshift平移后），它一直是通过的，同时观察变量区的原图像傅里叶变换后的数据矩阵发现，最大数据小于并接近100，所以第五个滤波器设定为100半径，但是发现，滤波后图像的小a还是有一定的模糊，这是因为巴特沃斯滤波器不是理想滤波器，在截止频率处存在一定的过度带，所以小尺寸的物体可能会有模糊。
二维巴特沃斯滤波器1. 简介二维巴特沃斯滤波器是一种常用的图像处理方法，用于对图像进行频域滤波。
它基于巴特沃斯滤波器的原理，在频域中对图像进行平滑或增强。
本文将详细介绍二维巴特沃斯滤波器的原理、实现步骤以及应用案例。
2. 巴特沃斯滤波器原理巴特沃斯滤波器是一种频率域滤波器，通过调整截止频率和阶数来控制信号的频率响应。
其传输函数可以表示为：H (u,v )=11+(D (u,v )D 0)2n其中，D (u,v ) 是图像中每个点到中心点的距离，D 0 是截止频率，n 是阶数。
当 n 取不同值时，巴特沃斯滤波器可以实现不同程度的平滑或增强效果。
当 n >1 时，增加了阶数可以使得低频信号更加平坦；当 n <1 时，减小了阶数可以使得低频信号更加突出。
3. 实现步骤二维巴特沃斯滤波器的实现步骤如下：步骤 1：读取图像首先，需要从文件中读取待处理的图像。
可以使用 C 语言中的图像处理库，如 OpenCV ，来实现图像读取功能。
#include <opencv2/opencv.hpp>int main() {// 读取图像cv::Mat image = cv::imread("input.jpg", cv::IMREAD_GRAYSCALE);// 其他处理步骤...return 0;}步骤 2：进行傅里叶变换将读取的图像进行傅里叶变换，得到频域表示。
可以使用 OpenCV 提供的函数dft 来实现傅里叶变换。
#include <opencv2/opencv.hpp>int main() {// 读取图像cv::Mat image = cv::imread("input.jpg", cv::IMREAD_GRAYSCALE);// 进行傅里叶变换cv::Mat frequencyDomain;cv::dft(image, frequencyDomain, cv::DFT_COMPLEX_OUTPUT);// 其他处理步骤...return 0;}步骤 3：生成巴特沃斯滤波器根据巴特沃斯滤波器的传输函数公式，可以生成巴特沃斯滤波器的频域表示。
巴特沃斯算法
巴特沃斯滤波器是一种经典的频域滤波器，由英国工程师斯蒂芬·巴特沃斯于1930年提出。
它通过调整滤波器的截止频率和阶数，可以实现对不同频率信号的增强或抑制。
在信号处理领域，巴特沃斯低通滤波算法是一种常用的数字滤波算法，能够有效地去除信号中高频成分，保留低频成分，常用于音频处理、图像处理等领域。
此外，根据具体的使用场景和需求，可以在巴特沃斯滤波器算法的基础上进行一些延伸和改进，例如自适应巴特沃斯滤波可以根据信号或图像的局部特征自动调整截止频率或阶数，以提高滤波效果；多尺度巴特沃斯滤波使用不同尺度的巴特沃斯滤波器来处理具有多个频率成分的信号或图像；巴特沃斯带通/带阻滤波设计滤波器以在指定频率范围内增强或抑制信号或图像的频率成分等。
使用巴特沃斯滤波器算法通常包括以下步骤：确定滤波类型，根据需求选择巴特沃斯滤波器的类型，如低通、高通、带通或带阻等；设计滤波器，确定滤波器的截止频率和阶数，并使用相应的库函数或工具进行滤波器设计；滤波操作，将待处理的信号或图像与设计好的巴特沃斯滤波器进行卷积或滤波操作；可选的后处理，根据需要，可以对滤波结果进行后处理，如调整对比度、增强细节等。
C++ 3阶巴特沃斯低通滤波算法在信号处理领域，巴特沃斯低通滤波算法是一种常用的数字滤波算法，它能够有效地去除信号中高频成分，保留低频成分，常用于音频处理、图像处理等领域。
C++作为一种高效的编程语言，能够很好地支持这一算法的实现。
本文将结合C++语言，深入探讨3阶巴特沃斯低通滤波算法的原理、实现和应用。
1. 巴特沃斯低通滤波算法概述巴特沃斯低通滤波器是一种能够通过滤波器将信号中高频成分抑制、低频成分保留的数字滤波器。
其传输函数具有一定的特点，采用巴特沃斯低通滤波器可以实现对信号的平滑处理，去除高频噪声，保留低频信号。
3阶巴特沃斯低通滤波器具有更加优化的特性，能够更好地滤除高频噪声，保留低频信号，因此在实际应用中具有广泛的价值。
2. 3阶巴特沃斯低通滤波算法原理3阶巴特沃斯低通滤波算法是建立在巴特沃斯低通滤波器基础上的改进版本，其核心原理是通过多级滤波器级联的方式，增强滤波效果，同时减少不必要的波纹和相位失真。
其数学模型和传输函数较为复杂，需要通过C++编程语言实现。
3. C++实现3阶巴特沃斯低通滤波算法在C++中实现3阶巴特沃斯低通滤波算法，需要充分利用C++语言的面向对象特性、模板编程等特点。
可以采用模块化的设计思路，将滤波器的设计、参数设置、滤波处理等功能进行封装，从而提高代码的可复用性和可维护性。
C++的性能优势也能够保证算法的高效性。
4. 应用案例分析3阶巴特沃斯低通滤波算法在信号处理领域具有广泛的应用，比如在音频去噪、图像平滑处理、信号恢复等方面均有重要作用。
通过具体的应用案例分析，可以更好地展现算法的效果和实用性，也有助于读者深入理解算法的具体应用场景。
5. 个人观点和总结作为一种经典的数字滤波算法，3阶巴特沃斯低通滤波算法在实际应用中能够发挥重要作用。
在C++语言中实现该算法，既能够充分发挥C++语言的优势，也能够更好地与实际应用结合，为信号处理领域的工程实践提供技术支持。
在未来的发展中，可以进一步优化算法的性能、扩展算法的适用范围，从而更好地满足不同领域的需求。
实验报告实验名称:图像处理姓名:刘强班级:电信1102学号:1404110128实验一图像变换实验——图像点运算、几何变换及正交变换一、实验条件PC机数字图像处理实验教学软件大量样图二、实验目的1、学习使用“数字图像处理实验教学软件系统”,能够进行图像处理方面的简单操作;2、熟悉图像点运算、几何变换及正交变换的基本原理,了解编程实现的具体步骤;3、观察图像的灰度直方图,明确直方图的作用与意义;4、观察图像点运算与几何变换的结果,比较不同参数条件下的变换效果;5、观察图像正交变换的结果,明确图像的空间频率分布情况。
三、实验原理1、图像灰度直方图、点运算与几何变换的基本原理及编程实现步骤图像灰度直方图就是数字图像处理中一个最简单、最有用的工具,它描述了一幅图像的灰度分布情况,为图像的相关处理操作提供了基本信息。
图像点运算就是一种简单而重要的处理技术,它能让用户改变图像数据占据的灰度范围。
点运算可以瞧作就是“从象素到象素”的复制操作,而这种复制操作就是通过灰度变换函数实现的。
如果输入图像为A(x,y),输出图像为B(x,y),则点运算可以表示为:B(x,y)＝f[A(x,y)]其中f(x)被称为灰度变换(Gray Scale Transformation,GST)函数,它描述了输入灰度值与输出灰度值之间的转换关系。
一旦灰度变换函数确定,该点运算就完全确定下来了。
另外,点运算处理将改变图像的灰度直方图分布。
点运算又被称为对比度增强、对比度拉伸或灰度变换。
点运算一般包括灰度的线性变换、阈值变换、窗口变换、灰度拉伸与均衡等。
图像几何变换就是图像的一种基本变换,通常包括图像镜像变换、图像转置、图像平移、图像缩放与图像旋转等,其理论基础主要就是一些矩阵运算,详细原理可以参考有关书籍。
实验系统提供了图像灰度直方图、点运算与几何变换相关内容的文字说明,用户在操作过程中可以参考。
下面以图像点运算中的阈值变换为例给出编程实现的程序流程图,如下:2、图像正交变换的基本原理及编程实现步骤数字图像的处理方法主要有空域法与频域法,点运算与几何变换属于空域法。
基于巴特沃斯低通滤波器的图像增强
陈传峰，朱长仁，宋洪芹
(国防科技大学电子科学与工程学院湖南长沙410073)
l 引言
在进行图像处理的过程中，获取原始图像后，首先需要对图像进行图像预处理，因为在获取图像的过程中，往往会发生图像失真，使所得图像与原图像有某种程度上的差别。
在许多情况下，人们难以确切了解引起图像降质的具体物理过程及其数学模型，但却能估计出使图像降质的一些可能原因，针对这些原因采取简单易行的方法，改善图像质量。
由于噪声、光照等原因，使图像质量不高，为了改善视觉效果或便于人、机器对图像的分析理解，一般都需要对图像进行增强处理，但这个过程并没有统一的标准。
图像增强一般不能增加原图像信息，只能针对一些成像条件，把弱信号突出出来，使一些信息更容易分辨。
图像增强的方法比较多，可以大概分为对比度增强，直方图增强、平滑和锐化［1］4大类，其中，直方图均衡是图像增强的经典方法，因为其有效性和简单易用性已成为图像增强最常用的方法，他又分为全局均衡和局部均衡2种。
全局的直方图均衡［2.3］是对整幅图像进行均衡，使其灰度分布均匀，让每一个灰度等级上的象素个数基本相等，算法简单、计算量小、容易实现，但对图像细节部分增强不够；局部直方图均衡［4.5］则可以增强图像内部细节信息，得到很好的增强效果。
Yin［6］通过对小波分解各个分量进行直方图均衡，然后重构得到处理后图像，S．M．Pi-zer［2］提出自适应直方图均衡算法，这也是一种局部均衡算法。
在对全局直方图均衡和局部直方图均衡详细分析的基础上，本文提出了一种基于Buterworth低通滤波的图像增强方法。
由于Buterworth低通滤波器在抑制噪声的同时，图像边缘模糊程度大大减小，且没有振铃效应。
基于以上特点，用Buterworth低通滤波器将低频分量和高频分量分离，低频分量进行均衡后，再将两部分融合，实现图像的增强。
2直方图均衡化
直方图均衡化［8］是灰度变换的一个重要应用，他是以累计分布函数变换为基础的直方图修正法。
一幅均匀量化的自然图像的灰度直方图通常在低值灰度区间上频率较大，这样的图像的较暗区域中的细节常常看不清楚。
为使图像变清晰，一个自然的想法是使图像的灰度动态范围变大，并且让灰度频率较小的灰度级经变换后其频率变的大一些，从而使象素个数增多，即应使变换后的图像灰度直方图在较大的动态范围内趋于均衡，产生一幅灰度级分布具有均匀概率密度的图像，扩展象素的取值动态范围。
传统的直方图均衡算法(HE)是利用全局图像信息，对给定的数字图像进行对比度增强处理，其处理过程如下：
(1)计算原图像的灰度级直方图；
(2)求原图像各灰度级的累积概率分布函数，并由此构造灰度转换函数；
(3)根据灰度转换函数将原图像所有象素灰度值映射到输出图像，得到增强图像。
这种直方图均衡方法限制了图像中某些局部区域的对比度拉伸力度，使某些细节与背景之间的对比度难以得
到有效增强，甚至出现蜕化。
3局部直方图均衡化算法
全局的直方图均衡使直方图在整个灰度级上的分布趋于均匀，取得一定的增强效果，但忽略了细节的信息甚至出现细节信息的丢失，局部直方图均衡则可以较好地解决这一问题，突出细节信息。
在一幅图像中，图像高频分量对应于图像的细节信息，而低频分量对图像的整体视觉效果有着决定性的影响。
对一幅图像进行分频［9］，得到高频分量和低频分量，然后对低频分量用传统算法进行直方图均衡，高频分量乘以一个加权系数，再将处理得到的结果进行融合，得到最终的处理结果，实现对原图像的增强。
通过实验，证明了其有效性。
整个处理过程分为以下几步：
3．1 分频
对一幅图像来说，能量主要集中在低频分量中，而噪声和图像细节信息主要集中在高频分量中。
在图像处理过程中，将两种分量分离开，进行不同的操作与处理，可以取得更佳的处理效果。
分频的目的就是将输入图像的高频分量和低频分量分开。
进行分频时，既可以选择高通滤波器又可以选择低通滤波器，都能达到分频的目的。
本文实验的过程中使用的是低通滤波器。
在诸多低通滤波器中，Buterworth低通滤波器没有"振铃"现象，而且能够提高图像的细节清晰度，所以，本文算法采用的是Buter-worth低通滤波器，文献[9]采用的是高斯低通滤波器。
Buterworth低通滤波器(BLPF)传递函数：
式中D。
为截止频率，，在D(u，V)=D。
处，；"为阶数，取正整数，用来控制衰减速度。
本文Buterworth低通滤波器截止频率D。
，一50，阶数n一2。
该2阶Buterworth低通滤波器如图1所示。
对一幅输入图像f(r，y)，可将其分解为高频分量fh(x，y)和低频分量f1(x，y)，如下式所示：
将图像f(x，y)通过2阶Buterworth低通滤波器g(x,y)，即f(x,y)与g(x,y)进行卷积，可得到图像的低频分量f(x,y)，如下：
其中(D。
为截止频率)在图像f(x，y)中减去低频分量f1(x，y)就得到高频分量：
3．2低频分量均衡化
图像的能量主要集中在低频分量f1(x，y)中，对图像的整体视觉效果起着决定性的作用，必然要对该分量进行直方图均衡。
对低频分量f(x，y)利用传统算法进行直方图均衡，得到t1(r，y)：
3．3 线性组合得输出图像
噪声存在于高频分量中，如果直接进行直方图均衡，必将使噪声得到增强，影响图像质量。
因此，将fh(x，y)乘上一个加权系数a得到th(x，y)：
a的取值因具体情况而定：若图像受噪声影响比较严重严重，a的取值应小于1，以免放大噪声；若图豫受噪声影响比较小，a的取值应大于l。
输出图像t(T，y)为高频分量fh(x，y)和低频分量f1(T，y)之和：
所以，输出图像f(x，y)为fh(x，y)和HE(f1(x，y))的线性组合，即：
t(x，y)就是希望得到的图像。
通过分频，输人图像被分解为高频分量和低频分量，对两分量进行不同的操作，再经过线性组合得到最终的增强图像，图像清晰度高，细节信息丰富，增强效果明显。
本文算法如下：
(1)对输入图像，f(x，y)进行快速傅里叶变换得F(u，v)；
(2)将F(u，v)作为Buterworth低通滤波器的输入，输出得到低频分量Fl(u，v)，用F(u，v)减去F1(u，v)的高频分量Fh(u，v)，并进行傅里叶反变换，得到fh(x，y)和f1(x，y)；
(3)利用HE算法对低频分量f1(x，y)进行直方图均衡；
(4)高频分量fh(x，y)乘以一个加权系数，再对两分量进行线性组合，融合得到输出图像。
4实验及结果分析本文以lena，flower，frog图像为例进行了实验。
图2中(a)为原始图像；(b)为灰度级较低的退化图像，图像的细节信息和整体的视觉效果都比较差(实验中以退化图像作为输入图像，检验算法的增强效果)；(c)为使用传统HE算法得到的处理结果；(d)为利用本文算法对原图进行增强得到的结果。
4．1 实验结果
实验结果如图2所示：
4．2 结果分析
(1)主观评价
图1(c)中，图像的细节信息并未得到有效增强。
这是因为全局直方图均衡利用的是全局信息，限制了对图像局部区域对比度的拉伸力度，本文算法则利用了局部图像信息，充分再现了图像细节，使细节信息更加突出，
并且视觉效果更好。
(2)客观评价
为了从数据上显示出本文算法的效果，采用了一个评价参量一灰度标准差Sd，定义式如下：
其中，I'表示增强后图像，I表示原图，N为图像象素数。
计算Sd时，首先计算I'中每一象素与I中对应象素的灰度值的差平方；然后求平均；最后再求算术平方根。
分别计算传统算法、本文算法增强后图像与原图间的灰度标准差，数据统计结果如表1所示：
从表1的统计数据来看，对退化图像进行图像增强，本文算法的处理结果较传统算法与原图更为接近，通过增强，更好地将退化图像恢复出原始图像的信息，并且还有了更为丰富的信息。
因而，灰度标准差的值比较大，但数据依然很好地表明了本文算法的有效性。
无论直观地从处理结果上还是从统计数据上来看，与传统图像增强算法相比，本文算法都取得了较好的实验效果，是一种有效的图像增强算法。
5 结语
为使图像处理后的效果得到进一步的提高，文献[10]采用将对比度增强算法与直方图均衡算法级联，使得图中的局部信息和全局信息在增强时都能被利用。
国外文献[2，ll，12]介绍采用多尺度自适应的直方图均衡(Multi-scale Adaptive Histogram Equali-zation)对图像进行处理以增强图像的对比度。
本文采用Buterworth低通滤波器将图像分为高频分量和低频分量，低频分量进行HE算法直方图均衡，然后再将高频分量和低频分量进行融合得到最终处理结果，这样可以突出图像的细节信息，使图像整体视觉效果得到改善，有利于后续图像处理工作的进行。
实验结果证明了本文算法的有效性。
