Method and device for generating atmospheric environment parameters of full orbit of aircraft

Abstract

The invention discloses a method for generating atmospheric environment parameters of an aircraft full orbit. The method comprises the following steps: acquiring an aircraft orbit information set and position coordinates of aircraft orbit sampling points; the aircraft orbit information set comprises aircraft orbit range information, aircraft orbit potential height information and latitude information; carrying out sampling processing on the aircraft orbit range information to obtain a discrete height value set and a discrete grid point set; obtaining a discrete height layer meteorological element value set; the discrete height layer meteorological element value set comprises a meteorological element value of each discrete height; and on the basis of the aircraft orbit information set, the discrete height layer meteorological element value set and the discrete grid point set, performing prediction processing on atmospheric environment parameters of the aircraft orbit sampling points to obtain atmospheric environment parameter prediction values of the aircraft orbit sampling points.

Description

Technical Field

[0001] The present invention relates to the fields of atmospheric environment and numerical analysis, and in particular to a method and device for generating atmospheric environment parameters of a full orbit of an aircraft. Background Art

[0002] When an aircraft is flying in the air, it is affected by various atmospheric environmental factors, among which the wind, density, temperature and other meteorological factors have a great impact on its flight accuracy. Since some aircraft fly long distances and high altitudes, the cumulative effect of various meteorological factors on the track, especially the high-altitude wind, can cause a large deviation in the aircraft trajectory.

[0003] In previous aircraft simulation tests, only the air wind observation data or atmospheric standard values ​​at the take-off point were used, lacking the support of environmental data such as air wind and air pressure, density, and temperature on the entire orbit. For long-range aircraft, the changes in the atmospheric environment of the entire orbit have a great impact on the accuracy of the flight trajectory, which will cause unacceptable position deviations of the aircraft. Therefore, it is necessary to provide atmospheric environment data for the flight trajectory of the aircraft to provide support for the assessment and guarantee of the impact of atmospheric environment changes on the entire flight trajectory. Summary of the invention

[0004] The present invention mainly solves the problem that the atmospheric environment changes of a long-range aircraft throughout its entire orbit have a great impact on the accuracy of its flight trajectory, which may cause unacceptable position deviations of the aircraft. The present invention discloses a method and device for generating atmospheric environment parameters of an aircraft throughout its entire orbit.

[0005] In a first aspect of an embodiment of the present invention, a method for generating atmospheric environment parameters of a full trajectory of an aircraft is disclosed, comprising:

[0006] S1, obtaining a set of spacecraft orbit information and position coordinates of a spacecraft orbit sampling point; the set of spacecraft orbit information includes spacecraft orbit range information, spacecraft orbit geopotential height information and latitude information; the spacecraft orbit range information includes height range information and horizontal range information;

[0007] S2, sampling and processing the spacecraft orbit range information to obtain a discrete height value set and a discrete grid point set;

[0008] S3, obtaining a discrete altitude layer meteorological element value set; the discrete altitude layer meteorological element value set includes a meteorological element value for each discrete altitude;

[0009] S4, based on the aircraft orbit information set, the discrete altitude layer meteorological element value set and the discrete grid point set, the atmospheric environment parameters of the aircraft orbit sampling points are predicted and processed to obtain the predicted values ​​of the atmospheric environment parameters of the aircraft orbit sampling points.

[0010] The sampling process of the spacecraft orbit range information to obtain a discrete grid point set includes:

[0011] S21, uniformly sampling the height range information to obtain a discrete height value set; the discrete height value set includes a plurality of discrete height values;

[0012] S22, performing two-dimensional uniform sampling processing on the horizontal range information to obtain a discrete plane coordinate set; the discrete plane coordinate set includes a plurality of discrete plane coordinate points;

[0013] S23, combining the discrete height value set and the discrete plane coordinate set to obtain a discrete grid point set; the discrete grid point set includes a plurality of discrete grid points; the position coordinates of the discrete grid points include discrete height values ​​and discrete plane coordinate points.

