CN114216575A - Ultra-short pulse reconstruction system and method based on BP neural network - Google Patents
Ultra-short pulse reconstruction system and method based on BP neural network Download PDFInfo
- Publication number
- CN114216575A CN114216575A CN202111383297.0A CN202111383297A CN114216575A CN 114216575 A CN114216575 A CN 114216575A CN 202111383297 A CN202111383297 A CN 202111383297A CN 114216575 A CN114216575 A CN 114216575A
- Authority
- CN
- China
- Prior art keywords
- pulse
- neural network
- data
- scanning spectrum
- ultrashort
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000013078 crystal Substances 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 22
- 238000010586 diagram Methods 0.000 claims abstract description 17
- 230000003287 optical effect Effects 0.000 claims abstract description 17
- 238000003062 neural network model Methods 0.000 claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims description 47
- 239000013598 vector Substances 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 6
- QBLDFAIABQKINO-UHFFFAOYSA-N barium borate Chemical group [Ba+2].[O-]B=O.[O-]B=O QBLDFAIABQKINO-UHFFFAOYSA-N 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 abstract description 5
- 238000012512 characterization method Methods 0.000 abstract description 3
- 238000013135 deep learning Methods 0.000 abstract description 2
- 230000003595 spectral effect Effects 0.000 abstract description 2
- 230000003111 delayed effect Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 6
- 210000002569 neuron Anatomy 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000009407 construction method and process Methods 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J11/00—Measuring the characteristics of individual optical pulses or of optical pulse trains
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/02—Details
- G01J3/0205—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
- G01J3/0208—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using focussing or collimating elements, e.g. lenses or mirrors; performing aberration correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/02—Details
- G01J3/0205—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
- G01J3/021—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using plane or convex mirrors, parallel phase plates, or particular reflectors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/42—Absorption spectrometry; Double beam spectrometry; Flicker spectrometry; Reflection spectrometry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/42—Absorption spectrometry; Double beam spectrometry; Flicker spectrometry; Reflection spectrometry
- G01J2003/423—Spectral arrangements using lasers, e.g. tunable
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
Abstract
本发明属于超短脉冲表征技术领域，涉及基于深度学习BP神经网络的超短脉冲重建系统及方法。一种基于BP神经网络的超短脉冲重建系统，包括双光束自相关光路单元、主控单元；所述的双光束自相关光路单元用于获取超短脉冲的延迟扫描光谱行迹图；所述的主控单元用于系统的控制及数据处理，包括控制模块和数据处理模块；所述的数据处理模块包括BP神经网络模型；所述BP神经网络模型对输入的扫描光谱行迹图进行预测，输出待测超短脉冲的幅度和相位。本发明提出了针对待测脉冲波长的不同而进行BBO晶体的切割角计算，提高了系统的适用性和实验的成功率。结合BP神经网络重建待测脉冲信息，极大的提高了超短脉冲重建效率。
The invention belongs to the technical field of ultrashort pulse characterization, and relates to an ultrashort pulse reconstruction system and method based on a deep learning BP neural network. An ultra-short pulse reconstruction system based on a BP neural network, comprising a double-beam autocorrelation optical path unit and a main control unit; the double-beam auto-correlation optical path unit is used to obtain a delayed scanning spectral track diagram of an ultra-short pulse; the The main control unit is used for system control and data processing, including a control module and a data processing module; the data processing module includes a BP neural network model; Measure the amplitude and phase of ultrashort pulses. The invention proposes to calculate the cutting angle of the BBO crystal according to the different wavelengths of the pulses to be measured, which improves the applicability of the system and the success rate of the experiment. Combined with the BP neural network to reconstruct the pulse information to be measured, the reconstruction efficiency of ultra-short pulses is greatly improved.
Description
Technical Field
The invention belongs to the technical field of ultrashort pulse characterization, and relates to an ultrashort pulse reconstruction system and method based on deep learning BP neural network
Background
Since the advent of ultrashort pulse laser, the ultrashort pulse laser has irreplaceable important application in numerous fields such as microscopic imaging, laser processing, medical treatment and the like, has attracted extensive attention and great interest of society, and particularly has very key significance in the research of the quantum field.
The ultrashort pulse laser has very short laser pulse width, generally in picosecond or femtosecond level, and when the pulse width is femtosecond level, the pulse width exceeds the electronic response limit and cannot be represented by a general detector, so that the method for accurately and efficiently representing ultrashort pulses is very critical. In the sixties of the twentieth century, the problem of ultrashort pulse characterization was solved to some extent by the proposal of the autocorrelation method, and in the following decades, the self-referenced spectral interferometry, the frequency-resolved optical switching method, and the like were successively proposed. Compared with the former two methods, the frequency-resolved optical switching method, especially the second harmonic frequency-resolved optical switching method, has the characteristics of high resolution, high accuracy, simple optical path and the like, and is a mainstream method for representing ultrashort pulses at present, but the method needs a long time to calculate and restore the required data result from the spectrum trace, cannot detect in real time, has low signal-to-noise ratio and poor performance, and is widely subject to defects.
Disclosure of Invention
Aiming at the defects of low resolution, low accuracy (autocorrelation method), complex light path (self-reference spectrum interference method), low efficiency (frequency-resolved optical switching method), poor performance (all) of low signal-to-noise ratio and the like in the prior art, the invention provides a system for representing ultrashort pulses by using a double-reflection second harmonic frequency-resolved optical switching method based on a BP neural network, preferably selects an optical part of the system consisting of the double-reflection autocorrelation light path and a second harmonic nonlinear crystal, and uses the BP neural network to replace the traditional inversion algorithm (such as PCGPA) to calculate and obtain amplitude phase data.
The technical scheme adopted by the invention for solving the technical problems is as follows: the ultrashort pulse reconstruction system based on the BP neural network is characterized in that: the system comprises a double-beam autocorrelation light path unit and a main control unit; the double-beam autocorrelation light path unit is used for acquiring a scanning spectrum trace diagram of the ultrashort pulse; the main control unit is used for controlling the system and processing data and comprises a control module and a data processing module; the data processing module comprises a BP neural network model; and the BP neural network model predicts the input scanning spectrum trace graph and outputs the amplitude and the phase of the ultrashort pulse to be detected.
Preferably, the training set of the BP neural network model is composed of simulated scanning spectrum trace graphs, and the same group of noiseless pulse amplitudes and phases correspond to n simulated scanning spectrum trace graphs with different signal-to-noise ratios; n is more than or equal to 1 and less than or equal to 30.
Preferably, the simulated scanning spectrum trace diagram is obtained by the following steps:
(1) collecting pulse amplitude data, and carrying out N-point sampling on each group of collected pulse amplitude data, wherein N is more than or equal to 100 and less than or equal to 500;
(2) selecting N data from N sampling points of the pulse amplitude data, and performing random number processing of different degrees;
(3) and (3) according to the formula, the pulse amplitude data processed in the step (2) and the pulse phase data generated by the function:
in the formula
The spectrum is a scanning spectrum of second harmonic, E (t) is a pulse to be detected, E (t-tau) is a time-shifting copy, namely a gate function, of the pulse to be detected, omega is the frequency of the pulse to be detected, and tau is delay;
and calculating to obtain data of the simulated scanning spectrum trace graph, and generating n simulated scanning spectrum trace graphs with different signal-to-noise ratios.
Preferably, the different degrees of random numeralization processing are: respectively taking different random numbers between 0 and 1 for the selected n pulse amplitude data; and, the selected random number is different from the corresponding original value.
Preferably, the control module is configured to control a delay system of the dual-beam autocorrelation optical path unit to delay.
Preferably, the control module is further configured to control the spectrometer to continuously scan the frequency-doubled pulse spectrum while the delay system generates the delay, so as to generate a delay-scanned spectrum trace pattern.
Preferably, the nonlinear crystal in the two-beam autocorrelation optical path unit is a barium metaborate crystal; the cutting angle of the barium metaborate crystal is a main section collinear phase matching frequency doubling angle, and the calculation method comprises the following steps: setting two beams of pulse to be measured as k1And k2The frequency multiplication pulse is k3The phase matching condition under the condition of the main section is to satisfy the wave vector superposition, i.e.
Decomposing the vector to obtain a formula:
since the three pulses satisfy the vector superposition, so
The square addition operation is performed on the formulas (2) and (3), and the result is substituted into
To obtain:
according to the frequency multiplication, omega can be known3＝2ω1＝2ω2Therefore, it is
Wherein n isp(θ3) The angle between the frequency-doubled light and the main axis is theta3The refractive index of the time-doubled light; n isp(θ1)、np(θ2) For the pulse to be measured and the main shaft to form an included angle theta1、θ2The refractive index of the pulse to be measured;
the optical refractive index of the frequency doubling light e is calculated according to a formula:
wherein n isp(θ3) Is that
Combining the formulas (5) and (6), and drawing by matlab to obtain theta3With theta1The change map of (1), theta corresponding to the vertex of the image3The value is the desired cutting angle.
The invention also provides an ultrashort pulse reconstruction method based on the BP neural network, which comprises the following steps:
(1) acquiring an ultrashort pulse slow scanning spectrum trace diagram to be detected;
(2) and inputting the scanning spectrum trace graph into a pre-trained BP neural network model, and outputting the amplitude and the phase of the ultrashort pulse to be detected.
The ultra-short pulse reconstruction system and method based on the BP neural network have the advantages that: on the basis of measuring ultrashort pulses by the traditional FROG method, double reflection improvement is carried out on an autocorrelation light path, and mutual interference among the light paths is reduced. The cutting angle calculation of the BBO crystal is carried out according to the difference of the pulse wavelengths to be detected, and the applicability of the system and the success rate of the experiment are improved. The BP neural network is combined to replace traditional algorithms such as PCGPA and the like to reconstruct pulse information to be detected from the FROG trace, and the neural network is trained by using training sets with different signal-to-noise ratios, so that the ultra-short pulse reconstruction efficiency is greatly improved, and the accuracy of reconstructing ultra-short pulses under the condition of low signal-to-noise ratio is improved.
Drawings
Fig. 1 is a schematic structural diagram of an ultra-short pulse reconstruction system based on a BP neural network according to an embodiment of the present invention;
FIG. 2 is a BBO crystal collinear phase matching calculation result when the pulse wavelength to be measured is 800 nm;
FIG. 3 is a schematic diagram of a BP neural network structure;
FIG. 4 is a schematic diagram of a random number process;
fig. 5 is a flowchart of an ultra-short pulse reconstruction method based on a BP neural network according to an embodiment of the present invention;
in fig. 1: 1. a femtosecond laser; 2-1. a first silvered mirror; 2-2 a second silvered mirror; 2-3 a third silvered mirror; 3. a femtosecond beam splitter; 4. a converging mirror; 5. a nonlinear crystal; 6. a diaphragm; 7. a spectrometer; 8. a stepping motor; 9. a computer; 10 external electrical pulse generator.
Detailed Description
In order to facilitate an understanding of the invention, the invention is described in more detail below with reference to the accompanying drawings and specific examples. Preferred embodiments of the present invention are shown in the drawings. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.
The structure of a dual-beam autocorrelation optical path system with a BBO crystal as a nonlinear crystal is shown in fig. 1, and includes: the device comprises a femtosecond laser 1, 5 silvered reflectors, a femtosecond beam splitter 3, a converging reflector 4, a nonlinear crystal 5, a data acquisition system, a delay system and a delay system controller.
The data acquisition system consists of a spectrometer 7 and a computer 9. The delay system consists of two second silvered reflectors 2-2 and a stepping motor 8, a delay system controller consists of an external electric pulse generator 10 and a control module built in a computer 9 and is used for controlling the delay system to generate displacement, and a BP neural network model is built in the computer 9.
The fixed optical path is composed of two first silvered mirrors 2-1. The emergent light of the fixed light path is incident on the converging mirror 4. The emergent light of the delay system enters the converging reflector 4 after passing through the third silvered reflector.
The nonlinear crystal 5 is used for generating frequency doubling pulse signals, and the diaphragm 6 is used for blocking other pulses except the frequency doubling pulses.
In order to facilitate movement, the light path part and the stepping motor of the embodiment are arranged on a 25-by-25-hole aluminum honeycomb plate, the stepping motor 8 is connected with an external pulse generator 10 through a data line, the external pulse generator 10 is connected with a computer 9 through an RS232 data bus, the spectrometer 7 is connected with the computer 9 through an instrument control line, and all instruments in the system are controlled by the computer 9.
The femtosecond laser device 1 outputs ultrashort pulses to be detected, and the femtosecond beam splitter 3 is 1: 1 splitting ratio, and dividing the ultrashort pulse to be detected into two input delay systems and a fixed light path. The fixed light path and the delay system are reflected twice to separate the incident light beam from the emergent light beam.
The data acquisition system acquires a scanning spectrum by a spectrometer 7, and generates a scanning spectrum trace pattern by a computer 9. The delay system generates optical path delay through the stepping motor 8, the delay system controller is composed of a control module in the computer 9 and an external electric pulse generator 10, the control module in the computer 9 writes a program for controlling delay amount by labview and controls the external electric pulse generator 10 to output pulse electric signals to the stepping motor 8.
In this embodiment, the nonlinear crystal 5 is a barium metaborate crystal (BBO crystal) having a large phase matching angle. The cutting angle of the BBO crystal is calculated according to the wavelength to be measured, and the calculation method comprises the following steps:
setting two beams of pulse to be measured as k1And k2The frequency multiplication pulse is k3The phase matching condition under the condition of the main section is that the wave vectors meet the vector superposition, namely
Decomposing the vector to obtain a formula:
since the three pulses satisfy the vector superposition, so
The square addition operation is performed on the formulas (2) and (3), and the result is substituted into
To obtain:
according to the frequency multiplication, omega can be known3＝2ω1＝2ω2Therefore, it is
Wherein n isp(θ3) The angle between the frequency-doubled light and the main axis is theta3The refractive index of the time-doubled light; n isp(θ1)、np(θ2) For the pulse to be measured and the main shaft to form an included angle theta1、θ2The refractive index of the pulse to be measured; in this embodiment, a class of phase matching ooes is adopted, the refractive index of the pulse o to be measured does not change with the included angle of the main shaft, and the refractive index of the frequency doubling e is calculated according to a formula:
wherein n isp(θ3) Is that
Combining the formulas (5) and (6), and drawing by matlab to obtain theta3With theta1The change map of (1), theta corresponding to the vertex of the image3The value is the desired cutting angle.
Compared with the prior art, the nonlinear crystal cut angle calculation method provided by the embodiment is more intuitive, and phase matching types such as ooe, eoe and the like can be calculated. The angle of the pulse to be measured entering the nonlinear crystal in a real experiment is difficult to control, and the calculation is based on the fact that the pulse to be measured can generate frequency doubling pulses when entering the nonlinear crystal at any angle, so that the feasibility of the experiment is improved.
In this embodiment, the wavelength of the femtosecond laser is 800nm, the cutting angle of the nonlinear crystal is a main-section collinear phase matching frequency doubling angle calculated to be 29.3 °, and the calculation result is shown in fig. 2.
In this embodiment, the working principle of the hardware system is as follows:
the femtosecond pulse to be detected is divided into light intensity of 1: 1, one beam of the same pulse is incident to a converging reflector after passing through a fixed light path, the other beam of the same pulse is reflected to the converging reflector through a silver-plated reflector after passing through a delay system, the two beams of pulses to be detected are converged to the center of a BBO crystal by the converging reflector, a frequency doubling pulse is generated between the two pulses after the nonlinear effect of the BBO crystal, and the pinhole diaphragm blocks other pulses except the frequency doubling pulse so as to ensure that only the frequency doubling pulse is incident to a probe of a spectrometer. The spectrometer transmits the collected frequency doubling pulse signals to a computer, and a scanning spectrum trace pattern is obtained through the processing of the computer.
In this embodiment, the software system selects a back-propagation (back-propagation) neural network structure, and compared with a multi-layer perceptron, the back-propagation optimization parameter weight is adopted by the BP neural network, so that the accuracy is higher and more reliable, and the convolutional neural network structure is simpler.
Setting the neuron number of a model input layer to be 49152, adding an input layer bias unit to be 1, setting the neuron number of a hidden layer to be 49152, adding a hidden layer bias unit to be 1, setting the neuron number of the hidden layer to be 49152 according to a picture matrix with a to-be-detected scanning spectrum trace graph of 128 × 3, setting the neuron number of an output layer to be 256 according to 128 points of amplitude data and phase data of two types of output pulses, and introducing a nonlinear factor by using a relu activation function to complete the neural network setting, as shown in fig. 3.
In this embodiment, the training set of the BP neural network model adopts a constructed simulated scanning spectrum trace atlas, and the construction method and process of the training set are as follows:
(A) placing a femtosecond laser with the wavelength of 800nm to be detected in front of the built double-beam autocorrelation light path system, and adjusting light path collimation to focus two beams of pulses to be detected on a BBO crystal to generate second harmonic.
(B) The step length displacement of a stepping motor is controlled by using a program compiled by labview, delay is realized, a spectrometer collects the scanning spectrum of each displacement, and a scanning spectrum tracing diagram is generated at a computer end.
(C) Collecting pulse amplitude data from documents and experiments, and performing 128-point sampling on each group of collected pulse amplitude data;
(D) selecting 1-30 data from 128 sampling points of the amplitude data to carry out random number processing with different degrees, specifically:
as shown in fig. 4a, the collected pulse amplitude data is a gaussian pulse, which is a continuous gaussian curve, and 128 points are equally divided on the abscissa, and the value of the middle point of each part represents the part, so as to obtain a scatter diagram shown in fig. 4b, which is composed of 128 points and is subjected to random number processing of different degrees, as shown in fig. 4c, 6 data points are subjected to random number processing, and as shown in fig. 4c, the data of 27 th, 45 th, 61 th, 70 th, 81 th, and 101 th points are given a value between 0 and 1 which is different from the original value, so as to obtain a scatter diagram shown in fig. 4c, which is the random number processing.
(E) And (3) the 6 processed pulse amplitude data and the pulse phase data generated by the function are processed according to a formula:
in the formula
The second harmonic scanning spectrum, E (t), the pulse to be measured, E (t-tau) the time-shifted copy of the pulse to be measured, i.e. the gate function, omega the frequency of the pulse to be measured, tau the delay.
And calculating to obtain data of the simulated scanning spectrum trace graph, and generating 6 simulated scanning spectrum trace graphs with different signal-to-noise ratios. FIG. 4d is a plot of a scatter plot as shown in FIG. 4 c.
Through the random number processing, the generated simulated scanning spectrum trace graph presents different signal-to-noise ratios, and a real spectrum trace graph can be simulated better.
In the training set in this embodiment, if 30 data are selected for randomization, the same group of noiseless pulse amplitudes and phases correspond to 30 simulated scanning spectrum trace patterns with different signal-to-noise ratios.
In this embodiment, a spectrometer is used to collect a delay scanning spectrum trace pattern of a frequency doubling pulse generated through an autocorrelation light path, and a PCGPA algorithm is used to perform inversion calculation to obtain amplitude pulse width phase information of the pulse to be measured, which is used as a test data set.
Inputting the constructed training data set into a BP neural network model, training and optimizing the model, wherein a loss function adopts cross entropy:
when loss is reduced to 10-3The expectation is considered to be reached.
And inputting the test data set and a part of the training data set into the trained model, predicting pulse information to be tested, and comparing the pulse information with the original data to test the performance of the BP neural network.
Through prediction and comparison, under the condition that the scanning spectrum trace diagram presents a high signal-to-noise ratio (is not subjected to random number processing), the reconstruction accuracy of the BP neural network is about 95.1% (95.3% of the traditional PCGPA), which is basically consistent with that of the traditional method, and under the condition that the scanning spectrum trace diagram presents a low signal-to-noise ratio (10 random number processing), the reconstruction accuracy of the BP neural network is 61.7% (40% of the traditional PCGPA), which is greatly improved compared with the traditional method. In the aspect of reconstruction speed, the time for reconstructing data from a group of scanning spectrum trace maps by the BP neural network is 0.5 second (about 30 minutes in the traditional PCGPA), which is greatly improved compared with the traditional method.
Embodiment 2 this embodiment provides an ultrashort pulse reconstruction method based on a BP neural network, and the flow is shown in fig. 5, including the following steps:
(1) acquiring a scanning spectrum trace diagram of the short wave pulse to be detected by using a built double-beam autocorrelation light path system which takes a BBO crystal as a nonlinear crystal;
(2) inputting the obtained scanning spectrum trace map into a built-in BP neural network model in a computer, and outputting corresponding amplitude and phase through model prediction.
Claims (8)
1. An ultrashort pulse reconstruction system based on a BP neural network is characterized in that: the system comprises a double-beam autocorrelation light path unit and a main control unit; the double-beam autocorrelation light path unit is used for acquiring a scanning spectrum trace diagram of the ultrashort pulse; the main control unit is used for controlling the system and processing data and comprises a control module and a data processing module; the data processing module comprises a BP neural network model; and the BP neural network model predicts the input scanning spectrum trace graph and outputs the amplitude and the phase of the ultrashort pulse to be detected.
2. The BP neural network-based ultrashort pulse reconstruction system of claim 1, wherein: the training set of the BP neural network model consists of simulated scanning spectrum trace graphs, and the same group of noiseless pulse amplitude and phase correspond to n simulated scanning spectrum trace graphs with different signal-to-noise ratios; n is more than or equal to 1 and less than or equal to 30.
3. The BP neural network-based ultrashort pulse reconstruction system of claim 2, wherein: the simulated scanning spectrum trace diagram is obtained by adopting the following steps:
(1) collecting pulse amplitude data, and carrying out N-point sampling on each group of collected pulse amplitude data, wherein N is more than or equal to 100 and less than or equal to 500;
(2) selecting N data from N sampling points of the pulse amplitude data, and performing random number processing of different degrees;
(3) and (3) according to the formula, the pulse amplitude data processed in the step (2) and the pulse phase data generated by the function:
in the formula
The spectrum is a scanning spectrum of second harmonic, E (t) is a pulse to be detected, E (t-tau) is a time-shifting copy, namely a gate function, of the pulse to be detected, omega is the frequency of the pulse to be detected, and tau is delay;
and calculating to obtain data of the simulated scanning spectrum trace graph, and generating n simulated scanning spectrum trace graphs with different signal-to-noise ratios.
4. The BP neural network-based ultrashort pulse reconstruction system of claim 3, wherein: the random number processing of different degrees is as follows: respectively taking different random numbers between 0 and 1 for the selected n pulse amplitude data; and, the selected random number is different from the corresponding original value.
5. The BP neural network-based ultrashort pulse reconstruction system of claim 1, wherein: the control module is used for controlling the delay system of the double-beam autocorrelation light path unit to delay.
6. The BP neural network-based ultrashort pulse reconstruction system of claim 5, wherein: the control module is also used for controlling the spectrometer to continuously scan the frequency doubling pulse spectrum while the delay system generates delay, and generating a delay scanning spectrum trace diagram.
7. The BP neural network-based ultrashort pulse reconstruction system of claim 1, wherein: the nonlinear crystal in the double-beam autocorrelation light path unit is a barium metaborate crystal; the cutting angle of the barium metaborate crystal is a main section collinear phase matching frequency doubling angle, and the calculation method comprises the following steps:
setting two beams of pulse to be measured as k1And k2The frequency multiplication pulse is k3The phase matching condition under the condition of the main section is to satisfy the wave vector superposition, i.e.
Decomposing the vector to obtain a formula:
since the three pulses satisfy the vector superposition, so
The square addition operation is performed on the formulas (2) and (3), and the result is substituted into
To obtain:
according to the frequency multiplication, omega can be known3＝2ω1＝2ω2Therefore, it is
Wherein n isp(θ3) The angle between the frequency-doubled light and the main axis is theta3The refractive index of the time-doubled light; n isp(θ1)、np(θ2) For the pulse to be measured and the main shaft to form an included angle theta1、θ2The refractive index of the pulse to be measured;
the optical refractive index of the frequency doubling light e is calculated according to a formula:
wherein n isp(θ3) Is that
Combining the formulas (5) and (6), and drawing by matlab to obtain theta3With theta1The change map of (1), theta corresponding to the vertex of the image3The value is the desired cutting angle.
8. An ultrashort pulse reconstruction method based on a BP neural network is characterized by comprising the following steps:
(1) acquiring a scanning spectrum trace diagram of the ultrashort pulse to be detected;
(2) and inputting the scanning spectrum trace graph into a pre-trained BP neural network model, and outputting the amplitude and the phase of the ultrashort pulse to be detected.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111383297.0A CN114216575B (en) | 2021-11-22 | 2021-11-22 | Ultrashort pulse reconstruction system and method based on BP neural network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111383297.0A CN114216575B (en) | 2021-11-22 | 2021-11-22 | Ultrashort pulse reconstruction system and method based on BP neural network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114216575A true CN114216575A (en) | 2022-03-22 |
| CN114216575B CN114216575B (en) | 2024-04-05 |
Family
ID=80697707
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111383297.0A Active CN114216575B (en) | 2021-11-22 | 2021-11-22 | Ultrashort pulse reconstruction system and method based on BP neural network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114216575B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5530544A (en) * | 1992-10-26 | 1996-06-25 | Sandia Corporation | Method and apparatus for measuring the intensity and phase of one or more ultrashort light pulses and for measuring optical properties of materials |
| US20020057435A1 (en) * | 2000-11-14 | 2002-05-16 | Georgia Tech Research Corporation | Electromagnetic wave analyzer |
| CN101526465A (en) * | 2009-04-22 | 2009-09-09 | 天津大学 | Quick multi-wavelength tissue optical parameter measuring device and trans-construction method |
| US20140237014A1 (en) * | 2013-02-20 | 2014-08-21 | Daniel James Kane | Method and apparatus for determining wave characteristics using interaction with a known wave |
| US20190170588A1 (en) * | 2016-08-10 | 2019-06-06 | Sphere Ultrafast Photonics, S.A. | Ultrashort laser pulse characterization and compression method |
| WO2021123481A1 (en) * | 2019-12-19 | 2021-06-24 | Universidad De Salamanca | Method and system for the temporal and spectral characterization of the amplitude and phase of ultrashort laser pulses |
| CN113063507A (en) * | 2021-03-26 | 2021-07-02 | 中国科学院物理研究所 | Prediction method of ultrashort pulse width based on convolutional neural network |
-
2021
- 2021-11-22 CN CN202111383297.0A patent/CN114216575B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5530544A (en) * | 1992-10-26 | 1996-06-25 | Sandia Corporation | Method and apparatus for measuring the intensity and phase of one or more ultrashort light pulses and for measuring optical properties of materials |
| US20020057435A1 (en) * | 2000-11-14 | 2002-05-16 | Georgia Tech Research Corporation | Electromagnetic wave analyzer |
| CN101526465A (en) * | 2009-04-22 | 2009-09-09 | 天津大学 | Quick multi-wavelength tissue optical parameter measuring device and trans-construction method |
| US20140237014A1 (en) * | 2013-02-20 | 2014-08-21 | Daniel James Kane | Method and apparatus for determining wave characteristics using interaction with a known wave |
| US20190170588A1 (en) * | 2016-08-10 | 2019-06-06 | Sphere Ultrafast Photonics, S.A. | Ultrashort laser pulse characterization and compression method |
| WO2021123481A1 (en) * | 2019-12-19 | 2021-06-24 | Universidad De Salamanca | Method and system for the temporal and spectral characterization of the amplitude and phase of ultrashort laser pulses |
| CN113063507A (en) * | 2021-03-26 | 2021-07-02 | 中国科学院物理研究所 | Prediction method of ultrashort pulse width based on convolutional neural network |
Non-Patent Citations (1)
| Title |
|---|
| TOM ZAHAVY.ET.AL: "Deep learning reconstruction of ultrashort pulses", OPTICA, vol. 5, no. 5, 18 May 2018 (2018-05-18) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114216575B (en) | 2024-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6504612B2 (en) | Electromagnetic wave analyzer | |
| US7387027B2 (en) | Characterization of materials with optically shaped acoustic waveforms | |
| US6219142B1 (en) | Method and apparatus for determining wave characteristics from wave phenomena | |
| JP4790560B2 (en) | Single terahertz wave time waveform measurement device | |
| CN101608999A (en) | Real-time observation single-beam dual-mode parameter adjustable Z scanning device and measurement method | |
| JP2017181259A (en) | Pulse light waveform measurement method and waveform measurement device | |
| CN101732035B (en) | Method for optical super resolution based on optical path encoding and coherent synthesis | |
| CN102230827A (en) | Mach-Zehnder point diffraction interferometer and method for reconstructing laser complex amplitudes | |
| CN109883350A (en) | A high-precision measurement system and measurement method for the internal topography of a special-shaped curved surface structure | |
| Shi et al. | Imaging consecutive targets through scattering medium and around corners beyond the optical memory effect using untrained network | |
| CN107655832A (en) | Based on the femtosecond Two-dimensional spectroscopy acquisition method changed in chirp and system | |
| CN114216575A (en) | Ultra-short pulse reconstruction system and method based on BP neural network | |
| CN109186769A (en) | A method of the measurement ellipse inclined rate of chirped pulse | |
| CN108344711A (en) | A kind of method and system improving terahertz pulse imaging resolution | |
| Xie et al. | Deep learning for estimation of Kirkpatrick–Baez mirror alignment errors | |
| CN118294016A (en) | Ultrafast single-pixel imaging enhancement system based on spatial optical filter | |
| Trofimov et al. | Highly effective method for temporal terahertz spectroscopy under the condition of random probe signals | |
| Grachev et al. | An open-source 3d-printed terahertz pulse time-domain holographic detection module for broadband beam inspection | |
| CN102636830B (en) | Strip phase aperture, 4f phase coherent nonlinear imaging system using the aperture and method for measuring nonlinear refractive index by the system | |
| US20230333007A1 (en) | Dispersion measurement device and dispersion measurement method | |
| Kale | Towards High Fidelity Quantum Computation and Simulation with Rydberg Atoms | |
| CN113514149A (en) | Optical Pulse Signal Processing System | |
| JP2021047130A (en) | Spectrometry device and method | |
| Brinks | Time-Resolved X-Ray Diffraction from Coherent Acoustic Phonons Analyzed by Combining Microscopic Modeling with Deep Learning Based Strain Retrieval | |
| CN220399276U (en) | Single-shot terahertz time-domain spectrometer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |