Zemax 技术文档

Lumerical & Speos

Ansys 仿真

产品问答

方案 > Zemax 技术文档 > 什么是点扩散函数（PSF）

什么是点扩散函数（PSF）

概要

这篇文章讲述了：

1.什么是点扩散函数？

2.点列图

3.快速傅里叶变换计算的点扩散函数（FFT PSF）

4.惠更斯算法计算的点扩散函数（Huygens PSF）

5.如何使用非序列模式下的透镜和探测器观察惠更斯积分

6.对于序列模式下的透镜组，分析点扩散函数时如何在惠更斯和快速傅里叶变换计算两者间进行选择

本文使用的附件请在以下链接下载：

http://customers.zemax.com/support/knowledgebase/Knowledgebase-Attachments/What-is-a-Point-Spread-Function/PSF_SampleFiles.aspx

什么是点扩散函数（PSF）？

PSF是一个物空间的点光源经过光学系统后的辐射照度分布。望远镜对遥远行星成像就是一个很好的PSF例子：由于行星距离我们非常远以至于我们在任何条件下都可以把它当做一个点。虽然发光光源是一个理想的点，但是像点并不是一个点。这是因为两个主要原因：首先，光学系统的像差会影响像点，使其在有限的区域内扩散；其次，即使光学系统没有像差，衍射效应同样会使像点发生扩散。

Zemax OpticStudio内置了三种计算PSF的方法：几何（无衍射效应）点列图、基于衍射效应的快速傅里叶变换（FFT）PSF和惠更斯（Huygens）PSF。这篇文章将讨论每种算法的理论基础，并给出了一些使用建议。

点列图：几何PSF

点列图是Zemax OpticStudio中最基础的分析功能之一。该功能会加载多条从物空间点光源发出的光线，并对所有光线穿过光学系统进行追迹，然后将以图表的形式显示光线相对于特定参考点的坐标图。

本文中使用的示例系统是一个焦距为50mm、F数为5的抛物面反射镜。物体放置于无穷远处。这是一个简单的牛顿望远镜系统，示例文件名为PSF_Newtonian.ZMX。光学系统如下所示：

两个视场下的点列图如下所示，其中一个视场为轴上视场，另一个为2°的离轴视场：

如图可见，点列图是一堆点的集合，其中每一个点代表一根光线。每根光线间没有任何相互作用。点列图非常高效的展示了望远镜系统的几何效应或光线像差。其中，离轴几何PSF图像（2°视场下点列图）非常形象的展示了系统的慧差和像散。然而对于轴上视场，点列图给出了完美像点的分析结果。这是光学系统成像的正确结果吗？为了回答这一问题，我们需要将点列图的结果即点的分布，与衍射极限相比较。

将艾里斑（Airy Disk）的参考椭圆添加到点列图中是一个快速比对几何像差和衍射效应的方法。打开点列图的设置列表，并勾选显示艾里斑：

现在，点列图将在几何点分布上的相应位置处，显示出艾里斑的直径：

在轴上视场上点的分布远小于艾里斑，而在离轴视场上点的分布则远大于艾里斑。这表明点列图只能有效合理的展示离轴视场的成像结果。想要同时计算轴上和离轴视场正确的PSF，则需要考虑衍射效应带来的影响。

通常来说，当光学系统的像差影响相比衍射效应大很多时，点列图是评价光学系统成像质量非常有效的方法。

FFT PST

快速傅里叶变换（FFT）算法广泛应用于许多光电系统的频率分析中。从概念上讲，FFT将一个空间上的分布分解为频域上的分布。对于傅里叶光学的详细讲解请参考书籍《Introduction to Fourier optics》，在用户手册的“物理光学传播（Physical Optics Propagation）”一章中也有对衍射理论的总结。两份资料都介绍了菲涅尔和夫琅禾费衍射理论。

大多数光学系统都满足基于夫琅禾费衍射理论的FFT PSF算法所必须的简化假设。其主要假设有：

1.F数足够大使得标量衍射理论成立

2.衍射PSF能量显著区域远小于光学系统出瞳到像面的距离

3.出瞳相对入瞳没有明显畸变。这表示入瞳上均匀分布的光线在出瞳上也应是合理的均匀分布

4.高采样率以满足准确模拟PSF的要求

5.像面上的主光线和法线夹角很小

光学系统的FFT PSF是通过如下方法计算所得：一组呈网格状分布的光线从点光源发出射向系统出瞳。OpticStudio会使用幅值和光程差计算每一根光线出瞳处波前网格上对应点的复振幅。经过适当缩放后该网格的快速傅里叶变换经过平方得到PSF的实部值。如果计算是多波长的，则PSF将进行累加。想要计算序列模式系统的FFT PST，请在Zemax OpticStudio的主菜单中选择分析选项卡（Analysis）>点扩散函数（PSF）> FFT PSF。示例文件牛顿望远镜系统的轴上视场的FFT PSF，如下图所示，其分析设置参数稍后我们会做讨论。

注意图中相似的艾里斑波形。牛顿望远系统的轴上视场像差为零，如上图所示这也是符合预期的结果。

为了生成如上图所示的图像，FFT PST的设置菜单应该如下图所示：