[0014] The method of predicting the atmospheric environment parameters of the aircraft orbit sampling points based on the aircraft orbit information set, the discrete altitude layer meteorological element value set and the discrete grid point set to obtain the predicted values ​​of the atmospheric environment parameters of the aircraft orbit sampling points includes:

[0015] S41, using a geometric height calculation model, calculating and processing the orbital potential height information and latitude information of the aircraft to obtain a geometric height information set; the geometric height information set includes a plurality of geometric height values;

[0016] S42, for each discrete grid point in the discrete grid point set, determining a geometric height value closest to a discrete height value of the discrete grid point, and updating the discrete height value of the discrete grid point using the geometric height value;

[0017] S43, using the standard altitude layer meteorological element calculation model, calculating and processing the discrete altitude layer meteorological element value set to obtain a discrete grid point meteorological element set; the discrete grid point meteorological element set includes the meteorological element value of each discrete grid point; the meteorological element value includes wind field, air pressure, temperature, and air density;

[0018] S44, searching for nine discrete grid points closest to the spacecraft orbit sampling point in the discrete grid point set;

[0019] S45, performing prediction calculation processing on the nine discrete grid points closest to the aircraft orbit sampling point and the aircraft orbit sampling point to obtain a predicted value of the atmospheric environment parameter of the aircraft orbit sampling point.

[0020] The expression of the geometric height calculation model is:

[0021]

[0022] Where: h is the calculated geometric height value, in meters; H is the orbital potential height of the spacecraft; Indicates latitude information. Indicates latitude The gravitational acceleration at time , in m / s 2 ; Indicates latitude The calibration value of the earth's radius at ; g np Represents the dimensionless standard gravitational acceleration, with a value of 9.80665m / s 2 .

[0023] The expression of the standard altitude layer meteorological element calculation model is:

[0024]

[0025] Among them, X(h0, x0, y0) is the predicted value of the wind field of the discrete grid point with the geometric height value h0 and the discrete plane coordinate point (x0, y0), h1, h2, h3 are the three discrete height values ​​closest to h0 in the discrete height value set, X1, X2, X3 are the wind field values ​​corresponding to the three heights h1, h2, h3 in the discrete height layer meteorological element value set; Y(h0, x0, y0) is the predicted value of the air pressure, temperature or air density of the discrete grid point with the geometric height value h0 and the discrete plane coordinate point (x0, y0), Y1, Y2 are the air pressure, temperature or air density corresponding to the two heights h1 and h2 in the discrete height layer meteorological element value set; the predicted values ​​of the wind field, air pressure, temperature and air density of the discrete grid point are the meteorological element values ​​of the discrete grid point in the discrete grid point meteorological element set. Among them, X1, X2, X3, Y1, Y2 are all obtained from the discrete height layer meteorological element value set.

[0026] The predictive calculation processing of the nine discrete grid points closest to the aircraft orbit sampling point and the aircraft orbit sampling point to obtain the predicted value of the atmospheric environment parameter of the aircraft orbit sampling point includes:

[0027] S451, performing boundary discrimination processing on the coordinate value of the aircraft trajectory sampling point and the nine discrete grid points closest to the aircraft trajectory sampling point to obtain a boundary discrimination result;

[0028] S452, if the boundary determination result is yes, determine the meteorological element value of the discrete grid point closest to the aircraft orbit sampling point in the nine discrete grid points closest to the aircraft orbit sampling point, as the predicted value of the atmospheric environment parameter of the aircraft orbit sampling point;

[0029] If the judgment result is no, the fusion prediction model is used to calculate and process the meteorological element values ​​of the nine discrete grid points closest to the aircraft orbit sampling point to obtain the predicted value of the atmospheric environment parameter of the aircraft orbit sampling point.

[0030] The boundary determination process includes:

[0031] For each of the nine discrete grid points, determine whether the inequality group is satisfied:

[0032] ∣x n -x i ∣≤a,

[0033] ∣y n -y i ∣≤a,

[0034] ∣(x n -x i )(y n -y i )∣≤a2,

[0035] The preset distance error threshold is a, the preset area constraint threshold is a2, (x n ,y n ) is the discrete plane coordinate of the spacecraft orbit sampling point, (x i ,y i ) is the discrete plane coordinate of the ith discrete grid point among the nine discrete grid points; |(x n -x i )(y n -y i )| is the plane distance between the ith discrete grid point and the spacecraft orbit sampling point;

[0036] If none of the nine discrete grid points satisfy the inequality group, the boundary determination result is confirmed to be no; if there is a point among the nine discrete grid points that satisfies the inequality group, the boundary determination result is confirmed to be yes.

[0037] The fusion prediction model comprises:

[0038]

[0039] Among them, f n A is the predicted value of the atmospheric environment parameter at the sampling point of the aircraft track, wherein the predicted value of the atmospheric environment parameter is the predicted value of one parameter among wind field, air pressure, temperature or air density; p,q is the weight of the numbers p and q, f p,qis the meteorological element value of the discrete grid points with serial numbers p and q, and the meteorological element value is one of wind field, air pressure, temperature or air density; p, p', q, q' are integers with a value range of -1 to 1, which are used to represent the serial numbers of the nine discrete grid points closest to the aircraft orbit sampling point, (x n ,y n ) is the discrete plane coordinate of the spacecraft orbit sampling point, (x p ,y p )、(x p′ ,y p′ )、(x q ,y q )、(x q′ ,y q′ ) are the discrete plane coordinates of the discrete grid points with corresponding serial numbers among the nine discrete grid points.

[0040] In a second aspect of an embodiment of the present invention, a device for generating atmospheric environment parameters of a full trajectory of an aircraft is disclosed, the device comprising:

[0041] A memory storing executable program code;

[0042] a processor coupled to the memory;

[0043] The processor calls the executable program code stored in the memory to execute the method for generating atmospheric environment parameters of the entire trajectory of the aircraft.

[0044] According to a third aspect of an embodiment of the present invention, a computer storable medium is disclosed, wherein the computer storable medium stores computer instructions, and when the computer instructions are called by a computer, they are used to execute the method for generating atmospheric environment parameters of the entire orbit of the aircraft.

[0045] The beneficial effects of the present invention are:

[0046] The present invention can be used for solving orbital meteorological data products of new long-range and long-distance aircraft, and can effectively reduce the position deviation of the aircraft. After measurement and calculation, after using this method for weather forecasting, the deviation can be reduced by up to 25%-30%.

[0047] The present invention is simple and practical, and can directly extract required atmospheric environment parameter data from a large amount of atmospheric environment data quickly, thereby effectively reducing the meteorological forecasting burden of the system.

[0048] The present invention can be used to calculate the atmospheric parameters of the flight orbits of various aircraft such as airplanes, and provide data support for full-orbit meteorological service assurance. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 It is a flow chart for implementing the method of the present invention. DETAILED DESCRIPTION

[0050] In order to better understand the content of the present invention, an embodiment is given here.

[0051] Figure 1 It is a flow chart for implementing the method of the present invention.

[0052] In a first aspect of the embodiment of the present application, a method for generating atmospheric environment parameters of a full trajectory of an aircraft is disclosed, comprising:

[0053] S1, obtaining a set of spacecraft orbit information; the set of spacecraft orbit information includes spacecraft orbit range information, spacecraft orbit potential height information and latitude information; the spacecraft orbit range information includes height range information and horizontal range information;

[0054] S2, sampling and processing the spacecraft orbit range information to obtain a discrete height value set and a discrete grid point set;

[0055] S3, obtaining a discrete altitude layer meteorological element value set; the discrete altitude layer meteorological element value set includes a meteorological element value for each discrete altitude;

[0056] S4, based on the aircraft orbit information set, the discrete altitude layer meteorological element value set and the discrete grid point set, the atmospheric environment parameters of the aircraft orbit sampling points are predicted and processed to obtain the predicted values ​​of the atmospheric environment parameters of the aircraft orbit sampling points.

