An atmospheric environment prediction and evaluation method and device for the entire orbit of an aircraft
By preprocessing the atmospheric environmental parameters of the aircraft orbit and weighted sum of the evaluation model, the accuracy of the atmospheric environmental prediction of the aircraft's full orbit is solved, and the aircraft's flight accuracy is improved.
Patent Information
- Application Number
- CN202411447474.0
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-16
- Publication Date
- 2025-07-25
- Estimated Expiration
- 2044-10-16
AI Technical Summary
The prior art is difficult to accurately and effectively evaluate the atmospheric environment of the aircraft's full orbit, resulting in deviations in the flight trajectory and affecting the flight accuracy.
By obtaining the predicted values and measured values of the atmospheric environmental parameters of the aircraft orbit, performing outlier value removal and normalization, the weighted sum is used using the position difference, wind field and environmental prediction evaluation models to be evaluated, the difference value calculation is calculated, and the optimal solution vector is extracted for evaluation.
It realizes effective evaluation of atmospheric environmental parameters of the entire orbit of the aircraft, and improves prediction accuracy and accuracy, especially the prediction effect of air pressure, temperature and air density.
Smart Images

Figure CN119311783B_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the fields of atmospheric environment and numerical analysis, and particularly relates to a method and device for predicting and evaluating the atmospheric environment of an aircraft's entire orbit. Background Art
[0002] When an aircraft flies in the air, it is affected by various atmospheric environment factors. Among them, 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's trajectory. Therefore, accurately predicting the atmospheric parameters on the aircraft's orbit is of great significance for improving the flight trajectory accuracy of the aircraft. At the same time, during the flight of the aircraft, the flight area involved is extensive, and the meteorological parameters change complexly and diversely. Effectively evaluating the predicted parameters of the atmospheric environment of the entire orbit of the aircraft is of great significance for improving the method for predicting the atmospheric environment of the entire orbit of the aircraft and enhancing the prediction accuracy. Therefore, how to accurately and effectively evaluate the predicted parameters of the atmospheric environment of the entire orbit of the aircraft is an urgent problem to be solved currently. Summary of the Invention
[0003] The present invention mainly solves the problem of how to accurately and effectively evaluate the predicted parameters of the atmospheric environment of the entire orbit of the aircraft, and discloses a method and device for predicting and evaluating the atmospheric environment of the entire orbit of the aircraft.
[0004] In the first aspect of the embodiments of the present application, a method for predicting and evaluating the atmospheric environment of an entire orbit of an aircraft is disclosed, including:
[0005] S1, obtaining a set of predicted values of atmospheric environment parameters and a set of measured values of atmospheric parameters for the aircraft orbit; the set of predicted values of atmospheric environment parameters includes the position coordinate information of several discrete points on the aircraft orbit and the corresponding sequence of predicted values of atmospheric environment parameters; the set of measured values of atmospheric parameters includes the position coordinate information of several discrete points on the aircraft orbit and the corresponding sequence of measured values of atmospheric environment parameters;
[0006] S2, respectively preprocessing the set of predicted values of atmospheric environment parameters and the set of measured values of atmospheric parameters to obtain a preprocessed set of predicted values of atmospheric environment parameters and a preprocessed set of measured values of atmospheric parameters;
[0007] S3, performing a prediction evaluation process on the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters to obtain a prediction evaluation value; the prediction evaluation value is used to evaluate the prediction accuracy of the set of predicted values of atmospheric environment parameters.
[0008] The predicted value sequence of the atmospheric environment parameters includes a predicted value sequence of the wind field, a predicted value sequence of the air pressure, a predicted value sequence of the temperature, and a predicted value sequence of the air density;
[0009] The measured value sequence of the atmospheric environment parameters includes a measured value sequence of the wind field, a measured value sequence of the air pressure, a measured value sequence of the temperature, and a measured value sequence of the air density.
[0010] The preprocessing of the predicted value set of the atmospheric environment parameters and the measured value set of the atmospheric parameters respectively to obtain a preprocessed predicted value set of the atmospheric environment parameters and a preprocessed measured value set of the atmospheric parameters includes:
[0011] S21, performing outlier rejection processing on the predicted value set of the atmospheric environment parameters and the measured value set of the atmospheric parameters respectively to obtain a first predicted value set of the atmospheric environment parameters and a first measured value set of the atmospheric parameters;
[0012] S22, performing normalization processing on the first predicted value set of the atmospheric environment parameters and the first measured value set of the atmospheric parameters respectively to obtain a preprocessed predicted value set of the atmospheric environment parameters and a preprocessed measured value set of the atmospheric parameters.
[0013] The prediction evaluation processing of the preprocessed predicted value set of the atmospheric environment parameters and the preprocessed measured value set of the atmospheric parameters to obtain a prediction evaluation value includes:
[0014] S31, extracting the position coordinate information of all discrete points in the preprocessed predicted value set of the atmospheric environment parameters and the preprocessed measured value set of the atmospheric parameters, and performing position difference evaluation processing on the position coordinate information of all discrete points to obtain a position difference evaluation value;
[0015] S32, performing wind field prediction evaluation processing on the preprocessed predicted value set of the atmospheric environment parameters and the preprocessed measured value set of the atmospheric parameters to obtain a wind field prediction evaluation value;
[0016] S33, performing environmental prediction evaluation processing on the preprocessed predicted value set of the atmospheric environment parameters and the preprocessed measured value set of the atmospheric parameters to obtain an environmental prediction evaluation value;
[0017] S34, using a preset weight vector to perform weighted summation processing on the position difference evaluation value, the wind field prediction evaluation value, and the environmental prediction evaluation value to obtain a prediction evaluation value.
[0018] The performing position difference evaluation processing on the position coordinate information of all discrete points to obtain a position difference evaluation value includes:
[0019] S311, obtaining the position coordinate information of all discrete points in the predicted value set of the atmospheric environment parameters and the measured value set of the atmospheric parameters;
[0020] S312. For each discrete point in the predicted value set of the atmospheric environment parameters, determine the corresponding discrete point in the measured value set of the atmospheric parameters, and calculate the fusion deviation value for the two discrete points to obtain the fusion deviation value of the discrete points in the predicted value set of the atmospheric environment parameters.
[0021] The calculation expression of the fusion deviation value is as follows:
[0022] β = |x0 - x1| / |x0 + x1| + |1 - y0 / y1| + |exp(z0 / z1) - χ| / |z0 + z1|,
[0023] c1 = arctan(|x0 - x1| / |y0 - y1|)
[0024] c2 = arctan(|x0 - x1| / |z0 - z1|)
[0025]
[0026] where χ is a preset calculation parameter, β is the position deviation value of the discrete point, (x0, y0, z0) is the position coordinate information of the discrete point in the predicted value set of the atmospheric environment parameters, (x1, y1, z1) is the position coordinate information of the discrete point corresponding to the position coordinate information (x0, y0, z0) in the measured value set of the atmospheric parameters; c1 and c2 are the angular deviation values of the discrete point; ε is the fusion deviation value of the discrete point.
[0027] S313. Calculate the average value of the fusion deviation values of all discrete points in the predicted value set of the atmospheric environment parameters to obtain the position difference evaluation value.
[0028] The environmental prediction and evaluation process for the preprocessed predicted value set of the atmospheric environment parameters and the preprocessed measured value set of the atmospheric parameters to obtain the environmental prediction and evaluation value includes:
[0029] S331. Obtain the predicted pressure value sequence, predicted temperature value sequence, and predicted air density value sequence in the predicted value sequence of the atmospheric environment parameters, and use the predicted pressure value sequence, predicted temperature value sequence, and predicted air density value sequence to construct an environment matrix to be evaluated.
[0030] S332. Obtain the measured pressure value sequence, measured temperature value sequence, and measured air density value sequence in the measured value sequence of the atmospheric environment parameters, and use the measured pressure value sequence, measured temperature value sequence, and measured air density value sequence to construct a standard environment matrix.
[0031] S333. Calculate the difference value between the environment matrix to be evaluated and the standard environment matrix to obtain the environmental prediction and evaluation value.
[0032] The calculation and processing of the difference value includes:
[0033] Subtract the to-be-evaluated environment matrix from the standard environment matrix to obtain a difference matrix;
[0034] Perform a first standard calculation on the difference matrix to obtain a first eigenmatrix;
[0035] The expression of the first standard calculation is:
[0036]
[0037] where m represents the row dimension of the difference matrix, and x ij represents the element in the i-th row and j-th column of the difference matrix, and z ij represents the element in the i-th row and j-th column of the first eigenmatrix;
[0038] Perform a second extraction on the first eigenmatrix to obtain an optimal solution vector;
[0039] The second extraction is to take out the largest number in each column to form an optimal solution vector z + ; the expression of the optimal solution vector is:
[0040]
[0041] where n represents the column dimension of the difference matrix;
[0042] Perform a first extraction on the first eigenmatrix to obtain a worst solution vector;
[0043] The first extraction is to take out the smallest number in each column to form a worst solution vector z - ; the expression of the worst solution vector is:
[0044]
[0045] Perform an evaluation calculation on the worst solution vector and the optimal solution vector to obtain an environmental prediction evaluation value; the expression of the evaluation calculation is:
[0046]
[0047] In the formula, ω j is the preset j-th importance weight; ω j is obtained by pre-setting or by calculating the variance value of each column of the first eigenmatrix.
[0048] In the second aspect of the embodiments of the present application, an atmospheric environment prediction and evaluation device for the entire orbit of an aircraft is disclosed. The device includes:
[0049] A memory storing executable program code;
[0050] A processor coupled to the memory;
[0051] The processor calls the executable program code stored in the memory to execute the atmospheric environment prediction and evaluation method for the entire orbit of the aircraft.
[0052] In the third aspect of the embodiments of the present application, a computer-readable storage medium is disclosed. The computer-readable storage medium stores computer instructions, which are used to execute the atmospheric environment prediction and evaluation method for the entire orbit of the aircraft when called by a computer.
[0053] In the fourth aspect of the embodiments of the present application, an information data processing terminal is disclosed. The information data processing terminal is used to implement the atmospheric environment prediction and evaluation method for the entire orbit of the aircraft.
[0054] The beneficial effects of the present invention are as follows:
[0055] The present invention realizes the effective evaluation of the atmospheric environment prediction parameters for the entire flight orbit of the aircraft, which is of great significance for improving the atmospheric environment prediction method for the entire orbit of the aircraft and enhancing the prediction accuracy.
[0056] During the prediction and evaluation process of the present invention, different evaluation models are proposed for three types of data respectively, and finally the evaluation results are fused to improve the accuracy of the evaluation results.
[0057] For the measured values of air pressure, temperature, and air density, the present invention effectively extracts the subtle difference features of various types of data by jointly constructing a matrix and extracting the features of the difference matrix, and realizes the effective evaluation of the prediction effects of the three parameters of air pressure, temperature, and air based on the subtle difference features. The present invention effectively evaluates the position error of the predicted value from two aspects: angle and position.
[0058] Starting from the data characteristics of the wind field, the present invention introduces image data to characterize the wind field difference value, and realizes the effective evaluation of the prediction accuracy of the wind field data through the value range distribution of pixel values. BRIEF DESCRIPTION OF THE DRAWINGS
[0059] Figure 1 It is a flowchart of the implementation of the method of the present invention. DETAILED DESCRIPTION
[0060] To better understand the content of the present invention, an embodiment is given here.
[0061] Figure 1 This is the implementation flowchart of the method of the present invention.
[0062] In the first aspect of the embodiments of the present application, an atmospheric environment prediction and evaluation method for the entire orbit of an aircraft is disclosed, including:
[0063] S1. Obtain a set of predicted values of atmospheric environment parameters and a set of measured values of atmospheric parameters for the aircraft orbit; the set of predicted values of atmospheric environment parameters includes position coordinate information and corresponding sequences of predicted values of atmospheric environment parameters at several discrete points on the aircraft orbit; the set of measured values of atmospheric parameters includes position coordinate information and corresponding sequences of measured values of atmospheric environment parameters at several discrete points on the aircraft orbit;
[0064] The discrete points in the set of measured values of atmospheric parameters and the set of predicted values of atmospheric environment parameters are in one-to-one correspondence, that is, each sequence of predicted values of atmospheric environment parameters has a corresponding sequence of measured values of atmospheric environment parameters at the same location in the set of predicted values of atmospheric environment parameters. Specifically, the sequences in the set of measured values of atmospheric parameters and the set of predicted values of atmospheric environment parameters can be established in one-to-one correspondence according to the sequence numbers in the set.
[0065] S2. Perform preprocessing on the set of predicted values of atmospheric environment parameters and the set of measured values of atmospheric parameters respectively to obtain a preprocessed set of predicted values of atmospheric environment parameters and a preprocessed set of measured values of atmospheric parameters;
[0066] S3. Perform prediction and evaluation processing on the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters to obtain a prediction and evaluation value; the prediction and evaluation value is used to evaluate the prediction accuracy of the set of predicted values of atmospheric environment parameters.
[0067] The sequence of predicted values of atmospheric environment parameters includes a sequence of predicted wind field values, a sequence of predicted pressure values, a sequence of predicted temperature values, and a sequence of predicted air density values;
[0068] The sequence of measured values of atmospheric environment parameters includes a sequence of measured wind field values, a sequence of measured pressure values, a sequence of measured temperature values, and a sequence of measured air density values;
[0069] The performing preprocessing on the set of predicted values of atmospheric environment parameters and the set of measured values of atmospheric parameters respectively to obtain a preprocessed set of predicted values of atmospheric environment parameters and a preprocessed set of measured values of atmospheric parameters includes:
[0070] S21. Perform outlier rejection processing on the set of predicted values of atmospheric environment parameters and the set of measured values of atmospheric parameters respectively to obtain a first set of predicted values of atmospheric environment parameters and a first set of measured values of atmospheric parameters;
[0071] S22. Normalize the first set of predicted values of atmospheric environment parameters and the first set of measured values of atmospheric parameters respectively to obtain a preprocessed set of predicted values of atmospheric environment parameters and a preprocessed set of measured values of atmospheric parameters;
[0072] Performing prediction and evaluation processing on the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters to obtain a prediction and evaluation value, including:
[0073] S31. Extract the position coordinate information of all discrete points in the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters, and perform position difference evaluation processing on the position information of all discrete points to obtain a position difference evaluation value;
[0074] S32. Perform wind field prediction and evaluation processing on the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters to obtain a wind field prediction and evaluation value;
[0075] S33. Perform environmental prediction and evaluation processing on the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters to obtain an environmental prediction and evaluation value;
[0076] S34. Use a preset weight vector to perform weighted summation processing on the position difference evaluation value, the wind field prediction and evaluation value, and the environmental prediction and evaluation value to obtain a prediction and evaluation value.
[0077] The set of predicted values of atmospheric environment parameters includes the position coordinate information of several discrete points on the aircraft orbit and the corresponding sequence of predicted values of atmospheric environment parameters; the set of measured values of atmospheric parameters includes the position coordinate information of several discrete points on the aircraft orbit and the corresponding sequence of measured values of atmospheric environment parameters;
[0078] Performing position difference evaluation processing on the position information of all discrete points to obtain a position difference evaluation value, including:
[0079] S311. Obtain the position coordinate information of all discrete points in the set of predicted values of atmospheric environment parameters and the set of measured values of atmospheric parameters;
[0080] S312. For each discrete point in the set of predicted values of atmospheric environment parameters, determine the corresponding discrete point in the set of measured values of atmospheric parameters, and calculate the fusion deviation value of the two discrete points to obtain the fusion deviation value of the discrete points in the set of predicted values of atmospheric environment parameters;
[0081] The calculation expression of the fusion deviation value is:
[0082] β = |x0 - x1| / |x0 + x1| + |1 - y0 / y1| + |exp(z0 / z1) - χ| / |z0 + z1|,
[0083] c1 = arctan(|x0 - x1| / |y0 - y1|)
[0084] c2 = arctan(|x0 - x1| / |z0 - z1|)
[0085]
[0086] where χ is a preset calculation parameter, β is the position deviation value of the discrete points, (x0, y0, z0) is the position coordinate information of the discrete points in the set of predicted values of the atmospheric environment parameters, (x1, y1, z1) is the position coordinate information of the discrete points in the set of measured values of the atmospheric parameters; c1 and c2 are the angular deviation values of the discrete points; ε is the fusion deviation value of the discrete points.
[0087] S313. Calculate the average value of the fusion deviation values of all discrete points in the set of predicted values of the atmospheric environment parameters to obtain the position difference evaluation value.
[0088] The environmental prediction and evaluation process for the preprocessed set of predicted values of the atmospheric environment parameters and the preprocessed set of measured values of the atmospheric parameters to obtain the environmental prediction and evaluation value includes:
[0089] S331. Obtain the predicted air pressure value sequence, predicted temperature value sequence, and predicted air density value sequence in the sequence of predicted values of the atmospheric environment parameters, and use the predicted air pressure value sequence, predicted temperature value sequence, and predicted air density value sequence to construct the environment matrix to be evaluated;
[0090] S332. Obtain the measured air pressure value sequence, measured temperature value sequence, and measured air density value sequence in the sequence of measured values of the atmospheric environment parameters, and use the measured air pressure value sequence, measured temperature value sequence, and measured air density value sequence to construct the standard environment matrix;
[0091] S333. Calculate the difference value between the environment matrix to be evaluated and the standard environment matrix to obtain the environmental prediction and evaluation value.
[0092] The difference value calculation process includes:
[0093] Subtract the standard environment matrix from the environment matrix to be evaluated to obtain the difference matrix;
[0094] Perform the first standard calculation process on the difference matrix to obtain the first eigenmatrix;
[0095] The expression of the first standard calculation process is:
[0096]
[0097] Among them, m represents the row dimension of the difference matrix, and x ij represents the element in the i-th row and j-th column of the said difference matrix, and z ij represents the element in the i-th row and j-th column of the first feature matrix;
[0098] Perform a second extraction process on the said first feature matrix to obtain an optimal solution vector;
[0099] The said second extraction process is to take out the largest number in each column to form the optimal solution vector z + ; The expression of the said optimal solution vector is:
[0100]
[0101] Among them, n represents the column dimension of the difference matrix;
[0102] Perform a first extraction process on the said first feature matrix to obtain a worst solution vector;
[0103] The said first extraction process is to take out the smallest number in each column to form the worst solution vector z - ; The expression of the said worst solution vector is:
[0104]
[0105] Perform an evaluation calculation process on the said worst solution vector and the optimal solution vector to obtain an environmental prediction evaluation value; The expression of the said evaluation calculation process is:
[0106]
[0107] In the formula, ω j is the preset j-th importance weight; ω j is obtained by pre-setting, or by calculating the variance value of each column of the first feature matrix.
[0108] The said process of performing wind field prediction evaluation on the preprocessed atmospheric environment parameter prediction value set and the preprocessed atmospheric parameter measurement value set to obtain a wind field prediction evaluation value includes:
[0109] S321. Obtain the wind field prediction value sequence in the said atmospheric environment parameter prediction value sequence and the wind field measurement value sequence in the said atmospheric environment parameter measurement value sequence;
[0110] S322. Perform an absolute value difference process on the said wind field prediction value sequence and the wind field measurement value sequence to obtain a wind field difference sequence;
[0111] S323. Obtain the position coordinate information of the discrete points corresponding to each element of the wind field difference sequence;
[0112] S324. Using the position coordinate information of the discrete points of each element as the image coordinates, and using the value of the element as the gray value of the pixel point corresponding to the image coordinates, construct the wind field image information; the wind field image information includes the gray value of the pixel point and the image coordinates of the pixel point;
[0113] S325. Uniformly divide the wind field image information into several sub-images;
[0114] S326. Perform an evaluation process on all sub-images to obtain a wind field prediction evaluation value.
[0115] The preset gray value range can be [1,10], [11,22], [23,40], [41,50], [50,255].
[0116] The expression of the evaluation process is:
[0117]
[0118] where c iv is the number of pixel points whose gray value belongs to the v-th preset gray value range in the i-th sub-image, p iv is the occurrence probability of the v-th gray value range in the i-th sub-image, V is the total number of preset gray scale value ranges, W and H are the length value and width value of the sub-image respectively, h i is the intermediate measurement value of the i-th sub-image, L is the number of sub-images, and u is the wind field prediction evaluation value.
[0119] Specifically, W and H can be respectively represented by the number of pixel points of the sub-image in the length direction and the width direction.
[0120] The position coordinate information of the discrete points corresponding to each element of the wind field difference sequence is the position coordinate information of the discrete points corresponding to the wind field prediction value corresponding to the element.
[0121] The preset weight vector can be (0.2, 0.3, 0.5).
[0122] The outlier rejection process is to judge whether each data is within the preset threshold, and delete the data that is not within the preset threshold from the set of atmospheric environment parameter prediction values or the set of atmospheric parameter measurement values.
[0123] The normalization process is to obtain the maximum value of each numerical sequence, including the predicted value sequence and the measured value sequence, for each type of parameter, and divide the numerical values in the numerical sequence by the maximum value to obtain the normalized numerical sequence.
[0124] Taking the position coordinate information of the discrete points of each element as the image coordinates means taking the planar coordinates (x, y) in the position coordinates as the image coordinates.
[0125] In the second aspect of the embodiments of the present application, an atmospheric environment prediction and evaluation device for the full orbit of an aircraft is disclosed. The device includes:
[0126] A memory storing executable program code;
[0127] A processor coupled to the memory;
[0128] The processor calls the executable program code stored in the memory to execute the above-mentioned atmospheric environment prediction and evaluation method for the full orbit of the aircraft.
[0129] In the third aspect of the embodiments of the present application, a computer-readable storage medium is disclosed. The computer-readable storage medium stores computer instructions, which are used to execute the above-mentioned atmospheric environment prediction and evaluation method for the full orbit of the aircraft when called by a computer.
[0130] In the fourth aspect of the embodiments of the present application, an information data processing terminal is disclosed. The information data processing terminal is used to implement the above-mentioned atmospheric environment prediction and evaluation method for the full orbit of the aircraft.
[0131] The above are only the 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 predicting and evaluating the atmospheric environment of an aircraft's entire orbit, characterized in that, Including: S1. Obtain a set of predicted values of atmospheric environment parameters and a set of measured values of atmospheric parameters for the flight vehicle orbit; the set of predicted values of atmospheric environment parameters includes the position coordinate information of several discrete points on the flight vehicle orbit and the corresponding sequence of predicted values of atmospheric environment parameters; the set of measured values of atmospheric parameters includes the position coordinate information of several discrete points on the flight vehicle orbit and the corresponding sequence of measured values of atmospheric environment parameters; S2. Perform preprocessing on the set of predicted values of atmospheric environment parameters and the set of measured values of atmospheric parameters respectively to obtain a preprocessed set of predicted values of atmospheric environment parameters and a preprocessed set of measured values of atmospheric parameters; S3. Perform prediction evaluation processing on the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters to obtain a prediction evaluation value; the prediction evaluation value is used to evaluate the prediction accuracy of the set of predicted values of atmospheric environment parameters; The sequence of predicted values of atmospheric environment parameters includes a sequence of predicted values of wind field, a sequence of predicted values of air pressure, a sequence of predicted values of temperature, and a sequence of predicted values of air density; The sequence of measured values of atmospheric environment parameters includes a sequence of measured values of wind field, a sequence of measured values of air pressure, a sequence of measured values of temperature, and a sequence of measured values of air density; The performing preprocessing on the set of predicted values of atmospheric environment parameters and the set of measured values of atmospheric parameters respectively to obtain a preprocessed set of predicted values of atmospheric environment parameters and a preprocessed set of measured values of atmospheric parameters includes: S21. Perform outlier rejection processing on the set of predicted values of atmospheric environment parameters and the set of measured values of atmospheric parameters respectively to obtain a first set of predicted values of atmospheric environment parameters and a first set of measured values of atmospheric parameters; S22. Perform normalization processing on the first set of predicted values of atmospheric environment parameters and the first set of measured values of atmospheric parameters respectively to obtain a preprocessed set of predicted values of atmospheric environment parameters and a preprocessed set of measured values of atmospheric parameters; The performing prediction evaluation processing on the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters to obtain a prediction evaluation value includes: S31. Extract the position coordinate information of all discrete points in the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters, and perform position difference evaluation processing on the position coordinate information of all discrete points to obtain a position difference evaluation value; S32. Perform wind field prediction evaluation processing on the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters to obtain a wind field prediction evaluation value; S33. Perform environmental prediction evaluation processing on the preprocessed set of predicted values of atmospheric environment parameters and the preprocessed set of measured values of atmospheric parameters to obtain an environmental prediction evaluation value; S34. Use a preset weight vector to perform weighted summation processing on the position difference evaluation value, the wind field prediction evaluation value, and the environmental prediction evaluation value to obtain a prediction evaluation value; The performing position difference evaluation processing on the position coordinate information of all discrete points to obtain a position difference evaluation value includes: S311. Obtain the position coordinate information of all discrete points in the set of predicted values of atmospheric environment parameters and the set of measured values of atmospheric parameters; S312. For each discrete point in the predicted value set of the atmospheric environment parameters, determine the corresponding discrete point in the measured value set of the atmospheric parameters, and calculate the fusion deviation value of the two discrete points to obtain the fusion deviation value of the discrete points in the predicted value set of the atmospheric environment parameters; The calculation expression of the fusion deviation value is: β = |x0 - x1| / |x0 + x1| + |1 - y0 / y1| + |exp(z0 / z1) - χ| / |z0 + z1|, c1 = arctan(|x0 - x1| / |y0 - y1|) c2 = arctan(|x0 - x1| / |z0 - z1|), where χ is a preset calculation parameter, β is the position deviation value of the discrete point, (x0, y0, z0) is the position coordinate information of the discrete point in the predicted value set of the atmospheric environment parameters, (x1, y1, z1) is the position coordinate information of the discrete point corresponding to the position coordinate information (x0, y0, z0) in the measured value set of the atmospheric parameters; c1 and c2 are the angular deviation values of the discrete point; ε is the fusion deviation value of the discrete point; S313. Calculate the mean value of the fusion deviation values of all discrete points in the predicted value set of the atmospheric environment parameters to obtain the position difference evaluation value.
2. The method for predicting and evaluating the atmospheric environment of the entire orbit of an aircraft according to claim 1, characterized in that, The environmental prediction and evaluation process for the preprocessed predicted value set of the atmospheric environment parameters and the preprocessed measured value set of the atmospheric parameters to obtain the environmental prediction and evaluation value includes: S331. Obtain the predicted value sequence of air pressure, the predicted value sequence of temperature, and the predicted value sequence of air density in the predicted value sequence of the atmospheric environment parameters, and use the predicted value sequence of air pressure, the predicted value sequence of temperature, and the predicted value sequence of air density to construct an environment matrix to be evaluated; S332. Obtain the measured value sequence of air pressure, the measured value sequence of temperature, and the measured value sequence of air density in the measured value sequence of the atmospheric environment parameters, and use the measured value sequence of air pressure, the measured value sequence of temperature, and the measured value sequence of air density to construct a standard environment matrix; S333. Calculate the difference value between the environment matrix to be evaluated and the standard environment matrix to obtain the environmental prediction and evaluation value.
3. The method for predicting and evaluating the atmospheric environment of the entire orbit of an aircraft according to claim 2, characterized in that The difference value calculation process includes: Subtract the standard environment matrix from the environment matrix to be evaluated to obtain a difference matrix; Perform a first standard calculation process on the difference matrix to obtain a first eigenmatrix; The expression of the first standard calculation process is: where m represents the row dimension of the difference matrix, and x ij represents the element in the i-th row and j-th column of the difference matrix, and z ij represents the element in the i-th row and j-th column of the first feature matrix; Perform a second extraction process on the first eigenmatrix to obtain an optimal solution vector; The second extraction process is to take out the largest number in each column to form the optimal solution vector z + ; The expression of the optimal solution vector is as follows: where n represents the column dimension of the difference matrix; Perform a first extraction process on the first eigenmatrix to obtain a worst solution vector; The first extraction process is to take out the smallest number in each column to form the worst solution vector z - ; The expression of the worst solution vector is as follows: Perform an evaluation calculation process on the worst solution vector and the optimal solution vector to obtain the environmental prediction and evaluation value; the expression of the evaluation calculation process is: where ω j is the preset j-th importance weight; ω j is obtained by presetting or by calculating the variance value of each column of the first feature matrix.
4. The atmospheric environment prediction and evaluation method for the full orbit of an aircraft according to claim 1, characterized in that, The wind field prediction and evaluation process for the preprocessed predicted value set of the atmospheric environment parameters and the preprocessed measured value set of the atmospheric parameters to obtain the wind field prediction and evaluation value includes: S321. Obtain the wind field prediction value sequence in the predicted value sequence of the atmospheric environment parameters and the wind field measurement value sequence in the measured value sequence of the atmospheric environment parameters; S322. Perform an absolute value difference operation on the wind field prediction value sequence and the wind field measurement value sequence to obtain a wind field difference sequence; S323. Obtain the position coordinate information of the discrete points corresponding to each element of the wind field difference sequence; S324. Use the position coordinate information of the discrete points of each element as the image coordinates, and use the value of the element as the gray value of the pixel point corresponding to the image coordinates to construct wind field image information; the wind field image information includes the gray value of the pixel point and the image coordinates of the pixel point; S325. Uniformly divide the wind field image information into a number of sub-images; S326. Perform an evaluation process on all sub-images to obtain a wind field prediction evaluation value; The expression of the evaluation process is: Among them, c iv is the number of pixel points in the i-th sub-image whose gray values belong to the v-th preset gray value range, and p iv is the occurrence probability of the v-th gray value range in the i-th sub-image. V is the total number of preset gray scale value ranges. W and H are the length value and width value of the sub-image respectively, and h i is the intermediate metric value of the i-th sub-image, L is the number of sub-images, and u is the wind field prediction evaluation value.
5. An atmospheric environment prediction and evaluation device for the entire orbit of an aircraft, characterized in that, The device includes: 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 atmospheric environment prediction and evaluation method for the entire orbit of the aircraft according to any one of claims 1 to 4.
6. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions, which are used to execute the atmospheric environment prediction and evaluation method for the entire orbit of the aircraft according to any one of claims 1 to 4 when called by a computer.
7. An information data processing terminal, characterized in that, The information data processing terminal is used to implement the atmospheric environment prediction and evaluation method for the entire orbit of the aircraft according to any one of claims 1 to 4.
Citation Information
Patent Citations
Prediction accuracy evaluation device and method of the same
JP2015132914A