Method for calculating GEO satellite beam pointing deviation based on maneuver measurement analysis

Abstract

This invention relates to a method for calculating the beam pointing deviation of GEO satellites based on mobile measurement analysis, belonging to the field of satellite communication technology. This invention utilizes power data from four calibration beams measured by a calibration station to calculate the elevation and roll deviation angles between the pointing angle from satellite space point S to calibration station position point B and the pointing angle from satellite space point S to the actual beam pointing center point C. Then, using the spatial position information of points A, B, and C, the elevation and roll deviation angles between the pointing angle from satellite space point S to calibration station position point B and the pointing angle from satellite space point S to the theoretical beam pointing center point A are calculated. Finally, the satellite beam pointing deviation is calculated. This invention is designed for GEO fixed-beam communication satellites, using measurement signals from a mobile calibration station to calculate the beam pointing deviation, and can be used to assist in satellite beam pointing correction. This invention supports mobile platform deployment, reducing the difficulty and cost of site construction compared to fixed-site deployment.

Description

Technical Field

[0001] This invention relates to the field of satellite communication technology, and in particular to a method for calculating the beam pointing deviation of GEO satellites based on maneuvering measurement and analysis. This method is applicable to fixed beams of geosynchronous orbit satellites and can assist satellites in adjusting the beam pointing to the theoretical position. Background Technology

[0002] For GEO fixed-beam communication satellites, the system needs to maintain beam pointing accuracy to facilitate network planning and generate precise planning schemes, assisting users in switching services based on coverage area. When the satellite coverage area is within the country's territory, fixed ground calibration stations can be built to measure and calculate the satellite antenna pointing deviation, implement beam pointing control, and maintain the stability of the satellite beam service range.

[0003] However, for GEO satellites with coverage areas overseas, constructing fixed calibration stations is difficult. There is a need to propose a method for calculating GEO satellite beam pointing deviation based on mobile measurement and analysis, enabling calibration stations to be deployed on mobile platforms. Currently, such a method does not exist in existing technology. Summary of the Invention

[0004] The purpose of this invention is to provide a method for calculating GEO satellite beam pointing deviation based on mobile measurement and analysis. Addressing the need for mobile deployment of satellite beam calibration stations, this invention proposes an innovative method for calculating GEO satellite beam pointing deviation, which supports the deployment of mobile platforms for calibration stations.

[0005] The technical solution adopted in this invention is as follows:

[0006] A method for calculating GEO satellite beam pointing deviation based on mobile measurement analysis, applied to beam calibration stations, includes the following steps:

[0007] Step 1: Obtain the longitude of the satellite's theoretical orbit determination point S. The longitude of the theoretical center point A of the satellite beam on the ground. and latitude Load the beam pointing deviation power design dataset for four cardinal directions (north, south, east, and west). The data in the dataset are... ,in For pitch angle deviation, For the roll angle deviation, For the beam power value of the East Standard, To calibrate the beam power value for the West standard, To calibrate the beam power value for the South Standard, The beam power value is for the North Standard calibration.

[0008] Step 2, based on the data Calculate the pointing deviation data after power normalization. ,in, This is the normalized value for the east-west calibration beam power. The normalized value of the north-south calibration beam power is calculated using the following formula:

[0009]

[0010] Step 3, generate the pointing deviation judgment matrix D:

[0011] ,

[0012] in, express , , That is, the 0th row contains 101 data points with a roll angle deviation of -0.5°, and the roll angle of each subsequent row increases by 0.01°; the 0th column contains 101 data points with a pitch angle deviation of -0.5°, and the pitch angle of each subsequent column increases by 0.01°.

[0013] Step 4: The beam calibration station receives calibration signals from the four calibration beams of the satellite and generates measurement data. ,in, The east, west, south, and north beam power values ​​measured at the beam calibration station. To measure time, To determine the longitude of the station, To determine the latitude of the calibration station; calculate the normalized value of the east-west beam power. and north and south beam normalization values ;

[0014] Step 5, using Calculate the pitch deviation line using the pointing deviation judgment matrix D;