[0057] The sampling process of the spacecraft orbit range information to obtain a discrete grid point set includes:

[0058] S21, uniformly sampling the height range information to obtain a discrete height value set; the discrete height value set includes a plurality of discrete height values;

[0059] S22, performing two-dimensional uniform sampling processing on the horizontal range information to obtain a discrete plane coordinate set; the discrete plane coordinate set includes a plurality of discrete plane coordinate points;

[0060] S23, combining the discrete height value set and the discrete plane coordinate set to obtain a discrete grid point set; the discrete grid point set includes a plurality of discrete grid points; the discrete grid points include discrete height values ​​and discrete plane coordinate points;

[0061] The method of predicting the atmospheric environment parameters of the aircraft orbit sampling points based on the aircraft orbit information set, the discrete altitude layer meteorological element value set and the discrete grid point set to obtain the predicted values ​​of the atmospheric environment parameters of the aircraft orbit sampling points includes:

[0062] S41, using a geometric height calculation model, calculating and processing the orbital potential height information and latitude information of the aircraft to obtain a geometric height information set; the geometric height information set includes a plurality of geometric height values;

[0063] S42, for each discrete grid point in the discrete grid point set, determining a geometric height value closest to a discrete height value of the discrete grid point, and updating the discrete height value of the discrete grid point using the geometric height value;

[0064] S43, using the standard altitude layer meteorological element calculation model, calculating and processing the discrete altitude layer meteorological element value set to obtain a discrete grid point meteorological element set; the discrete grid point meteorological element set includes the meteorological element value of each discrete grid point; the meteorological element value includes wind field, air pressure, temperature, and air density;

[0065] S44, searching for nine discrete grid points closest to the spacecraft orbit sampling point in the discrete grid point set;

[0066] S45, performing prediction calculation processing on the nine discrete grid points closest to the aircraft orbit sampling point and the aircraft orbit sampling point to obtain a predicted value of the atmospheric environment parameter of the aircraft orbit sampling point;

[0067] The expression of the geometric height calculation model is:

[0068]

[0069]

[0070] Where: h is the calculated geometric height value, in meters; H is the orbital potential height of the spacecraft, in gpm; Indicates latitude information. Indicates latitude The gravitational acceleration at time , in m / s 2 ; Indicates latitude The calibration value of the earth's radius at ; g np Represents the dimensionless standard gravitational acceleration, with a value of 9.80665m / s 2 .

[0071] The expression of the standard altitude layer meteorological element calculation model is:

[0072]

[0073] Among them, X(h0, x0, y0) is the predicted value of the wind field of the discrete grid point with the geometric height value h0 and the discrete plane coordinate point (x0, y0), h1, h2, h3 are the three discrete height values ​​closest to h0 in the discrete height value set, X1, X2, X3 are the wind field values ​​corresponding to the three heights h1, h2, h3 in the discrete height layer meteorological element value set; Y(h0, x0, y0) is the predicted value of the air pressure, temperature or air density of the discrete grid point with the geometric height value h0 and the discrete plane coordinate point (x0, y0), Y1, Y2 are the air pressure, temperature or air density corresponding to the two heights h1 and h2 in the discrete height layer meteorological element value set; the predicted values ​​of the wind field, air pressure, temperature and air density of the discrete grid point are the meteorological element values ​​of the discrete grid point in the discrete grid point meteorological element set. X1, X2, X3, Y1, Y2 are all obtained from the discrete height layer meteorological element value set.

[0074] The predictive calculation processing of the nine discrete grid points closest to the aircraft orbit sampling point and the aircraft orbit sampling point to obtain the predicted value of the atmospheric environment parameter of the aircraft orbit sampling point includes:

[0075] S451, performing boundary discrimination processing on the coordinate value of the aircraft trajectory sampling point and the nine discrete grid points closest to the aircraft trajectory sampling point to obtain a boundary discrimination result;

