CN117741685A - Wind speed inversion method, system, storage medium and equipment - Google Patents
Wind speed inversion method, system, storage medium and equipment Download PDFInfo
- Publication number
- CN117741685A CN117741685A CN202410190281.5A CN202410190281A CN117741685A CN 117741685 A CN117741685 A CN 117741685A CN 202410190281 A CN202410190281 A CN 202410190281A CN 117741685 A CN117741685 A CN 117741685A
- Authority
- CN
- China
- Prior art keywords
- power spectrum
- inversion
- total power
- wind speed
- time window
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000003860 storage Methods 0.000 title claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims abstract description 141
- 239000000758 substrate Substances 0.000 claims abstract description 30
- 238000005070 sampling Methods 0.000 claims abstract description 18
- 230000006870 function Effects 0.000 claims description 30
- 230000008569 process Effects 0.000 claims description 10
- 238000013500 data storage Methods 0.000 claims description 7
- 230000015654 memory Effects 0.000 claims description 7
- 238000006467 substitution reaction Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 abstract description 20
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000002592 echocardiography Methods 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 description 1
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Abstract
The invention discloses a wind speed inversion method, a system, a storage medium and equipment, wherein laser radar echo is sampled in a single bit mode, single bit sampling data is converted into frequency spectrums by adopting fast Fourier transform based on a square wave substrate in a corresponding total power spectrum peak frequency range, wind speed inversion is carried out, and a frequency interval calculated by a back inversion period is narrowed by the total power spectrum peak frequency of a front inversion period, so that the calculated amount required by wind speed inversion can be greatly reduced, and the consumption of calculation resources is reduced.
Description
Technical Field
The invention relates to a wind speed inversion method, a system, a storage medium and equipment, and belongs to the field of Doppler wind lidar.
Background
The accurate observation of wind has great application value in the fields of environment, military, aerospace, weather, and the like, the Doppler wind-measuring laser radar has the advantage that other wind-measuring means are difficult to compare in measurement range and precision, and the Doppler effect is the basic principle of measuring wind speed. The laser pulse emitted by the laser radar acts with particles in the atmosphere to generate backward scattering, the velocity component of the particles in the laser direction enables the backward scattering echo to generate Doppler frequency shift, and the wind speed can be obtained through inversion by extracting Doppler frequency shift information.
The extraction of effective Doppler frequency spectrum information is of great significance to the realization of wind speed inversion, and a laser radar generally converts a time domain signal into a frequency domain signal, so that Doppler frequency shift information is obtained for subsequent analysis and processing. The traditional method of time-frequency conversion is to perform fast fourier transform on the signal, and although the fast fourier transform has less calculation amount compared with the discrete fourier transform, the calculation amount is still quite huge, and a large amount of calculation resources are required to be consumed.
Disclosure of Invention
The invention provides a wind speed inversion method, a wind speed inversion system, a storage medium and wind speed inversion equipment, and solves the problems disclosed in the background art.
In order to solve the technical problems, the invention adopts the following technical scheme:
a method of wind speed inversion, comprising:
for T j,m,i The echo in the device is sampled by a single bit; wherein T is j,m,i An ith echo time window for an mth laser in the jth inversion period;
at the estimated T j,i In the corresponding total power spectrum peak frequency range, adopting the fast Fourier transform based on the square wave substrate to convert T j,m,i The single bit sampled data in the data storage unit is converted into a frequency spectrum; wherein T is j,i An ith echo time window, T, for the jth inversion period j,i The corresponding peak frequency range of the power spectrum is according to T j-1,i And T j-2,i Estimating the peak frequency of the total power spectrum in the reactor; t (T) j-1,i And T j-2,i The total power spectrum of the ith echo time window in the inversion period is accumulated in the ith echo time window of all lasers in the inversion period;
will T j,m,i The frequency spectrum within the spectrum is converted to a power spectrum,t is carried out according to the power spectrum in the ith echo time window of each laser in the jth inversion period j,i And (5) inverting the internal wind speed.
Converting 0's in the single bit sampled data to-1's is also included before performing the fast fourier transform based on the square wave substrate.
The fast fourier transform based on square wave substrates is:
adopting the function sq_exp%jx) Substitution function exp # -jx) Performing fast Fourier transform;
wherein, sq_exp # (jx)= sq_cos(x)+ jsq_sin(x) ；
Function of；
Function of；
In the method, in the process of the invention,xis an independent variable.
T is carried out according to the power spectrum in the ith echo time window of each laser in the jth inversion period j,i An internal wind speed inversion comprising:
obtaining T according to the power spectrum in the ith echo time window of each laser in the jth inversion period j,i A total power spectrum within;
according to T j,i The total power spectrum in the wind speed range is searched for T j,i Total power spectrum peak frequency in;
according to T j,i Calculating the frequency shift of the power spectrum according to the peak frequency of the total power spectrum in the power spectrum;
calculating T according to the power spectrum frequency shift and the laser wavelength j,i A velocity component of the wind velocity in the laser direction.
T j,i The peak frequency range of the corresponding total power spectrum isf j-1,i -df，f j-1,i +df]The method comprises the steps of carrying out a first treatment on the surface of the Wherein,f j-1,i is T j-1,i The peak frequency of the total power spectrum within,dfin order to adjust the parameters of the device,df=max{df min ，q|f j-1,i -f j-2,i |}，qis a constant value, and is used for the treatment of the skin,df min is thatdfIs set to be a minimum value of (c),f j-2,i is T j-2,i Peak frequencies of the total power spectrum within.
A wind speed inversion system, comprising:
sampling module for T j,m,i The echo in the device is sampled by a single bit; wherein T is j,m,i An ith echo time window for an mth laser in the jth inversion period;
fast Fourier transform module based on square wave substrate, in estimated T j,i In the corresponding total power spectrum peak frequency range, adopting the fast Fourier transform based on the square wave substrate to convert T j,m,i The single bit sampled data in the data storage unit is converted into a frequency spectrum; wherein T is j,i An ith echo time window, T, for the jth inversion period j,i The corresponding peak frequency range of the power spectrum is according to T j-1,i And T j-2,i Estimating the peak frequency of the total power spectrum in the reactor; t (T) j-1,i And T j-2,i The total power spectrum of the ith echo time window in the inversion period is accumulated in the ith echo time window of all lasers in the inversion period;
inversion module, T is j,m,i The frequency spectrum in the inversion period is converted into a power spectrum, and T is carried out according to the power spectrum in the ith echo time window of each laser in the jth inversion period j,i And (5) inverting the internal wind speed.
In the fast fourier transform module based on the square wave substrate, the fast fourier transform based on the square wave substrate is as follows:
adopting the function sq_exp%jx) Substitution function exp # -jx) Performing fast Fourier transform;
wherein, sq_exp # (jx)= sq_cos(x)+ jsq_sin(x) ；
Function of；
Function of；
In the method, in the process of the invention,xis an independent variable.
In the fast Fourier transform module based on square wave substrate, T j,i The peak frequency range of the corresponding total power spectrum isf j-1,i -df，f j-1,i +df]The method comprises the steps of carrying out a first treatment on the surface of the Wherein,f j-1,i is T j-1,i The peak frequency of the total power spectrum within,dfin order to adjust the parameters of the device,df=max{df min ，q|f j-1,i -f j-2,i |}，qis a constant value, and is used for the treatment of the skin,df min is thatdfIs set to be a minimum value of (c),f j-2,i is T j-2,i Peak frequencies of the total power spectrum within.
A computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform a wind speed inversion method.
A computer device comprising one or more processors, and one or more memories in which one or more programs are stored and configured to be executed by the one or more processors, the one or more programs comprising instructions for performing a wind speed inversion method.
The invention has the beneficial effects that: according to the method, the laser radar echo is subjected to single-bit sampling, the single-bit sampling data is converted into the frequency spectrum by adopting the fast Fourier transform based on the square wave substrate in the range of the peak frequency of the corresponding total power spectrum, so that wind speed inversion is performed, the frequency interval calculated by the inversion period after the inversion period is narrowed by the total power spectrum peak frequency of the inversion period before the wind speed inversion is performed, the calculated amount required by wind speed inversion can be greatly reduced, and the consumption of calculation resources is reduced.
Drawings
FIG. 1 is a flow chart of a wind speed inversion method.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for more clearly illustrating the technical aspects of the present invention, and are not intended to limit the scope of the present invention.
As shown in fig. 1, a wind speed inversion method includes the following steps:
step 1, for T j,m,i The echo in the device is sampled by a single bit; wherein T is j,m,i An ith echo time window for the mth laser in the jth inversion period.
Step 2, in the estimated T j,i In the corresponding total power spectrum peak frequency range, adopting the fast Fourier transform based on the square wave substrate to convert T j,m,i The single bit sampled data in the data storage unit is converted into a frequency spectrum; wherein T is j,i An ith echo time window, T, for the jth inversion period j,i The corresponding peak frequency range of the power spectrum is according to T j-1,i And T j-2,i Estimating the peak frequency of the total power spectrum in the reactor; t (T) j-1,i And T j-2,i The total power spectrum of the ith echo time window in the inversion period is accumulated in the ith echo time window of all lasers in the inversion period.
Step 3, T is taken as j,m,i The frequency spectrum in the inversion period is converted into a power spectrum, and T is carried out according to the power spectrum in the ith echo time window of each laser in the jth inversion period j,i And (5) inverting the internal wind speed.
According to the method, the laser radar echo is subjected to single-bit sampling, the single-bit sampling data is converted into the frequency spectrum by adopting the fast Fourier transform based on the square wave substrate in the range of the peak frequency of the corresponding total power spectrum, so that wind speed inversion is performed, the frequency interval calculated by the inversion period after the inversion period is narrowed by the total power spectrum peak frequency of the inversion period before the inversion period after the inversion period is shortened, the calculated amount required by wind speed inversion is greatly reduced, and the consumption of calculation resources is reduced.
After the Doppler wind-measuring laser radar emits a beam of laser, the laser acts with particles at different distances to generate echoes, and the farther the distance is, the later the time for the echoes to reach the laser radar is, so that the time for receiving the echoes after the beam of laser is emitted is divided into a plurality of echo time windows. And (3) transmitting a plurality of lasers in one inversion period, wherein the number of echo time windows of each laser is consistent, the time spans of the time windows are consistent, and carrying out wind speed inversion according to echoes in the echo time windows of different lasers.
In each echo time window, during inversion, the sampled data is required to be converted into a frequency spectrum, then the frequency spectrum is further converted into a power spectrum, and the total power spectrum with enough signal-to-noise ratio is obtained in an accumulation mode, so that the wind speed is inverted.
Conventional data sampling uses high precision sampling and employs a multi-bit sampling circuit, which can replace data with conventional high precision sampled data (typically 8-bit data), i.e., for T, with single-bit data in order to simplify subsequent computation and circuitry j,m,i The echo in the signal is sampled by a single bit, T j,m,i For the ith echo time window of the mth laser in the jth inversion period, the sampling circuit of the laser radar echo can be simplified into a single-bit sampling circuit by a multi-bit sampling circuit.
Since the single bit data has a value of 0 or 1, its mean value is 0.5, which may cause the low frequency part of the calculated spectrum to be affected. Thus, 0 in the single bit sampled data may be converted to-1 prior to spectral conversion, and then subsequent calculations may be performed with a data sequence consisting of 1 and-1.
Assuming that the laser radar echo corresponding to a certain echo time window is sampled and is subjected to proper preprocessing to obtainx(n)，n=0,1,…,N-1，N=2 L ，LIs a positive integer which is used for the preparation of the high-voltage power supply,x(n) Is the firstnThe number of pre-processed sample data,Nis the total number of the sample data after preprocessing in the echo time window, if the discrete Fourier transform pair is adoptedx(n) Performing time-frequency conversion to obtain frequency spectrum dataX(k) The method comprises the following steps:
；
in the method, in the process of the invention,k=0,1,…,N-1，X(k) Is plural, willX(k) Halving and obtaining a power spectrumP(k)。
Since the discrete Fourier transform has large calculation amount, the fast Fourier transform can be adopted for replacement, thereby saving a large amount of operation, but the fast Fourier transform still needs to be usedSuch high precision complex parameters still result in a huge computational effort, and therefore according to the euler formula:
exp(jx)= cos(x)+ jsq_sin(x)
due toTo simplify the calculation, let:
square wave cosine function；
Square wave sine function；
Wherein, the function sq_exp%jx)= sq_cos(x)+ jsq_sin(x)，xIs an independent variable.
Therefore, the invention adopts the function sq_exp%jx) Substitution of exponential function expjx) Performing a fast Fourier transform, i.e. single bit sampled data conversion using a fast Fourier transform based on a square wave substrate, suchThe term is replaced by 1 or-1, which can greatly reduce the computational complexity.
The fast Fourier transform based on the square wave substrate obtains that each frequency spectrum is a complex number, and the real part and the imaginary part of the complex number are calculated by square sum, so that the frequency spectrum can be converted into a power spectrum. Because the echo signal of a single laser pulse is weak, for each inversion period, the power spectrums obtained by a plurality of laser pulses need to be accumulated to obtain a total power spectrum with enough signal-to-noise ratio, and particularly, the ith echo time window of all the beams of laser in the jth inversion period can be used forAccumulating the power spectrums to obtain T j,i A total power spectrum within.
After the accumulation is finished, according to T j,i The total power spectrum in the wind speed range is searched for T j,i Peak frequency of total power spectrum in the spectrum according to T j,i The peak frequency of the total power spectrum in the wind speed is calculated, the power spectrum frequency shift corresponding to the wind speed is calculated, and T is calculated according to the power spectrum frequency shift and the laser wavelength j,i The velocity component of the wind velocity in the laser direction can be expressed as follows:
v=λ△f/2
in the method, in the process of the invention,vas a velocity component of the wind velocity in the laser direction,λthe wind speed range can be-30 m/s to +30 m/s for laser wavelength.
In order to further reduce the calculation amount required by the calculation of the power spectrum, after the total power spectrum with distance resolution is obtained through accumulation, the method of the invention estimates the possible range of the peak frequency of the total power spectrum in the same time window of the next period according to the peak frequency of the total power spectrum in the same time window of the last two inversion periods besides the calculation of the wind speed, so that in the following spectrum calculation, the calculated frequency interval is narrowed, i.e. the calculation is not performed any moreX(0) To the point ofX(N-1) all the calculations are performed once, but only the range covered part is calculated, so that unnecessary spectrum calculation is avoided and the calculation amount required for time-frequency conversion of the laser echo signal can be greatly reduced.
The process of determining the possible range of the estimated power spectrum peak frequency may be as follows:
definition T j,i An ith echo time window, T, for the jth inversion period j-1,i And T j-2,i An ith echo time window, T, for the jth-1 inversion period and the jth-2 inversion period, respectively j-1,i The peak frequency of the total power spectrum in the internal isf j-1,i ，T j-2,i The peak frequency of the total power spectrum in the internal isf j-2,i Then T j,i The peak frequency range of the corresponding total power spectrum isf j-1,i -df，f j-1,i +df]In the case of calculation within the range onlyX(k) Outside the range does not needCalculation, by which the frequency range of the spectrum calculation can be adaptively adjusted; wherein,dffor adjusting parameters, adjusting according to the wind measurement condition to ensure that the calculated amount is reduced and the peak frequency of the power spectrum is not covered;df=max{df min ，q|f j-1,i -f j-2,i |}，qis constant, typically between 1 and 3,df min is thatdfIs a minimum of (2).
Through the analysis, T is adopted j,i For example, the specific process of the wind speed inversion method may be as follows:
1) For T j,m,i The laser radar echo in the laser radar echo detection device performs single-bit sampling, and 0 in the single-bit sampling data is converted into-1.
2) If the current inversion period j is the first period or the second period, adopting the fast Fourier transform based on the square wave substrate to convert T j,m,i The single bit sampled data in the data storage unit is converted into a frequency spectrum;
if the current inversion period j is not the first period or the second period, according to T j-1,i And T j-2,i Total power spectrum peak frequency estimate T within j,i Corresponding power spectrum peak frequency range, at T j,i In the corresponding power spectrum peak frequency range, adopting the fast Fourier transform based on the square wave substrate to convert T j,m,i The single bit sampled data within is converted to a spectrum.
3) Will T j,m,i The frequency spectrum in the inversion period j is converted into a power spectrum, and the power spectrums in the ith echo time window of all lasers in the inversion period j are accumulated to obtain T j,i Total power spectrum in according to T j,i The total power spectrum in the wind speed range is searched for T j,i Peak frequencies of the total power spectrum within.
4) According to T j,i The peak frequency of the total power spectrum in the laser is calculated, the power spectrum frequency shift is calculated, and T is calculated according to the power spectrum frequency shift and the laser wavelength j,i A velocity component of the wind velocity in the laser direction.
Based on the same technical scheme, the invention also discloses a software system of the method, and a wind speed inversion system comprises:
sampling module for T j,m,i The echo in the device is sampled by a single bit; wherein T is j,m,i An ith echo time window for the mth laser in the jth inversion period.
Fast Fourier transform module based on square wave substrate, in estimated T j,i In the corresponding total power spectrum peak frequency range, adopting the fast Fourier transform based on the square wave substrate to convert T j,m,i The single bit sampled data in the data storage unit is converted into a frequency spectrum; wherein T is j,i An ith echo time window, T, for the jth inversion period j,i The corresponding peak frequency range of the power spectrum is according to T j-1,i And T j-2,i Estimating the peak frequency of the total power spectrum in the reactor; t (T) j-1,i And T j-2,i The total power spectrum of the ith echo time window in the inversion period is accumulated in the ith echo time window of all lasers in the inversion period.
In the fast fourier transform module based on the square wave substrate, the fast fourier transform based on the square wave substrate is as follows:
adopting the function sq_exp%jx) Function substitution function exp @jx) Performing fast Fourier transform;
wherein, sq_exp # (jx)= sq_cos(x)+ jsq_sin(x) ；
Function of；
Function of；
In the method, in the process of the invention,xis an independent variable.
In the fast Fourier transform module based on square wave substrate, T j,i The peak frequency range of the corresponding total power spectrum isf j-1,i -df，f j-1,i +df]The method comprises the steps of carrying out a first treatment on the surface of the Wherein,f j-1,i is T j-1,i The peak frequency of the total power spectrum within,dfin order to adjust the parameters of the device,df=max{df min ，q|f j-1,i -f j-2,i |}，qis a constant value, and is used for the treatment of the skin,df min is thatdfIs set to be a minimum value of (c),f j-2,i is T j-2,i Peak frequencies of the total power spectrum within.
Inversion module, T is j,m,i The frequency spectrum in the inversion period is converted into a power spectrum, and T is carried out according to the power spectrum in the ith echo time window of each laser in the jth inversion period j,i And (5) inverting the internal wind speed.
According to the system, the laser radar echo is subjected to single-bit sampling, the single-bit sampling data is converted into the frequency spectrum by adopting the fast Fourier transform based on the square wave substrate in the range of the peak frequency of the corresponding total power spectrum, so that wind speed inversion is performed, the frequency interval calculated by the inversion period after the inversion period is narrowed by the total power spectrum peak frequency of the inversion period before the inversion period after the inversion period is shortened, the calculated amount required by wind speed inversion is greatly reduced, and the consumption of calculation resources is reduced.
Based on the same technical solution, the present invention also discloses a computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform a wind speed inversion method.
Based on the same technical scheme, the invention also discloses a computer device, which comprises one or more processors and one or more memories, wherein one or more programs are stored in the one or more memories and are configured to be executed by the one or more processors, and the one or more programs comprise instructions for executing the wind speed inversion method.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing is illustrative of the present invention and is not to be construed as limiting thereof, but rather as providing for the use of additional embodiments and advantages of all such modifications, equivalents, improvements and similar to the present invention are intended to be included within the scope of the present invention as defined by the appended claims.
Claims (10)
1. A method of wind speed inversion comprising:
for T j,m,i The echo in the device is sampled by a single bit; wherein T is j,m,i An ith echo time window for an mth laser in the jth inversion period;
at the estimated T j,i In the corresponding total power spectrum peak frequency range, adopting the fast Fourier transform based on the square wave substrate to convert T j,m,i The single bit sampled data in the data storage unit is converted into a frequency spectrum; wherein T is j,i An ith echo time window, T, for the jth inversion period j,i The corresponding peak frequency range of the power spectrum is according to T j-1,i And T j-2,i Estimating the peak frequency of the total power spectrum in the reactor; t (T) j-1,i And T j-2,i The total power spectrum of the ith echo time window in the inversion period is accumulated in the ith echo time window of all lasers in the inversion period;
will T j,m,i The frequency spectrum in the inversion period is converted into a power spectrum, and T is carried out according to the power spectrum in the ith echo time window of each laser in the jth inversion period j,i And (5) inverting the internal wind speed.
2. The wind speed inversion method of claim 1 further comprising converting 0's in the single bit sampled data to-1's prior to performing the fast fourier transform based on the square wave substrate.
3. The wind speed inversion method of claim 1 wherein the fast fourier transform based on square wave substrate is:
adopting the function sq_exp%jx) Substitution function exp # -jx) Performing fast Fourier transform;
wherein, sq_exp # (jx)= sq_cos(x)+ jsq_sin(x)；
Function of；
Function of；
In the method, in the process of the invention,xis an independent variable.
4. The method according to claim 1, wherein T is performed based on a power spectrum in an ith echo time window of each laser in a jth inversion period j,i An internal wind speed inversion comprising:
obtaining T according to the power spectrum in the ith echo time window of each laser in the jth inversion period j,i A total power spectrum within;
according to T j,i The total power spectrum in the wind speed range is searched for T j,i Total power spectrum peak frequency in;
according to T j,i Calculating the frequency shift of the power spectrum according to the peak frequency of the total power spectrum in the power spectrum;
calculating T according to the power spectrum frequency shift and the laser wavelength j,i A velocity component of the wind velocity in the laser direction.
5. The method of wind speed inversion according to claim 1, wherein T j,i The peak frequency range of the corresponding total power spectrum isf j-1,i -df，f j-1,i +df]The method comprises the steps of carrying out a first treatment on the surface of the Wherein,f j-1,i is T j-1,i The peak frequency of the total power spectrum within,dfin order to adjust the parameters of the device,df=max{df min ，q|f j-1,i - f j-2,i |}，qis a constant value, and is used for the treatment of the skin,df min is thatdfIs set to be a minimum value of (c),f j-2,i is T j-2,i Peak frequencies of the total power spectrum within.
6. A wind speed inversion system, comprising:
sampling module for T j,m,i The echo in the device is sampled by a single bit; wherein T is j,m,i An ith echo time window for an mth laser in the jth inversion period;
fast Fourier transform module based on square wave substrate, in estimated T j,i In the corresponding total power spectrum peak frequency range, adopting the fast Fourier transform based on the square wave substrate to convert T j,m,i The single bit sampled data in the data storage unit is converted into a frequency spectrum; wherein T is j,i An ith echo time window, T, for the jth inversion period j,i The corresponding peak frequency range of the power spectrum is according to T j-1,i And T j-2,i Estimating the peak frequency of the total power spectrum in the reactor; t (T) j-1,i And T j-2,i The total power spectrum of the ith echo time window in the inversion period is accumulated in the ith echo time window of all lasers in the inversion period;
inversion module, T is j,m,i The frequency spectrum in the inversion period is converted into a power spectrum, and T is carried out according to the power spectrum in the ith echo time window of each laser in the jth inversion period j,i And (5) inverting the internal wind speed.
7. The wind speed inversion system of claim 6 wherein in the fast fourier transform based on square wave substrate module, the fast fourier transform based on square wave substrate is:
adopting the function sq_exp%jx) Substitution function exp # -jx) Performing fast Fourier transform;
wherein, sq_exp # (jx)= sq_cos(x)+ jsq_sin(x)；
Function of；
Function of；
In the method, in the process of the invention,xis an independent variable.
8. The wind speed inversion system of claim 6 wherein in the fast fourier transform module based on square wave substrates, T j,i The peak frequency range of the corresponding total power spectrum isf j-1,i -df，f j-1,i +df]The method comprises the steps of carrying out a first treatment on the surface of the Wherein,f j-1,i is T j-1,i The peak frequency of the total power spectrum within,dfin order to adjust the parameters of the device,df=max{df min ，q|f j-1,i - f j-2,i |}，qis a constant value, and is used for the treatment of the skin,df min is thatdfIs set to be a minimum value of (c),f j-2,i is T j-2,i Peak frequencies of the total power spectrum within.
9. A computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform any of the methods of claims 1-5.
10. A computer device, comprising:
one or more processors, and one or more memories, one or more programs stored in the one or more memories and configured to be executed by the one or more processors, the one or more programs comprising instructions for performing any of the methods of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410190281.5A CN117741685A (en) | 2024-02-21 | 2024-02-21 | Wind speed inversion method, system, storage medium and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410190281.5A CN117741685A (en) | 2024-02-21 | 2024-02-21 | Wind speed inversion method, system, storage medium and equipment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117741685A true CN117741685A (en) | 2024-03-22 |
Family
ID=90251244
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410190281.5A Pending CN117741685A (en) | 2024-02-21 | 2024-02-21 | Wind speed inversion method, system, storage medium and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117741685A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105022036A (en) * | 2015-08-26 | 2015-11-04 | 成都信息工程大学 | Wind profile radar wind speed determination method |
| CN107390229A (en) * | 2017-06-14 | 2017-11-24 | 南京牧镭激光科技有限公司 | A kind of processing method and its processing unit of anemometry laser radar signal |
| US20200064485A1 (en) * | 2016-11-28 | 2020-02-27 | Mitsubishi Electric Corporation | Data processing device, laser radar device, and wind measurement system |
| CN112505719A (en) * | 2020-10-15 | 2021-03-16 | 中国科学院西安光学精密机械研究所 | Laser wind measuring radar wind measuring method and system with Doppler frequency correction secondary compensation |
| CN114355387A (en) * | 2021-11-26 | 2022-04-15 | 山东师范大学 | Wind field inversion method and system based on eight-beam wind profile laser radar |
| CN117076833A (en) * | 2023-10-16 | 2023-11-17 | 深圳大学 | FPGA-based single-bit dimension-reduction fast Fourier transform method |
-
2024
- 2024-02-21 CN CN202410190281.5A patent/CN117741685A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105022036A (en) * | 2015-08-26 | 2015-11-04 | 成都信息工程大学 | Wind profile radar wind speed determination method |
| US20200064485A1 (en) * | 2016-11-28 | 2020-02-27 | Mitsubishi Electric Corporation | Data processing device, laser radar device, and wind measurement system |
| CN107390229A (en) * | 2017-06-14 | 2017-11-24 | 南京牧镭激光科技有限公司 | A kind of processing method and its processing unit of anemometry laser radar signal |
| CN112505719A (en) * | 2020-10-15 | 2021-03-16 | 中国科学院西安光学精密机械研究所 | Laser wind measuring radar wind measuring method and system with Doppler frequency correction secondary compensation |
| CN114355387A (en) * | 2021-11-26 | 2022-04-15 | 山东师范大学 | Wind field inversion method and system based on eight-beam wind profile laser radar |
| CN117076833A (en) * | 2023-10-16 | 2023-11-17 | 深圳大学 | FPGA-based single-bit dimension-reduction fast Fourier transform method |
Non-Patent Citations (1)
| Title |
|---|
| 李光福: "基于FPGA的激光雷达测风系统高速信号采集处理研究", 《中国优秀硕士学位论文全文数据库基础科学辑》, 15 February 2023 (2023-02-15), pages 5 - 18 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106772352B (en) | It is a kind of that Weak target detecting method is extended based on the PD radar of Hough and particle filter | |
| CN110673108B (en) | Airborne marine laser radar signal processing method based on iteration Klett | |
| CN109581275B (en) | Two-dimensional underwater DOA estimation method and device based on non-circular signal and three-dimensional orthogonal array | |
| CN111610503B (en) | Linear frequency modulation signal parameter estimation method based on improved LVD | |
| Miao et al. | A new parallel implementation for particle filters and its application to adaptive waveform design | |
| CN113640790A (en) | Wide-bandwidth pulse high-speed target detection method based on two-dimensional adaptive spectrum estimation | |
| CN103364769A (en) | Adaptive estimation method for spatially non-uniform sea clutter covariance matrix | |
| US11921242B2 (en) | Method and device for time-of-flight estimate | |
| Tian et al. | Underwater multi-target passive detection based on transient signals using adaptive empirical mode decomposition | |
| Tran et al. | Microwave radar imaging of rotating blades | |
| CN112285742B (en) | Method for estimating carrier-to-noise ratio of coherent wind lidar in frequency domain | |
| CN112328965B (en) | Method for DOA tracking of multiple mechanical signal sources using acoustic vector sensor array | |
| CN117741685A (en) | Wind speed inversion method, system, storage medium and equipment | |
| CN110658517B (en) | Dimensionality reduction sparse STAP method and device based on uncertain priori knowledge | |
| CN110850421A (en) | Underwater target detection method based on space-time adaptive processing of reverberation symmetric spectrum | |
| CN1609570A (en) | Weak signal detection under complex background and characteristic analysis system | |
| Xu et al. | Parameter estimation of underwater moving sources by using matched Wigner transform | |
| CN115345216A (en) | FMCW radar interference elimination method fusing prior information | |
| CN109581363B (en) | Small-size space debris detection and parameter extraction method based on incoherent scattering radar | |
| CN111736146A (en) | Bistatic pre-detection tracking method and device based on speed filtering | |
| CN116908834B (en) | Multi-target positioning tracking method and system based on millimeter wave radar | |
| Shen et al. | Wind turbine clutter mitigation for weather radar by an improved low-rank matrix recovery method | |
| Ju et al. | A full parameter estimation method of LFM pulse signal | |
| CN117309079B (en) | Ultrasonic flying time measuring method, device, equipment and medium based on time difference method | |
| Van der Merwe et al. | Realtime generation of K-Distributed sea clutter for hardware in the loop radar evaluation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |