Method and system for obtaining fine structure by improving height resolution of middle and upper atmospheric wind field

Abstract

The invention relates to radio physics and radar signal processing technologies, and particularly discloses a method and a system for obtaining a fine structure by improving the height resolution of a middle and upper atmospheric wind field, and the method comprises the steps: screening meteors participating in a wind field inversion program, extracting an echo signal, an elevation angle, an azimuth angle and a Doppler velocity of each screened meteor; the screened meteor echo signals are converted into a frequency spectrum, a distance resolution improvement multiple is set, the accurate group distance of the meteor echo signals is recalculated through the Capon algorithm, the effect of improving the meteor group distance resolution is achieved, the corresponding height is calculated, and according to the azimuth angle, the elevation angle and the Doppler velocity of each meteor, the Doppler velocity of each meteor is calculated. And substituting the height into a wind field inversion program to obtain a refined wind field with higher distance resolution. According to the method, extra equipment and hardware cost are not needed, the distance resolution of meteor detection can be greatly improved, and meanwhile a refined wind field result of a middle-upper atmosphere area can be obtained.

Description

Technical Field

[0001] The present invention relates to radio physics and radar signal processing technologies, and particularly to a method and system for obtaining fine structures by improving the height resolution of the middle and upper atmosphere wind field. Background Art

[0002] The middle and upper atmosphere includes the stratosphere, mesosphere and thermosphere, and the mesosphere and lower thermosphere are the middle and upper atmosphere regions 60 - 110 km away from the ground. Long-term observations of parameters such as the atmospheric wind field, temperature, and density in this region play an important role in understanding the changes and coupling relationships of the Earth's atmospheric spheres. Ground-based detection technologies, especially radio radars and lidars, have become important means for detecting the atmosphere in this region. Meteor radars, due to their advantages such as low cost and being unaffected by weather, are the most widely used detection methods. By detecting the Doppler drift velocity of meteor plasma wakes using radio radars, meteor echoes can be used as a medium to detect the mesosphere atmospheric wind field. Currently, wind-measuring meteor radars, i.e., all-sky meteor radars, have been widely applied.

[0003] Based on the distribution of underdense meteors detected by all-sky meteor radars, the height range of wind field inversion is generally 70 - 110 km, the time resolution is 1 hour, and limited by the range resolution of meteor radars, the height resolution of the wind field inversion is generally 2 km. Improving the range resolution of meteor radars in hardware will force the device's bandwidth to increase, introducing more interference from non-target frequencies, which is not conducive to meteor detection. Currently, there is a lack of effective methods to improve the range resolution. Summary of the Invention

[0004] To solve the problem that improving the range resolution of meteor radars in hardware will force the device's bandwidth to increase and introduce more interference from non-target frequencies, which is not conducive to meteor detection, the present invention provides a method and system for obtaining fine structures by improving the height resolution of the middle and upper atmosphere wind field. By steps such as screening meteors, converting signals, calculating covariance matrices, and applying the Capon algorithm, the height resolution of the middle and upper atmosphere wind field is improved.

[0005] According to one aspect of the specification of the present invention, a method for obtaining fine structures by improving the height resolution of the middle and upper atmosphere wind field is provided, including:

[0006] Screening out the meteors participating in the wind field inversion program, and extracting the echo signals, elevation angles, azimuth angles, and Doppler velocities of each screened meteor;

[0007] Converting the echo signal of each meteor from the time domain spectrum to the frequency spectrum, and intercepting the frequency matrix at the position where the meteor echo signal is located from the frequency spectrum;

[0008] Calculate the new range resolution after increasing the target multiple according to the original range resolution, calculate the covariance matrix of the meteor according to the frequency matrix, substitute the covariance matrix into the Capon algorithm to calculate the pseudo-spectrum at the new range resolution, obtain the accurate group distance of each meteor according to the pseudo-spectrum, and calculate the accurate height of each meteor according to the obtained accurate group distance;

[0009] Substitute the elevation angle, azimuth angle, Doppler velocity and accurate height of each meteor into the wind field inversion program to recalculate the wind field, and obtain the refined wind field at the new range resolution.

[0010] As a further technical solution, the meteors participating in the wind field inversion are screened out, and the echo signal of each meteor is extracted, including: tracing back the inversion process of the wind field, reversely calibrating the meteors participating in the wind field calculation according to the corresponding height and time information, and screening out the time-domain signal of the meteor echo.

[0011] As a further technical solution, the frequency matrix at the position of the meteor echo signal is intercepted from the frequency spectrum, including: preprocessing the time-domain spectrum of the meteor echo signal to initially determine the position of the signal, and intercepting a new frequency matrix.

[0012] As a further technical solution, calculating the covariance matrix of the meteor according to the frequency matrix, including:

[0013] Calculate the duration of each chip of the meteor radar modulation coding through the original range resolution inherent in the meteor radar detection device, and calculate the signal bandwidth through the chip duration;

[0014] According to the new frequency matrix, perform slicing in the range dimension, and list the frequency spectrum sequences under each range gate as sub-frequency matrices;

[0015] Calculate the frequency component interval of the sub-frequency matrix through the signal bandwidth and the length of the original frequency spectrum dimension;

[0016] Starting from the first frequency component interval, divide and reconstruct the frequency spectrum to obtain the frequency matrix;

[0017] Obtain the covariance matrix of each frequency component through the reconstructed frequency matrix.

[0018] As a further technical solution, substitute the covariance matrix into the Capon algorithm to calculate the pseudo-spectrum at the new range resolution, and obtain the accurate group distance of each meteor according to the pseudo-spectrum, including:

[0019] Set the target multiple for increasing the range resolution and then calculate the increased new range resolution;

[0020] Construct the steering matrix of the frequency matrix according to the increased target multiple;

[0021] Substitute the covariance matrix and the guiding matrix into the Capon algorithm to obtain the pseudo-spectrum of the new range resolution.

[0022] Analyze the pseudo-spectrum along the range dimension to obtain the accurate group range of the meteor echo signal.

[0023] As a further technical solution, calculate the accurate altitude of the meteor according to the accurate group range, including: calculate the accurate altitude of the meteor according to the accurate group range and the elevation angle of the meteor obtained from the initial wind field inversion program. By introducing the elevation angle of the meteor, the altitude of the meteor can be calculated more accurately, providing more accurate input data for subsequent wind field inversion.

[0024] As a further technical solution, during the process of screening the meteors participating in the wind field calculation, information such as the azimuth angle, elevation angle, and Doppler velocity of the meteors is obtained simultaneously. Substitute the accurate altitude into the wind field calculation program, and re-divide the altitude interval with the new range resolution to obtain a refined wind field.

[0025] According to one aspect of the specification of the present invention, the present invention provides a system for obtaining fine structure by improving the altitude resolution of the middle and upper atmosphere wind field, including:

[0026] The first main module is used to screen out the meteors participating in the calculation from the wind field inversion program, and extract their echo signals, elevation angles, azimuth angles, and Doppler velocities.

[0027] The second main module is used to convert the echo signal of each meteor from the time-domain spectrum to the frequency spectrum, and intercept the frequency matrix at the position where the meteor echo signal is located from the frequency spectrum.

[0028] The third main module is used to calculate the new range resolution after increasing the target multiple according to the original range resolution, calculate the covariance matrix of the meteor according to the frequency matrix, substitute the covariance matrix into the Capon algorithm to calculate the pseudo-spectrum at the new range resolution, obtain the accurate group range of each meteor according to the pseudo-spectrum, and calculate the accurate altitude of each meteor according to the obtained accurate group range.

[0029] The fourth main module is used to substitute the elevation angle, azimuth angle, Doppler velocity, and accurate altitude of each meteor into the wind field inversion program to recalculate the wind field, and obtain a refined wind field at the new range resolution. The second main module converts the meteor echo signal into a frequency spectrum and intercepts the frequency matrix, then the third main module extracts the accurate group range and accurate altitude of the signal and inputs them into the fourth main module, thereby quickly calculating the refined middle and upper atmosphere wind field.

