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

1. 项目背景与核心价值

在光学系统设计和成像质量评估中，点扩散函数（PSF）和调制传递函数（MTF）是两个最关键的量化指标。对于常规圆形孔径的光学系统，PSF和MTF的计算方法已经相当成熟。但当遇到矩形孔径这类特殊结构时，传统的计算模型就会遇到明显局限——这正是本项目要解决的核心问题。

我最近在开发一套红外成像系统时，就遇到了这个典型场景。由于探测器冷屏的机械限制，光学系统的出瞳必须设计为矩形。当时市面上主流的光学设计软件（如Zemax、CodeV）对矩形孔径的PSF/MTF计算支持有限，特别是当需要分析像差影响时，现有工具的计算精度和灵活性都难以满足需求。这促使我开发了这套专门针对矩形孔径的高级计算工具。

这套计算方法的价值主要体现在三个方面：

• 首先，它突破了传统圆形孔径的理论限制，可以准确反映矩形孔径带来的特殊光学效应（如非对称衍射模式）
• 其次，它支持像差场的自定义导入，能够分析像差与矩形孔径的耦合影响
• 最后，计算过程充分考虑了采样率和离散化误差，通过自适应算法确保数值精度

2. 理论基础与算法设计

2.1 矩形孔径的衍射特性

与圆形孔径的旋转对称性不同，矩形孔径的衍射图样具有明显的方向性差异。其夫琅禾费衍射场可以表示为两个一维sinc函数的乘积：

code复制PSF(x,y) ∝ sinc²(πa_x x/λz) · sinc²(πa_y y/λz)

其中a_x和a_y分别是矩形孔径在x和y方向的半宽度，λ为波长，z为传播距离。这种可分离特性是本算法的重要基础。

注意：当系统存在像差时，这个简单关系会被打破，需要采用更复杂的数值计算方法。

2.2 算法实现框架

整个计算流程采用模块化设计，主要包含以下核心组件：

1. 孔径建模模块

   • 支持任意长宽比的矩形孔径定义
   • 可叠加渐晕、遮蔽等实际约束条件
   • 提供参数化像差场输入接口

2. 波前传播模块

   • 采用角谱传播理论
   • 自动优化采样间隔以满足Nyquist准则
   • 支持多波长加权计算

3. PSF计算引擎

   • 基于快速傅里叶变换(FFT)实现
   • 包含场曲/畸变等几何校正
   • 提供矢量衍射计算选项

4. MTF推导模块

   • 通过PSF的