一种用于实现高斯光束整形为平顶光束的方法，本发明涉及非成像光学领域，尤其涉及一种用于将高斯光束整形为平顶光束的方法。本发明是要解决全局算法计算时间长，局部算法初始相位的设定对结果影响很大的问题，而提供了一种用于实现高斯光束整形为平顶光束的方法。（1）确定高斯光束的直径D1，平顶光束的直径D2，以及两个光学元件的间距L；（2）根据能量守恒计算平顶光束的光强；（3）根据光线追迹建立一一对应关系；（4）计算光线偏角；（5）确定整形元件的初始相位分布；（6）确定相位校正元件的初始相位分布；（7）局部算法优化计算相位分布。本发明应用于非成像光学领域。
1.一种用于实现高斯光束整形为平顶光束的方法，其特征在于用于实现高斯光束整形为平顶光束的方法按以下步骤实现：(1)确定高斯光束的直径D1，平顶光束的直径D2，以及两个光学元件的间距L：若激光高斯光束的半径为R1，则D1＝2R1；若需求的平顶光束的半径为R2，则D2＝2R2；为保证经过光学整形元件后的光线的偏角尽可能的小，取L≥D1；
(5)确定整形元件的初始相位分布：由光程函数与光线偏角间的关系，可以得到对上式积分，得到整形元件的初始相位分布；
(7)局部算法优化计算相位分布：利用光线追迹求得整形元件和相位校正元件的初始相位分布 和 后，将其作为局部算法的初始量，通过GS算法优化计算，最后求得优化后的相位分布，即完成了利用光线追迹法和局部算法求得整形元件和相位校正元件的相位，实现了高斯光束整形为平顶光束。
2.根据权利要求1所述的一种用于实现高斯光束整形为平顶光束的方法，其特征在于步骤(2)中：平顶光束的光强I2可由下式确定式中，I1为高斯光束的中心光强，ω1为高斯光束的束腰半径。
3.根据权利要求1所述的一种用于实现高斯光束整形为平顶光束的方法，其特征在于步骤(3)中根据光线追迹建立一一对应关系具体为：P1面上从距离光束中心r1处发出的光线入射到P2面上距离光束中心r2处，其中
4.根据权利要求1所述的一种用于实现高斯光束整形为平顶光束的方法，其特征在于步骤(6)中确定相位校正元件的初始相位分布具体为：经过整形元件调制后的光束在位相校正元件的前表面P2上的光强分布为平顶的，但是传播方向不再沿原光轴传播，需要通过位相校正元件的调制，使得整形后的光束沿着原光轴传播，由光程函数与光线偏角间的关系可以得到对上式积分，得到相位校正元件的初始相位分布。
5.根据权利要求1所述的一种用于实现高斯光束整形为平顶光束的方法，其特征在于步骤(7)中GS算法优化计算具体为：一、选取DOE的位相初始值 与输入面上的已知振幅分布A0(x,y)组成输入面复振幅分布U1(x,y)；
二、对U1(x,y)作光学正变换，在输出面得到再现像的场分布U′2(u,v)；
三、为使输出面的再现像不断逼近目标图像，保持再现像的位相分布不变，用目标图像的已知振幅分布B0(u,v)替代再现像的振幅分布|U′2(u,v)|；
四、对经过修正后的再现像的场分布U2(u,v)做光学逆变换，在输入面得到场分布U′1(x,y)；
五、返回第一步，以U′1(x,y)的位相部分取代 作下一次迭代。
一种用于实现高斯光束整形为平顶光束的方法
技术领域
[0001] 本发明涉及非成像光学领域，尤其涉及一种用于将高斯光束整形为平顶光束的方法。
背景技术
[0002] 自激光问世以来，其受到了广泛的应用。在有的应用领域对激光束的均匀性有较高要求，例如在激光光刻、激光加工、激光存储和记录、医疗领域等方面，因此需要通过特殊的设计方法将高斯光束整形为均匀分布的平顶光束。
[0003] 光束的整形方法有很多，利用二元光学元件、非球面透镜组、全息滤波器、振幅调制光栅、液晶空间光调制器、微透镜阵列整形、双折射透镜组、相位型光束整形均可将高斯光束整形为平顶光束。
[0004] 比较而言，非球面透镜组的方法结构简单，但是基本未考虑衍射效应。二元光学元件整形的理论基础（光波的衍射理论）相对来说更完整，目前一般采用迭代算法进行求解得到相位分布；迭代算法主要分为全局算法和局部算法，全局算法能有效找到最优值，但是计算时间长，局部算法计算时间短，但是算法中初始相位的设定对结果影响很大。
发明内容
[0005] 本发明是要解决全局算法计算时间长，局部算法初始相位的设定对结果影响很大的问题，而提供了一种用于实现高斯光束整形为平顶光束的方法。
[0006] 用于实现高斯光束整形为平顶光束的方法按以下步骤实现：
[0007] （1）确定高斯光束的直径D1，平顶光束的直径D2，以及两个光学元件的间距L：若激光高斯光束的半径为R1，则D1=2R1；若需求的平顶光束的半径为R2，则D2=2R2；为保证经过光学整形元件后的光线的偏角尽可能的小，一般取L≥D1；
[0008] （2）根据能量守恒计算平顶光束的光强；
[0009] （3）根据光线追迹建立一一对应关系；
[0010] （4）计算光线偏角：由几何关系，可以确定
[0011]
[0012] （5）确定整形元件的初始相位分布：由光程函数与光线偏角间的关系，可以得到[0013]
[0014] 对上式积分，得到整形元件的初始相位分布；
[0015] （6）确定相位校正元件的初始相位分布；
[0016] （7）局部算法优化计算相位分布：利用光线追迹求得整形元件和相位校正元件的初始相位分布 和 后，将其作为局部算法的初始量，通过GS算法优化计算，最后求得优化后的相位分布，即完成了利用光线追迹法和局部算法求得整形元件和相位校正元件的相位，实现了高斯光束整形为平顶光束。
[0017] 发明效果：
[0018] 本发明利用光线追迹法得到了整形元件和相位校正元件的初始相位，然后通过将初始相位代入局部算法中得到了整形元件和相位校正元件的优化相位分布，整形元件和相位校正元件的组合能够实现将高斯光束整形为平顶光束，整形后的平顶光束的均匀性理论上能接近100%，满足一些特殊应用。而且局部算法计算所用时间大约为全局算法计算所用时间的百分之一，同时把光线追迹得到的初始相位带入局部算法后，最后得到的优化相位分布更有利于加工。
附图说明
[0019] 图1是本发明流程图；
[0020] 图2是具体实施方式一中高斯光束整形为平顶光束的结构图，1为整形元件，2为相位校正元件；
[0021] 图3是具体实施方式一中高斯光束整形为平顶光束的原理图；
[0022] 图4是具体实施方式一中的光程函数与光线偏角间的关系图；
[0023] 图5是具体实施方式四中的高斯光束的波形图；
[0024] 图6是具体实施方式四中的平顶光束的波形图。
具体实施方式
[0025] 具体实施方式一：结合图1～图4说明本实施方式，本实施方式的一种用于实现高斯光束整形为平顶光束的方法按以下步骤实现：
[0026] （1）根据图2确定高斯光束的直径D1，平顶光束的直径D2，以及两个光学元件的间距L：如图3，若激光高斯光束的半径为R1，则D1=2R1；如图4，若需求的平顶光束的半径为R2，则D2=2R2；为保证经过光学整形元件后的光线的偏角尽可能的小，一般取L≥D1；
[0027] （2）根据能量守恒计算平顶光束的光强；
[0028] （3）根据光线追迹建立一一对应关系；
[0029] （4）计算光线偏角：由几何关系，可以确定
[0030]
[0031] （5）确定整形元件的初始相位分布：如图6，由光程函数与光线偏角间的关系并结合图5，可以得到
[0032]
[0033] 对上式积分，得到整形元件的初始相位分布；
[0034] （6）确定相位校正元件的初始相位分布；
[0035] （7）局部算法优化计算相位分布：利用光线追迹求得整形元件和相位校正元件的初始相位分布 和 后，将其作为局部算法的初始量，通过算法优化计算，最后求得优化后的相位分布。
[0036] 本实施方式效果：
[0037] 本实施方式利用光线追迹法得到了整形元件和相位校正元件的初始相位，然后通过将初始相位代入局部算法中得到了整形元件和相位校正元件的优化相位分布，整形元件和相位校正元件的组合能够实现将高斯光束整形为平顶光束，整形后的平顶光束的均匀性理论上能接近100%，满足一些特殊应用。而且局部算法计算所用时间大约为全局算法计算所用时间的百分之一，同时把光线追迹得到的初始相位带入局部算法后，最后得到的优化相位分布更有利于加工。
[0038] 具体实施方式二：本实施方式与具体实施方式一不同的是：步骤（2）中：平顶光束的光强I2可由下式确定
[0039]
[0040] 式中，I1为高斯光束的中心光强，ω1为高斯光束的束腰半径。
[0041] 其它步骤及参数与具体实施方式一相同。
[0042] 具体实施方式三：本实施方式与具体实施方式一或二不同的是：步骤（3）中根据光线追迹建立一一对应关系具体为：
[0043] P1面上从距离光束中心r1处发出的光线入射到P2面上距离光束中心r2处，其中0≤r1≤D1，D1≤r2≤D2，且满足关系
[0044]
[0045] 利用上式可以将r2写成关于r1的表达式。
[0046] 其它步骤及参数与具体实施方式一或二相同。
[0047] 具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：步骤（6）中确定相位校正元件的初始相位分布具体为：经过整形元件调制后的光束在位相校正元件的前表面P2上的光强分布为平顶的，但是传播方向不再沿原光轴传播，需要通过位相校正元件的调制，使得整形后的光束沿着原光轴传播，由光程函数与光线偏角间的关系可以得到[0048]
[0049] 对上式积分，得到相位校正元件的初始相位分布。
[0050] 其它步骤及参数与具体实施方式一至三之一相同。
[0051] 具体实施方式五：本实施方式与具体实施方式一至四之一不同的是：步骤（7）中GS算法优化计算具体为：
[0052] 一、选取DOE的位相初始值 与输入面（DOE所在平面）上的已知振幅分布A0(x,y)组成输入面复振幅分布U1(x,y)；
[0053] 二、对U1(x,y)作光学正变换，在输出面得到再现像的场分布U＇2(u,v)；
[0054] 三、为使输出面的再现像不断逼近目标图像，保持再现像的位相分布不变，用目标图像的已知振幅分布B0(u,v)替代再现像的振幅分布|U＇2(u,v)|；
[0055] 四、对经过修正后的再现像的场分布U2(u,v)做光学逆变换，在输入面得到场分布U＇1(x,y)；
[0056] 五、返回第一步，以U＇1(x,y)的位相部分取代 作下一次迭代。其它步骤及参数与具体实施方式一至四之一相同。