[0015] Step 6, using The rolling angle deviation line is calculated using the pointing deviation judgment matrix D;

[0016] Step 7: Calculate the intersection of the pitch deviation line and the roll deviation line. , and That is, the pitch deviation angle and roll deviation angle between the pointing angle from satellite space point S to calibration station position point B and the pointing angle from satellite space point S to the actual beam pointing center point C;

[0017] Step 8: Based on the position of the satellite's theoretical orbit determination point S The theoretical location of the satellite beam pointing to its center point A on the ground. He Biao School Station Location Calculate the elevation deviation angle of the pointing angle from satellite space point S to calibration station position point B relative to the pointing angle from satellite space point S to the theoretical beam pointing center point A. and rolling deviation angle ;

[0018] Step 9: Based on the calculation results of Steps 7 and 8, calculate the elevation deviation angle of the pointing angle from the satellite space point S to the actual beam pointing center point C relative to the pointing angle from the satellite space point S to the theoretical beam pointing center point A. and rolling deviation angle The calculation result is the satellite beam pointing deviation.

[0019] Furthermore, step 5 is specifically implemented as follows:

[0020] Step 5-1: Select the deviation range [-0.5°, 0.5°], where the power of the west calibration beam monotonically decreases from west to east, and the power of the east calibration beam monotonically increases. Pointing deviation judgment matrix. The beam normalization values ​​of each row are monotonically increasing;

[0021] Step 5-2, determine the matrix based on the pointing deviation. Take out the first one If we disregard the normalized power values ​​of the north-south calibrated beams in the row of data, then any data point in that row is... ,in For column number, from arrive The values ​​are monotonically increasing; search for the first value starting from the first column. Value greater than The column, column number is ,if Or there is no greater than If a column is not found, it indicates that the row has no intersection points; otherwise, the column number is calculated according to the following formula. Intersection of line data and pitch deviation line :

[0022]

[0023]

[0024] Step 5-3: Repeat step 5-2 to calculate the intersection points from row 0 to row 100, i.e., the pitch deviation points. The line connecting the intersection points forms the pitch angle deviation line.

[0025] Furthermore, step 6 is performed as follows:

[0026] Step 6-1: Select the deviation range [-0.5°, 0.5°], where the power of the north calibration beam decreases monotonically from north to south, and the power of the south calibration beam increases monotonically. (Pointing deviation judgment matrix) The normalized values ​​of the north and south beams in each column decrease monotonically.

[0027] Step 6-2, determine the matrix based on the pointing deviation. Take out the first one If we disregard the normalized power values ​​of the east-west calibrated beams in this column, then any data point in this column is... ,in For line numbers, from arrive The value is monotonically decreasing; search for the first value starting from the first row. Value less than The line number is ,if Or there is no less than If a row has no intersection, it means that the row has no intersection points; otherwise, calculate the nth row according to the following formula. Intersection of column data and rolling deviation line :

[0028]

[0029]

[0030] Step 6-3: Repeat step 6-2 to calculate the intersection points from column 0 to column 100, i.e., the rolling deviation points. The line connecting the intersection points forms the roll angle deviation line.

[0031] Furthermore, step 7 is specifically implemented as follows:

[0032] Step 7-1: Obtain the pitch deviation point from Step 5 and the roll deviation point from Step 6;

[0033] Step 7-2, obtain the first pitch deviation point coordinates and the pitch deviation point coordinates ;

[0034] Step 7-3, obtain the first One rolling deviation point coordinates and the One rolling deviation point coordinates ;

[0035] Step 7-4, determine the first pitch deviation point and the pitch deviation point Does the connection match the first? One rolling deviation point and the One rolling deviation point The lines intersect, and the following method is used to determine if they do:

[0036] Step 7-4-1: Determine whether the two line segments coincide on both the pitch deviation angle axis and the roll deviation angle axis. If not, there can be no intersection point. Proceed to step 7-5. If yes, proceed to step 7-4-2.

[0037] Step 7-4-2: If both of the following discriminants are true, then the pair intersects; otherwise, they do not intersect. The discriminants are as follows:

[0038]

[0039]

[0040] Step 7-5: If no intersection point exists, take the next rolling deviation point and repeat steps 7-3 to 7-5. If the two lines intersect, the calculation is complete. Proceed to step 8; the formula is as follows:

[0041]

[0042]

[0043] in, As an intermediate variable, the calculation method is as follows:

[0044] ;

[0045] Steps 7-6, if the first The and the first If the line connecting each pitch deviation point does not intersect with the line connecting each roll deviation point, then take the next pitch point and repeat steps 7-2 to 7-6.

[0046] Step 7-7: If the lines connecting each pitch deviation point and the lines connecting each roll deviation point do not intersect, the calculation accuracy is exceeded, and the satellite beam pointing deviation cannot be calculated.

[0047] Furthermore, step 8 is specifically implemented as follows:

[0048] Step 8-1, satellite spatial point S is The beam theory points to the center point A as The location of the marker station is point B. The unit for latitude and longitude is °. The altitude of the GEO satellite, For the Earth's radius, and The value range is [0, 60];

[0049] Step 8-2: Calculate the pointing angle from satellite space point S to calibration station location point B. The formula is as follows:

[0050]

[0051]

[0052] Step 8-3: Calculate the pointing angle from satellite spatial point S to the theoretical beam pointing center point A. The formula is as follows:

[0053]

[0054]

[0055] Step 8-4: Based on the results of steps 8-2 and 8-3, calculate... The formula is as follows:

[0056]

[0057] .

[0058] Furthermore, step 9 is specifically implemented as follows:

[0059] Step 9-1, calculate the beam pointing elevation deviation angle:

[0060]

[0061] Step 9-2, calculate the beam pointing roll deviation angle:

[0062]

[0063] Step 9-3: Obtain the satellite beam pointing deviation angle. .

[0064] Compared with the prior art, the present invention has the following advantages:

[0065] 1. This invention supports the deployment of satellite beam calibration stations on mobile platforms;

[0066] 2. This invention can assist in adjusting the satellite beam pointing, which helps maintain stable satellite coverage.

[0067] 3. The algorithm of this invention is simple, has low implementation complexity, and has value for widespread application. Attached Figure Description

[0068] Figure 1 This is a schematic diagram illustrating the principle of the satellite beam pointing deviation calculation method for mobile calibration stations.

[0069] Figure 2 This diagram illustrates the calculation of the actual beam pointing of a satellite relative to the elevation and roll angle deviations from the satellite's pointing direction at the calibration station, using measurement data from the calibration station.

[0070] Figure 3 This is a flowchart of a method for calculating the beam pointing deviation of GEO satellites based on motion measurement analysis. Detailed Implementation

[0071] The present invention will now be described in detail with reference to the accompanying drawings.

[0072] A method for calculating GEO satellite beam pointing deviation based on maneuvering measurement and analysis is proposed. This method uses power data from four calibration beams measured by a calibration station to calculate the pitch and roll deviation angles between the pointing angle from satellite space point S to calibration station position point B and the pointing angle from satellite space point S to the actual beam pointing center point C. Then, using the spatial position information of points A, B, and C, calculate the pitch and roll deviation angles of the pointing angle from satellite spatial point S to calibration station position point B relative to the pointing angle from satellite spatial point S to the theoretical beam pointing center point A. Finally, utilize and Calculate satellite beam pointing deviation .

[0073] Assume point A is the theoretical center of the satellite beam, point B is the location of the beam calibration station, point C is the actual center of the satellite beam, and point S is the satellite position. Figure 1 The image shows three points A, B, and C on the ground.

[0074] like Figure 3 As shown, the method specifically includes the following steps:

[0075] Step (1): Initialize the calibration system and obtain the longitude of the satellite's theoretical orbit determination point S. The longitude of the theoretical center point A of the satellite beam on the ground. and latitude Load the beam pointing deviation power design dataset for four cardinal directions from the database. ,in For pitch angle deviation, For the roll angle deviation, For the beam power value of the East Standard, To calibrate the beam power value for the West standard, To calibrate the beam power value for the South Standard, The beam power values ​​for the North Standard calibration are calculated in 0.01° increments, with pitch and roll angle errors within the range of [-0.5°, 0.5°]. The pointing deviation power design dataset contains a total of 10201 data sets.

[0076] Step (2): Based on the data Calculate the pointing deviation data after power normalization. ,in This is the normalized value for the east-west calibration beam power. The normalized power value for the north-south calibration beam is given, with a power normalization range of (-1, 1). The calculation formula is as follows:

[0077]

[0078] Step (3): Use the calculation results of step (2) to generate the pointing deviation judgment matrix. ,in, express , , That is, the 0th row contains 101 data points with a roll angle deviation of -0.5°, and the roll angle of each subsequent row increases by 0.01°; the 0th column contains 101 data points with a pitch angle deviation of -0.5°, and the pitch angle of each subsequent column increases by 0.01°.

[0079] Step (4): The mobile calibration station receives calibration signals transmitted by the four calibration beams of the satellite and generates measurement data. ,in To measure time, To determine the longitude of the station, To determine the latitude of the calibration station, the normalized value of the east-west beam power is calculated according to the formula in step (2). and north and south beam normalization values ;

[0080] Step (5): Using the information from step (4) After power normalization in step (3), the pitch deviation matrix D is used to calculate the pitch deviation line, which is composed of multiple points. The pitch deviation line includes the following steps:

[0081] Step (5-1): Within the selected deviation range [-0.5°, 0.5°], the power of the west calibration beam decreases monotonically from west to east, while the power of the east calibration beam increases monotonically. The pointing deviation judgment matrix... The beam normalization values ​​of each row are monotonically increasing;

[0082] Step (5-2): From the matrix Take out the first one row data, ,in For column number, from arrive The values ​​are monotonically increasing; start searching from the first column for the first value. Value greater than The column, column number is ,if or no greater than If the column is empty, it indicates that the row has no intersection; otherwise, the column is empty. Intersection of line data and pitch deviation line The calculation formula is as follows:

[0083]

[0084]

[0085] Step (5-3): Repeat step (5-2) to calculate the intersection points from row 0 to row 100, i.e. The line connecting the intersection points forms the pitch angle deviation line.

[0086] Step (6): Utilize the information from step (4) After power normalization in step (3), the roll angle deviation line is calculated from the deviation matrix D, which consists of multiple points. The roll angle deviation line includes the following steps:

[0087] Step (6-1): Within the selected deviation range [-0.5°, 0.5°], the power of the north calibration beam decreases monotonically from north to south, while the power of the south calibration beam increases monotonically. The pointing deviation judgment matrix... The normalized values ​​of the north and south beams in each column decrease monotonically.

[0088] Step (6-2): From the matrix Take out the first one Column data, ,in For column number, from arrive The value is monotonically decreasing, so we search for the first one starting from the first row. Value less than The line number is ,if or no less than If the line is true, it means that the line has no intersection; otherwise... Intersection of column data and rolling deviation line The calculation formula is as follows:

[0089]

[0090]

[0091] Step (6-3): Repeat step (6-2) to calculate the intersection points of columns 0 to 100, i.e. The line connecting the intersection points forms the roll angle deviation line.

[0092] Figure 2 This diagram illustrates the calculation of the elevation and roll deviations of the satellite's actual beam pointing relative to the beam pointing from the calibration station, using measurement data from the calibration station. The diagram shows the calculation process for the elevation and roll deviations within the range [-0.1, 0.1]. The elevation deviation line is formed by the intersection of the calibration station's east-west beam normalized values ​​with the east-west beam power normalized values ​​in each row of matrix D. The roll deviation line is formed by the intersection of the calibration station's north-south beam normalized values ​​with the north-south beam power normalized values ​​in each column of matrix D. The intersection of the elevation and roll deviation lines represents the elevation and roll deviations between the pointing angle from the satellite's spatial point S to the calibration station's position point B and the pointing angle from point S to the actual beam pointing center point C.

