CN105651399B - The weak phase measurement device of time domain phase recovery full optical fiber laser pulse and measurement method - Google Patents
The weak phase measurement device of time domain phase recovery full optical fiber laser pulse and measurement method Download PDFInfo
- Publication number
- CN105651399B CN105651399B CN201610028749.6A CN201610028749A CN105651399B CN 105651399 B CN105651399 B CN 105651399B CN 201610028749 A CN201610028749 A CN 201610028749A CN 105651399 B CN105651399 B CN 105651399B
- Authority
- CN
- China
- Prior art keywords
- phase
- laser pulse
- optical fiber
- complex amplitude
- modulator
- 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.)
- Active
Links
- 239000013307 optical fiber Substances 0.000 title claims abstract description 43
- 238000005259 measurement Methods 0.000 title claims abstract description 21
- 238000011084 recovery Methods 0.000 title claims abstract description 16
- 238000000691 measurement method Methods 0.000 title claims abstract description 8
- 239000000835 fiber Substances 0.000 claims abstract description 62
- 238000012360 testing method Methods 0.000 claims abstract description 44
- 238000000034 method Methods 0.000 claims description 23
- 238000009826 distribution Methods 0.000 claims description 16
- 230000003321 amplification Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 7
- 230000003595 spectral effect Effects 0.000 claims description 6
- 239000006185 dispersion Substances 0.000 claims description 5
- 210000001367 artery Anatomy 0.000 claims description 4
- 230000001934 delay Effects 0.000 claims description 4
- 210000003462 vein Anatomy 0.000 claims description 4
- 238000013101 initial test Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 abstract description 3
- 230000005540 biological transmission Effects 0.000 abstract description 2
- 230000010363 phase shift Effects 0.000 description 6
- 238000003754 machining Methods 0.000 description 5
- 230000004927 fusion Effects 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 241000239290 Araneae Species 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000013078 crystal Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 239000010437 gem Substances 0.000 description 1
- 229910001751 gemstone Inorganic materials 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000003698 laser cutting Methods 0.000 description 1
- 238000005542 laser surface treatment Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000008188 pellet Substances 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000004080 punching Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
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
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
Abstract
A kind of weak phase measurement device of time domain phase recovery full optical fiber laser pulse and measurement method, it is fiber optic splitter that device composition, which includes: along testing laser pulse input direction, testing laser pulse is divided into strong, weak two-beam by the fiber optic splitter, it is successively adjustable optic fibre delayer, high speed fibre phase-modulator, dispersive optical fiber and oscillograph along strong beam direction, it is successively high speed PIN photoelectric tube, arbitrary waveform generator, electric signal amplifier, the modulation input terminal of the output termination of the electric signal amplifier high speed fibre phase-modulator along weak beam direction.The present invention carries out phase-modulation and dispersive transmission to laser pulse using all optical fibre structure, has the characteristics that compact-sized, simple and flexible, and it is different from other laser pulse phase measurement modes, the present invention can measure with weak phase picosecond or nanosecond laser pulses, be suitable in the case of high repetition frequency and low-repetition-frequency.
Description
Technical field
The present invention relates to weak phase laser pulse, especially a kind of weak phase of time domain phase recovery full optical fiber laser pulse is surveyed
Measure device and measurement method.The present invention carries out phase-modulation to laser pulse using fibre-optic waveguide phase-modulator, then when passing through
The method of domain phase recovery obtains the waveform and phase of laser pulse.The present invention is suitable for the skin with weak phase, weak chirp
Second, perhaps nanosecond laser pulses can work in high repetition frequency or low-repetition-frequency.The present apparatus uses all -fiber
Structure can increase device stability and compactedness, to obtain reliable and stable laser pulse phase distribution.
Background technique
Since First laser manufactures so far, the application of laser and laser has gradually penetrated into each row of society
In each industry, especially for precision machining industry, the application of laser greatly improves machining accuracy.Wherein high-energy nanosecond or
Ps Laser Pulse studies (inertial confinement fusion), laser precision machining, laser cutting, laser thunder in laser physics
It reaches, the multiple fields such as Ultrafast spectrum, medicine, high-energy physics are widely used.For example, high-energy nanosecond laser pulses can be applied
In inertial confinement fusion, so that pellet reaches nuclear fusion condition, a large amount of energy is released, laser inertial confinement core is poly-
Change is expected to realize controllable nuclear fusion in future, fundamentally solves energy problem.High-energy picosecond pulse can be applied to laser
Processing and Laser Surface Treatment, fast laser pulse and being greatly improved without heat deposition and contactless characteristic during material effect
The processing characteristics such as controllability, machining accuracy and the surface smoothness of processing pattern, in metal, crystal, jewel, glass, high score
Excellent characteristic is shown in the processing cutting of the multiple materials such as sub- polymer even explosive, in auto industry, medical instrument, work
The Precision Machinings fields such as industry safety present wide application prospect.High-energy ps pulsed laser and ns pulsed laser and picosecond pulse laser at present
It is all made of master oscillator+casacade multi-amplifier structure.But in amplification process, due to nonlinear effect such as Self-phase modulation
Presence, laser pulse can accumulate nonlinear phase shift, to have an impact to rear class amplification.But unlike that the feelings of ultrashort pulse
The nonlinear phase shift that condition, nanosecond or picosecond pulse are accumulated in transmission process is typically small, and from spectrally, spectrum widening is not
Obviously, but these small nonlinear phase shifts can generate large effect in subsequent amplification process.Further, since there are some
The compensation way of nonlinear phase shift, such as the precompensation method of Direct Phase modulation, can be in its known nonlinear phase shift size
In the case where, nonlinear phase shift is compensated using precompensation means.This is just needed to these laser arteries and veins with weak phase
The phase of punching carries out precise measurement.There are many kinds of the modes of burst length phase (spectrum phase measurement) measurement at present, including
The modes such as FROG (frequency resolution optical switch), SPIDER (self-reference spectral interference) and time auto-correlation, but these sides
Formula is all based on auto-correlation/cross-correlation principle, carries out related operation or related to reference pulse progress by pulse itself
Operation.FROG is using being that non-linear process carries out related operation, and SPIDER scheduling algorithm is to carry out phase by way of interference
Close operation.But the structure of these modes is generally all more complicated, while related operation has to the intensity of measured pulse
Seek (FROG) perhaps lower (SPIDER etc.) of resolution ratio and these methods are for weak phase nanosecond or Ps Laser Pulse
The case where, it is often helpless.
Therefore the present invention propose using time domain phase recovery method combination Direct Phase modulate device come realize for
The high-precision phase measurement of weak phase nanosecond or Ps Laser Pulse.This method, which is suitable for measurement, has the narrow of weak phase distribution
High-precision phase measurement may be implemented in spectral laser pulse, the data calculation process of Simultaneous Iteration process.
Summary of the invention
It is an object of the invention to overcome the time phase measurement method of above-mentioned existing laser pulse not to be capable of measuring weak phase
The shortcomings that nanosecond or Ps Laser Pulse, propose a kind of weak phase measurement device of time domain phase recovery full optical fiber laser pulse and
Measurement method, this method carry out laser pulse by the way of time domain phase Iterative restoration, using Direct Phase modulating device
Phase-modulation, by the time waveform by the laser pulse by phase-modulation after dispersive medium as known conditions meter
Calculation obtains the waveform and phase distribution of primary laser pulse.Using the structure of all -fiber can be improved device flexibility with it is compact
Property, realize the high-acruracy survey of the time phase distribution to weak phase laser pulse.
Technical solution of the invention is as follows:
The weak phase measurement device of a kind of time domain phase recovery full optical fiber laser pulse, it is characterised in that its composition includes: edge
Testing laser pulse input direction is fiber optic splitter, which is divided into strong, weak two-beam, edge for testing laser pulse
Strong beam direction is successively adjustable optic fibre delayer, high speed fibre phase-modulator, dispersive optical fiber and oscillograph, along dim light Shu Fang
To being successively high speed PIN photoelectric tube, arbitrary waveform generator, electric signal amplifier, the output of the electric signal amplifier terminates institute
The modulation input terminal for the high speed fibre phase-modulator stated, the delay adjustment precision of the adjustable optic fibre delayer are 1ps, institute
The length for stating dispersive optical fiber meets the following conditions:
Wherein, β2For the second-order dispersion of the dispersive optical fiber, L is the length of dispersive optical fiber, and Δ ν is by high speed fibre
The spectral width of testing laser pulse after the phase-modulation of phase-modulator；
The modulation signal that the arbitrary waveform generator generates is single order Gaussian pulse, and the single order Gaussian pulse
Pulse width τ meets τ≤Δ T, and Δ T is the pulsewidth of testing laser pulse.
Using the weak phase measurement device of above-mentioned time domain phase recovery full optical fiber laser pulse to the weak phase of laser pulse to be measured
Measurement method, it is characterised in that method includes the following steps:
1. setting t0The modulation through the electric signal amplifier amplification output generated for the arbitrary waveform generator
Reach the height relative to testing laser pulse at the time of the center of electric signal reaches the high speed fibre phase-modulator
Initial relative time delay at the time of fast fibre optic phase modulator adjusts the delay time of the adjustable optic fibre delayer every time
It is Δ t that delay, which increases, every primary delay of adjustment, described one testing laser pulse strength I of oscillograph recordingm, successively obtain
I1、I2、┄、Im、┄I2n+1, the delay that the adjustable optic fibre delayer generates after the m times adjusting is t0+ m Δ t, the m times adjusting
Afterwards, testing laser pulse collects this by the oscillograph after the high speed fibre phase-modulator and dispersive optical fiber
Delay time (t0Laser pulse intensity I under+m Δ t)mAre as follows:
Im=| Am|2 (2)
Wherein, m=1,2,3 ... 2n+1, n are any positive integer and meet 2n+1 >=Δ T/ Δ t, AmFor testing laser pulse
By the light field complex amplitude after dispersive optical fiber (7)；
2. using time domain Phase Retrieve Algorithm to the laser pulse intensity ImData processing is carried out, testing laser is calculated
The phase distribution of pulse, the specific steps are as follows:
1) data initialization is arranged: i is current iteration number, and m is the serial number of the m times adjusting adjustable optic fibre delayer, enables i
The maximum value of=0, m are 2n+1；N0,2n+1(t) COMPLEX AMPLITUDE of the testing laser pulse generated for computer random, β are most
Small calculating error, K are maximum number of iterations；
2) i=i+1, m=0, the initial testing laser pulse complex amplitude N that current iteration calculates are enabledi,mIt (t) is i-1 iteration
Middle m=2n+1 corresponds to calculated light field complex amplitude, i.e. Ni,m(t)=Ni-1,2n+1(t)；
3) m=m+1 is enabled, currently to testing laser pulse complex amplitude E initial in the iterative calculation of mi,mIt (t) is m-1 iteration
The middle calculated light field complex amplitude of correspondence, i.e. Ei,m(t)=Ni,m-1(t)；
4) the light field complex amplitude after high speed fibre phase-modulator is calculated as follows
Wherein, V is the voltage magnitude of modulation signal, VπFor the half-wave voltage of high speed fibre phase-modulator,
(t-m Δ t) is the modulation signal of the time delays with m Δ t to B；
5) it is calculated as follows againLight field complex amplitude A after dispersive optical fiberi,m(t) and light field phase distributionIt is respectively as follows:
Wherein: F is Fourier transformation, F-1For inverse Fourier transform, ω is light field angular frequency；
6) laser pulse intensity by high speed fibre phase-modulator and dispersive optical fiber measured using oscillograph the m times
Im, instead of the light field complex amplitude that (4) formula is calculated, and retain phase invariant, obtain updated complex amplitude
7) by the light field complex amplitudeThe input terminal for inversely traveling to dispersive optical fiber, obtains updated incident light
Field complex amplitude
8) the complex amplitude N of high speed fibre phase-modulator input terminal testing laser pulse is calculated according to the following formulai,m(t):
Wherein, | φ (t-m Δ t) |maxThe phase tune generated on high speed fibre phase-modulator is loaded for modulation signal
System, conj (*) are function complex conjugate, and α is to prevent except null divisor；
9) as m < 2n+1, return step 3)；As m=2n+1, the error that current i-th iteration calculates is calculated as follows
Error are as follows:
If 10) Error < β, stop iterating to calculate, carry out next step 11)；If Error>=β, and i<K, then return
12) step 2) enters step if i=K；
11)Ni,2n+1(t) be testing laser pulse light field complex amplitude, phase distribution is according to light field complex amplitude Ni,2n+1
(t) obtaining the phase distribution of testing laser pulse in time using conventional phase unwrapping algorithm is Ni,2n+1(t)/|
Ni,2n+1(t)|；
12) as i=K, show that the maximum value 2n+1 of current m is not able to satisfy the requirement of computational accuracy, then enable n=n+1,
Increase laser pulse intensity ImMeasurement number, return step 1), continue to calculate.
The invention has the advantages that:
1. apparatus of the present invention use the structure of all-fiber, apparatus of the present invention are compact-sized, are conveniently adjusted.
2. the method restored using time phase, can measure the phase of weak phase laser pulse.
3. can flexibly control modulation signal using the method that Direct Phase is modulated.
4. can work in Gao Zhongying and low repetition.
Detailed description of the invention
Fig. 1 is the structural block diagram of the weak phase measurement device of time domain phase recovery full optical fiber laser pulse of the present invention.
Fig. 2 is the time phase restoration methods flow chart of the weak phase laser pulse of the present invention.
Specific embodiment
Below with reference to embodiment and attached drawing, the present invention will be further described, but protection model of the invention should not be limited with this
It encloses.
First referring to Fig. 1, Fig. 1 is the structure of the weak phase measurement device of time domain phase recovery full optical fiber laser pulse of the present invention
Block diagram.As seen from the figure, the composition of the weak phase measurement device of time domain phase recovery full optical fiber laser pulse of the present invention includes: along to be measured
Laser pulse In input direction is fiber optic splitter 1, which divides testing laser pulse In for strong, weak two-beam,
It is successively adjustable optic fibre delayer 5, high speed fibre phase-modulator 6, dispersive optical fiber 7 and oscillograph 8 along strong beam direction, along weak
Beam direction is successively high speed PIN photoelectric tube 2, arbitrary waveform generator 3, electric signal amplifier 4, the electric signal amplifier 4
The modulation input terminal of the output termination high speed fibre phase-modulator 6, the delay adjustment of the adjustable optic fibre delayer 5
Precision is 1ps, and the length of the dispersive optical fiber 7 meets the following conditions:
Wherein, β2For the second-order dispersion of dispersive optical fiber 7, L is the length of dispersive optical fiber 7, and Δ ν is by high speed fibre phase
The spectral width of testing laser pulse after 6 phase-modulation of modulator；
The modulation signal that the arbitrary waveform generator 3 generates is single order Gaussian pulse, should and single order Gaussian pulse
Pulse width meet τ≤Δ T, wherein τ be modulation signal pulse width, Δ T be testing laser pulse pulsewidth.
It is two parts that incident laser pulse, which first passes around 1 beam splitting of fiber optic splitter, and wherein splitting ratio is 10% and 90%.
10% end laser pulse, which is converted into electric signal as the triggering of arbitrary waveform generator 3 by high-speed photodetector PIN pipe 2, to be believed
Number.The arbitrary waveform generator 3 being triggered exports the electric arteries and veins that a pulse width is less than or equal to testing laser pulse width
Signal is rushed, the electric impulse signal is after the amplification of high-gain high-speed electric amplifier 4 as the modulation electricity of fibre optic phase modulator 6
Signal.The effect of optical fiber adjustable light delay 5 is then passed through in 90% end testing laser pulse, changes it between electrical modulation signal
Relative time delay.Laser pulse after time delays enters fibre optic phase modulator 6, obtains phase tune in time
System, laser pulse light spectrum broaden.Laser pulse after phase-modulation then enters the dispersion having compared with large dispersion
Optical fiber 7 is broadened by the laser pulse of phase-modulation in time, is detected the waveform of laser pulse by oscillograph 8, and its
Time waveform represents the spectral signature of laser pulse.
It can use time phase recovery by the collected Laser pulse time Wave data of apparatus of the present invention oscillograph 8
Method obtains the phase distribution of laser pulse, and referring to fig. 2, detailed process is as follows for the flow chart of time phase restoration methods:
Assuming that t0The tune through the amplification output of electric signal amplifier 4 generated for the arbitrary waveform generator 3
Reach the high speed relative to testing laser pulse at the time of the center of electric signal processed reaches on high speed fibre phase-modulator 6
Initial relative time delay at the time of fibre optic phase modulator 6 adjusts the delay time of the adjustable optic fibre delayer 5 every time
It is Δ t that delay, which increases, then the m times delay for adjusting the adjustable optic fibre delayer 5 is t0+ m Δ t, the m times testing laser pulse
The delay time is collected by the oscillograph 8 after the high speed fibre phase-modulator 6 and dispersive optical fiber 7
(t0Laser pulse intensity I under+m Δ t)mAre as follows:
Im=| Am|2 (9)
Wherein, m=1,2,3 ... 2n+1, n are any positive integer, AmIt is testing laser pulse by after dispersive optical fiber 7
Light field complex amplitude.
1) data initialization is carried out first, and enabling i is current iteration number, and m is the m times adjusting adjustable optic fibre delayer 5, i
=0, m=2n+1；Enable N0,2n+1It (t) is the complex amplitude of the testing laser pulse in current iteration calculating, N0,2n+1It (t) is random raw
At testing laser pulse COMPLEX AMPLITUDE；
2) i=i+1, m=0 are enabled, while initial testing laser pulse complex amplitude N in current iteration calculatingi,m(t) it is i-1 times
M=2n+1 corresponds to calculated light field complex amplitude, i.e. N in iterationi,m(t)=Ni-1,2n+1(t)；
3) m=m+1 is enabled, while currently to testing laser pulse complex amplitude E initial in the iterative calculation of mi,m(t) it is m-1 times
Calculated light field complex amplitude, i.e. E are corresponded in iterationi,m(t)=Ni,m-1(t)；
4) it calculates by the light field complex amplitude after high speed fibre phase-modulator 6Are as follows:
Wherein V is the voltage magnitude of modulation signal, VπFor the half-wave voltage of high speed fibre phase-modulator (6),
(t-m Δ t) is the modulation signal of the time delays with m Δ t to B；
5) it calculates againLight field complex amplitude A after dispersive optical fiber 7i,m(t) and light field phase distributionRespectively
Are as follows:
Wherein: F is Fourier transformation, F-1For inverse Fourier transform, ω is light field angular frequency；
6) it is strong that the obtained laser pulse by high speed fibre phase-modulator 6 and dispersive optical fiber 7 is measured according to oscillograph 8
Spend Im, instead of the light field complex amplitude that (11) formula is calculated, and retain phase invariant, obtain updated complex amplitudeIts mistake
Journey is as follows:
7) by updated light field complex amplitudeThe input terminal for inversely traveling to dispersive optical fiber 7, obtains updated incidence
Light field complex amplitude
8) the complex amplitude N of 6 input terminal testing laser pulse of high speed fibre phase-modulator is calculated according to the following formulai,m(t):
Wherein: | φ (t-m Δ t) |maxThe phase generated on high speed fibre phase-modulator 6 is loaded for modulation signal
Modulation, conj (*) are function complex conjugate, and α is to prevent except null divisor；
9) as m < 2n+1, return step 3)；Otherwise as m=2n+1, the error that current i-th iteration calculates is calculated
Error are as follows:
If 10) Error < β, stop iterating to calculate, carry out next step 11)；If otherwise Error>=β, and i<K, then
Return step 2), wherein β is minimum of computation error, and K is maximum number of iterations；If i=K, step (12) are carried out；
11) when iterative calculation terminates, the light field complex amplitude of testing laser pulse is when last iterative calculation terminates
Ni,2n+1(t), phase distribution is according to light field complex amplitude Ni,2n+1(t) testing laser is obtained using conventional phase unwrapping algorithm
The phase distribution of pulse in time is Ni,2n+1(t)/|Ni,2n+1(t)|；
12) as i=K, show that the maximum value 2n+1 of current m is not able to satisfy the requirement of computational accuracy, then enable n=n+1
Increase laser pulse intensity ImMeasurement number, return step 1), continue to calculate.
Practice have shown that the present invention uses the structure of all-fiber, apparatus of the present invention are compact-sized, are conveniently adjusted.Using when
Between phase recovery method, the phase of weak phase laser pulse can be measured.It can using the method that Direct Phase is modulated
Flexibly to control modulation signal.It can work in Gao Zhongying and low repetition.
Claims (2)
1. a kind of weak phase measurement device of time domain phase recovery full optical fiber laser pulse, it is characterised in that its constitute include: along to
Surveying laser pulse input direction is fiber optic splitter (1), which is divided into strong, weak two beam for testing laser pulse
Light is successively adjustable optic fibre delayer (5), high speed fibre phase-modulator (6), dispersive optical fiber (7) along strong beam direction and is shown
Wave device (8) is successively high speed PIN photoelectric tube (2), arbitrary waveform generator (3), electric signal amplifier (4) along weak beam direction,
The modulation input terminal of the output termination of the electric signal amplifier (4) the high speed fibre phase-modulator (6), described is adjustable
The delay adjustment precision of fiber delay time device (5) is 1ps, and the length of the dispersive optical fiber (7) meets the following conditions:
Wherein, β2For the second-order dispersion of the dispersive optical fiber (7), L is the length of dispersive optical fiber (7), and Δ ν is by high-speed light
The spectral width of testing laser pulse after the phase-modulation of fine phase-modulator (6)；
The modulation signal that the arbitrary waveform generator (3) generates is single order Gaussian pulse, and the single order Gaussian pulse
Pulse width τ meets τ≤Δ T, and Δ T is the pulsewidth of testing laser pulse.
2. using the weak phase measurement device of time domain phase recovery full optical fiber laser pulse described in claim 1 to testing laser arteries and veins
Rush the measurement method of weak phase, it is characterised in that method includes the following steps:
1. setting t0The modulation through electric signal amplifier (4) the amplification output generated for the arbitrary waveform generator (3)
Reach described relative to testing laser pulse at the time of the center of electric signal reaches high speed fibre phase-modulator (6)
Initial relative time delay at the time of high speed fibre phase-modulator (6) adjusts prolonging for adjustable optic fibre delayer (5) every time
When the time delay to increase be Δ t, every primary delay of adjustment, the oscillograph (8) records a testing laser pulse strength
Im, successively obtain I1、I2、┄、Im、┄I2n+1, the delay that the adjustable optic fibre delayer (5) generates after the m times adjusting is t0+
M Δ t, after the m times is adjusted, testing laser pulse quilt after the high speed fibre phase-modulator (6) and dispersive optical fiber (7)
The oscillograph (8) collects the delay time (t0Laser pulse intensity I under+m Δ t)mAre as follows:
Im=| Am|2 (2)
Wherein, m=1,2,3 ... 2n+1, n are any positive integer and meet 2n+1 >=Δ T/ Δ t, AmFor testing laser pulse process
Light field complex amplitude after dispersive optical fiber (7)；
2. using time domain Phase Retrieve Algorithm to the laser pulse intensity ImData processing is carried out, testing laser pulse is calculated
Phase distribution, the specific steps are as follows:
1) data initialization is arranged: i is current iteration number, and m is the m times serial number for adjusting adjustable optic fibre delayer (5), enables i
The maximum value of=0, m are 2n+1；N0,2n+1(t) COMPLEX AMPLITUDE of the testing laser pulse generated for computer random, β are most
Small calculating error, K are maximum number of iterations；
2) i=i+1, m=0, the initial testing laser pulse complex amplitude N that current iteration calculates are enabledi,mIt (t) is m in i-1 iteration
=2n+1 corresponds to calculated light field complex amplitude, i.e. Ni,m(t)=Ni-1,2n+1(t)；
3) m=m+1 is enabled, currently to testing laser pulse complex amplitude E initial in the iterative calculation of mi,mIt (t) is right in m-1 iteration
Answer calculated light field complex amplitude, i.e. Ei,m(t)=Ni,m-1(t)；
4) the light field complex amplitude after high speed fibre phase-modulator (6) is calculated as follows
Wherein, V is the voltage magnitude of modulation signal, VπFor the half-wave voltage of high speed fibre phase-modulator (6), B (t-m Δ t)
For the modulation signal of the time delays with m Δ t, j=(- 1)1/2For imaginary unit；
5) it is calculated as follows againLight field complex amplitude A after dispersive optical fiber (7)i,m(t) and light field phase distribution
It is respectively as follows:
Wherein: F is Fourier transformation, F-1For inverse Fourier transform, ω is light field angular frequency；
6) the laser arteries and veins by high speed fibre phase-modulator (6) and dispersive optical fiber (7) measured using oscillograph (8) the m times
Rush intensity Im, instead of the light field complex amplitude that (4) formula is calculated, and retain phase invariant, obtain updated complex amplitude
7) by the light field complex amplitudeThe input terminal for inversely traveling to dispersive optical fiber (7), obtains updated incident field
Complex amplitude
8) the complex amplitude N of high speed fibre phase-modulator (6) input terminal testing laser pulse is calculated according to the following formulai,m(t):
Wherein, | φ (t-m Δ t) | the phase-modulation generated on high speed fibre phase-modulator (6) is loaded for modulation signal,
Conj (*) is function complex conjugate, and α is to prevent except null divisor；
9) as m < 2n+1, return step 3)；As m=2n+1, the error that current i-th iteration calculates is calculated as follows
Error are as follows:
If 10) Error < β, stop iterating to calculate, carry out next step 11)；If Error>=β, and i<K, then return step
2) it, if i=K, enters step 12)；
11)Ni,2n+1(t) be testing laser pulse light field complex amplitude, phase distribution is according to light field complex amplitude Ni,2n+1(t)
Obtaining the phase distribution of testing laser pulse in time using conventional phase unwrapping algorithm is Ni,2n+1(t)/|Ni,2n+1
(t)|；
12) as i=K, show that the maximum value 2n+1 of current m is not able to satisfy the requirement of computational accuracy, then enable n=n+1, increase
Laser pulse intensity ImMeasurement number, return step 1), continue to calculate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610028749.6A CN105651399B (en) | 2016-01-15 | 2016-01-15 | The weak phase measurement device of time domain phase recovery full optical fiber laser pulse and measurement method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610028749.6A CN105651399B (en) | 2016-01-15 | 2016-01-15 | The weak phase measurement device of time domain phase recovery full optical fiber laser pulse and measurement method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105651399A CN105651399A (en) | 2016-06-08 |
| CN105651399B true CN105651399B (en) | 2019-04-19 |
Family
ID=56487626
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610028749.6A Active CN105651399B (en) | 2016-01-15 | 2016-01-15 | The weak phase measurement device of time domain phase recovery full optical fiber laser pulse and measurement method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105651399B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107490434B (en) * | 2017-07-24 | 2020-01-17 | 南开大学 | Method and device for high-speed measurement of multimode optical pulse cluster space-time spectrum information |
| CN108956432B (en) * | 2018-08-10 | 2020-08-07 | 武汉大学 | Flow type high-speed super-resolution imaging device and method based on structured light |
| CN109100304B (en) * | 2018-08-10 | 2020-12-01 | 武汉大学 | Single-pixel high-speed super-resolution imaging device and method based on time domain stretching |
| CN110864806A (en) * | 2019-11-15 | 2020-03-06 | 天津大学 | Real-time spectral measurement system based on grating-centered time domain spread dispersion Fourier transform |
| CN111340902B (en) * | 2019-12-11 | 2023-08-18 | 华中科技大学苏州脑空间信息研究院 | Optical phase modulation method and spatial light modulation method for irradiation of arbitrary position and shape |
| CN110987200B (en) * | 2019-12-17 | 2021-04-20 | 华中科技大学 | Method and device for measuring attosecond pulse |
| CN112510477A (en) * | 2020-12-08 | 2021-03-16 | 中国工程物理研究院激光聚变研究中心 | Device and method for improving stability of optical fiber laser pulse amplitude |
| CN113091922B (en) * | 2021-03-16 | 2023-01-20 | 中国科学院上海光学精密机械研究所 | Rapid convergence ultrashort pulse waveform and phase reconstruction method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4990762A (en) * | 1981-12-21 | 1991-02-05 | The United States Of America As Represented By The Secretary Of The Navy | Laser beam phase measurement and control system |
| CN103575408A (en) * | 2013-10-29 | 2014-02-12 | 桂林电子科技大学 | Coherent light detection system and method for optical fiber delay line phase control |
| CN103776550A (en) * | 2014-01-07 | 2014-05-07 | 浙江大学 | Supercontinuum pulse laser measurement device based on nonlinear nano material and method |
-
2016
- 2016-01-15 CN CN201610028749.6A patent/CN105651399B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4990762A (en) * | 1981-12-21 | 1991-02-05 | The United States Of America As Represented By The Secretary Of The Navy | Laser beam phase measurement and control system |
| CN103575408A (en) * | 2013-10-29 | 2014-02-12 | 桂林电子科技大学 | Coherent light detection system and method for optical fiber delay line phase control |
| CN103776550A (en) * | 2014-01-07 | 2014-05-07 | 浙江大学 | Supercontinuum pulse laser measurement device based on nonlinear nano material and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105651399A (en) | 2016-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105651399B (en) | The weak phase measurement device of time domain phase recovery full optical fiber laser pulse and measurement method | |
| US11143558B2 (en) | Method and system for measuring transient time width of ultrashort pulse | |
| McDonald | Negative group velocity | |
| US20160323036A1 (en) | System and method for removing far field limitations in microwave photonic arbitrary waveform generators | |
| CN103389164A (en) | Extreme ultraviolet attosecond pulse width measurement method and device | |
| Kou et al. | All-phase FFT based distance measurement in laser self-mixing interferometry | |
| CN103148940A (en) | Light asynchronous sampling signal measurement method and system | |
| JP6539930B2 (en) | Brillouin frequency shift distribution measurement system and Brillouin frequency shift distribution measurement method | |
| Baudin et al. | Fiber‐based measurement of temporal intensity and phase profiles of an optical telecommunication pulse through self‐phase modulation | |
| Rhodes et al. | Validating data analysis of broadband laser ranging | |
| Zolotovskii et al. | Modulation instability of pulsed radiation in an optical waveguide in the presence of the traveling refractive index wave | |
| EP0250917B1 (en) | Method for signal transmission and optical communications | |
| Wang et al. | Real‐Time and High‐Accuracy Microwave Frequency Identification Based on Ultra‐Wideband Optical Chirp Chain Transient SBS Effect | |
| CN104913853B (en) | Method and system for measuring ultra-short intense laser prepulsing | |
| CN108845330B (en) | Laser ranging method based on loop frequency measurement | |
| CN102937481B (en) | Femtosecond laser pulse complete reconstruction system and method | |
| Dorrer et al. | Single-shot high-resolution characterization of optical pulses by spectral phase diversity | |
| Kikuchi | Theory of sonogram characterization of optical pulses | |
| Fontaine et al. | Frequency-to-time assisted interferometry for polarization-diversified, single-shot, full-field waveform measurement | |
| CN109116275B (en) | Method for accurately measuring consistency of two magnetic fields at small distance | |
| CN108519567B (en) | Pulse magnetic field measurement method and corresponding device | |
| CN111693157B (en) | Time-frequency multiplexing-based large pump bandwidth ultrafast pulse time domain measurement method and system | |
| JP2022030729A (en) | Terahertz waveform detection device and terahertz waveform detection method | |
| Kurucz et al. | Dispersive Fourier transform based single-shot CEP drift measurement at arbitrary repetition rate | |
| Voskoboinik et al. | Fast and distributed dynamic sensing of strain using Sweep-Free Brillouin Optical Time-Domain analysis (SF-BOTDA) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |