矩形孔径光学系统的PSF与MTF计算原理及应用

1. 矩形孔径系统的基本概念与光学特性

在光学成像系统中，孔径的形状直接影响着系统的点扩散函数(PSF)和调制传递函数(MTF)。矩形孔径作为一种非圆形孔径，其光学特性与传统的圆形孔径有着显著差异。典型的矩形孔径常见于数字投影仪、扫描系统和某些特殊设计的相机中。

矩形孔径的数学描述可以用二维矩形函数表示：

code复制rect(x/a, y/b) = { 1, |x|≤a/2且|y|≤b/2
                 { 0, 其他情况

其中a和b分别表示矩形在两个正交方向上的边长。当a=b时，矩形退化为正方形，这是一种特殊的矩形孔径。

2. 矩形孔径系统的点扩散函数(PSF)计算

2.1 理论基础与推导过程

矩形孔径系统的PSF可以通过夫琅禾费衍射理论计算得到。根据傅里叶光学原理，出瞳函数的傅里叶变换就是系统的PSF。对于矩形孔径，其PSF表达式为：

code复制PSF(x,y) = |F{rect(x/a,y/b)}|² = (ab)² sinc²(πax/λf) sinc²(πby/λf)

其中λ是光波长，f是系统焦距，sinc函数定义为sinc(x)=sin(x)/x。这个结果表明，矩形孔径的PSF在x和y方向上是可分离的。

2.2 数值计算方法

在实际计算中，我们通常采用离散傅里叶变换(DFT)来数值计算PSF：

1. 首先在空间域创建矩形孔径的离散采样矩阵
2. 对矩阵进行二维FFT运算
3. 取变换结果的模平方得到PSF
4. 进行必要的归一化处理

Python实现示例：

python复制import numpy as np
import matplotlib.pyplot as plt

def rect_psf(a_ratio, b_ratio, size=512):
    # 创建矩形孔径
    aperture = np.zeros((size, size))
    a_start = int(size/2 - size*a_ratio/2)
    a_end = int(size/2 + size*a_ratio/2)
    b_start = int(size/2 - size*b_ratio/2)
    b_end = int(size/2 + size*b