[0076] S452, if the boundary determination result is yes, determine the meteorological element value of the discrete grid point closest to the aircraft orbit sampling point in the nine discrete grid points closest to the aircraft orbit sampling point, as the predicted value of the atmospheric environment parameter of the aircraft orbit sampling point;

[0077] If the judgment result is no, the fusion prediction model is used to calculate and process the meteorological element values ​​of the nine discrete grid points closest to the aircraft orbit sampling point to obtain the predicted value of the atmospheric environment parameter of the aircraft orbit sampling point.

[0078] The boundary determination process includes:

[0079] For each of the nine discrete grid points, determine whether the inequality group is satisfied:

[0080] ∣x n -x i ∣≤a,

[0081] ∣y n -y i ∣≤a,

[0082] ∣(x n -x i)(y n -y i )∣≤a2,

[0083] The preset distance error threshold is a, the preset area constraint threshold is a2, (x n ,y n ) is the discrete plane coordinate of the spacecraft orbit sampling point, (x i ,y i ) is the discrete plane coordinate of the ith discrete grid point among the nine discrete grid points; |(x n -x i )(y n -y i )| is the plane distance between the ith discrete grid point and the spacecraft orbit sampling point;

[0084] If none of the nine discrete grid points satisfy the inequality group, the boundary determination result is confirmed to be negative; if there is a point among the nine discrete grid points that satisfies the inequality group, the boundary determination result is confirmed to be positive.

[0085] The fusion prediction model comprises:

[0086]

[0087]

[0088] Among them, f n A is the predicted value of the atmospheric environment parameter at the sampling point of the aircraft track, wherein the predicted value of the atmospheric environment parameter is the predicted value of one parameter among wind field, air pressure, temperature or air density; p,q is the weight of the numbers p and q, f p,q is the meteorological element value of the discrete grid points with serial numbers p and q, and the meteorological element value is one of wind field, air pressure, temperature or air density; p, p', q, q' are integers with a value range of -1 to 1, which are used to represent the serial numbers of the nine discrete grid points closest to the aircraft orbit sampling point, (x n ,y n ) is the discrete plane coordinate of the spacecraft orbit sampling point, (x p ,y p )、(x p′ ,y p′ )、(x q ,y q )、(x q′ ,y q′ ) are the discrete plane coordinates of the discrete grid points with corresponding serial numbers among the nine discrete grid points.

[0089] Combined with the aircraft flight experiment process, the steps of the aircraft full track atmospheric environment generation method of the present invention are as follows:

[0090] Step 1: Read the three-dimensional atmospheric environment data of the area at takeoff, including geopotential height, temperature, zonal wind, radial wind, and density. The data coordinates are (longitude, latitude, and pressure). Taking the geopotential heights corresponding to the standard pressure layers in the air at the grid point (121E, 39N) as an example, they are (1000hPa, 97gpm), (925hPa, 742gpm), (850hPa, 1444gpm), (700hPa, 3010gpm), (500hPa, 5590gpm), (400hPa, 7210gpm), (300hPa, 9180gpm), and (250hPa, 10350gpm) ..., convert the potential height into geometric height. The geometric heights corresponding to the standard pressure layers in the air are (1000hPa, 97.3m), (925hPa, 744.1m), (850hPa, 1448.2m), (700hPa, 3019.5m), (500hPa, 5609.9m), (400hPa, 7237.6m), (300hPa, 9218.1m), and (250hPa, 10394.8m).

[0091] Step 2: Determine the standard geometric altitude layer, with an altitude range of 0 to 50km, and take points at intervals of 500m, namely 0, 500, 1000...45000, 50000. According to the standard altitude layer meteorological element calculation model, calculate the environmental values ​​of the zonal wind and radial wind of the standard geometric altitude layer; taking the radial wind of each standard pressure layer in the air at the grid point (121E, 39N) as an example, they are respectively (97.3m, 1.4m / s), (744.1m, 2.8m / s), (1448.2m, 8.9m / s), (3019.5m, 9.0m / s), (5609.9m, 16.9m / s), (7237.6m, 20.9m / s), (9218.1m, 2 2.9m / s), (10394.8m, 27m / s)..., according to the interpolation calculation shown in formula (4), the radial winds at the standard geometric altitude layers are (500m, 1.8m / s), (1000m, 4.5m / s), (1500m, 9.3m / s), (2000m, 11.1m / s), (2500m, 11.0m / s), (3000m, 9.1m / s), (3500m, 9.7m / s), (4000m, 10.8m / s)...

[0092] According to the standard altitude layer meteorological element calculation model, the values ​​of meteorological elements such as air pressure, temperature, and density of the standard geometric altitude layer are calculated. Taking the temperatures of the standard pressure layers in the air at the grid point (121E, 39N) as an example, they are (97.3m, 286.1K), (744.1m, 283.5K), (1448.2m, 280.5K), (3019.5m, 270.4K), (5609.9m, 254.1K), (7237.6m, 242.0K), (9218.1m, 224.2K), (10394.8m, 213 .4K)……, according to the interpolation calculation of the fusion prediction model, the radial winds at the standard geometric altitude layers are (500m, 284.5K), (1000m, 282.4K), (1500m, 280.2K), (2000m, 276.9K), (2500m, 273.7K), (3000m, 279.5K), (3500m, 267.3K), and (4000m, 264.1K).

[0093] At this point, the three-dimensional atmospheric environment data is converted from isobaric surface coordinates to geometric height coordinates.

[0094] Step 3: Read the spacecraft orbit sampling points. For the sampling points (x n ,y n ,h n ), when the sampling point coincides with the data coordinate point, directly extract the element data of the coordinate point; when the sampling point coincides with the data coordinate point (x n ,y n ) position coincides, use the standard geometric altitude layer meteorological elements, directly follow the second step to calculate the sampling point height h n Otherwise, select the nine grid points closest to the point (x i ,y j ,h n ), follow the instructions in the second step to calculate the atmospheric environmental element values ​​of the nine points respectively, and use the Lagrange interpolation method to calculate the values ​​of each meteorological element at the sampling point according to the location weights of the nine-point data.

[0095] The method for generating atmospheric environment parameters of the entire trajectory of the aircraft further includes:

[0096] At several moments, the predicted values ​​of the atmospheric environment parameters of the spacecraft orbit sampling points are obtained respectively according to step S4, and the meteorological element values ​​of the nearest nine discrete grid points obtained at each moment are used to construct the time series values ​​of the meteorological element values; all the time series values ​​of the meteorological element values ​​are mutually processed to obtain a cross-correlation matrix; the cross-correlation matrix is ​​subjected to eigenvalue decomposition processing to obtain eigenvectors; and the predicted values ​​of the atmospheric environment parameters of the spacecraft orbit sampling points at each moment are weighted summed using the eigenvectors to obtain the final predicted values ​​of the atmospheric environment parameters of the spacecraft orbit sampling points.

[0097] Perform eigenvalue decomposition on the cross-correlation matrix C and obtain:

[0098] C=VDV H ,

[0099] Among them, V is the eigenvector matrix, D is the eigenvalue matrix, and the diagonal elements of the matrix D are normalized and used as the weight vector.

[0100] In a second aspect of an embodiment of the present invention, a device for generating atmospheric environment parameters of a full trajectory of an aircraft is disclosed, the device comprising:

[0101] A memory storing executable program code;

[0102] a processor coupled to the memory;

[0103] The processor calls the executable program code stored in the memory to execute the method for generating atmospheric environment parameters of the entire trajectory of the aircraft.

[0104] According to a third aspect of an embodiment of the present invention, a computer storable medium is disclosed, wherein the computer storable medium stores computer instructions, and when the computer instructions are called by a computer, they are used to execute the method for generating atmospheric environment parameters of the entire orbit of the aircraft.

[0105] The above is only an embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.

Claims

