[go: up one dir, main page]

CN102662188A - Initializing method of mobile satellite communication antenna - Google Patents

Initializing method of mobile satellite communication antenna Download PDF

Info

Publication number
CN102662188A
CN102662188A CN2012101035437A CN201210103543A CN102662188A CN 102662188 A CN102662188 A CN 102662188A CN 2012101035437 A CN2012101035437 A CN 2012101035437A CN 201210103543 A CN201210103543 A CN 201210103543A CN 102662188 A CN102662188 A CN 102662188A
Authority
CN
China
Prior art keywords
antenna
mtd
satellite
angle
msubsup
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.)
Pending
Application number
CN2012101035437A
Other languages
Chinese (zh)
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.)
BEIJING STARNETO TECHNOLOGY Co Ltd
Original Assignee
BEIJING STARNETO TECHNOLOGY Co Ltd
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 BEIJING STARNETO TECHNOLOGY Co Ltd filed Critical BEIJING STARNETO TECHNOLOGY Co Ltd
Priority to CN2012101035437A priority Critical patent/CN102662188A/en
Publication of CN102662188A publication Critical patent/CN102662188A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/14Relay systems
    • H04B7/15Active relay systems
    • H04B7/185Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851Systems using a satellite or space-based relay
    • H04B7/18517Transmission equipment in earth stations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Radio Relay Systems (AREA)

Abstract

The invention relates to an initializing method of a mobile satellite communication antenna. According to the method, the mobile satellite communication antenna provided with a low-cost inertia navigation system can fast align a satellite and accurately estimate an initial course angle of the low-cost inertia navigation system. The method comprises the steps of: firstly calculating a satellite searching instruction by utilizing the initial course and gesture information of the low-cost inertia navigation system, then controlling a pitch axis and a polarization axis of the antenna to rotate according to calculated pitch instruction angle and polarization instruction angle; controlling an azimuth axis to scan between 0 and 360 degrees at a constant speed, and finding out an accurate course instruction angle through a satellite beacon recognition technology; accurately estimating the initial course angle error of the low-cost inertia navigation system by utilizing three-coordinate calculation and an iterative algorithm; finally calculating the satellite searching instruction of the antenna again according to the estimated inertia navigation course and gesture information, and controlling the antenna to rotate according to the instruction angle and accurately align the satellite. After finishing initialization of the antenna by using the method, the antenna can track the satellite in real time with high precision.

Description

