1.本发明涉及光学技术领域，特别涉及一种基于模式分解快速产生动态扭曲部分相干光的系统及方法。


背景技术：

2.众所周知，相位是光场的一个重要调控参量。自1993年simon等人提出扭曲相位后，人们通过模式分解、维格纳分布和张量表示等方法对扭曲部分相干光进行了大量的研究，证实了扭曲部分相干光在强度、相干度、偏振度和轨道角动量等性质上具有全新的统计特性，在光束整形、粒子捕获、提高成像分辨率和自由空间光通信等方面有广阔的应用前景。能够实际产生扭曲部分相干光是满足应用要求的基础前提，因此扭曲部分相干光的实际产生方法得到了一定程度的探索。
3.扭曲部分相干光不同于带有一般相位的光束，它的相位具有不可分离性，因此不能通过直接产生电场的方法来获得扭曲部分相干光，而需要通过获得其他易产生的光束进行转化或者叠加得到。光束转化为扭曲部分相干光的原理即先产生一个光源，然后经过一个光学系统将之转化为扭曲部分相干光。比较典型的一个实验是王海云等人根据williamson定理，设计出了一个3柱面镜转化系统，将借助空间光调制器提前产生的各向异性高斯谢尔模光束转化为了扭曲高斯谢尔模光束，若需要产生不同参数的光，则需要调整对应的光源和转化系统。模式叠加方法的原理则是产生对应基模并进行含权叠加得到扭曲部分相干光，根据基模的形式可以分为本征模、伪模、随机模等叠加方法。比较典型的一个实验是tian的伪模叠加方法实验，他们将扭曲高斯谢尔模光束分解为基模为伪模的含权叠加形式。实验中，空间光调制器加载对应全息图，通过叠加得到扭曲高斯谢尔模光束，随机相位起消除不需要的模式关联项的作用。相比于光束转化的方法，伪模叠加的方法优点在于不需要改变光学系统就能产生动态光束。
4.扭曲部分相干光已有的产生方法具有的耗时且不可调性等缺点，给光束的产生带来了较大的限制，很难满足扭曲部分相干光的实际使用要求，因此拥有一种能够快速产生动态可调的扭曲部分相干光的方法显得尤为重要。
5.目前，一共有三种产生扭曲部分相干光束的实验得到了报道。前两种方法基于williamson定理，提前设计好合适的方差矩阵和传输矩阵，然后实验产生对应的初始光场和传输系统，得到所需光场；第三种方法基于模式分解，合成含权基模，得到光场。
6.第一种方法，是friberg等人先通过一个声光偏转器以及滤波系统产生了各向异性高斯谢尔模型光束，再经过一个设计好的特殊位置摆放的由6个柱面镜组成的光学转换系统，最后转换为扭曲高斯谢尔模型光束；而后王海云等人在此理论和实验基础上进行了方法的改进，将扩束准直光照射在加载了全息图的空间光调制器上，直接产生各项异性高斯谢尔模光束且产生效果更好，并且设计了一个3柱面镜系统等效替代了原本的6柱面镜系统，最后实验产生了一个严格精确的扭曲高斯谢尔模型光束；第三种方法是tian等人将扭曲高斯谢尔模型光束分解成了含权伪模的采样叠加形式，并通过给各模式引入随机相位的
方法降低因模式之间的相关性所产生的额外的统计关联项，最终叠加得到了近似的扭曲高斯谢尔模型光束。
7.综上所述，基于williamson定理的两种方法中，一组初始光源和光学转化系统只能产生对应参数的扭曲部分相干光，不满足产生光束动态可调的实际使用要求。基于模式分解的方法中则需要大量引入随机相位的全息图，无法较快产生所需的图片。


技术实现要素：

