CN102540156B - Wind profile radar signal and data processing method based on two channels - Google Patents
Wind profile radar signal and data processing method based on two channels Download PDFInfo
- Publication number
- CN102540156B CN102540156B CN201010586765.XA CN201010586765A CN102540156B CN 102540156 B CN102540156 B CN 102540156B CN 201010586765 A CN201010586765 A CN 201010586765A CN 102540156 B CN102540156 B CN 102540156B
- Authority
- CN
- China
- Prior art keywords
- groups
- signal
- sigma
- fft
- wind profile
- 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.)
- Active
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a wind profile radar signal and data processing method based on two channels. The wind profile radar signal and data processing method is implemented by a signal processor and a data terminal, and comprises the specific steps that: the signal processor sets two groups of processing parameters; 2, the signal processor generates two groups of power spectrum data; 3, the data terminal processes the two groups of power spectrum data respectively; and 4, the data terminal selects a processing result. Therefore, wind profile radar signal and data processing based on the two channels is implemented. The method has the advantages that: the data output rate of a wind profile radar system is not changed; the new increased processing load is relatively small; the method is easily implemented; and the extraction capacity and the estimation precision of a wind profile radar to weak signals are improved effectively.
Description
Technical field
The present invention relates to a kind of radar signal and data processing method, particularly a kind of based on twin-channel Wind profile radar signal and data processing method.
Background technology
Current Wind profile radar signal and data processing method realize based on signal processor and data terminal, and the key step that signal is processed is: coherent accumulation, and by N
cIthe individual echoed signal of in succession returning is carried out coherent accumulation, wherein, and N
cIbe called coherent accumulation number; Spectral transformation, by N
fFTechoed signal after individual coherent accumulation is carried out fast fourier transform (FFT), obtains a power spectral density function, N
fFTbeing called FFT counts; Spectrum is average, by N
sPCindividual power spectral density function averages, and obtains a new average power spectral density function (being designated hereinafter simply as power spectrum), N
sPCbe called spectral average.Afterwards, data terminal is carried out Spectrum data processing to power spectrum, obtains echo power Pr, Doppler shift f
d, speed spectrum width σ and signal to noise ratio snr.
In above signal processing and Spectrum data processing, coherent accumulation is the committed step that wind profile radar is realized weak signal extraction, and its major function is to improve signal to noise ratio (S/N ratio), in theory, and N
cIinferior coherent accumulation can improve signal to noise ratio (S/N ratio) N
cIdoubly.The limit that coherent accumulation improves signal to noise ratio (S/N ratio) depends primarily on the auto-correlation time of return signal in the phase stability of Wind Profiler Radar System itself and atmosphere, when system phase degree of stability meets the demands, the number of times that carries out coherent accumulation is main relevant with auto-correlation time and the recurrent interval of atmosphere return signal.The wavelength that in atmosphere, the auto-correlation time of return signal selects with atmospheric condition and wind profile radar is relevant, for troposphere wind profiler radar (pattern-band), the auto-correlation time selects 0.15~0.45 second, and the PBL wind profile radar (L-band) is selected 0.05~0.15 second.
Because atmospheric condition is ever-changing, when coherent accumulation time that system adopts and return signal the auto-correlation time, difference was larger time, theory and practice all shows, can affect the effect of coherent accumulation, is unfavorable for the extraction of feeble signal.Although in Wind Profiler Radar System design, coherent accumulation number has certain selectable range, for carrying out in actual detection preferably, but most wind profile radar all take one group fixedly processing parameter automatically move, though part often changes processing parameter for the wind profile radar of studying for research needs, and can not automatically adapt to different atmospheric conditions.
Summary of the invention
The object of the present invention is to provide a kind of based on twin-channel Wind profile radar signal and data processing method, while solving traditional Wind profile radar signal processing, the coherent accumulation time existing and return signal auto-correlation time can not Auto-matchings, are unfavorable for the problem of weak signal extraction.
Based on twin-channel Wind profile radar signal and a data processing method, based on signal processor and data terminal, realize, its concrete steps are:
First step signal processor arranges two groups of processing parameters
For the auto-correlation time, troposphere wind profiler radar is selected respectively 0.15~0.225,0.3~0.45 second, and the PBL wind profile radar is selected respectively 0.05~0.075,0.1~0.15 second.Auto-correlation time and wave beam residence time are respectively:
τ=PRT×N
CI×N
FFT (1)
T
BEAM=PRT×N
CI×N
FFT×N
SPC (2)
In formula, τ is the auto-correlation time, T
bEAMfor wave beam residence time, PRT is the pulse repetition time, N
cIfor coherent accumulation number, N
fFTfor FFT counts, N
sPCfor spectral average.
Determine after wave beam residence time, according to select two groups of auto-correlation times, keep pulse repetition time, FFT to count constant, calculate two groups of different processing parameters.Wherein, the coherent accumulation number of first group of parameter is two times of second group of parameter, and spectral average is 1/2nd of second group of parameter.
Second step signal processor generates two groups of power spectrum datas
Signal processor, to identical echoed signal, according to two groups of processing parameters, carries out respectively coherent accumulation, spectral transformation, spectrum average treatment successively, generates two groups of power spectrum datas.Two groups of power spectrum datas have different not fuzzy speed and velocity resolution, are respectively:
V
MAX=λ/(4×PRT×N
CI) (3)
V
MIN=λ/(2×PRT×N
CI×N
FFT) (4)
In formula, V
mAXfor not fuzzy speed, V
mINfor velocity resolution, λ is radar wavelength.
The 3rd step data terminal is carried out respectively the processing of two groups of power spectrum datas
Data terminal, to two groups of power spectrum datas, is used identical processing parameter to carry out respectively square estimation, obtains echo power, Doppler shift, speed spectrum width and signal to noise ratio (S/N ratio).Formula is:
In formula, Pr is echo power, f
dfor Doppler shift, σ speed spectrum width, SNR is signal to noise ratio (S/N ratio), P
nfor noise power, p
ibe i point performance number, f
iit is i point Doppler shift value.
The 4th step data terminal is carried out result and is chosen
Data terminal is the result of two groups of power spectrum datas relatively, according to the high person's priority principle of signal to noise ratio (S/N ratio), chooses wind spectrum discrimination and the estimated result of high s/n ratio.Formula is:
SNR1≥SNR2,Pr=Pr1,f
d=f
d1,σ=σ1,SNR=SNR1
SNR1<SNR2,Pr=Pr2,f
d=f
d2,σ=σ2,SNR=SNR2 (9)
So far, realized based on twin-channel Wind profile radar signal and data processing.
This method adopts two kinds of different coherent accumulations and incoherent accumulation number, solves coherent accumulation and does not mate with real atmosphere situation, causes signal and data processing to extract difficult problem to weak echo signal.The advantage of the method is not change Wind Profiler Radar System output data rate, newly-increased processing load is less, realization is simple, can effectively strengthen wind profile radar to the extractability of feeble signal and estimated accuracy.
Embodiment
Based on twin-channel Wind profile radar signal and a data processing method, based on signal processor and data terminal, realize, its concrete steps are:
First step signal processor arranges two groups of processing parameters
For the auto-correlation time, troposphere wind profiler radar is selected respectively 0.2,0.4 second, and the PBL wind profile radar is selected respectively 0.065,0.13 second.Auto-correlation time and wave beam residence time are respectively:
τ=PRT×N
CI×N
FFT (1)
T
BEAM=PRT×N
CI×N
FFT×N
SPC (2)
In formula, τ is the auto-correlation time, T
bEAMfor wave beam residence time, PRT is the pulse repetition time, N
cIfor coherent accumulation number, N
fFTfor FFT counts, N
sPCfor spectral average.
Determine after wave beam residence time, according to select two groups of auto-correlation times, keep pulse repetition time, FFT to count constant, calculate two groups of different processing parameters.Wherein, the coherent accumulation number of first group of parameter is two times of second group of parameter, and spectral average is 1/2nd of second group of parameter.
Second step signal processor generates two groups of power spectrum datas
Signal processor, to identical echoed signal, according to two groups of processing parameters, carries out respectively coherent accumulation, spectral transformation, spectrum average treatment successively, generates two groups of power spectrum datas.Two groups of power spectrum datas have different not fuzzy speed and velocity resolution, are respectively:
V
MAX=λ/(4×PRT×N
CI) (3)
V
MIN=λ/(2×PRT×N
CI×N
FFT) (4)
In formula, V
mAXfor not fuzzy speed, V
mINfor velocity resolution, λ is radar wavelength.
The 3rd step data terminal is carried out respectively the processing of two groups of power spectrum datas
Data terminal, to two groups of power spectrum datas, is used identical processing parameter to carry out respectively square estimation, obtains echo power, Doppler shift, speed spectrum width and signal to noise ratio (S/N ratio).Formula is:
In formula, Pr is echo power, f
dfor Doppler shift, σ speed spectrum width, SNR is signal to noise ratio (S/N ratio), P
nfor noise power, p
ibe i point performance number, f
iit is i point Doppler shift value.
The 4th step data terminal is carried out result and is chosen
Data terminal is the result of two groups of power spectrum datas relatively, according to the high person's priority principle of signal to noise ratio (S/N ratio), chooses wind spectrum discrimination and the estimated result of high s/n ratio.Formula is:
SNR1≥SNR2,Pr=Pr1,f
d=f
d1,σ=σ1,SNR=SNR1
SNR1<SNR2,Pr=Pr2,f
d=f
d2,σ=σ2,SNR=SNR2 (9)
So far, realized based on twin-channel Wind profile radar signal and data processing.
Claims (1)
1. based on twin-channel Wind profile radar signal and a data processing method, based on signal processor and data terminal, realize, it is characterized in that concrete steps are:
First step signal processor arranges two groups of processing parameters
For the auto-correlation time, troposphere wind profiler radar is selected respectively 0.15~0.225,0.3~0.45 second, and the PBL wind profile radar is selected respectively 0.05~0.075,0.1~0.15 second; Auto-correlation time and wave beam residence time are respectively:
τ=PRT×N
CI×N
FFT (1)
T
BEAM=PRT×N
CI×N
FFT×N
SPC (2)
In formula, τ is the auto-correlation time, T
bEAMfor wave beam residence time, PRT is the pulse repetition time, N
cIfor coherent accumulation number, N
fFTfor FFT counts, N
sPCfor spectral average;
Determine after wave beam residence time, according to select two groups of auto-correlation times, keep pulse repetition time, FFT to count constant, calculate two groups of different processing parameters; Wherein, the coherent accumulation number of first group of parameter is two times of second group of parameter, and spectral average is 1/2nd of second group of parameter;
Second step signal processor generates two groups of power spectrum datas
Signal processor, to identical echoed signal, according to two groups of processing parameters, carries out respectively coherent accumulation, spectral transformation, spectrum average treatment successively, generates two groups of power spectrum datas; Two groups of power spectrum datas have different not fuzzy speed and velocity resolution, are respectively:
V
MAX=λ/(4×PRT×N
CI) (3)
V
MIN=λ/(2×PRT×N
CI×N
FFT) (4)
In formula, V
mAXfor not fuzzy speed, V
mINfor velocity resolution, λ is radar wavelength;
The 3rd step data terminal is carried out respectively the processing of two groups of power spectrum datas
Data terminal, to two groups of power spectrum datas, is used identical processing parameter to carry out respectively square estimation, obtains echo power, Doppler shift, speed spectrum width and signal to noise ratio (S/N ratio); Formula is:
In formula, Pr is echo power, f
dfor Doppler shift, σ speed spectrum width, SNR is signal to noise ratio (S/N ratio), P
nfor noise power, p
ibe i point performance number, f
iit is i point Doppler shift value;
The 4th step data terminal is carried out result and is chosen
Data terminal is the result of two groups of power spectrum datas relatively, according to the high person's priority principle of signal to noise ratio (S/N ratio), chooses wind spectrum discrimination and the estimated result of high s/n ratio; Formula is:
If SNR1>=SNR2, shellfish Pr=Pr1, f
d=f
d1, σ=σ 1, SNR=SNR1
If SNR1 < is SNR2, shellfish Pr=Pr2, f
d=f
d2, σ=σ 2, SNR=SNR2 (9)
So far, realized based on twin-channel Wind profile radar signal and data processing.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010586765.XA CN102540156B (en) | 2010-12-14 | 2010-12-14 | Wind profile radar signal and data processing method based on two channels |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010586765.XA CN102540156B (en) | 2010-12-14 | 2010-12-14 | Wind profile radar signal and data processing method based on two channels |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102540156A CN102540156A (en) | 2012-07-04 |
CN102540156B true CN102540156B (en) | 2014-03-12 |
Family
ID=46347491
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010586765.XA Active CN102540156B (en) | 2010-12-14 | 2010-12-14 | Wind profile radar signal and data processing method based on two channels |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102540156B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10706056B1 (en) | 2015-12-02 | 2020-07-07 | Palantir Technologies Inc. | Audit log report generator |
CN106443678B (en) * | 2016-08-31 | 2018-10-30 | 王�华 | Utilize the atmospheric duct monitoring method of wind profile radar and RASS |
CN109471093A (en) * | 2018-11-07 | 2019-03-15 | 中国人民解放军国防科技大学 | Single pulse radar sum and difference correlation target detection method and system |
CN109614579B (en) * | 2018-12-24 | 2022-05-17 | 雷象科技(北京)有限公司 | Phased array weather radar maximum value signal extraction method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0436048A1 (en) * | 1990-01-02 | 1991-07-10 | Max-Planck-Gesellschaft zur Förderung der Wissenschaften e.V. | Oblique spaced antenna method and system for measuring atmospheric wind fields |
US5568151A (en) * | 1995-05-17 | 1996-10-22 | Merritt; David A. | Statistical averaging method for wind profiler doppler spectra |
CN101881824A (en) * | 2009-05-05 | 2010-11-10 | 何平 | Objective and fast determination method of noise threshold of power spectrum density data |
-
2010
- 2010-12-14 CN CN201010586765.XA patent/CN102540156B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0436048A1 (en) * | 1990-01-02 | 1991-07-10 | Max-Planck-Gesellschaft zur Förderung der Wissenschaften e.V. | Oblique spaced antenna method and system for measuring atmospheric wind fields |
US5568151A (en) * | 1995-05-17 | 1996-10-22 | Merritt; David A. | Statistical averaging method for wind profiler doppler spectra |
CN101881824A (en) * | 2009-05-05 | 2010-11-10 | 何平 | Objective and fast determination method of noise threshold of power spectrum density data |
Non-Patent Citations (4)
Title |
---|
SUSANNE DRECHSEL et al..Three-Dimensional Wind Retrieval: Application of MUSCAT to Dual-Doppler Lidar.《Journal of Atmospheric and Oceanic Technology》.2009,第26卷(第3期), |
Three-Dimensional Wind Retrieval: Application of MUSCAT to Dual-Doppler Lidar;SUSANNE DRECHSEL et al.;《Journal of Atmospheric and Oceanic Technology》;20090331;第26卷(第3期);全文 * |
基于滤波器组的风廓线雷达信号处理技术;李广柱 等;《雷达科学与技术》;20071231(第6期);全文 * |
李广柱 等.基于滤波器组的风廓线雷达信号处理技术.《雷达科学与技术》.2007,(第6期), |
Also Published As
Publication number | Publication date |
---|---|
CN102540156A (en) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102169175B (en) | Method for sea clutter suppression and moving target detection based on FRFT (Fractional Fourier Transform) domain non-coherent accumulation | |
CN108693527A (en) | A kind of sea surface small target acquisition radar system | |
CN102540156B (en) | Wind profile radar signal and data processing method based on two channels | |
RU2594005C1 (en) | Method of processing radar signal in pulse-doppler radar set | |
CN101334469B (en) | Wind profile radar clutter suppression method based on fraction order Fourier transform | |
CN102426354A (en) | Broadband radar detection method based on weighted sequence statistics and multiple-pulse coherence accumulation | |
CN112363121B (en) | C & I interference suppression method based on echo preprocessing and coherent accumulation | |
CN106199549B (en) | A method of LFMCW radar signal-to-noise ratio is promoted using spectrum-subtraction | |
CN109116326B (en) | Self-adaptive radar sea clutter suppression method based on median estimation | |
CN106443626A (en) | Unmanned area target detection method | |
CN101881824A (en) | Objective and fast determination method of noise threshold of power spectrum density data | |
CN106019243A (en) | Inter-pulse initial phase third-power modulation and FRFT (fractional Fourier transform)-based DRFM (digital radio frequency memory) interference detection and resisting method | |
EP0336273B1 (en) | Pulse Doppler radar | |
CN109085568B (en) | Frequency modulation continuous wave multi-target detection method based on secondary frequency mixing | |
CN102645648B (en) | Pulse accumulating method for improving target detection performance of ship radar | |
CN106646422A (en) | Preprocessing system for reinforcing signal-to-noise ratio of Doppler frequency shift signal of coherent wind finding radar | |
Haghshenas et al. | Suppressing sidelobe levels in random phase modulated radar | |
Ganveer et al. | SAR implementation using LFM signal | |
Zhao et al. | Range-Doppler sidelobe and clutter suppression via time-range adaptive processing | |
CN110045346B (en) | Laden Fourier transform blind speed side lobe suppression method | |
Kabakchiev et al. | CFAR BI detector for Mariner targets in time domain for bistatic forward scattering radar | |
Kong et al. | Fast implementation of adaptive multi-pulse compression via dimensionality reduction technique | |
CN102819012B (en) | Phase-control gust profile radar signal processing method based on non-negative ICA (independent component analysis) | |
Zhao et al. | Adaptive radar detection for targets in compound-Gaussian clutter with inverse-gamma texture | |
RU2802367C1 (en) | Method for selecting moving targets at high pulse repetition rate of a probing linear-frequency-modulated signal with a small duty cycle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200102 Address after: 214000 28 Weiming Road, Binhu District, Wuxi City, Jiangsu Province Patentee after: Aerospace new weather Technology Co., Ltd Address before: 100854 box 203, box 45, No. 142, Beijing Patentee before: No.23 Institute of No.20 Academy, CASIC |