您当前的位置：
首页 >
文章列表页 >
基于CNN-SPGD算法的相干光通信像差校正方法研究
信息科学 | 更新时间：2022-04-26
|
基于CNN-SPGD算法的相干光通信像差校正方法研究
- Research on aberration correction method of coherent optical communication based on CNN-SPGD algorithm
- 光学精密工程
- 作者机构：
1.吉林大学 通信工程学院，吉林 长春 130012
2.中国科学院 长春光学精密机械与物理研究所，吉林 长春 130033
3.陕西科技大学 陕西人工智能联合实验室， 陕西 西安 710021
- 作者简介：
[ "刘 维（1986-），女，吉林延边人，副教授，硕士生导师，2009年、2015年于吉林大学分别获得学士、博士学位，2017年博士后流动站出站，主要从事基于自适应光学的相干探测地形测绘激光雷达大气补偿关键技术、自由空间光通信大气信道补偿技术、基于深度学习的目标识别技术等方面的研究。E-mail： jdlw@jlu.edu.cn" ]
刘 维（1986-），女，吉林延边人，副教授，硕士生导师，2009年、2015年于吉林大学分别获得学士、博士学位，2017年博士后流动站出站，主要从事基于自适应光学的相干探测地形测绘激光雷达大气补偿关键技术、自由空间光通信大气信道补偿技术、基于深度学习的目标识别技术等方面的研究。E-mail： jdlw@jlu.edu.cn
[ "石文孝（1960-），男，黑龙江哈尔滨人，教授，博士生导师，1983年于长春邮电学院获得学士学位，1991年于哈尔滨工业大学获得硕士学位，2006年于吉林大学获得博士学位，主要从事宽带无线移动通信网理论与技术、高速移动互联网理论与技术、无线光通信理论与技术等方面的研究。E-mail： swx@jlu.edu.cn" ]
石文孝（1960-），男，黑龙江哈尔滨人，教授，博士生导师，1983年于长春邮电学院获得学士学位，1991年于哈尔滨工业大学获得硕士学位，2006年于吉林大学获得博士学位，主要从事宽带无线移动通信网理论与技术、高速移动互联网理论与技术、无线光通信理论与技术等方面的研究。E-mail： swx@jlu.edu.cn
- 基金信息：国家自然科学基金资助项目(62001448)
- DOI：10.37188/OPE.20223006.0743 中图分类号： O436
纸质出版日期：2022-03-25，
收稿日期：2021-08-26，
修回日期：2021-10-11，
扫 描 看 全 文
引用本文
刘维,徐珺楠,金玳冉等.基于CNN-SPGD算法的相干光通信像差校正方法研究[J].光学精密工程,2022,30(06):743-754.
LIU Wei,XU Junnan,JIN Dairan,et al.Research on aberration correction method of coherent optical communication based on CNN-SPGD algorithm[J].Optics and Precision Engineering,2022,30(06):743-754.
刘维,徐珺楠,金玳冉等.基于CNN-SPGD算法的相干光通信像差校正方法研究[J].光学精密工程,2022,30(06):743-754. DOI： 10.37188/OPE.20223006.0743.
LIU Wei,XU Junnan,JIN Dairan,et al.Research on aberration correction method of coherent optical communication based on CNN-SPGD algorithm[J].Optics and Precision Engineering,2022,30(06):743-754. DOI： 10.37188/OPE.20223006.0743.
基于CNN-SPGD算法的相干光通信像差校正方法研究
基于CNN-SPGD算法的相干光通信像差校正方法研究 Research on aberration correction method of coherent optical communication based on CNN-SPGD algorithm 1 first-author 刘 维（1986-），女，吉林延边人，副教授，硕士生导师，2009年、2015年于吉林大学分别获得学士、博士学位，2017年博士后流动站出站，主要从事基于自适应光学的相干探测地形测绘激光雷达大气补偿关键技术、自由空间光通信大气信道补偿技术、基于深度学习的目标识别技术等方面的研究。E-mail： jdlw@jlu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021160&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021162&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021161&type= 22.01332855 32.00399780 刘 维 （1986-），女，吉林延边人，副教授，硕士生导师，2009年、2015年于吉林大学分别获得学士、博士学位，2017年博士后流动站出站，主要从事基于自适应光学的相干探测地形测绘激光雷达大气补偿关键技术、自由空间光通信大气信道补偿技术、基于深度学习的目标识别技术等方面的研究。E-mail： jdlw@jlu.edu.cn jdlw@jlu.edu.cn 1 1 1 corresp E-mail： swx@jlu.edu.cn E-mail： swx@jlu.edu.cn 石文孝（1960-），男，黑龙江哈尔滨人，教授，博士生导师，1983年于长春邮电学院获得学士学位，1991年于哈尔滨工业大学获得硕士学位，2006年于吉林大学获得博士学位，主要从事宽带无线移动通信网理论与技术、高速移动互联网理论与技术、无线光通信理论与技术等方面的研究。E-mail： swx@jlu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021163&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021165&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021164&type= 22.01332855 32.00400543 石文孝 （1960-），男，黑龙江哈尔滨人，教授，博士生导师，1983年于长春邮电学院获得学士学位，1991年于哈尔滨工业大学获得硕士学位，2006年于吉林大学获得博士学位，主要从事宽带无线移动通信网理论与技术、高速移动互联网理论与技术、无线光通信理论与技术等方面的研究。E-mail： swx@jlu.edu.cn swx@jlu.edu.cn 1 2 3 吉林大学 通信工程学院，吉林 长春 130012 吉林大学 通信工程学院，吉林 长春 130012 College of Communication Engineering， Jilin University， Changchun 130012， China College of Communication Engineering， Jilin University， Changchun 130012， China 中国科学院 长春光学精密机械与物理研究所，吉林 长春 130033 中国科学院 长春光学精密机械与物理研究所，吉林 长春 130033 Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China 陕西科技大学 陕西人工智能联合实验室， 陕西 西安 710021 陕西科技大学 陕西人工智能联合实验室， 陕西 西安 710021 Shaanxi Joint Laboratory of Artificial Intelligence，Shaanxi University of Science & Technology， Xi'an 710021， China Shaanxi Joint Laboratory of Artificial Intelligence，Shaanxi University of Science & Technology， Xi'an 710021， China 为了降低大气湍流对自由空间光通信系统传输性能的影响，建立了一套自适应光学校正系统来校正高阶像差。对该系统的混频效率，误码率和迭代次数等参数进行研究。首先，将图像的Zernike系数划分为8个大类，256个小类，通过CNN模型预测输入光斑的大类。对于每一小类的像差，取每一阶系数范围的中点作为标准点形成一个标准Zernike系数向量，按照其对应的标准电压作为初始校正电压进行校正。实验结果表明：随机并行梯度下降算法使混频效率达到0.80、0.85和0.9需要的迭代次数分别为110次、161次和280次；在使混频效率达到相同值的前提下，CNN-SPGD（Convolutional Neural Networks-Stochastic Parallelism Gradient Descent）算法所需的迭代次数分别为4次、37次和141次。而在相同的迭代次数下，CNN-SPGD算法在系统中的混频效率更高，误码率更低。CNN-SPGD算法与SPGD算法相比，像差校正速度更快。该CNN-SPGD算法可以大幅度减少传统自适应光学系统的迭代次数，满足激光通信的各种需求。 In order to reduce the influence of atmospheric turbulence on the transmission performance of free-space optical communication systems， an adaptive optical correction system is established to correct high-order aberrations. The mixing efficiency， bit error rate， and iteration times of the system are studied. First， the Zernike coefficient of the image is divided into 8 categories and 256 sub-categories， and the categories of input spots are predicted by a CNN model. For each small class of aberrations， the midpoint of the coefficient range of each order is taken as the standard point to form a standard Zernike coefficient vector， and its corresponding standard voltage is taken as the initial correction voltage. The experimental results show that the number of iterations required by the SPGD algorithm to achieve mixing efficiencies of 0.80， 0.85， and 0.9 is 110， 161， and 280， respectively. To reach the same mixing efficiencies， the number of iterations required by the CNN-SPGD algorithm is 4， 37， and 141. For the same number of iterations， the mixing efficiency of the CNN-SPGD algorithm in the system is higher and the bit error rate is lower. Moreover， the CNN-SPGD algorithm has a faster aberration correction speed compared to the SPGD algorithm. The CNN-SPGD algorithm can considerably reduce the number of iterations required by traditional adaptive optics systems and meet the various needs of laser communication. 自由空间光通信 无波前传感 混频效率 误码率 free space optical communication wavefront sensorless mixing efficiency bit error rate 1 引 言 1 自由空间光通信（Free Space Optical Communication，FSOC）是一种无线通信技术，具有抗电磁干扰能力强、无需许可证费用以及易于部署等优势 ［ 1 1 ］ 。根据通信体制的不同可以将FSOC技术划分为两种类型：非相干光通信和相干光通信。非相干光通信是采用强度调制/直接检测的方法实现通信的一种技术。这种通信方式在技术上较为成熟且成本较低。与非相干光通信系统相比，相干光通信系统在接收端采用的是相干检测的方式，主要具有调制/解调方式多样化、接收灵敏度更高以及无中继传输距离更远等优势，是FSOC技术未来的重要发展方向之一。 当相干FSOC系统在大气环境中通信时，通信性能会受到大气湍流的影响。自适应光学（Adaptive Optics，AO）技术已经被证明可以有效降低大气湍流对FSOC系统传输性能的影响 ［ 2 2 ］ 。刘超等应用127单元的AO系统对零差检测FSOC系统进行了波前畸变补偿实验，实验结果显示AO技术可以显著降低大气湍流对系统传输性能的影响 ［ 3 3 ］ 。李明等对海洋大气湍流影响下相干FSOC系统的通信性能进行了研究。实验结果显示，当应用闭环AO系统对波前畸变进行补偿时，相干FSOC系统的BER明显降低 ［ 4 4 ］ 。中国科学院长春光学精密机械与物理研究所的研究人员将基于349单元连续镜面变形镜的大规模高速AO系统应用于相干FSOC系统中，分析了相干FSOC系统的性能变化情况，计算了不同Greenwood频率下的混频效率和误码率并将其作为相干FSOC系统的性能指标。实验结果表明，相干FSOC系统的性能得到了显著提升 ［ 5 5 ］ 。 传统AO系统主要包括波前传感器、波前控制器和波前校正器 ［ 6 6 ］ 。除此之外，一种新型的AO技术——无波前传感AO技术也受到了越来越多的关注。李兆坤等提出了一种基于反向传播神经网络的波前检测方法，有效提升了AO系统的实时性 ［ 7 7 ］ 。马慧敏等提出了一种基于卷积神经网络（Convolutional Neural Networks，CNN）的无波前传感算法。卷积神经网络的作用是根据输入图像预测其对应的泽尼克系数。结果显示，经过畸变补偿后，光斑的斯特列尔比有了明显提高 ［ 8 8 ］ 。西安理工大学的研究人员开发了AO校正模块和偏振控制模块来缓解大气湍流对长距离相干FSOC系统性能的影响，并建立了基于外差差分二进制相移键控调制的相干FSOC系统，实现了10.2 km的远距离视频传输，其中AO校正模块采用了无波前传感技术 ［ 9 9 ］ 。 随机并行梯度下降算法（Stochastic Parallelism Gradient Descent，SPGD）是使用范围最广的无波前传感算法之一，容易实现且较为成熟，广泛应用于高分辨率成像领域。但其收敛速度较慢，实时性较差，这对相干FSOC系统的应用存在不利影响。自适应增益系数的Admspgd算法是一种提高收敛速度的有效算法 ［ 10 10 ］ ，此算法通过对增益系数的实时调整，改善了光纤激光组合系统的收敛速度、鲁棒性和有效带宽。朱丹等 ［ 11 11 ］ 对弹性反向算法和SPGD算法的像差校正能力比较后，得到了SPGD算法在复杂波前畸变的恢复能力较强的结论，证明了SPGD算法的巨大潜力。 为了提升SPGD算法的像差校正速度，本文将卷积神经网络与SPGD算法相结合，提出了CNN-SPGD算法。CNN用来对输入光斑进行分类，得到光斑的类型，进而得到对应的像差类型，之后由SPGD算法完成波前像差校正。实验的结果表明，在达到相同系统性能的前提下，CNN-SPGD算法相比于SPGD算法像差校正速度更快。 2 相干FSOC技术原理 1 无波前相干FSOC的系统框图如 图1 图1 所示 ［ 9 9 ］ ，经过大气湍流扰动的光波通过变形镜的闭环控制得以校正，其中CCD被用来读取光斑，PC通过算法得到相应的控制电压来控制变形镜进行形变，从而将畸变波前校正至正常波前，然后对信号光进行相干检测接收。 在接收端，信号光与本地振荡器生成的本振光进行混频。本振光偏振态较为稳定，而信号光经过光学系统和大气后偏振态可能会发生一些变化。在相干FSOC系统中可以采用偏振控制器保持信号光稳定的偏振态输出。信号光和本振光可以表示为如下形式 ［ 12 12 ］ ： <math display="block" id="M1"><msub><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mi>t</mi><mo>+</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021170&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021169&type= 41.14799881 4.91066647 ， <math display="block" id="M2"><msub><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mi>t</mi><mo>+</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021172&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021171&type= 47.58266449 4.91066647 ， 其中： <math id="M3"><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021174&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021173&type= 4.23333359 3.72533321 、 <math id="M4"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021176&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021175&type= 3.89466691 3.72533321 和 <math id="M5"><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021178&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021177&type= 3.72533321 4.65666676 分别表示信号光的振幅、角频率和初始相位； <math id="M6"><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021420&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021408&type= 5.84200001 3.72533321 ， <math id="M7"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021182&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021181&type= 5.50333309 3.72533321 和 <math id="M8"><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021400&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021399&type= 5.33400011 4.65666676 分别表示本振光的振幅、角频率和初始相位。 信号光和本振光在光电检测器上形成的光场可以用 <math id="M9"><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>+</mo><msub><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021138&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021185&type= 16.59466553 4.82600021 表示。根据检测器的平方定律可以得到光电检测器响应输出的电流为： <math display="block" id="M10"><mtable><mtr><mtd><mi>I</mi><mo>=</mo><mi>K</mi><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>+</mo><msub><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>=</mo><mi>K</mi><mo stretchy="false">{</mo><msubsup><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><msup><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mi>t</mi><mo>+</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><msubsup><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><msup><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mi>t</mi><mo>+</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo><mo>+</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">[</mo><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>+</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo><mi>t</mi><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>+</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>+</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">[</mo><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo><mi>t</mi><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo stretchy="false">}</mo><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021140&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021139&type= 108.79666138 18.11866570 其中， K 为比例系数。 式（3） 式（3） 中，电流 <math id="M11"><mi>I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021142&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021141&type= 1.26999998 2.79399991 可以看作由四个单项式的和组成：第一个单项式和第二个单项式表示直流项。第三个单项式中， <math id="M12"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>+</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021189&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021187&type= 14.13933277 3.80999994 的频率较高，一般的检测器无法响应。令 <math id="M13"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">F</mi></mrow></msub><mo>=</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021192&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021191&type= 23.62199974 3.72533321 ，当 <math id="M14"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021202&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021146&type= 4.65666676 3.72533321 低于检测器的截止频率 <math id="M15"><msub><mrow><mi>ω</mi></mrow><mrow><mi>C</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021198&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021196&type= 4.06400013 3.72533321 时，检测器的响应电流为 ［ 13 13 ］ ： <math display="block" id="M16"><mi>I</mi><mo>=</mo><mi>K</mi><mfenced open="{" close="}" separators="|"><mrow><mfrac><mrow><msubsup><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><msubsup><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo>+</mo><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">[</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">F</mi></mrow></msub><mi>t</mi><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021200&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021199&type= 82.71932983 9.65200043 . 由 式（4） 式（4） 可以看出，检测器的响应输出电流中包含直流和交流分量，所要传输的信息包含在交流分量中。相干接收机接收到的信号频率一般限制在以 <math id="M17"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021202&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021146&type= 4.65666676 3.72533321 为中心的带通范围内，可直接过滤掉不含信息的直流分量。直流分量被滤除后，可以用下式表示检测器响应电流的中频分量 ［ 12 12 ］ ： <math display="block" id="M18"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">F</mi></mrow></msub><mo>=</mo><mi>K</mi><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">[</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">F</mi></mrow></msub><mi>t</mi><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021205&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021204&type= 59.52066422 4.91066647 . 式（4） 式（4） 还可以用光功率的形式描述（将式中的光强度 <math id="M19"><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021406&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021405&type= 4.23333359 3.72533321 和 <math id="M20"><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021420&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021408&type= 5.84200001 3.72533321 用光功率 <math id="M21"><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021216&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021214&type= 3.80999994 3.72533321 和 <math id="M22"><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021223&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021218&type= 5.41866684 3.72533321 代替） ［ 14 14 ］ ： <math display="block" id="M23"><mtable columnalign="left"><mtr><mtd><mi>I</mi><mo>=</mo><mi>R</mi><mo stretchy="false">(</mo><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>+</mo><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mn mathvariant="normal">2</mn><mi>R</mi><mroot><mrow><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></mrow><mrow/></mroot><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">[</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">F</mi></mrow></msub><mi>t</mi><mo>+</mo></mtd></mtr><mtr><mtd><mo stretchy="false">(</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021221&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021224&type= 67.47933960 11.68400002 其中， <math id="M24"><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021240&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021238&type= 2.70933342 2.79399991 为检测器的响应度。根据 <math id="M25"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021226&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021242&type= 3.89466691 3.72533321 和 <math id="M26"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021228&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021244&type= 5.50333309 3.72533321 的关系，相干检测可以进一步分为零差检测和外差检测两种检测方式 ［ 13 13 ］ 。当 <math id="M27"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">F</mi></mrow></msub><mo>≠</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021247&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021229&type= 11.34533310 3.97933316 ，即 <math id="M28"><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>≠</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021249&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021248&type= 14.13933277 3.97933316 时，这种检测方式称之为外差检测；反之，称为零差检测。 3 SPGD算法 1 在AO系统中，求解波前校正器最优控制电压属于多维变量最优化问题。SPGD算法适用于多个控制变量的优化过程，其结构简单且算法复杂度低。采用SPGD算法进行像差校正的数学模型为 ［ 15 15 ］ ： <math display="block" id="M29"><mtable><mtr><mtd><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>J</mi><mo stretchy="false">[</mo><mi>u</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mtext> </mtext><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><mi>f</mi><mo stretchy="false">[</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>k</mi></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>,</mo><mi>u</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><mi>φ</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021253&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021251&type= 46.73600006 16.08666801 ， 其中： <math id="M30"><mi>J</mi><mo stretchy="false">[</mo><mi>u</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021256&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021232&type= 13.80066681 3.55599999 是表征性能指标 J 与控制变量 u （ r ）之间关系的目标函数，本文使用SR作为性能指标 J 评估SPGD算法对零差检测FSOC系统混频效率的影响； u （ r ）是波前校正器生成的补偿相位，作为控制变量；假设经过 k 次迭代校正后的残余波前像差为 <math id="M31"><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>k</mi></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021234&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021233&type= 8.72066593 4.82600021 ，则 k +1次迭代校正后的波前为： <math display="block" id="M32"><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>k</mi></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>+</mo><mi>u</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021260&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021235&type= 37.25333405 4.82600021 . 当 k =0时， <math id="M33"><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021263&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021237&type= 8.89000034 4.82600021 称为初始像差。 SPGD算法的计算流程可以描述为： 定义性能指标 J 的变化量： <math display="block" id="M34"><mi mathvariant="normal">Δ</mi><msup><mrow><mi>J</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msup><mo>=</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mi>J</mi></mrow><mrow><mo>+</mo></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>-</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mi>J</mi></mrow><mrow><mo>-</mo></mrow><mrow><mi>k</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021267&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021265&type= 29.71800041 3.97933316 ， <math display="block" id="M35"><mi mathvariant="normal">Δ</mi><msubsup><mrow><mi>J</mi></mrow><mrow><mo>±</mo></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><mi>J</mi><mo stretchy="false">[</mo><msup><mrow><mi>V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msup><mo>±</mo><mi mathvariant="normal">Δ</mi><msup><mrow><mi>V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msup><mo stretchy="false">]</mo><mo>-</mo><mi>J</mi><mo stretchy="false">[</mo><msup><mrow><mi>V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msup><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021285&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021283&type= 53.59400177 4.65666676 . 新的控制电压 <math id="M36"><mi>V</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021287&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021286&type= 3.04800010 2.87866688 可以通过下式获得： <math display="block" id="M37"><msup><mrow><mi>V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo>+</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msup><mo>=</mo><msup><mrow><mi>V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msup><mo>+</mo><mi>γ</mi><mi mathvariant="normal">Δ</mi><msup><mrow><mi mathvariant="bold-italic">V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msup><mi mathvariant="normal">Δ</mi><msup><mrow><mi>J</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021291&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021289&type= 42.67200089 3.89466691 ， 其中： <math id="M38"><mi>γ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021295&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021293&type= 1.94733346 3.64066648 是正的增益系数； k 代表迭代次数； <math id="M39"><mi mathvariant="normal">Δ</mi><msup><mrow><mi mathvariant="bold-italic">V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msup><mo>=</mo><mo stretchy="false">[</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>,</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>,</mo><mo>⋯</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mi>V</mi></mrow><mrow><mi>m</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msubsup><msup><mrow><mo stretchy="false">]</mo></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021297&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021296&type= 51.56200027 4.65666676 ，代表第 k 次迭代时生成的扰动电压； m 的值等于驱动器的个数。 SPGD算法的具体过程为： （1）设置增益系数 <math id="M40"><mi>γ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021561&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021549&type= 1.94733346 3.64066648 和控制电压 <math id="M41"><mi>V</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021274&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021273&type= 3.04800010 2.87866688 的初始值 <math id="M42"><msup><mrow><mi>V</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021304&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021303&type= 5.92666674 3.13266683 ； （2）生成随机扰动向量 <math id="M43"><mi mathvariant="normal">Δ</mi><msup><mrow><mi mathvariant="bold-italic">V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msup><mo>=</mo><mo stretchy="false">[</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>,</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>,</mo><mo>⋯</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mi>V</mi></mrow><mrow><mi>m</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo stretchy="false">)</mo></mrow></msubsup><msup><mrow><mo stretchy="false">]</mo></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021276&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021306&type= 51.64666748 4.65666676 ； （3）由 式（11） 式（11） 更新电压值，生成新的电压向量 <math id="M44"><msup><mrow><mi mathvariant="bold-italic">V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo>+</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021309&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021308&type= 9.39799976 3.13266683 ； （4）计算 k +1次迭代校正后性能指标 J 的值。如果迭代次数达到最大值 <math id="M45"><mi>μ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021576&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021590&type= 1.94733346 3.64066648 则停止迭代；否则返回（2）并继续生成扰动电压直至迭代次数达到最大值 <math id="M46"><mi>μ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021576&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021590&type= 1.94733346 3.64066648 。 4 变形镜模型 1 本文的实验中采用的波前校正器是标准分布的32单元连续镜面变形镜（Deformable Dirrors，DMs） ［ 15 15 ］ ，其驱动器分布如 图2 图2 所示。其中黑点表示波前校正器的中点，32个变形镜按照标号的位置均匀分布在波前校正器上，标号1~32。这32个变形镜就可以完美覆盖大圆圈中所示的波前校正器，其实物图如 图3 图3 ，为中国科学院长春光学精密机械与物理研究所研制的32单元连续镜面变形镜。 每一个变形镜由电机驱动，其电压大小可以线性控制变形镜的位置，使波前像差得到补偿。如 图4 图4 所示。 一般可以用高斯函数模型来模拟DMs的光学影响函数。DMs镜面坐标 <math id="M47"><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021339&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021337&type= 9.05933285 3.97933316 处，第 <math id="M48"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021342&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021322&type= 1.10066664 2.79399991 个驱动器生成的补偿相位可以表示为如下形式 ［ 16 16 ］ ： <math display="block" id="M49"><msub><mrow><mi>ψ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced open="[" close="]" separators="|"><mrow><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mi>ω</mi><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mroot><mrow><mo stretchy="false">(</mo><mi>x</mi><mo>-</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mo stretchy="false">(</mo><mi>y</mi><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot></mrow><mrow><mi>d</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mi>α</mi></mrow></msup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021344&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021326&type= 73.91400146 14.56266594 ， 其中： <math id="M50"><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021346&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021359&type= 2.79399991 3.72533321 和 <math id="M51"><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021348&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021363&type= 2.62466669 4.57200003 分别代表第 i 个驱动器在DMs镜面中的横坐标和纵坐标； <math id="M52"><mi>ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021368&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021366&type= 2.53999996 2.79399991 指相邻驱动器之间的交联值； d 代表驱动器的归一化间距； <math id="M53"><mi>α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021356&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021370&type= 1.94733346 2.79399991 是高斯函数指数。DM的某个子孔径的局部斜率就是其波前函数在该范围的局部导数，而电压与变形镜的斜率线性相关，所以通过相位共轭原理就可使变形镜的补偿相位与畸变波前相位相符。 在多个驱动器的共同作用下，变形镜会发生形变，产生的全部补偿相位为 ［ 17 17 ］ ： <math display="block" id="M54"><mi>u</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>j</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">32</mn></mrow></munderover></mstyle><mrow><msub><mrow><mi>V</mi></mrow><mrow><mi>j</mi></mrow></msub><msub><mrow><mi>ψ</mi></mrow><mrow><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021376&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021371&type= 37.42266846 8.97466660 ， 其中， <math id="M55"><msub><mrow><mi>V</mi></mrow><mrow><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021390&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021378&type= 3.72533321 3.80999994 是施加到第 j 个驱动器上的电压。 5 性能指标分析 1 5.1 混频效率 2 混频效率是评估相干FSOC系统性能的一项重要指标。信号光在与本振光混频之前，会受到诸如空间传输、大气湍流等多种因素的影响，这些影响最终会在混频效率这一性能指标上得到体现。因此，混频效率可以被认为是一个综合性的指标。对于采用零差检测的FSOC系统，混频效率定义为 ［ 17 17 ］ ： <math display="block" id="M56"><mi>η</mi><mo>=</mo><mfrac><mrow><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mo>∫</mo><mrow><mi>U</mi></mrow></msub><mrow><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><mi mathvariant="normal">Δ</mi><mi>φ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>U</mi></mrow></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msub><mo>∫</mo><mrow><mi>U</mi></mrow></msub><mrow><msubsup><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mi mathvariant="normal">d</mi><mi>U</mi><msub><mo>∫</mo><mrow><mi>U</mi></mrow></msub><mrow><msubsup><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mi mathvariant="normal">d</mi><mi>U</mi></mrow></mrow></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021393&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021381&type= 46.48199844 16.51000023 ， 其中： <math id="M57"><mi mathvariant="normal">Δ</mi><mi>φ</mi><mo>=</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021384&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021395&type= 23.45266724 4.65666676 ， <math id="M58"><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021403&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021398&type= 3.72533321 4.65666676 和 <math id="M59"><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021400&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021399&type= 5.33400011 4.65666676 分别表示信号光和本振光的相位； <math id="M60"><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021406&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021405&type= 4.23333359 3.72533321 和 <math id="M61"><msub><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">L</mi><mi mathvariant="normal">O</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021420&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021408&type= 5.84200001 3.72533321 分别代表信号光和本振光的振幅； <math id="M62"><mi>U</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021412&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021410&type= 3.04800010 2.79399991 代表光电检测器的入射场。 对于采用零差检测的相干FSOC系统，混频效率近似等于远场焦平面图像的斯特列尔比（Strehl Rate，SR），这一结论对于评估相干FSOC系统的性能很有意义，实验中将通过计算SR来分析畸变补偿前后零差检测FSOC系统混频效率的变化情况。SR与波前像差的均方根（Root Mean Square，RMS）之间存在如下近似关系 ［ 18 18 ］ ： <math display="block" id="M63"><mi>S</mi><mi>R</mi><mo>≈</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">(</mo><mo>-</mo><mi>R</mi><mi>M</mi><msup><mrow><mi>S</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=27021414&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021413&type= 30.98800087 4.14866638 ， 其中，RMS的计算公式为： <math display="block" id="M64"><mi>R</mi><mi>M</mi><mi>S</mi><mo>=</mo><mroot><mrow><mfenced open="<" close=">" separators="|"><mrow><msup><mrow><mfenced separators="|"><mrow><mi>φ</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>-</mo><mfenced open="<" close=">" separators="|"><mrow><mi>φ</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced></mrow><mrow/></mroot><mo>=</mo><mroot><mrow><mfenced open="<" close=">" separators="|"><mrow><msup><mrow><mi>φ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><mi>ρ</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo></mrow></mfenced><mo>-</mo><msup><mrow><mfenced open="<" close=">" separators="|"><mrow><mi>φ</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021431&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021415&type= 102.53133392 8.38199997 ， 其中： <math id="M65"><mfenced open="<" close=">" separators="|"><mrow><mo>⋅</mo></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021436&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021434&type= 4.14866638 3.47133350 代表取数学期望， <math id="M66"><mi>φ</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo>,</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021451&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021450&type= 11.59933376 3.97933316 表示入射光束的波前像差， <math id="M67"><mi>ρ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021442&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021439&type= 1.86266661 3.64066648 表示接收器内任意一点到接收器中心的极径， <math id="M68"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021444&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021455&type= 1.77800000 2.79399991 表示该点的极角。 5.2 误码率 2 相干FSOC系统的误码率（Bit error rate，BER）定义为 ［ 19 19 ］ ： <math display="block" id="M69"><mi>B</mi><mi>E</mi><mi>R</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mi mathvariant="normal">e</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">f</mi><mi mathvariant="normal">c</mi><mfenced separators="|"><mrow><mfrac><mrow><mi>Q</mi></mrow><mrow><mroot><mrow><mn mathvariant="normal">2</mn></mrow><mrow/></mroot></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021463&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021446&type= 31.07266808 10.66800022 ， 其中， <math id="M70"><mi mathvariant="normal">e</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">f</mi><mi mathvariant="normal">c</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021467&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021466&type= 11.43000031 3.55599999 为互补误差函数，具体表达式为： <math display="block" id="M71"><mi mathvariant="normal">e</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">f</mi><mi mathvariant="normal">c</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mroot><mrow><mi>π</mi></mrow><mrow/></mroot></mrow></mfrac><msubsup><mo>∫</mo><mrow><mi>x</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mo>-</mo><msup><mrow><mi>ξ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></msup><mi mathvariant="normal">d</mi><mi>ξ</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021470&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021468&type= 36.15266418 9.22866726 . Q 定义为： <math display="block" id="M72"><mi>Q</mi><mo>=</mo><mroot><mrow><mi>S</mi><mi>N</mi><mi>R</mi></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021474&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021483&type= 18.62666702 4.31799984 ， 其中，SNR是相干FSOC系统的信噪比。 对于一个采用BPSK调制的通信系统，接收端的光信号功率为 ［ 20 20 ］ ： <math display="block" id="M73"><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>=</mo><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">P</mi></mrow></msub><mi>h</mi><mi>v</mi><mi>B</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021489&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021485&type= 19.38866615 3.80999994 ， 其中： <math id="M74"><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">P</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021496&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021493&type= 4.31799984 3.80999994 为接收端单比特光子数； h 、 v 和 B 分别为普朗克常数、载波频率和数据比特率。 在没有大气湍流影响时，SNR定义为 ［ 20 20 ］ ： <math display="block" id="M75"><mi>S</mi><mi>N</mi><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn><mi>δ</mi><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></mrow><mrow><mi>h</mi><mi>v</mi><mi>B</mi></mrow></mfrac><mo>=</mo><mn mathvariant="normal">2</mn><mi>δ</mi><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">P</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021501&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021491&type= 34.79800034 7.95866632 ， 其中， <math id="M76"><mi>δ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021512&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021502&type= 1.94733346 2.79399991 是光电检测器的量子效率。 对于采用零差检测的通信系统， <math id="M77"><mi>S</mi><mi>N</mi><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021507&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021505&type= 7.61999989 2.87866688 和 <math id="M78"><mi>S</mi><mi>N</mi><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021522&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021519&type= 8.80533314 3.80999994 的关系为 ［ 17 17 ］ ： <math display="block" id="M79"><mi>S</mi><mi>N</mi><mi>R</mi><mo>=</mo><mn mathvariant="normal">2</mn><mi>S</mi><mi>N</mi><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>⋅</mo><mi>η</mi><mo>=</mo><mn mathvariant="normal">4</mn><mi>δ</mi><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">P</mi></mrow></msub><mi>η</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021532&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021528&type= 42.16400146 4.57200003 ， 其中， <math id="M80"><mi>η</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021545&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021535&type= 1.86266661 3.64066648 代表混频效率。 将式（ 17 17 ）、（ 19 19 ）和（22）联立可以得到： <math display="block" id="M81"><mi>B</mi><mi>E</mi><mi>R</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mi mathvariant="normal">e</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">f</mi><mi mathvariant="normal">c</mi><mo stretchy="false">(</mo><mroot><mrow><mn mathvariant="normal">2</mn><mi>δ</mi><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">P</mi></mrow></msub><mi>η</mi></mrow><mrow/></mroot><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021547&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021555&type= 39.20066452 7.87400007 . 进一步，根据混频效率和SR的近似关系，可以通过 式（23） 式（23） 得到BPSK调制的零差检测相干FSOC系统的误码率。 6 CNN-SPGD算法 1 卷积神经网络具备利用局部操作实现分级抽象表示的能力，因此非常适合计算机视觉处理。在一个典型的CNN模型中，一般包括卷积层、池化层和全连接层。CNN使用局部连接的方式解决了全连接的权重矩阵参数爆炸问题。此外，经过计算的特征图依然保持着原来的空间位置，同时通过权值共享的方式，使得网络的权重系数进一步减小。 本文提出的CNN-SPGD算法是将CNN模型与SPGD算法组合后进行波前像差校正。在算法中，像差被分为8大类，每一类都有唯一对应的初始电压，CNN负责选出当前像差的类别，在类别被确定的同时，控制器就会得到与之相应的初始电压，这里的初始电压只是大类电压的数学期望，而真正的最优电压还需要通过SPGD算法进一步获得。 CNN-SPGD算法的具体实现过程为： （1）设置SPGD算法中增益系数 <math id="M82"><mi>γ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021561&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021549&type= 1.94733346 3.64066648 的初始值； （2）CNN模型对接收到的光斑进行分类； （3）根据光斑的分类结果得到对应的像差类型，之后对电压进行选择，选择出最优电压作为SPGD算法的初始控制电压 <math id="M83"><msup><mrow><mi>V</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021565&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021563&type= 5.92666674 3.13266683 ； （4）生成随机扰动向量； （5）由 式（11） 式（11） 更新电压值，生成新的电压向量 <math id="M84"><msup><mrow><mi mathvariant="bold-italic">V</mi></mrow><mrow><mo stretchy="false">(</mo><mi>k</mi><mo>+</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021585&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021583&type= 9.39799976 3.13266683 ； （6）计算 k +1次迭代校正后性能指标 J 的值。如果迭代次数达到最大值 <math id="M85"><mi>μ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021576&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021590&type= 1.94733346 3.64066648 则停止迭代；否则返回（4）继续生成扰动电压直至迭代次数达到最大值 <math id="M86"><mi>μ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021576&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021590&type= 1.94733346 3.64066648 。 CNN-SPGD系统控制回路图和系统工作框图分别如 图5 图5 和 图6 图6 所示。 7 实验与性能分析 1 7.1 生成训练样本与电压选择 2 我们可以使用一系列Zernike多项式的线性组合表示波前像差 ［ 21 21 ］ 。在Zernike多项式中，第一项是活塞模式（常数），对成像质量不产生影响；第二项和第三项是倾斜像差，可以采用专用的倾斜校正系统进行校正，所以在本文的实验中不考虑前三阶像差，用第四阶到第十一阶Zernike多项式来表示初始像差，即： <math display="block" id="M87"><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><mi>φ</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">4</mn></mrow><mrow><mn mathvariant="normal">11</mn></mrow></munderover></mstyle><mrow><msub><mrow><mi>a</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>Z</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021626&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021624&type= 41.90999985 8.55133343 . 本文采用Nicolas Roddier提出的方法来获取初始Zernike系数 ［ 22 22 ］ ：首先构造独立的Karhunen-Loeve系数，之后转换为Zernike系数。参数 <math id="M88"><mi>D</mi><mo>/</mo><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021686&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021701&type= 7.28133297 4.14866638 设置为4，其中 <math id="M89"><mi>D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021636&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021634&type= 2.79399991 2.79399991 代表接收孔径， <math id="M90"><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021639&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021647&type= 2.62466669 3.72533321 代表大气相干长度。引入的初始像差Zernike系数如 表1 表1 所示。 Order Coefficient 4th -0.76 5th 0.85 6th 1.44 7th 0.14 8th -0.19 9th -0.11 10th 0.09 11th 0.15 相关的参数设置如下：驱动器的交联值 <math id="M91"><mi>ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021656&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021662&type= 2.53999996 2.79399991 设为0.2，归一化间距 d 为0.392，高斯函数指数 <math id="M92"><mi>α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021670&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021667&type= 1.94733346 2.79399991 为2；加载到变形镜驱动器上的初始控制电压 <math id="M93"><msup><mrow><mi mathvariant="bold-italic">V</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021674&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021688&type= 6.01133299 3.13266683 设置为 <math id="M94"><msup><mrow><mo stretchy="false">[</mo><mn mathvariant="normal">0,0</mn><mo>⋯</mo><mo>,</mo><mn mathvariant="normal">0</mn><mo stretchy="false">]</mo></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021678&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021676&type= 18.37266731 3.80999994 ；增益系数 <math id="M95"><mi>γ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021683&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021695&type= 1.94733346 3.64066648 的值设为1.2；随机扰动向量为32维的向量，每个元素之间相互独立，服从伯努利分布，实验中设置为 <math id="M96"><mfenced open="|" close="|" separators="|"><mrow><mi mathvariant="normal">Δ</mi><msub><mrow><mi>V</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfenced><mo>=</mo><mn mathvariant="normal">0.3</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021700&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021698&type= 17.94933319 4.91066647 。通过 式（15） 式（15） ，可以得到初始SR为0.33。 为确定各阶系数的范围，生成了多组Zernike系数（参数 <math id="M97"><mi>D</mi><mo>/</mo><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021686&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021701&type= 7.28133297 4.14866638 同样设置为4）。对各阶系数 <math id="M98"><msub><mrow><mi>a</mi></mrow><mrow><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021705&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021703&type= 2.62466669 3.72533321 进行了区间划分，因此，可以将输入像差划分为4×4×4×2×2=256种类型。 各阶Zernike多项式的系数范围呈对称分布。考虑到这一特点，对划分好的各阶Zernike多项式的系数区间作进一步处理，将对称区间归为同一类，如 表2 表2 所示。 Coefficient Range <math id="M99"><msub><mrow><mi>a</mi></mrow><mrow><mn>4</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021707&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021718&type= 2.62466669 3.21733332 <math id="M100"><mo stretchy="false">[</mo><mo>-</mo><mn>1.8</mn><mo>,</mo><mo>-</mo><mn>0.9</mn><mo stretchy="false">)</mo><mo>⋃</mo><mo stretchy="false">[</mo><mn>0.9,1.8</mn><mo stretchy="false">]</mo><mo>,</mo><mo stretchy="false">[</mo><mo>-</mo><mn>0.9,0</mn><mo stretchy="false">)</mo><mo>⋃</mo><mo stretchy="false">[</mo><mn>0,0.9</mn><mo stretchy="false">)</mo><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021728&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021733&type= 63.92333221 3.04800010 <math id="M101"><msub><mrow><mi>a</mi></mrow><mrow><mn>5</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021719&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021712&type= 2.62466669 3.21733332 <math id="M102"><mo stretchy="false">[</mo><mo>-</mo><mn>1.8</mn><mo>,</mo><mo>-</mo><mn>0.9</mn><mo stretchy="false">)</mo><mo>⋃</mo><mo stretchy="false">[</mo><mn>0.9,1.8</mn><mo stretchy="false">]</mo><mo>,</mo><mo stretchy="false">[</mo><mo>-</mo><mn>0.9,0</mn><mo stretchy="false">)</mo><mo>⋃</mo><mo stretchy="false">[</mo><mn>0,0.9</mn><mo stretchy="false">)</mo><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021728&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021733&type= 63.92333221 3.04800010 <math id="M103"><msub><mrow><mi>a</mi></mrow><mrow><mn>6</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021717&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021723&type= 2.62466669 3.21733332 <math id="M104"><mo stretchy="false">[</mo><mo>-</mo><mn>1.8</mn><mo>,</mo><mo>-</mo><mn>0.9</mn><mo stretchy="false">)</mo><mo>⋃</mo><mo stretchy="false">[</mo><mn>0.9,1.8</mn><mo stretchy="false">]</mo><mo>,</mo><mo stretchy="false">[</mo><mo>-</mo><mn>0.9,0</mn><mo stretchy="false">)</mo><mo>⋃</mo><mo stretchy="false">[</mo><mn>0,0.9</mn><mo stretchy="false">)</mo><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021728&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021733&type= 63.92333221 3.04800010 <math id="M105"><msub><mrow><mi>a</mi></mrow><mrow><mn>7</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021739&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021729&type= 2.62466669 3.21733332 <math id="M106"><mo stretchy="false">[</mo><mo>-</mo><mn>0.9,0</mn><mo stretchy="false">)</mo><mo>⋃</mo><mo stretchy="false">[</mo><mn>0,0.9</mn><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021744&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021753&type= 24.97666740 3.04800010 <math id="M107"><msub><mrow><mi>a</mi></mrow><mrow><mn>8</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021751&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021749&type= 2.62466669 3.21733332 <math id="M108"><mo stretchy="false">[</mo><mo>-</mo><mn>0.9,0</mn><mo stretchy="false">)</mo><mo>⋃</mo><mo stretchy="false">[</mo><mn>0,0.9</mn><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021744&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021753&type= 24.97666740 3.04800010 <math id="M109"><msub><mrow><mi>a</mi></mrow><mrow><mn>9</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021763&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021746&type= 2.62466669 3.21733332 <math id="M110"><mo stretchy="false">[</mo><mo>-</mo><mn>0.85,0.85</mn><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021778&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021762&type= 18.20333290 3.04800010 <math id="M111"><msub><mrow><mi>a</mi></mrow><mrow><mn>10</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021761&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021759&type= 3.55599999 3.21733332 <math id="M112"><mo stretchy="false">[</mo><mo>-</mo><mn>0.85,0.85</mn><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021778&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021762&type= 18.20333290 3.04800010 <math id="M113"><msub><mrow><mi>a</mi></mrow><mrow><mn>11</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021767&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021782&type= 3.55599999 3.21733332 <math id="M114"><mo stretchy="false">[</mo><mo>-</mo><mn>0.5,0.5</mn><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021783&type= 15.07066536 3.04800010 如 表2 表2 所示，各阶Zernike系数的区间经过重组后，像差分为8大类，其中 a 4 ， a 5 和 a 6 各有两大种，组合为8种，则光斑也被相应的划分为8个大类别；每一个大类别的像差中，又可以通过 表2 表2 细分为2×2×2×2×2=32个小类别，其中 a 4 ， a 5 ， a 6 ， a 7 和 a 8 各有两种，组合为32种。CNN模型的作用是预测出输入光斑属于8大类光斑中的哪一类，根据分类结果就可以得到相应的像差类型。对于每一大类像差，生成对应的10 000张焦平面光斑作为训练集样本，1 000张焦平面光斑作为验证集样本，1 000张焦平面光斑作为测试集样本，所以训练集样本包含80 000张图像，验证集样本包含8 000张图像，测试集样本包含8 000张图像，每张图像的分辨率设置为128×128像素，生成的图像为jpg格式。每张图像的命名规则为：类别_个数_jpg。例如0_156_jpg，指第0类光斑中（类别从0开始计数，0~7共8大类）的第156张图像，jpg表示生成图像的格式为jpg。 图7 图7 展示了部分类别光斑的图像样本。 对于每一小类别的波前像差，取每一阶系数范围的中点作为标准点形成一个“标准Zernike系数向量”，处于该区间的Zernike系数向量都可以归类为预先设定好的“标准Zernike系数向量” ［ 23 23 ］ 。按照这种方法，一共可以获得256组“标准Zernike系数向量”。由“标准Zernike系数向量”构成的波前像差称之为“标准像差”。选择初始电压的具体步骤为： （1）在像差校正过程开始前，预先用SPGD算法校正256个标准像差，使SR≥0.9时的电压，称之为“标准电压”。 （2）用CNN模型对输入光斑进行分类，根据光斑的分类结果得到相应的像差类型，每一大类像差中包含32个“标准像差”，自然也对应32个“标准电压”。 （3）逐一用32个“标准电压”生成补偿相位校正波前像差，选择使SR最大的“标准电压”作为SPGD算法的初始电压。 7.2 深度学习框架与CNN模型构建 2 实验中采用Keras完成CNN模型的搭建。Keras是一个高层神经网络应用程序编程接口。Keras的底层实现选择Tensorflow，CNN模型用Python语言编写。 实验中构建的CNN模型共计有5个卷积层，5个池化层，2个全连接层和一个输出层，其具体结构如 图8 图8 所示。模型中的Dropout是一种针对神经网络模型的正则化方法，主要是为了减轻模型的过拟合问题 ［ 23 23 ］ ，其基本思想为在对输入样本进行训练过程中，随机的忽略部分神经元，这样可以增强模型的泛化性，因为模型不会过于依赖某些局部的特征。 7.3 实验与算法对比 2 CNN模型构建完成后，在GPU上进行了训练。实验硬件环境如 表3 表3 所示。 CPU GPU Intel Core i7-4770 Nvidia Geforce GTX 1660 SUPER 训练网络时参数设置如下：训练轮数epochs=20；每一批次读取样本的数量batch_size=200。故在每一轮训练中，在训练集中运行的步数steps_per_epoch=80000/200=400步。训练效果可以通过调用evaluate_generator函数来评估。训练结果如 图9 图9 所示。经过20轮训练后，训练集的准确率为97.62%，测试集的准确率为91.46%。 CNN模型训练完成后，像差校正过程开始。由 表1 表1 中的Zernike系数生成的波前像差对应的焦平面图像如 图10 图10 所示。 初始波前像差以及CNN-SPGD算法进行畸变补偿后的残余波前像差随迭代次数变化如 图11 图11 所示；相应的焦平面上光斑图像如 图12 图12 所示。 由 图12 图12 可以看出，经过基于CNN-SPGD算法的无波前传感AO系统像差校正后，焦平面的光斑图像分布逐渐集中，说明波前像差已得到校正，光束质量得到改善，在500次迭代时得以收敛，这也验证了CNN-SPGD算法的有效性。为了更直观的对比CNN-SPGD算法相对于SPGD算法在像差校正速度方面的优势，我们将两种无波前传感算法对系统混频效率和误码率的影响曲线进行比较，结果如 图13 图13 所示。 1;2; 由 图13 图13 可见，在混频效率方面，SPGD算法使混频效率达到0.80、0.85和0.9需要的迭代次数分别为110次、161次和280次；在使混频效率达到相同值的前提下，CNN-SPGD算法所需的迭代次数分别为4次、37次和141次。为了对比CNN-SPGD算法和SPGD算法在时间复杂度上的关系，我们在1 660 s的GPU上进行了一百次重复实验，在实验中，光斑分类和电压选择所需的平均时间约为SPGD算法16次迭代的时间。也就是说，CNN-SPGD算法只需大约16次迭代的过程，这可以节省数十次甚至上百次的迭代运算，故CNN-SPGD算法相比于SPGD算法在校正速度上得到了提升。用两种算法对八种不同的初始像差进行像差校正，系统混频效率大于等于0.90时所需的迭代次数如 表4 表4 所示。在误码率方面，CNN-SPGD算法经过152次迭代后，便可以使BER达到 <math id="M115"><msup><mrow><mn mathvariant="normal">10</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">6</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021893&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021901&type= 6.60400009 3.04800010 ；而达到同样的误码率，SPGD算法需要334次迭代。即便考虑到光斑分类和电压选择的过程，CNN-SPGD算法的校正速度依然更快。 2 Experiment 2 Number of iterations SPGD CNN-SPGD 1 233 23 2 155 36 3 224 26 4 128 43 5 215 26 6 248 30 7 148 36 8 194 23 由于CNN模块分类后，每类初始控制电压对应的Zernike系数是该类像差的Zernike系数的期望，所以在该类像差的区间内，无论湍流强度如何，CNN模块都能在一次迭代过程中将像差补偿到很小的程度。在达到相同系统性能的前提下，考虑到光斑分类和电压选择的过程，CNN-SPGD算法与SPGD算法相比仍减少了迭代次数，像差校正速度得到了提升。 为了进一步验证CNN-SPGD算法在不同湍流强度的效果，以大气折射率结构常数作为强弱指标，根据Davis不等式生成了强度不同的三组湍流。其湍流强度和BER常数如 表5 表5 所示。 Turbulence intensity Initial BER BER after correction <math id="M116"><mn>4.17</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>17</mn></mrow></msup><msup><mrow><mi>m</mi></mrow><mrow><mo>-</mo><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021834&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021861&type= 23.53733444 2.62466669 <math id="M117"><mn>9.5</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>4</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021865&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021863&type= 13.97000027 2.62466669 <math id="M118"><mn>4.5</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>6</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021868&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021867&type= 13.97000027 2.62466669 <math id="M119"><mn>6.63</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>14</mn></mrow></msup><msup><mrow><mi>m</mi></mrow><mrow><mo>-</mo><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021836&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021870&type= 23.53733444 2.62466669 <math id="M120"><mn>5.0</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>3</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021883&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021837&type= 13.97000027 2.62466669 <math id="M121"><mn>6.8</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>6</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021885&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021884&type= 13.97000027 2.62466669 <math id="M122"><mn>3.24</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>13</mn></mrow></msup><msup><mrow><mi>m</mi></mrow><mrow><mo>-</mo><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021877&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021886&type= 23.53733444 2.62466669 <math id="M123"><mn>7.4</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>3</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021879&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021878&type= 13.97000027 2.62466669 <math id="M124"><mn>6.9</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>6</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021882&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021880&type= 13.97000027 2.62466669 分析 表5 表5 可知，未校正时的BER与湍流强度有关，湍流强度越强，初始BER越大，而经过CNN-SPGD算法的校正后，BER可以下降三个数量级左右，使系统性能得到提升，证明了算法的有效性。 8 结 论 1 本文针对SPGD算法对大气湍流造成的像差影响补偿的迭代次数较多，畸变补偿速度较慢的缺点，改进了传统自适应光学的校正机制，提出了CNN-SPGD优化方法。该方法通过训练CNN的模型得到像差的分类结果，再把CNN模型和SPGD算法相结合。通过CNN模型对输入光斑进行分类达到减少迭代次数的效果，再使用SPGD校正残余高阶像差，从而提高系统性能，改善校正结果。在自由空间光通信系统中，CNN-SPGD模型使误码率下降到 <math id="M125"><msup><mrow><mn mathvariant="normal">10</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">6</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021893&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27021901&type= 6.60400009 3.04800010 时的迭代次数相对传统SPGD可以降低约80%。而在混频效率达到0.9时，CNN-SPGD模型可以降低约85%的迭代次数。综上可以得出结论，在达到相同系统性能的前提下，即便CNN-SPGD算法有光斑分类和电压选择的过程，相比较于SPGD算法仍减少了迭代次数，有效提升了像差校正的速度，对FSOC的设计有一定的意义。在未来的工作中，我们将优化本系统的训练模型，同时更换更好的GPU来提高训练准确率，降低算法收敛时间。 参考文献 HONG Y Q ， HAN S K . Polarization-dependent SOA-based PolSK modulation for turbulence-robust FSO communication ［J］. Optics Express ， 2021 ， 29 （ 10 ）： 15587 - 15594 . doi: 10.1364/oe.421808 http://dx.doi.org/10.1364/oe.421808 Polarization-dependent SOA-based PolSK modulation for turbulence-robust FSO communication OSBORN J ， TOWNSON M J ， FARLEY O J D ， et al . Adaptive optics pre-compensated laser uplink to LEO and GEO ［J］. Optics Express ， 2021 ， 29 （ 4 ）： 6113 - 6132 . doi: 10.1364/oe.413013 http://dx.doi.org/10.1364/oe.413013 Adaptive optics pre-compensated laser uplink to LEO and GEO LIU C ， CHEN S ， LI X Y ， et al . Performance evaluation of adaptive optics for atmospheric coherent laser communications ［J］. Optics Express ， 2014 ， 22 （ 13 ）： 15554 - 15563 . doi: 10.1364/oe.22.015554 http://dx.doi.org/10.1364/oe.22.015554 Performance evaluation of adaptive optics for atmospheric coherent laser communications LI M ， CVIJETIE M . Coherent free space optics communications over the maritine atmosphere with use of adptive optics for beam wavefront correction ［J］. Applied Optics ， 2015 ， 54 （ 6 ）： 1453 - 1463 . doi: 10.1364/ao.54.001453 http://dx.doi.org/10.1364/ao.54.001453 Coherent free space optics communications over the maritine atmosphere with use of adptive optics for beam wavefront correction YANG L ， YAO K ， WANG J ， et al . Performance analysis of 349-element adaptive optics unit for a coherent free space optical communication system ［J］. Scientific reports ， 2019 ， 9 ： 13150 . doi: 10.1038/s41598-019-48338-3 http://dx.doi.org/10.1038/s41598-019-48338-3 Performance analysis of 349-element adaptive optics unit for a coherent free space optical communication system 高世杰 ， 盛磊 ， 吴志勇 ， 等 . 大气激光通信光斑图像的快速复原与实时检测 ［J］. 光学 精密工程 ， 2015 ， 23 （ 8 ）： 2393 - 2399 . doi: 10.3788/OPE.20152308.2393 http://dx.doi.org/10.3788/OPE.20152308.2393 大气激光通信光斑图像的快速复原与实时检测 GAO S J ， SHENG L ， WU Z Y ， et al . Rapid restoration and real-time detection on spot image of atmospheric laser communication ［J］. Opt. Precision Eng. ， 2015 ， 23 （ 8 ）： 2393 - 2399 . （in Chinese） . doi: 10.3788/OPE.20152308.2393 http://dx.doi.org/10.3788/OPE.20152308.2393 Rapid restoration and real-time detection on spot image of atmospheric laser communication LI Z K ， ZHAO X H . BP artificial neural network based wave front correction for sensor-less free space optics communication ［J］. Optics Communications ， 2017 ， 385 ： 219 - 228 . doi: 10.1016/j.optcom.2016.10.037 http://dx.doi.org/10.1016/j.optcom.2016.10.037 BP artificial neural network based wave front correction for sensor-less free space optics communication MA H M ， LIU H Q ， QIAO Y ， et al . Numerical study of adaptive optics compensation based on Convolutional Neural Networks ［J］. Optics Communications ， 2019 ， 433 ： 283 - 289 . doi: 10.1016/j.optcom.2018.10.036 http://dx.doi.org/10.1016/j.optcom.2018.10.036 Numerical study of adaptive optics compensation based on Convolutional Neural Networks 洪艳 . 基于相干检测的自由空间光通信系统研究 ［D］. 北京 ： 北京邮电大学 ， 2016 . 基于相干检测的自由空间光通信系统研究 HONG Y . Research on Free Space Optical Communication System Based on Coherent Detection ［D］. Beijing ： Beijing University of Posts and Telecommunications ， 2016 . （in Chinese） Research on Free Space Optical Communication System Based on Coherent Detection CHE D B ， LI Y Y ， WU Y H ， et al . Theory of AdmSPGD algorithm in fiber laser coherent synthesis ［J］. Optics Communications ， 2021 ， 492 ： 126953 . doi: 10.1016/j.optcom.2021.126953 http://dx.doi.org/10.1016/j.optcom.2021.126953 Theory of AdmSPGD algorithm in fiber laser coherent synthesis ZHU D ， WANG R Y ， ŽURAUSKAS M ， et al . Automated fast computational adaptive optics for optical coherence tomography based on a stochastic parallel gradient descent algorithm ［J］. Optics Express ， 2020 ， 28 （ 16 ）： 23306 . doi: 10.1364/oe.395523 http://dx.doi.org/10.1364/oe.395523 Automated fast computational adaptive optics for optical coherence tomography based on a stochastic parallel gradient descent algorithm 柯熙政 ， 邓莉君 . 无线光通信 ［M］. 北京 ： 科学出版社 ， 2016 . doi: 10.1007/978-3-658-10924-0_13 http://dx.doi.org/10.1007/978-3-658-10924-0_13 无线光通信 KE X Z ， DENG L J . Wireless Optical Communication ［M］. Beijing ： Science Press ， 2016 . （in Chinese） . doi: 10.1007/978-3-658-10924-0_13 http://dx.doi.org/10.1007/978-3-658-10924-0_13 Wireless Optical Communication 曹景太 . 基于自适应光学的相干自由空间光通信系统性能分析 ［D］. 长春 ： 吉林大学 ， 2017 . doi: 10.1364/oe.25.015299 http://dx.doi.org/10.1364/oe.25.015299 基于自适应光学的相干自由空间光通信系统性能分析 CAO J T . Performance Analysis of Coherent Free Space Optical Communication System Based on Adaptive Optics ［D］. Changchun ： Jilin University ， 2017 . （in Chinese） . doi: 10.1364/oe.25.015299 http://dx.doi.org/10.1364/oe.25.015299 Performance Analysis of Coherent Free Space Optical Communication System Based on Adaptive Optics 曹明华 ， 武鑫 ， 杨顺信 ， 等 . Log-normal湍流信道中超奈奎斯特传输系统的误码性能 ［J］. 光学 精密工程 ， 2020 ， 28 （ 2 ）： 465 - 473 . doi: 10.3788/OPE.20202802.0465 http://dx.doi.org/10.3788/OPE.20202802.0465 Log-normal湍流信道中超奈奎斯特传输系统的误码性能 CAO M H ， WU X ， YANG S X ， et al . BER performance of Faster-than-Nyquist communications under Log-normal turbulence channel ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 2 ）： 465 - 473 . （in Chinese） . doi: 10.3788/OPE.20202802.0465 http://dx.doi.org/10.3788/OPE.20202802.0465 BER performance of Faster-than-Nyquist communications under Log-normal turbulence channel CAO J T ， ZHAO X H ， LI Z K ， et al . Stochastic parallel gradient descent laser beam control algorithm for atmospheric compensation in free space optical communication ［J］. Optik ， 2014 ， 125 （ 20 ）： 6142 - 6147 . doi: 10.1016/j.ijleo.2014.06.127 http://dx.doi.org/10.1016/j.ijleo.2014.06.127 Stochastic parallel gradient descent laser beam control algorithm for atmospheric compensation in free space optical communication 李兆坤 . 自适应光学系统在自由空间光通信中的波前像差校正研究 ［D］. 长春 ： 吉林大学 ， 2017 . doi: 10.1016/j.optcom.2016.10.037 http://dx.doi.org/10.1016/j.optcom.2016.10.037 自适应光学系统在自由空间光通信中的波前像差校正研究 LI Z K . Research of Adaptive Optics for Wave Front Distortion Correction in Free Space Optics Communication ［D］. Changchun ： Jilin University ， 2017 . （in Chinese） . doi: 10.1016/j.optcom.2016.10.037 http://dx.doi.org/10.1016/j.optcom.2016.10.037 Research of Adaptive Optics for Wave Front Distortion Correction in Free Space Optics Communication 王婷婷 . 自由空间光通信中的波前畸变补偿方法研究 ［D］. 长春 ： 吉林大学 ， 2017 . doi: 10.5152/tjg.2021.20791 http://dx.doi.org/10.5152/tjg.2021.20791 自由空间光通信中的波前畸变补偿方法研究 WANG T T . Study on Compensation for Distorted Wavefront in Free Space Optical Communication ［D］. Changchun ： Jilin University ， 2017 . （in Chinese） . doi: 10.5152/tjg.2021.20791 http://dx.doi.org/10.5152/tjg.2021.20791 Study on Compensation for Distorted Wavefront in Free Space Optical Communication LI Z K ， CAO J T ， ZHAO X H ， et al . Combinational-deformable-mirror adaptive optics system for atmospheric compensation in free space communication ［J］. Optics Communications ， 2014 ， 320 ： 162 - 168 . doi: 10.1016/j.optcom.2014.01.042 http://dx.doi.org/10.1016/j.optcom.2014.01.042 Combinational-deformable-mirror adaptive optics system for atmospheric compensation in free space communication CAO J T ， ZHAO X H ， LIU W ， et al . BER analysis of coherent free-space optical communication systems with a focal-plane-based wavefront sensor ［J］. Journal of the Korean Physical Society ， 2018 ， 72 （ 5 ）： 555 - 560 . doi: 10.3938/jkps.72.555 http://dx.doi.org/10.3938/jkps.72.555 BER analysis of coherent free-space optical communication systems with a focal-plane-based wavefront sensor CAO J T ， ZHAO X H ， LIU W ， et al . Performance analysis of a coherent free space optical communication system based on experiment ［J］. Optics Express ， 2017 ， 25 （ 13 ）： 15299 . doi: 10.1364/oe.25.015299 http://dx.doi.org/10.1364/oe.25.015299 Performance analysis of a coherent free space optical communication system based on experiment NOLL R J . Zernike polynomials and atmospheric turbulence ［J］. Journal of the Optical Society of America ， 1976 ， 66 （ 3 ）： 207 . doi: 10.1364/josa.66.000207 http://dx.doi.org/10.1364/josa.66.000207 Zernike polynomials and atmospheric turbulence 王建立 ， 董玉磊 ， 姚凯男 ， 等 . 349单元自适应光学波前处理器 ［J］. 光学 精密工程 ， 2018 ， 26 （ 5 ）： 1007 - 1013 . doi: 10.3788/ope.20182605.1007 http://dx.doi.org/10.3788/ope.20182605.1007 349单元自适应光学波前处理器 WANG J L ， DONG Y L ， YAO K N ， et al . Three hundred and fourty-nine unit adaptive optical wavefront processor ［J］. Opt. Precision Eng. ， 2018 ， 26 （ 5 ）： 1007 - 1013 . （in Chinese） . doi: 10.3788/ope.20182605.1007 http://dx.doi.org/10.3788/ope.20182605.1007 Three hundred and fourty-nine unit adaptive optical wavefront processor GU H J ， LIU M Q ， LIU H Y ， et al . An algorithm combining convolutional neural networks with SPGD for SLAO in FSOC ［J］. Optics Communications ， 2020 ， 475 ： 126243 . doi: 10.1016/j.optcom.2020.126243 http://dx.doi.org/10.1016/j.optcom.2020.126243 An algorithm combining convolutional neural networks with SPGD for SLAO in FSOC 10.37188/OPE.20223006.0743 O436 A 1004-924X（2022）06-0743-12 刘维,徐珺楠,金玳冉等.基于CNN-SPGD算法的相干光通信像差校正方法研究[J].光学精密工程,2022,30(06):743-754. 刘维,徐珺楠,金玳冉等.基于CNN-SPGD算法的相干光通信像差校正方法研究[J].光学精密工程,2022,30(06):743-754. LIU Wei,XU Junnan,JIN Dairan,et al.Research on aberration correction method of coherent optical communication based on CNN-SPGD algorithm[J].Optics and Precision Engineering,2022,30(06):743-754. LIU Wei,XU Junnan,JIN Dairan,et al.Research on aberration correction method of coherent optical communication based on CNN-SPGD algorithm[J].Optics and Precision Engineering,2022,30(06):743-754. 国家自然科学基金青年科学基金项目（No.62001448） 2021-08-26 2021-10-11 2022-03-25 2022-04-26 光学精密工程 6 30
查看原文 
光学 精密工程 | 计算机视觉-光通信的新手段
0
浏览量
47
下载量
1
CSCD
Alert me when the article has been cited
文章被引用时，请邮件提醒。
提交
工具集
下载
- 网络PDF下载
- WORD下载
- Meta-XML下载
参考文献导出
- BibTeX下载
- Endnote下载
- RefMan 下载
- NoteExpress下载
- NoteFirst下载
收藏
添加至我的专辑
关联资源
分享给微信好友或者朋友圈
相关文章
基于单光子探测器的深空激光通信阵列
12.7 km水平信道空间激光通信闪烁抑制
Gamma-Gamma湍流信道下广义空时脉冲位置调制
可见光通信的间隔多脉冲位置调制
显示信息高速串行传输的差错控制
相关作者
刘向南
王建军
刘荣科
李晓亮
佟首峰
江伦
滕云杰
冯悦姝
相关机构
北京航空航天大学 电子信息工程学院
北京遥测技术研究所
长春理工大学 光电工程学院
长春理工大学 空间光电技术国家地方联合工程研究中心
吉林工程技术师范 学院量子信息技术交叉学科研究院
0
回到顶部