Initialization method of mobile satellite communication antenna
Technical Field
The application relates to an initialization method of a mobile satellite communication antenna, and relates to the fields of automatic control, inertial navigation and mobile satellite communication. The method can be widely applied to various mobile satellite communication devices, such as vehicle-mounted, ship-mounted and airborne mobile satellite communication antennas.
Background
The mobile satellite communication equipment can enable the satellite antenna to be always aligned to the geostationary satellite on a mobile carrier, realizes high-bandwidth real-time image and data transmission, and has very wide application prospect in departments of national defense, frontier defense, counter terrorism, emergency disaster relief, government and the like.
The mobile satellite communication antenna consists of an antenna feed system and a servo control system. The antenna feed system comprises an antenna surface, a duplexer, a filter, a down converter, an up converter, a waveguide and other components, and forms basic elements of the satellite communication system; the servo control system needs to drive the antenna surface to complete the functions of scanning, finding the satellite, stably tracking and the like, and the antenna surface is ensured to accurately point to the satellite in the moving process of the carrier.
The initialization process of the mobile satellite communication antenna is a process of completing initial satellite finding and determining the initial course of a carrier after the antenna is powered on, is a premise of ensuring that the antenna can track a satellite in real time and high precision under a dynamic state, and is also one of main core technologies of the mobile satellite communication antenna. Currently, the process of antenna initialization mainly refers to calculating the satellite-finding command angle (pitch angle, polarization angle, azimuth angle) of the antenna by means of position and attitude angle information given by a high-precision vehicle-mounted inertial navigation system, and controlling the antenna to align to a satellite. The scheme has high requirement on the inertial navigation system, in order to ensure that the antenna can lock the maximum satellite signal, the autonomous north-seeking precision of the inertial navigation system must reach more than 0.1 degrees, and the cost of the high-precision inertial navigation system far exceeds the cost of the mobile satellite communication antenna, which cannot be borne by most development units and users. In order to effectively reduce the cost, a low-precision vehicle-mounted inertial navigation system is required to be selected, and the low-precision inertial navigation system cannot automatically find north and output a correct course, so that an antenna cannot obtain an accurate satellite-aiming instruction angle. In order to solve the problem of initialization of a mobile satellite communication antenna based on low-precision inertial navigation, most of the current research units adopt the technology of 0-360-degree scanning of an antenna azimuth axis and satellite beacon signal peak recognition to realize antenna satellite finding, but the problem of calculating the initial course of an inertial navigation system by using satellite finding auxiliary information cannot be solved well, and particularly under the condition that a carrier provided with an antenna has a large pitch angle or a large roll angle, the course error of inertial navigation estimated by the existing method is large, so that the tracking performance of the mobile satellite communication antenna is greatly reduced.
Disclosure of Invention
Aiming at the technical problem, the method for initializing the mobile satellite communication antenna utilizes the antenna to scan by 0-360 degrees to search the satellite, and then estimates the initial course angle of the low-precision vehicle-mounted inertial navigation system through an algorithm of complete coordinate conversion and iterative progression. The method has good environmental adaptability, is suitable for the initial satellite finding process of the carrier under any attitude angle, and after the initialization is finished, the initial course precision of the inertial navigation system calculated by satellite finding auxiliary information can reach more than 0.1 degree. Compared with the existing initialization algorithm based on plane coordinate calculation, the method provided by the application can accurately estimate the initial course of the vehicle-mounted inertial navigation system in any attitude environment, so that the antenna can be ensured to accurately track the satellite after initialization is finished.
The antenna initialization method provided by the application has the design idea that: firstly, calculating a satellite searching command of an antenna by utilizing initial course and attitude information (course information contains larger errors) of a low-cost inertial navigation system, controlling a pitching axis and a polarization axis of the antenna to rotate according to a calculated pitching command angle and a calculated polarization command angle, then controlling an azimuth axis to scan between 0 and 360 degrees at a constant speed, and finding out an accurate course command angle by a satellite beacon identification technology; then accurately estimating the initial course angle error of the low-cost inertial navigation system by utilizing a three-coordinate calculation and iterative approximation algorithm; and finally, recalculating the satellite finding instruction of the antenna according to the correct inertial navigation course and attitude information, controlling the antenna to rotate according to the instruction angle and accurately aligning the antenna to the satellite. After the antenna initialization is completed by the method, the antenna can be ensured to track the satellite in real time with high precision.
The method comprises the following concrete implementation steps:
firstly, selecting a satellite to be tracked, and calculating a geographic coordinate system O-X of the satellite to be tracked at the location of an antenna according to formula 1tYtZtCoordinate value of (X)ta,Yta,Zta)T
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>X</mi> <mi>ta</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mi>ta</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mi>ta</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mo>-</mo> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>sin</mi> <mi>&Delta;&lambda;</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mi></mi> <mi>&Delta;&lambda;</mi> <mi>sin</mi> <msub> <mi>L</mi> <mi>a</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mi></mi> <mi>&Delta;&lambda;</mi> <mi>cos</mi> <msub> <mi>L</mi> <mi>a</mi> </msub> <mo>-</mo> <mi>R</mi> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> Formula 1
In the formula 1, R is the radius of the earth, and H is the height of a satellite to be tracked; λ ═ λa-λsIs the difference between the longitude of the location of the antenna and the longitude of the location of the satellite, whereaIs the longitude, λ, of the location of the antennasLongitude of the location of the satellite; l isaRepresenting the latitude value at the location of the antenna.
Secondly, after the antenna is electrified and reset, the beam center of the antenna is supposed to point to the geographic coordinate system OY of the antennatIf the axes are consistent, the pitch angle, azimuth angle, and polarization angle of the antenna beam required to point to the satellite can be expressed as:
Figure BSA00000698522400031
formula 2
In formula 2, θtγtRespectively representing the pitching, the azimuth and the polarization tracking command angle of the antenna under the geographic coordinate system.
Thirdly, calculating the coordinate system O-X of the antennaaYaZaPitching command angle theta of lower antenna pointing to satellitepAzimuth command angle
Figure BSA00000698522400033
Polarization command angle gammapThe calculation steps are as follows:
(1) calculating a rotation matrix of the antenna coordinate system (a system) when the geographic coordinate system (t system) is aligned with the satellite
Formula 3
(2) Calculating a rotation matrix from an antenna carrier coordinate system (system b, consistent with an antenna base coordinate system) to a geographic coordinate system
Figure BSA00000698522400036
Figure BSA00000698522400037
Formula 4
In formula 4, θb,γb
Figure BSA00000698522400038
And the pitch angle, the roll angle and the azimuth angle of the antenna carrier are respectively measured by the inertial navigation system. Because the inertial navigation system has lower precision and can not independently find the north
Figure BSA00000698522400039
With a large error from the true azimuth angle, θb,γbThe initial precision of the antenna mainly depends on the accelerometer, generally, the output precision of the antenna is high, the error is small compared with the real attitude angle, and the requirements of the satellite finding and tracking indexes of the antenna can be met.
(3) Solving a rotation matrix between an antenna carrier coordinate system (b system, consistent with an antenna base coordinate system) and an antenna coordinate system when the antenna is aligned with a satellite
Figure BSA00000698522400041
<math> <mrow> <msubsup> <mi>C</mi> <mi>b</mi> <mi>a</mi> </msubsup> <mo>=</mo> <msubsup> <mi>C</mi> <mi>b</mi> <mi>t</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>C</mi> <mi>t</mi> <mi>a</mi> </msubsup> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>c</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>12</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>13</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>22</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>23</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mn>31</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>32</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>33</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> Formula 5
(4) After the antenna is electrified and reset, the antenna coordinate system is superposed with the carrier coordinate system, and at the moment, the calculation instruction angle of the antenna aiming at the satellite to be tracked is as follows:
formula 6
In formula 6, θp,γp
Figure BSA00000698522400044
The antenna pitching, polarization and azimuth star finding command angles are calculated according to the low-precision inertial navigation output. After the inertial navigation system is initially powered on, thetab,γbThe precision is high through the initial alignment;
Figure BSA00000698522400045
by self-north-seeking, because the inertia device has lower precision,
Figure BSA00000698522400046
with large errors. Thus, the antenna satellite finding command angle θp,γp
Figure BSA00000698522400047
There are large errors, and the antenna still cannot be aligned with the satellite after the antenna executes the above instructions.
Fourthly, controlling the pitching of the antenna and the polarization motor according to the command angle thetapAnd gammapThe antenna azimuth control system controls the antenna to slowly rotate for a circle around the azimuth axis at a certain angular speed, monitors the amplitude of a satellite beacon signal received by the antenna all the time in the rotating process, and records the azimuth angle of the antenna relative to the antenna coordinate system zero position at the maximum moment of the satellite signal
Figure BSA00000698522400048
Fifthly, according to
Figure BSA00000698522400049
And calculating the initial course angle of the inertial navigation system by using the following calculation formula:
Figure BSA000006985224000410
formula 7
Wherein,
<math> <mrow> <msubsup> <mi>C</mi> <mi>b</mi> <mi>t</mi> </msubsup> <mo>=</mo> <msubsup> <mi>C</mi> <mi>b</mi> <mi>a</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>C</mi> <mi>a</mi> <mi>t</mi> </msubsup> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>11</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>12</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>c</mi> <mi>bt</mi> <mn>13</mn> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>21</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>22</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>23</mn> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>31</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>32</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>33</mn> </msubsup> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> formula 8
Figure BSA00000698522400052
Formula 9
In formula 9, θp,γpThe satellite finding command angle calculated for equation 6,
Figure BSA00000698522400053
recording the azimuth angle of the antenna at the maximum time of the obtained satellite signal for the fourth step;
Figure BSA00000698522400054
matrix calculated for equation 3
Figure BSA00000698522400055
The transposed matrix of (a), namely:
Figure BSA00000698522400056
sixthly, calculating the course angle of the inertial navigation obtained in the step five
Figure BSA00000698522400057
Substituting into step three, use
Figure BSA00000698522400058
Alternative procedure in formula 4
Figure BSA00000698522400059
Repeating the third step to the fifth step, and calculating to obtain the inertial navigation course angle
Figure BSA000006985224000510
Then the course angle
Figure BSA000006985224000511
Substituting into step three to replace that in step three formula 4
Figure BSA000006985224000512
Repeating the third step to the fifth step, and calculating to obtain the inertial navigation courseCornerBy such iteration and recursion, the initial course angle of inertial navigation can be calculated to beWherein n is more than or equal to 5.
Seven, handleSubstituting into the step three, calculating to obtain the final satellite finding command angle
Figure BSA000006985224000516
The antenna is controlled to rotate according to the command angle, so that the antenna can be accurately aligned with the satellite.
Drawings
Fig. 1 is a schematic diagram of a command angle required by an antenna to align with a satellite in a geographic coordinate system according to the present application.
Fig. 2 is a relationship between a geographical coordinate system and a carrier coordinate system to which the present application relates.
Fig. 3 is a flow chart of initialization of a mobile satellite communication antenna according to the present invention.
DETAILED DESCRIPTION OF EMBODIMENT (S) OF INVENTION
Suppose the satellite to be tracked is a satellite 5 of the middle star with a longitude of 110.5 degrees, and the longitude lambda of the place where the mobile satellite communication antenna is locatedaIs 116 degrees at latitude LaThe angle is 40 degrees, the radius R of the earth is 6378.17km, and the height H of the satellite is 36000 km; the pitch angle, roll angle and course angle of the antenna in the satellite-finding carrier coordinate system (the antenna carrier coordinate system is consistent with the antenna base coordinate system) are respectively assumed as follows: 0 °, 6 °, 0 °; because the low-cost inertial navigation system cannot accurately give a course angle, a pitch angle and a pitch angle can be accurately givenAnd if the roll angle is the following, the pitch angle, roll angle and course angle of the carrier measured by the inertial navigation system are respectively assumed to be: 0 degrees, 6 degrees and 10 degrees, and the course angle error of inertial navigation output is 10 degrees.
According to the formula 1, the coordinate values of the satellite in the geographic coordinate system can be calculated as follows:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>X</mi> <mi>ta</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mi>ta</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mi>ta</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mo>-</mo> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>sin</mi> <mi>&Delta;&lambda;</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mi></mi> <mi>&Delta;&lambda;</mi> <mi>sin</mi> <msub> <mi>L</mi> <mi>a</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mi></mi> <mi>&Delta;&lambda;</mi> <mi>cos</mi> <msub> <mi>L</mi> <mi>a</mi> </msub> <mo>-</mo> <mi>R</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mo>-</mo> <mn>4061.76</mn> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mn>27114.75</mn> </mtd> </mtr> <mtr> <mtd> <mn>25935.94</mn> </mtd> </mtr> </mtable> </mfenced> <mi>km</mi> </mrow> </math>
the instruction calculated according to equation 2 is:
Figure BSA00000698522400062
the satellite finding command angle calculated according to formula 6 is:
Figure BSA00000698522400063
the rotation angle of the antenna corresponding to the time when the satellite signal is maximum when the antenna scans for 0-360 degrees can be calculated according to the real positions of the satellite and the antenna
Figure BSA00000698522400064
Has a value of 194.68 degrees.
The calculation according to equation 7 of step five can be obtained:
Figure BSA00000698522400065
after iteration, it can be calculated that:
Figure BSA00000698522400072
Figure BSA00000698522400073
from the above calculation results, after 3 iterations, the error between the azimuth angle of the inertial navigation system and the true azimuth angle is better than 0.005 °, so that the estimated inertial navigation azimuth angle can be ensured to be consistent with the true azimuth angle under the condition that the iteration number n is greater than or equal to 5.