1. A method for generating atmospheric environment parameters of a full trajectory of an aircraft, characterized in that: include: S1, obtaining a set of spacecraft orbit information and position coordinates of a spacecraft orbit sampling point; the set of spacecraft orbit information includes spacecraft orbit range information, spacecraft orbit geopotential height information and latitude information; the spacecraft orbit range information includes height range information and horizontal range information; S2, sampling and processing the spacecraft orbit range information to obtain a discrete height value set and a discrete grid point set; S3, obtaining a discrete altitude layer meteorological element value set; the discrete altitude layer meteorological element value set includes a meteorological element value for each discrete altitude; S4, based on the aircraft orbit information set, the discrete altitude layer meteorological element value set and the discrete grid point set, the atmospheric environment parameters of the aircraft orbit sampling points are predicted and processed to obtain the predicted values ​​of the atmospheric environment parameters of the aircraft orbit sampling points.

2. The method for generating atmospheric environment parameters of a full trajectory of an aircraft according to claim 1, characterized in that: The sampling process of the spacecraft orbit range information to obtain a discrete grid point set includes: S21, uniformly sampling the height range information to obtain a discrete height value set; the discrete height value set includes a plurality of discrete height values; S22, performing two-dimensional uniform sampling processing on the horizontal range information to obtain a discrete plane coordinate set; the discrete plane coordinate set includes a plurality of discrete plane coordinate points; S23, combining the discrete height value set and the discrete plane coordinate set to obtain a discrete grid point set; the discrete grid point set includes a plurality of discrete grid points; the position coordinates of the discrete grid points include discrete height values ​​and discrete plane coordinate points.

3. The method for generating atmospheric environment parameters of a full trajectory of an aircraft as claimed in claim 2, characterized in that: The method of predicting the atmospheric environment parameters of the aircraft orbit sampling points based on the aircraft orbit information set, the discrete altitude layer meteorological element value set and the discrete grid point set to obtain the predicted values ​​of the atmospheric environment parameters of the aircraft orbit sampling points includes: S41, using a geometric height calculation model, calculating and processing the orbital potential height information and latitude information of the aircraft to obtain a geometric height information set; the geometric height information set includes a plurality of geometric height values; S42, for each discrete grid point in the discrete grid point set, determining a geometric height value closest to a discrete height value of the discrete grid point, and updating the discrete height value of the discrete grid point using the geometric height value; S43, using the standard altitude layer meteorological element calculation model, calculating and processing the discrete altitude layer meteorological element value set to obtain a discrete grid point meteorological element set; the discrete grid point meteorological element set includes the meteorological element value of each discrete grid point; the meteorological element value includes wind field, air pressure, temperature, and air density; S44, searching for nine discrete grid points closest to the spacecraft orbit sampling point in the discrete grid point set; S45, performing prediction calculation processing on the nine discrete grid points closest to the aircraft orbit sampling point and the aircraft orbit sampling point to obtain a predicted value of the atmospheric environment parameter of the aircraft orbit sampling point.

4. The method for generating atmospheric environment parameters of a full trajectory of an aircraft as claimed in claim 3, characterized in that: The expression of the geometric height calculation model is: Where: h is the calculated geometric height value, in meters; H is the orbital potential height of the spacecraft; Indicates latitude information. Indicates latitude The gravitational acceleration at time , in m / s 2 ; Indicates latitude The calibration value of the earth's radius at ; g np Represents the dimensionless standard gravitational acceleration, with a value of 9.80665m / s 2 .

5. The method for generating atmospheric environment parameters of a full trajectory of an aircraft as claimed in claim 4, characterized in that: The expression of the standard altitude layer meteorological element calculation model is: Among them, X(h0, x0, y0) is the predicted value of the wind field of the discrete grid point with a geometric height value of h0 and a discrete plane coordinate point of (x0, y0), h1, h2, h3 are the three discrete height values ​​closest to h0 in the discrete height value set, X1, X2, X3 are the wind field values ​​corresponding to the three heights h1, h2, h3 in the discrete height layer meteorological element value set; Y(h0, x0, y0) is the predicted value of the air pressure, temperature or air density of the discrete grid point with a geometric height value of h0 and a discrete plane coordinate point of (x0, y0), Y1, Y2 are the air pressure, temperature or air density corresponding to the two heights h1 and h2 in the discrete height layer meteorological element value set; the predicted values ​​of the wind field, air pressure, temperature and air density of the discrete grid point are the meteorological element values ​​of the discrete grid point in the discrete grid point meteorological element set.

