图像频率域低通滤波处理程序设计

1 基本原理简介

1.1 MATLAB 简介

MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品

语言和交互式环境，主要包括MATLAB和Simulink两大部分。

MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程

大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。

创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处

MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并

己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。

1.2傅立叶变换基本原理

傅里叶变换能将满足一定条件的某个函数表示成三角函数（正弦和/或余弦函数）或者它们的积分的线性组合。在不同的研究领域，傅里叶变换具有多种不同的变体形式，如连续傅里叶变换和离散傅里叶变换。最初傅里叶分析是作为热过程的解析分析的工具被提出的[1]。傅里叶变换属于谐波分析。傅里叶变换的逆变换容易求出，而且形式与正变换非常类似。正弦基函数是微分运算的本征函数，从而使得线性微分方程的求解可以转化为常系数的代数方程的求解。在线性时不变的物理系统内，频率是个不变的性质，从而系统对于复杂激励的响应可以通过组合其对不同频率正弦信号的响应来获取。卷积定理指出：傅里叶变换可以化复杂的卷积运算为简单的乘积运算，从而提供了计算卷积的一种简单手段。离散形式的傅里叶变换可以利用数字计算机快速的实现（其算法称为快速傅里叶变换算法（FFT））。

MATLAB中提供的变换函数

（1） fft2:用于计算二维快速傅立叶变换，语句格式：

B=fft2(I,m,n) 按指定的点数计算m,返回矩阵B的大小为m×n,不写默认为原图像大小。（2）fftn:用于计算n维快速傅立叶变换

（3）fftshift:用于将变换后的图像频谱中心从矩阵的原点移到矩阵地中心，语法格式：B=fftshift(I)

(4)ifft2:用于计算图像的二维傅立叶反变换，语法格式

B=ifft2(i)

(5)ifftn:用于计算n维傅立叶变换

快速卷积实验：傅立叶变换一个重要特性是可以实现快速卷积。

设A为M×N矩阵，B为P×Q的矩阵，快速卷积方法如下：

*对A和B补0，使其大小都为(M+P-1)×(N+Q-1)

*利用fft2对矩阵A和B进行二维变换

*将两个FFT结果相乘，利用ifft2对得到的乘积进行傅立叶反变换

1.3滤波电路原理及基本滤波器电路介绍

滤波的基本概念

滤波是信号处理中的一个重要概念。滤波分经典滤波和现代滤波。

经典滤波的概念，是根据富立叶分析和变换提出的一个工程概念。根据高等数学理论，任何一个满足一定条件的信号，都可以被看成是由无限个正弦波叠加而成。换句话说，就是工程信号是不同频率的正弦波线性叠加而成的，组成信号的不同频率的正弦波叫做信号的频率成分或叫做谐波成分。只允许一定频率范围内的信号成分正常通过，而阻止另一部分频率成分通过的电路，叫做经典滤波器或滤波电路。

实际上，任何一个电子系统都具有自己的频带宽度（对信号最高频率的限制），频率特性反映出了电子系统的这个基本特点。而滤波器，则是根据电路参数对电路频带宽度的影响而设计出来的工程应用电路。

用模拟电子电路对模拟信号进行滤波，其基本原理就是利用电路的频率特性实现对信号中频率成分的选择。根据频率滤波时，是把信号看成是由不同频率正弦波叠加而成的模拟信号，通过选择不同的频率成分来实现信号滤波。

当允许信号中较高频率的成分通过滤波器时，这种滤波器叫做高通滤波器。

当允许信号中较低频率的成分通过滤波器时，这种滤波器叫做低通滤波器。

当只允许信号中某个频率范围内的成分通过滤波器时，这种滤波器叫做带通滤波器。理想滤波器的行为特性通常用幅度-频率特性图描述，也叫做滤波器电路的幅频特性。理想滤波器的幅频特性如图所示。图中，w1和w2叫做滤波器的截止频率。

滤波器频率响应特性的幅频特性图

对于滤波器，增益幅度不为零的频率范围叫做通频带，简称通带，增益幅度为零的频率范围叫做阻带。例如对于LP，从-w1当w1之间，叫做LP的通带，其他频率部分叫做阻带。通带所表示的是能够通过滤波器而不会产生衰减的信号频率成分，阻带所表示的是被滤波器衰减掉的信号频率成分。通带内信号所获得的增益，叫做通带增益，阻带中信号所得到的衰减，叫做阻带衰减。在工程实际中，一般使用dB作为滤波器的幅度增益单位。

低通滤波器的基本电路特点是，只允许低于截止频率的信号通过。二阶低通

Butterworth滤波电路如下：下图是用运算放大器设计的二阶低通Butterworth滤波电路。

2 课程设计的基本步骤

2.1 检查图片格式

info=imfinfo('3.jpg')

info =

Filename: '3.jpg'

FileModDate: '24-Jun-2012 15:00:32'

FileSize: 19289

Format: 'jpg'

FormatVersion: ''

Width: 256

Height: 256

BitDepth: 24

ColorType: 'truecolor'

FormatSignature: ''

NumberOfSamples: 3

CodingMethod: 'Huffman'