Crosstalk correction of array fibers applied to beam quality measurement

应用于光束质量测量的阵列光纤串扰校正 Crosstalk correction of array fibers applied to beam quality measurement 1 2 3 first-author 罗　杰（1996-），女，安徽池州人，博士研究生，2018年于广西大学获得学士学位，主要从事激光大气传输及参数测量技术的研究。E-mail： jiebao517@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246391&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246405&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246393&type= 22.01332855 32.00399780 罗　杰 （1996-），女，安徽池州人，博士研究生，2018年于广西大学获得学士学位，主要从事激光大气传输及参数测量技术的研究。E-mail： jiebao517@163.com jiebao517@163.com 2 3 corresp E-mail： laqin@aiofm.ac.cn E-mail： laqin@aiofm.ac.cn 秦来安（1983-），男，河南南阳人，博士，副研究员，2006年于郑州大学获得学士学位，2011年于中国科学院合肥物质科学研究院获得博士学位，主要从事激光大气传输及参数测量技术的研究。E-mail： laqin@aiofm.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246395&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246400&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246398&type= 22.01333618 32.00399780 秦来安 （1983-），男，河南南阳人，博士，副研究员，2006年于郑州大学获得学士学位，2011年于中国科学院合肥物质科学研究院获得博士学位，主要从事激光大气传输及参数测量技术的研究。E-mail： laqin@aiofm.ac.cn laqin@aiofm.ac.cn 2 3 1 2 3 2 3 中国科学技术大学 环境科学与光电技术学院，安徽 合肥 230026 中国科学技术大学 环境科学与光电技术学院，安徽 合肥 230026 School of Environmental Science and Optoelectronic Technology， University of Science and Technology of China， Hefei 230026， China School of Environmental Science and Optoelectronic Technology， University of Science and Technology of China， Hefei 230026， China 中国科学院合肥物质科学研究院安徽光学精密机械研究所 中国科学院大气光学重点实验室，安徽 合肥 230031 中国科学院合肥物质科学研究院安徽光学精密机械研究所 中国科学院大气光学重点实验室，安徽 合肥 230031 Key Laboratory of Atmospheric Optics，Anhui Institute of Optics and Fine Mechanics， Hefei Institutes of Physical Science，Chinese Academy of Sciences，Hefei 230031，China Key Laboratory of Atmospheric Optics，Anhui Institute of Optics and Fine Mechanics， Hefei Institutes of Physical Science，Chinese Academy of Sciences，Hefei 230031，China 先进激光技术安徽省实验室，安徽 合肥 230037 先进激光技术安徽省实验室，安徽 合肥 230037 Advanced Laser Technology Laboratory of Anhui Province， Hefei 230037， China Advanced Laser Technology Laboratory of Anhui Province， Hefei 230037， China 为降低光束质量测量系统中阵列光纤输出串扰导致的测量误差，本文提出了硬件设计与算法构建相结合的串扰校正方法，并就此展开了相关的原理研究与方法验证。首先，基于光束质量测量要求和光纤传光原理分析了阵列光纤的串扰影响；其次，结合朗伯散射原理和实验结果验证了校正串扰的原理：使用朗伯体降低阵列光纤输出光发散角差异，并建立统一的弥散光斑串扰模型，再采用反卷积算法复原到靶光斑；再次，介绍了针对性的反卷积算法原理，并就相关参数的设计展开了讨论；最后，对真实光斑、光纤输出的未校正光斑及已校正光斑进行对比分析，实验验证了校正方法的可行性。实验结果表明：与未校正光斑相比，校正后光斑强度分布的相对均方根误差由36.06%降至4.67%，桶中功率的相对均方根误差由7.79%降低至0.73%，86.5%桶中功率所在束宽的测量相对误差由10.83%降至3.46%，结合校正算法的图像处理和参数计算总时间约为8 s。 To reduce the error in a beam quality measurement system， which is caused by the output crosstalk of array fibers， a correction method combining hardware design and algorithm construction is proposed， and related principal research and method verification are conducted. First， the crosstalk effect of the array fibers is analyzed based on beam quality measurement requirements and the fiber transmission principle. Second， the principle of crosstalk correction is verified by combining the Lambertian scattering principle and experimental results： the Lambertian scattering principle is used to reduce the divergence angle difference of the output light of fibers， establishing a unified diffuse model to correct spot crosstalk， and the deconvolution algorithm is then used to restore the target spot. Third， the principle of a targeted deconvolution algorithm is introduced， and the design of the relevant algorithm parameters is discussed. Finally， a comparison and analysis of the real spot， the uncorrected spot， and the corrected spot output produced by the fiber are conducted， and the experiment is found to verify the feasibility and effectiveness of the correction method. The experimental results reveal that compared to the uncorrected spot， the relative root mean square error of the intensity distribution of the corrected spot is reduced from 36.06% to 4.67%， and the relative root mean square error of the power in the bucket is reduced from 7.79% to 0.73%. The relative error of the beam width where the power in the bucket is 86.5% is reduced from 10.83% to 3.46%. The total time of image processing and parameter calculation using correction algorithm is about 8 s. 光束质量测量 阵列光纤 串扰校正 反卷积 beam quality measurement array fiber crosstalk correction deconvolution 1 引　言 1 光束质量研究在激光器设计、激光传能、激光通信、激光探测等多领域都具有重要意义，其典型测量方法为取样光束，传输剖面光斑至探测元件，对响应值分布进行复原和校正得到强度分布，进一步计算相关参数完成质量评价 ［ 1 1 - 2 2 1-2 ］ 。相关的研究重点在于取样、传输和探测系统等硬件的设计，以及复原算法和参数计算的优化 ［ 3 3 - 6 6 3-6 ］ 。 近年来，阵列光纤由于突出的综合特性在光束质量测量系统中得到多方面应用。作为取样元件，光纤芯径小，分辨率高，且对于数值孔径（Numerical Aperture，NA）对应角度内入射的光可以全反射低损通过，实现多角度下的高一致性取样 ［ 7 7 - 8 8 7-8 ］ ；作为衰减元件，发散输出的光纤配合适当距离放置的探测器可以实现10 3 ~10 4 倍的衰减 ［ 9 9 ］ ；作为传光元件，高透过率的柔性光纤通过接收端稀疏、输出端密集的排布方式能够对大面积光斑实现缩束，进而使用外径毫米级的电感耦合器件（Charge Coupled Device，CCD）靶面在无需考虑成像畸变的前提下测量分米级的光束截面 ［ 10 10 - 11 11 10-11 ］ 。光束质量的有效测量基于光斑分布的准确获取，因此要求阵列光纤具有高一致性的传输效率和输出光场，以实现探测元件对各光纤输出光能的响应度一致，一般要求传输效率和响应度的非均匀性不超过6% ［ 12 12 - 14 14 12-14 ］ 。目前不同NA和透光率的光纤都已有成熟的制造工艺，可以保证传输效率的可控 ［ 15 15 ］ ；而光纤的输出光场较为复杂，一般表现为发散光束，易相互串扰；且发散角度和远场光斑分布受到光束入射角度 ［ 16 16 ］ 、光纤NA ［ 17 17 ］ 、弯曲曲率 ［ 18 18 ］ 等多方面因素的综合影响，导致串扰难以单纯使用算法进行预测和校正。增大单元间距可以有效降低串扰但限制了系统的取样分辨率；另外使用CCD相机对阵列光纤的输出端面进行成像，而非使用探测器直接采集发散光，在原理上可以忽略发散光的串扰。但在实际应用中发现，阵列柔性光纤在实现转向、缩束时难以达到弯曲程度的一致，导致各光纤输出发散角的差异，进一步造成相机拍摄的各光纤输出的单元光斑发生程度不一的弥散，在像面上相互串扰。 本文针对应用于光束质量测量的阵列光纤，结合硬件设计和算法处理提出一类串扰校正方法：在光纤输出端增设透射朗伯体，使用CCD相机拍摄朗伯体发光面，降低光纤输出发散角差异的影响；再配套相应反卷积算法，复原弥散光斑，最终测量出阵列光纤的真实输出，进而实现后续的光束质量测量。本文介绍了阵列光纤的串扰原理与影响，并基于朗伯散射原理和反卷积理论进行了串扰校正的仿真分析与实验研究：首先实验验证了所选材料的朗伯特性和作用；并以朗伯散射所成高斯光斑作为卷积核，建立相应反卷积算法，实现串扰的校正；再使用光线追迹的仿真方法选取最佳算法参数；最后对比分析了未经光纤传输的真实光斑、光纤输出未校正串扰的光斑及已校正串扰的复原光斑，实验验证了校正方法可以有效降低串扰引起的测量误差。 2 测量原理 1 2.1　阵列光纤串扰影响 2 为避免阵列光纤对截面光斑取样占空比及耦合效率过低，此应用一般使用芯径和可接收孔径角更大的多模光纤，其光传输示意图如 图1 图1 。输出光束的最大发散角 <math id="M1"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi><mo>,</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246438&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246429&type= 8.89000034 3.72533321 处于 <math id="M2"><mo stretchy="false">[</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">A</mi></mrow></msub><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246410&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246408&type= 13.63133430 4.48733330 ， <math id="M3"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246414&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246421&type= 3.47133350 3.72533321 为光束输入角度， <math id="M4"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">A</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246483&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246481&type= 4.99533367 3.72533321 为光纤NA所对应的半孔径角，远场光斑分布有高斯型、近平顶型和近环形等 ［ 19 19 ］ ， <math id="M5"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi><mo>,</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246428&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246434&type= 8.89000034 3.72533321 和远场光斑分布都会随着光纤弯曲应力的变化而变化。阵列中，各光纤所受应力和弯曲程度的差异导致各输出光束的 <math id="M6"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi><mo>,</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246438&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246429&type= 8.89000034 3.72533321 在 <math id="M7"><mo stretchy="false">[</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">A</mi></mrow></msub><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246448&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246431&type= 13.63133430 4.48733330 内分布，即使在光束直入射的情况下，弯曲光纤的输出发散角度仍有可能达到 <math id="M8"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">A</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246483&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246481&type= 4.99533367 3.72533321 。拍摄时，一方面发散角的差异导致了相机视场光阑对各光纤输出光所限的入光量不同，造成拍摄失真；另一方面各光纤输出光的焦面无法统一，离焦产生的弥散斑在像面上相互串扰，且弥散程度不一，弥散量不易确定，在池化处理中统计所拍图像的强度分布时，存在难以单纯使用算法消除的误差 ［ 20 20 ］ 。 2.2　朗伯散射特性 2 针对2.1节介绍的应用阵列光纤时产生的串扰，提出在光纤输出端面增设密接透射朗伯体的方法，使阵列光纤输出的不同发散角光束都转化为遵循 式（1） 式（1） 分布的朗伯辐射光 ［ 21 21 ］ 。 <math display="block" id="M9"><msub><mrow><mi>I</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo>=</mo><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">N</mi></mrow></msub><mo>⋅</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246467&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246465&type= 19.55800056 3.72533321 ， 其中， <math id="M10"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">N</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246470&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246458&type= 2.87866688 3.72533321 为入射光透过朗伯体后在法线方向的输出光强， <math id="M11"><msub><mrow><mi>I</mi></mrow><mrow><mi>θ</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246474&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246472&type= 2.37066650 3.72533321 则为和法线成 <math id="M12"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246493&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246477&type= 1.77800000 2.79399991 角方向的光强，由 式（1） 式（1） 可知朗伯辐射光的分布与入射角度无关。当朗伯体端面与光纤输出端面密接时，即使光纤输出光的发散角不同，但其尚未在空间中自由发散就传输至朗伯体中，并转为与光纤输出发散角无关的朗伯辐射光，从而可以消除光纤输出发散角差异的影响。 式（1） 式（1） 为理想朗伯体的特性，而在实际应用中，还需要对所用材料进行朗伯散射特性的实验验证。本方法所用材料由透过型体散射玻璃加工而成，材料表面进行精密抛光，采用边缘压紧的安装方式与光纤输出端面保持密接。为进行朗伯特性验证，测量不同入射角度下材料的发光面光斑，如分布保持稳定，即满足消除光纤发散角差异的要求。考虑到此类应用中常用光纤NA一般不大于0.6，其对应 <math id="M13"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">A</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246483&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246481&type= 4.99533367 3.72533321 和 <math id="M14"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi><mo>,</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246497&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246495&type= 8.89000034 3.72533321 为36.8°，因此使用相对角度在［ - 40°，40°］范围内变化的1 064 nm激光束模拟光纤输出的不同角度发散光，入射至厚度为0.12 mm的材料的同一位置，使用CCD相机拍摄材料发光面光斑，获得光斑总模拟数字单元（Analog Digital Unit，ADU）及光斑分布随入射角度的变化分别如 图2 图2 和 图3 图3 所示。 图2 图2 中，统计光斑的总ADU，发现其波动范围不超过4.00%。 图3 图3 中，光斑剖线基本一致，各角度下的剖线与正入射之间的相对均方根误差（Root Mean Square Error，RMSE ）不超过平均值的9.00%，可认为所选材料的发光面光斑分布基本不受入射角度的影响，具有所需朗伯特性。 2.3　测量原理验证 2 为验证并量化阵列光纤的输出串扰及增设密接朗伯体后的效果，使用CCD相机拍摄单纤的输出端面，实验系统如 图4 图4 所示，此时不增设朗伯体。由于S型弯曲为光纤阵列应用中较为常见的弯曲方式，因此在实验中将光纤穿入管道中实现S型弯曲；换用管道以改变曲率半径 <math id="M15"><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246544&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246535&type= 2.70933342 2.79399991 ， <math id="M16"><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246544&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246535&type= 2.70933342 2.79399991 的取值分别为纤芯半径（ R 0 ）的120倍、90倍和30倍，以模拟柔性光纤在自由状态下不可避免的排布差异；光源及拍摄环境不变。实验结果如 图5 图5 （a）所示，对比不同 <math id="M17"><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246540&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246538&type= 2.70933342 2.79399991 对应的相机所摄的输出端面，其存在不同程度的弥散，叠加三幅图像可以发现其相互之间发生串扰，对每个光纤端面划分网格并统计网格之外的总ADU即串扰至相邻光纤的能量，分别占光纤总输出的10.43%（120倍）、12.84%（90倍）和17.76%（30倍），验证了直接拍摄光纤端面时，光纤排布差异将导致光纤之间的相互串扰不一致，且随着光纤的逐渐拉紧，串扰变得更加严重。为验证朗伯体的作用，在此实验系统中增设与光纤输出端面密接的朗伯体，使用相机拍摄朗伯体发光表面，改变 <math id="M18"><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246544&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246535&type= 2.70933342 2.79399991 后获得光斑再进行叠加得到 图5 图5 （b），统计得到串扰能量占每根光纤总输出的49.48%（120倍）、48.72%（90倍）和49.07%（30倍）。朗伯体的散射使得串扰增加，但大幅降低了串扰的差异性，使所摄光斑遵循统一的弥散模型，为消除串扰的算法处理建立前提。 3 校正算法 1 3.1　反卷积校正原理 2 从 图3 图3 可知，朗伯体发光面光斑呈高斯分布，可将光斑看作经高斯模糊的卷积处理结果，卷积核可通过实验获得。当经朗伯体作用后的输出光斑和卷积核都已知时，通过反卷积处理可解得阵列光纤的真实输出。具体求解如下：其中 <math id="M19"><mi mathvariant="bold-italic">I</mi><mo stretchy="false">(</mo><mi>m</mi><mo>×</mo><mi>n</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246588&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246586&type= 14.81666565 3.55599999 表示朗伯体相对输入即光纤真实输出； <math id="M20"><mi mathvariant="bold-italic">O</mi><mo stretchy="false">(</mo><mi>m</mi><mo>×</mo><mi>n</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246593&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246590&type= 16.08666801 3.55599999 表示朗伯体相对输出， <math id="M21"><mi>m</mi><mo>×</mo><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246580&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246578&type= 9.56733322 2.79399991 代表 <math id="M22"><mi>m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247014&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247013&type= 2.79399991 2.79399991 行 <math id="M23"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247016&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247015&type= 1.86266661 2.79399991 列的阵列光纤。 （1）求解高斯型卷积核 <math id="M24"><mi>G</mi><mo stretchy="false">(</mo><mi>l</mi><mo>×</mo><mi>l</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246622&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246620&type= 13.80066681 3.55599999 ， <math id="M25"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 代表卷积核尺寸，由于后续将使用中心对齐的same模式，一般将 <math id="M26"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 设置为奇数，考虑到实际拍摄画面的幅度， <math id="M27"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 一般不超过（ <math id="M28"><mn mathvariant="normal">2</mn><mi>m</mi><mo>-</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246642&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246640&type= 11.43000031 2.79399991 ）和（ <math id="M29"><mn mathvariant="normal">2</mn><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246661&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246652&type= 10.49866676 2.79399991 ）： <math display="block" id="M30"><mi>G</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi>π</mi><msup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">(</mo><mo>-</mo><mo stretchy="false">(</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo><mo>/</mo><mn mathvariant="normal">2</mn><msup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246673&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246656&type= 57.99666977 8.12800026 ， 其中，由于卷积核只表征单元强度，因此其积分为1，只求解尺度参数 <math id="M31"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 ，求解方式为获取实验光斑进行二维高斯拟合。 （2）获取输出矩阵 <math id="M32"><mi mathvariant="bold-italic">O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246808&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246806&type= 2.96333337 2.79399991 ，对实验所摄光斑需要进行池化处理，将分辨率由像素数转为光纤数；池化网格的定位需根据光纤位置预先测量。 （3）获取卷积关系，使用 <math id="M33"><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246932&type= 2.79399991 2.79399991 对 <math id="M34"><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246836&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246835&type= 1.43933344 2.79399991 进行遍历，当 <math id="M35"><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246932&type= 2.79399991 2.79399991 的尺寸 <math id="M36"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 不足以访问 <math id="M37"><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246836&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246835&type= 1.43933344 2.79399991 中的每个元素时，对 <math id="M38"><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246932&type= 2.79399991 2.79399991 进行补零即可：构建 <math id="M39"><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mi>m</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246733&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246716&type= 14.22399998 3.55599999 行 <math id="M40"><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246736&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246723&type= 13.29266548 3.55599999 列的零矩阵，将 <math id="M41"><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246932&type= 2.79399991 2.79399991 以中心位置对齐的same模式填入其中，得到矩阵 <math id="M42"><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246753&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246747&type= 3.72533321 3.72533321 。 <math id="M43"><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246836&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246835&type= 1.43933344 2.79399991 、 <math id="M44"><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246753&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246747&type= 3.72533321 3.72533321 和 <math id="M45"><mi mathvariant="bold-italic">O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246808&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246806&type= 2.96333337 2.79399991 存在以下卷积关系： <math display="block" id="M46"><mtable><mtr><mtd><mi mathvariant="bold-italic">O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="bold-italic">I</mi><mo>·</mo><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo>=</mo></mtd></mtr><mtr><mtd><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>m</mi></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>j</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>n</mi></mrow></munderover></mstyle><mrow><mo stretchy="false">[</mo><mi mathvariant="bold-italic">I</mi><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo><mo>⋅</mo><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mi>m</mi><mo>+</mo><mi>x</mi><mo>-</mo><mi>i</mi><mo>,</mo><mi>n</mi><mo>+</mo><mi>y</mi><mo>-</mo><mi>j</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow></mrow></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246771&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246757&type= 64.00800323 14.81666565 . （4）反卷积求解输入矩阵 <math id="M47"><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246836&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246835&type= 1.43933344 2.79399991 。为使用矩阵求逆的方法方便快捷地实现反卷积，构建矩阵表达 式（4） 式（4） 以求解 式（3） 式（3） ： <math display="block" id="M48"><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">G</mi><mo stretchy="false">(</mo><mi>m</mi><mi>n</mi><mo>×</mo><mi>m</mi><mi>n</mi><mo stretchy="false">)</mo><mo>·</mo><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">I</mi><mo stretchy="false">(</mo><mi>m</mi><mi>n</mi><mo>×</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="bold-italic">O</mi><mi mathvariant="bold-italic">O</mi><mo stretchy="false">(</mo><mi>m</mi><mi>n</mi><mo>×</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246794&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246776&type= 68.07199860 3.55599999 ， 其中， <math id="M49"><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246829&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246820&type= 2.87866688 2.79399991 和 <math id="M50"><mi mathvariant="bold-italic">O</mi><mi mathvariant="bold-italic">O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246802&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246800&type= 5.84200001 2.79399991 分别为 <math id="M51"><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246836&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246835&type= 1.43933344 2.79399991 和 <math id="M52"><mi mathvariant="bold-italic">O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246808&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246806&type= 2.96333337 2.79399991 按行变形后得到的一维矩阵，据 式（3） 式（3） 求解出卷积核组合矩阵 <math id="M53"><mi mathvariant="bold-italic">G</mi><mi mathvariant="bold-italic">G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246812&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246810&type= 5.41866684 2.79399991 ，可得具体表达式如下： <math display="block" id="M54"><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mn mathvariant="normal">1,1</mn><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd/><mtd><mo>⋮</mo></mtd><mtd/><mtd><mo>⋮</mo></mtd><mtd/><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mn mathvariant="normal">2</mn><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mn mathvariant="normal">1,2</mn><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd/><mtd><mo>⋮</mo></mtd><mtd/><mtd><mo>⋮</mo></mtd><mtd/><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mi>m</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mi>m</mi><mo>-</mo><mn mathvariant="normal">1,1</mn><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd/><mtd><mo>⋮</mo></mtd><mtd/><mtd><mo>⋮</mo></mtd><mtd/><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mi>m</mi><mo>-</mo><mn mathvariant="normal">1,2</mn><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mi>m</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mn mathvariant="normal">2</mn><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mtd><mtd><mo>…</mo></mtd><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">z</mi></mrow></msub><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mtd></mtr></mtable></mrow></mfenced><mo>·</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1,1</mn><mo stretchy="false">)</mo></mrow></msub></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mrow></msub></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msub></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">O</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1,1</mn><mo stretchy="false">)</mo></mrow></msub></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">O</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mrow></msub></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">O</mi></mrow><mrow><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msub></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">O</mi></mrow><mrow><mo stretchy="false">(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246823&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246814&type= 154.68600464 37.93066406 . 求解 <math id="M55"><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">I</mi><mo>=</mo><mi mathvariant="bold-italic">G</mi><msup><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mo>·</mo><mi mathvariant="bold-italic">O</mi><mi mathvariant="bold-italic">O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246818&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246816&type= 23.36800003 3.04800010 ，变形 <math id="M56"><mi mathvariant="bold-italic">I</mi><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246829&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246820&type= 2.87866688 2.79399991 则可得 <math id="M57"><mi>m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247014&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247013&type= 2.79399991 2.79399991 行 <math id="M58"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247016&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247015&type= 1.86266661 2.79399991 列的输入 <math id="M59"><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246836&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246835&type= 1.43933344 2.79399991 。 3.2　参数设计 2 本算法的主要参数设计围绕卷积核 <math id="M60"><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246932&type= 2.79399991 2.79399991 展开，主要参数为高斯分布的尺度参数 <math id="M61"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 和卷积核尺寸 <math id="M62"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 ，二者共同决定了校正的精准度。 3 . 2 . 1　高斯分布尺度参数 <math id="M63"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 3 在朗伯体的材料确定之后， <math id="M64"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 一般由材料厚度决定，进行单光纤与材料的组合实验及相应仿真，得到不同厚度材料的 <math id="M65"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 的实验值与仿真值如 表1 表1 所示。 Thickness /mm Experimental value Simulation value Relative error 0.12 12.22 11.43 6.46% 0.25 23.81 22.07 7.31% 0.38 34.15 32.88 3.72% 0.50 46.95 45.51 3.07% 表1 表1 显示朗伯体越厚，尺度参数 <math id="M68"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 越大，则发光面的弥散光斑越大。原因是所用材料为体散射材料，其散射效果可视为多片薄的面散射材料叠加，因此弥散程度与材料厚度成正相关。在实际应用中一般使用预实验所获取的 <math id="M69"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 作为卷积核 <math id="M70"><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246932&type= 2.79399991 2.79399991 的尺度参数而非仿真值，而本文中希望采用建模仿真方法讨论参数选取的规律，因此还在 表1 表1 中对比了 <math id="M71"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 的实验值与仿真值，发现其误差不超过实验值的8%，认为此应用下的光传输仿真具有可信度。获取尺度参数 <math id="M72"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 及卷积核 <math id="M73"><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246932&type= 2.79399991 2.79399991 后，仿真了阵列光纤的输出在不同厚度材料下的校正效果，结果如 图6 图6 所示，仿真条件如 表2 表2 所示。 Source file format File size/kB Boundary condition seq， txt 35，157 <math id="M74"><mi>m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246922&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246940&type= 2.45533323 2.37066650 ， <math id="M75"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246924&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246923&type= 1.60866666 2.37066650 =20 图6 图6 （a）、（d）、（g）和（j）表示阵列光源，其中心间距为光源直径1.5倍，对400单元的光源强度整体分布赋值为高斯分布，模拟光束质量测量中高频出现的高斯光斑经过阵列光纤后的输出光。 图6 图6 （b）、（e）、（h）和（k）为材料发光面未经过反卷积处理的池化光斑，随着材料厚度的增加，单纤输出光斑愈加弥散，其构成的整体光斑也逐渐散开，这会造成光束质量参数测量的误差，因此需要进行反卷积校正。由 图6 图6 （c）、（f）、（i）和（l）可知，随着弥散程度的增加，算法的稳定性有所下降，主要体现在背景噪声愈加明显。在 图6 图6 （l）中，当材料厚度为0.50 mm时，背景噪声的起伏已达光斑分布峰值的4.3%和平均值的50.4%，对于光斑分布统计和光束质量参数计算的影响已无法忽视，因此在满足朗伯特性的前提下，此应用一般选取最薄的材料。由于2.2节已验证0.12 mm材料的朗伯特性的材料，结合此节的讨论及工艺限制，应用此现有最薄材料进行后续实验，此时参考 表1 表1 ，对应 <math id="M76"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 为12.22。 3 . 2 . 2　卷积核尺寸 <math id="M77"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 3 在确定了尺度参数 <math id="M78"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 即材料厚度之后，对卷积核尺寸 <math id="M79"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 进行讨论， <math id="M80"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 一般取值为 <math id="M81"><mo stretchy="false">[</mo><mn mathvariant="normal">3,2</mn><mi>m</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246973&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246956&type= 18.37266731 3.55599999 中的奇数，光纤行（列）数 <math id="M82"><mi>m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247014&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247013&type= 2.79399991 2.79399991 在此节中取20， <math id="M83"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 只取奇数是利于反卷积处理时的same模式的中心对齐。对0.12 mm材料发光面的光斑进行不同大小卷积核的反卷积处理，为便于展示校正后光斑与真实光斑的对比，使用两类光斑的相对RMSE来表征其差异，相对RMSE越接近于零则差异越小，校正效果越优。分别统计了光斑强度分布的相对RMSE及光斑桶中功率（Power in the Bucket， PIB）曲线的相对RMSE，前者可以展示光斑实时显示的准确度，后者以相应半径所在圆内功率的差异显示计算光束质量参数的准确度，两者共同表征了校正效果的优劣；为了讨论卷积核尺寸的选取规律，还增加了卷积核的PIB作为参照，用以显示相应大小的卷积核内所占功率，如 图7 图7 所示。 由 图7 图7 可以发现：随着 <math id="M88"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 的增加，两类RMSE逐渐降低，表明校正有效性的提升；在卷积核 l 在［17，39］的区间内，光斑强度分布的RMSE（红色点线）在7.09%到7.28%之间变化，起伏不超过0.2%；PIB的RMSE（蓝色点线）在0.13%到0.15%之间变化，起伏不超过0.02%。此应用场景下，综合考虑校正算法的实时运行时间与校正效果， <math id="M89"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 可取17。同时对比卷积核的PIB曲线可以注意到， <math id="M90"><mi>l</mi><mo>=</mo><mn mathvariant="normal">17</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246995&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246994&type= 9.73666668 2.79399991 所对应的PIB为0.993 0，而 <math id="M91"><mi>l</mi><mo>=</mo><mn mathvariant="normal">16</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246997&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246996&type= 9.73666668 2.79399991 所对应的PIB为0.989 7，因此当推广至其它的卷积核尺度参数 <math id="M92"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30246998&type= 1.69333339 2.79399991 与光纤行（列）数 <math id="M93"><mi>m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247014&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247013&type= 2.79399991 2.79399991 的应用时，可以选取卷积核的PIB为99%所对应的就近整数半径为最佳卷积核尺寸 <math id="M94"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247002&type= 1.18533325 2.79399991 。 4 实验验证 1 4.1　实验平台 2 为了测试应用于光束质量测量的阵列光纤串扰校正的效果，使用了自主研发的光束质量测量系统，如 图8 图8 所示，该测量系统利用光纤阵列两端面不同疏密度的锥形排列，对截面光斑进行了40倍缩束以供相机进行无畸变拍摄。实验采用1 064 nm光源入射系统的接收端面即光纤疏排端面，由20行20列的矩形光纤阵列取样并传输至输出端面即光纤密排端面，经朗伯散射、衰减、滤光等处理后被CCD相机拍摄，所摄光斑经序列文件矩阵化、背景噪音过滤、光纤对应网格池化、校正算法处理、插值扩束、功率复原等实时显示为光斑强度分布，并可以开展后续的光束质量参数的计算。 4.2　实验结果 2 在本实验中，为了验证校正方法的效果，使用 图8 图8 所示平台展开了已校正串扰的光斑测量实验；并设计了未校正串扰的对比实验，实验条件为：实验平台拆卸朗伯体，相机直接拍摄输出端面，且不对光斑进行反卷积处理；同时还使用漫反射法测量结果作为真实光斑，实验条件为：接收端面同一处换置漫反射屏，使用大视场CCD相机拍摄屏上成像光斑。分别将未校正和已校正的光斑与真实光斑进行对比分析，得到的处理过程与结果如 图9 图9 所示，PIB曲线对比如 图10 图10 所示。 从 图9 图9 的强度分布测量中可以发现，未校正光斑（d）在强度分布上与真实光斑（a）存在着较为明显的差距，量化两者的数值矩阵并计算（d）与（a）的相对RMSE为36.06%，而经过朗伯散射与反卷积的已校正光斑（h）则与（a）基本吻合，其相对RMSE降低至4.67%。对 图9 图9 的强度分布进行光束质量参数的计算，得到如 图10 图10 所示的分析结果。在典型束宽计算采用的PIB=86.5%处 ［ 22 22 ］ ，未校正光斑（d）的对应束宽为60.80 mm，与真实光斑（a）束宽54.86 mm相差了5.94 mm，差值占（a）束宽的10.83%；而已校正光斑（h）的束宽为56.76 mm，与（a）差值为1.90 mm，差值占真实光斑束宽的3.46%，将束宽测量的误差降低了7.37%；另外未校正光斑（d）与（a）的PIB曲线的相对RMSE为7.79%，而已校正光斑（h）的为0.73%，认为校正方法可以有效降低阵列光纤串扰带来的测量误差。光斑（h）与（a）之间的差距除了校正方法本身存在的误差外，还包含了系统自带误差，如光纤阵列透过率差异未完全标定补偿等。但考虑到总误差已基本达到光束质量测量的要求 ［ 11 11 - 13 13 11-13 ］ ，可以认为校正方法能够较高精度地复原到靶光斑，进而可以顺利展开后续的光束质量测量。同时光束质量测量还要求算法的快速运行，当阵列行列数 <math id="M95"><mi>m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247014&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247013&type= 2.79399991 2.79399991 和 <math id="M96"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247016&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247015&type= 1.86266661 2.79399991 由10增至50时，本校正算法的运行时间由175 ms增至398 ms。考虑到CCD相机拍摄的100帧、 <math id="M97"><mn mathvariant="normal">1</mn><mtext> </mtext><mn mathvariant="normal">000</mn><mo>×</mo><mn mathvariant="normal">1</mn><mtext> </mtext><mn mathvariant="normal">000</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247018&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=30247017&type= 21.25133324 2.87866688 像素的seq源文件进行光斑复原处理及质量参数计算的时间一般不低于8 000 ms，校正算法运行时间一般不超过总运行时间的5%，对总测量速度影响较小。 5 结　论 1 本文基于光束质量测量系统，研究了阵列光纤的输出串扰，并提出了一类串扰校正方法。结果表明：阵列光纤的输出串扰造成了光束质量测量的误差，主要表现在单元光纤输出值与测量值响应度不一致，及其造成的光斑分布和PIB曲线偏差；针对这一问题设计了增设朗伯体并结合反卷积处理的串扰校正方法，并展开了相关参数设计的研究：一般使用相对较薄的朗伯体，可以有效降低反卷积处理后的背景噪声；同时考虑卷积核尺寸时，一般选取卷积核PIB为99%所对应的尺寸，最后实验验证了校正方法的有效性。在光斑分布的测量中，校正处理使光斑分布测量的相对RMSE得到了21.67%的降低；在PIB的测量中则得到了3.41%的降低；另外在PIB=86.5%处，串扰影响使得相应束宽的测量产生了10.83%的相对误差，校正处理使得误差降低了7.37%，而校正运算时间一般不超过总运算时间的5%；综上认为校正方法快速、有效降低了阵列光纤串扰对光束质量测量造成的误差，为此类测量提供了可用的优化方向。 参考文献 JORGE K C ， RIVA R ， RODRIGUES N A S ， et al . Scattered light imaging method （SLIM） for characterization of arbitrary laser beam intensity profiles ［J］. Applied Optics ， 2014 ， 53 （ 20 ）： 4555 - 4564 . doi: 10.1364/ao.53.004555 http://dx.doi.org/10.1364/ao.53.004555 Scattered light imaging method （SLIM） for characterization of arbitrary laser beam intensity profiles 李晓明 ， 沈学举 ， 李刚 ， 等 . 基于光斑轮廓特征的激光光束质量简易评价方法 ［J］. 中国激光 ， 2015 ， 42 （ 5 ）： 0502001 . doi: 10.3788/cjl201542.0502001 http://dx.doi.org/10.3788/cjl201542.0502001 基于光斑轮廓特征的激光光束质量简易评价方法 LI X M ， SHEN X J ， LI G ， et al . Evaluating specifications of laser beam quality based on beam profile ［J］. Chinese Journal of Lasers ， 2015 ， 42 （ 5 ）： 0502001 . （in Chinese） . doi: 10.3788/cjl201542.0502001 http://dx.doi.org/10.3788/cjl201542.0502001 Evaluating specifications of laser beam quality based on beam profile SUMACHEV K E ， SHARKOV V V ， SAVIKIN A P ， et al . Measurement of 194 µm YAP： Tm 3+ laser beam quality using a CCD camera ［J］. Journal of Optical Technology ， 2020 ， 87 （ 5 ）： 262 . doi: 10.1364/jot.87.000262 http://dx.doi.org/10.1364/jot.87.000262 Measurement of 194 µm YAP： Tm3+ laser beam quality using a CCD camera 冯国斌 ， 杨鹏翎 ， 王振宝 ， 等 . 光电/量热复合式近红外高能激光光斑探测器 ［J］. 光学 精密工程 ， 2013 ， 21 （ 6 ）： 1417 - 1424 . doi: 10.3788/OPE.20132106.1417 http://dx.doi.org/10.3788/OPE.20132106.1417 光电/量热复合式近红外高能激光光斑探测器 FENG G B ， YANG P L ， WANG Z B ， et al . Photoelectric-calorimetric compound beam profile detector for near-infrared high energy laser ［J］. Opt. Precision Eng. ， 2013 ， 21 （ 6 ）： 1417 - 1424 . （in Chinese） . doi: 10.3788/OPE.20132106.1417 http://dx.doi.org/10.3788/OPE.20132106.1417 Photoelectric-calorimetric compound beam profile detector for near-infrared high energy laser 王艳茹 ， 王建忠 ， 冉铮惠 ， 等 . 高能激光光束质量β因子的影响因素分析 ［J］. 中国光学 ， 2021 ， 14 （ 2 ）： 353 - 360 . doi: 10.37188/CO.2020-0137 http://dx.doi.org/10.37188/CO.2020-0137 高能激光光束质量β因子的影响因素分析 WANG Y R ， WANG J Z ， RAN Z H ， et al . Analysis of effects on the beam quality β factor of high power laser ［J］. Chinese Optics ， 2021 ， 14 （ 2 ）： 353 - 360 . （in Chinese） . doi: 10.37188/CO.2020-0137 http://dx.doi.org/10.37188/CO.2020-0137 Analysis of effects on the beam quality β factor of high power laser 吴晗平 . 激光光束质量的评价与应用分析 ［J］. 光学 精密工程 ， 2000 ， 8 （ 2 ）： 128 - 132 . doi: 10.3321/j.issn:1004-924X.2000.02.003 http://dx.doi.org/10.3321/j.issn:1004-924X.2000.02.003 激光光束质量的评价与应用分析 WU H P . Evaluation and applied analysis of laser beam quality ［J］. Opt. Precision Eng. ， 2000 ， 8 （ 2 ）： 128 - 132 . （in Chinese） . doi: 10.3321/j.issn:1004-924X.2000.02.003 http://dx.doi.org/10.3321/j.issn:1004-924X.2000.02.003 Evaluation and applied analysis of laser beam quality 罗杰 ， 秦来安 ， 侯再红 ， 等 . 激光光斑分布测量系统中光纤传光特性 ［J］. 光学学报 ， 2021 ， 41 （ 12 ）： 1214002 . doi: 10.3788/aos202141.1214002 http://dx.doi.org/10.3788/aos202141.1214002 激光光斑分布测量系统中光纤传光特性 LUO J ， QIN L A ， HOU Z H ， et al . Light-transmitting characteristics of optical fiber in laser spot distribution measurement system ［J］. Acta Optica Sinica ， 2021 ， 41 （ 12 ）： 1214002 . （in Chinese） . doi: 10.3788/aos202141.1214002 http://dx.doi.org/10.3788/aos202141.1214002 Light-transmitting characteristics of optical fiber in laser spot distribution measurement system 管雯璐 ， 谭逢富 ， 侯再红 ， 等 . 高功率密度激光宽角度阵列探测技术研究 ［J］. 光学学报 ， 2022 ， 42 （ 2 ）： 159 - 166 . 高功率密度激光宽角度阵列探测技术研究 GUAN W L ， TAN F F ， HOU Z H ， et al . Wide angle array detection technology for high power density laser ［J］. Acta Optica Sinica ， 2022 ， 42 （ 2 ）： 159 - 166 . （in Chinese） Wide angle array detection technology for high power density laser 吴勇 ， 杨鹏翎 ， 陈绍武 ， 等 . 基于光纤阵列的激光光斑测量方法 ［J］. 强激光与粒子束 ， 2013 ， 25 （ 2 ）： 276 - 280 . doi: 10.3788/HPLPB20132502.0276 http://dx.doi.org/10.3788/HPLPB20132502.0276 基于光纤阵列的激光光斑测量方法 WU Y ， YANG P L ， CHEN S W ， et al . Laser intensity distribution measuring method based on fiber array ［J］. High Power Laser and Particle Beams ， 2013 ， 25 （ 2 ）： 276 - 280 . （in Chinese） . doi: 10.3788/HPLPB20132502.0276 http://dx.doi.org/10.3788/HPLPB20132502.0276 Laser intensity distribution measuring method based on fiber array 黄伟 ， 马松山 ， 李晓芹 ， 等 . 复合式激光远场光斑分布定量测量技术研究 ［J］. 激光与红外 ， 2018 ， 48 （ 8 ）： 1031 - 1035 . doi: 10.3969/j.issn.1001-5078.2018.08.017 http://dx.doi.org/10.3969/j.issn.1001-5078.2018.08.017 复合式激光远场光斑分布定量测量技术研究 HUANG W ， MA S S ， LI X Q ， et al . Study on quantitative compound measurement of far field spot distribution of laser ［J］. Laser & Infrared ， 2018 ， 48 （ 8 ）： 1031 - 1035 . （in Chinese） . doi: 10.3969/j.issn.1001-5078.2018.08.017 http://dx.doi.org/10.3969/j.issn.1001-5078.2018.08.017 Study on quantitative compound measurement of far field spot distribution of laser GU D F ， WINKER B ， CHEN Y C ， et al . High energy laser target board apparatus ， US9423298 ［P］. 2016-8-23 . High energy laser target board apparatus PANG M ， RONG J ， ZHOU S ， et al . Array detector for high energy laser based on diffuse transmission sampling ［J］. The Review of Scientific Instruments ， 2014 ， 85 （ 1 ）： 013105 . doi: 10.1063/1.4861390 http://dx.doi.org/10.1063/1.4861390 Array detector for high energy laser based on diffuse transmission sampling 王振宝 ， 冯国斌 ， 杨鹏翎 ， 等 . 近红外激光光斑功率密度时空分布探测器 ［J］. 红外与激光工程 ， 2011 ， 40 （ 5 ）： 935 - 938 . doi: 10.3969/j.issn.1007-2276.2011.05.032 http://dx.doi.org/10.3969/j.issn.1007-2276.2011.05.032 近红外激光光斑功率密度时空分布探测器 WANG Z B ， FENG G B ， YANG P L ， et al . Detector array for measuring spatio-temporal distribution of near infrared laser power density ［J］. Infrared and Laser Engineering ， 2011 ， 40 （ 5 ）： 935 - 938 . （in Chinese） . doi: 10.3969/j.issn.1007-2276.2011.05.032 http://dx.doi.org/10.3969/j.issn.1007-2276.2011.05.032 Detector array for measuring spatio-temporal distribution of near infrared laser power density 陈洪耀 ， 张黎明 ， 杨本永 ， 等 . 漫射板反射式大直径激光功率-能量测量方法 ［J］. 强激光与粒子束 ， 2009 ， 21 （ 10 ）： 1441 - 1445 . 漫射板反射式大直径激光功率-能量测量方法 CHEN H Y ， ZHANG L M ， YANG B Y ， et al . Diffuser reflection method for measuring power/energy of far-field laser spot ［J］. High Power Laser and Particle Beams ， 2009 ， 21 （ 10 ）： 1441 - 1445 . （in Chinese） Diffuser reflection method for measuring power/energy of far-field laser spot 刘星洋 ， 翟尚礼 ， 潘望 ， 等 . 红外传能光纤的传能特性研究进展 ［J］. 激光与红外 ， 2020 ， 50 （ 8 ）： 907 - 913 . doi: 10.3969/j.issn.1001-5078.2020.08.002 http://dx.doi.org/10.3969/j.issn.1001-5078.2020.08.002 红外传能光纤的传能特性研究进展 LIU X Y ， ZHAI S L ， PAN W ， et al . Research progress on energy transfer properties of infrared transmitting energy fiber ［J］. Laser & Infrared ， 2020 ， 50 （ 8 ）： 907 - 913 . （in Chinese） . doi: 10.3969/j.issn.1001-5078.2020.08.002 http://dx.doi.org/10.3969/j.issn.1001-5078.2020.08.002 Research progress on energy transfer properties of infrared transmitting energy fiber 张美 ， 魏福利 ， 邱孟通 ， 等 . 不同入射角激光束在大芯径SI-PDF光纤中的传输损耗 ［J］. 光学 精密工程 ， 2009 ， 17 （ 10 ）： 2406 - 2410 . doi: 10.3321/j.issn:1004-924X.2009.10.009 http://dx.doi.org/10.3321/j.issn:1004-924X.2009.10.009 不同入射角激光束在大芯径SI-PDF光纤中的传输损耗 ZHANG M ， WEI F L ， QIU M T ， et al . Transmission attenuation of laser beams with different incident angles in large-core SI-PDF optical fibers ［J］. Opt. Precision Eng. ， 2009 ， 17 （ 10 ）： 2406 - 2410 . （in Chinese） . doi: 10.3321/j.issn:1004-924X.2009.10.009 http://dx.doi.org/10.3321/j.issn:1004-924X.2009.10.009 Transmission attenuation of laser beams with different incident angles in large-core SI-PDF optical fibers 周德春 ， 于凤霞 ， 谭芳 ， 等 . 大数值孔径传像光纤的制备及其光学性能研究 ［J］. 激光与光电子学进展 ， 2010 ， 47 （ 12 ）： 120605 . doi: 10.3788/LOP47.120605 http://dx.doi.org/10.3788/LOP47.120605 大数值孔径传像光纤的制备及其光学性能研究 ZHOU D C ， YU F X ， TAN F ， et al . Preparation and optical performance of large numerical aperture image transfer fiber ［J］. Laser & Optoelectronics Progress ， 2010 ， 47 （ 12 ）： 120605 . （in Chinese） . doi: 10.3788/LOP47.120605 http://dx.doi.org/10.3788/LOP47.120605 Preparation and optical performance of large numerical aperture image transfer fiber 谭中伟 ， 杨婧雅 ， 刘艳 ， 等 . 基于卷积神经网络的增敏型光纤弯曲传感器 ［J］. 光学 精密工程 ， 2020 ， 28 （ 7 ）： 1454 - 1461 . doi: 10.37188/OPE.20202807.1454 http://dx.doi.org/10.37188/OPE.20202807.1454 基于卷积神经网络的增敏型光纤弯曲传感器 TAN Z W ， YANG J Y ， LIU Y ， et al . Enhanced fiber optic bending senor based on convolutional neural network ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 7 ）： 1454 - 1461 . （in Chinese） . doi: 10.37188/OPE.20202807.1454 http://dx.doi.org/10.37188/OPE.20202807.1454 Enhanced fiber optic bending senor based on convolutional neural network GAMBLING W A ， PAYNE D N ， MATSUMURA H . Mode conversion coefficients in optical fibers ［J］. Applied Optics ， 1975 ， 14 （ 7 ）： 1538 - 1542 . doi: 10.1364/ao.14.001538 http://dx.doi.org/10.1364/ao.14.001538 Mode conversion coefficients in optical fibers VELSINK M C ， LYU Z P ， PINKSE P W H ， et al . Comparison of round- and square-core fibers for sensing， imaging， and spectroscopy ［J］. Optics Express ， 2021 ， 29 （ 5 ）： 6523 . doi: 10.1364/oe.417021 http://dx.doi.org/10.1364/oe.417021 Comparison of round- and square-core fibers for sensing， imaging， and spectroscopy 郁道银 ， 谈恒英 . 工程光学 ［M］. 3版 . 北京 ： 机械工业出版社 ， 2011 . 工程光学 YU D Y ， TAN H Y . Engineering optics ［M］. Beijing ： China Machine Press ， 2011 . （in Chinese） Engineering optics 张雨秋 ， 黄良金 ， 常琦 ， 等 . 光纤激光基模光束的β因子 ［J］. 物理学报 ， 2021 ， 70 （ 20 ）： 79 - 85 . doi: 10.7498/aps.70.20210281 http://dx.doi.org/10.7498/aps.70.20210281 光纤激光基模光束的β因子 ZHANG Y Q ， HUANG L J ， CHANG Q ， et al . Β factor of fundamental mode of fiber laser beam ［J］. Acta Physica Sinica ， 2021 ， 70 （ 20 ）： 79 - 85 . （in Chinese） . doi: 10.7498/aps.70.20210281 http://dx.doi.org/10.7498/aps.70.20210281 Β factor of fundamental mode of fiber laser beam 10.37188/OPE.20223012.1418 TN206 TN253 A 1004-924X（2022）12-1418-11 罗杰,秦来安,侯再红等.应用于光束质量测量的阵列光纤串扰校正[J].光学精密工程,2022,30(12):1418-1428. 罗杰,秦来安,侯再红等.应用于光束质量测量的阵列光纤串扰校正[J].光学精密工程,2022,30(12):1418-1428. LUO Jie,QIN Laian,HOU Zaihong,et al.Crosstalk correction of array fibers applied to beam quality measurement[J].Optics and Precision Engineering,2022,30(12):1418-1428. LUO Jie,QIN Laian,HOU Zaihong,et al.Crosstalk correction of array fibers applied to beam quality measurement[J].Optics and Precision Engineering,2022,30(12):1418-1428. 国家自然科学基金项目（No.41875033） 2022-01-12 2022-02-21 2022-06-25 2022-06-25 光学精密工程 12 30