搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于光场一阶关联的时域成像

张瑞雪 李洪国 李宗国

基于光场一阶关联的时域成像

张瑞雪, 李洪国, 李宗国
PDF
HTML
导出引用
导出核心图
  • 与通常利用二阶强度关联测量实现时域鬼成像不同, 本文利用时域热光源借助干涉仪通过一阶关联实现时域成像. 基于空域光束的近轴衍射和时域窄带脉冲在色散介质中色散之间的空间-时间二象性, 在时域脉冲响应函数的基础上得到了表征一阶关联时域成像的强度表达式, 分析研究了光源脉冲宽度和相干时间对成像可见度和分辨率的影响. 结果一方面表明基于热光场一阶关联的时域成像在不需要额外色散补偿或消除条件下可以实现时域物体信号的再现, 另一方面表明当光源脉冲宽度一定时, 成像可见度随光源脉冲相干时间的增加而增加, 但是成像分辨率逐渐降低, 其中当光源脉冲宽度约为100 ps, 相干时间约为0.5 ps时, 间隔为20 ps, 宽度为8 ps的时域矩形波型物体的成像质量(兼顾可见度和分辨率)较好. 该结果对于基于热光一阶关联的时域成像在时序信号测量中的应用具有重要意义.
      通信作者: 李宗国, zgli@tjut.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11604243)和天津市自然科学基金(批准号: 16JCQNJC01600)资助的课题.
    [1]

    Padgett M J, Boyd R W 2017 Phil. Trans. R. Soc. A 375 20160233

    [2]

    Pittman T B, Shih Y H, Strekalov D V, Sergienko A V 1995 Phys. Rev. A 52 R3429

    [3]

    Bennink R S, Bentley S J, Boyd R W 2002 Phys. Rev. Lett. 89 113601

    [4]

    Gatti A, Brambilla E, Bache M, Lugiato L A 2004 Phys. Rev. A 70 013802

    [5]

    Cheng J, Han S 2004 Phys. Rev. Lett. 92 093903

    [6]

    Cao D Z, Xiong J, Wang K G 2005 Phys. Rev. A 71 013801

    [7]

    Valencia A, Scarcelli G, D’ Angelo M, Shih Y H 2005 Phys. Rev. Lett. 94 063601

    [8]

    Ferri F, Magatti D, Gatti A, Bache M, Brambilla E, Lugiato L A 2005 Phys. Rev. Lett. 94 183602

    [9]

    Cai Y, Zhu S Y 2005 Phys. Rev. E 71 056607

    [10]

    Zhang D, Zhai Y H, Wu L A, Chen X H 2005 Opt. Lett. 30 2354

    [11]

    Cai Y, Wang F 2007 Opt. Lett. 32 205

    [12]

    Liu X F, Chen X H, Yao X R, Yu W K, Zhai G J, Wu L A 2014 Opt. Lett. 39 2314

    [13]

    Sun B, Edgar M P, Bowman R, Vittert L E, Welsh S, Bowman A, Padgett M J 2013 Science 340 844

    [14]

    Bromberg Y, Katz O, Silberberg Y 2009 Phys. Rev. A 79 053840

    [15]

    Shapiro J H 2008 Phys. Rev. A 78 061802

    [16]

    Zhao C Q, Gong W L, Chen M L, Li E R, Wang H, Xu W D, Han S S 2012 Appl. Phys. Lett. 101 141123

    [17]

    Hong Y, Li E R, Gong W L, Han S S 2015 Opt. Express 23 14541

    [18]

    Chen M, Li E, Gong W L, Bo Z, Xu X, Zhao C, Shen X, Xu W, Han S S 2013 Opt. Photonics J. 3 83

    [19]

    Li S, Cropp F, Kabra K, Lane T J, Wetzstein G, Musumeci P, Ratner D 2018 Phys. Rev. Lett. 121 114801

    [20]

    Cheng J 2009 Opt. Express 17 7916

    [21]

    Cheng J, Lin J 2013 Phys. Rev. A 87 043810

    [22]

    Cao D Z, Xiong J, Zhang S H, Lin L F, Gao L, Wang K G 2008 Appl. Phys. Lett. 92 201102

    [23]

    Chan K W C, O’ Sullivan M N, Boyd R W 2010 Opt. Express 18 5562

    [24]

    Zhang D J, Li H G, Zhao Q L, Wang S, Wang H B, Xiong J, Wang K G 2015 Phys. Rev. A 92 013823

    [25]

    Li H G, Zhang D J, Xu D J, Zhao Q L, Wang S, Wang H B, Xiong J, Wang K G 2015 Phys. Rev. A 92 043816

    [26]

    Katz O, Bromberg Y, Silberberg Y 2009 Appl. Phys. Lett. 95 131110

    [27]

    仲亚军, 刘娇, 梁文强, 赵生妹 2015 物理学报 64 014202

    Zhong Y J, Liu J, Liang W Q, Zhao S M 2015 Acta Phys. Sin. 64 014202

    [28]

    Gao C, Wang X, Wang Z, Li Z, Du G, Chang F, Yao Z 2017 Phys. Rev. A 96 023838

    [29]

    Cao D H, Li Q H, Zhuang X C, Ren H, Zhang S H, Song X B 2018 Chin. Phys. B 27 123401

    [30]

    Yang H, Wu H, Wang H B, Cao D H, Zhang S H, Xiong J, Wang K 2018 Phys. Rev. A 98 053853

    [31]

    Salem R, Foster M A, Gaeta A L 2013 Adv. Opt. Photonics 5 274

    [32]

    Foster M A, Salem R, Geraghty D F, Turner-Foster A C, Lipson M, Gaeta A L 2008 Nature 456 81

    [33]

    Schröder J, Wang F, Clarke A, Ryckeboer E, Pelusi M , Roelens M A, Eggleton B J 2010 Opt. Commun. 283 2611

    [34]

    Fridman M, Farsi A, Okawachi Y, Gaeta A L 2012 Nature 481 62

    [35]

    Ryczkowski P, Barbier M, Friberg A T, Dudley J M, Genty G 2016 Nat. Photonics 10 167

    [36]

    Shirai T, Setälä T, Friberg A T 2010 J. Opt. Soc. Am. B 27 2549

    [37]

    Setälä T, Shirai T, Friberg A T 2010 Phys. Rev. A 82 043813

    [38]

    Chen Z, Li H, Li Y, Shi J, Zeng G 2013 Opt. Eng. 52 076103

    [39]

    Gao L, Zhang S H, Xiong J, Gan S, Feng L J, Cao D Z, Wang K G 2009 Phys. Rev. A 80 021806

    [40]

    Vabre L, Dubois A, Boccara A C 2002 Opt. Lett. 27 530

    [41]

    Kolner B H 1994 IEEE J. Quant. Electron. 30 1951

    [42]

    Cai Y, Zhu S 2004 Opt. Lett. 29 2716

    [43]

    Qu L, Bai Y, Nan S, Shen Q, Li H, Fu X 2018 Opt. Laser Technol. 104 197

  • 图 1  基于光场一阶关联的时域成像装置示意图

    Fig. 1.  Schematic illustration for first-order temporal imaging.

    图 2  光源相干时间为0.5 ps, 光源脉冲宽度宽度分别为5, 10, 50和100 ps条件下的关联像

    Fig. 2.  The correlated images with different pulse widths of the source 5, 10, 50 and 100 ps, for the source’s coherence time 0.5 ps

    图 3  光源相干时间Tc为0.5 ps时, 成像可见度随光源脉冲宽度的变化曲线

    Fig. 3.  The imaging visibility as a function of the source’s pulse width for the source’s coherence time 0.5 ps

    图 4  光源脉冲宽度为T0 = 100 ps, 光源相干时间Tc分别为0.5, 1.5, 5和10 ps条件下的关联像

    Fig. 4.  The correlated images with different source’s coherence time 0.5, 1.5, 5 and 10 ps, for the source’s pulse width 100 ps

    图 5  光源脉冲宽度为100 ps时, 成像可见度随光源相干时间的变化曲线

    Fig. 5.  The imaging visibility as a function of the source coherence time for the source’s pulse width 100 ps.

  • [1]

    Padgett M J, Boyd R W 2017 Phil. Trans. R. Soc. A 375 20160233

    [2]

    Pittman T B, Shih Y H, Strekalov D V, Sergienko A V 1995 Phys. Rev. A 52 R3429

    [3]

    Bennink R S, Bentley S J, Boyd R W 2002 Phys. Rev. Lett. 89 113601

    [4]

    Gatti A, Brambilla E, Bache M, Lugiato L A 2004 Phys. Rev. A 70 013802

    [5]

    Cheng J, Han S 2004 Phys. Rev. Lett. 92 093903

    [6]

    Cao D Z, Xiong J, Wang K G 2005 Phys. Rev. A 71 013801

    [7]

    Valencia A, Scarcelli G, D’ Angelo M, Shih Y H 2005 Phys. Rev. Lett. 94 063601

    [8]

    Ferri F, Magatti D, Gatti A, Bache M, Brambilla E, Lugiato L A 2005 Phys. Rev. Lett. 94 183602

    [9]

    Cai Y, Zhu S Y 2005 Phys. Rev. E 71 056607

    [10]

    Zhang D, Zhai Y H, Wu L A, Chen X H 2005 Opt. Lett. 30 2354

    [11]

    Cai Y, Wang F 2007 Opt. Lett. 32 205

    [12]

    Liu X F, Chen X H, Yao X R, Yu W K, Zhai G J, Wu L A 2014 Opt. Lett. 39 2314

    [13]

    Sun B, Edgar M P, Bowman R, Vittert L E, Welsh S, Bowman A, Padgett M J 2013 Science 340 844

    [14]

    Bromberg Y, Katz O, Silberberg Y 2009 Phys. Rev. A 79 053840

    [15]

    Shapiro J H 2008 Phys. Rev. A 78 061802

    [16]

    Zhao C Q, Gong W L, Chen M L, Li E R, Wang H, Xu W D, Han S S 2012 Appl. Phys. Lett. 101 141123

    [17]

    Hong Y, Li E R, Gong W L, Han S S 2015 Opt. Express 23 14541

    [18]

    Chen M, Li E, Gong W L, Bo Z, Xu X, Zhao C, Shen X, Xu W, Han S S 2013 Opt. Photonics J. 3 83

    [19]

    Li S, Cropp F, Kabra K, Lane T J, Wetzstein G, Musumeci P, Ratner D 2018 Phys. Rev. Lett. 121 114801

    [20]

    Cheng J 2009 Opt. Express 17 7916

    [21]

    Cheng J, Lin J 2013 Phys. Rev. A 87 043810

    [22]

    Cao D Z, Xiong J, Zhang S H, Lin L F, Gao L, Wang K G 2008 Appl. Phys. Lett. 92 201102

    [23]

    Chan K W C, O’ Sullivan M N, Boyd R W 2010 Opt. Express 18 5562

    [24]

    Zhang D J, Li H G, Zhao Q L, Wang S, Wang H B, Xiong J, Wang K G 2015 Phys. Rev. A 92 013823

    [25]

    Li H G, Zhang D J, Xu D J, Zhao Q L, Wang S, Wang H B, Xiong J, Wang K G 2015 Phys. Rev. A 92 043816

    [26]

    Katz O, Bromberg Y, Silberberg Y 2009 Appl. Phys. Lett. 95 131110

    [27]

    仲亚军, 刘娇, 梁文强, 赵生妹 2015 物理学报 64 014202

    Zhong Y J, Liu J, Liang W Q, Zhao S M 2015 Acta Phys. Sin. 64 014202

    [28]

    Gao C, Wang X, Wang Z, Li Z, Du G, Chang F, Yao Z 2017 Phys. Rev. A 96 023838

    [29]

    Cao D H, Li Q H, Zhuang X C, Ren H, Zhang S H, Song X B 2018 Chin. Phys. B 27 123401

    [30]

    Yang H, Wu H, Wang H B, Cao D H, Zhang S H, Xiong J, Wang K 2018 Phys. Rev. A 98 053853

    [31]

    Salem R, Foster M A, Gaeta A L 2013 Adv. Opt. Photonics 5 274

    [32]

    Foster M A, Salem R, Geraghty D F, Turner-Foster A C, Lipson M, Gaeta A L 2008 Nature 456 81

    [33]

    Schröder J, Wang F, Clarke A, Ryckeboer E, Pelusi M , Roelens M A, Eggleton B J 2010 Opt. Commun. 283 2611

    [34]

    Fridman M, Farsi A, Okawachi Y, Gaeta A L 2012 Nature 481 62

    [35]

    Ryczkowski P, Barbier M, Friberg A T, Dudley J M, Genty G 2016 Nat. Photonics 10 167

    [36]

    Shirai T, Setälä T, Friberg A T 2010 J. Opt. Soc. Am. B 27 2549

    [37]

    Setälä T, Shirai T, Friberg A T 2010 Phys. Rev. A 82 043813

    [38]

    Chen Z, Li H, Li Y, Shi J, Zeng G 2013 Opt. Eng. 52 076103

    [39]

    Gao L, Zhang S H, Xiong J, Gan S, Feng L J, Cao D Z, Wang K G 2009 Phys. Rev. A 80 021806

    [40]

    Vabre L, Dubois A, Boccara A C 2002 Opt. Lett. 27 530

    [41]

    Kolner B H 1994 IEEE J. Quant. Electron. 30 1951

    [42]

    Cai Y, Zhu S 2004 Opt. Lett. 29 2716

    [43]

    Qu L, Bai Y, Nan S, Shen Q, Li H, Fu X 2018 Opt. Laser Technol. 104 197

  • [1] 刘雪峰, 姚旭日, 李明飞, 俞文凯, 陈希浩, 孙志斌, 吴令安, 翟光杰. 强度涨落在热光鬼成像中的作用. 物理学报, 2013, 62(18): 184205. doi: 10.7498/aps.62.184205
    [2] 李龙珍, 姚旭日, 刘雪峰, 俞文凯, 翟光杰. 基于压缩感知超分辨鬼成像. 物理学报, 2014, 63(22): 224201. doi: 10.7498/aps.63.224201
    [3] 周阳, 张红伟, 钟菲, 郭树旭. 基于自适应阈值方法实现迭代降噪鬼成像. 物理学报, 2018, 67(24): 244201. doi: 10.7498/aps.67.20181240
    [4] 李曙光, 刘晓东, 侯蓝田. 一种晶体光纤基模色散特性的矢量法分析. 物理学报, 2004, 53(6): 1873-1879. doi: 10.7498/aps.53.1873
    [5] 仲亚军, 刘娇, 梁文强, 赵生妹. 针对多散斑图的差分压缩鬼成像方案研究. 物理学报, 2015, 64(1): 014202. doi: 10.7498/aps.64.014202
    [6] 杨磊, 刘楠楠, 李小英. Sagnac光纤环制备并分离简并关联光子对的实验研究. 物理学报, 2016, 65(19): 194202. doi: 10.7498/aps.65.194202
    [7] 孙健明, 于洁, 郭霞生, 章东. 基于分数导数研究高强度聚焦超声的非线性声场. 物理学报, 2013, 62(5): 054301. doi: 10.7498/aps.62.054301
    [8] 乔耀军, 韩庆生, 李蔚. 基于全光时域分数阶傅里叶变换的光脉冲最小损伤传输新方法. 物理学报, 2011, 60(1): 014219. doi: 10.7498/aps.60.014219
    [9] 魏东宾, 赵兴涛, 苑金辉, 孟 佳, 王海云, 周桂耀, 侯蓝田. 一种新型的多包层光子晶体光纤的分析方法. 物理学报, 2008, 57(5): 3011-3015. doi: 10.7498/aps.57.3011
    [10] 陈翔, 张心贲, 祝贤, 程兰, 彭景刚, 戴能利, 李海清, 李进延. 色散补偿光子晶体光纤结构参数对其色散的影响. 物理学报, 2013, 62(4): 044222. doi: 10.7498/aps.62.044222
    [11] 李林栗, 冯国英, 杨浩, 周国瑞, 周昊, 周寿桓, 朱启华, 王建军. 纳米光纤的色散特性及其超连续谱产生. 物理学报, 2009, 58(10): 7005-7011. doi: 10.7498/aps.58.7005
    [12] 李政颖, 孙文丰, 李子墨, 王洪海. 基于色散补偿光纤的高速光纤光栅解调方法. 物理学报, 2015, 64(23): 234207. doi: 10.7498/aps.64.234207
    [13] 耿易星, 李荣凤, 赵研英, 王大辉, 卢海洋, 颜学庆. 色散对双晶交叉偏振滤波输出特性的影响. 物理学报, 2017, 66(4): 040601. doi: 10.7498/aps.66.040601
    [14] 周康, 黎华, 万文坚, 李子平, 曹俊诚. 太赫兹量子级联激光器频率梳的色散. 物理学报, 2019, 68(10): 109501. doi: 10.7498/aps.68.20190217
    [15] 任国斌, 王 智, 娄淑琴, 简水生. 高折射率芯Bragg光纤的色散特性研究. 物理学报, 2004, 53(6): 1862-1867. doi: 10.7498/aps.53.1862
    [16] 李曙光, 刘晓东, 侯蓝田. 光子晶体光纤色散补偿特性的数值研究. 物理学报, 2004, 53(6): 1880-1886. doi: 10.7498/aps.53.1880
    [17] 赵兴涛, 侯蓝田, 刘兆伦, 王 伟, 魏红彦, 马景瑞. 改进的全矢量有效折射率方法分析光子晶体光纤的色散特性. 物理学报, 2007, 56(4): 2275-2280. doi: 10.7498/aps.56.2275
    [18] 甘琛利, 聂志强, 李 岭, 姜 彤, 沈磊剑, 李沛哲, 宋建平, 张彦鹏, 卢克清. 倒V形四能级亚飞秒极化拍的三光子吸收和色散. 物理学报, 2008, 57(1): 243-251. doi: 10.7498/aps.57.243
    [19] 张小民, 王建军, 许党朋, 张锐, 林宏焕, 邓颖, 耿远超, 余晓秋, 黄小东. 色散对高能激光光纤前端FM-AM效应的影响. 物理学报, 2010, 59(3): 1857-1862. doi: 10.7498/aps.59.1857
    [20] 赵岩, 施伟华, 姜跃进. 中心外缺陷对带隙型光子晶体光纤色散特性的影响. 物理学报, 2010, 59(9): 6279-6283. doi: 10.7498/aps.59.6279
  • 引用本文:
    Citation:
计量
  • 文章访问数:  276
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-02-01
  • 修回日期:  2019-03-05
  • 上网日期:  2019-06-05
  • 刊出日期:  2019-05-01

基于光场一阶关联的时域成像

  • 天津理工大学理学院, 天津 300384
  • 通信作者: 李宗国, zgli@tjut.edu.cn
    基金项目: 国家自然科学基金(批准号: 11604243)和天津市自然科学基金(批准号: 16JCQNJC01600)资助的课题.

摘要: 与通常利用二阶强度关联测量实现时域鬼成像不同, 本文利用时域热光源借助干涉仪通过一阶关联实现时域成像. 基于空域光束的近轴衍射和时域窄带脉冲在色散介质中色散之间的空间-时间二象性, 在时域脉冲响应函数的基础上得到了表征一阶关联时域成像的强度表达式, 分析研究了光源脉冲宽度和相干时间对成像可见度和分辨率的影响. 结果一方面表明基于热光场一阶关联的时域成像在不需要额外色散补偿或消除条件下可以实现时域物体信号的再现, 另一方面表明当光源脉冲宽度一定时, 成像可见度随光源脉冲相干时间的增加而增加, 但是成像分辨率逐渐降低, 其中当光源脉冲宽度约为100 ps, 相干时间约为0.5 ps时, 间隔为20 ps, 宽度为8 ps的时域矩形波型物体的成像质量(兼顾可见度和分辨率)较好. 该结果对于基于热光一阶关联的时域成像在时序信号测量中的应用具有重要意义.

English Abstract

    • 鬼成像, 又称关联成像, 是一种基于光的涨落关联获得物体图像的新型成像技术. 具有涨落关联性质的光源发出的光被分束器分成两路, 在一条光路(物光路)中, 透过物体的光被一个不具有空间分辨的桶探测器接收, 在另一条光路(参考光路)中光经过一段自由传播后直接被一个具有空间分辨的探测器接收, 其中单个探测器都不能够得到物体的像, 但是对两个探测器接收到的信号进行关联测量却可重建物体的图像信息[1]. 国内外关于鬼成像的研究逐步走向应用, 从最初的纠缠光源[2]到热光源[3-12]、空间光调制器或投影仪调制光源[13], 从双臂鬼成像到单臂鬼成像、计算鬼成像[14,15], 从室内近距离成像到室外远距离遥感成像应用[16-18]、从光子鬼成像到电子鬼成像[19]. 另外在鬼成像的应用方面还有关于鬼成像质量影响因素的相关研究, 如大气湍流对成像的影响[20,21]、高阶关联鬼成像的可见度及信噪比[22,23]、多波长光源对信噪比的影响[24]、基于光学参量放大的鬼成像[25]、压缩感知鬼成像[26,27]、计算鬼成像质量的优化[28]等. 最近, Cao等[29]研究了基于热光源的分数阶鬼成像, Yang等[30]首次从概率密度角度理论解释了热光的鬼成像, 并进行了实验验证. 利用空域中光的近轴衍射和时域中脉冲信号在色散介质中色散之间的空间-时间二象性, 很多空域系统的时域对应因此得到了快速发展, 比如超快光信号处理[31]、超快示波器[32,33]、时域隐身[34]等.

      最近, Ryczkowski等[35]实现了基于二阶关联的时域鬼成像[36-38], 为利用关联方法进行时域信号的探测开辟了新方向. 此前Gao等[39]基于物理上描述鬼成像的强度关联(二阶场关联)函数可以分解为一阶场关联, 提出并实现了基于热光场一阶关联的空域关联成像. 而目前关于一阶关联的时域成像还没有相关研究, 本文基于空间-时间二象性, 设计并研究基于热光场一阶关联的时域成像, 促进利用光场关联实现时域信号探测的应用. 这种成像模式采用 Mach-Zehnder干涉仪通过物光场和参考光场的关联获得物体的图像信息, 装置类似于基于一阶场关联的全息成像[39], 但区别在于一阶关联全息成像中是记录物体的干涉条纹, 然后通过再现获得物体的图像. 另外这种成像模式与基于一阶场关联的光学相干层析(OCT)成像的相同点是都利用了低相干光源的时间相干性, 但区别在于 OCT成像装置是基于迈克耳孙干涉仪[40]. 本文第2部分给出实现一阶关联时域成像的装置图并进行理论分析研究, 第3部分给出基于一阶关联实现时域物体成像的结果, 并分析光源的脉冲宽度和相干时间对成像可见度和分辨率的影响, 最后进行 总结.

    • 为了便于后面得到描述一阶时域成像的关联函数, 首先分析给出时域脉冲响应函数. 在空域中, 描述近轴光场衍射的波方程为

      $ \frac{{{\partial ^2}E}}{{\partial {x^2}}} + \frac{{{\partial ^2}E}}{{\partial {y^2}}} + 2{\rm i}k\frac{{\partial E}}{{\partial z}} = 0. $

      该波方程的通解形式上可写为

      $ \begin{split} & E\left( {x,y,z} \right) \\ =& \iint{E\left( {{x}_{0}},{{y}_{0}},0 \right)h\left( x,y;{{x}_{0}},{{y}_{0}} \right){\rm d}{{x}_{0}}{\rm d}{{y}_{0}},} \end{split} $

      其中$ E\left( {x_0, y_0, 0} \right) $$ E\left( {x, y, z} \right) $分别表示位置$ \left( {x, y, z} \right) $和初始位置$ \left( {x_0, y_0, 0} \right) $处的光场, $ h\left( {x, y; x_0, y_0} \right) $是描述光学系统的脉冲响应函数, 其表达式为

      $ \begin{split} & h\left( x,y;{{x}_{0}},{{y}_{0}} \right)\\ =& \frac{k}{{\rm i}2{\text{π}} z}\exp \left\{\!\frac{{\rm i}k}{2z}\left[ {{\left( x-{{x}_{0}} \right)}^{2}}+{{\left( y -{{y}_{0}} \right)}^{2}} \right] \right\}. \end{split} $

      上面的解也可以通过菲涅耳衍射理论得到. 这表明近轴波方程与菲涅耳衍射理论存在内在联系. 相应地, 在一阶色散近似条件下[41], 描述窄带光脉冲在色散介质中色散的微分方程为

      $ \frac{\partial^{2}E}{\partial \tau^{2}}+\frac{2{\rm i}}{\beta}\frac{\partial E}{\partial l} = 0, $

      其中$ E $表示与在色散介质中传播的时间$ \tau $和距离$ l $有关的光场分布, $ \beta $是光在色散介质中的群速度色散系数. 根据以上方程, 描述脉冲色散的微分方程与描述近轴光波空间衍射的波方程形式相同. 因此时域脉冲在色散介质中的色散可看作空域中光波近轴衍射的时间对应, 即空间-时间二象性. (4)式的解与近轴衍射波方程解的数学形式相似, 可写为

      $ E\left( {\tau} \right) = \int E\left( {\tau_0} \right)h\left( {\tau, \tau_0} \right){\rm d}\tau_0, $

      其中脉冲响应函数$ h\left( {\tau, \tau_0} \right) $表征在色散介质中系统在$ \tau $时刻对初始时刻$ \tau_0 $输入脉冲的响应, 其形式为[41]

      $ h\left( {\tau ,{\tau _0}} \right) = \sqrt {\frac{1}{{2{\text{π}} {\rm i}\varPhi }}} \exp \left[ {\frac{\rm i}{{2\varPhi }}{{\left( {\tau - {\tau _0}} \right)}^2}} \right], $

      其中$ \varPhi = \beta l $表示长度为$ l $的色散介质的群延迟色散参量. 基于上述二像性, 也可写出其他光学系统的脉冲响应函数. 因此类比空域物体的脉冲响应函数, 时域物体的脉冲响应函数为

      $ h\left( {\tau, \tau_{0}} \right) = P\left( {\tau_0} \right)\delta\left( {\tau-\tau_0} \right), $

      其中$ P\left( {\tau} \right) $是描述物体时序结构的函数, 可以通过电光调制器调制产生. 同理, 时域透镜的脉冲响应函数可写为

      $ h\left( {\tau, \tau_{0}} \right) = \exp\left( {-\frac{{\rm i}\tau^{2}}{2\varPhi_{\rm f}}} \right)\delta\left( {\tau-\tau_{0}} \right), $

      其中$ \varPhi_{\rm f} $表征时域透镜的“焦时间”, 类似于空间透镜的焦距. 时域透镜用于对光脉冲施以二次相位调制, 该调制也可通过电光调制器实现.

    • 采用类似于空域中借助 Mach-Zehnder干涉仪实现一阶关联成像的思想, 本文设计的一阶关联时域成像装置示意图如图1所示. 具有光场涨落关联性质的时序光源$ {\rm TS} $发出的光在色散介质如光纤中传播一段距离后被光纤分束器$ {\rm FBS}_1 $分为强度相等的两束光, 一束光(物光束, $ {\rm Path }\;1 $)在色散介质中传播, 依次经过时域物体$ {\rm TO} $和时域透镜$ {\rm TL}_1 $后到达光纤分束器$ {\rm FBS}_2 $; 另一束光(参考光束, $ {\rm Path }\;2 $)在色散介质中经过时域透镜$ {\rm TL}_2 $后也到达光纤分束器$ {\rm FBS}_2 $, 两束光在光纤分束器$ {\rm FBS}_2 $中汇合后被探测器$ {\rm Detector} $接收, 探测器采集得到的信号通过数据采集卡将信号数据存储至计算机.

      图  1  基于光场一阶关联的时域成像装置示意图

      Figure 1.  Schematic illustration for first-order temporal imaging.

      根据(6)式—(8)式, 并考虑脉冲响应函数的级联性质, 可得物光路$ {\rm Path}\; 1 $和参考光路$ {\rm Path }\;2 $系统的脉冲响应函数分别为:

      $\begin{split}{h_1}\left( {\tau ,{\tau _0}} \right) = &\frac{1}{{2{\text{π}} {\rm i}\sqrt {{\varPhi _0}{\varPhi _{{\rm f_1}}}} }}\int {P\left( {{\tau ^\prime }} \right)}\\&\times\exp \left[ {\frac{\rm i}{{2{\varPhi _0}}}{{\left( {{\tau ^\prime } - {\tau _0}} \right)}^2} - \frac{{{\rm i}\tau {\tau ^\prime }}}{{{\varPhi _{{\rm f_1}}}}}} \right]{\rm d}{\tau ^\prime },\end{split}$

      $ \begin{split} {{h_2}\left( {\tau ,{\tau _0}} \right)}= & \frac{1}{{2{\text{π}} {\rm i}\sqrt {{\varPhi _1}{\varPhi _2}} }}\exp \left( {\frac{{{\rm i}\tau _0^2}}{{2{\varPhi _1}}} + \frac{{{\rm i}{\tau ^2}}}{{2{\varPhi _2}}}} \right)\\ & \times \int {\exp } \left[ {\frac{{{\rm i}\tau _2^{\prime 2}}}{2}\left( { \!-\! \frac{1}{{{\varPhi _{{\rm f_2}}}}} \!+\! \frac{1}{{{\varPhi _1}}} \!+\! \frac{1}{{{\varPhi _2}}}} \right)} \right.\\ &\left. { - {\rm i}\tau _2^\prime \left( {\frac{{{\tau _0}}}{{{\varPhi _1}}} + \frac{\tau }{{{\varPhi _2}}}} \right)} \right]{\rm d}\tau _2^\prime , \end{split} $

      其中$ \varPhi_0 $表示光源与时域物体之间色散介质的群延迟色散参量, $ \varPhi_{\rm f_1} $$ \varPhi_{\rm f_2} $分别表征时域透镜$ {\rm TL}_1 $$ {\rm TL}_2 $的“焦时间”, 时域透镜$ {\rm TL}_1 $与时域物体$ {\rm TO} $, 时域透镜$ {\rm TL}_1 $与光纤分束器$ {\rm FBS}_2 $之间色散介质的群延迟色散参量也均等于$ \varPhi_{\rm f_1} $, 时域透镜$ {\rm TL}_2 $与光源$ {\rm TS} $, 时域透镜$ {\rm TL}_2 $与光纤分束器$ {\rm FBS}_2 $之间色散介质的群延迟色散参量分别为$ \varPhi_{1} $$ \varPhi_{2} $. 为了便于讨论, 假设时序热光源是时域非相干光的, 其关联性质满足

      $ \left\langle {E_{\rm s}^*\left( {{\tau _0}} \right){E_{\rm s}}\left( {\tau _0^\prime } \right)} \right\rangle = {I_{\rm s}}\delta \left( {{\tau _0} - \tau _0^\prime } \right), $

      其中$ E $表示与在色散介质中传播的时间$ \tau $和距离$ l $有关的光场分布, $ \beta $是光在色散介质中的群速度色散系数. 根据以上方程, 描述脉冲色散的微分方程与描述近轴光波空间衍射的波方程形式相同. 因此时域脉冲在色散介质中的色散可看作空域中光波近轴衍射的时间对应, 即空间-时间二象性. (4)式的解与近轴衍射波方程解的数学形式相似, 可写为

      $\begin{split}\left\langle E_{1}^{*}(\tau) E_{2}(\tau)\right\rangle= &\iint h_{1}^{*}\left(\tau, \tau_{0}\right) h_{2}\left(\tau, \tau_{0}^{\prime}\right)\\&\times\left\langle E_{\rm s}^{*}\left(\tau_{0}\right) E_{\rm s}\left(\tau_{0}^{\prime}\right)\right\rangle {\rm d} \tau_{0} {\rm d} \tau_{0}^{\prime}\end{split}.$

      将(9)式—(11)式代入(12)式, 探测面处的一阶场关联可进一步写为

      $\begin{split} & \left\langle {E_1^*(\tau ){E_2}(\tau )} \right\rangle \\ =\;& {I_{\rm s}}\int {h_1^*} \left( {\tau ,{\tau _0}} \right){h_2}\left( {\tau ,{\tau _0}} \right){\rm d}{\tau _0} \\ =\;& \frac{{{I_{\rm s}}}}{{2{\text{π}} }}\sqrt {\frac{\rm i}{{2{\text{π}} {\varPhi _0}{\varPhi _{\rm R}}{\varPhi _{{\rm f_1}}}}}} \exp \left[ {\frac{{{\rm i}{\tau ^2}}}{{2{\varPhi _{\rm R}}}}\left( {1 - \frac{{{\varPhi _1}}}{{{\varPhi _{{\rm f_2}}}}}} \right)} \right] \\ & \times \iint\exp\left( {\frac{{ - {\rm i}{\tau ^{\prime 2}}}}{{2{\varPhi _0}}} + \frac{{{\rm i}\tau {\tau ^\prime }}}{{{\varPhi _{{\rm f_1}}}}}} \right) \exp \bigg[ - \frac{{{\rm i}\tau _0^2}}{{2{\varPhi _0}}} \\ & + \frac{{{\rm i}\tau _0^2}}{{2{\varPhi _{\rm R}}}}\left( {1 - \frac{{{\varPhi _2}}}{{{\varPhi _{{{\rm f}_2}}}}}} \right) + \frac{{{\rm i}{\tau _0}{\tau ^\prime }}}{{{\varPhi _0}}} - \frac{{{\rm i}\tau {\tau _0}}}{{{\varPhi _{\rm R}}}} \bigg]\\ & \times {P^*}\left( {{\tau ^\prime }} \right){\rm d}{\tau _0}{\rm d}{\tau ^\prime }, \end{split} $

      其中$ \varPhi_{\rm R} = \varPhi_1+\varPhi_2-\varPhi_1\varPhi_2/\varPhi_{\rm f_2} $. 当成像条件$ \dfrac{1}{\varPhi_{\rm f_2}} = \dfrac{1}{\varPhi_2}+\dfrac{1}{\varPhi_1-\varPhi_0} $满足时, 上式化简为

      $\begin{split} & \left\langle E_{1}^{*}\left( \tau \right){{E}_{2}}\left( \tau \right) \right\rangle \\ =\; & {{I}_{\rm s}}\sqrt{\frac{{{\varPhi }_{1}}-{{\varPhi }_{0}}}{{\rm i}2{\text{π}} {{\varPhi }_{2}}{{\varPhi }_{{{\rm f}_{1}}}}}}{{P}^{*}}\left( -\tau \frac{{{\varPhi }_{1}}-{{\varPhi }_{0}}}{{{\varPhi }_{2}}} \right) \\ & \times \exp \left[ {\rm i}{{\tau }^{2}}\frac{{{\varPhi }_{1}}-{{\varPhi }_{0}}}{{{\varPhi }_{2}}}\left( \frac{1}{2{{\varPhi }_{{{\rm f}_{2}}}}}-\frac{1}{{{\varPhi }_{{{\rm f}_{1}}}}} \right) \right].\end{split}$

      (14)式表明时域物体的像会受到二次相位因子的影响. 分析可知, 当各部分色散介质的群延迟色散参量之间关系满足$ \varPhi_1-\varPhi_0 = \varPhi_2 = \varPhi_{\rm f_1} = 2\varPhi_{\rm f_2} $时, 二次相位因子项可以被消掉, 这样(14)式最后化简为

      $ \left\langle {E_1^*\left( \tau \right){E_2}\left( \tau \right)} \right\rangle = {I_{\rm s}}\sqrt {\frac{1}{{{\rm i}2{\text{π}} {\varPhi _2}}}} {P^*}\left( { - \tau } \right). $

      由此可知, 通过一阶场关联可再现时域物体信号的像, 即可以实现时域成像, 不需要额外的色散补偿或消除.

    • 前面讨论了利用时序非相干热光源实现一阶关联时域成像, 考虑到实际光源发出的脉冲宽度和相干时间是有限的, 因此利用更为普遍的热光源模型即高斯-谢尔模型, 分析成像可见度和分辨率. 高斯-谢尔型时序光源的一阶场关联满足

      $ \left\langle {E_{\rm s}^*\left( {{\tau _0}} \right){E_{\rm s}}\left( {\tau _0^\prime } \right)} \right\rangle = {I_{\rm s}}\exp \left[ { - \frac{{\tau _0^2 + {\tau ^\prime }_0^2}}{{T_0^2}} - \frac{{{{\left( {{\tau _0} - {\tau ^\prime }_0} \right)}^2}}}{{2T_{\rm c}^2}}} \right], $

      其中$ T_0 $$ T_{\rm c} $分别表示光源发出脉冲的宽度和相干时间. 探测面处的强度是物光路和参考光路的光场叠加的模方, 即

      $\begin{split}\langle I(\tau)\rangle= \;&{\left\langle\left|E_{1}(\tau)\right|^{2}\right\rangle+\left\langle\left|E_{2}(\tau)\right|^{2}\right\rangle} \\ &+{\left\langle E_{1}^{*}(\tau) E_{2}(\tau)\right\rangle+\left\langle E_{1}(\tau) E_{2}^{*}(\tau)\right\rangle}\end{split}.$

      将(5)式、(9)式、(10)式、(12)式和(16)式代入(17)式, 假设时域物体距离光源很近, 即$ \varPhi_0 = 0 $, 并考虑成像条件和消去二次相位因子条件可得探测面处一阶关联时域成像的强度表达式为

      $ \begin{split} \langle I(\tau )\rangle = &\frac{{{I_{\rm s}}}}{{2{\text{π}} {\varPhi _2}}}\iint{P^*}\left( {{\tau _0}} \right)P\left( {\tau _0^\prime } \right)\exp \left( { - \frac{{{\rm i}\tau _0^\prime \tau }}{{{\varPhi _2}}} + \frac{{{\rm i}{\tau _0}\tau }}{{{\varPhi _2}}}} \right)\\ &\times\exp \left[ { - \frac{{\tau _0^2 + \tau _0'^2}}{{T_0^2}} - \frac{{{{\left( {{\tau _0} - \tau _0^\prime } \right)}^2}}}{{2T_{\rm c}^2}}} \right]{\rm d}{\tau _0}{\rm d}\tau _0^\prime + {I_{\rm s}}\\ &\times\exp \!\left(\!{ - \frac{{2{\tau ^2}}}{{T_0^2}}} \right)\!\!+\! {I_{\rm s}}\sqrt {\frac{2}{{{\text{π}}{\rm i}{\varPhi _2}}}}\!\int \!\!{\exp}\!\left( {\frac{{{\rm i}\tau {\tau _0}}}{{{\varPhi _2}}} \!+\! \frac{{\rm i{\tau ^2}}}{{{\varPhi _2}}}} \right)\\ &\times{P^*}\left( {{\tau _0}} \right)\exp \left[ { - \frac{{\tau _0^2 + {\tau ^2}}}{{T_0^2}} - \frac{{{{\left( {{\tau _0} + \tau } \right)}^2}}}{{2T_{\rm c}^2}}} \right]{\rm d}{\tau _0}. \end{split} $

      假设时域物体为具有一定间隔的两个矩形波信号, 矩形波间隔和宽度分别为20和8 ps, 选取色散介质的群延迟色散参量为$\varPhi_0 = 0, \ \varPhi_1 = \varPhi_2 = $$\varPhi_{{\rm f}_1} = 2\varPhi_{{\rm f}_2} = 50\; {\rm ps}^2 $. 首先分析保持光源脉冲相干时间不变, 改变光源脉冲宽度的成像结果. 根据(18)式, 模拟得到了光源脉冲相干时间为$ T_{\rm c} $ = 0.5 ps, 脉冲宽度$ T_0 $分别为5, 10, 50和100 ps条件下的成像结果, 如图2所示. 根据时域物体成像质量的可见度$ V=\dfrac{{{\left\langle I \right\rangle }_{\rm {max}}}-{{\left\langle I \right\rangle }_{\rm {min}}}}{{{\left\langle I \right\rangle }_{\rm {max}}}+{{\left\langle I \right\rangle }_{\rm {min}}}} $, 图3给出了光源相干时间$ T_{\rm c}= $0.5 ps 条件下, 成像可见度随光源脉冲宽度的变化曲线. 由图2可以看出, 随着脉冲宽度的增加, 时域物体的图像可见度逐渐降低但分辨率逐渐变好, 即时域物体的矩形波逐渐显现. 由图3可以看出, 在保持光源脉冲相干时间不变时, 随光源脉冲宽度的增加, 成像可见度逐渐减小, 这与空间域二阶鬼成像的结果类似[42]. 这是因为当光源脉冲宽度大于时域物体信号持续时间时, 场关联才能完整展现物体的信息.

      图  2  光源相干时间为0.5 ps, 光源脉冲宽度宽度分别为5, 10, 50和100 ps条件下的关联像

      Figure 2.  The correlated images with different pulse widths of the source 5, 10, 50 and 100 ps, for the source’s coherence time 0.5 ps

      图  3  光源相干时间Tc为0.5 ps时, 成像可见度随光源脉冲宽度的变化曲线

      Figure 3.  The imaging visibility as a function of the source’s pulse width for the source’s coherence time 0.5 ps

      此外, Qu等[43]的研究结果表明, 在二阶时域鬼成像中时域物体图像的保真度(表示像与原物体信号的相似度)随脉冲宽度的增加而增加. 图2图3表明随着脉冲宽度的增加, 一阶时域成像也有类似的结果, 即图像与原物体信号的相似度逐渐增加. 一阶时域成像的优点在于成像质量除了与光源的脉冲宽度和相干时间两个参量有关以外, 还与物光和参考光的强度比有关, 相同条件下成像可见度比二阶时域鬼成像的高[39]. 但是一阶时域成像对物光信号不执行桶探测, 需要响应较快的探测器, 而二阶时域鬼成像在这方面具有优势, 因为二阶时域鬼成像对物光执行桶探测只需慢响应探测器即可.

      最后分析当保持光源脉冲宽度不变, 改变光源相干时间的成像结果. 根据(18)式, 模拟得到了光源脉冲宽度为$ T_0 $ = 100 ps, 光源相干时间$ T_{\rm c} $分别为0.5, 1.5, 5和10 ps条件下的成像结果, 如图4所示. 由图4可以看出, 随着相干时间的增加, 时域物体的像逐渐失真, 即成像的分辨率(物体像的半高宽)逐渐变差. 由图5可以看出, 随着光源相干时间的增加, 成像可见度逐渐增加最终基本保持不变. 光源的相干时间类似于空域鬼成像中光源的横向关联长度(散斑尺寸), 光源的相干时间越短, 成像系统对于时域物体的时序结构细节分辨能力越好, 但可见度越差. 这就说明基于一阶关联的时域成像中, 可见度和分辨率二者有相互矛盾的关系. 因此需要适当选取光源的脉冲宽度和相干时间以获得兼顾可见度和分辨率的理想成像结果. 由以上分析可知, 当光源的相干时间为0.5 ps、脉冲宽度为$ T_0 $ = 100 ps 时, 间隔为20 ps、宽度为8 ps 的时域矩形波型物体的成像质量(兼顾可见度和分辨率)较好.

      图  4  光源脉冲宽度为T0 = 100 ps, 光源相干时间Tc分别为0.5, 1.5, 5和10 ps条件下的关联像

      Figure 4.  The correlated images with different source’s coherence time 0.5, 1.5, 5 and 10 ps, for the source’s pulse width 100 ps

      图  5  光源脉冲宽度为100 ps时, 成像可见度随光源相干时间的变化曲线

      Figure 5.  The imaging visibility as a function of the source coherence time for the source’s pulse width 100 ps.

    • 本文根据空间光束近轴衍射和光脉冲在色散介质中色散现象之间的空间-时间二象性, 在时域脉冲响应函数基础上, 研究了基于一阶场关联的时域成像, 并以高斯-谢尔型热光源为例, 模拟再现了时域物体信号的成像结果, 分析了光源脉冲宽度和相干时间对成像可见度和分辨率的影响. 结果表明, 利用一阶场关联在不需要色散补偿或消除的条件下可以实现时域物体信号的成像, 并且随着光源脉冲宽度的增加, 像的可见度逐渐增加, 而随着光源脉冲相干时间的增加, 成像分辨率逐渐变差, 但是可见度逐渐增加. 该结果对于场关联在时域信号测量中的应用具有促进作用.

参考文献 (43)

目录

    /

    返回文章
    返回