CN102401896A - Method for calculating Doppler signal spectrum - Google Patents

Method for calculating Doppler signal spectrum Download PDF

Info

Publication number
CN102401896A
CN102401896A CN2011102496366A CN201110249636A CN102401896A CN 102401896 A CN102401896 A CN 102401896A CN 2011102496366 A CN2011102496366 A CN 2011102496366A CN 201110249636 A CN201110249636 A CN 201110249636A CN 102401896 A CN102401896 A CN 102401896A
Authority
CN
China
Prior art keywords
frequency
antennas
doppler
doppler signal
signal
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.)
Granted
Application number
CN2011102496366A
Other languages
Chinese (zh)
Other versions
CN102401896B (en
Inventor
姚斌
王磊
李苏红
贺雄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
HUNAN XIANGYI RAILROAD LOCOMOTIVE ELECTRICAL EQUIPMENT CO Ltd
Original Assignee
HUNAN XIANGYI RAILROAD LOCOMOTIVE ELECTRICAL EQUIPMENT CO Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by HUNAN XIANGYI RAILROAD LOCOMOTIVE ELECTRICAL EQUIPMENT CO Ltd filed Critical HUNAN XIANGYI RAILROAD LOCOMOTIVE ELECTRICAL EQUIPMENT CO Ltd
Priority to CN 201110249636 priority Critical patent/CN102401896B/en
Publication of CN102401896A publication Critical patent/CN102401896A/en
Application granted granted Critical
Publication of CN102401896B publication Critical patent/CN102401896B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a method for calculating Doppler signal spectrum, which sequentially comprises the following steps that: two antennas with different launch angles are arranged in a speed measurement radar system so as to acquire the Doppler signal; the original signals are subjected to band-pass or low-pass filtering and then FFT (Fast Fourier Transform Algorithm) transform to acquire the corresponding spectrums; the two spectrums are synthetized, and points of discontinuity of the two spectrums are removed; the strength and signal to noise ratio of the signal are assessed, and the proximal frequency points are searched in the spectrum synthesis; and the spectral line of the crossed spectral line is gathered to obtain the frequency estimate as well as smooth and calculate the parameters. The method disclosed by the invention can effectively avoid the error of calculation of the single Doppler signal spectrum caused by roadbed function, so that the relative running speed of a railway locomotive with high accuracy can be acquired.

Description

A kind of method that is used to calculate Doppler signal frequency spectrum
Technical field
The present invention relates to a kind of method that is used to calculate Doppler signal frequency spectrum, especially relate to a kind of method that in railway locomotive velocity radar system, is used for the signal spectrum of computing machine vehicle speed.
Background technology
The frequency spectrum Doppler technology is widely used in auto industry as a kind of method of basis; And for rail vehicles; The rapid variation of considering track ground must cause the variation of Doppler frequency spectrum, and the center of Doppler frequency spectrum also can change simultaneously, and causes that the correction factor of speed is also changing; To this problem, we have introduced a kind of computing method of the railway locomotive travelling speed being carried out valuation according to the intersection frequency of double antenna frequency spectrum.
Summary of the invention
The purpose of this invention is to provide a kind of method that is used to calculate Doppler signal frequency spectrum; It can be according to the Doppler signal amplitude, centre frequency, the bandwidth that estimate; Utilize the singularity of Doppler frequency spectrum, calculate required parameter exactly, to overcome the various defectives that prior art exists.
The so-called a kind of method that is used to calculate Doppler signal frequency spectrum of the present invention is characterized in that it experiences following step successively:
In the velocity radar system, install emission angle be respectively 40 ° with 50 ° and two antennas that highly are 500mm are set, their distinguish radiation millimeter wave, with transmitted wave and reflection echo difference frequency, the signal that obtains is Doppler signal; Be with logical or LPF to original Doppler signal, should the filtering high frequency interference with noise in high fdrequency component, necessary anti-aliasing filter before also will accomplishing AD and sampling is to satisfy nyquist sampling theorem; Data to after the sampling are carried out the FFT conversion to obtain corresponding frequency spectrum; The frequency spectrum of two antennas is synthesized and removes catastrophe point; Intensity, signal to noise ratio (S/N ratio) to Doppler signal are assessed, immediate frequency during the search frequency spectrum is synthetic; The intersection frequency that searches is gathered the heart to obtain Frequency Estimation; Confirm filter range according to frequency change rate repeatedly, be with logical or LPF further to reduce parameter error; Calculate, two corresponding Doppler's expression formulas of antenna do
Figure BDA0000086643920000022
Add relational expression by two formulas arrangement gained
Figure BDA0000086643920000023
With F ι d - F d 2 + F d = F k - - - 4 , Can obtain F from relational expression 3 and 4 d(i.e. the Doppler frequency of 40 ° of antennas) also has
Figure BDA0000086643920000025
The i.e. Doppler frequency of 50 ° of antennas, substitution relational expression again 1 and 2 just can obtain corresponding railway locomotive travelling speed V, the f in the relational expression 1 0Be the transmission frequency of 40 ° of antennas, in the relational expression 2
Figure BDA0000086643920000026
Be the transmission frequency of 50 ° of antennas, the C in two formulas then is the light velocity 3 * 10 8M/s.
The invention has the beneficial effects as follows: when measuring the railway locomotive travelling speed, can avoid the influence of different surface state effectively, significantly improve the accuracy of parameter estimation, thereby obtain high-precision measurement result spectrum analysis according to these computing method.
Description of drawings
Fig. 1 is radar microwave radiation areas figure;
Fig. 2 is the vertical view of Fig. 1;
Fig. 3 is the spectrogram under the two microwave antenna ideal situations of vehicle-mounted velocity radar.
Among the figure: the frequency spectrum that 1, obtains under the 2 one multidirectional same sex reflection cases, the frequency spectrum that obtains under 3, the 4 one multidirectional different in nature reflection cases, F KThe intersection point of one frequency spectrum.
Embodiment
To measure the railway locomotive travelling speed is example, Vehicular radar system 500mm highly be provided with emission angle be respectively 40 ° with two antennas of 50 °, there is covering to a certain degree in the millimeter wave that two antennas are launched on reflecting surface.A kind of method that is used to calculate Doppler signal frequency spectrum provided by the invention experiences following step successively: with the transmitted wave of two antennas and reflection echo difference frequency to obtain Doppler signal; Be with logical or LPF to original Doppler signal, should the filtering high frequency interference and noise in high fdrequency component, also to accomplish necessary anti-aliasing filter before the sampling, to satisfy nyquist sampling theorem; Data to after the sampling are carried out the FFT conversion to obtain corresponding frequency spectrum; The frequency spectrum of two antennas is synthesized and removes catastrophe point; Intensity, signal to noise ratio (S/N ratio) to Doppler signal are assessed, immediate frequency during the search frequency spectrum is synthetic; The Frequency Estimation that obtains that the spectral line of the intersection frequency that searches is gathered the heart; Confirm filter range according to frequency change rate repeatedly, be with logical or LPF with further minimizing parameter error; Calculate, two corresponding Doppler's expression formulas of antenna are:
With
Figure BDA0000086643920000032
Add relational expression by two formulas arrangement gained
Figure BDA0000086643920000033
Also have F ι d - F d 2 + F d = F k , Earlier try to achieve F from latter two relational expression dWith
Figure BDA0000086643920000035
Here it is we to obtain in isotropy (ground function U I=1) ideal and free from error Doppler frequency under the ideally-reflecting face, preceding two relational expressions of substitution just can obtain the corresponding relative travelling speed V of railway locomotive again.

Claims (1)

1. method that is used to calculate Doppler signal frequency spectrum is characterized in that it experiences following step successively:
1. in the velocity radar system, install emission angle be respectively 40 ° with 50 ° and two antennas that highly are 500mm are set, their distinguish radiation millimeter wave, with transmitted wave and reflection echo difference frequency, the signal that obtains is Doppler signal;
2. be with logical or LPF to original Doppler signal, should the filtering high frequency interference with noise in high fdrequency component, necessary anti-aliasing filter before also will accomplishing AD and sampling is to satisfy nyquist sampling theorem;
3. the data after the sampling are carried out the FFT conversion to obtain corresponding frequency spectrum;
4. the frequency spectrum of two antennas is synthesized and remove catastrophe point;
5. intensity, the signal to noise ratio (S/N ratio) of Doppler signal are assessed immediate frequency during the search frequency spectrum is synthetic;
6. the spectral line to the intersection frequency that searches gathers the heart to obtain Frequency Estimation;
7. confirm filter range based on frequency change rate repeatedly, be with logical or LPF further to reduce parameter error;
8. calculate, two corresponding Doppler's expression formulas of antenna do
Figure FDA0000086643910000011
Figure FDA0000086643910000012
Add relational expression by two formulas arrangement gained Also have F ι d - F d 2 + F d = F k - - - 4 , Can obtain F from relational expression 3 and 4 dPromptly the Doppler frequency of 40 ° of antennas also has
Figure FDA0000086643910000021
The i.e. Doppler frequency of 50 ° of antennas, substitution relational expression again 1 and 2 just can obtain corresponding rail locomotive travelling speed V, the f in the relational expression 1 0Be the transmission frequency of 40 ° of antennas, in the relational expression 2
Figure FDA0000086643910000022
Be the transmission frequency of 50 ° of antennas, the C in two formulas then is the light velocity 3 * 10 8M/s.
CN 201110249636 2011-08-26 2011-08-26 Method for calculating Doppler signal spectrum Active CN102401896B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110249636 CN102401896B (en) 2011-08-26 2011-08-26 Method for calculating Doppler signal spectrum

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110249636 CN102401896B (en) 2011-08-26 2011-08-26 Method for calculating Doppler signal spectrum

Publications (2)

Publication Number Publication Date
CN102401896A true CN102401896A (en) 2012-04-04
CN102401896B CN102401896B (en) 2013-04-24

Family

ID=45884300

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110249636 Active CN102401896B (en) 2011-08-26 2011-08-26 Method for calculating Doppler signal spectrum

Country Status (1)

Country Link
CN (1) CN102401896B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104635231A (en) * 2015-02-12 2015-05-20 西安电子科技大学 Locomotive radar speed measuring method based on adaptive Kalman filter
CN104678367A (en) * 2013-11-28 2015-06-03 中国航空工业集团公司雷华电子技术研究所 Velocity deception interference resisting method for airborne pulse Doppler radar
CN105334511A (en) * 2015-11-03 2016-02-17 周青松 Railroad locomotive velocity measurement method based on extraction of Doppler signal frequency
CN105372649A (en) * 2015-11-03 2016-03-02 黄中瑞 Railroad locomotive speed measurement method based on adaptive double-path sampling filtering
CN112557420A (en) * 2019-09-26 2021-03-26 深圳市雷立科技有限公司 Probe of microwave moisture meter

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5252980A (en) * 1992-07-23 1993-10-12 The United States Of America As Represented By The Secretary Of The Air Force Target location system
DE4406404C2 (en) * 1994-02-26 2000-01-20 Daimler Chrysler Ag Method and arrangement for determining the position and movement of an object
CN1257695A (en) * 1999-12-30 2000-06-28 复旦大学 Dual-ultrasonic Doppler method for measuring blood flow speed
KR20050103400A (en) * 2004-04-26 2005-10-31 이시동 Radio position finder with radio-direction finder
CN101320086A (en) * 2008-06-27 2008-12-10 北京航空航天大学 Echo signal processing equipment and method of Doppler speed measuring laser radar

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5252980A (en) * 1992-07-23 1993-10-12 The United States Of America As Represented By The Secretary Of The Air Force Target location system
DE4406404C2 (en) * 1994-02-26 2000-01-20 Daimler Chrysler Ag Method and arrangement for determining the position and movement of an object
CN1257695A (en) * 1999-12-30 2000-06-28 复旦大学 Dual-ultrasonic Doppler method for measuring blood flow speed
KR20050103400A (en) * 2004-04-26 2005-10-31 이시동 Radio position finder with radio-direction finder
CN101320086A (en) * 2008-06-27 2008-12-10 北京航空航天大学 Echo signal processing equipment and method of Doppler speed measuring laser radar

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104678367A (en) * 2013-11-28 2015-06-03 中国航空工业集团公司雷华电子技术研究所 Velocity deception interference resisting method for airborne pulse Doppler radar
CN104678367B (en) * 2013-11-28 2017-02-22 中国航空工业集团公司雷华电子技术研究所 Velocity deception interference resisting method for airborne pulse Doppler radar
CN104635231A (en) * 2015-02-12 2015-05-20 西安电子科技大学 Locomotive radar speed measuring method based on adaptive Kalman filter
CN105334511A (en) * 2015-11-03 2016-02-17 周青松 Railroad locomotive velocity measurement method based on extraction of Doppler signal frequency
CN105372649A (en) * 2015-11-03 2016-03-02 黄中瑞 Railroad locomotive speed measurement method based on adaptive double-path sampling filtering
CN112557420A (en) * 2019-09-26 2021-03-26 深圳市雷立科技有限公司 Probe of microwave moisture meter

Also Published As

Publication number Publication date
CN102401896B (en) 2013-04-24

Similar Documents

Publication Publication Date Title
CN107683422B (en) Vehicle radar system
CN102401896B (en) Method for calculating Doppler signal spectrum
EP2584373B1 (en) Radar device
Song et al. The velocity and range detection using the 2D-FFT scheme for automotive radars
CN101738611B (en) Underwater acoustic target signal detection and identification method
CN106199582B (en) Method applied to indoor two-frequency CW radar human body locating and tracking
CN102798856B (en) Small-wavelet-packet-based 24GHz LFMCW (Linear Frequency Modulation Continuous Wave) radar ranging method
CN103403572B (en) Vehicle inclination detection device
WO2014106907A1 (en) Radar device
CN101363913B (en) Signal processing method of fan-shaped beam circular cone scanning microwave scatterometer
JP2017166985A (en) Rader system and target object detection method
CN107783121B (en) Unmanned automobile anti-collision radar system signal processing system and method based on combined waveform
CN109596856B (en) Flow velocity measuring method applied to radar wave flow velocity meter and radar wave flow velocity meter
US8884814B2 (en) Processing method for FMCW radar signal with dual pulse repetition frequency
CN112224446B (en) High-speed coaxial dual-rotor blade tip distance measuring method based on phase distance measuring principle
JP7075925B2 (en) Speed measuring device and speed measuring method
CN113805166A (en) Target tracking and ranging method and system of radar level meter
CN111929685A (en) Radar detection method and device based on linear frequency modulation continuous wave, radar device and motor vehicle
CN106501796B (en) A kind of locomotive velocity measuring method, apparatus and system
JP2009014405A (en) In-vehicle radar apparatus
CN113812932A (en) Millimeter wave radar-based in-vehicle life target detection and positioning technology
CN109324327A (en) It is a kind of based on fmcw radar precision ranging and the method and system to test the speed
JP3773779B2 (en) Radar signal processing device
JP2000081480A (en) Fmcw radar, storage medium, and vehicle control device
CN105334511A (en) Railroad locomotive velocity measurement method based on extraction of Doppler signal frequency

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant