CN110686702A - Pulsar photon flow simulation system and method based on light field high-order correlation - Google Patents
Pulsar photon flow simulation system and method based on light field high-order correlation Download PDFInfo
- Publication number
- CN110686702A CN110686702A CN201910924850.3A CN201910924850A CN110686702A CN 110686702 A CN110686702 A CN 110686702A CN 201910924850 A CN201910924850 A CN 201910924850A CN 110686702 A CN110686702 A CN 110686702A
- Authority
- CN
- China
- Prior art keywords
- time
- photon
- detector
- sequence
- time difference
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- Manufacturing & Machinery (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention provides a pulsar photon flow simulation system and method based on light field high-order correlation, wherein the system comprises: the photon flow simulation subsystem is used for generating a photon time sequence, normalizing the light intensity distribution of the detection surface as space probability distribution, dividing detection surface regions, summing region probabilities to obtain the probability of photon collapse to each region, sequentially collapsing photon flow according to the region probabilities, marking photons collapsed to the detector region, extracting the time of marked photons from the photon time sequence, and outputting two paths of detector time sequences to the coincidence algorithm subsystem; and the coincidence algorithm subsystem is used for identifying the coincidence photon pairs, performing algorithm coincidence on the time sequences of the two paths of detectors, recording the time difference of the coincidence photon pairs and outputting a time difference distribution curve. The simulation system can obtain the time difference distribution curve consistent with the thermo-optic coincidence measurement result, and realizes the simulation of the pulsar X-ray intensity correlation interferometry system based on coincidence counting.
Description
Technical Field
The invention relates to the field of pulsar observation and navigation, in particular to a pulsar photon flow simulation system and method based on light field high-order association, and particularly relates to a photon flow simulation and data conformity processing method.
Background
Pulsar is a neutron star rotating at high speed and periodically emits a pulse signal, called "lighthouse" in the universe. The high-precision measurement of the pulsar information is beneficial to establishing a satellite navigation network and realizing the high-precision autonomous navigation of the spacecraft. In the existing pulsar astronomical observation means, the observation of a radio array is limited by wavelength, and the resolution can reach 10 milli-angular seconds at most; the visible light wave band observation is limited by the wavelength and the caliber of the detector, and the resolution ratio is difficult to realize the improvement of the order of magnitude precision. Considering that the wavelength of X-rays is several orders of magnitude shorter than that of visible light, observing pulsar in an X-ray band is expected to realize high-precision measurement in the micro-angle second order.
The interferometry is an important method for astronomical observation, the conventional interferometry acquires information by utilizing the first-order correlation characteristic of a light field, and for an X-ray system, the processing precision of an optical device is required to reach the nanometer level, which has great difficulty in engineering realization. The pulsar radiation field is considered to be a thermal light field, has random fluctuation and obeys Gaussian statistics, and the information is acquired by utilizing the high-order correlation of pulsar X-rays, so that the method is a feasible way for realizing high-resolution observation of pulsars.
Pulsar is more than several thousand years of light from earth, and the flux of light reaching the satellite detector is extremely low (e.g., the brightest Crab pulsar, which also has a flux of only 1.54photon/cm2S) can only be detected with single photons. The existing pulsar simulation system is used for simulating and measuring the relative distance between satellites, determining the phase information of pulsar and the like by recording the photon time sequence reaching a detector, thereby being beneficial to realizing high-precision pulsar navigation. However, the existing pulsar simulation system does not consider the high-order correlation characteristic of the pulsar X-ray in photon flow simulation and data signal processing, and cannot be applied to the development of a pulsar interferometry system based on the high-order correlation.
In prior art 1 (friendship, commander, a high-precision X-ray pulsar signal simulation method, CN104567937A) and prior art 2 (huijiejun, shoufei, an X-ray pulsar signal simulation source based on a random single photon emission mechanism, CN103697908A), signal simulation of an X-ray pulsar is performed based on a single photon emission mechanism, after initial photon arrival time is determined, photon phase is calculated according to a pulsar phase prediction model, then an inverse function method is used to recur to obtain the next photon phase, and the next photon arrival time can be solved by substituting the phase into the phase prediction model. The process in which photons arrive at the detector and are recorded is the poisson process, and the photon arrival time is a random number.
In prior art 3 (zhanhua, scholanzi, a frequency domain weighted phase comparison method for X-ray pulsar photon sequences, CN105300386A), two arrays of photon arrival time sequences received by two spacecrafts are sampled at equal intervals to obtain photon intensity sequences, transformed into the frequency domain by FFT to extract phase information to obtain phase differences, and the ratio of the phase differences with respect to frequency is weighted by energy to obtain normalized delay, so that the estimation accuracy of the phase can be effectively improved.
In the prior art 4 (grandson peak, Fanghai swallow, X-ray pulse two-star photon sequence simulation method, CN108981750A) and the prior art 5 (grandson peak, Fanghai swallow, on-orbit X-ray pulsar timing model construction method, CN107144274A), the arrival time of X-ray pulsar photons is preprocessed, and when the fixed time of the near-earth orbit spacecraft is converted into the coordinate of a solar system centroid coordinate system, the response delay of an instrument, Einstein delay of relativistic effect, Shapiro delay of space-time bending caused by a large-mass celestial body, clock correction, year-round parallax correction and the like are considered, so that the estimation precision of the photon arrival time is improved.
In prior art 6 (zhanhua, xu tie, a constellation orientation simulation system and method based on X-ray pulsar, CN103674020B), the significance of photon arrival time to X-ray pulsar navigation is elucidated, and by the optical travel time difference of X-ray pulsar signals propagated between satellites, the relative distance between satellites can be measured, and the included angle between the satellite baselines in the constellation can be obtained, so as to realize satellite navigation and formation satellite orientation.
In prior art 7 (liu jian macro, a cloud start, a coincidence measurement system and method, CN105450215A), pulse shaping is performed first, and then a carry chain latch module is used to perform carry chain latch and encoding on multiple pulse signals; correcting errors of the multi-path codes, decoding and converting into time information; and calculating the relative time value, comparing the relative time value with the coincidence gate width, and judging whether the coincidence is successful.
In summary, in the prior art 1 and the prior art 2, the photon time sequence of the X-ray pulsar is simulated by simulation and experiment means, and in the prior art 3, phase difference information is extracted from the two photon time sequences, but the prior art does not utilize the high-order correlation characteristic of the light field to simulate the photon time sequence, does not consider the wave optics theory, and does not relate to the quantum theory of photon probability collapse. The prior art 4, 5, 6 is concerned with on-track observation correction and does not involve photon flow generation. The prior art 7 is mainly to improve the coincidence measurement hardware system to improve the coincidence accuracy, and does not relate to a data coincidence processing method. Therefore, the prior art does not consider the high-order correlation characteristic of the pulsar as a thermal light source, and the intensity-correlated interferometry simulation cannot be performed.
According to the quantum optical principle, the thermal light source has a beam bunching effect, when the radiant luminous flux is extremely low and the coherence time is extremely short, the intensity correlation between two points (two detector positions) in space can be directly obtained by adopting intensity correlation measurement (also called coincidence measurement) based on coincidence counting, and second-order interference fringes are obtained by scanning the detector on the space position, so that the angle information of the pulsar is obtained. Therefore, the core of the pulsar X-ray interferometry simulation based on light field high-order correlation lies in generating photon flow with thermo-optic bunching effect, and realizing second-order interferometry through data coincidence processing.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a pulsar photon flow simulation system and method based on light field high-order association, which are not only suitable for pulsar X-ray sources, but also suitable for other thermal light sources and have universality in different wave bands. The data accord with the processing algorithm, help to utilize the satellite carried single photon detector to carry on the pulsar interference measurement.
In order to achieve the above object, the technical solution of the present invention is as follows:
a pulsar photon flow simulation system and method based on light field high-order association, the simulation system includes photon flow simulation subsystem, accords with two parts of the algorithm subsystem; the photon flow simulation subsystem comprises a photon time collapse module, a photon space collapse module and a detector time distribution module; the photon time collapse module comprises a time probability submodule and a time collapse submodule; the photon space collapse module comprises a space light intensity sub-module, a space probability sub-module and a region collapse sub-module; the coincidence algorithm subsystem comprises a photon pair time difference searching module and a time difference channel counting module;
the time probability submodule generates time probability distribution which is used as the probability of collapse of each time point of the photon on a time axis and outputs the time probability distribution to the time collapse submodule;
the time collapse submodule collapses photon flow to a time axis according to time probability distribution, and after all photon collapse is finished, the photon collapse submodule is numbered again according to the time sequence to generate a photon time sequence and outputs the photon time sequence to the space probability submodule and the detector time distribution module;
the space light intensity submodule is given different circular complex Gaussian random initial phases aiming at the Gaussian statistical property of the thermo-optic by a transfer function method of classical wave optics to obtain different detection surface light intensity distributions, and the detection surface light intensity distributions are output to the space probability submodule;
the space probability submodule is used for normalizing the light intensity distribution of the detection surface to obtain space probability distribution which is used as the probability of collapse of each position of the photon on the detection surface, calculating the time difference of adjacent photons by utilizing a photon time sequence, distributing a space probability distribution to each photon according to the time difference condition, and outputting the space probability distribution of the photon flow to the area collapse submodule;
the region collapse submodule divides regions where photons are collapsed on a detection surface into a first detector aperture region, a second detector aperture region, …, an Nth detector aperture region and an outer detector region, wherein N is the number of detectors on the detection surface, the probabilities of coverage ranges of the divided regions are summed in spatial probability distribution respectively to obtain the region probabilities of the photons collapsed in the regions, photon flows are collapsed sequentially according to the region probabilities, photons collapsed to the first detector aperture, the second detector aperture, …, the Nth detector aperture and the outer detector are respectively marked as 1, 2, …, N and 0, and a photon marking sequence is output to the detector time distribution module;
the detector time allocation module allocates photon time according to the photon time sequence and the photon marking sequence, selects two detectors from the first detector to the Nth detector, extracts the photon time marked as one of the detectors by using the photon time sequence, arranges the photon time from small to large, and generates the time sequence { A of the detector 11，A2，A3，A4…, extracting the photon times marked as another detector, and arranging them from small to large to generate a detector 2 time series { B1，B2，B3，B4…, outputting the detector 1 time series and the detector 2 time series to the photon pair time difference lookup module;
the time probability submodule is characterized in that: under the thermo-optic simulation condition, the time probability distribution is uniform distribution;
the spatial probability submodule is characterized in that: setting the coherent time fluctuation as a Gaussian random number (Gaussian wave packet describing the random fluctuation of the width) with the mean value equal to the coherent time, wherein the time difference condition means that two adjacent photons with the time difference in the coherent time fluctuation range share the same spatial probability distribution, otherwise, different spatial probability distributions are used;
the photon pair time difference searching module uses a detector 1 time sequence { A1，A2，A3，A4… } for detector 2, the time series is { B }1，B2，B3，B4… }, will A1As a starting time tstart1Find and tstart1Adjacent next time value Bj(Bj>tstart1And min | Bj-tstart1L) as the stop time tstop1Recording time difference Δ t1＝tstop1-tstart1Completing a search cycle; then subsequently find and tstop1Adjacent next time value Ai(Ai>tstop1And min | Ai-tstop1L) as a new start time tstart2And circulating in the period to output a time difference sequence (delta t) to the time difference channel statistical module1，Δt2，Δt3，Δt4，…}；
The time difference channel counting module divides the time difference coordinate into a plurality of channels at equal intervals, counts the time difference number of each channel in the time difference sequence to obtain the time difference channel count, and finally obtains the time difference distribution curve.
A pulsar photon flow simulation method based on light field high-order correlation comprises the following steps:
s1, generating time probability distribution;
s2, collapsing the photon flow according to the time probability distribution to obtain a photon time sequence;
s3, generating the light intensity distribution of the detection surface;
s4, normalizing the light intensity distribution of the detection surface to obtain spatial probability distribution;
s5, calculating the time difference of adjacent photons by utilizing the photon time sequence, and distributing spatial probability distribution to each photon according to the time difference condition;
s6, dividing the detection surface area into a first detector aperture area, a second detector aperture area, a …, an Nth detector aperture area and an outer detector area, and performing probability summation on the spatial probability distribution of each photon according to the divided areas to obtain area probability;
s7, collapsing the photon flow according to the regional probability to obtain a photon marking sequence;
s8, selecting two detectors from the first detector to the Nth detector, recording the detectors as a detector 1 and a detector 2, and distributing photon time according to the photon time sequence and the photon mark sequence to obtain a detector 1 time sequence and a detector 2 time sequence;
s9, searching the time difference of the corresponding photon pair from the time sequence of the detector 1 and the time sequence of the detector 2, and recording the time difference into the time difference sequence;
and S10, performing channel statistics on the time difference sequence to obtain time difference channel counts, and finally obtaining a time difference distribution curve.
The invention has the technical effects that:
(1) simulating the high-order correlation characteristic of a light field according to a quantum optical probability wave collapse principle, generating photon flow with a thermo-optic bunching effect, obtaining the arrival time sequence of correlated photons of two detectors, and realizing pulsar X-ray signal detection simulation based on intensity correlation;
(2) designing an equivalent coincidence algorithm according to a coincidence circuit principle, carrying out data coincidence processing on the arrival time sequence of the associated photons to obtain a photon arrival time difference distribution curve, and realizing the intensity correlation characteristic value extraction simulation of the pulsar X-ray;
(3) intensity correlation characteristic values corresponding to detection of different spatial positions are obtained by changing the distance between detectors, and second-order interference fringes are fitted to realize intensity correlation interferometry of the pulsar angle information.
Drawings
FIG. 1 is a block diagram of a pulsar photon flow simulation system based on light field high order correlation according to the present invention.
Fig. 2 is a flowchart example of the procedure of the light field high-order correlation-based pulsar photon flow simulation method of the present invention.
Fig. 3 is a schematic diagram of the photon flow simulation subsystem of the present invention.
Fig. 4 is an exemplary graph of a 10-sample temporal probability distribution.
Fig. 5 is an exemplary graph of photon time collapse for 10 samples.
FIG. 6 is a schematic diagram of the area collapsing submodule of the present invention.
Fig. 7 is an exemplary 10-dimensional spatial probability distribution diagram, wherein (a) is the spatial probability distribution of the entire detection area, (b) is the spatial probability distribution of the region of the detector 1, and (c) is the spatial probability distribution of the region of the detector 2.
FIG. 8 is a schematic diagram of a compliant algorithm subsystem of the present invention.
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.
Fig. 1 is a module frame diagram of a light field high-order correlation-based pulsar photon flow simulation system of the present invention, which comprises a photon flow simulation subsystem 1 and a coincidence algorithm subsystem 2; the photon flow simulation subsystem 1 comprises a photon time collapse module 11, a photon space collapse module 12 and a detector time distribution module 13; the photon time collapsing module 11 comprises a time probability submodule 111 and a time collapsing submodule 112; the photon space collapse module 12 comprises a space light intensity sub-module 121, a space probability sub-module 122 and a region collapse sub-module 123; the coincidence algorithm subsystem 2 comprises a photon pair time difference searching module 21 and a time difference channel counting module 22; the time probability submodule 111 is configured to generate time probability distribution, which is used as the probability of each point of photon collapse on the time axis, and output the time probability distribution to the time collapse submodule; the time collapsing submodule 112 is configured to generate a photon time sequence, sequentially "randomly roll dice" a specified number of photons onto a time axis according to the time probability distribution, extract the time of each photon from the time axis after all the photons are collapsed, sort the time from small to large to generate a photon time sequence, and output the photon time sequence to the space probability submodule and the detector time allocation module; the spatial light intensity submodule 121 is configured to generate detection surface light intensity distribution, assign different random initial phases to characteristics of thermo-light by a classical wave optics method, obtain a plurality of different detection surface light intensity distributions, and output the detection surface light intensity distributions to the spatial probability submodule; the spatial probability submodule 122 is configured to generate spatial probability distribution, normalize the light intensity distribution of the plurality of different detection surfaces to obtain corresponding spatial probability distribution, and calculate adjacent photons according to the photon time sequenceAccording to the time difference, distributing a space probability distribution to each photon according to the time difference condition, and outputting the space probability distribution to the region collapse submodule; the area collapsing submodule 123 is configured to generate a photon marking sequence, the area where the photons are collapsed and arrive on the detection surface is divided into a first detector aperture area, a second detector aperture area, …, an nth detector aperture area, and an area outside the detector, the area probability of each area is the sum of the probabilities of the coverage areas of the corresponding areas in the spatial probability distribution, the photons are sequentially collapsed by "randomly rolling dice" according to the respective area probabilities, the photons collapsed to the first detector aperture, the second detector aperture, …, the nth detector aperture, and the outside of the detector are respectively marked as 1, 2, …, N, and 0, where N is the number of detectors on the detection surface, and the photon marking sequence is output to the detector time allocation module; the detector time allocation module 13 is configured to allocate photon time according to the photon time sequence and the photon mark sequence, select two detectors from the first detector to the nth detector, extract the photon time marked as one of the detectors by using the photon time sequence, arrange the photon time from small to large, and generate a time sequence { a of the detector 11，A2，A3，A4…, extracting the photon times marked as another detector, and arranging them from small to large to generate a detector 2 time series { B1，B2，B3，B4…, outputting the detector 1 time series and the detector 2 time series to the photon pair time difference lookup module; the photon pair time difference searching module 21 is used for searching and recording the time difference conforming to the photon pair so as to obtain the time sequence { A ] of the detector 11，A2，A3，A4… } for detector 2, the time series is { B }1，B2，B3，B4… }, in A1As a starting time tstart1Find and tstart1Adjacent next time value Bj(Bj>tstart1And min | Bj-tstart1L) as the stop time tstop1Recording time difference Δ t1＝tstop1-tstart1To complete a searchAnd (5) finding a period. Then subsequently find and tstop1Adjacent next time value Ai(Ai>tstop1And min | Ai-tstop1L) as a new start time tstart2And circulating in the period to output a time difference sequence (delta t) to the time difference channel statistical module1，Δt2，Δt3，Δt4… }; the time difference channel counting module 22 is configured to perform channel counting on the time difference sequence, divide the time difference coordinates into a plurality of channels at equal intervals, count the time difference channels according to the time difference number of each channel in the time difference sequence, and finally obtain a time difference distribution curve.
FIG. 2 is a flowchart of a procedure of the light field high-order correlation-based pulsar photon flow simulation method of the present invention, in which the main adjustable input parameters include the number of samples N and the number of photons per sample NpA coherence time tcThe detector spacing dx; generating N time probability distributions, each sampling NpCollapsing each photon to a time axis according to time probability distribution to obtain N groups of photon time sequences; in order to ensure that each photon can be at least distributed to a different space probability distribution, the original phase of the complex Gaussian random of the thermo-optic circular type is changed through a transfer function method of classical optics and a series of fixed parameters to generate the NxNpThe light intensity distribution of the amplitude detection surface is normalized to obtain NxNpAmplitude spatial probability distribution; calculating NxNpTime difference of adjacent photons of a photon, if the time difference is at coherence time tcWithin the fluctuation, two adjacent photons share the same spatial probability distribution, if the time difference is not in the coherence time tcChanging a different spatial probability distribution in the fluctuation to ensure the coherence of the photons in the coherent time fluctuation; calculating the probability of the photons appearing in the detector 1 region or the probability of the photons appearing in the detector 2 region through spatial probability distribution, wherein each photon is collapsed according to the region probability, and the photon time sequence is distributed to two detectors to obtain N groups of detector 1 time sequences and N groups of detector 2 time sequences; each sampling is carried out by utilizing a coincidence algorithm to accord two detector time sequences to obtain N groups of time difference sequences, and N groups of time difference sequences are obtainedAnd combining the time difference sequences, and performing channel statistics to obtain a time difference distribution curve.
FIG. 3 is a schematic diagram of the photon flow simulation subsystem according to the present invention, which includes a time probability sub-module 111, a time collapse sub-module 112, a space probability sub-module 122, a region collapse sub-module 123, and a detector time allocation module 13; the time probability submodule 111 generates a time probability distribution with uniform distribution under a thermo-optic simulation condition; the time collapse submodule 112 is a time sequence formed by collapsing photon groups on a time axis according to time probability distribution, and rearranges the time sequence; the spatial probability submodule 122 normalizes the light intensity distribution of the detection surface to obtain spatial probability distribution, sets the coherent time fluctuation as a gaussian random number with the mean value equal to the coherent time, and shares the same spatial probability distribution with adjacent photons of the time difference in the coherent time fluctuation, otherwise uses different spatial probability distributions; the region collapse submodule 123 calculates the probability of the photon appearing in the detector 1 region or the detector 2 region according to the spatial probability distribution, and each photon collapses according to the region probability; the detector time distribution module 13 extracts the time corresponding to the photon appearing in the detector 1 from the photon time sequence, distributes the time corresponding to the photon appearing in the detector 1 to the detector 1 time sequence, extracts the time corresponding to the photon appearing in the detector 2, and distributes the time corresponding to the photon appearing in the detector 2 to the detector 2 time sequence;
FIG. 4 is an exemplary graph of a 10-sample temporal probability distribution with 10 time-axis grids5Grid spacing is set to 10- 10s, summing the time probabilities on all the grids to 1; the temporal probability distribution is a uniform distribution;
FIG. 5 is an exemplary 10-sample photon time collapse diagram with 10 time-axis grids5Grid spacing is set to 10- 10s, the number of photons per sampling is set to be 500; each photon is randomly collapsed onto a time axis according to the time probability distribution, and when the number of photons is large, the photon distribution approaches to the time probability distribution (uniform distribution);
FIG. 6 is a schematic diagram of the area collapse sub-module of the present invention, wherein the area probability is the probability of the detector aperture coverage area summed as the probability of a photon appearing in the detector area in the spatial probability distribution, and the probability distribution of the detector aperture coverage area changes as the detector moves across the detection plane; in order to avoid the phenomenon that probability superposition possibly caused by overlapping of the positions of two detectors is larger than 1, each detector is respectively allocated with the same spatial probability distribution, the total probability sum of the two spatial probability distributions is 1, and the probability sum of one spatial probability distribution is 0.5;
FIG. 7 is an exemplary graph of 10 spatial probability distributions, in the case of one-dimensional space, the X-ray wavelength is 1nm, the light source size is 20 μm, the distance from the light source to the detection surface is 120m, the detection surface size is 0.1m, the detector aperture is 6mm, and the detector pitch is 12 mm; the spatial probability distribution is normalized detection surface light intensity distribution which can be obtained by a transfer function method of classical optics, and different random initial phases are selected to obtain different detection surface light intensity distributions; fig. 7(a) is the spatial probability distribution of the entire detection plane, fig. 7(b) is the spatial probability distribution of the region of the detector 1, with the center of the detector located at the position of 0mm of the detection plane x, fig. 7(c) is the spatial probability distribution of the region of the detector 2, with the center of the detector located at the position of 12mm of the detection plane x;
FIG. 8 is a schematic diagram of a subsystem conforming to an algorithm according to the present invention, which includes a photon pair time difference lookup module 21 and a time difference channel statistics module 22; the photon pair time difference searching module 21 sets the time sequence of the input detector 1 as { A }1，A2，A3，A4…, the input probe 2 time series is { B }1，B2，B3，B4…, generally based on probe 1 time series, will be a1As a starting time tstart1Find and tstart1Adjacent next time value Bj(Bj>tstart1And min | Bj-tstart1L) as the stop time tstop1Recording time difference Δ t1＝tstop1-tstart1Completing a search cycle; then subsequently find and tstop1Adjacent next time value Ai(Ai>tstop1And min | Ai-tstop1L) as a new start time tstart2And circulating in the period to output a time difference sequence (delta t) to the time difference channel statistical module1，Δt2，Δt3，Δt4… }; the time difference channel counting module 22 is configured to perform channel counting on the time difference sequence, divide the time difference coordinates into a plurality of channels at equal intervals, count the time difference channels according to the time difference number of each channel in the time difference sequence, and finally obtain a time difference distribution curve.
Claims (8)
1. A pulsar photon flow simulation system based on light field high order correlation is characterized by comprising:
the photon flow simulation subsystem is used for generating a photon time sequence, distributing photon time according to the principle of probability wave collapse in quantum optics and outputting two paths of detector time sequences to the coincidence algorithm subsystem;
and the coincidence algorithm subsystem is used for identifying the time difference of the coincidence photon pairs, performing algorithm coincidence on the two paths of detector time sequences and outputting a time difference distribution curve.
2. The light field higher order correlation based pulsar photon flow simulation system according to claim 1, wherein the photon flow simulation subsystem comprises:
the photon time collapse module is used for generating a photon time sequence and outputting the photon time sequence to the photon space collapse module and the detector time distribution module;
the photon space collapse module is used for generating a photon marking sequence, sequentially marking the photons collapsed to the first detector aperture, the second detector aperture, …, the Nth detector aperture and the outside of the detector as 1, 2, …, N and 0, wherein N is the number of detectors on a detection surface, and outputting the photon marking sequence to the detector time distribution module;
a detector time allocation module for allocating photons according to the photon time sequence and the photon mark sequenceOptionally selecting two detectors from the first detector to the Nth detector, extracting the photon time marked as one of the detectors by utilizing the photon time sequence, and arranging the photon time from small to large to generate a time sequence { A ] of the detector 11，A2，A3，A4…, extracting the photon times marked as another detector, and arranging them from small to large to generate a detector 2 time series { B1，B2，B3，B4…, the detector 1 time series and the detector 2 time series are output to the photon pair time difference lookup module.
3. The light field higher order correlation based pulsar photon flow simulation system according to claim 2, wherein the photon time collapsing module comprises:
the time probability submodule is used for generating time probability distribution which is used as the probability of photon collapse at each point on a time axis and outputting the time probability distribution to the time collapse submodule;
and the time collapsing submodule is used for generating a photon time sequence, sequentially collapsing the designated number of photons to a time axis by 'randomly rolling dice' according to the time probability distribution, extracting the time of each photon from the time axis after all the photons are collapsed, sequencing the time from small to large to generate the photon time sequence, and outputting the photon time sequence to the space probability submodule and the detector time distribution module.
4. The light field higher order correlation based pulsar photon flow simulation system of claim 3, wherein the temporal probability sub-module is configured to have a uniform temporal probability distribution under thermo-optic simulation conditions.
5. The light field higher order correlation based pulsar photon flow simulation system of claim 2, wherein the photon space collapsing module comprises:
the space light intensity submodule is used for generating detection surface light intensity distribution, endows different random initial phases aiming at the characteristics of thermo-light through a classical wave optics method to obtain a plurality of different detection surface light intensity distributions, and outputs the detection surface light intensity distributions to the space probability submodule;
the space probability submodule is used for generating space probability distribution, normalizing the light intensity distribution of the plurality of different detection surfaces to obtain corresponding space probability distribution, calculating the time difference of adjacent photons according to the photon time sequence, distributing one space probability distribution to each photon according to the time difference condition, and outputting the space probability distribution to the area collapse submodule;
the region collapse submodule is used for generating a photon marking sequence, regions which are reached by photon collapse are divided into a first detector aperture region, a second detector aperture region, …, an Nth detector aperture region and a detector outer region on a detection surface, the region probability of each region is the probability sum of the coverage range of the corresponding region in the space probability distribution, photons are sequentially rolled on dice randomly according to the respective region probability to collapse, the photons which are collapsed to the first detector aperture, the second detector aperture, …, the Nth detector aperture and the detector outer region are sequentially marked as 1, 2, …, N and 0, and the photon marking sequence is output to the detector time distribution module.
6. The light field higher order correlation based pulsar photon flow simulation system according to claim 5, wherein the spatial probability sub-module: setting the coherent time fluctuation as a Gaussian random number surrounding the coherent time, wherein the time difference condition means that two time-adjacent photons with the time difference in the coherent time fluctuation share the same spatial probability distribution, otherwise, different spatial probability distributions are used.
7. The light field higher order correlation based pulsar photon flow simulation system according to claim 1, wherein said algorithm-compliant subsystem comprises:
a photon pair time difference searching module for searching and recording the time difference conforming to the photon pair to detect the time sequence { A ] of the detector 11，A2，A3，A4… } reference, probe 2 time sequenceColumn is { B1，B2，B3，B4… }, in A1As a starting time tstart1Find and tstart1Adjacent next time value Bj(Bj>tstart1And min | Bj-tstart1L) as the stop time tstop1Recording time difference Δ t1＝tstop1-tstart1And completing a search cycle. Then subsequently find and tstop1Adjacent next time value Ai(Ai>tstop1And min | Ai-tstop1L) as a new start time tstart2And circulating in the period to output a time difference sequence (delta t) to the time difference channel statistical module1，Δt2，Δt3，Δt4，…}；
And the time difference channel counting module is used for carrying out channel counting on the time difference sequence, dividing the time difference coordinate into a plurality of channels at equal intervals, counting the time difference number of each channel in the time difference sequence to obtain time difference channel counting, and finally obtaining a time difference distribution curve.
8. A pulsar photon flow simulation method based on light field high-order correlation is characterized by comprising the following steps:
s1, generating time probability distribution;
s2, collapsing the photon flow according to the time probability distribution to obtain a photon time sequence;
s3, generating the light intensity distribution of the detection surface;
s4, normalizing the light intensity distribution of the detection surface to obtain spatial probability distribution;
s5, calculating the time difference of adjacent photons by utilizing the photon time sequence, and distributing spatial probability distribution to each photon according to the time difference condition;
s6, dividing the detection surface area into a first detector aperture area, a second detector aperture area, a …, an Nth detector aperture area and an outer detector area, and performing probability summation on the spatial probability distribution of each photon according to the divided areas to obtain area probability;
s7, collapsing the photon flow according to the regional probability to obtain a photon marking sequence;
s8, selecting two detectors from the first detector to the Nth detector, recording the detectors as a detector 1 and a detector 2, and distributing photon time according to the photon time sequence and the photon mark sequence to obtain a detector 1 time sequence and a detector 2 time sequence;
s9, searching the time difference of the corresponding photon pair from the time sequence of the detector 1 and the time sequence of the detector 2, and recording the time difference into the time difference sequence;
and S10, performing channel statistics on the time difference sequence to obtain time difference channel counts, and finally obtaining a time difference distribution curve.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910924850.3A CN110686702B (en) | 2019-09-27 | 2019-09-27 | Pulsar photon flow simulation system and method based on light field high-order correlation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910924850.3A CN110686702B (en) | 2019-09-27 | 2019-09-27 | Pulsar photon flow simulation system and method based on light field high-order correlation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110686702A true CN110686702A (en) | 2020-01-14 |
| CN110686702B CN110686702B (en) | 2022-10-11 |
Family
ID=69110796
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910924850.3A Active CN110686702B (en) | 2019-09-27 | 2019-09-27 | Pulsar photon flow simulation system and method based on light field high-order correlation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110686702B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111649735A (en) * | 2020-06-12 | 2020-09-11 | 中国空间技术研究院 | Pulsar signal noise reduction method based on photon probability |
| CN111666142A (en) * | 2020-06-12 | 2020-09-15 | 中国空间技术研究院 | Fast grouping photon intensity correlation method |
| CN112229397A (en) * | 2020-09-11 | 2021-01-15 | 中国科学院上海光学精密机械研究所 | Satellite angular position intensity correlation measurement system and method based on spatial modulation |
| CN112325874A (en) * | 2020-10-21 | 2021-02-05 | 中国科学院上海光学精密机械研究所 | System and method for measuring strength correlation of star angular position based on snake-shaped light path |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5512743A (en) * | 1994-01-25 | 1996-04-30 | The United States Of America As Represented By The Secretary Of The Navy | Space-based asteroid detection and monitoring system |
| US5583774A (en) * | 1994-06-16 | 1996-12-10 | Litton Systems, Inc. | Assured-integrity monitored-extrapolation navigation apparatus |
| CN103697908A (en) * | 2013-11-28 | 2014-04-02 | 中国航天科技集团公司第五研究院第五一三研究所 | X-ray pulsar signal simulation source based on random single-photon emission mechanism |
| CN103759726A (en) * | 2014-01-03 | 2014-04-30 | 西安电子科技大学 | Fast cyclostationary possion signal simulation method and hardware system thereof |
| CN104567937A (en) * | 2015-01-27 | 2015-04-29 | 中国空间技术研究院 | High-precision X-ray pulsar signal simulation method |
| CN105093236A (en) * | 2015-08-06 | 2015-11-25 | 中国科学院上海光学精密机械研究所 | Cross track direction scaling resampling method based on down-looking synthetic aperture laser imaging radar |
| CN106773615A (en) * | 2016-12-21 | 2017-05-31 | 山东航天电子技术研究所 | A kind of pulsar photon timing method based on APD detectors |
| CN108981750A (en) * | 2018-07-10 | 2018-12-11 | 西安电子科技大学 | X-ray pulse double star photon sequence emulation mode |
| CN108981749A (en) * | 2018-07-05 | 2018-12-11 | 西安电子科技大学 | X-ray pulsar angular position measurement ground simulation verifies system and its application method |
| CN109683188A (en) * | 2019-01-15 | 2019-04-26 | 南昌大学 | A kind of Qi Liankefu event and gamma events meet imaging device and method |
| CN110132286A (en) * | 2019-03-04 | 2019-08-16 | 杭州电子科技大学 | Consider the X-ray pulsar navigation method of spacecraft dynamic effect and system deviation |
| CN110186464A (en) * | 2019-05-30 | 2019-08-30 | 西安电子科技大学 | A kind of X-ray pulsar navigation TOA estimation method based on Bayesian Estimation |
-
2019
- 2019-09-27 CN CN201910924850.3A patent/CN110686702B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5512743A (en) * | 1994-01-25 | 1996-04-30 | The United States Of America As Represented By The Secretary Of The Navy | Space-based asteroid detection and monitoring system |
| US5583774A (en) * | 1994-06-16 | 1996-12-10 | Litton Systems, Inc. | Assured-integrity monitored-extrapolation navigation apparatus |
| CN103697908A (en) * | 2013-11-28 | 2014-04-02 | 中国航天科技集团公司第五研究院第五一三研究所 | X-ray pulsar signal simulation source based on random single-photon emission mechanism |
| CN103759726A (en) * | 2014-01-03 | 2014-04-30 | 西安电子科技大学 | Fast cyclostationary possion signal simulation method and hardware system thereof |
| CN104567937A (en) * | 2015-01-27 | 2015-04-29 | 中国空间技术研究院 | High-precision X-ray pulsar signal simulation method |
| CN105093236A (en) * | 2015-08-06 | 2015-11-25 | 中国科学院上海光学精密机械研究所 | Cross track direction scaling resampling method based on down-looking synthetic aperture laser imaging radar |
| CN106773615A (en) * | 2016-12-21 | 2017-05-31 | 山东航天电子技术研究所 | A kind of pulsar photon timing method based on APD detectors |
| CN108981749A (en) * | 2018-07-05 | 2018-12-11 | 西安电子科技大学 | X-ray pulsar angular position measurement ground simulation verifies system and its application method |
| CN108981750A (en) * | 2018-07-10 | 2018-12-11 | 西安电子科技大学 | X-ray pulse double star photon sequence emulation mode |
| CN109683188A (en) * | 2019-01-15 | 2019-04-26 | 南昌大学 | A kind of Qi Liankefu event and gamma events meet imaging device and method |
| CN110132286A (en) * | 2019-03-04 | 2019-08-16 | 杭州电子科技大学 | Consider the X-ray pulsar navigation method of spacecraft dynamic effect and system deviation |
| CN110186464A (en) * | 2019-05-30 | 2019-08-30 | 西安电子科技大学 | A kind of X-ray pulsar navigation TOA estimation method based on Bayesian Estimation |
Non-Patent Citations (2)
| Title |
|---|
| LING-QING ZHANG.ETC: "Time-evolution of X-ray emission nonuniformity in X-ray laser medium", 《IEEE》 * |
| 宋佳凝等: "基于相位差测量的脉冲星时间相对导航方法", 《宇航学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111649735A (en) * | 2020-06-12 | 2020-09-11 | 中国空间技术研究院 | Pulsar signal noise reduction method based on photon probability |
| CN111666142A (en) * | 2020-06-12 | 2020-09-15 | 中国空间技术研究院 | Fast grouping photon intensity correlation method |
| CN111649735B (en) * | 2020-06-12 | 2021-11-16 | 中国空间技术研究院 | Pulsar signal noise reduction method based on photon probability |
| CN112229397A (en) * | 2020-09-11 | 2021-01-15 | 中国科学院上海光学精密机械研究所 | Satellite angular position intensity correlation measurement system and method based on spatial modulation |
| CN112229397B (en) * | 2020-09-11 | 2022-08-30 | 中国科学院上海光学精密机械研究所 | Satellite angular position intensity correlation measurement system and method based on spatial modulation |
| CN112325874A (en) * | 2020-10-21 | 2021-02-05 | 中国科学院上海光学精密机械研究所 | System and method for measuring strength correlation of star angular position based on snake-shaped light path |
| CN112325874B (en) * | 2020-10-21 | 2023-02-14 | 中国科学院上海光学精密机械研究所 | System and method for measuring strength correlation of star angular position based on snake-shaped light path |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110686702B (en) | 2022-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110686702B (en) | Pulsar photon flow simulation system and method based on light field high-order correlation | |
| Schröder | Radio detection of cosmic-ray air showers and high-energy neutrinos | |
| Shi et al. | Low-frequency E-field Detection Array (LFEDA)—Construction and preliminary results | |
| Casandjian et al. | A revised catalogue of EGRET-ray sources | |
| Morales et al. | Reionization and Cosmology with 21-cm Fluctuations | |
| Asaoka et al. | Cherenkov τ shower earth-skimming method for PeV–EeV ντ observation with Ashra | |
| Guo et al. | Extended reduced-order surrogate models for scalar-tensor gravity in the strong field and applications to binary pulsars and gravitational waves | |
| Shirasaki et al. | Cross-correlation of the extragalactic gamma-ray background with luminous red galaxies | |
| Abbiendi et al. | A study of muon-electron elastic scattering in a test beam | |
| Mostafá et al. | The high-altitude water cherenkov observatory | |
| Aab et al. | The Pierre auger observatory and its upgrade | |
| Jhansi et al. | The angular resolution of GRAPES-3 EAS array after improved timing and shower front curvature correction based on age and size | |
| Romero-Shaw et al. | Inferring interference: Identifying a perturbing tertiary with eccentric gravitational wave burst timing | |
| Liu et al. | Dark-siren cosmology with Decihertz gravitational-wave detectors | |
| Epimakhov | Exploring cosmic ray origins with ground-based EAS arrays Tunka and HiSCORE | |
| Wild et al. | Triangulating the height of cosmic noise absorption: A method for estimating the characteristic energy of precipitating electrons | |
| Wang et al. | Angle measurement of pulsars based on spatially modulated X-ray intensity correlation | |
| Timashkov et al. | Ground level enhancement of December 13, 2006 observed by means of muon hodoscope | |
| Hervé | Measurement of the mass composition of the highest energy cosmic rays with the Pierre Auger observatory. | |
| Lundeen | Searches for Beyond the Standard Model Phenomena with the HAWC Detector | |
| Tuttle et al. | Horizontal electric fields from flow of auroral O+(2P) ions at sub-second resolution | |
| Whelan | A magnetic spectrometer analysis method for ultra high energy cosmic ray data. | |
| Blaess | Extracting a less model dependent cosmic ray composition from Xmax distributions | |
| He et al. | Simulation and experimental comparison of the performance of four-corner-readout plastic scintillator muon-detector system | |
| Bogomilov et al. | Neutrino factory near detector |
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 |