DFDI仪器多普勒干涉条纹对比度分析方法与流程

技术领域：

本发明涉及光学技术领域，尤其涉及一种dfdi仪器多普勒干涉条纹对比度分析方法，即基于dfdi仪器的工作原理，分析多普勒干涉条纹对比度的方法。

背景技术：
：

色散固定光程差干涉仪(dispersedfixed-delayinterferometer，以下简称为dfdi)是一种新型的视向速度法实现技术，主要用于系外行星的探测。当行星绕其主恒星旋转时，由于行星对恒星的牵引而导致恒星相对于地面观测者的视向速度产生周期性的变化。

dfdi仪器主要由固定延迟干涉仪和中低分辨率后色散器件组成，因此当恒星视向速度产生周期性变化时，恒星吸收谱线经过仪器获取的多普勒干涉条纹会出现周期性的相位变化。因而通过测量dfdi仪器多普勒干涉条纹的周期性相位变化，来探测恒星视向速度的变化，进而探测恒星周围是否存在行星。

dfdi仪器通过有效结合干涉仪和光谱仪的优点，使用中低分辨率色散器件来实现与高精度阶梯光栅等同的探测精度，有效提高仪器透过率、大幅减小仪器体积的同时，降低了仪器对环境影响的敏感度，且具有优良的性价比，是传统阶梯光栅光谱仪的有力补充。

测量多普勒干涉条纹的相位变化，首先应准确提取吸收线对应的干涉条纹，其次应准确解析某一时刻的相位。这就需要获取条纹清晰的二维干涉条纹，也就是说要求所采集的多普勒干涉条纹对比度比较高，否则很难准确定位吸收线位置、且解析的相位精度往往比较低、进而探测的视向速度准确度也比较低。

干涉条纹对比度是dfdi仪器的关键因素。而干涉条纹的对比度与dfdi系统各项关键参数紧密相关，如光栅分辨率、固定光程差、吸收线特性。即在不同光栅分辨率，或在不同的固定光程差，或吸收线特性不同，均会导致干涉条纹的对比度不同。而通过实验逐步调节系统参数来获取对比度较高的干涉条纹，往往需要投入大量的时间和精力，且即使调整出较清晰的干涉条纹，也很难确定当前参数是最佳取值范围，有可能是局部最佳，因而仪器性能无法充分发挥。

技术实现要素：
：

本发明提供一种dfdi仪器多普勒干涉条纹对比度分析方法，为了使dfdi仪器生成对比度较高的干涉条纹，通过dfdi仪器多普勒干涉条纹对比度分析模型，可系统分析仪器各项关键参数对干涉条纹对比度的影响，从而有助于确定各项参数的最佳取值范围。

为实现上述目的，本发明采用的技术方案为：dfdi仪器多普勒干涉条纹对比度分析方法，其特征在于：所述的步骤为：

步骤1)：根据光程差d和光栅分辨率gr确定中间参数g1，其中δk0由gr决定，即在吸收线中心波长λa处δk0＝1/(λa·gr)；

步骤2)：根据光程差d，光栅分辨率gr，在吸收线中心波长λa处δk0＝1/(λa·gr)，以及恒星吸收线长半高宽δλa确定中间参数g2；

步骤3)：根据光栅分辨率gr，恒星吸收线的吸收强度a，吸收线长半高宽δλa以及步骤1)获取的中间参数g1和步骤2)获取的中间参数g2，建立dfdi仪器多普勒干涉条纹对比度分析模型；

步骤4)：根据应用需求4个参数光程差d、光栅分辨率gr、吸收线的吸收强度a、吸收线长半高宽δλa均设定为固定参数值，利用步骤3)建立的干涉条纹对比度分析模型，获取当前参数值组合下对应多普勒干涉条纹对比度；

步骤5)：根据应用需求4个参数光程差d、光栅分辨率gr、吸收线的吸收强度a、吸收线长半高宽δλa中，任意一个参数设定为一定取值范围，其余三个参数设定为固定参数值，利用步骤3)建立的干涉条纹对比度分析模型，获取该任意一个参数在其取值范围内对应的多普勒干涉条纹对比度一维分布；

