CN108594444B - 基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法 - Google Patents

基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法 Download PDF

Info

Publication number
CN108594444B
CN108594444B CN201810314144.2A CN201810314144A CN108594444B CN 108594444 B CN108594444 B CN 108594444B CN 201810314144 A CN201810314144 A CN 201810314144A CN 108594444 B CN108594444 B CN 108594444B
Authority
CN
China
Prior art keywords
film
mathieu
cone mirror
modulation
amplitude modulation
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
CN201810314144.2A
Other languages
English (en)
Other versions
CN108594444A (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.)
Zhejiang Normal University CJNU
Original Assignee
Zhejiang Normal University CJNU
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 Zhejiang Normal University CJNU filed Critical Zhejiang Normal University CJNU
Priority to CN201810314144.2A priority Critical patent/CN108594444B/zh
Publication of CN108594444A publication Critical patent/CN108594444A/zh
Application granted granted Critical
Publication of CN108594444B publication Critical patent/CN108594444B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/09Beam shaping, e.g. changing the cross-sectional area, not otherwise provided for
    • G02B27/0938Using specific optical elements

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
  • Diffracting Gratings Or Hologram Optical Elements (AREA)

Abstract

基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法,是一种精确、高效地产生一簇无衍射Mathieu光束的方法。该方法基于角向Mathieu函数的定义,首先构造具有角向Mathieu函数分布的振幅分布灰度图,再将灰度图颜色黑白反转,用胶片输出仪打印在高分辨率胶片上,打印好的胶片经曝光处理后,得到具有角向Mathieu函数分布的振幅调制胶片。将制作的胶片紧贴锥镜入瞳面放置,入射的平行光束经胶片调制振幅、锥镜调制相位后,在锥镜后产生Mathieu光束。本发明实验装置结构简单,可以方便、精确地产生Mathieu光束。

Description

