CN111238779B - Dfdi仪器多普勒干涉条纹对比度分析方法 - Google Patents

Dfdi仪器多普勒干涉条纹对比度分析方法 Download PDF

Info

Publication number
CN111238779B
CN111238779B CN202010065235.4A CN202010065235A CN111238779B CN 111238779 B CN111238779 B CN 111238779B CN 202010065235 A CN202010065235 A CN 202010065235A CN 111238779 B CN111238779 B CN 111238779B
Authority
CN
China
Prior art keywords
interference fringe
fringe contrast
parameters
absorption
instrument
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010065235.4A
Other languages
English (en)
Other versions
CN111238779A (zh
Inventor
吴银花
魏儒义
吴慎将
聂亮
赵保银
张发强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi depule precision optics Co.,Ltd.
Original Assignee
Xian Technological University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Technological University filed Critical Xian Technological University
Priority to CN202010065235.4A priority Critical patent/CN111238779B/zh
Publication of CN111238779A publication Critical patent/CN111238779A/zh
Application granted granted Critical
Publication of CN111238779B publication Critical patent/CN111238779B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M11/00Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
    • G01M11/02Testing optical properties
    • G01M11/0242Testing optical properties by measuring geometrical properties or aberrations
    • G01M11/0271Testing optical properties by measuring geometrical properties or aberrations by using interferometric methods

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明涉及一种DFDI仪器多普勒干涉条纹对比度分析方法。通过DFDI仪器多普勒干涉条纹对比度分析模型,可系统分析仪器各项关键参数对干涉条纹对比度的影响,从而有助于确定各项参数的最佳取值范围。本发明采用的步骤为:1、根据光程差d和光栅分辨率gr,确定中间参数g1和中间参数g2;2、建立DFDI仪器多普勒干涉条纹对比度分析模型;3、利用干涉条纹对比度分析模型,获取当前参数值组合下对应多普勒干涉条纹对比度;4、利用干涉条纹对比度分析模型,获取该任意一个参数在其取值范围内对应的多普勒干涉条纹对比度一维分布;5、利用干涉条纹对比度分析模型,可获取该任意两个参数在其取值范围内对应的多普勒干涉条纹对比度二维分布。

Description

DFDI仪器多普勒干涉条纹对比度分析方法
技术领域:
本发明涉及光学技术领域,尤其涉及一种DFDI仪器多普勒干涉条纹对比度分析方法,即基于DFDI仪器的工作原理,分析多普勒干涉条纹对比度的方法。
背景技术:
色散固定光程差干涉仪(Dispersed Fixed-Delay Interferometer,以下简称为DFDI)是一种新型的视向速度法实现技术,主要用于系外行星的探测。当行星绕其主恒星旋转时,由于行星对恒星的牵引而导致恒星相对于地面观测者的视向速度产生周期性的变化。
DFDI仪器主要由固定延迟干涉仪和中低分辨率后色散器件组成,因此当恒星视向速度产生周期性变化时,恒星吸收谱线经过仪器获取的多普勒干涉条纹会出现周期性的相位变化。因而通过测量DFDI仪器多普勒干涉条纹的周期性相位变化,来探测恒星视向速度的变化,进而探测恒星周围是否存在行星。
DFDI仪器通过有效结合干涉仪和光谱仪的优点,使用中低分辨率色散器件来实现与高精度阶梯光栅等同的探测精度,有效提高仪器透过率、大幅减小仪器体积的同时,降低了仪器对环境影响的敏感度,且具有优良的性价比,是传统阶梯光栅光谱仪的有力补充。
测量多普勒干涉条纹的相位变化,首先应准确提取吸收线对应的干涉条纹,其次应准确解析某一时刻的相位。这就需要获取条纹清晰的二维干涉条纹,也就是说要求所采集的多普勒干涉条纹对比度比较高,否则很难准确定位吸收线位置、且解析的相位精度往往比较低、进而探测的视向速度准确度也比较低。
干涉条纹对比度是DFDI仪器的关键因素。而干涉条纹的对比度与DFDI系统各项关键参数紧密相关,如光栅分辨率、固定光程差、吸收线特性。即在不同光栅分辨率,或在不同的固定光程差,或吸收线特性不同,均会导致干涉条纹的对比度不同。而通过实验逐步调节系统参数来获取对比度较高的干涉条纹,往往需要投入大量的时间和精力,且即使调整出较清晰的干涉条纹,也很难确定当前参数是最佳取值范围,有可能是局部最佳,因而仪器性能无法充分发挥。
发明内容:
本发明提供一种DFDI仪器多普勒干涉条纹对比度分析方法,为了使DFDI仪器生成对比度较高的干涉条纹,通过DFDI仪器多普勒干涉条纹对比度分析模型,可系统分析仪器各项关键参数对干涉条纹对比度的影响,从而有助于确定各项参数的最佳取值范围。
为实现上述目的,本发明采用的技术方案为:DFDI仪器多普勒干涉条纹对比度分析方法,其特征在于:所述的步骤为:
步骤1):根据光程差d和光栅分辨率gr确定中间参数g1,其中Δk0由gr决定,即在吸收线中心波长λa处Δk0=1/(λa·gr);
Figure BDA0002375781020000021
步骤2):根据光程差d,光栅分辨率gr,在吸收线中心波长λa处Δk0=1/(λa·gr),以及恒星吸收线长半高宽Δλa
Figure BDA0002375781020000034
确定中间参数g2;
Figure BDA0002375781020000031
步骤3):根据光栅分辨率gr,恒星吸收线的吸收强度A,吸收线长半高宽Δλa
Figure BDA0002375781020000033
以及步骤1)获取的中间参数g1和步骤2)获取的中间参数g2,建立DFDI仪器多普勒干涉条纹对比度分析模型;
Figure BDA0002375781020000032
步骤4):根据应用需求4个参数光程差d、光栅分辨率gr、吸收线的吸收强度A、吸收线长半高宽Δλa
Figure BDA0002375781020000035
均设定为固定参数值,利用步骤3)建立的干涉条纹对比度分析模型,获取当前参数值组合下对应多普勒干涉条纹对比度;
步骤5):根据应用需求4个参数光程差d、光栅分辨率gr、吸收线的吸收强度A、吸收线长半高宽Δλa
Figure BDA0002375781020000036
中,任意一个参数设定为一定取值范围,其余三个参数设定为固定参数值,利用步骤3)建立的干涉条纹对比度分析模型,获取该任意一个参数在其取值范围内对应的多普勒干涉条纹对比度一维分布;
步骤6):根据应用需求4个参数光程差d、光栅分辨率gr、吸收线的吸收强度A、吸收线长半高宽Δλa
Figure BDA0002375781020000037
中,任意两个参数设定为一定取值范围,其余两个参数设定为固定参数值,利用步骤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
Figure BDA0002375781020000052
确定中间参数g2,如公式(13)所示。
步骤3:根据光栅分辨率gr(在吸收线中心波长λa处Δk0=1/(λa·gr)),恒星吸收线的吸收强度A,吸收线长半高宽Δλa
Figure BDA0002375781020000051
以及步骤1获取的中间参数g1和步骤2获取的中间参数g2,建立DFDI仪器多普勒干涉条纹对比度分析模型,如公式(11)所示。
步骤4:根据应用需求4个参数光程差d、光栅分辨率gr(在吸收线中心波长λa处Δk0=1/(λa·gr))、吸收线的吸收强度A、吸收线长半高宽Δλa
Figure BDA0002375781020000053
均设定为固定参数值,利用步骤3建立的干涉条纹对比度分析模型,可获取当前参数值组合下对应多普勒干涉条纹对比度。
步骤5:根据应用需求4个参数光程差d、光栅分辨率gr(在吸收线中心波长λa处Δk0=1/(λa·gr))、吸收线的吸收强度A、吸收线长半高宽Δλa
Figure BDA0002375781020000054
中,任意一个参数设定为一定取值范围,其余三个参数设定为固定参数值,利用步骤3建立的干涉条纹对比度分析模型,可获取该任意一个参数在其取值范围内对应的多普勒干涉条纹对比度一维分布。
步骤6:根据应用需求4个参数光程差d、光栅分辨率gr(在吸收线中心波长λa处Δk0=1/(λa·gr))、吸收线的吸收强度A、吸收线长半高宽Δλa
Figure BDA0002375781020000064
中,任意两个参数设定为一定取值范围,其余两个参数设定为固定参数值,利用步骤3建立的干涉条纹对比度分析模型,可获取该任意两个参数在其取值范围内对应的多普勒干涉条纹对比度二维分布。
以下对本发明方案的步骤进行更详细地论述:
光的干涉是光的波动性的基本特征,是指两束或多束光在空间相遇时,在重叠区内形成稳定的强弱强度分布的现象。两束光干涉后强度I如公式(1)所示,其中I1和I2是两束光的光强,θ是两束光振动方向间的夹角,
Figure BDA0002375781020000065
是两束光的相位差。
Figure BDA0002375781020000061
当干涉的两束光振动方向相同,即夹角θ为0,又假设两束光的光强I1和I2大小一样,即I1=I2=I0,则(1)式可简化为公式(2)。
Figure BDA0002375781020000062
干涉条纹对比度是用来度量干涉场中某一点附近的条纹的清晰程度。干涉条纹对比度γ定义如公式(3)所示,干涉条纹光强的最大值Imax和最小值Imin之差与最大值Imax和最小值Imin之和的比值。
Figure BDA0002375781020000063
从公式(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/λaa为中心波长),吸收线波数半高宽为
Figure BDA0002375781020000074
(Δλa为波长半高宽),Δk0(Δλ0)是影响艾里斑大小的关键因素,由光栅分辨率gr决定。
S(k)=S1+S2 (4)
Figure BDA0002375781020000071
Figure BDA0002375781020000072
当k=ka时,即在吸收线中心波长位置处,S2又可简化成公式(7),因而此时干涉光谱S也可简化为公式(8)所示。
Figure BDA0002375781020000073
Figure BDA0002375781020000081
公式(8)中,当cos(2πdk)=1时S为最大值Smax,当cos(2πdk)=-1时S为最小值Smin,即Smax和Smin分别如公式(9)和(10)所示。
Figure BDA0002375781020000082
Figure BDA0002375781020000083
将Smax和Smin分别带入公式(3)中Imax和Imin,便可得DFDI仪器多普勒干涉条纹对比度分析模型,如公式(11)所示,其中g1和g2见公式(12)、(13)。
Figure BDA0002375781020000084
Figure BDA0002375781020000085
Figure BDA0002375781020000086
由DFDI仪器多普勒干涉条纹对比度分析模型可知,经干涉、色散后获取的恒星多普勒干涉条纹对比度γcod_absorb值小于1,且主要由4个参数决定,分别是光程差d,光栅分辨率gr(Δk0=k/gr,在吸收线中心波长λa处Δk0=1/(λa·gr)),吸收线的吸收强度A,吸收线长半高宽
Figure BDA0002375781020000087
实施例:
为了更好地说明本发明涉及的DFDI仪器多普勒干涉条纹对比度分析方法,用给定参数值分析DFDI仪器多普勒干涉条纹对比度。
1)假设4个参数均设定为固定参数值,如光程差d为10mm、光栅分辨率gr为15000(在吸收线中心波长λa=800nm处,Δk0=1/(λa·gr)=1/12mm)、吸收线的吸收强度A为0.95、吸收线波长半高宽Δλa为0.02nm
Figure BDA0002375781020000091
通过利用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
Figure BDA0002375781020000092
通过利用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
Figure BDA0002375781020000093
通过利用DFDI仪器多普勒干涉条纹对比度分析模型,可获取当前参数值组合下对应多普勒干涉条纹对比度γcod_absorb值二维分布,如图2所示。
上述实施例仅例示性说明本发明的原理及其功效,以及部分运用的实施例,对于本领域的普通技术人员来说,在不脱离本发明创造构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。

Claims (1)

1.DFDI仪器多普勒干涉条纹对比度分析方法,其特征在于:方法步骤为:
步骤1):根据光程差d和光栅分辨率gr确定中间参数g1,其中Δk0由gr决定,即在吸收线中心波长λa处Δk0=1/(λa·gr);
Figure FDA0003097633440000011
步骤2):根据光程差d,光栅分辨率gr,在吸收线中心波长λa处Δk0=1/(λa·gr),以及恒星吸收线波长半高宽Δλa
Figure FDA0003097633440000012
确定中间参数g2;ΔKa为吸收线波数半高宽;
Figure FDA0003097633440000013
步骤3):根据光栅分辨率gr,恒星吸收线的吸收强度A,吸收线波长半高宽Δλa,以及步骤1)获取的中间参数g1和步骤2)获取的中间参数g2,建立DFDI仪器多普勒干涉条纹对比度分析模型;
Figure FDA0003097633440000014
步骤4):根据应用,需求4个参数:光程差d、光栅分辨率gr、吸收线的吸收强度A、吸收线波长半高宽Δλa,均设定为固定参数值,利用步骤3)建立的干涉条纹对比度分析模型,获取当前参数值组合下对应多普勒干涉条纹对比度;
步骤5):根据应用,需求4个参数光程差d、光栅分辨率gr、吸收线的吸收强度A、吸收线波长半高宽Δλa中,任意一个参数设定为一定取值范围,其余三个参数设定为固定参数值,利用步骤3)建立的干涉条纹对比度分析模型,获取该任意一个参数在其取值范围内对应的多普勒干涉条纹对比度一维分布;
步骤6):根据应用,需求4个参数光程差d、光栅分辨率gr、吸收线的吸收强度A、吸收线波长半高宽Δλa中,任意两个参数设定为一定取值范围,其余两个参数设定为固定参数值,利用步骤3)建立的干涉条纹对比度分析模型,获取该任意两个参数在其取值范围内对应的多普勒干涉条纹对比度二维分布。
CN202010065235.4A 2020-01-20 2020-01-20 Dfdi仪器多普勒干涉条纹对比度分析方法 Active CN111238779B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010065235.4A CN111238779B (zh) 2020-01-20 2020-01-20 Dfdi仪器多普勒干涉条纹对比度分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010065235.4A CN111238779B (zh) 2020-01-20 2020-01-20 Dfdi仪器多普勒干涉条纹对比度分析方法

