CN113507579B - Inter-satellite laser link capturing method and system for space gravitational wave detection - Google Patents

Inter-satellite laser link capturing method and system for space gravitational wave detection Download PDF

Info

Publication number
CN113507579B
CN113507579B CN202110891730.5A CN202110891730A CN113507579B CN 113507579 B CN113507579 B CN 113507579B CN 202110891730 A CN202110891730 A CN 202110891730A CN 113507579 B CN113507579 B CN 113507579B
Authority
CN
China
Prior art keywords
scanning
laser
capturing
ccd
angle
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
CN202110891730.5A
Other languages
Chinese (zh)
Other versions
CN113507579A (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 Jiaotong University
Original Assignee
Shanghai Jiaotong University
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 Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN202110891730.5A priority Critical patent/CN113507579B/en
Publication of CN113507579A publication Critical patent/CN113507579A/en
Application granted granted Critical
Publication of CN113507579B publication Critical patent/CN113507579B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/70SSIS architectures; Circuits associated therewith
    • H04N25/71Charge-coupled device [CCD] sensors; Charge-transfer registers specially adapted for CCD sensors
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/11Arrangements specific to free-space transmission, i.e. transmission through air or vacuum
    • H04B10/118Arrangements specific to free-space transmission, i.e. transmission through air or vacuum specially adapted for satellite communication

Abstract

The invention provides an inter-satellite laser link capturing method and system for space gravitational wave detection, comprising the following steps: step 1: scanning a target terminal area through laser, adjusting the visual axes of a transmitting end and a detecting end according to response signals of a detection target, and establishing a communication link after the visual axes are aligned; step 2: setting laser capturing parameters including a capturing area, a laser beam divergence angle and a scanning mode; step 3: constructing a Simulink model to capture laser, wherein the Simulink model comprises CCD output simulation and CCD angle calculation; step 4: setting an initial alignment error, calculating a CCD capturing error, and adjusting capturing precision. According to the invention, through analyzing the working principle of laser capturing of space gravitational wave detection, a laser capturing mode and corresponding capturing parameters are designed, and a Simulink simulation model of laser capturing is built on the basis, so that the simulation model is packaged, and the application of practical engineering is facilitated.

Description

Inter-satellite laser link capturing method and system for space gravitational wave detection
Technical Field
The invention relates to the technical field of laser communication, in particular to an inter-satellite laser link capturing method and system for space gravitational wave detection.
Background
In the space gravitational wave detection, in order to realize scientific measurement of gravitational wave signals, a million kilometer laser link between satellites is required to be constructed, laser capturing is the first stage of the link establishment process, the capturing process is to continuously adjust the visual axis directions of two laser communication terminals, realize mutual alignment, and finally realize closed loop of an optical path between the two terminals.
In the laser capturing process, a Charge Coupled Device (CCD) is used as a capturing detector, and because the CCD field of view is smaller, a star Sensor (STR) is adopted to perform initial alignment, and CCD laser capturing is performed after the initial alignment. Because of errors of the sensor and the actuating mechanism, an initial aiming error exists between the two satellites after initial alignment, in space gravitational wave detection, the distance between the two satellites is in the order of million kilometers, the actual position of the satellite is greatly deviated from the reference position due to the initial aiming error, and energy is lost in the transmission process of laser, so that a proper capturing scheme and capturing parameters are required to be determined for realizing rapid and accurate laser capturing, and the capturing uncertainty area, the scanning step length, the scanning mode and the selection of laser beam divergence angle are included. The invention designs a space gravitational wave detection inter-satellite laser link capturing method based on a Matlab/Simulink platform, and establishes a Simulink model through analysis of a capturing principle to simulate an inter-satellite link capturing process in space gravitational wave detection.
Patent document CN101567721a (application number: CN 200910071922.0) discloses a beam capturing scanning method for quickly establishing a laser link between a relay star and a user star, which comprises the following steps: step one: the relay star terminal sends out a beacon beam to a scanning point of the area to be scanned; step two: the relay star terminal sends out a beacon beam to the next scanning point in the area to be scanned at time intervals; step three: when the beacon light beam in the first step or the second step reaches the area to be scanned, the satellite terminal of the user sends out a return light beam; step four: the relay star terminal judges whether the return light beam in the step three is received or not, if the judgment result is negative, the step two is executed; if the result is yes, the relay star terminal stops sending out the beacon beam, and at the moment, the user star is positioned at the previous scanning point, and the relay star captures the user star successfully. However, this patent suffers from energy loss during laser transmission during scanning, thus causing positional deviation.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide an inter-satellite laser link capturing method and system for space gravitational wave detection.
The inter-satellite laser link capturing method for space gravitational wave detection provided by the invention comprises the following steps:
step 1: scanning a target terminal area through laser, adjusting the visual axes of a transmitting end and a detecting end according to response signals of a detection target, and establishing a communication link after the visual axes are aligned;
step 2: setting laser capturing parameters including a capturing area, a laser beam divergence angle and a scanning mode;
step 3: constructing a Simulink model to capture laser, wherein the Simulink model comprises CCD output simulation and CCD angle calculation;
step 4: setting an initial alignment error, calculating a CCD capturing error, and adjusting capturing precision.
Preferably, after the initial alignment of the execution mechanisms of the transmitting end and the detecting end, the receiving view field of the detecting end is larger than the capturing area, and the beam divergence angle of the laser beam of the transmitting end is smaller than the capturing area.
Preferably, the scanning mode is a spiral scanning mode, the equal-pitch equal-linear speed spiral scanning is adopted according to the probability distribution and the shape characteristics of the capturing area, and the area where the target terminal exists is scanned from the center position with the highest probability of occurrence of the target;
the polar equation for helical scanning is:
Figure BDA0003196164710000021
Figure BDA0003196164710000022
wherein θ is the polar angle, ρ is the polar diameter, I θ For scanning step length, θ b Is the beam divergence angle of the light beam;
the scan trajectory curves for azimuth and elevation directions are:
Figure BDA0003196164710000023
Figure BDA0003196164710000024
wherein x is the azimuth direction and y is the pitch direction;
in the polar coordinate system, scanning is started from a central point, and the polar angle θ (i) and the polar diameter ρ (i) of each corresponding scanning point are expressed as:
Figure BDA0003196164710000025
Figure BDA0003196164710000026
when the scanning stay points are converted into a rectangular coordinate system, the coordinates of each scanning stay point are as follows:
Figure BDA0003196164710000031
Figure BDA0003196164710000032
where i denotes the ith scan point in the scan path, i is 1,2,3, …, N scan
Figure BDA0003196164710000033
Wherein N is scan For maximum number of scanning positions, θ FOU Representing the half-angle size of the capture area.
Preferably, the CCD output simulation includes:
calculating the offset on the CCD according to the incidence angle of the laser, and simulating the position of the light spot on the CCD according to the offset;
simulating the light spot intensity of each pixel on the image by adopting Gaussian distribution, and outputting the light spot intensity as a simulated light spot image;
the components of the incident laser light in azimuth and pitch directions are radians of each pixel in x, y directions.
Preferably, the CCD angle calculation includes:
the light spot centroid is calculated according to the light intensity of each pixel collected on the photosensitive surface of the CCD camera, the beacon light spot image on the detector consists of m multiplied by n pixels, and the position coordinate of each pixel is (x i ,y j ) The corresponding gray value is p ij (x, y), m and n are the number of transverse and longitudinal pixels of the beacon spot image, and i and j are serial numbers in the x and y directions respectively;
the calculation formula of the centroid (X, Y) of the light spot is as follows:
Figure BDA0003196164710000034
Figure BDA0003196164710000035
the invention provides an inter-satellite laser link capturing system for space gravitational wave detection, which comprises the following components:
module M1: scanning a target terminal area through laser, adjusting the visual axes of a transmitting end and a detecting end according to response signals of a detection target, and establishing a communication link after the visual axes are aligned;
module M2: setting laser capturing parameters including a capturing area, a laser beam divergence angle and a scanning mode;
module M3: constructing a Simulink model to capture laser, wherein the Simulink model comprises CCD output simulation and CCD angle calculation;
module M4: setting an initial alignment error, calculating a CCD capturing error, and adjusting capturing precision.
Preferably, after the initial alignment of the execution mechanisms of the transmitting end and the detecting end, the receiving view field of the detecting end is larger than the capturing area, and the beam divergence angle of the laser beam of the transmitting end is smaller than the capturing area.
Preferably, the scanning mode is a spiral scanning mode, the equal-pitch equal-linear speed spiral scanning is adopted according to the probability distribution and the shape characteristics of the capturing area, and the area where the target terminal exists is scanned from the center position with the highest probability of occurrence of the target;
the polar equation for helical scanning is:
Figure BDA0003196164710000041
Figure BDA0003196164710000042
wherein θ is the polar angle, ρ is the polar diameter, I θ For scanning step length, θ b Is a light beamIs a beam divergence angle of (2);
the scan trajectory curves for azimuth and elevation directions are:
Figure BDA0003196164710000043
Figure BDA0003196164710000044
wherein x is the azimuth direction and y is the pitch direction;
in the polar coordinate system, scanning is started from a central point, and the polar angle θ (i) and the polar diameter ρ (i) of each corresponding scanning point are expressed as:
Figure BDA0003196164710000045
Figure BDA0003196164710000046
when the scanning stay points are converted into a rectangular coordinate system, the coordinates of each scanning stay point are as follows:
Figure BDA0003196164710000047
/>
Figure BDA0003196164710000048
where i denotes the ith scan point in the scan path, i is 1,2,3, …, N scan
Figure BDA0003196164710000049
Wherein N is scan For maximum number of scanning positions, θ FOU Representing the half-angle size of the capture area.
Preferably, the CCD output simulation includes:
calculating the offset on the CCD according to the incidence angle of the laser, and simulating the position of the light spot on the CCD according to the offset;
simulating the light spot intensity of each pixel on the image by adopting Gaussian distribution, and outputting the light spot intensity as a simulated light spot image;
the components of the incident laser light in azimuth and pitch directions are radians of each pixel in x, y directions.
Preferably, the CCD angle calculation includes:
the light spot centroid is calculated according to the light intensity of each pixel collected on the photosensitive surface of the CCD camera, the beacon light spot image on the detector consists of m multiplied by n pixels, and the position coordinate of each pixel is (x i ,y j ) The corresponding gray value is p ij (x, y), m and n are the number of transverse and longitudinal pixels of the beacon spot image, and i and j are serial numbers in the x and y directions respectively;
the calculation formula of the centroid (X, Y) of the light spot is as follows:
Figure BDA0003196164710000051
Figure BDA0003196164710000052
compared with the prior art, the invention has the following beneficial effects:
according to the invention, through analyzing the working principle of laser capturing of space gravitational wave detection, a laser capturing mode and corresponding capturing parameters are designed, and a Simulink simulation model of laser capturing is built on the basis, so that the simulation model is packaged, and the application of practical engineering is facilitated.
Drawings
Other features, objects and advantages of the present invention will become more apparent upon reading of the detailed description of non-limiting embodiments, given with reference to the accompanying drawings in which:
FIG. 1 is a flow chart of a spatial laser acquisition method;
FIG. 2 is a schematic diagram of a space laser captured Simulink simulation model;
fig. 3 is a schematic view of a helical scan.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the present invention, but are not intended to limit the invention in any way. It should be noted that variations and modifications could be made by those skilled in the art without departing from the inventive concept. These are all within the scope of the present invention.
Examples:
referring to fig. 1, the process of the spatial laser capturing method according to the present invention includes the following steps:
(1) Determining a capture scheme and a capture mode
The capture scheme employed is a "laser scanning capture scheme". The scheme adopts laser to scan the possible area of the target terminal, and detects the response signal of the target, continuously adjusts the visual axes of the two terminals, and gradually and accurately aligns and establishes a communication link.
The capturing mode is a staring/scanning capturing mode, after the execution mechanisms of the transmitting end and the detecting end are initially aligned, the receiving view field of the detecting end is larger than the capturing uncertainty area, and the beam divergence angle of the laser beam of the transmitting end is smaller than the capturing uncertainty area.
(2) Determining acquisition parameters
The size of the capture uncertainty area is determined by a plurality of factors, mainly by the total pointing error of the visual axis of the transmitting terminal in space and the capture probability, and the uncertainty area theta is selected in simulation FOU =15.2μrad;
The beam divergence angle of the laser is given by the ideal full width of the emission beam reduced by the gesture shake of the local spacecraft and the pointing shake of the telescope, and the beam divergence angle is selected as theta in the simulation b =1.43μrad;
The scanning mode is a spiral scanning mode, the capturing efficiency is high, the method is easy to realize, and according to the characteristics of the probability distribution and the uncertain area shape of the target, the area where the target terminal possibly exists is scanned by adopting equal-pitch equal-linear speed spiral scanning from the center position with the highest probability of occurrence of the target; as shown in fig. 2, a spiral scanning schematic diagram is shown, wherein an outer circle represents an uncertain region, a spiral circle represents a laser scanning path, and the size of the inner circle represents a beam divergence angle of laser;
the polar equation for helical scanning is:
Figure BDA0003196164710000061
wherein θ is the polar angle, ρ is the polar diameter, I θ For scanning step length, also pitch of helix, angle θ with beam divergence of beam b In connection, to ensure complete coverage of the target uncertainty region, it is common to take:
Figure BDA0003196164710000062
the scan trajectory curves for azimuth and elevation directions are then:
Figure BDA0003196164710000063
x represents the azimuth direction and y represents the pitch direction;
the maximum number of scan positions N required is due to the fact that sufficient overlap must be provided between the scan positions scan It can be estimated that:
Figure BDA0003196164710000064
wherein θ FOU Representing the half angle size of the uncertainty region;
in the polar coordinate system, the polar angle θ (i) and the polar diameter ρ (i) of each corresponding scanning point, starting from the center point, can be expressed as:
Figure BDA0003196164710000071
when the scanning stay points are converted into a rectangular coordinate system, the coordinates of each scanning stay point are as follows:
Figure BDA0003196164710000072
where i denotes the ith scan point in the scan path, i is 1,2,3, …, N scan
(3) Building Simulink model
As shown in fig. 3, the Simulink model mainly includes 3 parts:
and a scanning module: the input is the initial aiming error, and the output is the position of the receiving end in the uncertain region when receiving the signal of the transmitting end, and the position is expressed by the angle offset.
CCD analog output module: calculating the current laser incidence angle according to the output of the scanning module, calculating the offset on the CCD to simulate the position of a light spot on the CCD, simulating the light intensity of the light spot of each pixel on the image by adopting Gaussian distribution, wherein p in the output represents the light intensity of the pixel point on the CCD, dx represents the radian represented by each pixel in the x direction, dy represents the radian represented by each pixel in the y direction, and yaw and pitch represent the angles of the incident laser in the yaw and pitch directions at the current moment respectively.
CCD angle calculation module: calculating to obtain a spot centroid according to the intensity of each pixel acquired on a photosensitive surface of a CCD camera by adopting a centroid algorithm, and assuming that a beacon spot image on a detector consists of m multiplied by n pixels, the position coordinate of each pixel is (x i ,y j ) The corresponding gray value is p ij (x, y). The spot centroid (X, Y) calculation formula can be expressed as:
Figure BDA0003196164710000073
Figure BDA0003196164710000074
the outputs yaw_c and pitch_c in the block represent the offset of the current incident laser calculated by the CCD in yaw and pitch directions.
(4) Setting initial error, calculating CCD capture error
After the Simulink model is built, setting an initial alignment error, wherein the initial alignment error obeys Gaussian distribution, selecting a random module in Simulink as an initial error module, and selecting a variance of 5.184 multiplied by 10 in the embodiment -11 The obtained CCD capturing error is 3.206 ×10 in the pitching direction -7 rad, heading direction 6.86×10 -7 rad;
The error can ensure that the incident laser can enter the sensor visual field range with higher precision after the actuating mechanism adjusts the incident angle of the laser, namely the capturing is considered to be successful.
The invention provides an inter-satellite laser link capturing system for space gravitational wave detection, which comprises the following components: module M1: scanning a target terminal area through laser, adjusting the visual axes of a transmitting end and a detecting end according to response signals of a detection target, and establishing a communication link after the visual axes are aligned; module M2: setting laser capturing parameters including a capturing area, a laser beam divergence angle and a scanning mode; module M3: constructing a Simulink model to capture laser, wherein the Simulink model comprises CCD output simulation and CCD angle calculation; module M4: setting an initial alignment error, calculating a CCD capturing error, and adjusting capturing precision.
After the execution mechanisms of the transmitting end and the detecting end are initially aligned, the receiving view field of the detecting end is larger than the capturing area, and the beam divergence angle of the laser beam of the transmitting end is smaller than the capturing area. The scanning mode is a spiral scanning mode, the equal-pitch equal-linear speed spiral scanning is adopted according to the probability distribution and the shape characteristics of the capturing area, and the area where the target terminal exists is scanned from the center position with the highest probability of occurrence of the target; the polar equation for helical scanning is:
Figure BDA0003196164710000081
Figure BDA0003196164710000082
wherein θ is the polar angle, ρ is the polar diameter, I θ For scanning step length, θ b Is the beam divergence angle of the light beam;
the scan trajectory curves for azimuth and elevation directions are:
Figure BDA0003196164710000083
Figure BDA0003196164710000084
wherein x is the azimuth direction and y is the pitch direction;
in the polar coordinate system, scanning is started from a central point, and the polar angle θ (i) and the polar diameter ρ (i) of each corresponding scanning point are expressed as:
Figure BDA0003196164710000085
Figure BDA0003196164710000086
when the scanning stay points are converted into a rectangular coordinate system, the coordinates of each scanning stay point are as follows:
Figure BDA0003196164710000091
Figure BDA0003196164710000092
where i denotes the ith scan point in the scan path, i is 1,2,3, …, N scan
Figure BDA0003196164710000093
Wherein N is scan For maximum number of scanning positions, θ FOU Representing the half-angle size of the capture area.
The CCD output simulation includes: calculating the offset on the CCD according to the incidence angle of the laser, and simulating the position of the light spot on the CCD according to the offset; simulating the light spot intensity of each pixel on the image by adopting Gaussian distribution, and outputting the light spot intensity as a simulated light spot image; the components of the incident laser light in azimuth and pitch directions are radians of each pixel in x, y directions.
The CCD angle calculation includes: the centroid algorithm is adopted, the centroid of the light spot is calculated according to the light intensity of each pixel collected on the photosensitive surface of the CCD camera, the beacon light spot image on the detector consists of m multiplied by n pixels, and the position coordinate of each pixel is (x i ,y j ) The corresponding gray value is p ij (x, y), m and n are the number of transverse and longitudinal pixels of the beacon spot image, and i and j are serial numbers in the x and y directions respectively; the calculation formula of the centroid (X, Y) of the light spot is as follows:
Figure BDA0003196164710000094
Figure BDA0003196164710000095
those skilled in the art will appreciate that the systems, apparatus, and their respective modules provided herein may be implemented entirely by logic programming of method steps such that the systems, apparatus, and their respective modules are implemented as logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers, etc., in addition to the systems, apparatus, and their respective modules being implemented as pure computer readable program code. Therefore, the system, the apparatus, and the respective modules thereof provided by the present invention may be regarded as one hardware component, and the modules included therein for implementing various programs may also be regarded as structures within the hardware component; modules for implementing various functions may also be regarded as being either software programs for implementing the methods or structures within hardware components.
The foregoing describes specific embodiments of the present invention. It is to be understood that the invention is not limited to the particular embodiments described above, and that various changes or modifications may be made by those skilled in the art within the scope of the appended claims without affecting the spirit of the invention. The embodiments of the present application and features in the embodiments may be combined with each other arbitrarily without conflict.

Claims (2)

1. An inter-satellite laser link acquisition method for spatial gravitational wave detection, comprising:
step 1: scanning a target terminal area through laser, adjusting the visual axes of a transmitting end and a detecting end according to response signals of a detection target, and establishing a communication link after the visual axes are aligned;
step 2: setting laser capturing parameters including a capturing area, a laser beam divergence angle and a scanning mode;
step 3: constructing a Simulink model for laser capturing, wherein the Simulink model comprises a scanning module, CCD output simulation and CCD angle calculation;
step 4: setting an initial alignment error, calculating a CCD capturing error, and adjusting capturing precision;
after the execution mechanisms of the transmitting end and the detecting end are initially aligned, the receiving view field of the detecting end is larger than the capturing area, and the beam divergence angle of the laser beam of the transmitting end is smaller than the capturing area;
the scanning mode is a spiral scanning mode, the equal-pitch equal-linear speed spiral scanning is adopted according to the probability distribution and the shape characteristics of the capturing area, and the area where the target terminal exists is scanned from the center position with the highest probability of occurrence of the target;
the polar equation for helical scanning is:
Figure QLYQS_1
Figure QLYQS_2
wherein (1)>
Figure QLYQS_3
For polar angle>
Figure QLYQS_4
Is of the polar diameter>
Figure QLYQS_5
For the scanning step size +.>
Figure QLYQS_6
Is the beam divergence angle of the light beam;
the scan trajectory curves for azimuth and elevation directions are:
Figure QLYQS_7
Figure QLYQS_8
wherein x is the azimuth direction and y is the pitch direction;
in the polar coordinate system, scanning is started from a central point, and the polar angle of each corresponding scanning point
Figure QLYQS_9
And polar diameter->
Figure QLYQS_10
Expressed as:
Figure QLYQS_11
Figure QLYQS_12
conversion to straightWhen in the angular coordinate system, the coordinates of each scanning resident point are as follows:
Figure QLYQS_13
Figure QLYQS_14
wherein i represents the ith scanning point in the scanning path, i is 1,2,3, …,>
Figure QLYQS_15
;
Figure QLYQS_16
wherein (1)>
Figure QLYQS_17
For the maximum number of scanning positions->
Figure QLYQS_18
Representing the half-angle size of the capture area;
the scanning module includes: the input is initial alignment error, the output is the position of the receiving end in an uncertain region when receiving a signal of the transmitting end, and the position is expressed by angle offset;
the CCD output simulation includes:
calculating a current laser incidence angle according to the output of the scanning module, calculating an offset on the CCD according to the laser incidence angle, and simulating the position of a light spot on the CCD according to the offset;
simulating the light spot intensity of each pixel on the image by adopting Gaussian distribution, and outputting the light spot intensity as a simulated light spot image;
the components of the incident laser in azimuth and pitching directions are radians of each pixel in the x-axis direction and the y-axis direction;
the CCD angle calculation includes:
the centroid of the light spot is calculated according to the light intensity of each pixel collected on the photosensitive surface of the CCD camera,
the beacon spot image on the detector consists of m×n pixels, and the position coordinate of each pixel is
Figure QLYQS_19
The corresponding gray value is +.>
Figure QLYQS_20
M and n are the number of horizontal and vertical pixels of the beacon spot image, and i and j are serial numbers in the x and y directions respectively;
the calculation formula of the centroid (X, Y) of the light spot is as follows:
Figure QLYQS_21
Figure QLYQS_22
2. an inter-satellite laser link acquisition system for spatial gravitational wave detection, comprising:
module M1: scanning a target terminal area through laser, adjusting the visual axes of a transmitting end and a detecting end according to response signals of a detection target, and establishing a communication link after the visual axes are aligned;
module M2: setting laser capturing parameters including a capturing area, a laser beam divergence angle and a scanning mode;
module M3: constructing a Simulink model for laser capturing, wherein the Simulink model comprises a scanning module, CCD output simulation and CCD angle calculation;
module M4: setting an initial alignment error, calculating a CCD capturing error, and adjusting capturing precision;
after the execution mechanisms of the transmitting end and the detecting end are initially aligned, the receiving view field of the detecting end is larger than the capturing area, and the beam divergence angle of the laser beam of the transmitting end is smaller than the capturing area;
the scanning mode is a spiral scanning mode, the equal-pitch equal-linear speed spiral scanning is adopted according to the probability distribution and the shape characteristics of the capturing area, and the area where the target terminal exists is scanned from the center position with the highest probability of occurrence of the target;
the polar equation for helical scanning is:
Figure QLYQS_23
Figure QLYQS_24
wherein (1)>
Figure QLYQS_25
For polar angle>
Figure QLYQS_26
Is of the polar diameter>
Figure QLYQS_27
For the scanning step size +.>
Figure QLYQS_28
Is the beam divergence angle of the light beam;
the scan trajectory curves for azimuth and elevation directions are:
Figure QLYQS_29
Figure QLYQS_30
wherein x is the azimuth direction and y is the pitch direction;
in the polar coordinate system, scanning is started from a central point, and the polar angle of each corresponding scanning point
Figure QLYQS_31
And polar diameter->
Figure QLYQS_32
Expressed as:
Figure QLYQS_33
Figure QLYQS_34
when the scanning stay points are converted into a rectangular coordinate system, the coordinates of each scanning stay point are as follows: />
Figure QLYQS_35
Figure QLYQS_36
Wherein i represents the ith scanning point in the scanning path, i is 1,2,3, …,>
Figure QLYQS_37
Figure QLYQS_38
wherein (1)>
Figure QLYQS_39
For the maximum number of scanning positions->
Figure QLYQS_40
Representing the half-angle size of the capture area;
the scanning module includes: the input is initial alignment error, the output is the position of the receiving end in an uncertain region when receiving a signal of the transmitting end, and the position is expressed by angle offset;
the CCD output simulation includes:
calculating a current laser incidence angle according to the output of the scanning module, calculating an offset on the CCD according to the laser incidence angle, and simulating the position of a light spot on the CCD according to the offset;
simulating the light spot intensity of each pixel on the image by adopting Gaussian distribution, and outputting the light spot intensity as a simulated light spot image;
the components of the incident laser in azimuth and pitching directions are radians of each pixel in the x-axis direction and the y-axis direction;
the CCD angle calculation includes:
the centroid of the light spot is calculated according to the light intensity of each pixel collected on the photosensitive surface of the CCD camera,
the beacon spot image on the detector consists of m×n pixels, and the position coordinate of each pixel is
Figure QLYQS_41
The corresponding gray value is +.>
Figure QLYQS_42
M and n are the number of horizontal and vertical pixels of the beacon spot image, and i and j are serial numbers in the x and y directions respectively;
the calculation formula of the centroid (X, Y) of the light spot is as follows:
Figure QLYQS_43
Figure QLYQS_44
/>
CN202110891730.5A 2021-08-04 2021-08-04 Inter-satellite laser link capturing method and system for space gravitational wave detection Active CN113507579B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110891730.5A CN113507579B (en) 2021-08-04 2021-08-04 Inter-satellite laser link capturing method and system for space gravitational wave detection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110891730.5A CN113507579B (en) 2021-08-04 2021-08-04 Inter-satellite laser link capturing method and system for space gravitational wave detection

Publications (2)

Publication Number Publication Date
CN113507579A CN113507579A (en) 2021-10-15
CN113507579B true CN113507579B (en) 2023-05-12

Family

ID=78015005

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110891730.5A Active CN113507579B (en) 2021-08-04 2021-08-04 Inter-satellite laser link capturing method and system for space gravitational wave detection

Country Status (1)

Country Link
CN (1) CN113507579B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115065410B (en) * 2022-07-06 2023-09-15 北京理工大学 Method for scanning and capturing high-precision laser links between satellites
CN115941035B (en) * 2022-12-06 2023-06-02 威海激光通信先进技术研究院 Quick calibration method for capturing view field of laser communication terminal based on two-dimensional scanning light beam
CN117590484B (en) * 2024-01-12 2024-03-19 中国民用航空总局第二研究所 Airport runway foreign matter positioning method and system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101567721A (en) * 2009-04-29 2009-10-28 哈尔滨工业大学 Beam-capturing scanning method for rapidly establishing laser link between relay satellite and user satellite
CN111970058A (en) * 2020-09-22 2020-11-20 长春理工大学 Beacon-free rapid acquisition method for laser communication between satellites

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101207440B (en) * 2007-12-19 2011-08-31 哈尔滨工业大学 Capture performance test method of space optical communication system
CN102075243A (en) * 2010-12-28 2011-05-25 哈尔滨工业大学 Error detection device and control method for laser communication link beams
CN112953631A (en) * 2018-12-29 2021-06-11 长沙天仪空间科技研究院有限公司 Laser capturing and communication system and method based on satellite
CN112436892A (en) * 2020-12-30 2021-03-02 中国电子科技集团公司第三十四研究所 Wireless optical communication automatic scanning, capturing and tracking system and method thereof

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101567721A (en) * 2009-04-29 2009-10-28 哈尔滨工业大学 Beam-capturing scanning method for rapidly establishing laser link between relay satellite and user satellite
CN111970058A (en) * 2020-09-22 2020-11-20 长春理工大学 Beacon-free rapid acquisition method for laser communication between satellites

Also Published As

Publication number Publication date
CN113507579A (en) 2021-10-15

Similar Documents

Publication Publication Date Title
CN113507579B (en) Inter-satellite laser link capturing method and system for space gravitational wave detection
CN113296128B (en) System and method for establishing high-capture-rate low-orbit inter-satellite laser communication link
CN101738161B (en) Equipment and method for measuring six-dimensional pose of moving object
CN102324962B (en) Method for capturing, tracking and processing satellite optical communication
CN111950873B (en) Satellite real-time guiding task planning method and system based on deep reinforcement learning
US20070247612A1 (en) System and method for improving lidar data fidelity using pixel-aligned lidar/electro-optic data
US20220237738A1 (en) Information processing device, information processing method, information processing program, image processing device, and image processing system for associating position information with captured images
CN105827310B (en) A kind of optical antenna for multipoint laser communication based on wide-angle beam expanding lens
CN1825787A (en) Method for establishing laser link between low orbit safellite and ground station
CN109543292B (en) Simulation method of inter-satellite link signal capturing process based on orbit prediction
CN113486447B (en) Inter-satellite solar time prediction method and device and electronic equipment
CN105043417A (en) Multi-target continuous imaging drift angle compensation method
CN111669224A (en) Inter-satellite laser communication aiming deviation on-orbit measurement and correction method
CN115096293A (en) Multi-robot cooperation detection map construction method and device and mine sweeping robot
CN110967005B (en) Imaging method and imaging system for on-orbit geometric calibration through star observation
CN112874819B (en) Attitude planning method for mutually establishing laser links after spacecraft orbit entering
CN112000124A (en) Unmanned aerial vehicle inspection method applied to power grid
Burisch et al. Heliostat calibration using attached cameras and artificial targets
US20190302806A1 (en) Image capturing system for shape measurement of structure, on-board controller
CN110233666A (en) A kind of optical antenna for multipoint laser communication based on biconvex reflection microscope group
CN104393932B (en) Real-time correction method for optical axis of telescope of quantum communication ground station
CN201285280Y (en) Equipment for measuring six-dimension pose of moving object
CN106405855A (en) Vortex light beam real-time alignment system and method having OAM (Orbital Angular Momentum)
CN116338716B (en) Multi-target association method of air-ground unmanned system based on azimuth topological structure
CN115065410B (en) Method for scanning and capturing high-precision laser links between satellites

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