基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法
技术领域
本发明涉及结构光束的产生、调控及变换技术,是一种利用胶片和锥镜精确产生Mathieu光束的方法。
背景技术
很多科学实验,比如光与物质相互作用,需要形态稳定不变的光束。但衍射是最普遍的光学现象。光束传输过程中的衍射效应,导致光束展宽及形变。为克服光波的衍射,近年来引入和产生的无衍射结构光束,是构建传输不变光束的一种重要手段。无衍射结构光束是波动方程的特解,这类光束的一个共同特点是,光束分布与传输距离无关。因此无衍射结构光束能在光与物质相互作用的许多研究方面,发挥着特殊的作用。1987年,在自由空间产生的第一种具有实用价值的无衍射光束是Bessel光束。
在实际的科学研究中,不同类型的科学实验需要不同光学形态的无衍射光束。通过激光调制技术,实现光束分布及其光束拓扑特性的多样化,是激光技术的一个重要研究方向。2000年,在椭圆柱坐标下,通过求解波动方程,人们找到了一种新的传输不变解——Mathieu函数解[Opt.Lett.25(20),1493(2000)]。这表明,满足Mathieu函数分布的光束,也是一种无衍射光束。相对于Bessel光束,Mathieu光束的光学形态更为丰富(事实上,Bessel光束只是Mathieu光束的一个特例),但实验产生也复杂的多。
在无衍射Mathieu光束的产生方面,Marcelino Anguiano-Morales等利用锥镜和有一定张角的透明扇形相结合,产生了光学形态接近Mathieu函数的近似Mathieu光束[Opt.Commun.281,401(2008)]。但近似Mathieu光束,并不具有理想的无衍射传输特性。
在实际的使用中,唯有精确产生Mathieu光束,才能将其更好地用于实际的科学研究。利用Durnin装置,可以精确地产生Mathieu光束,见引文[Opt.Commun.222,75(2003)]和[光学学报,37(2):0205002(2017)]。但因为Durnin装置只能利用透射圆环部分的少量光能量产生Mathieu光束,对光能的利用率较低。
发明内容
本发明提供了一种高效、精确产生Mathieu光束的方法。通过具有角向Mathieu函数振幅分布的胶片作为振幅调制元件,合并锥镜的相位调制特性产生Mathieu光束。
本发明利用胶片输出仪(型号:Propalette 8000)制作具有角向Mathieu函数振幅分布的振幅调制元件。首先通过角向Mathieu函数,计算得到具有角向Mathieu函数的振幅分布灰度图,再将灰度图颜色黑白反转,用胶片输出仪打印在高分辨率感光胶片上。打印好的胶片经过照相底片冲洗工艺,得到振幅调制元件——振幅调制胶片。
胶片制作好后,将胶片紧贴锥镜入射面固定,入射平行光经胶片振幅调制和锥镜相位调制后,经菲涅尔衍射后产生Mathieu光束。
本发明产生Mathieu光束的理论过程如下:
根据奇、偶模及光束相对于
Figure BSA0000162117410000021
Figure BSA0000162117410000022
的对称性,Mathieu光束可分为四大类。分别为:第一、二类偶型Mathieu光束,可写为ce2m(η,q)Je2m(ξ,q)和ce2m+1(η,q)Je2m+1(ξ,q),m=0,1,2,...,ce和Je分别为偶型角向、径向马丢函数;第一、二类奇型Mathieu光束,可写为se2m+2(η,q)Jo2m+2(ξ,q)和se2m+1(η,q)Jo2m+1(ξ,q),m=1,2,3,...,其中,se和Jo分别为奇型角向、径向Mathieu函数。式中q是无量纲椭圆参量。显然,Mathieu光束的阶次m和无量纲椭圆参量q共同决定了Mathieu光束的分布。
以第一类偶阶角向Mathieu函数为例,第一类偶阶角向Mathieu函数的傅里叶展开形式为
Figure BSA0000162117410000023
其中,
Figure BSA0000162117410000024
为径向波数,kz是轴向波数,h为椭圆参数。
设加工的振幅调制胶片具有式(1)所示的透过率函数。
在利用锥镜产生Mathieu光束的方案中,径向波数kr取决于锥镜。由
Figure BSA0000162117410000025
可知参量kr和h共同决定了q。给定kr的情况下,h的大小决定了q;给定h的情况下,kr的大小决定了q。
T(r)为锥镜变化函数
Figure BSA0000162117410000026
式中,k是波数,n为锥镜折射率,θ0为锥镜底角(锥镜入射平面和出射锥面间的夹角),R为锥镜入瞳半径。
假定振幅为E0平行光先经过振幅调制元件调制后,再经锥镜调制相位后,锥镜后的光场复振幅分布为
Figure BSA0000162117410000031
锥镜后的光场的传输可写为菲涅尔衍射的形式:
Figure BSA0000162117410000032
式中,
Figure BSA0000162117410000033
为初始光场(场源),r,
Figure BSA0000162117410000034
分别为场源的径向距离和方位角。U(ρ,θ,z)为轴向距离为z处的光场分布,ρ,θ分别为观察面场的径向距离和方位角。
在推导式(4)的过程中,下面的关系式(5)将被使用
Figure BSA0000162117410000035
式中,C2m=2πA0/ce2m(0;q)ce2m(π/2;q)为比例常数,不影响光束分布。
利用式(5),方程(4)可推导为如下的形式
Figure BSA0000162117410000036
上式中,有
Figure BSA0000162117410000037
Figure BSA0000162117410000038
g(r)=Je2m(ξ,qz)ce2m(η,qz)r。当k→∞,稳相法能用来求解形如∫g(r)exp[ikf(r)]的积分表达式。根据
Figure BSA0000162117410000039
可求得稳相点为r0=(n-1)θ0z。当r=r0∈[0,R],方程(6)可写为
Figure BSA00001621174100000310
(7)式即为锥镜后的光场振幅分布。(7)式表明,锥镜后的光束可以被看作无衍射Mathieu光束。
以上理论推导结果证明,利用胶片振幅调制合并锥镜相位调制的方法,可精确地产生无衍射Mathieu光束。
与过去利用Durnin装置产生Mathieu光束的方法相比,本发明利用胶片振幅调制和锥镜相位调制的方法,将胶片紧贴锥镜入射面固定,对光能的利用率更高,可以高效、精确地产生Mathieu光束。另外,通过加工任意类型、任意阶次、任意椭圆参数的角向Mathieu函数分布的胶片,可以方便灵活、精确地产生相应的任意Mathieu光束。
附图说明
图1为本发明通过胶片振幅调制和锥镜相位调制精准产生Mathieu光束的示意图
具体实施方式
以下结合附图与实施例对本发明作进一步的说明。
图1为本发明通过胶片振幅调制和锥镜相位调制精准产生Mathieu光束的实验装置图,实验装置包括:激光束(1)、扩束准直系统(2、3)、胶片(4)、锥镜(5)、科学CCD相机(6)。
实验实施过程如下:
1、通过角向Mathieu函数计算,得到具有角向Mathieu函数的振幅分布图,制作具有该振幅分布图的胶片(4)。胶片(4)的具体制作步骤:通过角向Mathieu函数,计算得到具有角向Mathieu函数的振幅分布灰度图;将灰度图颜色进行黑白反转;用胶片输出仪将颜色反转后的灰度图输出在高分辨率感光胶片上;打印好的胶片经过照相底片冲洗工艺,得到振幅调制元件——振幅调制胶片。
2、激光束(1)经过扩束准直系统(2、3组成)对光束进行扩束准直后,入射至胶片(4)上。
3、胶片(4)对入射的平行光进行振幅调制。
4、振幅调制后的光束经锥镜(5)相位调制后,在锥镜(5)的焦平面后用科学CCD相机(6)记录Mathieu光束。
5、加工任意类型、任意阶次、任意椭圆参数的角向Mathieu函数分布的胶片(4),即可得到相应的不同类型、不同阶次、不同椭圆系的Mathieu光束。