6. The method for generating atmospheric environment parameters of a full trajectory of an aircraft as claimed in claim 4, characterized in that: The predictive calculation processing of the nine discrete grid points closest to the aircraft orbit sampling point and the aircraft orbit sampling point to obtain the predicted value of the atmospheric environment parameter of the aircraft orbit sampling point includes: S451, performing boundary discrimination processing on the coordinate value of the aircraft trajectory sampling point and the nine discrete grid points closest to the aircraft trajectory sampling point to obtain a boundary discrimination result; S452, if the boundary determination result is yes, determine the meteorological element value of the discrete grid point closest to the aircraft orbit sampling point in the nine discrete grid points closest to the aircraft orbit sampling point, as the predicted value of the atmospheric environment parameter of the aircraft orbit sampling point; If the judgment result is no, the fusion prediction model is used to calculate and process the meteorological element values ​​of the nine discrete grid points closest to the aircraft orbit sampling point to obtain the predicted value of the atmospheric environment parameter of the aircraft orbit sampling point.

7. The method for generating atmospheric environment parameters of a full trajectory of an aircraft as claimed in claim 6, characterized in that: The boundary determination process includes: For each of the nine discrete grid points, determine whether the inequality group is satisfied: ∣x n -x i ∣≤a, ∣and n -and i ∣≤a, ∣(x n -x i )(and n -and i )∣≤a2, The preset distance error threshold is a, the preset area constraint threshold is a2, (x n ,y n ) is the discrete plane coordinate of the spacecraft orbit sampling point, (x i ,y i ) is the discrete plane coordinate of the ith discrete grid point among the nine discrete grid points; |(x n -x i )(y n -y i )| is the plane distance between the ith discrete grid point and the spacecraft orbit sampling point; If none of the nine discrete grid points satisfy the inequality group, the boundary determination result is confirmed to be no; if there is a point among the nine discrete grid points that satisfies the inequality group, the boundary determination result is confirmed to be yes.

8. The method for generating atmospheric environment parameters of a full trajectory of an aircraft as claimed in claim 6, characterized in that: The fusion prediction model comprises: Among them, f n A is the predicted value of the atmospheric environment parameter at the sampling point of the aircraft track, wherein the predicted value of the atmospheric environment parameter is the predicted value of one parameter among wind field, air pressure, temperature or air density; p,q is the weight of the numbers p and q, f p,q is the meteorological element value of the discrete grid points with serial numbers p and q, and the meteorological element value is one of wind field, air pressure, temperature or air density; p, p', q, q' are integers with a value range of -1 to 1, which are used to represent the serial numbers of the nine discrete grid points closest to the aircraft orbit sampling point, (x n ,y n ) is the discrete plane coordinate of the spacecraft orbit sampling point, (x p ,y p )、(x p′ ,y p′ )、(x q ,y q )、(x q′ ,y q′ ) are the discrete plane coordinates of the discrete grid points with corresponding serial numbers among the nine discrete grid points.

9. A device for generating atmospheric environment parameters of a full trajectory of an aircraft, characterized in that: The device comprises: A memory storing executable program code; a processor coupled to the memory; The processor calls the executable program code stored in the memory to execute the method for generating atmospheric environment parameters of the entire trajectory of an aircraft according to any one of claims 1 to 8.

10. A computer storable medium, characterized in that: The computer storable medium stores computer instructions, and when the computer instructions are called by a computer, they are used to execute the method for generating atmospheric environment parameters of the entire trajectory of an aircraft as described in any one of claims 1 to 8.

Images