步骤6)：根据应用需求4个参数光程差d、光栅分辨率gr、吸收线的吸收强度a、吸收线长半高宽δλa中，任意两个参数设定为一定取值范围，其余两个参数设定为固定参数值，利用步骤3)建立的干涉条纹对比度分析模型，获取该任意两个参数在其取值范围内对应的多普勒干涉条纹对比度二维分布。

本发明根据dfdi仪器工作原理，建立了多普勒干涉条纹对比度分析模型，模型中涉及dfdi仪器的各项关键参数，且可根据应用需求设定较宽参数范围，并通过利用对比度模型的分析，可确定dfdi仪器各项关键参数的最佳取值范围，从而不仅有助于仪器光路搭建实验，也有助于仪器设计。与现有技术相比，本发明具有的优点和效果如下：

1)本发明根据应用需求设定不同的参数值组合，dfdi仪器多普勒干涉条纹对比度分析模型便可给出当前参数值组合下的多普勒干涉条纹对比度值，这种方式相比通过实验调光路、观测干涉条纹便捷、快速、有效；

2)本发明若将其中的任一参数设定为某个取值范围，其余参数设定为固定参数值，dfdi仪器多普勒干涉条纹对比度分析模型将用曲线的形式显示该任一参数在取值范围内对应的干涉条纹对比度分布情况，进而可有效分析这一参数对干涉条纹对比度的影响；

3)本发明若将其中的任意两个参数设定为某个取值范围，其余参数设定为固定参数值，dfdi仪器多普勒干涉条纹对比度分析模型将用二维图形显示这两个参数在其取值范围内对应的干涉条纹对比度分布情况，进而不仅可有效分析参数对干涉条纹对比度的影响，还可分析两个参数之间相互制约的关系。

附图说明：

图1为对比度随光程差变化的一维分布；

图2为对比度随光程差和光栅分辨率变化的二维分布。

具体实施方式：

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

本发明提供的dfdi仪器多普勒干涉条纹对比度分析方法具体步骤如下：

步骤1：根据光程差d和光栅分辨率gr(在吸收线中心波长λa处δk0＝1/(λa·gr))，确定中间参数g1，如公式(12)所示。

步骤2：根据光程差d，光栅分辨率gr(在吸收线中心波长λa处δk0＝1/(λa·gr))，以及恒星吸收线长半高宽δλa确定中间参数g2，如公式(13)所示。

步骤3：根据光栅分辨率gr(在吸收线中心波长λa处δk0＝1/(λa·gr))，恒星吸收线的吸收强度a，吸收线长半高宽δλa以及步骤1获取的中间参数g1和步骤2获取的中间参数g2，建立dfdi仪器多普勒干涉条纹对比度分析模型，如公式(11)所示。

步骤4：根据应用需求4个参数光程差d、光栅分辨率gr(在吸收线中心波长λa处δk0＝1/(λa·gr))、吸收线的吸收强度a、吸收线长半高宽δλa均设定为固定参数值，利用步骤3建立的干涉条纹对比度分析模型，可获取当前参数值组合下对应多普勒干涉条纹对比度。

步骤5：根据应用需求4个参数光程差d、光栅分辨率gr(在吸收线中心波长λa处δk0＝1/(λa·gr))、吸收线的吸收强度a、吸收线长半高宽δλa中，任意一个参数设定为一定取值范围，其余三个参数设定为固定参数值，利用步骤3建立的干涉条纹对比度分析模型，可获取该任意一个参数在其取值范围内对应的多普勒干涉条纹对比度一维分布。

步骤6：根据应用需求4个参数光程差d、光栅分辨率gr(在吸收线中心波长λa处δk0＝1/(λa·gr))、吸收线的吸收强度a、吸收线长半高宽δλa中，任意两个参数设定为一定取值范围，其余两个参数设定为固定参数值，利用步骤3建立的干涉条纹对比度分析模型，可获取该任意两个参数在其取值范围内对应的多普勒干涉条纹对比度二维分布。

以下对本发明方案的步骤进行更详细地论述：

光的干涉是光的波动性的基本特征，是指两束或多束光在空间相遇时，在重叠区内形成稳定的强弱强度分布的现象。两束光干涉后强度i如公式(1)所示，其中i1和i2是两束光的光强，θ是两束光振动方向间的夹角，是两束光的相位差。

当干涉的两束光振动方向相同，即夹角θ为0，又假设两束光的光强i1和i2大小一样，即i1＝i2＝i0，则(1)式可简化为公式(2)。