Claims (3)

1.基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法,其特征在于:利用Propalette 8000胶片输出仪制作具有角向Mathieu函数分布的振幅调制胶片;将胶片紧贴锥镜入瞳面,入射激光束经胶片振幅调制和锥镜调制相位,在锥镜后精确产生Mathieu光束。
2.根据权利要求1所述的基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法,其特征在于该方法的具体步骤如下:
(1)、构造具有角向Mathieu函数分布的振幅分布灰度图;
(2)、将灰度图颜色黑白反转,输入胶片输出仪打印在高分辨率感光胶片上,经曝光处理后得到振幅调制胶片;
(3)、用制作的胶片调制入射激光束的振幅分布;
(4)、将胶片紧贴锥镜入瞳面放置,胶片振幅调制后的光束经锥镜相位调制,在锥镜的焦平面后用科学CCD相机记录产生的Mathieu光束。
3.根据权利要求1所述的基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法,其特征在于:该方法可以制作任意类型、任意阶次、任意椭圆参数的角向Mathieu函数分布的振幅调制胶片,精确地产生相应的任意Mathieu光束。
CN201810314144.2A 2018-03-28 2018-03-28 基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法 Active CN108594444B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810314144.2A CN108594444B (zh) 2018-03-28 2018-03-28 基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810314144.2A CN108594444B (zh) 2018-03-28 2018-03-28 基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法

Publications (2)

Publication Number Publication Date
CN108594444A CN108594444A (zh) 2018-09-28
CN108594444B true CN108594444B (zh) 2020-12-04

Family

ID=63621436

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810314144.2A Active CN108594444B (zh) 2018-03-28 2018-03-28 基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法

Country Status (1)

