CN112229397B - Satellite angular position intensity correlation measurement system and method based on spatial modulation - Google Patents

Satellite angular position intensity correlation measurement system and method based on spatial modulation Download PDF

Info

Publication number
CN112229397B
CN112229397B CN202010952170.5A CN202010952170A CN112229397B CN 112229397 B CN112229397 B CN 112229397B CN 202010952170 A CN202010952170 A CN 202010952170A CN 112229397 B CN112229397 B CN 112229397B
Authority
CN
China
Prior art keywords
pixel detector
metal hole
hole screen
light
light path
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
Application number
CN202010952170.5A
Other languages
Chinese (zh)
Other versions
CN112229397A (en
Inventor
尹少齐
喻虹
谈志杰
韩申生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Optics and Fine Mechanics of CAS
Original Assignee
Shanghai Institute of Optics and Fine Mechanics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Optics and Fine Mechanics of CAS filed Critical Shanghai Institute of Optics and Fine Mechanics of CAS
Priority to CN202010952170.5A priority Critical patent/CN112229397B/en
Publication of CN112229397A publication Critical patent/CN112229397A/en
Application granted granted Critical
Publication of CN112229397B publication Critical patent/CN112229397B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C1/00Measuring angles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations

Abstract

A system and a method for measuring intensity correlation of star body angular positions based on spatial modulation comprise a pulsar to be measured, a first metal hole screen, a first single-pixel detector, a first fixed platform, a second metal hole screen, a second single-pixel detector, a second fixed platform and a computer, wherein light emitted by a pulsar radiation source passes through the first metal hole screen and the first single-pixel detector to serve as a light path A, and light emitted by the pulsar radiation source passes through the second metal hole screen and the second single-pixel detector to serve as a light path B. The invention relates two light paths, wherein the light path A is fixed and the light path B rotates. When the second-order correlation degree of the two light paths is maximum, the angle information of the two metal hole screens at the moment is recorded, and the field angle information to be measured of the two light paths is finally obtained by utilizing the geometric relation of the star angular position intensity correlation measurement system based on space modulation. The method is based on the high-order correlation characteristic of the light field, and has the advantage of realizing high-precision angle position measurement on the target star to be measured.

Description

Satellite angular position intensity correlation measurement system and method based on spatial modulation
Technical Field
The invention relates to the field of navigation of pulsar, in particular to a satellite angular position intensity correlation measurement system and method based on spatial modulation.
Background
Pulsar is a kind of neutron star, has the characteristic of sending periodic pulses, and can be used for deep space autonomous navigation. In particular to a millisecond pulsar which has a stable period, the frequency range of emitted electromagnetic waves is the range of microwaves and X rays, and the detection limit resolution ratio by utilizing the X rays is high and is several orders of magnitude higher than that of visible light. In consideration of adverse influence factors such as turbulence and cosmic noise in measurement, the second-order intensity correlation measurement method based on statistical optics is adopted for processing, so that the noise influence is expected to be overcome, and high-resolution pulsar observation is realized.
The full width at half maximum (FWHM) of an interference fringe for intensity-correlated interferometry appears as the transverse coherence dimension of the speckle, only within which the speckle can be correlated, with a spatial coherence dimension of
Figure BDA0002677362420000011
For a thermal light source of finite size, the greater the distance from the light source to the scattering medium, the greater the spatial coherence scale; the smaller the source size, the larger the spatial coherence scale. The X-rays emitted by the millisecond pulsar have a shorter coherence time and energy resolution
Figure BDA0002677362420000012
The time resolution of the detector is estimated to be ps or fs magnitude, while the time resolution of the existing detector can reach ns magnitude and does not meet the requirement.
In prior art 1 (yangting high, gaoyiping, a universal pulsar autonomous navigation measurement model construction method, CN104316048B), pulsar navigation is realized by giving a theoretical relationship between the time when a pulsar pulse reaches an aircraft and the pulsar emission time.
Prior art 2 (zhao shenmei, liang weng, cheng super, dong xiao liang, a quantum correlation imaging method based on angular position entanglement, CN104407485A) is implemented by encoding imaged object information into angular slits of different angles and loading the different angular slits onto a light beam. Because a monotonous functional relation exists between the coincidence count and the angle seam angle, the information of the object can be obtained on the reference light path according to the coincidence count value, and the imaging of the object is obtained.
Prior art 3 (chenjun, a material wave correlation imaging generation method and apparatus thereof, CN103412304A) obtains an image of an object by performing bucket measurement on a signal material wave after the object is projected and performing spatial correlation with the reference material wave.
Prior art 4 (chen zhipeng, quantum imaging method and quantum imaging system, CN102087411A), obtains an image of the object by performing multiple correlation measurements on the object and performing signal processing on the data of the correlation measurements by using a compressive sensing algorithm.
In the prior art 5 (Zhanhua, Shiga, Liu Meng, Guo Yi pan, ChaoKai, Chengming, Wei Shi Yong, X-ray pulsar space navigation ground simulation positioning system and positioning method, CN110068339A), the coordinate position of the X-ray radiation source is calculated according to the position relationship between the X-ray radiation source and the detector in the coordinate system.
In prior art 6 (chenxihao, montouying, boring, wuweiqi, paying, Zhang, Xiuyiemin, super-resolution associated imaging system and imaging method based on low-pass filtering, CN107219638A), according to a threshold, an area array data signal obtained from each time sequence point is sequentially processed through a reference arm spatial filter or a detection arm spatial filter for low-pass filtering operation, and finally the two groups of area array data obtained from the reference arm optical path and the detection arm optical path are processed according to the traditional thermo-optical associated imaging principle and method to realize super-resolution associated imaging of an object to be imaged.
In prior art 7 (yangting high, gaoyiping, pulsar-based navigation constellation time synchronization and directional parameter measurement method, CN104316056B), in case of applying inter-satellite link navigation technology, the pulsar space-time reference frame is used as a reference standard to measure the navigation constellation overall rotation and realize the navigation constellation time synchronization. A single-station observation method is provided by selecting a navigation satellite to observe pulsar; the same pulsar is observed by selecting two navigation satellites, so that a poor observation method is provided.
Prior art 8 (surname bin, a two-dimensional compression ghost imaging system and method based on coincidence measurement, CN103323396B) and prior art 9 (shuncky, liuxuefeng, surname bin, zhazhijie, an entanglement imaging system and method based on dual-compression coincidence measurement, CN103308189A), coincidence measurement is performed on the total light intensity of the object arm light path and the reference arm light path through a coincidence measurement circuit, a coincidence measurement value is output, and an algorithm module reconstructs spatial correlation coefficient distribution by using a compression sensing algorithm according to a measurement matrix and the measurement value.
In summary, the above prior art relates to the improvement of the imaging quality of an object, or the calculation of the coordinate position of a radiation source through the position relationship between the radiation source and a detector in a coordinate system. The invention aims to realize the detection of the field angles to be detected of the two light paths, obtains the second-order correlation degree of the two light paths through correlation calculation, does not relate to object information in the light paths, searches the angle information of the two metal hole screens corresponding to the maximum second-order correlation degree, and completes the detection of the field angles to be detected of the two light paths by utilizing the geometric relation of a star body angular position intensity correlation measurement system based on spatial modulation. From the aspect of the method, the invention provides a new scheme for the high-precision angular position measurement of the pulsar, the scheme comprises field angles to be measured of two light paths, angle information of two metal hole screens and position information of two single-pixel detectors, and the method has universal applicability to the measurement of different field angles to be measured of a light source. According to the scheme, the direct detection of the angle to be detected is converted into the searching of the angle information of the two metal hole screens corresponding to the two optical paths when the second-order correlation degree of the two optical paths is the maximum, so that the detection of the field angle to be detected of the two optical paths is realized, the measurement precision of the field angle to be detected is improved, and the detection precision of the field angle to be detected of the two optical paths is improved by controlling the rotation precision of the two metal hole screens.
Disclosure of Invention
The technical problem to be solved by the present invention is to overcome the defects of the prior art, and provide a system and a method for measuring intensity correlation of star angular positions based on spatial modulation, wherein a light intensity sequence acquired by two optical paths is subjected to correlation operation by fixing the first metal aperture screen and the first single-pixel detector of the optical path a, rotating the second metal aperture screen and the second single-pixel detector of the optical path B, so as to obtain angle information of the first metal aperture screen and the second metal aperture screen corresponding to the second-order correlation peak value when the second-order correlation peak value is maximum, and finally obtain the field angles to be measured of the two optical paths.
The technical solution of the invention is as follows:
a star angular position intensity correlation measurement system based on spatial modulation comprises a first metal hole screen, a first single-pixel detector, a first fixed platform, a second metal hole screen, a second single-pixel detector, a second fixed platform and a computer, wherein light emitted by a pulsar to be measured passes through a first goldThe aperture screen and the first single-pixel detector are used as a light path A, and light emitted by the pulsar to be detected passes through the second metal aperture screen and the second single-pixel detector to be used as a light path B. The distance from the first metal hole screen to the first single-pixel detector is equal to the distance from the second metal hole screen to the second single-pixel detector, and the distance is d. First metal hole screen with first single pixel detector all fix first fixed platform on, second metal hole screen with second single pixel detector all fix second fixed platform on, second fixed platform have the rotation function. The first metal hole screen and the second metal hole screen have the same spatial structure distribution, M 2 N 2 Is the transverse symmetry axis of the second metal hole screen. The input end of the computer is connected with the output ends of the first single-pixel detector and the second single-pixel detector, and the computer is provided with a program for performing correlation operation on the collected light intensity sequence.
The method for measuring the intensity correlation of the angular position of the star body by utilizing the intensity correlation measurement of the angular position of the star body based on the spatial modulation is characterized by comprising the following steps of:
<1> selecting said single pixel detector, satisfying the following formula:
Figure BDA0002677362420000041
in the formula, the transverse diameter of the single-pixel detector is l, lambda is the wavelength of the ray emitted by the pulsar to be detected, d is the distance from the metal hole screen to the single-pixel detector, and a is the transverse length of the metal hole screen;
<2>the first single-pixel detector (3) and the second single-pixel detector (5) are exposed for k times at the same time to respectively obtain the initial angle theta 0 Initial light intensity sequence of lower light path A
Figure BDA0002677362420000042
And light path B initial light intensity sequence
Figure BDA0002677362420000043
And the light path B is started to generate a light intensity sequence
Figure BDA0002677362420000044
And light path A initial light intensity sequence
Figure BDA0002677362420000045
Performing correlation operation to obtain intensity correlation distribution in the correlated imaging sequence under the initial angle
Figure BDA0002677362420000046
<3>Calculating the m times, wherein m is 0, N is the total times of rotation, and the light intensity sequence of the light path B
Figure BDA0002677362420000047
Average value of (2)
Figure BDA0002677362420000048
Calculating the light intensity sequence of the light path A
Figure BDA0002677362420000049
Average value of (2)
Figure BDA00026773624200000410
Calculating the average value of the intensity correlation distribution of the light intensity sequence A and the light intensity sequence B
Figure BDA0002677362420000051
<4> calculating the second-order strength correlation value, the formula is as follows:
Figure BDA0002677362420000052
<5>winding the second fixed platform (8) around the transverse central axis M of the second metal hole screen (4) 2 N 2 After rotating by the next angle thetaThe first single-pixel detector (3) and the second single-pixel detector (5) are exposed for k times, and after the mth rotation is recorded, theta is recorded m =θ 0 + m θ, m-th light intensity sequence of optical path A
Figure BDA0002677362420000053
And light path B mth light intensity sequence
Figure BDA0002677362420000054
And the light intensity sequence of the mth time of the light path B
Figure BDA0002677362420000055
And light path A mth light intensity sequence
Figure BDA0002677362420000056
Performing correlation operation to obtain intensity correlation distribution in the correlated imaging sequence under the mth angle
Figure BDA0002677362420000057
<6>Repeating the steps<3>、<4>、<5>Obtaining a plurality of groups of second-order intensity correlation values until the peak value of the second-order intensity correlation values appears, and recording the angle xi of the first metal aperture screen under a Cartesian coordinate system at the moment 1 And angle xi of the second metal aperture screen in a Cartesian coordinate system 2
And (7) through the geometric relationship of the star angular position intensity correlation measurement system based on spatial modulation, the to-be-measured field angles of the two light paths meet the following relationship:
Δξ=|ξ 12 |
wherein, Δ ξ represents the field angle of the two beams of light emitted by the pulsar to be measured.
Compared with the prior art, the invention has the following technical effects:
1) the invention relates two light paths, wherein the light path A is fixed and the light path B rotates. When the second-order correlation degree of the two light paths is maximum, the angle information of the two metal hole screens at the moment is recorded, and the field angle information to be measured of the two light paths is finally obtained by utilizing the geometric relation of the star angular position intensity correlation measurement system based on space modulation.
2) The device has universal applicability to measurement of different field angles of the light source, and has the characteristics of reasonable structural design, few components and convenience in implementation.
3) The direct detection of the angle is changed into the search of the angle of the two metal hole screens corresponding to the maximum second-order correlation degree, so that the detection of the field angle to be detected of the two light paths is realized, and the measurement precision of the field angle to be detected is improved.
Drawings
FIG. 1 is a schematic diagram of the optical path of the star angular position intensity correlation measurement system based on spatial modulation, in which:
1: pulsar to be measured, 2: first metal aperture screen, 3: first single-pixel detector, 4: second metal aperture screen, 5: second single-pixel detector, 6: and (4) a computer.
Fig. 2 is a positional relationship of a single-pixel detector with respect to a metal aperture screen. Wherein (a) is the position relationship between the first metal aperture screen and the first single-pixel detector, M 1 N 1 Is the transverse symmetry axis of the first metal aperture screen, O 1 The first metal hole screen and the first single-pixel detector are both arranged on the first fixed platform and are projection points of the first single-pixel detector relative to the first metal hole screen. (b) For the positional relationship of the second metal aperture screen and the second single-pixel detector, M 2 N 2 Is the transverse symmetry axis of the second metal aperture screen, O 2 For the projection point of second single pixel detector for second metal hole screen, second metal hole screen and second single pixel detector all settle on second fixed platform, and second fixed platform possesses rotation function, in the picture:
2: first metal aperture screen, 3: first single-pixel detector, 4: second metal aperture screen, 5: second single-pixel detector, 7: first fixed platform, 8: second fixed platform
FIG. 3 is a diagram of a simulated light field distribution of a light field at a metal aperture screen that has freely propagated to a plane in which a single pixel detector is locatedThe positions of (a) are outlined in black boxes in the figure. Wherein (a) is the condition that the relative angle between the included angle of the optical path A and the first metal hole screen and the included angle of the optical path B and the second metal hole screen is 0 degree, and (B) is the relative angle between the included angle of the optical path A and the first metal hole screen and the included angle of the optical path B and the second metal hole screen is (5 multiplied by 10) -6 ) DEG, in the case of a magnetic field.
Fig. 4 is a correlation diagram for two-way light based on fig. 3. Wherein (a) is the condition that the relative angle between the included angle of the optical path A and the first metal hole screen and the included angle of the optical path B and the second metal hole screen is 0 degree, and (B) is the relative angle between the included angle of the optical path A and the first metal hole screen and the included angle of the optical path B and the second metal hole screen is (5 multiplied by 10) -6 ) DEG, in the case of a magnetic field.
Fig. 5 is a schematic geometric relationship diagram of the system and method for measuring intensity correlation of star angular positions based on spatial modulation according to the present invention, which omits the positions of two single-pixel detectors,
in the figure: o: pulsar to be measured, EF: first metal aperture screen, CD: a second metal aperture screen.
Detailed Description
For better understanding of the objects, technical solutions and advantages of the present invention, the following description of the present invention with reference to the accompanying drawings is provided for further description, but should not be construed to limit the scope of the present invention.
The invention relates to a star angular position intensity correlation measurement system based on spatial modulation, which comprises a pulsar 1 to be measured, a first metal hole screen 2, a first single-pixel detector 3, a second metal hole screen 4, a second single-pixel detector 5 and a computer 6, as shown in figure 1. The distance from the first metal hole screen 2 to the first single-pixel detector 3 is equal to the distance from the second metal hole screen 4 to the second single-pixel detector 5. The spatial structure distribution of the first metal hole screen 2 and the second metal hole screen 4 is the same; the light intensity information of the first metal hole screen 2 is received and recorded by the first single-pixel detector 3, and the light intensity information of the second metal hole screen 4 is received and recorded by the second single-pixel detector 5. The computer 6 is connected with the output ends of the two single-pixel detectors and is provided with a program for performing correlation operation on the collected light intensity sequences.
The position relationship between the metal aperture screen and the single-pixel detector is shown in FIG. 2, wherein (a) is the position relationship between the first metal aperture screen and the first single-pixel detector, and M is 1 N 1 Is the transverse symmetry axis of the first metal aperture screen, O 1 The first metal hole screen and the first single-pixel detector are both fixed on the first fixing platform, and the first single-pixel detector is a projection point of the first single-pixel detector relative to the first metal hole screen. (b) For the positional relationship of the second metal aperture screen and the second single-pixel detector, M 2 N 2 Is the transverse symmetry axis of the second metal aperture screen, O 2 The projection point of the second single-pixel detector is relative to the second metal hole screen, the second metal hole screen and the second single-pixel detector are both fixed on a second fixed platform, and the second fixed platform has a rotating function;
the method for measuring the intensity correlation of the angular position of the star body based on the spatial modulation, which is used for carrying out experiments on X-rays, comprises the following steps:
<1> selecting a single pixel detector, satisfying the following formula:
Figure BDA0002677362420000081
in the formula, l is the transverse diameter of the single-pixel detector, lambda is the wavelength of a ray emitted by the pulsar (1) to be detected, d is the distance from the metal hole screen to the single-pixel detector, and a is the transverse length of the metal hole screen;
in this embodiment, the distance between the first single-pixel detector 3 and the first metal aperture screen 2 and the distance d between the second single-pixel detector 5 and the second metal aperture screen 4 are 10 meters; the transverse length a of the first metal hole screen 2 and the second metal hole screen 4 is 1 mm;
<2>the first single-pixel detector 3 and the second single-pixel detector 5 are exposed k times at the same time to respectively obtain the initial angle theta 0 Initial light intensity sequence of lower light path A
Figure BDA0002677362420000082
And light path B initial light intensity sequence
Figure BDA0002677362420000083
And the light path B is started to generate a light intensity sequence
Figure BDA0002677362420000084
And light path A initial light intensity sequence
Figure BDA0002677362420000085
Performing correlation operation to obtain intensity correlation distribution in the correlated imaging sequence under the initial angle
Figure BDA0002677362420000086
<3>Calculating the m times, wherein m is 0, N is the total times of rotation, and the light intensity sequence of the light path B
Figure BDA0002677362420000087
Average value of (2)
Figure BDA0002677362420000088
Calculating the light intensity sequence of the light path A
Figure BDA0002677362420000089
Average value of (2)
Figure BDA00026773624200000810
Calculating the average value of the intensity correlation distribution of the light intensity sequence of the light path A and the light intensity sequence of the light path B
Figure BDA00026773624200000811
<4> calculating the second-order intensity correlation value, the formula is as follows:
Figure BDA00026773624200000812
<5>shielding the second stationary platform 8 around the second metal aperture4 transverse central axis M 2 N 2 After rotating a certain angle theta, exposing the first single-pixel detector 3 and the second single-pixel detector 5 for k times, and respectively recording theta after the mth rotation m =θ 0 + m θ, m-th light intensity sequence of optical path A
Figure BDA0002677362420000091
And light path B mth light intensity sequence
Figure BDA0002677362420000092
And the light intensity sequence of the mth time of the light path B
Figure BDA0002677362420000093
And light path A mth light intensity sequence
Figure BDA0002677362420000094
Performing correlation operation to obtain intensity correlation distribution in the correlated imaging sequence under the mth angle
Figure BDA0002677362420000095
<6>Repeating the steps<3>、<4>、<5>And obtaining a plurality of groups of second-order intensity correlation values until the peak value of the second-order intensity correlation values appears. In the present embodiment, the angle θ corresponding to the peak of the second-order intensity correlation value N At 0 deg. relative angle theta m Ranging from- (1X 10) -5 ) At an angle of 0 DEG, per rotation (1X 10) -7 ) Degree. The left diagram in FIG. 3 is the relative angle theta between the included angle between the light path A and the first metal hole screen 2 and the included angle between the light path B and the second metal hole screen 4 m Is- (5X 10) -6 ) In the case of DEG, the right diagram is a relative angle theta between an included angle between the light path A and the first metal hole screen 2 and an included angle between the light path B and the second metal hole screen 4 m The 0 deg. case, where the single pixel detector is located, is outlined in black in the figure. Fig. 4 is a correlation diagram of two optical fields obtained based on fig. 3, and the central peak point of the right diagram is the position of the single-pixel detector relative to the metal aperture screen. Moving over at angular positionsIn the process, the position of the single-pixel detector is unchanged, and the second-order correlation peak point is not in the same position as the right image as seen from the left image, so that the correlation light intensity values of the two light paths are different at different angular positions. When the second-order intensity correlation value is recorded to be at the peak position, the angle xi of the first metal hole screen 2 in a Cartesian coordinate system 1 And the angle xi of the second metal aperture screen 4 in a Cartesian coordinate system 2
<7> the size of the field angle to be measured of the two optical paths is determined by a geometric relational graph of the star angular position intensity correlation measurement system based on the spatial modulation as shown in figure 5. In a geometric relational graph of a star angular position intensity correlation measurement system based on spatial modulation, auxiliary lines C ' D ' and O ' A are made, and the following relations are satisfied:
C'D'//CD,O'A//OB
the included angle between the light path A and the first metal hole screen 2 and the auxiliary line C ' D ' is OAD '. When the relative angle theta of the first metal hole screen 2 and the second metal hole screen 4 m When the angle is 0 degree, the included angle between the optical path A of the pulsar to be measured and the first metal hole screen 2 and the included angle between the optical path B of the pulsar to be measured and the second metal hole screen 4 are equal to be alpha, and the angle information of the first metal hole screen 2 under the Cartesian coordinate system is xi 1 And the angle information of the second metal hole screen 4 in the Cartesian coordinate system is ξ 2 The following relationship is satisfied:
α=γ+Δξ,Δξ=|ξ 12 |
according to the geometrical relationship, the opening angles to be measured of the two optical paths satisfy the following relationship:
∠AOB=∠O'AO=Δξ,Δξ=|ξ 12 |
therefore, the flare angle [ AOB ] of two beams of light emitted by the pulsar to be detected is equal to the angle [ xi ] of the first metal hole screen 2 in a Cartesian coordinate system 1 And the angle xi of the second metal hole screen 4 in a Cartesian coordinate system 2 The difference of (a).