[0030] According to one aspect of the specification of the present invention, the present invention provides an electronic device, including a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the steps of the method for obtaining a fine structure by enhancing the height resolution of the middle and upper atmospheric wind field are implemented.

[0031] According to one aspect of the specification of the present invention, the present invention provides a computer-readable storage medium that stores a computer program. When the computer program is executed by a processor, the steps of the method for obtaining a fine structure by enhancing the height resolution of the middle and upper atmospheric wind field are implemented.

[0032] The beneficial technical effects of the present invention are as follows: Without changing the hardware parameters of the meteor radar, by performing signal processing and analysis on the frequency spectrum of the meteor echo signal, the frequency matrix at the position of the meteor echo signal for each detection is intercepted, reducing the amount of arithmetic processing. The precise group distance of the meteor echo signal can be obtained according to the required magnification factor, and then the precise height of the meteor can be obtained. Through signal processing, the range resolution of the meteor echo is improved, and at the same time, the height interval of the wind field becomes smaller, so that a middle and upper atmospheric wind field with a smaller height resolution can be obtained. The reduction of the height resolution can reveal the fine structure of the wind field and obtain more detailed wind field change information. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings used in the description of the embodiments or the prior art. Obviously, the following drawings are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can also be obtained based on these drawings without creative efforts.

[0034] Figure 1 It is a schematic diagram of the overall process of a method for obtaining a fine structure by enhancing the height resolution of a middle and upper atmospheric wind field according to the present invention;

[0035] Figure 2 It is an overall diagram of the time-domain signal at the original range resolution provided by an embodiment of the present invention;

[0036] Figure 3 It is a detailed diagram of the region where the signal is located at the original range resolution provided by an embodiment of the present invention;

[0037] Figure 4 It is an overall diagram of the time-domain signal after enhancing the range resolution provided by an embodiment of the present invention;

[0038] Figure 5 It is a detailed diagram of the region where the signal is located after enhancing the range resolution provided by an embodiment of the present invention;

[0039] Figure 6This is the local wind field map within the altitude range of 84 - 90 km with the original altitude resolution of 2 km provided by the embodiments of the present invention;

[0040] Figure 7 This is the local wind field map within the altitude range of 84 km - 90 km with the resolution of 500 m after the altitude resolution is increased by 4 times provided by the embodiments of the present invention;

[0041] Figure 8 This is the system block diagram for obtaining the fine structure by enhancing the altitude resolution of the middle and upper atmosphere wind field provided by the embodiments of the present invention. Detailed implementation manners

[0042] Next, the technical solutions in the embodiments of the present invention will be described clearly and completely in conjunction with the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0043] It should be noted that, without conflict, the embodiments in the present invention and the features in the embodiments may be combined with each other.

[0044] Next, the present invention will be further described in conjunction with specific embodiments, but it is not a limitation of the present invention.

[0045] This embodiment provides a method for obtaining the fine structure by enhancing the altitude resolution of the middle and upper atmosphere wind field. The echo signals of the meteors participating in the wind field inversion are screened out from the meteor echo signal data within one day, the screened signals are converted into frequency spectra, a distance resolution enhancement multiple is set, the group distance of the meteor echo signals is recalculated using the Capon algorithm to achieve the effect of enhancing the meteor group distance resolution and calculate the corresponding altitude, and according to the azimuth, elevation angle and Doppler velocity of each meteor, the altitude is substituted into the wind field inversion program to obtain a refined wind field with a higher distance resolution. The specific steps are as follows:

[0046] S1. Screen out the meteors participating in the wind field inversion program, and extract the echo signals, elevation angles, azimuth angles and Doppler velocities of each screened meteor;

[0047] S2. Convert the echo signal of each meteor from the time domain spectrum to the frequency spectrum, and intercept the frequency matrix at the position of the meteor echo signal from the frequency spectrum;