这里的采样是指对射向系统入瞳处的网格分布光线的采样率。Zemax OpticStudio会自动设置显示网格的尺寸比采样网格大一倍，并将超出入瞳区域的数据用零代替。因此PSF的分析结果显示的点要比网格采样点多4倍。如果系统像差很小，则我们感兴趣的区域将集中在图表的中心。为了减少显示这些近似等于零的数据，我们选择显示网格为较小的尺寸。

在分析参数设置中可以使用多种网格尺寸来显示同样的PSF数据，可以设置如下图所示参数：

显示网格的尺寸为128x128，视场编号为2，类型为对数，显示为伪彩色图。其分析结果如下图所示：

惠更斯PSF

从概念上讲，惠更斯PSF通过将点列图上每一根光线转换为一个小的平面波。回想一下，光线可以模拟一个小区域的平面波，因此光线在各向同性介质中与波前垂直。由光线模拟的平面波拥有振幅、相位和传播方向这三个参数。因此像面上的总辐照度分布可以考虑为所有追迹光线所代表的平面波的相干叠加。对所有光线的积分可以直接得到衍射PSF。

Zemax OpticStudio中大多数衍射效应是基于标量衍射理论的（F数不会很小），但惠更斯算法在“使用偏振（Use Polarization）”开启时可以考虑电场的矢量特性。所有基于惠更斯算法的分析都全面考虑了偏振向量和偏振相位引入的像差。这一算法通过分别计算偏振电场的Ex,Ey,Ez分量，然后进行非相干叠加。和计算其他相位差一样，惠更斯算法会计算电场中每一个正交分量引入的偏振相位差。

实际上所有光学系统都满足进行惠更斯PSF计算的前提假设：

采样率足够高以准确模拟PSF

惠更斯PSF不是基于FFT进行计算的。虽然计算速度方面惠更斯PSF通常要比FFT PSF慢，但是在不满足FFT PSF前提假设的情况下惠更斯PSF要更加准确。使用FFT PSF假设有问题从而需要使用惠更斯PSF的情况有：

1.像面与主光线的垂直面存在明显倾斜

2.出瞳相对于入瞳存在明显的畸变

光学系统的惠更斯PSF的计算如下。从点光源出发追迹网格分布的光线并传播到像面。使用光线的振幅、坐标、方向余弦和光程差，计算每根光线在像空间网格上每点处的平面波的复振幅。将像空间每个网格点上所有光线的复振幅相干叠加。则像空间每个网格点上的强度为叠加复振幅的平方。如果计算中考虑多个波长，则对不同波长的PSF结果进行非相干叠加。

如果您想要计算序列模式下系统的惠更斯PSF，在Zemax OpticStudio主菜单中选择分析选项卡（Analysis）> PSF > 惠更斯PSF（Huygens PSF）。惠更斯PSF同样可以计算非序列元件（NSC）系统，我们将在后文中详细介绍。需要注意的是FFT PSF不能应用于非序列系统。

由用户定义的惠更斯PSF的关键设置参数是光瞳采样、像面采样和像面采样间距。这些参数可以在惠更斯PSF分析参数设置表中输入。打开参数设置表并如下图输入参数：

像面采样间距为像面上每个网格点之间的距离，其单位为毫米。计算PSF的总体区域的尺寸是像面采样和像面采样间距之积。

如下图所示为示例牛顿望远镜系统轴上视场的惠更斯PSF：

离轴视场如下图所示（将视场编号改为2）：

光线数量和像点越多，PSF的分辨率和计算精度越高。但这会带来计算时间的增加。

观察惠更斯积分

观察积分过程的一种方法是一次观察一根光线的相干叠加结果。这可以通过在Zemax rOpticStudio的非序列元件功能中设置相干探测器来完成。示例文件为HPSF_Integration.ZMX。

该示例文件包括一个椭圆光源、一个单透镜和一个矩形探测器物体。光源在一个圆形区域内产生多条随机光线。所有光线沿Z轴平行出射。光源模拟一束准直光或无穷远的点光源。光线输出条数设为20，分析光线的条数设为1。这将使系统在同一时刻只分析一条光线，我们会在后文中详细讨论。透镜为一个简单的单透镜，它用来将平行光很好的聚焦在探测器上。探测器属性设置为吸收，像素数量为120x120。

需要注意的是，探测器的第11个参数“PSF波长模式（PSF Wave#）”设为1。

这一设置可以让探测器计算相干惠更斯PSF的积分。每一根射到探测器上的光线都被转换为局部平面波，该平面波照射探测器上的所有像素。随后平面波在每一像素上的相干振幅，叠加到已探测的相干振幅上。这允许我们可根据需求每次只追迹一条光线，所以我们可以看到每一根单独光线的叠加效果。

每次积分一根光线

接下来让我们通过示例观察这一积分过程，打开示例文件HPSF_Integration，选择分析选项卡（Analysis）- 光线追迹（Trace Rays）- 光线追迹控制（Ray Trace Control）。选中自动更新（Auto Update）。点击“追迹（Trace）”。由于光源只定义了一根分析光线，因此一根随机的光线被追击并且探测器更新当前结果。再次点击“追迹”，但注意不先点击“清空探测器（Clear Detectors）”来追击第二根光线。此时两根光线的追迹结果将显示出两个传播方向有一定夹角的平面波的干涉结果，探测器将显示条纹图案。由于光线的选取是随机的，因此条纹图案每次也不相同，因此探测器结果可能与下图并不完全相同。