基于黑白光栅的宽波段大视场激光告警

Laser warning with broadband and large field of view based on black-white grating 基于黑白光栅的宽波段大视场激光告警 1 2 3 first-author 杨雪梅（1999-），女，甘肃天水人，硕士研究生，2020年于中北大学大学获得学士学位，主要从事光电信息检测方面的研究。E-mail：yxm123_1234@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037390&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037395&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037393&type= 22.01333237 32.00400162 杨雪梅 （1999-），女，甘肃天水人，硕士研究生，2020年于中北大学大学获得学士学位，主要从事光电信息检测方面的研究。E-mail： yxm123_1234@163.com yxm123_1234@163.com 1 2 3 corresp E-mail： zhangrui@nuc.edu.cn E-mail： zhangrui@nuc.edu.cn 张 瑞（1987—），男，山西长治人，博士，教授，硕士生导师，主要从事光电探测、光谱（偏振）成像、光电目标识别、光学高精度位移测量等方面的研究。E-mail：zhangrui@nuc.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037397&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037400&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037399&type= 22.01333618 32.00400162 张 瑞 （1987—），男，山西长治人，博士，教授，硕士生导师，主要从事光电探测、光谱（偏振）成像、光电目标识别、光学高精度位移测量等方面的研究。E-mail： zhangrui@nuc.edu.cn zhangrui@nuc.edu.cn 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 中北大学 山西省光电信息与仪器工程技术研究中心，山西 太原 030051 中北大学 山西省光电信息与仪器工程技术研究中心，山西 太原 030051 Engineering and Technology Research Center of Shanxi Provincial for Optical-Electric Information and Instrument， North University of China， Taiyuan 030051， China Engineering and Technology Research Center of Shanxi Provincial for Optical-Electric Information and Instrument， North University of China， Taiyuan 030051， China 中北大学 前沿交叉科学研究院，山西 太原 030051 中北大学 前沿交叉科学研究院，山西 太原 030051 Academy for Advanced Interdisciplinary Research， North University of China， Taiyuan 030051， China Academy for Advanced Interdisciplinary Research， North University of China， Taiyuan 030051， China 中北大学 南通智能光机电研究院，江苏 南通 710065 中北大学 南通智能光机电研究院，江苏 南通 710065 Nantong Institute of Intelligent Opto-Mechatronics， North University of China， Nantong 710065， China Nantong Institute of Intelligent Opto-Mechatronics， North University of China， Nantong 710065， China To realized laser-warning systems for high-precision detection with a broadband and large field of view， a wide field two-dimensional laser-warning detection system was presented based on broadband black-white grating diffraction. In addition， an optical system comprising a broadband black-white grating and large field-of-view lens was designed. A five-group seven-element combined lens structure was applied to achieve the large field-of-view detection with an azimuth angle of 95° and a pitch angle of 72°. The lens aperture was 2.0， and the pupil diameter was 36 mm. These values ensured that the system exhibited a considerable amount of light transmission.Furthermore， a broadband black-white grating with a grating constant of 2.5 μm （and duty cycle of 0.6） was designed and processed， with broadband detection performed from 0.5-1.7 μm. Moreover， a multi-narrowband filter is added to the system to filter out background-light interference. Finally， testing with lasers of typical wavelengths （532， 1 064， and 1 550 nm）， the angle error achieved by the system is better than 0.7°， the laser center wavelength detection accuracy is better than 8 nm， thus verifying the feasibility of the method experimentally. 为了实现激光告警机宽光谱、大视场的高精度探测，提出了一种基于宽波段光栅衍射的大视场二维激光告警探测系统。设计了由宽波段光栅和大视场镜头组成的光学系统，通过采用五组七片式组合透镜结构，实现了水平方向95°、垂直方向72°的大视场探测。镜头的光圈为2.0，入瞳直径为36 mm，保证系统具有很好的通光量。设计加工了光栅常数为2.5 μm（占空比为0.6）的宽波段黑白光栅，实现了0.5~1.7 μm的宽波段探测。系统中加入多窄带滤光片以滤除背景光的干扰。最后，用典型波长（532，1 064，1 550 nm）的激光器进行测试，系统探测得到的角度测量误差优于0.7°、激光中心波长的探测精度优于8 nm，通过实验验证了该方法的可行性。 激光告警 光学设计 宽波段黑白光栅 激光波长 高精度二维角度 laser warning optical design broadband black-white grating laser wavelength high-precision two-dimensional angle 1 Introduction 1 The detection of the parameters of incoming laser is associated with an important and constantly expanding group of security， military， aircraft， and satellite safety systems ［ 1 1 ］ . Over the last century， the development of laser detection technology has become one of the most debated issues surrounding the protection of strategic centers against various laser sources， such as continuous-wave lasers and pulsed lasers ［ 2 2 ］ . Therefore， a laser detection system has been designed to identify a threat， collect information accurately， and issue warnings promptly. Furthermore， this system can combine excellent electronic devices with advanced optical technologies ［ 2 2 ］ . Recently proposed laser detection systems can be classified into three categories： coherent recognition， scatter recognition， and spectrum recognition ［ 3 3 - 5 5 3-5 ］ . Spectrum recognition can also be categorized depending on the involvement of imaging. Primarily， the imaging laser detection system consists of a wide-angle fisheye lens and a charge-coupled device （CCD） photodetector ［ 6 6 - 7 7 6-7 ］ . The most significant advantage of this method is that full air space can be monitored with high orientation accuracy， but without scanning ［ 8 8 ］ . However， this detection system cannot identify the laser wavelength， which is detrimental to determining the type of a laser threat ［ 8 8 ］ . Another methodology involves the application of quadrant photodetectors or discrete point detectors located within a selected geometry ［ 9 9 ］ . For sensors of this type， the ratios of received powers （and related photocurrents） accommodate the evaluation of the incoming radiation for its angle ［ 9 9 ］ . However， in most circumstances， such solutions suffer from a limited field of view， low angular resolution， susceptibility to the background， and heavy reliance on beam irradiance distribution ［ 10 10 ］ . Furthermore， despite being perfectly Gaussian in the early stages， a laser beam is affected by refractive index spatial inhomogeneities （turbulence） after travelling hundreds of meters in the atmosphere. This phenomenon can distort the wave front and lead to irradiance inhomogeneities ［ 10 10 ］ . Recently， a novel method was proposed to recognize the incident angle and laser wavelength in a large field of view （FOV） ［ 10 10 ］ . In this system， a large field-of-view lens is combined with a sinusoidal amplitude grating ［ 11 11 ］ . First， the system's directional resolution is determined by the angular accuracy of the hardware characteristics exhibited by imaging devices， such as the FOV of the lens and the resolution of imaging ［ 12 12 ］ . Subsequently， a lens is designed to cover this wavelength range， which is significant for improving the capability of the imaging laser detection system ［ 12 12 ］ . The objective of this article is to provide a lens and grating for an imaging laser detection system ［ 13 13 ］ that can supposedly detect rays in both the visible and NIR spectra （wavelength of laser threats） ［ 14 14 ］ . The inverse-telephoto objective lens was applied as a first-order design to reduce the angle of the incident ray ［ 14 14 ］ . Finally， the results of the proposed structure were examined using ZEMAX optical design software ［ 15 15 ］ . Based on the free mirror proposed in the Netherlands last year， the laser-warning system has a detection accuracy of more than 1°， which falls far short of the laser detection requirement ［ 15 15 ］ . Here， a large field-of-view lens and wide-spectrum grating were taken as the primary originals to achieve high-precision two-dimensional laser warnings. 2 The diagram of the laser warning 1 The laser-warning system adopted a modular structure， as shown in Figure 1 Figure 1 . Therefore， it is advantageous for comprehensively detecting the incoming laser owing to its wavelength， azimuth， pitch angle， and pulse width. Figure 2 Figure 2 shows the optical path diagram of the detection of the incoming laser based on broadband grating. The system consists mainly of broadband black-white grating， a lens， a multi-narrowband filter， and a visible and near-infrared area array detector （CCD）， which is placed on the focal plane of the lens. This diagram shows that the CCD imager records the spots once the lens detects the laser beam. Subsequently， the direction and wavelength of the incoming laser beam can be determined using the coordinates of the imaged spots through an image processing circuit system. A black-white grating diffracts the incoming ray while the lens is imaged on the area array detector. Figure 3 Figure 3 （a） shows the projection of the incoming laser beam along the xz -plane， which is used to measure the azimuth angle α . The projection in the yz -plane is used to measure the pitch angle γ ， as shown in Figure 3 Figure 3 （b）. After entering the optical system， the incoming beam is diffracted by a broadband black-white grating before passing through a broadband， large field of view， achromatic lens to converge the diffracted beam on the CCD at the focal plane and image， thus forming a light spot. As shown in Figure 3 Figure 3 （a）， the position of the 0 diffraction order spectrum is denoted by x 0 ， that of the 1 diffraction order spectrum is indicated by x + ， that of the - 1 diffraction order spectrum is indicated by x - ， the 0 diffraction order angle is represented by α ， and the 1 and - 1 diffraction order angles are expressed as β + and β - ， respectively， which can be obtained using Equation （1） Equation （1） . <math display="block" id="M1"><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><msub><mrow><mi>β</mi></mrow><mrow><mo>+</mo></mrow></msub><mo>=</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>α</mi><mo>+</mo><mfrac><mrow><mi>λ</mi></mrow><mrow><mi>d</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><msub><mrow><mi>β</mi></mrow><mrow><mo>-</mo></mrow></msub><mo>=</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>α</mi><mo>-</mo><mfrac><mrow><mi>λ</mi></mrow><mrow><mi>d</mi></mrow></mfrac></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037423&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037421&type= 29.88733101 18.79599953 . The CCD was placed at the focal plane of the lens， and the spot information， as imaged on the CCD， was processed. Thus， the positions of x 0 ， x + ， and x - can be obtained as follows： <math display="block" id="M2"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mi>f</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mi>α</mi></mtd></mtr><mtr><mtd><msub><mrow><mi>x</mi></mrow><mrow><mo>+</mo></mrow></msub><mo>=</mo><mi>f</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><msub><mrow><mi>β</mi></mrow><mrow><mo>+</mo></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>x</mi></mrow><mrow><mo>-</mo></mrow></msub><mo>=</mo><mi>f</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><msub><mrow><mi>β</mi></mrow><mrow><mo>-</mo></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037428&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037425&type= 21.25133324 16.08666801 . Considering Equations （ 1 1 ） and （ 2 2 ）： <math display="block" id="M3"><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mi>α</mi><mo>=</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfrac><mrow><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>f</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mi>λ</mi><mo>=</mo><mi>d</mi><mfenced open="[" close="]" separators="|"><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfrac><mrow><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>f</mi></mrow></mfrac><mo stretchy="false">)</mo><mo>+</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mfenced separators="|"><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfrac><mrow><msub><mrow><mi>x</mi></mrow><mrow><mo>+</mo></mrow></msub></mrow><mrow><mi>f</mi></mrow></mfrac></mrow></mfenced></mrow></mfenced></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037432&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037430&type= 63.58466721 22.35199928 . The angle of interest， γ， （ Figure 2 Figure 2 （b）） can be inferred using the same principle： <math display="block" id="M4"><mi>γ</mi><mo>=</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfrac><mrow><mi>y</mi></mrow><mrow><mi>f</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037437&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037435&type= 19.21933365 8.80533314 ， where d represents the grating period， and f refers to the focal length of the lens. 3 Optical design 1 3.1　Optical lens design 2 The laser-warning system studied here can measure the beam for its wavelength， azimuth angle （ α ）， and pitch angle （ γ ）. It can detect wavelengths of up to 0.5–1.7 μm， and the laser center wavelength detection accuracy is better than 8 nm. Meanwhile， the azimuth can attain 95°， whereas the pitch angle can attain 72°. Ultimately， the angular error is greater than 0.7°. The CCD chosen for this study was supplied by the Shanxi Guohui Optoelectronics Technology Co.， Ltd. Figure 4 Figure 4 shows some indices of the CCD that can detect the spectrum from visible to NIR wavelengths. This is effective for broadening the detection range. After the laser beam enters the general optical imaging system， a phase difference arises owing to the variation in lens thickness， thus achieving the focusing and imaging function of the lens. Moreover， there are certain requirements regarding the shape and thickness of the lens to improve its functioning. This study aimed to simultaneously detect laser beams with a broadband and a large field of view. In general， it is difficult to achieve this without an effective solution. Thus， a negative meniscus lens is required to reduce the maximum laser-beam angle and steer the beam away from the entire hemisphere to the end of the lens. Furthermore， multiple negative elements were used in the front group of the lenses to further reduce the angle of the laser beam. Therefore， a broadband and wide-field achromatic lens was designed. The ZEMAX optical automatic design software was applied according to the required parameters of the first-order model， and the merit function was used to perform the primary optimization. Thus， a high-resolution lens produces excellent performance in the context of visible and infrared spectra. Figure 5 Figure 5 illustrates the optical layout of the designed lens， which has a diagonal field of view of up to 150°. By contrast， Figure 6 Figure 6 shows the corresponding diffused spots of laser beams with different wavelengths and directions. Table 1 Table 1 shows the size of the diffused spots corresponding to different spectra. The size of the full-spectrum spot is smaller than 30 μm， and the spot radius meeting our requirements is smaller than the size of the two pixels of the detector. The root mean square （RMS） is a requisite parameter to evaluate an optical system's performance. According to this diagram， the maximum RMS is 10.1 μm， which means that the system can achieve the resolution required to provide high-quality images in the visible and near-infrared wavelengths range. Wavelength/nm Spot diameter/μm 500 29.3 532 25.8 632 19.6 1 064 20.6 1 550 23.7 1 700 23.2 Modulation transfer function （MTF）， one of the most comprehensive standards for evaluating lens quality， illustrates the modulation of an optical system from an object to an image as a function of spatial frequency， ranging between 0 and 1. Specific spatial frequency and aperture parameters were set to obtain the MTF images. The spatial frequency was set to 30 pairs/mm. The MTF curves are shown in Figure 7 Figure 7 . The closer the MTF value is to 1， the better the lens performance. Figure 7 Figure 7 shows the MTF characteristics of the lens for visible and infrared wavelengths. It can be observed from this figure that the MTF value of the lens invariably exceeded 0.4， indicating that the overall characteristics of the lens were satisfactory. A broadband， large field of view， achromatic lens was constructed based on this design， as shown in Figure 8 Figure 8 . The lens has an aperture of 2.0 and a pupil diameter of 36 mm. 3.2　Broadband black-white grating design 2 In line with these requirements， a broadband black-white grating of 32 mm × 32 mm and grating period of 2.5 μm was designed. The duty cycle of the designed grating is 0.6. A schematic of the grating structure is presented in Figure 9 Figure 9 . Compared with the traditional sinusoidal amplitude grating， this black-white grating has easy processing. The method discussed was implemented in the simulation. A multi-narrowband filter was used to filter out the interference of background light， and the filter was replaced according to the detection requirements. Without background-light interference， it is possible to add filters. Figure 10 Figure 10 shows the transmittance curve of the filter used in the experiment. MATLAB software was used to simulate the intensity of the laser beam at different wavelengths and incident angles. In the visible and infrared spectra， the laser beam wavelength was set to 532， 632， 1 064， and 1 550 nm for simulation. The simulation results obtained when the incident angle was 0° are shown in Figure 11 Figure 11 . For the field of view of 95°， both the 0 1 and - 1 diffraction order spots could be detected when the laser was incident at the maximum wavelength. The simulation results obtained for the incident angle of 48° are shown in Figure 12 Figure 12 ； a minimum of two spots are observed on the detector when the maximum wavelength of the laser is incident. According to the simulation diagram， the normalized intensity of ±1 diffraction order spots is in 50% excess， while the diffracted 0 order and ±1 order light intensities starkly differ. Consequently， the 0 and ±1 order spots can be effectively distinguished when the laser diffraction spot is imaged on the CCD. After processing the electrical signal on the CCD using a high-speed signal-processing circuit， information， such as the wavelength and direction， of the incident laser beam can be obtained. 4 Simulation and experiment 1 The broadband black-white grating period used here was 2.5 μm. The size of the CCD is 9.6 mm × 7.68 mm， the number of pixels is 640 × 512， and the area of each pixel is 15 μm × 15 μm. The results are shown in Figures 13 13 and 14 14 . The theoretical wavelength resolution of the incident laser beam is lower than 8 nm， while the angle resolution can attain 0.18°. In addition， this design can detect an azimuth of up to 95° and a pitch angle of up to 72° for an incident laser beam. However， given the experimental environment and algorithm errors， the angle error could still be maintained at a level higher than 0.7°. Multiple ultra-narrowband filters were used to filter the light source into monochromatic light for the simplicity of the experiment. The beam wavelengths were set to 532， 1 064， and 1 550 nm. The experimental equipment and results are shown in Figure 15 Figure 15 . As shown in Figure 15 Figure 15 （b）， 15（c）， and 15（d）， a comparison was performed for the spot positions after beam diffraction with different wavelengths at the same incident angle. Figure 15 Figure 15 （d） and 15（d） present a comparison of the spot positions after beam diffraction for the same laser wavelength but different incident angles. Regardless， the detector can detect the zero and one-diffraction order spots of the incident beam. Finally， after beam diffraction， the spot position information was processed by the DSP algorithm and high-speed signal processing circuit to obtain the angle， wavelength， and other details of the incident laser beam. A portion of the experimental data is presented in Table 2 Table 2 . The detection errors of the angle and wavelength can be expressed as follows： <math display="block" id="M5"><mtable columnalign="left"><mtr><mtd><mo>∆</mo><mi>α</mi><mo>=</mo><mfenced open="|" close="|" separators="|"><mrow><msup><mrow><mi>α</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>-</mo><mi>α</mi></mrow></mfenced></mtd></mtr><mtr><mtd><mo>∆</mo><mi>γ</mi><mo>=</mo><mfenced open="|" close="|" separators="|"><mrow><msup><mrow><mi>γ</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>-</mo><mi>γ</mi></mrow></mfenced></mtd></mtr><mtr><mtd><mo>∆</mo><mi>λ</mi><mo>=</mo><mfenced open="|" close="|" separators="|"><mrow><msup><mrow><mi>λ</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>-</mo><mi>λ</mi></mrow></mfenced></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037552&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037549&type= 20.99733353 16.42533302 ， 2 Wavelength/nm 2 Actual value/（°） 2 Measured value/（°） 2 Error/（°） 2 Wavelength and error/nm α γ α ' γ ' Δ α Δ γ λ ' Δ λ 2 532 39.0 22 39.6 22.3 0.6 0.3 536 4 -16 -4 -16.2 -3.4 0.2 0.6 536 4 2 1 064 10 25 10.1 25.0 0.1 0.0 1 066 2 20 -8 20.1 -8.1 0.1 0.1 1 060 4 2 1 550 -4 0 -3.9 0.6 0.1 0.6 1 557 7 -29 7 -29.1 7.4 0.1 0.4 1 557 7 where α and γ denote the azimuth and pitch angles of the turntable， respectively， and α ' and γ ' represent the measured actual azimuth and pitch angles， respectively. In addition， Δ α and Δ γ represent the azimuth and pitch angle errors， respectively. λ represents the central wavelength of the laser， λ ' represents the measured wavelength， and Δ λ represents the central laser wavelength error. It can be concluded from Table 2 Table 2 that the maximum angular error is 0.6° and the maximum central laser wavelength error is 7 nm. We can draw the following important conclusions from these values： the optical angle error is better than 0.7°. Simultaneously， the laser center-wavelength detection accuracy was better than 8 nm. Thus， high-precision laser warnings can be obtained. 5 Conclusion 1 Here， we designed an optical system with a broadband， large field of view， achromatic optical lens and broadband black and white grating to achieve high-precision detection of the two-dimensional angle of the incoming laser. According to the simulation and experimental results， the maximum measuring azimuth of the system could attain 95°， and the pitch angle could attain 72°. In addition， the detectable spectral range was expanded to 0.5–1.7 μm. Compared with the existing laser-warning system， this design achieved a large field of view， broadband， and high-precision detection simultaneously. The error angle was better than 0.7°， and the central laser wavelength detection accuracy was better than 8 nm. These effectively improved the reliability of laser warnings. In the future， the lens and grating designs will be optimized to further expand the detection field of view and improve detection accuracy. References XU C Y ， CHENG H B ， FENG Y P ， et al . Large field distributed aperture laser semiactive angle measurement system design with imaging fiber bundles ［J］. Applied Optics ， 2016 ， 55 （ 25 ）： 7000 - 7008 . doi: 10.1364/ao.55.007000 http://dx.doi.org/10.1364/ao.55.007000 Large field distributed aperture laser semiactive angle measurement system design with imaging fiber bundles YU J ， SHEN Z X ， WANG Z S . Compact dual band/dual FOV infrared imaging system with freeform prism ［J］. Optics Letters ， 2021 ， 46 （ 4 ）： 829 - 832 . doi: 10.1364/ol.412091 http://dx.doi.org/10.1364/ol.412091 Compact dual band/dual FOV infrared imaging system with freeform prism NEJAD S M ， SHESHKELANI N R . Design and performance analysis of a fisheye-based optical head for an imaging laser detecting system ［J］. Iranian Journal of Science and Technology， Transactions A： Science ， 2020 ， 44 （ 5 ）： 1595 - 1604 . doi: 10.1007/s40995-020-00966-5 http://dx.doi.org/10.1007/s40995-020-00966-5 Design and performance analysis of a fisheye-based optical head for an imaging laser detecting system EL-SHERIF A F ， AYOUB H S ， EL-SHARKAWY Y H ， et al . The design and implementation of photoacoustic based laser warning receiver for harsh environments ［J］. Optics & Laser Technology ， 2018 ， 98 ： 385 - 396 . doi: 10.1016/j.optlastec.2017.06.030 http://dx.doi.org/10.1016/j.optlastec.2017.06.030 The design and implementation of photoacoustic based laser warning receiver for harsh environments ZANDI M ， SUGDEN K ， BENTON D M . Low-cost laser detection system with a 360-deg horizontal field of view ［J］. Optical Engineering ， 2021 ， 60 （ 2 ）： 027106 . doi: 10.1117/1.oe.60.2.027106 http://dx.doi.org/10.1117/1.oe.60.2.027106 Low-cost laser detection system with a 360-deg horizontal field of view MIKŠ A ， POKORNÝ P . Calculation of a lens system with one or two aspherical surfaces having corrected spherical aberration ［J］. Journal of the Optical Society of America A， Optics， Image Science， and Vision ， 2020 ， 37 （ 9 ）： 1390 - 1397 . doi: 10.1364/josaa.399361 http://dx.doi.org/10.1364/josaa.399361 Calculation of a lens system with one or two aspherical surfaces having corrected spherical aberration YANG W J ， ZHANG X H ， MA H W ， et al . Laser beams-based localization methods for boom-type roadheader using underground camera non-uniform blur model ［J］. IEEE Access ， 2020 ， 8 ： 190327 - 190341 . doi: 10.1109/access.2020.3032368 http://dx.doi.org/10.1109/access.2020.3032368 Laser beams-based localization methods for boom-type roadheader using underground camera non-uniform blur model MOHAMMADNEJAD S ， AGHAEI F . Noise characteristics improvement of submicron InP/InGaAs avalanche photodiode for laser detection system ［J］. Optics Communications ， 2020 ， 455 ： 124561 . doi: 10.1016/j.optcom.2019.124561 http://dx.doi.org/10.1016/j.optcom.2019.124561 Noise characteristics improvement of submicron InP/InGaAs avalanche photodiode for laser detection system LIU S ， NIU Y X . Design and analysis of miniature high-resolution omnidirectional gaze optical imaging system ［J］. Optical Review ， 2017 ， 24 （ 3 ）： 334 - 344 . doi: 10.1007/s10043-017-0324-0 http://dx.doi.org/10.1007/s10043-017-0324-0 Design and analysis of miniature high-resolution omnidirectional gaze optical imaging system 杨皓琨 ， 赵冬娥 ， 刘吉 ， 等 . 基于衍射光栅的光谱检测系统的光学设计 ［J］. 科技视界 ， 2016 （ 4 ）： 161 - 162 . doi: 10.3969/j.issn.2095-2457.2016.04.116 http://dx.doi.org/10.3969/j.issn.2095-2457.2016.04.116 基于衍射光栅的光谱检测系统的光学设计 YANG H K ， ZHAO D E ， LIU J ， et al . Design of spectrum detecting system based on diffraction grating ［J］. Science ＆ Technology Vision ， 2016 （ 4 ）： 161 - 162 . （in Chinese） . doi: 10.3969/j.issn.2095-2457.2016.04.116 http://dx.doi.org/10.3969/j.issn.2095-2457.2016.04.116 Design of spectrum detecting system based on diffraction grating 李晓 ， 张瑞 ， 王志斌 ， 等 . 二维激光告警光学系统设计 ［J］. 红外与激光工程 ， 2015 ， 44 （ 6 ）： 1806 - 1810 . doi: 10.3969/j.issn.1007-2276.2015.06.023 http://dx.doi.org/10.3969/j.issn.1007-2276.2015.06.023 二维激光告警光学系统设计 LI X ， ZHANG R ， WANG ZH B ， et al . Design of two-dimensional laser warning optical system ［J］. Infrared and Laser Engineering ， 2015 ， 44 （ 6 ）： 1806 - 1810 . （in Chinese） . doi: 10.3969/j.issn.1007-2276.2015.06.023 http://dx.doi.org/10.3969/j.issn.1007-2276.2015.06.023 Design of two-dimensional laser warning optical system 曹俊卿 ， 赵冬娥 ， 王志斌 ， 等 . 一种光栅型激光告警系统设计与实现 ［J］. 激光与红外 ， 2015 ， 45 （ 1 ）： 5 - 8 . doi: 10.3969/j.issn.1001-5078.2015.01.001 http://dx.doi.org/10.3969/j.issn.1001-5078.2015.01.001 一种光栅型激光告警系统设计与实现 CAO J Q ， ZHAO D E ， WANG ZH B ， et al . Design and implementation of a grating laser warning system ［J］. Laser ＆ Infrared ， 2015 ， 45 （ 1 ）： 5 - 8 . （in Chinese） . doi: 10.3969/j.issn.1001-5078.2015.01.001 http://dx.doi.org/10.3969/j.issn.1001-5078.2015.01.001 Design and implementation of a grating laser warning system WANG L ， DONG H J ， ZHANG W A ， et al . Optimum design of wide angle laser detecting system based on fish eye lens and sinusoidal amplitude grating ［J］. Optics Communications ， 2014 ， 310 ： 173 - 178 . doi: 10.1016/j.optcom.2013.07.084 http://dx.doi.org/10.1016/j.optcom.2013.07.084 Optimum design of wide angle laser detecting system based on fish eye lens and sinusoidal amplitude grating GUTIÉRREZ C E ， SABRA A . Freeform lens design for scattering data with general radiant fields ［J］. Archive for Rational Mechanics and Analysis ， 2018 ， 228 （ 2 ）： 341 - 399 . doi: 10.1007/s00205-017-1196-y http://dx.doi.org/10.1007/s00205-017-1196-y Freeform lens design for scattering data with general radiant fields WOJTANOWSKI J ， JAKUBASZEK M ， ZYGMUNT M . Freeform mirror design for novel laser warning receivers and laser angle of incidence sensors ［J］. Sensors （Basel， Switzerland） ， 2020 ， 20 （ 9 ）： 2569 . doi: 10.3390/s20092569 http://dx.doi.org/10.3390/s20092569 Freeform mirror design for novel laser warning receivers and laser angle of incidence sensors 10.37188/OPE.20223019.2295 TH254 A 1004-924X（2022）19-2295-10 杨雪梅,张瑞,石金等.基于黑白光栅的宽波段大视场激光告警[J].光学精密工程,2022,30(19):2295-2304. 杨雪梅,张瑞,石金等.基于黑白光栅的宽波段大视场激光告警[J].光学精密工程,2022,30(19):2295-2304. YANG Xuemei,ZHANG Rui,SHI Jin,et al.Laser warning with broadband and large field of view based on black-white grating[J].Optics and Precision Engineering,2022,30(19):2295-2304. YANG Xuemei,ZHANG Rui,SHI Jin,et al.Laser warning with broadband and large field of view based on black-white grating[J].Optics and Precision Engineering,2022,30(19):2295-2304. 国家自然科学基金资助项目（No.62105302）；山西省自然科学基金资助项目（No.20210302124269， No.201901D111158， No.201901D111145）；山西省创新项目（No.2019L0530） 2022-04-09 2022-06-07 2022-10-10 2022-10-24 光学精密工程 19 30