[0048] S3. Calculate the new range resolution after increasing the target multiple according to the original range resolution, calculate the covariance matrix of meteors according to the frequency matrix, substitute the covariance matrix into the Capon algorithm to calculate the pseudo-spectrum at the new range resolution, obtain the precise group distance of each meteor according to the pseudo-spectrum, and calculate the precise height of each meteor according to the obtained precise group distance;

[0049] S4. Substitute the elevation angle, azimuth angle, Doppler velocity, and precise height of each meteor into the wind field inversion program to recalculate the wind field, and obtain the refined wind field at the new range resolution.

[0050] In this embodiment, the method for screening meteors in step S1 is as follows: Trace back the inversion process of the wind field, reverse-calibrate the meteors participating in the wind field calculation according to the corresponding height and time information, and screen out the time-domain signals of meteor echoes.

[0051] It should be noted that during the calculation of the original wind field, the azimuth angle, elevation angle, and Doppler velocity of each meteor have been calculated, and the azimuth angle, elevation angle, and Doppler velocity of the screened meteors can be directly referenced in the subsequent steps.

[0052] In step S2, for the echo signal of each meteor, convert the time-domain spectrum to a frequency spectrum. The time-domain spectrum of the screened meteor echo signal is converted to a frequency spectrum through Fourier transform , and the forms and relationships of the time-domain spectrum and the frequency spectrum are as follows:

[0053] (1),

[0054] (2),

[0055] (3),

[0056] It should be noted that in (1) and (2), m represents the distance dimension value of the signal spectrum, n represents the time / frequency dimension value of the signal spectrum, and a and b are the values in the matrix forms of the time-domain spectrum and the frequency spectrum respectively.

[0057] In step S2, set the required multiple of range resolution improvement, calculate the covariance matrix corresponding to the frequency spectrum in step 2, and use the Capon algorithm to calculate the pseudo-spectrum at the new range resolution to obtain the precise group distance of each meteor after the range resolution is improved. The specific method is as follows:

[0058] 1) The inherent original range resolution d0 of the meteor radar detection device is:

[0059] (4),

[0060] Where c is the speed of light, The signal bandwidth B for detection is:

[0061] (5),

[0062] 2) Through the time domain spectrum of meteor echo signal Preprocessing, preliminarily determine the location of the signal, and extract a new frequency matrix , which can be expressed as follows:

[0063] (6),

[0064] Among them, k and v represent the starting and ending values in the distance dimension of the original time domain spectrum, g and u represent the starting and ending values in the frequency dimension of the original time domain spectrum, that is, the newly intercepted frequency matrix is (vk) rows and (ug) columns;

[0065] 3) Frequency matrix Using the distance dimension to split, the frequency spectrum sequence under each distance gate can be expressed as:

[0066] (7),

[0067] in, is the frequency spectrum sequence under a single range gate, 、 are the starting frequency and ending frequency of the frequency matrix respectively, is the spectral component of the starting frequency in the frequency spectrum sequence, is the spectral component of the termination frequency in the frequency spectrum sequence;

[0068] The frequency component interval of the signal frequency spectrum sequence with bandwidth B is: (8),

[0069] Where n is the length (value) of the frequency dimension of the original frequency spectrum. is the frequency component interval of the submatrix, starting from the first frequency component, the frequency spectrum is divided and the frequency matrix is reconstructed:

[0070] (9),

[0071] Where N is the number of subbands, M is equal to ug, , Therefore The first frequency component divided by the frequency interval, the rest to They are the second to the last frequency components.

[0072] 4) Set the distance resolution improvement factor S, then the improved resolution is , and construct the frequency matrix according to the improvement factor S of the guiding matrix:

[0073] (10),

[0074] where the phase slope vector is:

[0075] (11),

[0076] is a mathematical expression, e is the natural logarithm, j is the imaginary unit, , ... are all phase slopes;

[0077] 5) The covariance matrix of each frequency component in the frequency spectrum can be expressed as:

[0078] (12),

[0079] where, represents of the conjugate transpose.

[0080] 6) Substitute the covariance matrix and the guiding matrix into the Capon algorithm to obtain the distance - dimension pseudo - spectrum with the resolution improved to :

[0081] (13),

[0082] In the above formula, is of the inverse matrix, C is the guiding matrix, is the transpose matrix of C.

[0083] By performing signal analysis on the pseudo - spectrum along the distance dimension, the accurate group distance of the meteor echo signal can be obtained .

[0084] Performing signal analysis on the pseudo - spectrum along the distance dimension specifically includes multiplying the number of rows in the distance dimension of the meteor echo signal in the pseudo - spectrum by the distance resolution to obtain the accurate group distance of the corresponding meteor echo signal.

[0085] In step S4, through the accurate group distance Calculate the exact altitude of the meteor , and the specific method is as follows:

[0086] (14),

[0087] It should be noted that the elevation angle of the meteor can be obtained together with other parameters of the meteor during the backtracking process of S1.

[0088] Moreover, in step S4, substitute the exact altitude of the meteor into the calculation program of the wind field to recalculate the wind field, and re-divide the altitude with as the resolution to obtain a refined wind field with a higher range resolution.

[0089] It should be noted that the wind field inversion program designed in this embodiment adopts a mature calculation program on the market.

[0090] Specifically, during implementation, the original range resolution is 2 km. In this implementation, the magnification factor is 4, and the improved range resolution should be 500 m.

[0091] As Figure 2-5 shown, in one implementation, it is a comparison diagram of the results of improving the resolution of the echo signal of a certain meteor in this application. Among them Figure 2 shows the time-domain signal diagram of this meteor under the original range resolution, Figure 3 is the detail diagram of the area where the meteor echo signal is located. It can be seen that the original range resolution in the "detection range" dimension is 2 km; Figure 4 is the time-domain signal diagram of this meteor after improving the range resolution, Figure 5 is the detail diagram of the area where the meteor echo signal is located. It can be seen that the improved range resolution in the "detection range" dimension is 500 m. From the comparison, it can be seen that the range resolution of the signal after processing has increased by 4 times.

[0092] As Figure 6-7 shown, it is a comparison diagram of the wind field results before and after improving the resolution. In order to clearly show the difference in altitude resolution, Figure 6-7 the pictures shown are all partial enlarged diagrams of the wind field details in a partial altitude range. As Figure 6 shown, it is a local wind field diagram in the altitude range of 84 - 90 km with the original altitude resolution of 2 km. It can be seen from the figure that at 14:00, there is a wind shear at altitudes of 84 km and 86 km, and the fine structure of this shear cannot be obtained with a 2-km altitude resolution. Similarly, at 15:00, there is a wind shear at altitudes of 84 km and 86 km, and the fine structure of this shear cannot be obtained with a 2-km altitude resolution; As Figure 7As shown, it is a local wind field map in the altitude range of 84 km - 90 km with a resolution of 500 m after the height resolution is increased by 4 times. It can be seen from the figure that at 14:00 and 15:00, the fine structure of the wind field gradually trimming at 84 km - 87 km can be clearly obtained. Compared with Figure 6 , the fine structure and change details of the wind shear in the wind field at high resolution can be clearly seen.

[0093] Based on the same inventive concept as the foregoing embodiments, the embodiments of the present invention further provide a system for obtaining the fine structure by improving the height resolution of the middle and upper atmosphere wind field, including:

[0094] The first main module is used to screen out the meteors participating in the calculation from the wind field inversion program, and extract their echo signals, elevation angles, azimuth angles, and Doppler velocities;

[0095] The second main module is used to convert the echo signal of each meteor from the time domain spectrum to the frequency spectrum, and intercept the frequency matrix at the position where the meteor echo signal is located from the frequency spectrum;

