矩形孔径光学系统PSF与MTF计算原理与实践

1. 项目概述：矩形孔径系统的光学特性分析

在光学成像系统设计中，点扩散函数（PSF）和调制传递函数（MTF）是评估系统性能的两个核心指标。当光学系统的孔径形状从传统的圆形变为矩形时，其成像特性会发生显著变化。这种矩形孔径设计常见于特殊用途的光学仪器、扫描系统和部分空间受限的光学装置中。

矩形孔径系统的PSF表现为二维sinc函数形式，这与圆形孔径的Airy斑图案有本质区别。MTF计算则需要考虑矩形对称性带来的方向性差异——系统在x和y方向上的分辨率表现会因矩形长宽比不同而产生分化。理解这种差异对光学工程师至关重要，特别是在设计条形扫描系统、光刻设备或特殊传感器时。

2. 矩形孔径的光学原理基础

2.1 衍射理论与孔径函数

任何光学系统的成像过程都可以用衍射理论来描述。当光波通过有限大小的孔径时，会发生衍射现象。对于矩形孔径，其透射函数可以表示为：

python复制# 矩形孔径函数示例
def rect_aperture(x, y, a, b):
    return np.where((np.abs(x) <= a/2) & (np.abs(y) <= b/2), 1, 0)

其中a和b分别代表矩形在两个正交方向上的尺寸。这个二元函数在傅里叶光学中起着关键作用，因为它直接决定了系统的衍射特性。

2.2 从孔径到PSF的转换

根据标量衍射理论，光学系统的PSF是孔径函数的傅里叶变换的模平方。对于矩形孔径，这个变换可以得到解析解：

code复制PSF(x,y) ∝ sinc²(πax/λf) * sinc²(πby/λf)

其中λ是波长，f是系统焦距，sinc(x) = sin(x)/x。这个结果表明，矩形孔径的PSF在x和y方向上具有不同的展宽特性，直接影响了系统的分辨率表现。

关键提示：当矩形长宽比a/b增大时，PSF在长边方向会变得更窄，这意味着该方向上的理论分辨率更高，但同时旁瓣效应也会更明显。

3. 高级PSF计算方法与实践

3.1 数值计算实现步骤

实际工程中，我们通常采用离散傅里叶变换(DFT)来计算PSF。以下是典型实现流程：

1. 定义采样网格：
   • 根据奈奎斯特准则确定采样频率
   • 创建N×N的二维坐标网格（通常N取256或512）