A method and device for generating atmospheric environment parameters of a full orbit of an aircraft

Abstract

The present invention discloses a method for generating atmospheric environment parameters of a full orbit of an aircraft, comprising: obtaining an aircraft orbit information set and position coordinates of an aircraft orbit sampling point; the aircraft orbit information set comprises aircraft orbit range information, aircraft orbit potential height information and latitude information; sampling and processing the aircraft orbit range information to obtain a discrete height value set and a discrete grid point set; obtaining a discrete altitude layer meteorological element value set; the discrete altitude layer meteorological element value set comprises a meteorological element value of each discrete altitude; based on the aircraft orbit information set, the discrete altitude layer meteorological element value set and the discrete grid point set, predicting and processing the atmospheric environment parameters of the aircraft orbit sampling point to obtain predicted values ​​of the atmospheric environment parameters of the aircraft orbit sampling point.

Description

Technical Field

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

[0002] When an aircraft flies in the air, it is affected by various atmospheric environment factors, among which meteorological elements such as air wind, density, and temperature have a great impact on its flight accuracy. Due to the long flight distance and high flight altitude of some aircraft, the cumulative effect of various meteorological elements on the orbit, especially the high-altitude wind, will 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, air pressure, density, and temperature on the entire orbital plane. For long-range aircraft, the atmospheric environment changes in the entire orbit have a great impact on the flight trajectory accuracy, which will cause unacceptable position deviations of the aircraft. Therefore, it is necessary to provide atmospheric environment data for the flight orbit of the aircraft to support the assessment and guarantee of the impact of atmospheric environment changes in the entire flight orbit. Summary of the Invention

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

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

[0006] S1, obtaining a set of aircraft orbit information and the position coordinates of the sampling points on the aircraft orbit; the set of aircraft orbit information includes aircraft orbit range information, aircraft orbit geopotential height information, and latitude information; the aircraft orbit range information includes altitude range information and horizontal range information;

[0007] S2, performing sampling processing on the aircraft orbit range information to obtain a set of discrete altitude values and a set of discrete grid points;

[0008] S3, obtaining a set of meteorological element values for discrete altitude layers; the set of meteorological element values for discrete altitude layers includes the meteorological element values for each discrete altitude;

[0009] S4, based on the set of aircraft orbit information, the set of meteorological element values for discrete altitude layers, and the set of discrete grid points, performing prediction processing on the atmospheric environment parameters of the sampling points on the aircraft orbit to obtain the predicted values of the atmospheric environment parameters of the sampling points on the aircraft orbit.

[0010] Performing sampling processing on the aircraft orbit range information to obtain a discrete grid point set, including:

[0011] S21, performing uniform sampling processing on the altitude range information to obtain a discrete altitude value set; the discrete altitude value set includes a number of discrete altitude 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 number of discrete plane coordinate points;

[0013] S23, performing combination processing on the discrete altitude value set and the discrete plane coordinate set to obtain a discrete grid point set; the discrete grid point set includes a number of discrete grid points; the position coordinates of the discrete grid points include discrete altitude values and discrete plane coordinate points.

[0014] Based on the aircraft orbit information set, the discrete altitude layer meteorological element value set, and the discrete grid point set, performing prediction processing on the atmospheric environment parameters of the aircraft orbit sampling points to obtain the predicted values of the atmospheric environment parameters of the aircraft orbit sampling points, including:

[0015] S41, using a geometric altitude calculation model to perform calculation processing on the geopotential altitude information and latitude information of the aircraft orbit to obtain a geometric altitude information set; the geometric altitude information set includes a number of geometric altitude values;

[0016] S42, for each discrete grid point in the discrete grid point set, determining the geometric altitude value closest to the discrete altitude value of the discrete grid point, and using the geometric altitude value to perform update processing on the discrete altitude value of the discrete grid point;

[0017] S43, using a standard altitude layer meteorological element calculation model to perform calculation processing on 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 values of each discrete grid point; the meteorological element values include wind field, air pressure, temperature, and air density;

[0018] S44, in the discrete grid point set, finding the nine discrete grid points closest to the aircraft orbit sampling point;