[0096] The third main module is used to calculate the new range resolution after increasing the target multiple based on the original range resolution, calculate the covariance matrix of the meteors according to the frequency matrix, substitute the covariance matrix into the Capon algorithm to calculate the pseudo-spectrum at the new range resolution, obtain the exact group distance of each meteor according to the pseudo-spectrum, and calculate the exact height of each meteor according to the obtained exact group distance;

[0097] The fourth main module is used to substitute the elevation angle, azimuth angle, Doppler velocity, and exact height of each meteor into the wind field inversion program to recalculate the wind field, and obtain the refined wind field at the new range resolution. The second main module converts the meteor echo signal into a frequency spectrum and intercepts the frequency matrix, then the third main module extracts the exact group distance and exact height of the signal and inputs them into the fourth main module, thereby quickly calculating the refined middle and upper atmosphere wind field.

[0098] In the application of this embodiment, the echo signals of the meteors participating in the wind field inversion are screened out from the first main module in the meteor echo signal data within a day, and then the echo signals are transmitted to the second main module. The second main module converts the echo signals from the time domain spectrum to the frequency spectrum, and extracts the frequency matrix at the position where the meteor echo signals are located from the frequency spectrum. Then, in the third main module, a distance resolution enhancement multiple, such as 4 times, is set to calculate the new distance resolution. Next, the covariance matrix of the meteors is calculated based on the frequency matrix, and then the covariance matrix is substituted into the Capon algorithm to calculate the pseudo-spectrum at the new distance resolution, and the precise group distance of each meteor is obtained based on the pseudo-spectrum. The precise altitude of each meteor is calculated based on the obtained precise group distance. Finally, the precise altitude of each meteor calculated from the obtained precise group distance is substituted into the fourth main module, and the elevation angle, azimuth angle, and Doppler velocity of each meteor calculated in the wind field inversion program are substituted, so as to obtain a refined wind field with a higher distance resolution.

[0099] Based on the same inventive concept as the foregoing embodiment, an embodiment of the present invention further provides an electronic device, including a memory and a processor. The memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions to implement a method for obtaining a fine structure by enhancing the altitude resolution of the middle and upper atmosphere wind field as proposed in the above embodiment.

[0100] An embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored. When the program is executed by a processor, it is used to more intuitively observe the results in computer tomography, expanding the application range of obtaining a fine structure by enhancing the altitude resolution of the middle and upper atmosphere wind field. The storage medium can be any non-volatile storage device such as a hard disk, a solid-state drive, a flash drive, an optical disc, etc., for storing computer program codes and necessary data files. The stored computer program includes: a data acquisition module, a photovoltaic panel identification module.

[0101] The above is only a preferred embodiment of the present invention, and it does not limit the embodiments and protection scope of the present invention. For those skilled in the art, it should be realized that all equivalent substitutions and obvious changes made by using the content of the specification of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for obtaining fine structures by enhancing the height resolution of the middle and upper atmosphere wind field, characterized in that, Including: Screening out the meteors participating in the wind field inversion program, and extracting the echo signal, elevation angle, azimuth angle, and Doppler velocity of each screened meteor; Converting the echo signal of each meteor from the time domain spectrum to the frequency spectrum, and intercepting the frequency matrix at the position where the meteor echo signal is located from the frequency spectrum; Calculating the new range resolution after increasing the target multiple according to the original range resolution, calculating the covariance matrix of the meteor according to the frequency matrix, substituting the covariance matrix into the Capon algorithm to calculate the pseudo-spectrum at the new range resolution, obtaining the precise group distance of each meteor according to the pseudo-spectrum, and calculating the precise height of each meteor according to the obtained precise group distance; Substituting the elevation angle, azimuth angle, Doppler velocity, and precise height of each meteor into the wind field inversion program to recalculate the wind field, and obtaining the refined wind field at the new range resolution.

