CN108519671B - Closed-loop correction control method for phase translation error of splicing telescope system - Google Patents
Closed-loop correction control method for phase translation error of splicing telescope system Download PDFInfo
- Publication number
- CN108519671B CN108519671B CN201810356033.8A CN201810356033A CN108519671B CN 108519671 B CN108519671 B CN 108519671B CN 201810356033 A CN201810356033 A CN 201810356033A CN 108519671 B CN108519671 B CN 108519671B
- Authority
- CN
- China
- Prior art keywords
- mirror
- sub
- central sub
- closed
- error
- 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
Images
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B26/00—Optical devices or arrangements for the control of light using movable or deformable optical elements
- G02B26/06—Optical devices or arrangements for the control of light using movable or deformable optical elements for controlling the phase of light
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B23/00—Telescopes, e.g. binoculars; Periscopes; Instruments for viewing the inside of hollow bodies; Viewfinders; Optical aiming or sighting devices
- G02B23/02—Telescopes, e.g. binoculars; Periscopes; Instruments for viewing the inside of hollow bodies; Viewfinders; Optical aiming or sighting devices involving prisms or mirrors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Astronomy & Astrophysics (AREA)
- Telescopes (AREA)
- Mechanical Light Control Or Optical Switches (AREA)
Abstract
The invention relates to a closed-loop correction control method for phase translation errors of a splicing telescope system, which comprises the following steps: after coarse phase sharing of the sparse aperture telescope is achieved, two-dimensional dispersion interference fringes between each pair of sub-mirrors are collected in real time, a target function value representing the size of a phase translation error and a variable representing the positive and negative of the phase translation error are calculated, then the step length of closed-loop iterative correction and the direction of closed-loop correction are selected according to the target function value and the positive and negative variable respectively to conduct real-time closed-loop iterative correction on the phase translation error between each pair of sub-mirrors, and finally the phase translation error is corrected within a set error range. The method does not need to carry out any physical modification on the existing equipment, is simple to realize, does not need to carry out any physical calibration in advance, does not have multiple complex data operations, and has good stability and robustness for the closed-loop correction of the phase translation error.
Description
Technical Field
The invention relates to a closed-loop correction control method for phase translation errors of a splicing telescope system.
Background
The telescope with high spatial resolution is used for better observing celestial bodies and universes, which are dreams of astronomysts and also targets pursued by astronomical instrument builders. The spatial resolution of the telescope is proportional to its aperture, i.e. the larger the aperture, the higher the resolution of the telescope. At present, the construction of a single-mirror telescope with more than 8 meters is limited by factors such as manufacturing process, manufacturing cost and the like, and if the spatial resolution of the telescope is to be improved, a splicing mirror technology can be adopted to construct a telescope with a larger aperture.
Besides astronomical application, the spliced telescope has good application prospects in the fields of aerospace remote sensing and the like, the research on the aspect of earth observation is also beneficial to improving the accuracy of forecasting on earth climate, space climate, solar activity time and the like, the influence of geological disasters, severe climate environments, solar activity change and the like on the earth and human is reduced, and the spliced telescope has great social value on the development of national economy and the improvement of quality and level of life of people.
The key of the spliced telescope for obtaining the image close to the diffraction limit resolution is to ensure that the light beams from the sub-mirrors are coherently superposed in an in-situ manner on a focal plane in the whole field of view, and the in-situ error must be stably controlled within one tenth of a wavelength range, which is a precondition for realizing stable Fizeau interference imaging. Therefore, the translational error control technique in the splicing telescope has become one of the hot spots of research in the related art. Currently, the proposed translational error detection and closed-loop control techniques mainly include the following:
1. phase difference method
2. Closed-loop common-phase control method based on rectangular pyramid sensor
The southern european astronomical observatory studies the application of the rectangular pyramid sensor in the co-phase detection. In the related experiment, besides measuring the inclination and the high-order aberration, the rectangular pyramid wavefront sensor can also measure the phase shift error between each sub-mirror. Experimental results show that the method has high measurement precision, the translation error after closed-loop correction reaches 5.7nm, but the method can only measure the optical path difference in a wavelength range (see: optical sensor for segmented calibration, Optics Letters, 2005, 30(19): 2572-. When using dual wavelength lambda1And λ2When closed-loop common-phase correction is carried out, the correctable maximum translation error does not exceed lambda1λ2/[4(λ1-λ2)]Therefore, the method can usually correct the translation error only in a range of a few micrometers.
3. Dispersion fringe method
Dispersion fringe method was proposed by Fang Shi et al, university of California, USA, for the detection of co-phase between the sub-mirrors of a Keck telescope (see: Experimental verification of discrete segmenting as a segmented phase detecting technique using the Keck telescope, applied optics Vol.43, Issue 23, pp.4474-4481 (2004)). Simulation calculation and experimental results show that the method is only used in a coarse common-phase detection stage, and when the absolute translation error is less than half wavelength, the method fails and needs other detection methods.
The Chinese patent application with the application number of 200810000577.7 provides a two-dimensional dispersion fringe analysis method for absolute distance measurement, and the method has the advantages of large measurement range and high measurement precision. However, in this method, the main peak position corresponding to each wavelength when the absolute distance is zero needs to be calibrated, which is difficult to achieve in practical use because: firstly, it is difficult to control the absolute distance between two sub-mirrors to be zero, and other detection means must be used; secondly, the calibration optical path and the actual measurement optical path are usually two different optical paths, and then, or the existence of factors such as temperature change, external vibration, atmospheric turbulence and the like can cause that the position of a main peak corresponding to each wavelength when the calibration absolute distance is zero has a larger deviation with the position of a main peak corresponding to each wavelength in an actual system, thereby finally causing the failure of the whole dispersion fringe analysis method.
The method needs to obtain a nonlinear relation between a peak value ratio and a translation error in advance through measurement or simulation, namely, the nonlinear relation between the peak value ratio of interference fringes and the translation error is calibrated, but external environment changes can cause uncertain random aberration to be generated in an actual system, so that the relation between the peak value ratio of the interference fringes and the translation error in the actual system is not matched with a calibrated value, the measurement result is influenced, and the correction effect is not good finally.
In view of the above, there is a need for an improved method for correcting phase shift errors of a spliced telescope.
Disclosure of Invention
In order to solve the problems in the prior art, the invention aims to provide a closed-loop correction control method for the phase translation error of a splicing telescope system, so as to overcome the defects of the existing phase translation error detection and correction control technology, realize the direct closed-loop correction of the phase translation error in the splicing telescope system within a millimeter-scale phase translation error range, control the final correction error within a tenth wavelength range, and do not need to carry out any preliminary calibration, and do not increase the complexity and the cost in engineering.
The invention relates to a closed-loop correction control method for phase translation errors of a splicing telescope system, wherein the splicing telescope system comprises a plurality of sub-mirrors, and the method comprises the following steps of:
s1, setting a sub-mirror close to the center of the telescope field of view as a central sub-mirror, and collecting a plurality of two-dimensional dispersion interference fringes between n non-central sub-mirrors and the central sub-mirror;
step S2, calculating the objective function value J between the ith non-central sub-mirror and the central sub-mirror according to the two-dimensional dispersion interference fringe between the ith non-central sub-mirror and the central sub-mirroriI is 1 … n; if Ji＜ασ/λ0Then the step ends, otherwise step S3 is executed, where σ is the set phase shift error closed loop correction residual, λ0For the center wavelength of the broadband light, α is a function of the value J characterizing the objective function and a center wavelength λ0A proportionality coefficient of approximate linear relation between phase translation errors between the non-central sub-mirror and the central sub-mirror within the range;
step S3, calculating to obtain a variable S representing the positive and negative of the phase shift error between the ith non-central sub-mirror and the central sub-mirrori；
Step S4, if the objective function value J between the ith non-central sub-mirror and the central sub-mirroriWhen the position is larger than α, the i-th non-central sub-mirror is made to generate a translation Pi＝-SPSiWherein S isPFor the first correction step length, otherwise make the i-th non-central sub-mirror generate the translation amount Pi＝-SLSiWherein S isLIs a second correction step length;
step S5, the current closed-loop correction process is completed, and when the next closed-loop correction process is entered, the step S1 is executed.
In the above method for controlling closed-loop correction of phase shift error of a splicing telescope system, the step S2 includes:
step S21, according to the two-dimensional dispersion interference fringe between the ith non-central sub-mirror and the central sub-mirror, establishing a coordinate system with the dispersion direction as the x axis and the baseline direction of the pair of sub-mirrors as the y axis, wherein the two-dimensional dispersion interference fringe comprises N one-dimensional sub-fringes along the y axis direction;
step S22, obtaining the second peak intensity I of the jth one-dimensional sub-stripe2(j) And a third peak intensity I3(j) Wherein, in the step (A),j＝1…N；
step S23, calculating the peak ratio r (j) of the jth one-dimensional sub-stripe according to formula (1):
R(j)＝I2(j)/I3(j)-1，j＝1…N (1)，
step S24, calculating the objective function value J between the ith non-central sub-mirror and the central sub-mirror according to the formula (2)i：
Wherein β represents the order of the center-to-center distance, which is a positive integer greater than 1.
In the above method for controlling closed-loop correction of phase shift error of a splicing telescope system, the step S3 includes:
step S31, according to the two-dimensional dispersion interference fringe between the ith non-central sub-mirror and the central sub-mirror, establishing a coordinate system with the dispersion direction as the x axis and the baseline direction of the pair of sub-mirrors as the y axis, wherein the two-dimensional dispersion interference fringe comprises N one-dimensional sub-fringes along the y axis direction;
step S32, obtaining the first peak intensity I of the jth one-dimensional sub-stripe1(j) And the first peak coordinate y1(j) Second peak intensity I2(j) And a second peak position y2(j)，j＝1…N；
Step S33, calculating a peak ratio q (j) of positive and negative information of the translation error of the jth one-dimensional sub-stripe according to formula (3):
Q(j)＝[I2(j)/I1(j)]·sign[y1(j)-y2(j)]，j＝1…N (3)，
wherein sign [ ] is a sign function;
step S34, calculating and obtaining a variable S representing the positive and negative of the phase translation error between the ith non-central sub-mirror and the central sub-mirror according to the formula (4)i：
Wherein sign [ ] is a sign function, and T is a positive integer.
In the closed-loop correction control method for the phase shift error of the splicing telescope system, the first correction step length SPHas a value range of [ sigma, lambda%0]。
In the above closed-loop correction control method for the phase shift error of the splicing telescope system, the second correction step SLHas a value range of [1, sigma ]]。
By adopting the technical scheme, after the coarse phase sharing is realized in the sparse aperture telescope, the two-dimensional dispersion interference fringes between each pair of sub-mirrors are collected in real time, the objective function value representing the size of the phase translation error and the variables representing the positive and negative of the phase translation error are calculated, then the step length of closed-loop iterative correction and the direction of closed-loop correction are selected according to the objective function value and the positive and negative variables respectively to carry out real-time closed-loop iterative correction on the phase translation error between each pair of sub-mirrors, and finally the phase translation error is corrected within a set error range. Compared with the prior art, the method improves the real-time performance of the phase translation error correction, and does not need to calibrate in advance in the whole process of detecting the closed-loop correction from the translation error, thereby improving the precision and the stability of the phase translation error closed-loop correction.
Drawings
FIG. 1 is a schematic diagram of the construction of a splicing telescope system according to the invention;
FIG. 2 is a diagram showing the comparison between the two-dimensional dispersion interference fringe light intensity distribution and the one-dimensional sub-fringe light intensity distribution;
FIG. 3 is a diagram illustrating the variation of objective function value with phase shift error in the present invention
Fig. 4 is a schematic diagram of the variation of the phase shift error with iterative correction control in the present invention.
Detailed Description
The preferred embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
The invention relates to a closed-loop correction control method for a phase translation error of a splicing telescope system, which is used for closed-loop correction of the phase translation error of the splicing telescope system shown in figure 1.
As shown in fig. 1, the split telescope system includes: a splicing telescope 10, a collimator lens 3, a first beam splitter 4, a second beam splitter 5, an imaging system 6, a wavefront detector 7, a translational error detector 8, and a wavefront controller 9, wherein,
the splicing telescope 10 receives peripheral incident light (such as starlight or broadband light), and specifically includes: a sub-mirror 1 and a plurality of sub-mirrors 2 each having a driver 20; in this embodiment, the telescope 10 is a cassegrain-type reflector telescope; the number of the sub-mirrors 2 is two, and the sub-mirrors are used for receiving and reflecting incident light; the plurality of drivers 20 mounted on the back of the sub-mirror 2 can realize six-degree-of-freedom motion and have millimeter-scale stroke and nanometer-scale displacement precision; the secondary mirror 1 is used for receiving the incident light reflected by the secondary mirror 2 and reflecting and converging the incident light in front of the secondary mirror 2 again;
the collimating mirror 3 is used for receiving the incident light reflected again by the secondary mirror 1 and generating parallel (or basically parallel) emergent light;
the first beam splitter 4 receives the emergent light from the collimating mirror 3, transmits a first split beam to the wavefront detector 7, and transmits a second split beam to the second beam splitter 5;
the second beam splitter 5 receives the second split beam, transmits the third split beam to the translational error detector 8, and transmits the fourth split beam to the imaging system 6 for imaging;
the wavefront detector 7 detects the first split beam and obtains system aberration data (including tilt and higher-order aberration) on the transmission path of each sub-mirror 2; in the present embodiment, the wavefront sensor 7 may be a hartmann-shack wavefront sensor, a pyramid wavefront sensor, or an interferometer;
the translational error detector 8 detects the third split beam and collects dispersion interference fringes between each pair of sub-mirrors 2;
the wavefront controller 9 is respectively connected with the wavefront detector 7, the translational error detector 8 and the driver of the sub-mirror 2, and firstly receives and outputs a first driving voltage to the driver 20 of the sub-mirror 2 according to the systematic offset error data of the sub-mirror 2 so as to drive the sub-mirror 2 to generate corresponding translation, inclination and rotation, thereby correcting the systematic offset error of the sub-mirror 2; then, the chromatic dispersion interference fringes between the sub-mirrors 2 collected by the translational error detector 8 are received, the magnitude and the direction of the phase translational error between the sub-mirrors 2 are calculated, and a second driving voltage is output to the driver 20 of the sub-mirror 2 to drive the sub-mirror 2 to generate translation, so that the compensation correction is performed on the phase translational error (the content of the part is the closed loop correction control method for the phase translational error of the splicing telescope system of the present invention which will be described in detail below).
Specifically, the method of the present invention comprises: after the systematic maladjustment errors except the phase shift error are corrected by the wavefront detector 7 and the wavefront controller 9 to make the sub-mirror 2 in the rough co-phase state (it should be understood that the co-phase errors mainly include the shift, tilt and defocus aberrations, and the rough co-phase, i.e. the tilt and defocus aberrations are corrected, the rough co-phase implementation process is relatively easy, no specific sensor is needed, and the rough co-phase can be implemented generally through a system far field image), the following steps are performed:
step S1, setting the sub-mirror 2 close to the center of the telescope field of view as a central sub-mirror (also called as a central sub-aperture), and collecting a plurality of two-dimensional dispersion interference fringes between other non-central sub-mirrors and the central sub-mirror (a two-dimensional dispersion interference fringe is formed between each pair of sub-mirrors 2) by using the translational error detector 8;
step S2, calculating an objective function value J between the ith non-central sub-mirror and the central sub-mirror according to the two-dimensional dispersion interference fringe between the ith non-central sub-mirror and the central sub-mirror by using the wavefront controller 9iI is 1 … n, n is the number of non-central sub-mirrors; if Ji＜ασ/λ0Then the process ends (this time, it means that the phase shift error between the ith non-central sub-mirror and the central sub-mirror is within the set phase shift error closed loop correction residual error range, so that it is not necessary to correct the phase shift error), otherwise (i.e. if J is performed)i≥ασ/λ0) Step S3 is executed (in this case, it is explained that the phase shift error between the i-th non-central sub-mirror and the central sub-mirror exceeds the set phase shift error closed-loop correction residual error range, so that it is necessary to correct the phase shift errorDifference correction), where σ is the set phase shift error closed loop correction residual, λ0The central wavelength of the broadband light (usually sigma ≦ λ) is used for collecting the two-dimensional dispersive interference fringes010) α characterizing the value of the objective function J and a central wavelength λ0The proportionality coefficient of approximate linear relation between phase shift errors between the non-central sub-mirror and the central sub-mirror in the range is when | | < lambda0When the phase shift error is equal to the target function value J, the approximate linear relation is satisfied, namely J is approximately equal to α/lambda0；
It should be noted that the above-mentioned proportionality coefficient α is mainly related to system detector parameters, and can be obtained by fitting in advance (for example, 0.25 is taken in this embodiment);
step S3, calculating to obtain a variable S representing the positive and negative of the phase shift error between the ith non-central sub-mirror and the central sub-mirrori(i.e., the direction characterizing the phase shift error); in particular, the phase shift error between the ith non-central sub-mirror and the central sub-mirroriWhen > 0, SiWhen 1 is equal toiWhen < 0, Si＝-1；
Step S4, if the objective function value J between the ith non-central sub-mirror and the central sub-mirrori(> α), a corresponding drive voltage (i.e., the second drive voltage described above) is applied by the wavefront controller 9 to the i-th non-central sub-mirror driver 20 to cause the i-th non-central sub-mirror to translate by an amount Pi＝-SPSi(unit is nm) in which SPIs a first correction step length with a value range of [ sigma, lambda0]σ is the set phase shift error closed loop correction residual, SPGenerally taken as0/3 (for example, in the present embodiment, the measurement wavelength of the broadband light used is 550nm to 650nm, and the center wavelength λ thereof is0Is 600nm, and therefore, the first correction step SPTaking the central wavelength λ0One third of (i.e., get 200), otherwise (i.e., if J)iα) to generate a translation P for the ith non-central sub-mirror by applying a corresponding drive voltage to the driver 20 for the ith non-central sub-mirror via the wavefront controller 9i＝-SLSi(in nm), wherein,SLis the second correction step length with the value range of [1, sigma ]]σ is the phase shift error closed loop correction residual, SLGenerally, σ/2 (e.g., σ ═ λ in the present embodiment)0/30，SL＝σ/2)；
It is noted that JiThe common phase error is about one wavelength when the common phase error is α, so as to determine whether the common phase error is more than one wavelength, if so, the maximum correction step can be lambda0(ii) a In addition, the first correction step SPFor correction steps with co-phase error greater than one wavelength, and a second correction step SLThe correction step length is adopted when the common-phase error is smaller than one wavelength, and the two step lengths are adopted to ensure the speed of correction iteration and simultaneously give consideration to the correction precision;
step S5, the current closed-loop detection and correction process is completed, and when the next closed-loop detection and correction process is entered, the process returns to step S1.
The step S2 specifically includes the following steps:
step S21, according to the two-dimensional dispersion interference fringe between the ith non-central sub-mirror and the central sub-mirror, establishing a coordinate system with the dispersion direction as the x axis and the baseline direction of the pair of sub-mirrors as the y axis, wherein the two-dimensional dispersion interference fringe comprises N one-dimensional sub-fringes along the y axis direction;
step S22, obtaining the second peak intensity I of the jth one-dimensional sub-stripe2(j) And a third peak intensity I3(j) (as shown in fig. 2), wherein j ═ 1 … N;
step S23, calculating the peak ratio r (j) of the jth one-dimensional sub-stripe according to formula (1):
R(j)＝I2(j)/I3(j)-1，j＝1…N (1)，
wherein [ I ]2(j)/I3(j)-1]It holds true that not less than 0 is constant, so that the objective function value JiIs 0 and the minimum value is independent of system parameters and is applicable to any system;
step S24, calculating the objective function value J between the ith non-central sub-mirror and the central sub-mirror according to the formula (2)i：
Where β represents the order of the center distance, which is a positive integer greater than 1 (for example, in this embodiment, 2 is taken, that is, formula (2) represents calculating the center distance of 2 orders).
The step S3 specifically includes:
step S31, according to the two-dimensional dispersion interference fringe between the ith non-central sub-mirror and the central sub-mirror, establishing a coordinate system with the dispersion direction as the x axis and the baseline direction of the pair of sub-mirrors as the y axis, wherein the two-dimensional dispersion interference fringe comprises N one-dimensional sub-fringes along the y axis direction;
step S32, obtaining the first peak intensity I of the jth one-dimensional sub-stripe1(j) And the first peak coordinate y1(j) Second peak intensity I2(j) And a second peak position y2(j)，j＝1…N；
Step S33, calculating a peak ratio q (j) of positive and negative information of the translation error of the jth one-dimensional sub-stripe according to formula (3):
Q(j)＝[I2(j)/I1(j)]·sign[y1(j)-y2(j)]，j＝1…N (3)，
wherein sign [ ] is a sign function;
step S34, calculating and obtaining a variable S representing the positive and negative of the phase translation error between the ith non-central sub-mirror and the central sub-mirror according to the formula (4)i：
Wherein sign [ ] is a sign function, T is a positive integer, and the value range is [5, N/10] (for example, 10 is taken in this embodiment).
In this example, the initial translational phase error between the two sub-mirrors is 5 μm, and the measurement wavelength used is 550nm-650nm (central wavelength λ)0600nm) under a bandwidth light condition. FIG. 3 is a correlation between an objective function J and a phase shift errorWhen the phase shift error is in a central wavelength range, the least square method is used for fitting to obtain an approximate linear relation between the two, wherein the approximate linear relation is as follows: j ≈ 0.25/lambda0. FIG. 4 shows the variation of the phase shift error with the number of iterative corrections, with the initial phase shift error of 5000nm and the final set closed-loop correction residual error of λ0And 30 (namely 20nm), the phase shift error becomes 15nm after 28 closed loop iterative corrections, and the set correction precision is achieved.
In summary, compared with the prior art, the invention has the following advantages:
1. the invention provides a method for directly judging the positive and negative signs of the phase shift error from two-dimensional dispersion interference fringe data, which effectively solves the 2 pi ambiguity problem and sign judgment in the phase shift error detection;
2. compared with the existing two-dimensional dispersion fringe analysis method and the phase translation error method based on the far field similarity, the method does not need to carry out any physical calibration in advance, and has better stability and robustness for the closed-loop correction of the phase translation error;
3. the invention has no multi-dimensional complex data operation, and improves the real-time performance of the phase shift error closed-loop correction;
4. the invention can be realized in the existing telescope common-phase device without any physical modification on the existing equipment, and is simple to realize.
The above embodiments are merely preferred embodiments of the present invention, which are not intended to limit the scope of the present invention, and various changes may be made in the above embodiments of the present invention. All simple and equivalent changes and modifications made according to the claims and the content of the specification of the present application fall within the scope of the claims of the present patent application. The invention has not been described in detail in order to avoid obscuring the invention.
Claims (4)
1. A closed loop correction control method for phase shift errors in a spliced telescope system comprising a plurality of sub-mirrors, the method comprising, after placing the sub-mirrors in a coarse common phase state, performing the steps of:
s1, setting a sub-mirror close to the center of the telescope field of view as a central sub-mirror, and collecting a plurality of two-dimensional dispersion interference fringes between n non-central sub-mirrors and the central sub-mirror;
step S2, calculating the objective function value J between the ith non-central sub-mirror and the central sub-mirror according to the two-dimensional dispersion interference fringe between the ith non-central sub-mirror and the central sub-mirroriI is 1 … n; if Ji＜ασ/λ0Then the step ends, otherwise step S3 is executed, where σ is the set phase shift error closed loop correction residual, λ0For the center wavelength of the broadband light, α is a function of the value J characterizing the objective function and a center wavelength λ0A proportionality coefficient of approximate linear relation between phase translation errors between the non-central sub-mirror and the central sub-mirror within the range;
the step S2 includes:
step S21, according to the two-dimensional dispersion interference fringe between the ith non-central sub-mirror and the central sub-mirror, establishing a coordinate system with the dispersion direction as the x axis and the baseline direction of the pair of sub-mirrors as the y axis, wherein the two-dimensional dispersion interference fringe comprises N one-dimensional sub-fringes along the y axis direction;
step S22, obtaining the second peak intensity I of the jth one-dimensional sub-stripe2(j) And a third peak intensity I3(j) Wherein j is 1 … N;
step S23, calculating the peak ratio r (j) of the jth one-dimensional sub-stripe according to formula (1):
R(j)＝I2(j)/I3(j)-1，j＝1…N (1)，
step S24, calculating the objective function value J between the ith non-central sub-mirror and the central sub-mirror according to the formula (2)i：
Wherein, beta represents the order of the center distance, and is a positive integer greater than 1;
step S3, calculating to obtain the positive and negative characteristics of the phase shift error between the ith non-central sub-mirror and the central sub-mirrorVariable S ofi；
Step S4, if the objective function value J between the ith non-central sub-mirror and the central sub-mirroriWhen the position is larger than α, the i-th non-central sub-mirror is made to generate a translation Pi＝-SPSiWherein S isPFor the first correction step length, otherwise make the i-th non-central sub-mirror generate the translation amount Pi＝-SLSiWherein S isLIs a second correction step length;
step S5, the current closed-loop correction process is completed, and when the next closed-loop correction process is entered, the step S1 is executed.
2. The closed-loop correction control method for the phase shift error of the splicing telescope system as recited in claim 1, wherein the step S3 comprises:
step S31, according to the two-dimensional dispersion interference fringe between the ith non-central sub-mirror and the central sub-mirror, establishing a coordinate system with the dispersion direction as the x axis and the baseline direction of the pair of sub-mirrors as the y axis, wherein the two-dimensional dispersion interference fringe comprises N one-dimensional sub-fringes along the y axis direction;
step S32, obtaining the first peak intensity I of the jth one-dimensional sub-stripe1(j) And the first peak coordinate y1(j) Second peak intensity I2(j) And a second peak position y2(j)，j＝1…N；
Step S33, calculating a peak ratio q (j) of positive and negative information of the translation error of the jth one-dimensional sub-stripe according to formula (3):
Q(j)＝[I2(j)/I1(j)]·sign[y1(j)-y2(j)]，j＝1…N (3)，
wherein sign [ ] is a sign function;
step S34, calculating and obtaining a variable S representing the positive and negative of the phase translation error between the ith non-central sub-mirror and the central sub-mirror according to the formula (4)i：
Wherein sign [ ] is a sign function, and T is a positive integer.
3. The closed-loop correction control method for phase shift error of splicing telescope system according to claim 1, wherein the first correction step SPHas a value range of [ sigma, lambda%0]。
4. The method for closed-loop correction control of phase shift error in a splicing telescope system as recited in claim 1, wherein the second correction step S is performedLHas a value range of [1, sigma ]]。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810356033.8A CN108519671B (en) | 2018-04-19 | 2018-04-19 | Closed-loop correction control method for phase translation error of splicing telescope system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810356033.8A CN108519671B (en) | 2018-04-19 | 2018-04-19 | Closed-loop correction control method for phase translation error of splicing telescope system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108519671A CN108519671A (en) | 2018-09-11 |
| CN108519671B true CN108519671B (en) | 2020-09-18 |
Family
ID=63429773
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810356033.8A Active CN108519671B (en) | 2018-04-19 | 2018-04-19 | Closed-loop correction control method for phase translation error of splicing telescope system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108519671B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109709668A (en) * | 2019-02-02 | 2019-05-03 | 北京空间机电研究所 | A kind of automatic phase control unit and adjusting method |
| CN110109230B (en) * | 2019-05-24 | 2020-07-28 | 西安交通大学 | Intelligent splicing assembly method for aspheric surface complex curved surface workpiece |
| CN111207910B (en) * | 2020-01-16 | 2021-05-25 | 中国科学院长春光学精密机械与物理研究所 | Spliced mirror common-phase error correction method based on dispersion fringe slope analysis |
| CN113917686A (en) * | 2021-10-19 | 2022-01-11 | 中国科学院光电技术研究所 | Image-based splicing diffraction telescope splicing error parallel correction method |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6649895B1 (en) * | 2000-07-14 | 2003-11-18 | Adaptive Optics Associates, Inc. | Dispersed Hartmann sensor and method for mirror segment alignment and phasing |
| CN100573038C (en) * | 2008-01-23 | 2009-12-23 | 北京理工大学 | The two-dimension chromatic dispersion fringe analysis method that is used for absolute distance measurement |
| US9200966B2 (en) * | 2013-07-24 | 2015-12-01 | Raytheon Company | Dual field of view telescope |
| CN105115427B (en) * | 2015-09-16 | 2017-12-26 | 中国科学院光电技术研究所 | A kind of method for improving dispersion striped sensor common phase system noise resisting ability |
| CN105700128B (en) * | 2016-05-03 | 2018-06-19 | 中国科学院上海天文台 | A kind of spliced telescope common phase level control device and control method |
| CN107450176B (en) * | 2017-09-21 | 2019-10-29 | 中国科学院上海天文台 | A kind of space sparse aperture telescope common phase control device and control method |
-
2018
- 2018-04-19 CN CN201810356033.8A patent/CN108519671B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108519671A (en) | 2018-09-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108519671B (en) | Closed-loop correction control method for phase translation error of splicing telescope system | |
| CN105700128B (en) | A kind of spliced telescope common phase level control device and control method | |
| Ellerbroek et al. | Inverse problems in astronomical adaptive optics | |
| CN107450176B (en) | A kind of space sparse aperture telescope common phase control device and control method | |
| CN107894326B (en) | Splicing main mirror common-phase error detection method based on multi-wavelength phase modulation | |
| CN106768876B (en) | Space solar telescope wavefront sensing methods based on asterism hot spot | |
| Lacour et al. | Reaching micro-arcsecond astrometry with long baseline optical interferometry-Application to the GRAVITY instrument | |
| CN113639676B (en) | Two-dimensional photoelectric auto-collimation method and device based on wavefront measurement and correction | |
| Lee et al. | Fine optical alignment correction of astronomical spectrographs via in-situ full-field moment-based wavefront sensing | |
| CN107797264B (en) | Co-phase adjusting device of synthetic aperture telescope | |
| CN105466576A (en) | Atmospheric turbulence height and angle anisoplanatism wavefront error synchronization measurement device and synchronization measurement method | |
| Acton et al. | Multi-field alignment of the james webb space telescope | |
| CN102589472B (en) | Method for highly precisely eliminating adjustment error in spherical surface shape interference detection | |
| CN111829671A (en) | High-resolution wavefront detection device and wavefront restoration method | |
| CN111399210B (en) | Coarse alignment method, device, equipment and storage medium for large-caliber large-field telescope | |
| Wang et al. | Vision ray model based stereo deflectometry for the measurement of the specular surface | |
| Eves et al. | Phase shifting angle interferometer | |
| Dhabal et al. | Optics alignment of a balloon-borne far-infrared interferometer BETTII | |
| Quirrenbach et al. | The PRIMA astrometric planet search project | |
| Sun et al. | Simultaneous measurement of displacement, pitch, and yaw angles of a cavity output mirror based on phase measuring deflectometry | |
| Krist et al. | End-to-end coronagraphic modeling including a low-order wavefront sensor | |
| Tang et al. | High-precision phase shift method for heavy-load reference mirrors based on nano-precision grating sensor monitoring | |
| CN115046476B (en) | Sparse aperture telescope edge sensing method | |
| Kotani et al. | Long-baseline optical fiber interferometer instruments and science | |
| Li et al. | Imaging performance evaluation and phasing error correction of sparse aperture telescope based on small satellites formation |
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 |