Claims (2)

1. A method for initializing a mobile satellite communication antenna is characterized in that: the method comprises the following implementation steps:
the first step is as follows: selecting a satellite to be tracked, and calculating the coordinate value of the satellite to be tracked in the antenna geographic coordinate system;
the second step is that: calculating a pitch angle, an azimuth angle and a polarization tracking instruction angle of the geographic coordinate system pointing to the satellite to be rotated;
the third step: calculating a pitching instruction angle, an azimuth instruction angle and a polarization instruction angle of the antenna beam pointing to the satellite to be rotated;
the fourth step: calculating to obtain an initial course angle of an inertial navigation system installed on the antenna by combining an antenna 0-360-degree scanning technology and an iterative recursion algorithm according to the pitching instruction angle, the azimuth instruction angle and the polarization instruction angle obtained in the third step;
the fifth step: recalculating an antenna satellite finding instruction angle according to the initial course angle of the inertial navigation system obtained in the fourth step;
and a sixth step: and controlling the antenna to align the satellite according to the satellite searching command angle obtained in the fifth step.
2. A method for initializing a mobile satellite communication antenna is characterized in that the method comprises the following implementation steps:
the first step is as follows: selecting a satellite to be tracked, and calculating a geographic coordinate system O-X of the satellite to be tracked at the location of the antenna according to the formula 1tYtZtCoordinate value of (X)ta,Yta,Zta)T
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>X</mi> <mi>ta</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mi>ta</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mi>ta</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mo>-</mo> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>sin</mi> <mi>&Delta;&lambda;</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mi></mi> <mi>&Delta;&lambda;</mi> <mi>sin</mi> <msub> <mi>L</mi> <mi>a</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mi></mi> <mi>&Delta;&lambda;</mi> <mi>cos</mi> <msub> <mi>L</mi> <mi>a</mi> </msub> <mo>-</mo> <mi>R</mi> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> Formula 1
In the formula 1, R is the radius of the earth, and H is the height of a satellite to be tracked; λ ═ λasIs the difference between the longitude of the location of the antenna and the longitude of the location of the satellite, whereaIs the longitude, λ, of the location of the antennasLongitude of the location of the satellite; l isaRepresenting the latitude value at the location of the antenna.
The second step is that: after the antenna is powered on and reset, the beam center of the antenna is supposed to point to the geographic coordinate system OY of the antennatIf the axes are consistent, the pitch angle, azimuth angle, and polarization angle of the antenna beam required to point to the satellite can be expressed as:
Figure FSA00000698522300021
formula 2
In formula 2, θt
Figure FSA00000698522300022
γtRespectively representing the pitching, the azimuth and the polarization tracking command angle of the antenna under the geographic coordinate system. The third step: calculating the antenna coordinate system O-XaYaZaPitching command angle theta of lower antenna pointing to satellitepAzimuth command angle
Figure FSA00000698522300023
Polarization command angle gammapThe calculation steps are as follows:
(1) calculating a rotation matrix of the antenna coordinate system (a system) when the geographic coordinate system (t system) is aligned with the satellite
Figure FSA00000698522300024
Figure FSA00000698522300025
Formula 3
(2) Calculating a rotation matrix from an antenna carrier coordinate system (system b, consistent with an antenna base coordinate system) to a geographic coordinate system
Figure FSA00000698522300026
Figure FSA00000698522300027
Formula 4
In formula 4, θb,γb
Figure FSA00000698522300028
And the pitch angle, the roll angle and the azimuth angle of the antenna carrier are respectively measured by the inertial navigation system. Because the inertial navigation system has lower precision and can not independently find the north
Figure FSA00000698522300029
With a large error from the true azimuth angle, θb,γbThe initial precision of the antenna mainly depends on the accelerometer, generally, the output precision of the antenna is high, the error is small compared with the real attitude angle, and the requirements of the satellite finding and tracking indexes of the antenna can be met.
(3) Solving antenna carrier seatRotation matrix between the system of coordinates (b, corresponding to the antenna base coordinate system) to the antenna coordinate system when the antenna is aligned with the satellite
Figure FSA00000698522300031
<math> <mrow> <msubsup> <mi>C</mi> <mi>b</mi> <mi>a</mi> </msubsup> <mo>=</mo> <msubsup> <mi>C</mi> <mi>b</mi> <mi>t</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>C</mi> <mi>t</mi> <mi>a</mi> </msubsup> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>c</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>12</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>13</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>22</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>23</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mn>31</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>32</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>33</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> Formula 5
(4) After the antenna is electrified and reset, the antenna coordinate system is superposed with the carrier coordinate system, and at the moment, the calculation instruction angle of the antenna aiming at the satellite to be tracked is as follows:
formula 6
In formula 6, θp,γpThe antenna pitching, polarization and azimuth star finding command angles are calculated according to the low-precision inertial navigation output. After the inertial navigation system is initially powered on, thetab,γbThe precision is high through the initial alignment;
Figure FSA00000698522300035
by self-north-seeking, because the inertia device has lower precision,
Figure FSA00000698522300036
with large errors. Thus, the antenna satellite finding command angle θp,γp
Figure FSA00000698522300037
There are large errors, and the antenna still cannot be aligned with the satellite after the antenna executes the above instructions.
The fourth step: controlling the pitching and polarizing motors of the antenna according to the command angle thetapAnd gammapThe antenna azimuth control system controls the antenna to slowly rotate for a circle around the azimuth axis at a certain angular speed, monitors the amplitude of a satellite beacon signal received by the antenna all the time in the rotating process, and records the azimuth angle of the antenna relative to the antenna coordinate system zero position at the maximum moment of the satellite signal
Figure FSA00000698522300038
The fifth step: according to
Figure FSA00000698522300039
And calculating the initial course angle of the inertial navigation system by using the following calculation formula:
formula 7
Wherein,
<math> <mrow> <msubsup> <mi>C</mi> <mi>b</mi> <mi>t</mi> </msubsup> <mo>=</mo> <msubsup> <mi>C</mi> <mi>b</mi> <mi>a</mi> </msubsup> <mo>&times;</mo> <msubsup> <mi>C</mi> <mi>a</mi> <mi>t</mi> </msubsup> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>11</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>12</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>c</mi> <mi>bt</mi> <mn>13</mn> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>21</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>22</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>23</mn> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>31</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>32</mn> </msubsup> </mtd> <mtd> <msubsup> <mi>C</mi> <mi>bt</mi> <mn>33</mn> </msubsup> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> formula 8
Figure FSA00000698522300042
Formula 9
In formula 9, θp,γpThe satellite finding command angle calculated for equation 6,
Figure FSA00000698522300043
recording the azimuth angle of the antenna at the maximum time of the obtained satellite signal for the fourth step;
Figure FSA00000698522300044
matrix calculated for equation 3The transposed matrix of (a), namely:
Figure FSA00000698522300046
and a sixth step: calculating the course angle of the inertial navigation obtained in the step fiveSubstituting into step three, useAlternative procedure in formula 4Repeating the third step to the fifth step, and calculating to obtain the inertial navigation course angle
Figure FSA000006985223000410
Then the course angle
Figure FSA000006985223000411
Substituting into step three to replace that in step three formula 4
Figure FSA000006985223000412
Repeating the third step to the fifth step, and calculating to obtain the inertial navigation course angle
Figure FSA000006985223000413
By such iteration and recursion, the initial course angle of inertial navigation can be calculated to be
Figure FSA000006985223000414
Wherein n is more than or equal to 5.
The seventh step: handle
Figure FSA000006985223000415
Substituting into the step three, calculating to obtain the final satellite finding command angle
Figure FSA000006985223000416
The antenna is controlled to rotate according to the command angle, so that the antenna can be accurately aligned with the satellite.
CN2012101035437A 2012-04-11 2012-04-11 Initializing method of mobile satellite communication antenna Pending CN102662188A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012101035437A CN102662188A (en) 2012-04-11 2012-04-11 Initializing method of mobile satellite communication antenna

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012101035437A CN102662188A (en) 2012-04-11 2012-04-11 Initializing method of mobile satellite communication antenna

Publications (1)

Publication Number Publication Date
CN102662188A true CN102662188A (en) 2012-09-12

Family

ID=46771715

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012101035437A Pending CN102662188A (en) 2012-04-11 2012-04-11 Initializing method of mobile satellite communication antenna

Country Status (1)

Country Link
CN (1) CN102662188A (en)

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104124528A (en) * 2014-05-05 2014-10-29 北京星网卫通科技开发有限公司 Inertia/GNSS (Global Navigation Satellite System)/satellite beacon based integrated communication on the move antenna stabilization tracking method
CN104125003A (en) * 2013-04-25 2014-10-29 成都国恒空间技术工程有限公司 Rapid satellite locking system of portable satellite receiver
CN104913790A (en) * 2015-05-28 2015-09-16 北京航天控制仪器研究所 Heading drift error closed-loop compensation method applied to communication-in-moving inertial navigation system
CN105021191A (en) * 2015-07-07 2015-11-04 中国人民解放军第二炮兵工程大学 Low-cost method for estimating antenna attitude of Satcom on the move measurement and control system
CN105116430A (en) * 2015-08-21 2015-12-02 北京航天控制仪器研究所 SOTM pseudo course sea ship dynamic satellite searching method based on Kalman filtering
CN105789894A (en) * 2016-04-11 2016-07-20 中国人民解放军国防科学技术大学 UHF antenna tracking method of satellite in sun-synchronous orbit
CN106092096A (en) * 2016-06-03 2016-11-09 上海航天控制技术研究所 In high-precision orbital emulation, the satellite position based on iterative approach method determines method
CN107579759A (en) * 2017-09-19 2018-01-12 清华大学 Method and device for stabilizing antenna beam in UAV satellite communication system
CN109149110A (en) * 2018-08-31 2019-01-04 捷信(浙江)通信技术有限公司 A kind of satellite Dynamic Tracking and antenna equipment
CN109302223A (en) * 2018-09-12 2019-02-01 上海无线电设备研究所 The antenna selecting method of group-net communication between multiple high dynamic carriers
CN109443385A (en) * 2018-11-13 2019-03-08 中国兵器装备集团自动化研究所 A kind of inertial navigation installation error automatic calibration method of antenna for satellite communication in motion
CN109582045A (en) * 2019-01-08 2019-04-05 北京慧清科技有限公司 The Initial Alignment Method of antenna when a kind of carrier inclined
CN110441797A (en) * 2019-08-07 2019-11-12 中国海洋大学 Highly reliable Beidou RDSS antenna automatic stabilisation tracking based on carrier-to-noise ratio
CN110764119A (en) * 2019-11-07 2020-02-07 中国人民解放军火箭军工程大学 Satellite antenna autonomous measurement and control method and system independent of satellite navigation information
CN110808447A (en) * 2019-10-24 2020-02-18 迪泰(浙江)通信技术有限公司 Shipborne satellite antenna system based on triaxial dynamic tracking technology
CN111064002A (en) * 2018-10-16 2020-04-24 成都空间矩阵科技有限公司 Servo control method for low-profile satellite communication antenna
CN111337055A (en) * 2020-05-07 2020-06-26 成都国卫通信技术有限公司 Calibration method for satellite mobile communication antenna inertial navigation
CN111864347A (en) * 2020-06-24 2020-10-30 宁波大学 A polarization dynamic matching method for VICTS antenna
CN111912404A (en) * 2019-05-09 2020-11-10 西安京东天鸿科技有限公司 Output attitude correction system and method for flight equipment
CN112350766A (en) * 2020-10-26 2021-02-09 中国电子科技集团公司第五十四研究所 Angle control system and method for antenna pointing to low-orbit communication satellite
CN111864348B (en) * 2020-06-24 2021-03-09 宁波大学 Initial satellite finding method of VICTS antenna
CN113849003A (en) * 2021-10-13 2021-12-28 西安尹纳数智能科技有限公司 Control method for motion isolation of communication-in-motion antenna
CN114583451A (en) * 2022-04-28 2022-06-03 成都迅翼卫通科技有限公司 Satellite tracking method, device and system for communication in motion and computer equipment
CN115987371A (en) * 2022-12-05 2023-04-18 迪泰(浙江)通信技术有限公司 Satellite loss identification method and device for satellite mobile terminal antenna

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020092350A1 (en) * 2000-11-28 2002-07-18 Business Arts Inc. Gravity gradiometry
US20020116125A1 (en) * 2000-12-23 2002-08-22 Ching-Fang Lin Positioning & data integrating method and system thereof
CN101916118A (en) * 2010-07-13 2010-12-15 北京爱科迪信息通讯技术有限公司 Method and system for seeking stars by using satellite antenna
US20110172850A1 (en) * 2009-09-14 2011-07-14 Israel Aerospace Industries Ltd. Infantry robotic porter system and methods useful in conjunction therewith
CN102486377A (en) * 2009-11-17 2012-06-06 哈尔滨工程大学 A Method for Attitude Acquisition of Initial Heading of Fiber Optic Gyro Strapdown Inertial Navigation System

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020092350A1 (en) * 2000-11-28 2002-07-18 Business Arts Inc. Gravity gradiometry
US20020116125A1 (en) * 2000-12-23 2002-08-22 Ching-Fang Lin Positioning & data integrating method and system thereof
US20110172850A1 (en) * 2009-09-14 2011-07-14 Israel Aerospace Industries Ltd. Infantry robotic porter system and methods useful in conjunction therewith
CN102486377A (en) * 2009-11-17 2012-06-06 哈尔滨工程大学 A Method for Attitude Acquisition of Initial Heading of Fiber Optic Gyro Strapdown Inertial Navigation System
CN101916118A (en) * 2010-07-13 2010-12-15 北京爱科迪信息通讯技术有限公司 Method and system for seeking stars by using satellite antenna

Cited By (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104125003A (en) * 2013-04-25 2014-10-29 成都国恒空间技术工程有限公司 Rapid satellite locking system of portable satellite receiver
CN104125003B (en) * 2013-04-25 2018-05-29 成都国恒空间技术工程有限公司 A kind of portable satellite receiver satellite quick lock systems
CN104124528A (en) * 2014-05-05 2014-10-29 北京星网卫通科技开发有限公司 Inertia/GNSS (Global Navigation Satellite System)/satellite beacon based integrated communication on the move antenna stabilization tracking method
CN104124528B (en) * 2014-05-05 2016-03-02 北京星网卫通科技开发有限公司 Exceedingly high line stabilization tracking in a kind of inertia/GNSS/ satellite beacon combined moving
CN104913790B (en) * 2015-05-28 2017-11-28 北京航天控制仪器研究所 A kind of inertial navigation system heading effect error closed loop compensation method applied to communication in moving
CN104913790A (en) * 2015-05-28 2015-09-16 北京航天控制仪器研究所 Heading drift error closed-loop compensation method applied to communication-in-moving inertial navigation system
CN105021191A (en) * 2015-07-07 2015-11-04 中国人民解放军第二炮兵工程大学 Low-cost method for estimating antenna attitude of Satcom on the move measurement and control system
CN105116430B (en) * 2015-08-21 2017-06-27 北京航天万达高科技有限公司 The sea pool state based on Kalman filtering for the pseudo- course of communication in moving searches star method
CN105116430A (en) * 2015-08-21 2015-12-02 北京航天控制仪器研究所 SOTM pseudo course sea ship dynamic satellite searching method based on Kalman filtering
CN105789894A (en) * 2016-04-11 2016-07-20 中国人民解放军国防科学技术大学 UHF antenna tracking method of satellite in sun-synchronous orbit
CN106092096A (en) * 2016-06-03 2016-11-09 上海航天控制技术研究所 In high-precision orbital emulation, the satellite position based on iterative approach method determines method
CN107579759A (en) * 2017-09-19 2018-01-12 清华大学 Method and device for stabilizing antenna beam in UAV satellite communication system
CN109149110A (en) * 2018-08-31 2019-01-04 捷信(浙江)通信技术有限公司 A kind of satellite Dynamic Tracking and antenna equipment
CN109149110B (en) * 2018-08-31 2021-09-03 捷信(浙江)通信技术有限公司 Satellite dynamic tracking method and antenna equipment
CN109302223B (en) * 2018-09-12 2021-11-02 上海无线电设备研究所 Antenna selection method for networking communication among multiple high dynamic carriers
CN109302223A (en) * 2018-09-12 2019-02-01 上海无线电设备研究所 The antenna selecting method of group-net communication between multiple high dynamic carriers
CN111064002A (en) * 2018-10-16 2020-04-24 成都空间矩阵科技有限公司 Servo control method for low-profile satellite communication antenna
CN111064002B (en) * 2018-10-16 2021-04-23 正成集团科技有限公司 Servo control method for low-profile satellite communication antenna
CN109443385A (en) * 2018-11-13 2019-03-08 中国兵器装备集团自动化研究所 A kind of inertial navigation installation error automatic calibration method of antenna for satellite communication in motion
CN109443385B (en) * 2018-11-13 2022-07-29 中国兵器装备集团自动化研究所有限公司 Inertial navigation installation error automatic calibration method of communication-in-moving antenna
CN109582045A (en) * 2019-01-08 2019-04-05 北京慧清科技有限公司 The Initial Alignment Method of antenna when a kind of carrier inclined
CN111912404B (en) * 2019-05-09 2024-04-05 西安京东天鸿科技有限公司 Output attitude correction system and method for flight equipment
CN111912404A (en) * 2019-05-09 2020-11-10 西安京东天鸿科技有限公司 Output attitude correction system and method for flight equipment
CN110441797A (en) * 2019-08-07 2019-11-12 中国海洋大学 Highly reliable Beidou RDSS antenna automatic stabilisation tracking based on carrier-to-noise ratio
CN110808447A (en) * 2019-10-24 2020-02-18 迪泰(浙江)通信技术有限公司 Shipborne satellite antenna system based on triaxial dynamic tracking technology
CN110764119A (en) * 2019-11-07 2020-02-07 中国人民解放军火箭军工程大学 Satellite antenna autonomous measurement and control method and system independent of satellite navigation information
CN111337055A (en) * 2020-05-07 2020-06-26 成都国卫通信技术有限公司 Calibration method for satellite mobile communication antenna inertial navigation
CN111337055B (en) * 2020-05-07 2023-06-02 成都国卫通信技术有限公司 Calibration method for satellite mobile communication antenna inertial navigation
CN111864348B (en) * 2020-06-24 2021-03-09 宁波大学 Initial satellite finding method of VICTS antenna
CN111864347A (en) * 2020-06-24 2020-10-30 宁波大学 A polarization dynamic matching method for VICTS antenna
CN112350766A (en) * 2020-10-26 2021-02-09 中国电子科技集团公司第五十四研究所 Angle control system and method for antenna pointing to low-orbit communication satellite
CN113849003A (en) * 2021-10-13 2021-12-28 西安尹纳数智能科技有限公司 Control method for motion isolation of communication-in-motion antenna
CN113849003B (en) * 2021-10-13 2024-04-26 复远芯(上海)科技有限公司 Control method for motion isolation of communication-in-motion antenna
CN114583451A (en) * 2022-04-28 2022-06-03 成都迅翼卫通科技有限公司 Satellite tracking method, device and system for communication in motion and computer equipment
CN114583451B (en) * 2022-04-28 2022-07-22 成都迅翼卫通科技有限公司 Satellite tracking method, device and system for communication in motion and computer equipment
CN115987371A (en) * 2022-12-05 2023-04-18 迪泰(浙江)通信技术有限公司 Satellite loss identification method and device for satellite mobile terminal antenna
CN115987371B (en) * 2022-12-05 2023-12-05 迪泰(浙江)通信技术有限公司 Satellite loss identification method and device for satellite mobile terminal antenna

Similar Documents

Publication Publication Date Title
CN102662188A (en) Initializing method of mobile satellite communication antenna
US10355351B2 (en) Antenna array pointing direction estimation and control
CN102508275B (en) Multiple-antenna GPS(Global Positioning System)/GF-INS (Gyroscope-Free-Inertial Navigation System) depth combination attitude determining method
CN104124528B (en) Exceedingly high line stabilization tracking in a kind of inertia/GNSS/ satellite beacon combined moving
CN101893440B (en) Celestial autonomous navigation method based on star sensors
CN108508918B (en) A high-precision real-time ground pointing control method for a static-orbit remote sensing satellite data transmission antenna
CN110926468B (en) Communication-in-motion antenna multi-platform navigation attitude determination method based on transfer alignment
CN103633417B (en) Airborne antenna based on strapdown attitude tenacious tracking high accuracy points to tracking
CN113237456B (en) Method for measuring initial installation angle of communication-in-motion antenna
CN108759819A (en) A kind of polarization navigation real-time location method based on omnimax polarization degree information
CN103759727B (en) Navigation and positioning method based on sky polarized light distribution mode
CN111864348B (en) Initial satellite finding method of VICTS antenna
CN101762273A (en) Autonomous optical navigation method for soft landing for deep space probe
CN106410410A (en) Satellite capturing and tracking method for VSAT (Very Small Aperture Terminal) antenna system with physical level platform
CN105116430B (en) The sea pool state based on Kalman filtering for the pseudo- course of communication in moving searches star method
CN103913163A (en) Method for calculating A-E-C triaxial antenna coordinate of shipborne satellite communication earth station
CN107576939A (en) A kind of single beacon distance-measuring and positioning method based on virtual ranging beacon
CN104880192A (en) Carrier course angle calculation method based on polarization compass
CN103972654A (en) Communication-in-moving antenna satellite alignment tracking device under shielding of rotor wings of helicopter
CN106005455B (en) A kind of two axis Towed bird systems for being directed toward control based on geographic coordinate system
CN103985952A (en) Shipborne A-E-C triaxial satellite communication antenna polarization deviation angle real-time correction method
CN106329120A (en) Sitcom On-the-Move (SOTM) low-cost measurement and control method
CN103471614A (en) Transfer alignment method in polar region based on inverse coordinate system
CN113108788B (en) A long-endurance inertial navigation/astronomical global integrated navigation method
CN107830873A (en) A kind of high-precision vehicle positioning and orienting method combined based on uniaxial horizontal rotation inertial navigation with odometer

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120912