2. The method for obtaining fine structure by enhancing the height resolution of the middle and upper atmosphere wind field according to claim 1, wherein: The screening out of the meteors participating in the wind field inversion program and extracting the echo signal of each meteor includes: tracing back the inversion process of the wind field, inversely calibrating the meteors participating in the wind field calculation according to the corresponding height and time information, and screening out the time domain signal of the meteor echo.

3. The method for obtaining fine structure by enhancing height resolution of the middle and upper atmosphere wind field according to claim 1, characterized in that: Intercepting the frequency matrix at the position where the meteor echo signal is located from the frequency spectrum includes: preprocessing the time domain spectrum of the meteor echo signal, initially determining the position of the signal, and intercepting a new frequency matrix.

4. The method for obtaining the fine structure by enhancing the height resolution of the middle and upper atmosphere wind field according to claim 3, characterized in that: Calculating the covariance matrix of the meteor according to the frequency matrix includes: Calculating the duration of each chip of the meteor radar modulation coding through the original range resolution inherent in the meteor radar detection device, and calculating the signal bandwidth through the chip duration; Splitting according to the distance dimension based on the new frequency matrix, and listing the frequency spectrum sequences under each range gate as sub-frequency matrices; Calculating the frequency component interval of the sub-frequency matrix through the signal bandwidth and the length of the original frequency spectrum dimension; Starting from the first frequency component interval, dividing and reconstructing the frequency spectrum to obtain the frequency matrix; Obtaining the covariance matrix of each frequency component through the reconstructed frequency matrix.

5. A method for obtaining fine structures by enhancing the height resolution of the middle and upper atmospheric wind field according to claim 4, characterized in that: Substituting the covariance matrix into the Capon algorithm to calculate the pseudo-spectrum at the new range resolution, and obtaining the precise group distance of each meteor according to the pseudo-spectrum includes: Setting the target multiple for increasing the range resolution and then calculating the increased new range resolution; Constructing the steering matrix of the frequency matrix according to the increased target multiple; Substituting the covariance matrix and the steering matrix into the Capon algorithm to obtain the pseudo-spectrum at the new range resolution; Analyzing the pseudo-spectrum along the distance dimension to obtain the precise group distance of the meteor echo signal.

6. The method for obtaining fine structure by enhancing height resolution of the middle and upper atmosphere wind field according to claim 5, characterized in that: Calculating the precise height of the meteor according to the precise group distance includes: calculating the precise height of the meteor according to the precise group distance and the elevation angle of the meteor obtained from the initial wind field inversion program.

7. A system for obtaining fine structures by enhancing the height resolution of the middle and upper atmosphere wind field, characterized in that, Including: The first main module is used to screen out the meteors participating in the calculation from the wind field inversion program, and extract their echo signals, elevation angles, azimuth angles, and Doppler velocities; The second main module is used to convert the echo signal of each meteor from the time domain spectrum to the frequency spectrum, and intercept the frequency matrix at the position where the meteor echo signal is located from the frequency spectrum; The third main module is used to calculate the new range resolution after increasing the target multiple according to the original range resolution, calculate the covariance matrix of meteors according to the frequency matrix, substitute the covariance matrix into the Capon algorithm to calculate the pseudospectrum at the new range resolution, obtain the precise group distance of each meteor according to the pseudospectrum, and calculate the precise height of each meteor according to the obtained precise group distance; The fourth main module is used to substitute the elevation angle, azimuth angle, Doppler velocity and precise height of each meteor into the wind field inversion program to recalculate the wind field, and obtain the refined wind field at the new range resolution; the meteor echo signal is converted into a frequency spectrum and the frequency matrix is intercepted through the second main module, and then the precise group distance and precise height of the signal are extracted by the third main module and input into the fourth main module, so as to quickly calculate the refined middle and upper atmosphere wind field.

8. An electronic device, comprising a memory and a processor, the memory storing a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method for obtaining the fine structure by improving the height resolution of the middle and upper atmosphere wind field according to any one of claims 1 to 6.

9. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the steps of the method for obtaining the fine structure by improving the height resolution of the middle and upper atmosphere wind field according to any one of claims 1 to 6.

Images