Claims (2)

1. A star angular position intensity correlation measurement system based on spatial modulation is characterized by comprising a first metal hole screen (2), a first single-pixel detector (3), a second metal hole screen (4), a second single-pixel detector (5), a computer (6), a first fixed platform (7) and a second fixed platform (8);
light emitted by the pulsar (1) to be detected passes through the first metal hole screen (2) and the first single-pixel detector (3) to be used as a light path A;
light emitted by the pulsar (1) to be detected passes through the second metal hole screen (4) and the second single-pixel detector (5) to be used as a light path B;
the distance from the first metal hole screen (2) to the first single-pixel detector (3) is equal to the distance from the second metal hole screen (4) to the second single-pixel detector (5), and the distance is d; the first metal hole screen (2) and the first single-pixel detector (3) are fixed on the first fixing platform (7), and the second metal hole screen (4) and the second single-pixel detector (5) are fixed on the second fixing platform (8);
the first metal hole screen (2) and the second metal hole screen (4) have the same spatial distribution;
the transverse central axis of the second metal hole screen (4) is M 2 N 2
The input end of the computer (6) is connected with the output ends of the first single-pixel detector (3) and the second single-pixel detector (5), the computer (6) is provided with a program for performing correlation operation on the acquired light intensity sequence, and the method specifically comprises the following steps:
<1> selecting a single pixel detector, satisfying the following formula:
Figure FDA0003708412410000011
in the formula, l is the transverse diameter of the single-pixel detector, lambda is the wavelength of a ray emitted by the pulsar (1) to be detected, d is the distance from the metal hole screen to the single-pixel detector, and a is the transverse length of the metal hole screen;
adjusting the first metal hole screen (2) and the first single-pixel detector (3) to be coaxial with the pulsar (1) to be detected, adjusting the second metal hole screen (4) and the second single-pixel detector (5) to be coaxial with the pulsar (1) to be detected, wherein the first metal hole screen (2), the first single-pixel detector (3) and the second metal hole screen (4) and the second single-pixel detector (5) in the light path A are not coaxial;
<2>the first single-pixel detector (3) and the second single-pixel detector (5) are exposed for k times at the same time to respectively obtain the initial angle theta 0 Initial light intensity sequence of lower light path A
Figure FDA0003708412410000022
And initial light intensity sequence of light path B
Figure FDA0003708412410000023
And the light path B is started to generate a light intensity sequence
Figure FDA0003708412410000024
And light path A initial light intensity sequence
Figure FDA0003708412410000025
Performing correlation operation to obtain intensity correlation distribution in the correlated imaging sequence under the initial angle
Figure FDA0003708412410000026
<3>Calculating the m times, wherein m is 0, N is the total times of rotation, and the light intensity sequence of the light path B
Figure FDA0003708412410000027
Average value of (2)
Figure FDA0003708412410000028
Calculating the light intensity sequence of the light path A
Figure FDA0003708412410000029
Average value of (2)
Figure FDA00037084124100000210
Calculating the average value of the intensity correlation distribution of the light intensity sequence of the light path A and the light intensity sequence of the light path B
Figure FDA00037084124100000211
<4> calculating the second-order intensity correlation value, the formula is as follows:
Figure FDA0003708412410000021
<5>winding the second fixed platform (8) around the transverse central axis M of the second metal hole screen (4) 2 N 2 After rotating the next angle theta, k times of exposure is carried out on the first single-pixel detector (3) and the second single-pixel detector (5), and after the mth rotation is recorded respectively, the theta is recorded m =θ 0 + m θ, m-th light intensity sequence of optical path A
Figure FDA00037084124100000212
And light path B mth light intensity sequence
Figure FDA00037084124100000213
And the light path B is subjected to the mth light intensity sequence
Figure FDA00037084124100000214
And light path A mth light intensity sequence
Figure FDA00037084124100000215
Performing correlation operation to obtain intensity correlation distribution in the correlated imaging sequence under the mth angle
Figure FDA00037084124100000216
<6>Repeating the steps<3>、<4>、<5>Obtaining a plurality of groups of second-order strength correlation values until the peak value of the second-order strength correlation values appears, and recording the angle xi of the first metal hole screen under the Cartesian coordinate system at the moment 1 And angle xi of the second metal aperture screen in a Cartesian coordinate system 2
And (7) through the geometric relationship of the star angular position intensity correlation measurement system based on spatial modulation, the to-be-measured field angles of the two light paths meet the following relationship:
Δξ=|ξ 12 |
wherein, Δ ξ represents the field angle of the two beams of light emitted by the pulsar to be measured.
2. The system according to claim 1, wherein the pulsar (1) to be measured emits X-rays.
CN202010952170.5A 2020-09-11 2020-09-11 Satellite angular position intensity correlation measurement system and method based on spatial modulation Active CN112229397B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010952170.5A CN112229397B (en) 2020-09-11 2020-09-11 Satellite angular position intensity correlation measurement system and method based on spatial modulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010952170.5A CN112229397B (en) 2020-09-11 2020-09-11 Satellite angular position intensity correlation measurement system and method based on spatial modulation