8.本发明要解决的技术问题是提供一种简单可行、需要模式数少、产生光束耗时短的基于模式分解快速产生动态扭曲部分相干光的系统及方法。
9.为了解决上述问题，本发明提供了一种基于模式分解快速产生动态扭曲部分相干光的系统，其包括：
10.激光器，用于产生激光；
11.第一计算机，用于将扭曲部分相干光束分解为含权的正交完备的本征基模的叠加形式；并根据扭曲部分相干光分解的本征基模叠加形式、要产生的扭曲部分相干光的具体参数大小以及精确度的要求来确定所需要的本征基模的模式及其权重大小；并加载包含各个本征基模的动态全息图至空间光调制器；
12.空间光调制器，用于接收激光并对激光进行调制，得到动态扭曲部分相干光。
13.作为本发明的进一步改进，还包括振幅过滤系统，用于对从所述空间光调制器出射的动态扭曲部分相干光进行过滤，只保留第一衍射级。
14.作为本发明的进一步改进，所述振幅过滤系统为4f系统，包括依次设置的第一透镜、光阑和第二透镜，从所述空间光调制器出射的动态扭曲部分相干光依次经过所述第一透镜、光阑和第二透镜。
15.作为本发明的进一步改进，在所述动态全息图中，各个本征基模的单个全息图按序排列，且显示时间与对应的权重成比例，所需要的本征基模及对应的权重比根据要产生的扭曲部分相干光的具体参数来确定。
16.作为本发明的进一步改进，还包括电荷耦合元件ccd和第二计算机，所述电荷耦合元件ccd和第二计算机连接，所述电荷耦合元件ccd用于接收所述动态扭曲部分相干光，所述第二计算机用于计算并验证所述动态扭曲部分相干光是否与理论一致。
17.作为本发明的进一步改进，还包括扩束镜，设置于所述激光器和空间光调制器之间，用于将所述激光器产生的激光变成扩束准直光。
18.作为本发明的进一步改进，还包括反射平面镜，所述反射平面镜设置于所述激光器和扩束镜之间，用于加长光路。
19.本发明还提供了一种基于模式分解快速产生动态扭曲部分相干光的方法，应用于如上述任一所述的基于模式分解快速产生动态扭曲部分相干光的系统，其包括以下步骤：
20.s1、所述第一计算机将扭曲部分相干光束分解为含权的正交完备的本征基模的叠加形式；
21.s2、所述第一计算机根据扭曲部分相干光分解的本征基模叠加形式、要产生的扭曲部分相干光的具体参数大小以及精确度的要求来确定所需要的本征基模的模式及其权重大小；
22.s2、所述空间光调制器接收所述激光器产生的激光，所述第一计算机加载包含各个本征基模的动态全息图至空间光调制器，利用所述空间光调制器对激光进行调制，得到动态扭曲部分相干光。
23.作为本发明的进一步改进，步骤s1包括：
24.s11、根据模式分解理论以及交叉谱密度函数构建方法，任意满足非负正定性的光束的交叉谱密度可以表示为非负权重函数λ(v)和模式函数φ(ρ，v)对参量v的积分形式：
25.w(ρ，ρ
′
)＝∫λ(v)φ
*
(ρ，v)φ(ρ
′
，v)d2v
ꢀꢀꢀ
(1)
26.s12、将λ(v)视为一组由带有正系数的狄拉克函数组成，即：
[0027][0028]
将式(2)的连续积分形式表示为模式的离散叠加形式：
[0029][0030]
其中，n代表模式的阶数，φ(ρ，vn)是模式函数，λ(vn)是对应的权重函数；
[0031]
s13、将扭曲部分相干光的交叉谱密度表示为：
[0032]
w(ρ，ρ
′
)＝a exp[f(σ，δ，ρ，ρ
′
)-ikμ(ρ
×
ρ
′
)
⊥
]
ꢀꢀꢀ
(4)
[0033]
其中，σ、δ和μ分别表示束腰宽度、相干长度和扭曲因子，f(σ，δ，ρ，ρ
′
)表示对应扭曲部分相干光的具体关联结构；
[0034]
s14、当f(σ，δ，ρ，ρ
′
)确定时，结合式(3)和式(4)，确定满足正交完备性的模式函数φ(ρ，vn)及其权重函数λ(vn)的具体形式，将对应的扭曲部分相干光进行本征模式分解，表示为：
[0035][0036]
作为本发明的进一步改进，步骤s2包括：
[0037]
s21、给定一个模式范围；
[0038]
s22、在给定的范围中根据式(5)中确定的权重函数计算各本征基模的权重值，同时进行归一化得到归一化权重值分布，按从大到小的顺序排列作为选取模式的优先级；
[0039]
s23、根据精度要求确定一个归一化权重值的标准，与上一步确定的优先级结合，选取归一化权重值在该标准上的本征基模，由此确定具体需要的本征基模及其权重。
[0040]
本发明的有益效果：
[0041]
本发明基于模式分解快速产生动态扭曲部分相干光的系统及方法基于比例确定的权重大小选取模式，需要的模式数少，产生光束耗时更短；只需要更改加载在空间光调制器中的全息图，就可以实现扭曲部分相干光束的动态调控；模式之间相互正交，不会产生额外的交叉项。
[0042]
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂，以下特举较佳实施例，并配合附图，详细说明如下。
附图说明
[0043]
图1是本发明优选实施例中基于模式分解快速产生动态扭曲部分相干光的系统的示意图一；
[0044]
图2是本发明优选实施例中基于模式分解快速产生动态扭曲部分相干光的系统的示意图二；
[0045]
图3是本发明优选实施例中基于模式分解快速产生动态扭曲部分相干光的系统在理论和实验中产生的动态扭曲部分相干光的强度、相干度分布及拟合情况图；
[0046]
图4是本发明优选实施例中基于模式分解快速产生动态扭曲部分相干光的系统在理论和实验中产生的动态扭曲部分相干光的不同距离的传输旋转情况图。
[0047]
标记说明：1、激光器；2、反射平面镜；3、扩束镜；4、空间光调制器；5、第一透镜；6、光阑；7、第二透镜；8、电荷耦合器件ccd；9、第一计算机；10、第二计算机；11、柱透镜；12、第三透镜。
具体实施方式
[0048]
下面结合附图和具体实施例对本发明作进一步说明，以使本领域的技术人员可以更好地理解本发明并能予以实施，但所举实施例不作为对本发明的限定。
[0049]
如图1所示，本发明优选实施例中的基于模式分解快速产生动态扭曲部分相干光的系统的示意图一，该系统包括：
[0050]
激光器1，用于产生激光；可选的，激光器1为氦氖激光器。
[0051]
第一计算机9，用于将扭曲部分相干光束分解为含权的正交完备的本征基模的叠加形式；并根据扭曲部分相干光分解的本征基模叠加形式、要产生的扭曲部分相干光的具体参数大小以及精确度的要求来确定所需要的本征基模的模式及其权重大小；并加载包含各个本征基模的动态全息图至空间光调制器；
[0052]
空间光调制器4(slm)，用于接收激光并对激光进行调制，得到动态扭曲部分相干光。
[0053]
进一步地，还包括振幅过滤系统，用于对从所述空间光调制器出射的动态扭曲部分相干光进行过滤，只保留第一衍射级。
[0054]
具体地，所述振幅过滤系统为4f系统，包括依次设置的第一透镜5、光阑6和第二透镜7，从所述空间光调制器4出射的动态扭曲部分相干光依次经过所述第一透镜5、光阑6和第二透镜7。光束经过4f系统后得到翻转平移的本征模光束，经过叠加合成得到所需的扭曲部分相干光。
[0055]
在一些实施例中，还包括扩束镜3，设置于激光器1和空间光调制器4之间，用于将激光器1产生的激光变成扩束准直光。
[0056]
可选地，还包括反射平面镜2，反射平面镜2设置于激光器1和扩束镜3之间，用于加长光路。
[0057]
在所述动态全息图中，各个本征基模的单个全息图按序排列，且显示时间与对应的权重成比例，所需要的本征基模及对应的权重比根据要产生的扭曲部分相干光的具体参数来确定。如束腰、相干长度、扭曲相位等。
[0058]
在一些实施例中，还包括电荷耦合元件ccd8和第二计算机10，所述电荷耦合元件
ccd8和第二计算机10连接，所述电荷耦合元件ccd8用于接收所述动态扭曲部分相干光，所述第二计算机10用于计算并验证所述动态扭曲部分相干光是否与理论一致。
[0059]
本发明优选实施例还公开了一种基于模式分解快速产生动态扭曲部分相干光的方法，应用于上述基于模式分解快速产生动态扭曲部分相干光的系统，其包括以下步骤：
[0060]
s1、所述第一计算机将扭曲部分相干光束分解为含权的正交完备的本征基模的叠加形式；
[0061]
s2、所述第一计算机根据扭曲部分相干光分解的本征基模叠加形式、要产生的扭曲部分相干光的具体参数大小以及精确度的要求来确定所需要的本征基模的模式及其权重大小；
[0062]
s2、所述空间光调制器接收所述激光器产生的激光，所述第一计算机加载包含各个本征基模的动态全息图至空间光调制器，利用所述空间光调制器对激光进行调制，得到动态扭曲部分相干光。
[0063]
进一步地，步骤s1包括：
[0064]
s11、根据模式分解理论以及gori等人提出的交叉谱密度函数构建方法，任意满足非负正定性的光束的交叉谱密度(csd)可以表示为非负权重函数λ(ν)和模式函数φ(ρ,ν)对参量ν的积分形式：
[0065]
w(ρ,ρ
′
)＝∫λ(ν)φ
*
(ρ,ν)φ(ρ
′
,ν)d2v
ꢀꢀꢀꢀꢀꢀ
(1)
[0066]
s12、将λ(ν)视为一组由带有正系数的狄拉克函数组成，即：
[0067][0068]
将式(2)的连续积分形式表示为模式的离散叠加形式：
[0069][0070]
其中，n代表模式的阶数，φ(ρ,νn)是模式函数，λ(νn)是对应的权重函数；
[0071]
s13、将扭曲部分相干光的交叉谱密度表示为：
[0072]
w(ρ,ρ
′
)＝aexp[f(σ,δ,ρ,ρ
′
)-ikμ(ρ
×
ρ
′
)
⊥
]
ꢀꢀꢀꢀꢀ
(4)
[0073]
其中，σ、δ和μ分别表示束腰宽度、相干长度和扭曲因子，f(σ,δ,ρ,ρ
′
)表示对应扭曲部分相干光的具体关联结构；
[0074]
s14、当f(σ,δ,ρ,ρ
′
)确定时，结合式(3)和式(4)，确定满足正交完备性的模式函数φ(ρ,νn)及其权重函数λ(νn)的具体形式，将对应的扭曲部分相干光进行本征模式分解，表示为：
[0075][0076]
进一步地，步骤s2包括：
[0077]
s21、给定一个模式范围；其中，虽然无法取遍所有的模式，但是可以给定一个的模式范围，若继续扩大该范围，强度几乎不会发生变化，则可以认为已经取到所需要的模式。
[0078]
s22、在给定的范围中根据式(5)中确定的权重函数计算各本征基模的权重值，同时进行归一化得到归一化权重值分布，按从大到小的顺序排列作为选取模式的优先级；使得只需要更少的模式数就能产生同等精确度的扭曲部分相干光。
[0079]
s23、根据精度要求确定一个归一化权重值的标准，与上一步确定的优先级结合，选取归一化权重值在该标准上的本征基模，由此确定具体需要的本征基模及其权重。
[0080]
由于扭曲部分相干光束的交叉谱密度可以表示式(5)，激光经过空间光调制器4后同时刻产生了加载的全息图对应的本征基模φ
n*
(ρ)φn(ρ
′
)，由于其播放时间与权重成比例，故在阶数n的全息图的持续时间段内，产生的光束的csd为其他模式同理，当整个动态全息图播放一轮，得到的光束的csd可以表示为：
[0081][0082]
即为所需要的扭曲部分相干光束的近似。
[0083]
如果需要产生不同参数的扭曲部分相干光束，如束腰宽度、相干长度和扭曲因子等，只要通过计算好所需的模式及其权重，就能在仅仅改变加载的动态全息图而不改变光学系统的情况下实现扭曲部分相干光束的动态产生。同时该方法只需要扩大或缩小模式范围(增加或减少叠加模式)就能对产生的扭曲部分相干光实现精确度的调整。且由于该方法是在正交完备基中按照权重大小来选取模式的，需要的模式少，实验产生更快速，仅仅受限于slm的刷新频率；且没有伪模及随机模模式之间的交叉关联项，实验产生的光束稳定。
[0084]
如图2所示，本发明优选实施例中的基于模式分解快速产生动态扭曲部分相干光的系统的示意图二，其基于本征模分解形式产生扭曲高斯谢尔模(tgsm)光束。
[0085]
根据上式(6)以及扭曲tgsm的csd解析式：
[0086][0087]
可以将tgsm表示为本征基模为拉盖尔高斯光束的含权叠加，即
[0088][0089]
为权重函数，是归一化拉盖尔形式的模式函数，表示为：
[0090][0091][0092]
其中ξ、t、w分别是由σ、δ和μ共同确定的参数：
[0093][0094][0095]
[0096]
且a＝1/(4σ2),b＝1/(2δ2),u＝kμ。
[0097]
根据式(8)、(9)、(10)可知，可以用不同含权的拉盖尔模式组叠加产生不同参数的tgsm光束，我们选取要产生的tgsm光束参数为：σ＝0.5mm,δ＝0.28mm,μ＝0.001mm-1
。根据式(9)中的和式(10)，可以得到各模式的归一化本征值谱选取了一个标准值为0.01。归一化本征值大于该标准的模式选用，小于该标准的模式忽略，最终确定的模式为：阶数n＝0，拓扑核数为m＝-2～25；n＝1，m＝0～11，一共40个模式。
[0098]
根据确定的模式及其权重使用matlab运行出动态全息图，动态全息图中各模式全息图按照n＝0，m＝-2～25；n＝1，m＝0～11的顺序排列，且播放时间与权重成比例，将动态全息图加载到slm上，使用的slm型号是bq-slm1024，像素数为1024
×
768，像素大小为18μm
×
18μm，在图2中将波长为632nm的扩束准直光入射到slm上，再经过两个焦距为15cm的透镜组成的振幅过滤系统后，在振幅过滤系统后焦面产生不同模式的拉盖尔高斯谢尔模光束。
[0099]
为了验证实验产生的光束是σ＝0.5mm,δ＝0.28mm,μ＝0.001mm-1
的tgsm光束，使用电荷耦合器件ccd8在4f系统后焦面接收了实验产生的各模式的拉盖尔高斯光束的强度照片，用于验证束腰宽度和相干长度和理论设计是否一致；该实施例在上述实施例的基础上加上了一个焦距为15cm的柱面镜11和焦距为20cm的第三透镜12，根据传输旋转情况验证扭曲因子。使用的ccd型号是gs3-u3-28s5m-c，像素为1920
×
1440，像素大小为4.54μm。
[0100]
其中束腰宽度可以由光束的强度来验证，根据式(7)，理论的强度分布均由束腰宽度唯一确定：
[0101][0102]
根据式(8)，实验的强度分布为各个模式强度的含权叠加，亦即ccd接收照片强度的叠加：
[0103][0104]
相干长度可以由权重谱、多张照片的强度分布以及某一点强度(选取了中心点)来表示，其过程如下：
[0105]
根据式(7)可得，中心点与其他点的理论相干度分布可以由相干长度唯一确定：
[0106][0107]
而根据式(8)，实验合成的扭曲部分相干光束的相干度可以表示为：
[0108][0109]
而由于本征模式之间相互正交，故可以简化为：
[0110][0111]
即相干长度可以由权重谱、强度分布以及中心点强度来表示。
[0112]
扭曲相位是造成光束传输旋转的唯一因素，因此在波长一定的情况下，扭曲因子和传输旋转情况一一对应，根据tgsm光束的张量传输定律，我们用matlab模拟了tgsm光束在图二中的理论传输旋转情况，并和实验旋转结果进行了拟合。
[0113]
图3是本发明优选实施例中基于模式分解快速产生动态扭曲部分相干光的系统在理论和实验中产生的动态扭曲部分相干光的强度、相干度分布及拟合情况图；其中，第一行及第二行分别是强度拟合和相干度的模方拟合，验证了实验产生光束的束腰宽度及相干长度和理论设计一致。
[0114]
图4是本发明优选实施例中基于模式分解快速产生动态扭曲部分相干光的系统在理论和实验中产生的动态扭曲部分相干光的不同距离的传输旋转情况图。其中，第一行和第二行分别是理论和实验的结果，验证了实验产生光束的扭曲因子和理论设计一致。
[0115]
以上实施例仅是为充分说明本发明而所举的较佳的实施例，本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换，均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。