[0093] Step (7): Based on the pitch angle deviation line and roll angle deviation line from steps (5) and (6), calculate the intersection of the two lines. The intersection point represents the pitch and roll deviation angles between the pointing angle from satellite space point S to calibration station position point B and the pointing angle from satellite space point S to the actual beam pointing center point C. The calculation steps are as follows:

[0094] Step (7-1): Obtain the pitch deviation point from step (5) and the roll deviation point from step (6);

[0095] Step (7-2): Obtain the first pitch deviation point coordinate and the pitch deviation point coordinate ;

[0096] Step (7-3): Obtain the first One rolling deviation point coordinate and the One rolling deviation point coordinate ;

[0097] Step (7-4): Determine the... pitch deviation point and the pitch deviation point Is the connection with the first One rolling deviation point and the One rolling deviation point The process for determining if lines intersect is as follows:

[0098] Step (7-4-1): First, determine whether the two line segments coincide on the pitch deviation angle axis and the roll deviation angle axis. If the condition is not met, there is no intersection point. Proceed to step (7-5). If they coincide, there may be an intersection point. Proceed to step (7-4-2).

[0099] Step (7-4-2): If both of the following discriminants are true, they intersect; otherwise, they do not intersect. Proceed to step (7-5). The discriminants are as follows:

[0100]

[0101]

[0102] Step (7-5): If there is no intersection point, take the next rolling deviation point and repeat the calculation process from step (7-3) to step (7-5). If the two lines intersect, the result can be calculated. After completion, proceed to step (8), the formula is as follows:

[0103]

[0104]

[0105]

[0106] Step (7-6): If the first The and the first If the line connecting each pitch deviation point does not intersect with the line connecting each roll deviation point, then take the next pitch point and repeat the calculation process from step (7-2) to step (7-6).

[0107] Step (7-7): If the lines connecting each pitch deviation point and the lines connecting each roll deviation point do not intersect, the calculation accuracy is exceeded and the satellite beam pointing deviation cannot be calculated. Then proceed to step (10).

[0108] Step (8): Determine the orbital position based on the satellite's theoretical position. The theoretical center position of the satellite beam on the ground He Biao School Station Location Calculate the pitch and roll deviation angles of the pointing angle from satellite space point S to calibration station position point B relative to the pointing angle from satellite space point S to the theoretical beam pointing center point A. The calculation steps are as follows:

[0109] Step (8-1): Satellite position point S Beam theory points to center point A Location of the marker station: Point B The unit for latitude and longitude is °. The altitude of the GEO satellite, For the Earth's radius, =42164Km, =6371km, and The value range is [0, 60];

[0110] Step (8-2): Calculate the pointing angle from satellite space point S to calibration station location point B. ,in Pi is represented by the following formula:

[0111]

[0112]

[0113] Step (8-3): Calculate the pointing angle from satellite spatial point S to the theoretical beam pointing center point A. The formula is as follows:

[0114]

[0115]

[0116] Step (8-4): Calculate based on the results of steps (8-2) and (8-3). The formula is as follows:

[0117]

[0118]

[0119] Step (9): Based on the calculation results of Steps (7) and (8), calculate the pitch deviation angle and roll deviation angle of the pointing angle from the satellite space point S to the actual beam pointing center point C relative to the pointing angle from the satellite space point S to the theoretical beam pointing center point A. The calculation result is the satellite beam pointing deviation. The calculation steps are as follows:

[0120] Step (9-1): Calculate the beam pointing pitch deviation angle, using the following formula:

[0121]

[0122] Step (9-2): Calculate the beam pointing roll deviation angle, using the following formula:

[0123]

[0124] Step (9-3): Based on steps (9-2) and (9-3), obtain the satellite beam pointing deviation angle. Proceed to step (10).

[0125] Step (10): Output the calculation results to complete the calculation process.