Publications (2)

Publication Number Publication Date
CN112229397A CN112229397A (en) 2021-01-15
CN112229397B true CN112229397B (en) 2022-08-30

Family

ID=74116155

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010952170.5A Active CN112229397B (en) 2020-09-11 2020-09-11 Satellite angular position intensity correlation measurement system and method based on spatial modulation

Country Status (1)

Country Link
CN (1) CN112229397B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102087411A (en) * 2010-12-02 2011-06-08 上海电机学院 Quantum imaging method and quantum imaging system
CN103308189A (en) * 2013-06-28 2013-09-18 中国科学院空间科学与应用研究中心 Entanglement imaging system and method based on dual-compression coincidence measurements
CN206192364U (en) * 2016-11-09 2017-05-24 中国科学院西安光学精密机械研究所 Device is verified on VLBI measurement system and ground based on X ray
CN108827988A (en) * 2018-06-12 2018-11-16 中国科学院上海光学精密机械研究所 Non-local modulation of X-ray diffraction imaging device and method based on light field High order correletion
CN110686702A (en) * 2019-09-27 2020-01-14 中国科学院上海光学精密机械研究所 Pulsar photon flow simulation system and method based on light field high-order correlation
CN110793528A (en) * 2019-09-27 2020-02-14 西安空间无线电技术研究所 Low-orbit satellite-based anchoring-based Beidou navigation constellation autonomous orbit determination method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102087411A (en) * 2010-12-02 2011-06-08 上海电机学院 Quantum imaging method and quantum imaging system
CN103308189A (en) * 2013-06-28 2013-09-18 中国科学院空间科学与应用研究中心 Entanglement imaging system and method based on dual-compression coincidence measurements
CN206192364U (en) * 2016-11-09 2017-05-24 中国科学院西安光学精密机械研究所 Device is verified on VLBI measurement system and ground based on X ray
CN108827988A (en) * 2018-06-12 2018-11-16 中国科学院上海光学精密机械研究所 Non-local modulation of X-ray diffraction imaging device and method based on light field High order correletion
CN110686702A (en) * 2019-09-27 2020-01-14 中国科学院上海光学精密机械研究所 Pulsar photon flow simulation system and method based on light field high-order correlation
CN110793528A (en) * 2019-09-27 2020-02-14 西安空间无线电技术研究所 Low-orbit satellite-based anchoring-based Beidou navigation constellation autonomous orbit determination method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Angular position measurement of pulsars based on X-ray intensity correlation;Su T , et al;《 Optik》;20181231;第8-11页 *
X 射线强度关联干涉测量能谱展宽校正;杨善初等;《光学学报》;20191031;第39卷(第10期);第1-8页 *

