HCIPy中双星成像模拟的关键参数设置解析

HCIPy中双星成像模拟的关键参数设置解析

hcipy A framework for performing optical propagation simulations, meant for high contrast imaging, in Python. 项目地址: https://gitcode.com/gh_mirrors/hc/hcipy

概述

在使用HCIPy光学仿真库进行双星系统成像模拟时，正确设置传播参数至关重要。本文将通过一个实际案例，详细分析在模拟两颗相距20角秒的恒星成像时需要注意的关键技术细节。

问题背景

在光学系统仿真中，我们经常需要模拟双星系统的成像效果。一个常见的需求是模拟两颗角距离为20角秒的恒星在望远镜焦平面上的成像情况。根据理论计算，使用3.910米的焦距时，两颗星在焦平面上的物理距离应为0.379毫米。

关键参数分析

1. 基本光学参数

• 望远镜口径：0.356米
• 中心遮挡比：0.13/0.356
• 焦距：3.910米
• 工作波长：550nm
• 像素尺寸：3.45μm

2. 仿真设置要点

在初始仿真中，用户创建了光瞳网格和遮挡光阑，并设置了两个点源的波前：

pos_star_1 = np.array([np.radians(star_separation/2/3600),0.0])
wf_1 = Wavefront(shimm_aperture * np.exp(2j * np.pi * np.dot(pos_star_1, pupil_grid.coords) / wavelength), wavelength)

3. 传播过程的关键问题

初始仿真结果与理论值不符的核心原因是在创建FraunhoferPropagator时未指定焦距参数。正确的传播器创建方式应为：

propagator = FraunhoferPropagator(pupil_grid, focal_grid, focal_length=focallength_t)

技术原理详解

1. 角秒与弧度转换

20角秒转换为弧度：

20" = 20/3600度 = (20/3600)×(π/180)弧度 ≈ 9.696×10⁻⁵弧度

2. 理论像面距离计算

在像面上，点源的位移Δx与角位移θ的关系为：

Δx = f × tanθ ≈ f × θ (小角度近似)

因此：

Δx = 3.910m × 9.696×10⁻⁵rad ≈ 0.379mm

3. HCIPy传播器工作原理

FraunhoferPropagator实现的是夫琅禾费衍射公式：

U(x,y) ∝ ∫∫ P(ξ,η) exp[-i2π(xξ+yη)/(λf)] dξdη

其中f是系统焦距。如果不指定焦距，HCIPy会默认使用1米的焦距，导致计算结果与预期不符。

完整解决方案

正确的仿真流程应包括：

1. 正确定义光学系统参数
2. 创建光瞳网格和光阑
3. 设置点源位置（转换为弧度）
4. 创建传播器时明确指定焦距
5. 执行传播计算

修正后的传播器创建代码：

propagator = FraunhoferPropagator(pupil_grid, focal_grid, focal_length=focallength_t)

结论

在HCIPy中进行光学传播仿真时，必须明确指定所有关键光学参数，特别是系统焦距。忽略这一参数会导致计算结果与理论预期严重偏离。通过正确设置FraunhoferPropagator的焦距参数，可以确保双星系统的成像模拟结果与理论计算一致。

这一案例也提醒我们，在进行光学仿真时，必须仔细检查每个物理参数的设置，确保它们与实际系统参数完全对应，才能获得可靠的仿真结果。

hcipy A framework for performing optical propagation simulations, meant for high contrast imaging, in Python. 项目地址: https://gitcode.com/gh_mirrors/hc/hcipy