CN106663913B - System for adjusting the phase of a large amount of laser sources - Google Patents
System for adjusting the phase of a large amount of laser sources Download PDFInfo
- Publication number
- CN106663913B CN106663913B CN201580041454.6A CN201580041454A CN106663913B CN 106663913 B CN106663913 B CN 106663913B CN 201580041454 A CN201580041454 A CN 201580041454A CN 106663913 B CN106663913 B CN 106663913B
- Authority
- CN
- China
- Prior art keywords
- phase
- plane
- light beam
- laser source
- adjusting
- 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
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B6/00—Light guides; Structural details of arrangements comprising light guides and other optical elements, e.g. couplings
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/05—Construction or shape of optical resonators; Accommodation of active medium therein; Shape of active medium
- H01S3/08—Construction or shape of optical resonators or components thereof
- H01S3/08004—Construction or shape of optical resonators or components thereof incorporating a dispersive element, e.g. a prism for wavelength selection
- H01S3/08009—Construction or shape of optical resonators or components thereof incorporating a dispersive element, e.g. a prism for wavelength selection using a diffraction grating
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/10—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
- H01S3/10053—Phase control
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/10—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
- H01S3/102—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating by controlling the active medium, e.g. by controlling the processes or apparatus for excitation
- H01S3/1022—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating by controlling the active medium, e.g. by controlling the processes or apparatus for excitation by controlling the optical pumping
- H01S3/1024—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating by controlling the active medium, e.g. by controlling the processes or apparatus for excitation by controlling the optical pumping for pulse generation
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/10—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
- H01S3/13—Stabilisation of laser output parameters, e.g. frequency or amplitude
- H01S3/1305—Feedback control systems
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/10—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
- H01S3/13—Stabilisation of laser output parameters, e.g. frequency or amplitude
- H01S3/1307—Stabilisation of the phase
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/23—Arrangements of two or more lasers not provided for in groups H01S3/02 - H01S3/22, e.g. tandem arrangements of separate active media
- H01S3/2383—Parallel arrangements
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/23—Arrangements of two or more lasers not provided for in groups H01S3/02 - H01S3/22, e.g. tandem arrangements of separate active media
- H01S3/2308—Amplifier arrangements, e.g. MOPA
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Optics & Photonics (AREA)
- Plasma & Fusion (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
- Lasers (AREA)
- Diffracting Gratings Or Hologram Optical Elements (AREA)
Abstract
The present invention relates to for adjusting the system with the phase of the laser source periodically configured, the system includes: for by the beam collimation from light source and the device being directed on the composed diffraction optical element (1) with periodic phase network, wherein incidence angle is different between each light beam, and incidence angle is determined according to network cycle；For controlling the device of the phase in the source based on the back reaction signal from composite light beam；The device (5) sampled for a part (12) to composite light beam；Fourier lense (6) on the path of the beam section, is provided with composed diffraction optical element (1) in its object plane；Fourier lense (6) as the detector array (7) in plane, be able to detect intensity distribution；For calculating the device (8) of back reaction signal based on these intensity distributions.
Description
Technical field
The field of the invention is the field of the optics coherence tomography in a large amount of fundamental laser light sources.
Background technique
Height in the case of generation and/or ultrashort pulse source of the optics coherence tomography of laser source especially suitable for high power laser source
The generation in energy laser source (for example, pulse width is less than 1 picosecond).
The acquisition in high power (or high-energy) and high brightness laser source is currently limited by the flux stabilized of gain material
Property.A kind of scheme for this problem is that amplification is concurrently distributed on multiple gain medias.This needs to be situated between from each gain
The same phase of laser beam of matter output, so that it is guaranteed that the best optics coherence tomography of whole laser beams.Therefore, it is necessary to dynamically compensate for due to
It is introduced in a large amount of M laser beams by the propagation of the component for the gain media (for example, fiber amplifier) being connected in parallel
Delay.Once locking phase, M newborn laser beam constructive interference, and thus form M times of the brightness that brightness is basic amplifier
Source, and keep its pencil quality (for example, in the case where single mode optical fiber, being limited to diffraction).Therefore, setting and transmitter one
Phaselocked loop more than sample is a problem.
Framework for adjusting the phase of laser source can classify according to multiple standards.The first be light beam spatially
The mode of synthesis or superposition.Thus, Liang Ge race can be distinguished:
Tile aperture synthetic: M laser beam is collimated and has the parallel direction of propagation.This synthesis model is radar wave
The optical equivalence object of beam shaping antenna.Aperture tile with strong main lobe and parasitic secondary lobe.
Fill aperture synthetic: by using polarizer or diffraction optical element (DOE) by M beam combination near field.
The advantages of filling aperture synthetic method is its efficiency, because not having secondary lobe in far field in this case.
Property followed by error signal and will allow to offset in phase between laser source it and it is excellent
Change the processing of their relevant addition.Four kinds of methods for being essentially available for the optics coherence tomography of laser beam are distinguished, according to being included in
Information content classification in negative-feedback signal:
The method for method of referred to as " climbing the mountain ": a part by extracting synthesis energy simply forms error signal, should
Part is maximized by changing the phase in the M channel (light beam) to be synthesized.The technology based on M-1 dimension based on
The optimization algorithm of gradient.In this case, complexity is Processing Algorithm, and the error signal as invariant signal is very simple
It is single and inexpensive.The disadvantages of the method are as follows the bandwidth of loop, changes 1/M.Therefore, this method is more suitable for a small amount of conjunction
At light beam, usually less than 10.
The method of referred to as OHD (optical heterodyne detection) method.In the method, by each transmitter relative to reference beam
Phase measurement composition error signal be vector signal；Each channel uses a detector.It is conciliate by heterodyne mixing
It adjusts and concurrently carries out M measurement.The disadvantages of this method are as follows:
RF component is utilized, which increase the costs in each channel；
By referring to light beam；
Error signal that is measuring before synthesis and not guaranteeing most preferably synthesis quality: it is flat that it cannot compensate phase measurement
Phase fluctuation between face and synthesis plane.Therefore, it is necessary to calibrate to system.
The method of referred to as LOCSET or synchronous multidither method.For hill climbing method, this method uses synthesis energy
A part is used as error signal, but in this case, the contribution from each channel passes through the RF at its proprietary frequency
Modulation carries out " frequency marker " to each channel to identify.Then each light is obtained by using the heterodyne mixing of reference beam
The error signal of beam.This method is advantageous, because it only needs a detector, and the availability of fast phase modulation device
Allow to imagine a large amount of channel.On the other hand, this method needs a large amount of RF in negative feedback loop (frequency mixer, modulator etc.)
Thus component increase the cost in each channel of system significantly.In this case, by with identical frequency in time
Each light beam is sequentially modulated, similar signal is obtained, but there is negative effect to the bandwidth of system.
Directly measure the method for the phase between transmitter, wherein error signal is to interfere with each other from light beam to be synthesized
Or the mapping with the phase of the interference pattern extraction of reference beam interference.This direct interferometric method is common: being passed through
Single image is recorded to obtain all phases by sensor array, and is therefore fully suitable for a large amount of transmitters.It is used
The cost of imager will be divided by the quantity in channel, and be not therefore conclusive.On the other hand, the bandwidth of system can be limited
In used sensor, especially in infrared ray.However, this is not to limit substantially.Finally, for OHD method, synthesizing
Phase is measured before；It cannot compensate phase measurement plane and synthesize the phase fluctuation between plane, therefore cannot be guaranteed best
Synthesis quality.Therefore, it is necessary to calibrate to system.
Following table summarizes the prior art of optics coherence tomography technology.The shortcomings that cell of grey indicates every kind of method.
Therefore, currently without the existing framework of the optics coherence tomography for laser beam, meet following condition: loop simultaneously
Bandwidth > 1kHz, number of beams may be 100,1000 or even higher, operate without calibration (error signal in synthesis plane)
And low cost.
Summary of the invention
Context of the invention is related to the system that wherein laser beam uses diffraction optical element (DOE) spatially to synthesize.
Original use of the system according to the present invention based on the diffraction element also allows to produce other than providing the space combination of light beam
The error signal of raw innovation, the error signal allow to compensate the phase difference between laser source.The error signal is closed by diffraction
It is calculated at the intensity of the higher-order diffraction of element.This error signal allows to meet all above-mentioned conditions.
More specifically, subject of the present invention is a kind of for adjusting with λ0Centered on phase co-wavelength M laser
The system of the phase in source, the M laser source have periodic space configuration, and M is the integer greater than 2, which includes:
M beam collimation for will generate from source and it is directed to the synthesis diffraction optics member with periodic phase grating
Device on part, wherein incidence angle is different between each light beam, these incidence angles are determined according to the period of grating；And
For controlling the device of the phase in the source based on the negative-feedback signal generated from synthesis light beam.
It is mainly characterized by, which includes:
For extracting the device of a part of synthesis light beam；
Fourier lense on the path of the part of synthesis light beam, the fourier lense have object plane and as flat
Face, wherein synthesis diffraction optical element is in its object plane；
Fourier lense as the detector array in plane, be able to detect the intensity point of the part of synthesis light beam
Cloth；
For calculating the device of negative-feedback signal based on these intensity distributions.
It actually obtained vector error signal (its size is provided by the quantity of measured higher-order diffraction), but not
It needs using RF component.Further, since optimization is not directed to locking phase (referring to OHD and direct interferometry technology), but (by making
The intensity of higher-order minimizes) directly against synthesis intensity, so system does not require to calibrate in its principle.
Therefore it is combined with following advantages:
For LOCSET and hill climbing method, error signal is generated in synthesis plane, therefore does not need to calibrate.
Error signal is made of one group of nonredundancy measurement, is enabled and is generated negative-feedback letter by simple processing operation
Number.
The system cost of each channel (light beam) is relatively low, because it is not related to any RF element, and each channel
Only need a detector.
The system is compatible with a large amount of channels, and has the bandwidth greater than 1kHz.
A feature according to the present invention, the device for calculating negative-feedback signal include for by square as defined below
The inverse matrix of battle array calculates the device of the product of the intensity distribution detected in the plane of detector array: the matrix is by leading to
It crosses the coefficient for over one period obtaining the phase unwrapping for synthesizing diffraction optical element at Fourier space to limit, if M
Be odd number then the matrix size be (2M-1) × (2M-1), the size of the matrix is 2M × 2M if M is even number.
In general, M > 100.
Preferably, source is arranged to the configuration of one-dimensional or two-dimensional space.
According to a preferred embodiment of the invention, the light beam generated from laser source has same exit plane, then is
System includes having object plane and another fourier lense as plane, and the exit plane of laser source is located in object plane, and synthesis is spread out
Optical element is penetrated to be located at as in plane.
Detailed description of the invention
It is following as non-limiting example and the specific descriptions that are provided referring to appended attached drawing by reading, of the invention its
Its feature and advantage will be apparent, in the drawing:
Fig. 1 schematically show use diffraction optical element as from " M to 1 " (in the example in the drawings for 5 to
1) device of light beam is synthesized；
Fig. 2 schematically shows the systems of illustrative adjustment phase place according to the present invention；
Fig. 3 schematically shows the intensity distribution of diffraction optical element, which is used as from 1 to M (attached
In the example of figure, from 1 to 5) beam splitter (Fig. 3 a), or be used as from M to 1 (in the example in the drawings, from 5 to 1) synthesize
The device of light beam, wherein intensity distribution passes through diffraction optical element diffraction.
Identical element appended drawing reference having the same in each figure.
Specific embodiment
Context of the invention is such system, is primarily based on and uses diffraction optical element 1 (or DOE) as being used for
The device of each laser beam 10 is synthesized, as shown in the example of figure 1.Laser beam 10 is entered with the angle limited by the space periodic of DOE
It penetrates on DOE 1.Once light beam is phase locked, and has optimum phase distribution (being arranged by DOE), then all light beams are in DOE
Constructive interference at the rank 0 (=main rank) of 11b, and the destructive interference at higher-order 11a.
The concept is that the intensity distribution 11a of diffraction is as error signal at the higher-order of DOE, so as to excellent
It is combined to, as shown in the example of figure 2.In our method, for LOCSET and climb the mountain technology, chain end (that is,
After DOE) measured error signal, to allow to consider all interference suffered by light beam.It will be noted that if the example of Fig. 2
The linear arrangement of laser beam is shown, and therefore illustrates one-dimensional DOE, then the solution proposed is equally applicable to swash
The two-dimensional arrangement and two dimension DOE of light beam.
System according to the present invention for adjusting the phase of M laser source is described referring to Fig. 2.M laser source has
With λ0Centered on phase co-wavelength.These laser sources can be clock；Pulse width might be less that 10-12s。
The system includes:
M phase-modulator: there is a modulator 4 in the output end of each laser source.
DOE 1 is synthesized, with the phase grating with the predetermined space period, positioned at the picture plane of fourier lense 14
In: the M light beam generated from modulator is directed on DOE 1 by the fourier lense 14.Each light beam is by the sky of DOE
Between period restrictions specific incidence angle irradiate DOE.
For extracting the device of a part 12 of synthesis light beam 11, high reflectivity mirror 5 (for example, extracting 1%) can be
Or polarizing cube beam splitter.It is preferably chosen extraction < 1/M.The another part for synthesizing light beam forms the output beam 13 of system.
The second fourier lense 6 in object plane, wherein synthesis DOE 1 is located in object plane.
The matrix of detector 7 in the picture plane (=plane B) of the second fourier lense 6, is able to detect by DOE 1
The intensity distribution 11b, 11a of a part of the order of diffraction of the light beam of synthesis.
For calculating the device 8 of negative-feedback signal from these distributions in the plane of detector array.These computing devices 8
It is connected to M phase-modulator 4, to control them.
M light beam can be adopted is directed into DOE 1 in various manners.In the upstream of DOE, which is included, for example:
Same master oscillator 2 is connected to " 1 arrives M " coupler 3, to generate M laser source；
Possibly, it is respectively connected to the M amplifier 9 of phase-modulator 4.
The exit plane (=plane A) of M laser beam (being generated by amplifier or modulator) is with spacing P's shown in Fig. 2
Space periodicity configuration is located in the object plane of fourier lense 14.
According to an alternative, M laser source has collimation lens associated with each source, and with periodic
Angle and space configuration directly position, so that light beam synthesizes DOE with the specific incidence angle irradiation limited by the space periodic of DOE.
The device 8 for calculating negative-feedback signal is considered now.Therefore, solved the problems, such as by these computing devices be:
The variable of problem is the spatial distribution of the electromagnetic field formed by the superposition of the electromagnetic field generated from each laser source.
It is assumed that the intensity distribution of electromagnetic field is in two independent planes it is known that the exit plane in source is (flat in Fig. 2
Face A) in uniform (or measurement) distribution IAThe distribution I measured later is synthesized with by DOE (in the plane B of Fig. 2)B。
Purpose is to calculate the phase in plane A and BWithDistribution so that from A to B digital communication electromagnetic fieldIt provides
This problem is similar to the survey of the phase aberrations of the image fault from the intensity for example encountered in astronomy
The problem of amount.Illustrative methods for solving this problem can be found in the literature.Following publication can be quoted:
R.G.Paxman and J.R.Fienup " is misaligned sensing and image reconstruction (Optical using the optics of phase difference
Misalignment sensing and image reconstruction using phase diversity) ",
J.Opt.Soc.Am.A 5,914-923 (1988) or J.N.Cederquist, J.R.Fienup, C.C.Wackerman,
S.R.Robinson and D.Kryskowski, " Wave-front phase of Fourier's ionization meter estimates (Wave-front phase
Estimation from Fourier intensity measurements) ", J.Opt.Soc.Am.A 6,1020-1026
(1989) or R.G.Paxman, T.J.Schulz, J.R.Fienup, " by using phase difference Combined estimator object and picture
Poor (Joint estimation of object and aberrations by using phase diversity) "
J.Opt.Soc.Am.A 9,1072-1085 (1992).
The major defect of such method is that they use digital fourier transformation (for during seeking solution
It is related to the calculating of the type of the light propagation from A to B), it may relate to the long calculating time (usually more much longer than 1s).
In a system in accordance with the invention, the main simplification of phase calculation is simply multiplied due to passing through via known matrix
And realize the magnetic distribution from the magnetic distribution Calculation Plane B of plane A.
Specifically, for the M laser beam considered in A, the field distribution in A can be write as according to the parity of M:
If M=2N+1, or:
If M=2N.
Otherwise:
Wherein:
ω is that light beam (it is assumed that Gauss) is with a tight waist in plane A；
P is the period of position of the light beam in plane A；
αkIt is in the weighting coefficient in terms of the amplitude of k-th of light beam (in this case, for example,-N when for M=2N+1
K between+N, and k, α when for M=2N between-N+1 and+Nk=1；Otherwise, αk=0)；
For the light phase of k-th of light beam；
δ is unit pulse (Dirac) function and * is convolution operator.
Field in the plane of DOE passes through EA(x) Fourier transformation and with the phase of DOE in terms of transfer functionIt is multiplied to obtain:
Pass through EDOE(u) Fourier transformation obtains the field for propagating to measurement plane (plane B) again:
Further, since the phase (passing through its structure) of DOE is the periodic function for being 1/P (u is considered in far field) in the period,
SoIt can be write as the form of its Fourier space:
And therefore:
Calculate the synthesis DOE for M light beam to be synthesized to one.On the contrary, by the single light beam of identical DOE diffraction by base
Its intensity is generated on this with identical order (≈ I1) M light beam (referred to as main beam) and compared with low-intensity I2(I2< < I1)
Infinite multiple higher-orders.In other words, it means that coefficient ckFollowing relationship:
K when for M=2N+1 between-N and+N and k when for M=2N between-N+1 and+N, | ck|2≈ 1/M,
Otherwise | ck|2< < 1/M.
Therefore, for following situations, above-mentioned EB(x) the item c of expression formulak+hαhIt can have the value that can not ignore:
For M=2N+1, [- N≤k+h≤+ N] ∪ [- N≤h≤+ N], or
For M=2N, [- N+1≤k+h≤+ N] ∪ [- N+1≤h≤+ N].
Or:
For M=2N+1, k ∈ { -2N+1 ... ,+2N }, or
For M=2N, k ∈ { -2N+1 ... ,+2N }.
Generally, for M=2N+1, EB(x) therefore expression formula is considered for being truncated from -2N to+2N
Index k and h be accurate:
Or for M=2N, -2N+1 to+2N:
Then the expression formula of for M=2N+1 the case where, matrix product are identified and write out as follows:
Or:
Wherein, EA, k(and EB,k, respectively) and it is field EA(x) (and EB(x), complex weighting factor respectively), so that
Near x=kP:
HDOECoefficient c for the Fourier space being launched into over one period by the phase of DOEkThe matrix of restriction.Cause
This, by constructing DOE, which is a priori known.As it appears from the above, in fact, for synthesis M laser beam calculating DOE,
Wherein odd number M is equal to 2N+1, and 2M-1 coefficient is needed in Fourier space, equal to realize:
For being equal to the even number M of 2N, then need 2M coefficient to include
This non-limiting selection that 2M-1 coefficient in the case where M=2N+1 is also shown in Fig. 3 a and 3b, because
To only have the order of diffraction (i.e. for k ∈ { -2N ..., 2N }, in plane B near x=kP) of order -2N to+2N to have for meter
Significant intensity.
Therefore, M light beam from plane A to the light propagation of plane B via if M is odd number size be (2M-1) × (2M-
1), the matrix (matrix H that size is 2M × 2M if M is even numberDOE) and vector EASimple product calculate.In usual side
Fourier transformation used in method is by the matrix HDOEInstead of.Therefore, it is detected in plane B from the array of photoelectric detector
Magnetic distribution calculate magnetic distribution in plane A, pass through the matrix HDOEInverse matrix detected in plane B
The simple product of intensity distribution realize.
According to the magnetic distribution calculated in plane A in this way, phase is calculated in a usual manner.Electromagnetic field in plane A
This simplified calculate of distribution substantially accelerates phase calculation algorithm (=calculating negative-feedback signal), such as iteration or search are most
Big type, and make even for thousands of light beams, these devices for being used to calculate negative-feedback signal can also be implemented in real time.
It is calculated as illustrative iterative phase, the calculating described in following discloses text can be quoted:
- J.Markham and J.A.Conchello, " parametric blind deconvolution: for estimating image and fuzzy robust simultaneously
Method (Parametric blind deconvolution:a robust method for the simultaneous
Estimation of image and blur) ", J.Opt.Soc.Am.A 16 (10), 2377-2391 (1999)；
- J.R.Fienup " phase retrieval algorithm: compares (Phase retrieval algorithms:a
Comparison) ", (15) Appl.Opt.21,2758-2769 (1982).
In this example, synthesis DOE is by transmitting come work, but the system according to the present invention is protected when using reflection DOE
It holds effectively.
Claims (7)
1. one kind has for adjusting with λ0Centered on phase co-wavelength M laser source phase system, the M laser
Source has periodic space configuration, and M is the integer greater than 2, which includes:
It is used for the M beam collimation generated from source and is directed to the synthesis diffraction optical element with periodic phase grating
(1) device on, wherein incidence angle θ2kDifferent between each light beam, these incidence angles are determined according to the period of grating；With
And
For controlling the device of the phase in the source based on the negative-feedback signal generated from synthesis light beam；
For extracting the device (5) of a part (12) of synthesis light beam；
Synthesis light beam the part path on fourier lense (6), the fourier lense (6) have object plane and
As plane, wherein synthesis diffraction optical element (1) is in its object plane；
Fourier lense (6) as the detector array (7) in plane, be able to detect the intensity of the part of synthesis light beam
Distribution；
It is characterized in that, the system includes:
For calculating the device (8) of negative-feedback signal based on these intensity distributions comprising for passing through square as defined below
The inverse matrix of battle array calculates the device of the product of the intensity distribution detected in the plane of detector array: matrix by by
The coefficient for obtaining the phase unwrapping for synthesizing diffraction optical element at Fourier space in a cycle limits, if M is odd
The size of several then matrixes is (2M-1) × (2M-1), and the size of matrix is 2M × 2M if M is even number.
2. according to claim 1 have for adjusting with λ0Centered on the phase of M laser source of phase co-wavelength be
System, which is characterized in that laser source is pulsed.
3. according to claim 2 have for adjusting with λ0Centered on the phase of M laser source of phase co-wavelength be
System, which is characterized in that pulse width is less than 10-12Second.
4. having according to claim 1 described in any one of -3 for adjusting with λ0Centered on phase co-wavelength M laser
The system of the phase in source, which is characterized in that M > 100.
5. according to claim 1 have for adjusting with λ0Centered on the phase of M laser source of phase co-wavelength be
System, which is characterized in that extraction section < 1/M.
6. according to claim 1 have for adjusting with λ0Centered on the phase of M laser source of phase co-wavelength be
System, which is characterized in that source is arranged to the configuration of one-dimensional or two-dimensional space.
7. according to claim 1 have for adjusting with λ0Centered on the phase of M laser source of phase co-wavelength be
System, which is characterized in that when the light beam generated from laser source has same exit plane, the system comprises another Fourier
Lens (14), another described fourier lense (14) have the object plane and synthesis diffraction light where the exit plane of laser source
Learn the picture plane where element (1).
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1401222 | 2014-05-28 | ||
| FR1401222A FR3021761B1 (en) | 2014-05-28 | 2014-05-28 | SYSTEM FOR CONDITIONING A LARGE NUMBER OF LASER SOURCES |
| PCT/EP2015/061524 WO2015181130A1 (en) | 2014-05-28 | 2015-05-26 | System for phasing a large number of laser sources |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106663913A CN106663913A (en) | 2017-05-10 |
| CN106663913B true CN106663913B (en) | 2019-03-15 |
Family
ID=51610152
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201580041454.6A Active CN106663913B (en) | 2014-05-28 | 2015-05-26 | System for adjusting the phase of a large amount of laser sources |
Country Status (10)
| Country | Link |
|---|---|
| US (1) | US9812840B2 (en) |
| EP (1) | EP3149813B1 (en) |
| JP (1) | JP6538084B2 (en) |
| KR (1) | KR102272867B1 (en) |
| CN (1) | CN106663913B (en) |
| ES (1) | ES2704111T3 (en) |
| FR (1) | FR3021761B1 (en) |
| HU (1) | HUE042506T2 (en) |
| IL (1) | IL249158B (en) |
| WO (1) | WO2015181130A1 (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105680317A (en) * | 2016-01-12 | 2016-06-15 | 上海理工大学 | Broadband spectrum beam-splitting beam-combining based wavelength-adjustable light source establishing method |
| US10444526B2 (en) * | 2016-08-01 | 2019-10-15 | The Regents Of The University Of California | Optical pulse combiner comprising diffractive optical elements |
| KR102620450B1 (en) * | 2017-11-07 | 2024-01-03 | 시반 어드밴스드 테크놀러지스 엘티디. | Optical phased array dynamic beam shaping with noise correction |
| FR3076958B1 (en) * | 2018-01-18 | 2021-12-03 | Compagnie Ind Des Lasers Cilas | METHOD AND SYSTEM FOR ADJUSTING THE PROFILE OF A LASER WAVEFRONT |
| DE102018211971A1 (en) | 2018-07-18 | 2020-01-23 | Trumpf Laser Gmbh | Device, laser system and method for combining coherent laser beams |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4649351A (en) * | 1984-10-19 | 1987-03-10 | Massachusetts Institute Of Technology | Apparatus and method for coherently adding laser beams |
| CN101159365A (en) * | 2007-09-27 | 2008-04-09 | 中国人民解放军空军工程大学 | Multi-path optical fiber laser coherent beam combination device and coherent beam combination method |
| CN101251653A (en) * | 2008-04-03 | 2008-08-27 | 哈尔滨工业大学 | Method for preparation of two-dimension optical high order diffraction beam splitter |
| CN101459315A (en) * | 2009-01-04 | 2009-06-17 | 中国人民解放军国防科学技术大学 | Pulse laser based on wavelength-division multiplex coherent synthesizing phase control |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5042922A (en) * | 1986-05-20 | 1991-08-27 | Hughes Aircraft Company | Method for improvidng the spatial resolution in an integrated adaptive optics apparatus |
| US5090795A (en) * | 1987-10-22 | 1992-02-25 | Hughes Aircraft Company | Integrated adaptive optics apparatus |
| JP3840794B2 (en) * | 1998-04-13 | 2006-11-01 | 富士ゼロックス株式会社 | Laser drive device |
| WO2006063474A1 (en) * | 2004-12-16 | 2006-06-22 | Vectronix Ag | Not temperature stabilized pulsed laser diode and all fibre power amplifier |
| US7440174B2 (en) * | 2006-02-24 | 2008-10-21 | Northrop Grumman Corporation | Coherent fiber diffractive optical element beam combiner |
| US7924894B2 (en) * | 2008-01-18 | 2011-04-12 | Northrop Grumman Systems Corporation | Digital piston error control for high-power laser system employing diffractive optical element beam combiner |
| US7756169B2 (en) * | 2008-01-23 | 2010-07-13 | Northrop Grumman Systems Corporation | Diffractive method for control of piston error in coherent phased arrays |
| US7821900B2 (en) * | 2008-05-15 | 2010-10-26 | Northrop Grumman Systems Corporation | Diffractive optical element and method of designing the same |
| WO2013155533A1 (en) * | 2012-04-13 | 2013-10-17 | California Institute Of Technology | Integrated light source independent linewidth reduction of lasers using electro-optical feedback techniques |
| US8792526B2 (en) * | 2012-05-24 | 2014-07-29 | Northrop Grumman Systems Corporation | Coherent laser array control system and method |
-
2014
- 2014-05-28 FR FR1401222A patent/FR3021761B1/en not_active Expired - Fee Related
-
2015
- 2015-05-26 KR KR1020167036621A patent/KR102272867B1/en active IP Right Grant
- 2015-05-26 HU HUE15725588A patent/HUE042506T2/en unknown
- 2015-05-26 ES ES15725588T patent/ES2704111T3/en active Active
- 2015-05-26 WO PCT/EP2015/061524 patent/WO2015181130A1/en active Application Filing
- 2015-05-26 EP EP15725588.6A patent/EP3149813B1/en active Active
- 2015-05-26 CN CN201580041454.6A patent/CN106663913B/en active Active
- 2015-05-26 JP JP2016569820A patent/JP6538084B2/en active Active
- 2015-05-26 US US15/313,542 patent/US9812840B2/en active Active
-
2016
- 2016-11-23 IL IL249158A patent/IL249158B/en active IP Right Grant
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4649351A (en) * | 1984-10-19 | 1987-03-10 | Massachusetts Institute Of Technology | Apparatus and method for coherently adding laser beams |
| CN101159365A (en) * | 2007-09-27 | 2008-04-09 | 中国人民解放军空军工程大学 | Multi-path optical fiber laser coherent beam combination device and coherent beam combination method |
| CN101251653A (en) * | 2008-04-03 | 2008-08-27 | 哈尔滨工业大学 | Method for preparation of two-dimension optical high order diffraction beam splitter |
| CN101459315A (en) * | 2009-01-04 | 2009-06-17 | 中国人民解放军国防科学技术大学 | Pulse laser based on wavelength-division multiplex coherent synthesizing phase control |
Also Published As
| Publication number | Publication date |
|---|---|
| IL249158A0 (en) | 2017-01-31 |
| IL249158B (en) | 2020-10-29 |
| KR102272867B1 (en) | 2021-07-02 |
| ES2704111T3 (en) | 2019-03-14 |
| EP3149813A1 (en) | 2017-04-05 |
| KR20170012445A (en) | 2017-02-02 |
| HUE042506T2 (en) | 2019-07-29 |
| FR3021761B1 (en) | 2018-02-09 |
| JP6538084B2 (en) | 2019-07-03 |
| CN106663913A (en) | 2017-05-10 |
| JP2017517153A (en) | 2017-06-22 |
| WO2015181130A1 (en) | 2015-12-03 |
| FR3021761A1 (en) | 2015-12-04 |
| US20170201063A1 (en) | 2017-07-13 |
| US9812840B2 (en) | 2017-11-07 |
| EP3149813B1 (en) | 2018-10-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106663913B (en) | System for adjusting the phase of a large amount of laser sources | |
| US10911142B2 (en) | Distributed array for direction and frequency finding | |
| CN103472457B (en) | Three-dimensional imaging system and method for calculating correlation flight time by means of sparse aperture compression | |
| CN103472455B (en) | Four-dimensional spectral imaging system and method for calculating correlation flight time by means of sparse aperture compression | |
| CN103472456B (en) | Active imaging system and method based on sparse aperture compressing calculation correlation | |
| Pen et al. | The GMRT EoR experiment: limits on polarized sky brightness at 150 MHz | |
| Blanchard et al. | Coherent optical beam forming with passive millimeter-wave arrays | |
| US8681047B2 (en) | System and method to form coherent wavefronts for arbitrarily distributed phased arrays | |
| CN101975947B (en) | Two-dimensional radial imaging method for mirror image synthetic aperture | |
| US20140192161A1 (en) | Three-dimensional reconstruction of a millimeter-wave scene by optical up-conversion and cross-correlation detection | |
| Molaei et al. | Development of fast Fourier-compatible image reconstruction for 3D near-field bistatic microwave imaging with dynamic metasurface antennas | |
| Chen et al. | Astrometry in two-photon interferometry using an earth rotation fringe scan | |
| CN108613742B (en) | Double-channel large-aperture static interference imaging spectrometer | |
| Ryan et al. | Instantaneous microwave-photonic spatial-spectral channelization via k-space imaging | |
| Hill et al. | Photon noise correlations in millimeter-wave telescopes | |
| Liu et al. | Video-rate terahertz interferometric and synthetic aperture imaging | |
| Denolle et al. | Real-time focal-plane wavefront sensing for compact imaging phased-array telescopes: numerical and experimental demonstration | |
| DeSantis et al. | Multi-beam solid-state LiDAR using star-coupler-based optical phased arrays | |
| US10267684B1 (en) | Device, system, and technique for characterizing a focal plane array modulation transfer function | |
| Rahmer et al. | Aircraft avoidance for laser propagation at the Large Binocular Telescope Observatory: life under a busy airspace | |
| McKean et al. | Calibration of LOFAR | |
| Yoo | Low-Mass Planar Photonic Imaging Sensor | |
| Murakowski et al. | RF k-space Tomography | |
| CN102305932A (en) | Moving target imaging method for Fresnel telescope imaging laser radar | |
| Laviada et al. | Phase retrieval technique for submillimetre‐wave frequency scanning‐based radar system |
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 |