[0126] This invention targets GEO fixed-beam communication satellites, utilizing measurement signals from mobile calibration stations to calculate beam pointing deviation, which can be used to assist in satellite beam pointing correction. This invention supports mobile platform deployment and can be used to realize satellite beam pointing deviation calculation based on measurement analysis from mobile calibration stations, reducing the difficulty and cost of site construction compared to fixed-site deployment.

[0127] The present invention has been described in detail above with reference to the accompanying drawings. However, those skilled in the art should understand that the specification is for interpreting the claims, and the scope of protection of the present invention shall be determined by the claims. Any modifications made based on the present invention shall be within the scope of protection claimed.

Claims

1. A method for calculating GEO satellite beam pointing deviation based on mobile measurement analysis, applied to beam calibration stations, characterized in that... Includes the following steps: Step 1: Obtain the longitude of the satellite's theoretical orbit determination point S. The longitude of the theoretical center point A of the satellite beam on the ground. and latitude Load the beam pointing deviation power design dataset for four cardinal directions (north, south, east, and west). The data in the dataset are... ,in For pitch angle deviation, For the roll angle deviation, For the beam power value of the East Standard, To calibrate the beam power value for the West standard, To calibrate the beam power value for the South Standard, The beam power value is for the North Standard calibration. Step 2, based on the data Calculate the pointing deviation data after power normalization. ,in, This is the normalized value for the east-west calibration beam power. The normalized value of the north-south calibration beam power is calculated using the following formula: Step 3, generate the pointing deviation judgment matrix D: ， in, express , , That is, the 0th row contains 101 data points with a roll angle deviation of -0.5°, and the roll angle of each subsequent row increases by 0.01°; the 0th column contains 101 data points with a pitch angle deviation of -0.5°, and the pitch angle of each subsequent column increases by 0.01°. Step 4: The beam calibration station receives calibration signals from the four calibration beams of the satellite and generates measurement data. ,in, The east, west, south, and north beam power values ​​measured at the beam calibration station. To measure time, To determine the longitude of the station, To determine the latitude of the calibration station; calculate the normalized power values ​​of the east and west calibration beams. and the normalized value of the north-south calibration beam power ; Step 5, using Calculate the pitch deviation line using the pointing deviation judgment matrix D; Step 6, using The rolling angle deviation line is calculated using the pointing deviation judgment matrix D; Step 7: Calculate the intersection of the pitch deviation line and the roll deviation line. , and That is, the pitch deviation angle and roll deviation angle between the pointing angle from satellite space point S to calibration station position point B and the pointing angle from satellite space point S to the actual beam pointing center point C; Step 8: Based on the position of the satellite's theoretical orbit determination point S The theoretical location of the satellite beam pointing to its center point A on the ground. He Biao School Station Location Calculate the elevation deviation angle of the pointing angle from satellite space point S to calibration station position point B relative to the pointing angle from satellite space point S to the theoretical beam pointing center point A. and rolling deviation angle ; Step 9: Based on the calculation results of Steps 7 and 8, calculate the elevation deviation angle of the pointing angle from the satellite space point S to the actual beam pointing center point C relative to the pointing angle from the satellite space point S to the theoretical beam pointing center point A. and rolling deviation angle The calculation result is the satellite beam pointing deviation.

2. The method for calculating GEO satellite beam pointing deviation based on maneuver measurement analysis according to claim 1, characterized in that, The specific method for step 5 is as follows: Step 5-1: Select the deviation range [-0.5°, 0.5°], where the power of the west calibration beam monotonically decreases from west to east, and the power of the east calibration beam monotonically increases. Pointing deviation judgment matrix. The beam normalization values ​​of each row are monotonically increasing; Step 5-2, determine the matrix based on the pointing deviation. Take out the first one If we disregard the normalized power values ​​of the north-south calibrated beams in the row of data, then any data point in that row is... ,in For column number, from arrive The values ​​are monotonically increasing; search for the first value starting from the first column. Value greater than The column, column number is ,if Or there is no greater than If a column is not found, it indicates that the row has no intersection points; otherwise, the column number is calculated according to the following formula. Intersection of line data and pitch deviation line : Step 5-3: Repeat step 5-2 to calculate the intersection points from row 0 to row 100, i.e., the pitch deviation points. The line connecting the intersection points forms the pitch angle deviation line.

