matlab由图像得到复振幅程序,基于Matlab的空间滤波实验的计算机仿真

35

基于Matlab 的空间滤波实验的计算机仿真

张奇辉，王 洪，蓝发超

(华南理工大学物理科学与技术学院，广东 广州 510640)

摘 要：利用阿贝-波特实验装置和空间滤波系统，从改变频谱入手改造一幅光学图形进行光学信息处理。在

此基础上，通过Matlab 环境编写程序完成阿贝-波特实验的物理模型的构建并进行计算机模拟实验。

关键词：计算机模拟；Matlab ；空间滤波

中图分类号：TP391.9 文献标识码：A 文章编号：1003-7551(2008)01-0035-04

1 引言

在工程设计领域中，人们通过对研究对象建立模型，用计算机程序实现系统的运行和得到运行结果，寻找出最优方案，然后再予以物理实现，这就是计算机仿真科学。在计算机日益普及的今天，计算机仿真技术作为虚拟实验手段已经成为计算机应用的一个重要分支。它是继理论分析和实物实验之后，认识客观规律性的新型手段。作为科学计算软件，Matlab 的特点是使用方便、输入便捷、运算功能齐全，而且有大量的函数可供使用。因此本文提出基于Matlab 软件，通过在频谱面上设置滤波器对空间频谱的处理，实现对阿贝-波特

实验装置和空间滤波系统的模拟。为了实现仿真实验操作的方便，本文设计出了图形用户可操作界面(GUI )

。 2 空间滤波原理

根据阿贝成像原理，相干光学成像过程可分为两步：第一步称为分频过程，即从物平面到光源的共轭像平面或曰频谱面，由输入的物作为衍射屏对照射光波产生夫琅和费衍射；第二步称为合频或频谱综合过程，即从频谱面到输入物的共轭像平面，被分解的频谱成分经进一步的衍射后再次叠加形成输入物的共轭像。按照傅里叶变换理论，两步成像过程实际上是光学系统对携带输入物信息的二维光场的复振幅分布进行的两次傅里叶变换过程。

以图1所示4f 成像系统为例，此时输入平面O(即物平面)位于透镜1L 的前焦平面，输出平面I(即像平面)位于透镜1L 的后焦平面。透镜1L 和2L 分别起分频(傅里叶变换)和合频(逆傅里叶变换)作用。设输入图像的复振幅分布为),(y x g ，透镜1L 后焦平面T(即频谱面)上的复振幅分布为),(ηξG ，按照傅里叶光学理论，当1L 的孔径无限大时，函数),(ηξG 即等于),(y x g 的傅里叶变换，而),(y x g 为),(ηξG 的傅里叶逆变换，即

(,)(,)exp i2()d d x y x y G f f g x y f x f y x y π∞

?∞

??=?+??∫∫ (1) (,)(,)exp[i2()]d d g x y G x y ξηπμνμν∞∞

?∞?∞=+∫∫ (2)

其中/f μξλ=，/f νηλ=，表示光场(,)G ξη的空间频率。设(',')g x y 为透镜2L 后焦平面I(输出平面)上的复振幅分布，同样，当2L 的孔径无限大时，(',')g x y 就等于的傅里叶变换：

(',')(,)exp[i2('')]d d g x y G x y ξηπμνμν∞∞?∞?∞=

+∫∫ (3)

可以得 (',')(,)g x y g x y ∝?? (4)

即输出图像是输入图像的倒置，且在几何上相似。

现如果在频谱平面T 上插入一衍射屏，并设其复振幅透射率函数为(,)T ξη，则透过衍射屏的光场复振幅分布为：

'(,)(,)(,)G G T ξηξηξη= (5)

经过透镜2L 变换后，在像平面上便得到经过衍射屏调制后的输出图像的光场复振幅分布：

111'(','){'(,)}{(,)}*{(,)}g x y F G F G F T ξηξηξη???== (6) * 收稿日期：2008-01-10