Also Published As

Publication number Publication date
CN112229397A (en) 2021-01-15

Similar Documents

Publication Publication Date Title
Maggio et al. How does a dark compact object ringdown?
KR101666688B1 (en) Device for nondestructively examining composite structure and nondestructive examination method
Bérujon et al. Two-dimensional x-ray beam phase sensing
Berujon et al. X-ray pulse wavefront metrology using speckle tracking
CN107862661A (en) A kind of optical coherence tomography system method for correcting image
JP3449791B2 (en) Imaging method and imaging apparatus
Namikawa et al. Detecting black-hole binary clustering via the second-generation gravitational-wave detectors
CN102809494A (en) Knife-edge measuring method of modulation transfer function of digital X-ray imaging system
JP3205160B2 (en) Magnetic field measuring method and charged particle beam device using the same
Mihailescu et al. Standoff 3D gamma-ray imaging
CN112229397B (en) Satellite angular position intensity correlation measurement system and method based on spatial modulation
JP2019502900A5 (en)
Craciunescu et al. The maximum likelihood reconstruction method for JET neutron tomography
KR20080089404A (en) Characterization of three-dimensional distribution of defects by x-ray topography
US11604152B2 (en) Fast industrial computed tomography for large objects
CN111025202B (en) Scanning type three-dimensional magnetic field detection method and device
CN112325874B (en) System and method for measuring strength correlation of star angular position based on snake-shaped light path
Mooley et al. Optical measurement of superluminal motion in the neutron-star merger GW170817
US20230252663A1 (en) X-ray digital image correlation
JP7437337B2 (en) Internal state imaging device and internal state imaging method
Mokhnachevskaya et al. Differential Cherenkov detectors: experimental technique and results
Wilke et al. Two-dimensional time-dependent imaging utilizing tomographic concepts
Jiang et al. Microwave Staring Correlated Imaging Based on Quasi-Stationary Platform with Motion Measurement Errors
JP3380921B2 (en) How to measure strain in crystals
CN115436944A (en) Radar imaging focusing method and device, computing equipment and storage medium

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