Country Link
CN (1) CN108594444B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112882243A (zh) * 2021-02-02 2021-06-01 浙江师范大学 一种基于稳相法构建椭圆螺旋马蒂厄涡旋光束的方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1997027509A1 (en) * 1996-01-26 1997-07-31 California Institute Of Technology Optical system using a radial harmonic pupil filter for generating collimated beams

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2821883C2 (de) * 1978-05-19 1980-07-17 Ibm Deutschland Gmbh, 7000 Stuttgart Vorrichtung zur Materialbearbeitung
EP0525801B1 (en) * 1991-07-31 1998-12-30 Nec Corporation Method and arrangement of generating a non-diffractive beam at a location which is remote from optical element and application thereof
WO2007014773A1 (en) * 2005-08-02 2007-02-08 Carl Zeiss Laser Optics Gmbh Optical system for attenuating and imaging an optical beam for a subsequent intensity measurement
DE102007032181B4 (de) * 2007-07-11 2012-02-16 Karlsruher Institut für Technologie Optische Anordnung und ihre Verwendung
US9791684B2 (en) * 2010-01-06 2017-10-17 Ecole polytechnique fédérale de Lausanne (EPFL) Optical coherence microscopy system having a filter for suppressing a specular light contribution
CN101975991B (zh) * 2010-09-30 2012-05-23 浙江大学 基于振幅调制与偏振及相位相结合扩展焦深的方法和装置
CN102393254B (zh) * 2011-08-23 2013-08-14 中国计量学院 一种拉格朗日高斯光束古依相位参量测量方法
CN102778820A (zh) * 2012-08-17 2012-11-14 杭州士兰明芯科技有限公司 基于空间光调制器的无掩膜图形曝光系统
CN103645562B (zh) * 2013-12-20 2015-10-21 中南大学 一种同时调制振幅和相位的光束整形方法
CN103760673A (zh) * 2014-01-06 2014-04-30 华侨大学 一种产生近似无衍射零阶Mathieu beam的光学系统
CN103792663B (zh) * 2014-01-17 2015-11-18 北京空间机电研究所 一种生成螺旋式贝塞尔光束的光学系统及生成方法
WO2015117115A1 (en) * 2014-02-03 2015-08-06 President And Fellows Of Harvard College Three-dimensional super-resolution fluorescence imaging using airy beams and other techniques
WO2016073785A1 (en) * 2014-11-05 2016-05-12 The Regents Of The University Of Colorado 3d imaging, ranging, and/or tracking using active illumination and point spread function engineering
US10047001B2 (en) * 2014-12-04 2018-08-14 Corning Incorporated Glass cutting systems and methods using non-diffracting laser beams
CN104777612B (zh) * 2015-04-28 2017-01-25 中南大学 一种对光束的振幅和相位同时整形的迭代方法
CN105182543B (zh) * 2015-07-23 2017-07-14 浙江师范大学 一种产生Pearcey光束簇的方法
CN105353507A (zh) * 2015-10-21 2016-02-24 昌吉市创惠信息科技有限公司 一种医疗光探测器的椭圆型光学衰减片设计方法
CN105607266B (zh) * 2016-01-06 2017-11-14 北京理工大学 三维矢量光束的生成方法与装置
WO2018029678A1 (en) * 2016-08-07 2018-02-15 Ramot At Tel-Aviv University Ltd. Method and system for imaging internal medium
CN107238575B (zh) * 2017-06-08 2019-12-20 深圳大学 一种基于完美涡旋光激发spr的光声显微系统
CN107621701B (zh) * 2017-09-07 2023-08-25 苏州大学 产生双指数贝塞尔高斯光束的方法及系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1997027509A1 (en) * 1996-01-26 1997-07-31 California Institute Of Technology Optical system using a radial harmonic pupil filter for generating collimated beams

Also Published As

Publication number Publication date
CN108594444A (zh) 2018-09-28

Similar Documents

Publication Publication Date Title
Yu et al. Ultrahigh enhancement of light focusing through disordered media controlled by mega-pixel modes
Liu et al. 4K-DMDNet: diffraction model-driven network for 4K computer-generated holography
Chen et al. Phase hologram optimization with bandwidth constraint strategy for speckle-free optical reconstruction
CN107942523A (zh) 一种基于光强传输测量计算的相位恢复系统
CN108594444B (zh) 基于胶片振幅调制和锥镜相位调制产生Mathieu光束的方法
Zhang et al. Ten-mega-pixel snapshot compressive imaging with a hybrid coded aperture
Khonina et al. Caustics of non-paraxial perfect optical vortices generated by toroidal vortex lenses
Christmas et al. Displays based on dynamic phase-only holography
Makowski et al. Miniature holographic projector with cloud computing capability
Gao et al. Compact and lightweight panoramic annular lens for computer vision tasks
CN105607393A (zh) 利用相位型空间光调制器的可变焦投影显示系统
Tudor et al. Multilevel Spiral Axicon for High-Order Bessel–Gauss Beams Generation
CN109613712A (zh) 一种采用轴棱锥相位分布和方位角随机调制产生无衍射散斑的方法及其实现装置
Kim et al. Accelerated synthesis of wide-viewing angle polygon computer-generated holograms using the interocular affine similarity of three-dimensional scenes
Davis et al. Programmable zoom lens system with two spatial light modulators: limits imposed by the spatial resolution
CN112882243A (zh) 一种基于稳相法构建椭圆螺旋马蒂厄涡旋光束的方法
KR100664872B1 (ko) 1차원 홀로그램을 이용한 영상 표시장치
CN112327557A (zh) 一种自加速光生成方法及装置
CN207586561U (zh) 一种基于光强传输测量计算的相位恢复系统
Kotlyar et al. Geometric progression of optical vortices
Junker et al. Scalar wave-optical reconstruction of plenoptic camera images
Petrov et al. Multi-view 3D integral imaging systems using projectors and mobile devices
CN109257524B (zh) 一种基于Wigner分布函数的全聚焦扫描成像方法
Khorin et al. Composite Diffraction-Free Beam Formation Based on Iteratively Calculated Primitives
TWI817692B (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