Publications (2)

Publication Number Publication Date
CN111238779A CN111238779A (zh) 2020-06-05
CN111238779B true CN111238779B (zh) 2021-07-13

Family

ID=70871377

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010065235.4A Active CN111238779B (zh) 2020-01-20 2020-01-20 Dfdi仪器多普勒干涉条纹对比度分析方法

Country Status (1)

Country Link
CN (1) CN111238779B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004093460A (ja) * 2002-09-02 2004-03-25 Olympus Corp レンズ検査方法
CN110399646A (zh) * 2019-07-01 2019-11-01 西安工业大学 一种用于系外行星探测的dfdi仪器模型建立方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5394317B2 (ja) * 2010-05-17 2014-01-22 富士フイルム株式会社 回転対称非球面形状測定装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004093460A (ja) * 2002-09-02 2004-03-25 Olympus Corp レンズ検査方法
CN110399646A (zh) * 2019-07-01 2019-11-01 西安工业大学 一种用于系外行星探测的dfdi仪器模型建立方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
干涉仪中部分偏振光干涉的理论分析;孙柏忠;《长春光学精密机械学院学报》;19931231;第16卷(第2期);全文 *

Also Published As

Publication number Publication date
CN111238779A (zh) 2020-06-05

Similar Documents

Publication Publication Date Title
US7206073B2 (en) Dispersed fourier transform spectrometer
CN103063304B (zh) 色散剪切像面干涉超光谱成像装置及方法
CN104266755B (zh) 一种提高光通量的光谱测量方法及系统
KR101066856B1 (ko) 주파수 차를 구비한 멀티래터럴 간섭계를 포함하는 파면 분석 방법
CN104655290A (zh) 斐索型双波长激光调谐移相干涉测试装置及其测试方法
US6924899B2 (en) System for measuring wavefront tilt in optical systems and method of calibrating wavefront sensors
CN105738073A (zh) 一种在空间频率域进行像素响应函数测量的方法
CN111238779B (zh) Dfdi仪器多普勒干涉条纹对比度分析方法
CN112556842B (zh) 一种双波段高光谱分辨率闪电高速成像仪
Cassaing et al. Optimized fringe tracker for the VLTI/PRIMA instrument
CN110399646A (zh) 一种用于系外行星探测的dfdi仪器模型建立方法
CN106526690A (zh) 视向速度测量的系外行星天基高精度探测系统及方法
WO2018146456A1 (en) Compact interferometer
CN102252765B (zh) 一种扩束准直系统波面像差检测装置及方法
CN111238644B (zh) Dfdi仪器干涉光谱的白光干扰去除方法
US10571442B2 (en) Sagnac fourier spectrometer (SAFOS)
Cai et al. Spatial heterodyne spectrometer based on the Mach–Zehnder interferometer
CN103047928B (zh) 一种对移相干涉仪随机误差模态的测试方法
CN109286809B (zh) 一种图像传感器全阵列像素响应函数测量方法
Ko et al. Analysis of high-resolution spectra from a hybrid interferometric/dispersive spectrometer
Semel A proposal for a solar magnetograph
CN109844474B (zh) 用于分析高清晰度电磁波的方法和设备
Baldi et al. Experimental image dataset for validation of the noise-induced bias that affects Digital Image Correlation
Saetiew et al. Real-time assessment of spectrometer alignment using modulation transfer function (MTF) measurement
CN103822716B (zh) 一种高调制度傅立叶变换光谱获取方法及光谱仪

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210910

Address after: 719000 room 3014, Gaoke building, Kaiyuan Avenue, high tech Industrial Park, Yulin City, Shaanxi Province

Patentee after: Shaanxi depule precision optics Co.,Ltd.

Address before: 710032 No. 2 Xuefu Middle Road, Weiyang District, Xi'an City, Shaanxi Province

Patentee before: XI'AN TECHNOLOGICAL University