3. The method for calculating GEO satellite beam pointing deviation based on maneuver measurement analysis according to claim 2, characterized in that, The specific method for step 6 is as follows: Step 6-1: Select the deviation range [-0.5°, 0.5°], where the power of the north calibration beam decreases monotonically from north to south, and the power of the south calibration beam increases monotonically. (Pointing deviation judgment matrix) The normalized values ​​of the north and south beams in each column decrease monotonically. Step 6-2, determine the matrix based on the pointing deviation. Take out the first one If we disregard the normalized power values ​​of the east-west calibrated beams in this column, then any data point in this column is... ,in For line numbers, from arrive The value is monotonically decreasing; search for the first value starting from the first row. Value less than The line number is ,if Or there is no less than If a row has no intersection, it means that the row has no intersection points; otherwise, calculate the nth row according to the following formula. Intersection of column data and rolling deviation line : Step 6-3: Repeat step 6-2 to calculate the intersection points from column 0 to column 100, i.e., the rolling deviation points. The line connecting the intersection points forms the roll angle deviation line.

4. The method for calculating GEO satellite beam pointing deviation based on maneuver measurement analysis according to claim 3, characterized in that, The specific method for step 7 is as follows: Step 7-1: Obtain the pitch deviation point from Step 5 and the roll deviation point from Step 6; Step 7-2, obtain the first pitch deviation point coordinates and the pitch deviation point coordinates ; Step 7-3, obtain the first One rolling deviation point coordinates and the One rolling deviation point coordinates ; Step 7-4, determine the first pitch deviation point and the pitch deviation point Does the connection match the first? One rolling deviation point and the One rolling deviation point The lines intersect, and the following method is used to determine if they do: Step 7-4-1: Determine whether the two line segments coincide on both the pitch deviation angle axis and the roll deviation angle axis. If not, there can be no intersection point. Proceed to step 7-5. If yes, proceed to step 7-4-2. Step 7-4-2: If both of the following discriminants are true, then the pair intersects; otherwise, they do not intersect. The discriminants are as follows: Step 7-5: If no intersection point exists, take the next rolling deviation point and repeat steps 7-3 to 7-5. If the two lines intersect, the calculation is complete. Proceed to step 8; the formula is as follows: in, As an intermediate variable, the calculation method is as follows: ； Steps 7-6, if the first The and the first If the line connecting each pitch deviation point does not intersect with the line connecting each roll deviation point, then take the next pitch point and repeat steps 7-2 to 7-6. Step 7-7: If the lines connecting each pitch deviation point and the lines connecting each roll deviation point do not intersect, the calculation accuracy is exceeded, and the satellite beam pointing deviation cannot be calculated.

5. The method for calculating GEO satellite beam pointing deviation based on maneuver measurement analysis according to claim 4, characterized in that, The specific method for step 8 is as follows: Step 8-1, satellite spatial point S is The beam theory points to the center point A as The location of the marker station is point B. The unit for latitude and longitude is °. The altitude of the GEO satellite, For the Earth's radius, and The value range is [0, 60]; Step 8-2: Calculate the pointing angle from satellite space point S to calibration station location point B. The formula is as follows: Step 8-3: Calculate the pointing angle from satellite spatial point S to the theoretical beam pointing center point A. The formula is as follows: Step 8-4: Based on the results of steps 8-2 and 8-3, calculate... The formula is as follows: 。 6. The method for calculating GEO satellite beam pointing deviation based on maneuver measurement analysis according to claim 5, characterized in that, The specific method for step 9 is as follows: Step 9-1, calculate the beam pointing elevation deviation angle: Step 9-2, calculate the beam pointing roll deviation angle: Step 9-3: Obtain the satellite beam pointing deviation angle. .

Images