[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 the predicted values of the atmospheric environment parameters of the aircraft orbit sampling point.

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

[0021]

[0022] where: h is the calculated geometric height value, with the unit of m; H is the geopotential height of the aircraft orbit; represents the latitude information, represents the latitude of when the acceleration due to gravity, with the unit of m / s 2 ; represents the calibrated value of the Earth's radius at the latitude of ; g np represents the dimensionless standard acceleration due to gravity, with a value of 9.80665 m / s 2 .

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

[0024]

[0025] where X(h0, x0, y0) is the predicted wind field value of the discrete grid point with the geometric height value of h0 and the discrete plane coordinate points of (x0, y0), h 1 , h 2 , h 3 are the three discrete height values closest to h0 in the discrete height value set, X 1 , X 2 , X 3 are the wind field values corresponding to h 1 , h 2 , h 3 at the three heights respectively in the discrete altitude 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 of h0 and the discrete plane coordinate points of (x0, y0), Y 1 , Y 2 are the air pressure, temperature or air density corresponding to h 1 , h 2 at the two heights respectively in the discrete altitude 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, X 1 , X 2 , X 3 , Y 1 , Y 2 are all obtained from the discrete altitude layer meteorological element value set.

[0026] 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 the predicted values of the atmospheric environment parameters of the aircraft orbit sampling point, including:

[0027] S451. Perform boundary discrimination processing on the coordinate values of the aircraft orbit sampling points and the nine discrete grid points closest to the aircraft orbit sampling points respectively to obtain boundary discrimination results;

[0028] S452. If the boundary discrimination result is yes, determine the meteorological element value of the discrete grid point with the closest planar distance to the aircraft orbit sampling point among the nine discrete grid points closest to the aircraft orbit sampling point as the predicted value of the atmospheric environment parameters of the aircraft orbit sampling point;

[0029] If the discrimination result is no, use the fusion prediction model to perform calculation processing on 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 parameters of the aircraft orbit sampling point.

[0030] The boundary discrimination processing includes:

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

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

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

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

[0035] where 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 aircraft orbit sampling point, (x i , y i ) is the discrete plane coordinate of the i-th discrete grid point among the nine discrete grid points; ∣(x n - x i )(y n - y i )∣ is the planar distance between the i-th discrete grid point and the aircraft orbit sampling point;

[0036] If none of the nine discrete grid points satisfy the inequality group, confirm that the boundary discrimination result is no; if there are points among the nine discrete grid points that satisfy the inequality group, confirm that the boundary discrimination result is yes.

[0037] The fusion prediction model includes:

[0038]

[0039] Among them, f n is the predicted value of the atmospheric environment parameter at the flight vehicle orbit sampling point, and the predicted value of the atmospheric environment parameter is the predicted value of one parameter among wind field, air pressure, temperature or air density; A p,q are the weights with serial numbers p and q, f p,q is the meteorological element value at 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, and the value range is -1 to 1, which are used to represent the serial numbers of the nine discrete grid points closest to the flight vehicle orbit sampling point, (x n , y n ) is the discrete plane coordinate of the flight vehicle 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 corresponding to the corresponding serial numbers among the nine discrete grid points.

[0040] In the second aspect of the embodiments of the present invention, an apparatus for generating atmospheric environment parameters of the entire orbit of a flight vehicle is disclosed, and the apparatus includes:

[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 and executes the method for generating atmospheric environment parameters of the entire orbit of the flight vehicle.

[0044] In the third aspect of the embodiments of the present invention, a computer-readable storage medium is disclosed, and the computer-readable storage medium stores computer instructions, and when the computer instructions are called by the computer, they are used to execute the method for generating atmospheric environment parameters of the entire orbit of the flight vehicle.

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

[0046] The present invention can be used for solving the orbit meteorological data products of new long-distance flight vehicles such as new long-range ones, and can effectively reduce the position deviation of the flight vehicle. After measurement, after using this method for meteorological forecasting, the deviation reduction can reach 25%-30%.

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

[0048] The present invention can be used for calculating atmospheric parameters of the flight orbits of various aircraft and other aircraft, and provides data support for the meteorological service guarantee of the entire orbit. Description of the Drawings

[0049] Figure 1 It is a flowchart of the implementation of the method of the present invention. Detailed Embodiments

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

[0051] Figure 1 It is a flowchart of the implementation of the method of the present invention.

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

[0053] S1, obtaining a set of aircraft orbit information; the set of aircraft orbit information includes aircraft orbit range information, aircraft orbit geopotential height information, and latitude information; the aircraft orbit range information includes altitude range information and horizontal range information;

[0054] S2, performing sampling processing on the aircraft orbit range information to obtain a set of discrete altitude values and a set of discrete grid points;

[0055] S3, obtaining a set of meteorological element values at discrete altitude levels; the set of meteorological element values at discrete altitude levels includes meteorological element values at each discrete altitude;

[0056] S4, based on the set of aircraft orbit information, the set of meteorological element values at discrete altitude levels, and the set of discrete grid points, performing prediction processing on the atmospheric environment parameters of the aircraft orbit sampling points to obtain predicted values of the atmospheric environment parameters of the aircraft orbit sampling points.

[0057] The performing sampling processing on the aircraft orbit range information to obtain a set of discrete grid points includes:

[0058] S21, performing uniform sampling processing on the altitude range information to obtain a set of discrete altitude values; the set of discrete altitude values includes several discrete altitude values;

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

[0060] S23. Combine and process 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 number of discrete grid points; the discrete grid points include discrete height values and discrete plane coordinate points;

[0061] Predict and process 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, including:

[0062] S41. Use the geometric height calculation model to calculate and process the geopotential height information and latitude information of the aircraft orbit to obtain a geometric height information set; the geometric height information set includes a number of geometric height values;

[0063] S42. For each discrete grid point in the discrete grid point set, determine the geometric height value closest to the discrete height value of the discrete grid point, and use the geometric height value to update the discrete height value of the discrete grid point;

[0064] S43. Use the standard altitude layer meteorological element calculation model to calculate and process 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 values of each discrete grid point; the meteorological element values include wind field, air pressure, temperature, and air density;

[0065] S44. In the discrete grid point set, find the nine discrete grid points closest to the aircraft orbit sampling point;

[0066] S45. Perform prediction calculation processing on the nine discrete grid points closest to the aircraft orbit sampling point and the aircraft orbit sampling point to obtain the predicted values of the atmospheric environment parameters 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, unit: m; H is the geopotential height of the aircraft orbit, unit: gpm; represents the latitude information, represents the latitude of when the gravitational acceleration, unit: m / s 2 ; represents the latitude of Calibration value of the Earth's radius at that time; g np Represents the dimensionless standard acceleration of gravity, with a value of 9.80665 m / s 2 .

[0071] The expression of the calculation model of the meteorological elements at the standard altitude layer is as follows:

[0072]

[0073] Among them, X(h0, x0, y0) is the predicted value of the wind field at the discrete grid point with the geometric height value of h0 and the discrete plane coordinate point (x0, y0), h 1 , h 2 , h 3 are the three discrete height values closest to h0 in the set of discrete height values, X 1 , X 2 , X 3 are the wind field values corresponding to h 1 , h 2 , h 3 respectively at the three heights in the set of meteorological element values at the discrete height layer; Y(h0, x0, y0) is the predicted value of the air pressure, temperature or air density at the discrete grid point with the geometric height value of h0 and the discrete plane coordinate point (x0, y0), Y 1 , Y 2 are the air pressure, temperature or air density corresponding to h 1 , h 2 respectively at the two heights in the set of meteorological element values at the discrete height layer; the predicted values of the wind field, air pressure, temperature and air density at the discrete grid point are the meteorological element values of the discrete grid point in the set of meteorological elements at the discrete grid point. X 1 , X 2 , X 3 , Y 1 , Y 2 are all obtained from the set of meteorological element values at the discrete height layer.

[0074] Performing predictive calculation processing on the nine discrete grid points closest to the aircraft orbit sampling point and the aircraft orbit sampling point to obtain the predicted values of the atmospheric environment parameters at the aircraft orbit sampling point, including:

[0075] S451. Performing boundary discrimination processing on the coordinate value of the aircraft orbit sampling point and the nine discrete grid points closest to the aircraft orbit sampling point respectively to obtain the boundary discrimination results;

[0076] S452. If the boundary discrimination result is yes, determine the meteorological element value of the discrete grid point closest to the aircraft orbit sampling point in terms of the plane distance among 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 discrimination result is no, use the fusion prediction model 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 discrimination process includes:

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

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

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

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

[0083] where 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 aircraft orbit sampling point, (x i ,y i ) is the discrete plane coordinate of the i-th discrete grid point among the nine discrete grid points; ∣(x n -x i )(y n -y i )∣ is the plane distance between the i-th discrete grid point and the aircraft orbit sampling point;

[0084] If none of the nine discrete grid points satisfy the inequality group, confirm that the boundary discrimination result is no; if there are points among the nine discrete grid points that satisfy the inequality group, confirm that the boundary discrimination result is yes.

[0085] The fusion prediction model includes:

[0086]

[0087]

[0088] where, fn is the predicted value of the atmospheric environment parameters at the aircraft orbit sampling point. The predicted value of the atmospheric environment parameters is the predicted value of one of the parameters of wind field, air pressure, temperature or air density; A p,q are the weights with serial numbers p and q, f p,q are the meteorological element values at the discrete grid points with serial numbers p and q. The meteorological element values are one of wind field, air pressure, temperature or air density; p, p′, q, q′ are integers, and the value range is -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 aircraft 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 corresponding to the serial numbers among the nine discrete grid points.

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

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

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

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

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

[0094] Step 3: Read in the sampling points of the aircraft orbit. For the sampling point (x n , y n , h n ), when the sampling point coincides with the data coordinate point, directly extract the element data of this coordinate point; when the sampling point and the data coordinate point (xn , y n ) coincide in position. Using the meteorological elements of the standard geometric height layer, directly calculate the height h of the sampling point according to the description in the second step n of the atmospheric environment element value; otherwise, select the nine grid points closest to this point (x i , y j , h n ). According to the description in the second step, calculate the atmospheric environment element values of the nine points respectively, and use the Lagrange interpolation method to calculate the meteorological element values of the sampling point according to the position weights of the nine-point data.

[0095] The method for generating the atmospheric environment parameters of the full orbit of the aircraft described above further includes:

[0096] At several moments, obtain the predicted values of the atmospheric environment parameters of the aircraft orbit sampling points according to step S4 respectively. Use the meteorological element values of the nine closest discrete grid points obtained at each moment to construct the time series values of the meteorological element values; perform cross-correlation processing on all the time series values of the meteorological element values to obtain the cross-correlation matrix; perform eigenvalue decomposition processing on the cross-correlation matrix to obtain the eigenvectors; use the eigenvectors to perform weighted summation on the predicted values of the atmospheric environment parameters of the aircraft orbit sampling points at each moment to obtain the final predicted values of the atmospheric environment parameters of the aircraft orbit sampling points.

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

[0098] C = VDV H ,

[0099] where 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 the second aspect of the embodiments of the present invention, an apparatus for generating atmospheric environment parameters of the full orbit of an aircraft is disclosed. The apparatus includes:

[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 and executes the method for generating the atmospheric environment parameters of the full orbit of the aircraft described above.

[0104] In the third aspect of the embodiments of the present invention, a computer-readable storage medium is disclosed. The computer-readable storage medium stores computer instructions, and when the computer instructions are called by the computer, they are used to execute the method for generating the atmospheric environment parameters of the full orbit of the aircraft described above.

[0105] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, various changes and modifications can be made to the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included within 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, predicting the atmospheric environment parameters of the aircraft orbit sampling point to obtain the predicted value of the atmospheric environment parameters of the aircraft orbit sampling point; 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.

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 1, 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 .

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 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.

5. The method for generating atmospheric environment parameters of a full trajectory of an aircraft as claimed in claim 3, 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.

6. The method for generating atmospheric environment parameters of a full trajectory of an aircraft as claimed in claim 5, 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.

7. The method for generating atmospheric environment parameters of a full trajectory of an aircraft as claimed in claim 5, 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, and 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.

8. 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 7.

9. 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 7.

Images