干涉条纹对比度是用来度量干涉场中某一点附近的条纹的清晰程度。干涉条纹对比度γ定义如公式(3)所示，干涉条纹光强的最大值imax和最小值imin之差与最大值imax和最小值imin之和的比值。

从公式(3)可以看出，当干涉光强的极小值imin＝0时，γ＝1，两光束完全相干，条纹最清晰；当imax＝imin时，γ＝0，两光束完全不相干；当imax≠imin≠0时，0<γ<1，两光束部分相干，条纹清晰度介于上述两种情况之间。

公式(2)表示的是理想的单色光干涉条纹，按公式(3)计算可得出其干涉条纹对比度为1，即此时两光束完全相干。

根据dfdi仪器的工作原理，dfdi仪器是由干涉仪和光谱仪两部分组成，干涉后形成的干涉条纹还要经过光谱仪进行后色散，但由于后色散模块存在光栅模糊效应，因此探测器所接收的干涉条纹不再是公式(2)所示的理想干涉条纹，而是经过光栅的线扩散函数调制后的干涉条纹。即此时的干涉条纹s由均匀连续光的干涉条纹s1和摩尔条纹s2两部分组成，如公式(4)～(6)所示，其中d表示光程差，k表示波数，吸收线的吸收强度为a，吸收线中心波数为ka＝1/λa(λa为中心波长)，吸收线波数半高宽为(δλa为波长半高宽)，δk0(δλ0)是影响艾里斑大小的关键因素，由光栅分辨率gr决定。

s(k)＝s1+s2(4)

当k＝ka时，即在吸收线中心波长位置处，s2又可简化成公式(7)，因而此时干涉光谱s也可简化为公式(8)所示。

公式(8)中，当cos(2πdk)＝1时s为最大值smax，当cos(2πdk)＝-1时s为最小值smin，即smax和smin分别如公式(9)和(10)所示。

将smax和smin分别带入公式(3)中imax和imin，便可得dfdi仪器多普勒干涉条纹对比度分析模型，如公式(11)所示，其中g1和g2见公式(12)、(13)。

由dfdi仪器多普勒干涉条纹对比度分析模型可知，经干涉、色散后获取的恒星多普勒干涉条纹对比度γcod_absorb值小于1，且主要由4个参数决定，分别是光程差d，光栅分辨率gr(δk0＝k/gr，在吸收线中心波长λa处δk0＝1/(λa·gr))，吸收线的吸收强度a，吸收线长半高宽

实施例：

为了更好地说明本发明涉及的dfdi仪器多普勒干涉条纹对比度分析方法，用给定参数值分析dfdi仪器多普勒干涉条纹对比度。

1)假设4个参数均设定为固定参数值，如光程差d为10mm、光栅分辨率gr为15000(在吸收线中心波长λa＝800nm处，δk0＝1/(λa·gr)＝1/12mm)、吸收线的吸收强度a为0.95、吸收线波长半高宽δλa为0.02nm通过利用dfdi仪器多普勒干涉条纹对比度分析模型，可获取当前参数值组合下对应多普勒干涉条纹对比度γcod_absorb值为0.2424。

2)假设4个参数中，光程差d取值范围为0～50mm、光栅分辨率gr为15000(在吸收线中心波长λa＝800nm处，δk0＝1/(λa·gr)＝1/12mm)、吸收线的吸收强度a为0.95、吸收线波长半高宽δλa为0.02nm通过利用dfdi仪器多普勒干涉条纹对比度分析模型，可获取当前参数值组合下对应多普勒干涉条纹对比度γcod_absorb值一维分布，如图1所示。

3)假设4个参数中，光程差d取值范围为0～50mm、光栅分辨率gr取值范围为5000～30000(在吸收线中心波长λa＝800nm处，δk0＝1/(λa·gr)取值范围为1/24mm～1/4mm)、吸收线的吸收强度a为0.95、吸收线波长半高宽δλa为0.02nm通过利用dfdi仪器多普勒干涉条纹对比度分析模型，可获取当前参数值组合下对应多普勒干涉条纹对比度γcod_absorb值二维分布，如图2所示。

上述实施例仅例示性说明本发明的原理及其功效，以及部分运用的实施例，对于本领域的普通技术人员来说，在不脱离本发明创造构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。