CN102707275A - Digital processing method of altimeter of linear frequency modulation continuous wave radar - Google Patents

Digital processing method of altimeter of linear frequency modulation continuous wave radar Download PDF

Info

Publication number
CN102707275A
CN102707275A CN2012101667795A CN201210166779A CN102707275A CN 102707275 A CN102707275 A CN 102707275A CN 2012101667795 A CN2012101667795 A CN 2012101667795A CN 201210166779 A CN201210166779 A CN 201210166779A CN 102707275 A CN102707275 A CN 102707275A
Authority
CN
China
Prior art keywords
frequency
signal
amplitude
slope
spectral line
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
CN2012101667795A
Other languages
Chinese (zh)
Other versions
CN102707275B (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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN 201210166779 priority Critical patent/CN102707275B/en
Publication of CN102707275A publication Critical patent/CN102707275A/en
Application granted granted Critical
Publication of CN102707275B publication Critical patent/CN102707275B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a digital processing method of an altimeter of a linear frequency modulation continuous wave radar, which belongs to the field of the radio altimeters. The digital processing method comprises the following steps: step 1, obtaining beating signals from an orthogonal mixer arranged at the front end of the altimeter of the radar, and carrying out A/D (analog-to-digital) sampling to obtain a digital signal sequence; step 2, performing Fourier transform on the digital signal sequence, and calculating the amplitude of a spectral line; step 3, judging according to the amplitude of the spectral line of an amplitude spectrum, and measuring the beating frequency of the digital signal sequence; step 4, refining the frequency spectrum, calculating the slope correction of saw-tooth wave, and correcting the slope of the voltage signals of the saw-tooth wave; and step 5, calculating the high value of the current altimeter according to the slope of the voltage signals of the saw-tooth wave. The FFT-CZT method adopted by the invention is insensitive to noises and can perfectly work at lower signal-to-noise ratio, but the measurement precision is reduced. Because of having no limit on the properties of the targets, the digital processing method can be applied to the measurement and the analysis of point targets or surface targets.

Description

A kind of digitalized processing method of linear frequency modulation continuous wave radar altimeter
Technical field
The present invention relates to a kind of digitalized processing method of linear frequency modulation continuous wave radar altimeter, belong to the radio altimeter field.
Background technology
The linear frequency modulation continuous wave radar altimeter; It is the LFMCW radar altimeter; It is a kind of height finder operator's section that obtains height value through echoed signal and the Beat Signal that transmits; Radiation power is little, measuring accuracy is high, equipment is simple relatively, good electron is resisted advantages such as performance and low probability of intercept owing to have, and has been widely used in various aircraft such as unmanned plane at present.LFMCW radar altimeter early adopts the mode of average frequency counting to obtain the elevation information in the Beat Signal; Afterwards; Arrived the eighties mid-term, because the especially application of microprocessor in height indicator of digital technology makes CW with frequency modulation system radar altimeter obtain remarkable progress.Be different from previous time domain processing mode; Begin to adopt Digital Signal Processing to realize extraction for the Beat Signal that obtains to elevation information; Wherein the most frequently used algorithm is FFT, but because FFT is the spectrum analysis of on the whole unit circle in Z territory, carrying out, so its frequency resolution is f s/ N, wherein f sBe SF, N is the effective length of signal x (n), wants to improve frequency resolution and then realizes that high altimetry precision need increase the intercepted length to signal x (n), but calculated amount is increased greatly.Simultaneously; When adopting traditional FFT that signal is carried out spectrum analysis; Can receive the influence of aliasing effect, quantization error, leakage effect and fence effect, " fence effect " of FFT especially, the distance spectrum that directly adopts FFT and obtained has fixed distance Δ R (range resolution of radar); Thereby produce the range error of Δ R/2, this relative error that makes radar altimeter under low clearance, measure is bigger.
Summary of the invention
To the problem that exists in the prior art, the present invention proposes a kind of digitalized processing method of linear frequency modulation continuous wave radar altimeter, in permanent beat system LFMCW radar altimeter system, uses the FFT-CZT method and improves measuring accuracy.The FFT-CZT method is a kind of method that both combine with FFT and Chirp-Z; Basic ideas are to utilize FFT that the beat frequency spectrum of echoed signal is carried out bigness scale earlier; The frequency band that uses the Chirp-Z transfer pair to set then carries out local refinement; Reduce the SI on the beat frequency spectrum, reach the purpose that improves measuring accuracy.The calculated amount of the method is compared with the FFT method of direct increase sampling number and is wanted much less.Replace traditional frequency discriminator that Beat Signal is handled with the FFT-CZT method; Calculate the output height value after the tracking attitude of proofreading and correct through bigness scale search and to the cycle of modulation sawtooth wave, simulation result shows that this kind signal processing method can realize high-acruracy survey.
The present invention proposes a kind of digitalized processing method of linear frequency modulation continuous wave radar altimeter, mainly comprises following step:
Step 1: obtain Beat Signal from the front end orthogonal mixer of radar altimeter, the difference frequency ratings is f P, Beat Signal is carried out the A/D sampling, obtain digital signal sequences x (n);
Step 2: digital signal sequences is carried out Fast Fourier Transform (FFT) and calculates the amplitude of each bar spectral line;
Obtain the Fourier transform of digital signal sequences according to the discrete Fourier transformation formula, the transformation relation formula is following:
X ( k ) = Σ n = 0 N - 1 x ( n ) W N kn , k = 0,1,2 , . . . , N - 1 x ( n ) = 1 N Σ k = 0 N - 1 X ( k ) W N - kn , n = 0,1,2 , . . . , N - 1 - - - ( 1 )
Wherein: x (n) is a digital signal sequences, and this sequence length is N, and X (k) is the plural expression-form of the Fourier transform of digital signal sequences x (n), W NBe the Fourier transform factor,
Figure BDA00001682909200022
N and k represent the position of the middle each point of Fourier transform sequence X (k) of digital signal sequences x (n) and x (n) sequence respectively; The Fourier transform sequence X (k) that obtains digital signal sequences x (n) afterwards; Through Fourier transform sequence X (k) delivery is obtained amplitude spectrum; And the amplitude of each bar spectral line in the calculating amplitude spectrum, the amplitude threshold value is set simultaneously;
Step 3: adjudicate and measure the difference frequency of digital signal sequences x (n) according to amplitude spectrum spectral line amplitude, if difference frequency is at [f P-Δ f P, f P+ Δ f P] in the scope, write down spectral line positions in Fourier transform sequence X (k) that all surpass the amplitude threshold value, get into step 4; Otherwise the slope of the sawtooth voltage signal that the adjustment saw-toothed wave generator produces is specially: in amplitude spectrum, surpass the corresponding frequency of the spectral line of amplitude threshold value less than f P-Δ f P, then increase the slope of sawtooth wave, in amplitude spectrum, surpass the corresponding frequency of the spectral line of amplitude threshold value greater than f P+ Δ f P, then reduce the slope of sawtooth wave, return step 1 and search for again, until the spectral line difference frequency that in amplitude spectrum, surpasses the amplitude threshold value can search at least one spectral line in this specified hunting zone till;
Step 4: through Chirp-Z transfer pair frequency band [f P-Δ f P, f P+ Δ f P] carry out the frequency spectrum refinement, calculate sawtooth wave slope correction, the slope of the sawtooth voltage signal that current saw-toothed wave generator is produced is revised, and makes difference frequency stable to difference frequency ratings f P
Difference frequency is at specified hunting zone [f P-Δ f P, f P+ Δ f P] in, do the frequency spectrum refinement through this frequency band of Chirp-Z transfer pair and analyze, obtain the exact value of difference frequency; Chirp-Z transformation relation formula is:
X ( k ) = Σ n = 0 N - 1 x ( n ) z k - n - - - ( 2 )
z k=AW -k,0≤k≤M-1(3)
Wherein: x (n) is a digital signal sequences, and sequence length is N, and X (k) is the plural expression-form of the transform of digital signal sequences x (n), and M is that the frequency sampling in the frequency band of being analyzed is counted, and A is the starting point of Chirp-Z conversion, and W is the ratio between X (k) each point, z kBe the sampled point of Chirp-Z conversion, k representes the position of X (k) sequence each point, and n representes the position of each point among the digital signal sequences x (n), and A and W are expressed as with polar form:
Figure BDA00001682909200031
A wherein 0Vector radius for the sampling starting point; W 0Be the inside still bandy parameter of control sampling wheel profile; The transform path of Chirp-Z conversion is one section circular arc on the unit circle, θ 0Phase angle for starting point;
Figure BDA00001682909200032
Be the angular frequency rate variance of adjacent spots, following relation arranged with difference frequency:
Figure BDA00001682909200033
F wherein 1Be the initial frequency of refinement analysis frequency band, f 2For stopping frequency, f 1=f P-Δ f P, f 2=f P+ Δ f P, by (5) formula, in this refinement frequency band, frequency resolution does
Figure BDA00001682909200034
Obtain finding the maximum spectral line of amplitude in its amplitude spectrum behind the X (k), the k value that spectral line is corresponding, the exact value that obtains difference frequency is Δ f,
Figure BDA00001682909200035
Ratings f with Δ f and difference frequency PCompare, if the exact value Δ f of difference frequency equals f P, or the exact value Δ f of difference frequency and ratings f PThe absolute value of difference less than frequency resolution, the slope of record sawtooth voltage signal this moment; If the exact value Δ f of difference frequency is not equal to f P, the exact value Δ f and the ratings f of while difference frequency PThe absolute value of difference more than or equal to frequency resolution, the correction that then calculates the sawtooth voltage signal slope does The slope of the sawtooth voltage signal that the adjustment saw-toothed wave generator produces carries out step 4 again and proofreaies and correct, until meeting the demands record sawtooth voltage signal slope γ 0
Step 5: the height value that calculates current height indicator according to the sawtooth voltage signal slope:
The sawtooth voltage signal γ that utilizes step 4 to obtain 0, go out the height value h of height indicator through computes 0:
h 0 = f p γ 0 · c 2 - - - ( 6 )
Wherein, c is a radio wave propagation speed.
The invention has the advantages that:
(1) the present invention proposes a kind of digitalized processing method of linear frequency modulation continuous wave radar altimeter, and the FFT-CZT method of employing is the basis with FFT, and Fast Fourier Transform (FFT) is quite ripe, and a large amount of general, special-purpose software and hardware supports is arranged;
(2) the present invention proposes a kind of digitalized processing method of linear frequency modulation continuous wave radar altimeter, and the FFT-CZT method that the present invention adopts still can be worked under lower signal to noise ratio (S/N ratio) to insensitive for noise preferably, and just measuring accuracy can descend to some extent;
(3) the present invention proposes a kind of digitalized processing method of linear frequency modulation continuous wave radar altimeter, and to the not restriction of characteristic of target, therefore both can be used for point target also can be used for appearance target Measurement and analysis.
Description of drawings
Fig. 1: the process flow diagram of the digitalized processing method of a kind of linear frequency modulation continuous wave radar altimeter that the present invention proposes;
Fig. 2: the signal flow diagram of the digitalized processing method of a kind of linear frequency modulation continuous wave radar altimeter that the present invention proposes.
Among the figure: the 1-saw-toothed wave generator; The 2-radio frequency source; 3-power amplifier A; The 4-circulator;
The 5-antenna; The 6-receiving front-end offsets processing module; 7-power amplifier B;
The 8-orthogonal mixer; The 9-low-frequency amplifier; The 10-correcting circuit.
Embodiment
To combine accompanying drawing that the present invention is done further detailed description below.
The present invention proposes a kind of digitalized processing method of linear frequency modulation continuous wave radar altimeter, and is as shown in Figure 1, comprises following step:
Step 1: obtain Beat Signal from the front end orthogonal mixer of radar altimeter, and Beat Signal is carried out obtaining digital signal sequences x (n) after A/D sampling, the quantification;
Obtain Beat Signal x (t) from the front end orthogonal mixer of radar altimeter, the difference frequency ratings of this Beat Signal is f p, for continuous permanent Beat Signal, carry random noise signal, this Beat Signal is carried out the A/D sampling, SF is f s, after over-sampling, obtain digital signal sequences x (n).
The procurement process of described Beat Signal: as shown in Figure 2; The sawtooth voltage signal controlling radio frequency source 2 that is produced by saw-toothed wave generator 1 produces the linear frequency modulation continuous wave signal; This linear frequency modulation continuous wave signal a part of signal after power amplifier A3 power amplification is fed to antenna 5 through circulator 4; By antenna 5 signal radiation is gone out; The signal that receives by antenna 5 again; Signal inputs to receiving front-end through circulator 4 and offsets processing module 6, offsets processing with another part signal from power amplifier A3, after offseting processing, carries out power amplification through power amplifier B7 again; Signal after amplifying gets into orthogonal mixer 8 and carries out Frequency mixing processing; The signal of this Frequency mixing processing is that signal and the radio frequency source 2 after power amplifier B7 power amplification produces the linear frequency modulation continuous wave signal, exports Beat Signal after orthogonal mixer 8 Frequency mixing processing, and this Beat Signal is amplified into the digital signal processing link through low-frequency amplifier 9 and handles; Be specially in order through the Fourier transform of step 2, the search of step 3 Beat Signal frequency, the processes such as tracking processing of step 4, produce the slope that correction signals are controlled the sawtooth voltage signal that saw-toothed wave generator produces through control correcting circuit 10 at last.
Step 2: digital signal sequences is carried out Fast Fourier Transform (FFT) and calculates the amplitude of each bar spectral line;
Obtain the Fourier transform of digital signal sequences according to the discrete Fourier transformation formula, the transformation relation formula is following
X ( k ) = Σ n = 0 N - 1 x ( n ) W N kn , k = 0,1,2 , . . . , N - 1 x ( n ) = 1 N Σ k = 0 N - 1 X ( k ) W N - kn , n = 0,1,2 , . . . , N - 1 - - - ( 1 )
Wherein: x (n) is a digital signal sequences, and this sequence length is N, and X (k) is the plural expression-form of the Fourier transform of digital signal sequences x (n), W NBe the Fourier transform factor,
Figure BDA00001682909200052
N and k represent the position of each point among digital signal sequences x (n) and its Fourier transform X (k) respectively.The Fourier transform sequence X (k) that obtains digital signal sequences x (n) through Fourier transform sequence X (k) delivery is obtained amplitude spectrum, and is calculated the amplitude of each bar spectral line in the amplitude spectrum afterwards, and the amplitude threshold value is set simultaneously.
Step 3: adjudicate and measure the difference frequency of digital signal sequences x (n) according to amplitude spectrum spectral line amplitude, if difference frequency is at difference frequency ratings f PNear, promptly at [f P-Δ f P, f P+ Δ f P] in the scope, then get into step 4, get into tracking mode; Otherwise the slope of the sawtooth voltage signal that the adjustment saw-toothed wave generator produces returns step 1; Difference frequency ratings f PDesign parameter by system provides, generally much smaller than SF f s, Δ f PFor the difference frequency float value, be generally 20%f P
Make comparisons the amplitude threshold value of the amplitude of every spectral line in the amplitude spectrum and setting, write down the position of spectral line in Fourier transform sequence X (k) that all surpass the amplitude threshold value, these spectral lines that surpass the amplitude threshold value are concentrated in the i.e. [f of one section narrower frequency band P-Δ f P, f P+ Δ f P] in the scope, if be in the spectral line in this scope, the slope of the sawtooth voltage signal that the adjustment saw-toothed wave generator produces returns step 1, searches for next time.According to the character of Fast Fourier Transform (FFT), N bar spectral line is to be evenly distributed in 0 on the frequency axis, f sBetween, f sBe the SF of A/D sampling in the step 1, therefore the position by each bar spectral line place just can obtain the corresponding difference frequency of each bar spectral line.If be not in [f P-Δ f P, f P+ Δ f P] spectral line of scope, two kinds of situation are then arranged, in amplitude spectrum, surpass the corresponding frequency of the spectral line of amplitude threshold value less than f P-Δ f P, then increase the slope of sawtooth wave; In amplitude spectrum, surpass the corresponding frequency of the spectral line of amplitude threshold value greater than f P+ Δ f P, then reduce the slope of sawtooth wave; Return step 1 and search for again, until the spectral line difference frequency that in amplitude spectrum, surpasses the amplitude threshold value at this specified hunting zone [f P-Δ f P, f P+ Δ f P] interior can institute searching till at least one spectral line.
Step 4: be in [f through the Chirp-Z transfer pair P-Δ f P, f P+ Δ f P] spectral line frequency band in the scope carries out the frequency spectrum refinement, calculates sawtooth wave slope correction, the sawtooth voltage signal that current saw-toothed wave generator is produced is revised and is made difference frequency stable to difference frequency ratings f P
Search procedure finishes the back and gets into tracing process.This moment, difference frequency was at specified hunting zone [f P-Δ f P, f P+ Δ f P] in, do the frequency spectrum refinement through this frequency band of Chirp-Z transfer pair and analyze, obtain the exact value of difference frequency.Chirp-Z transformation relation formula is following:
X ( k ) = Σ n = 0 N - 1 x ( n ) z k - n - - - ( 2 )
z k=AW -k,0≤k≤M-1(3)
Wherein: x (n) is a digital signal sequences, and sequence length is N, and X (k) is the plural expression-form of the transform of digital signal sequences x (n).M is that the frequency sampling in the frequency band of being analyzed is counted.A and W are plural arbitrarily, and A is the starting point of Chirp-Z conversion, and W is the ratio between X (k) each point.z kBe the sampled point of Chirp-Z conversion, k representes the position of each point in X (k) sequence, and n representes the position of each point among the digital signal sequences x (n).A and W are expressed as with polar form
Wherein: A 0Vector radius for the sampling starting point; W 0Be the inside still bandy parameter of control sampling wheel profile; Get A among the present invention 0=W 0=1, this moment, the transform path of Chirp-Z conversion was one section circular arc on the unit circle.θ 0Phase angle for starting point;
Figure BDA00001682909200063
Be the angular frequency rate variance of adjacent spots, following relation arranged with difference frequency,
Figure BDA00001682909200064
In the formula, f 1Be the initial frequency of refinement analysis frequency band, f 2For stopping frequency, in the present invention, f 1=f P-Δ f P, f 2=f P+ Δ f PCan know by (5) formula; In this refinement frequency band, frequency resolution is
Obtain finding the maximum spectral line of amplitude in its amplitude spectrum behind the X (k), the k value that spectral line is corresponding, the exact value that obtains difference frequency is Δ f,
Figure BDA00001682909200066
Ratings f with Δ f and difference frequency PCompare, if the exact value Δ f of difference frequency equals f P, or the exact value Δ f of difference frequency and ratings f PThe absolute value of difference less than frequency resolution, the slope of record sawtooth voltage signal this moment; If the exact value Δ f of difference frequency is not equal to f P, the exact value Δ f and the ratings f of while difference frequency PThe absolute value of difference more than or equal to frequency resolution, the correction that then calculates the sawtooth voltage signal slope does
Figure BDA00001682909200067
The slope of the sawtooth voltage signal that the adjustment saw-toothed wave generator produces carries out step 4 again and proofreaies and correct, until meeting the demands record sawtooth voltage signal slope γ 0
Step 5: the height value that calculates current height indicator according to the sawtooth voltage signal slope.
The sawtooth voltage signal γ that utilizes step 4 to obtain 0, can calculate the height value h of height indicator through following formula 0
h 0 = f p γ 0 · c 2 - - - ( 6 )
Wherein, c is a radio wave propagation speed.

Claims (3)

1. the digitalized processing method of a linear frequency modulation continuous wave radar altimeter is characterized in that: comprise following step:
Step 1: obtain Beat Signal from the front end orthogonal mixer of radar altimeter, the difference frequency ratings is f P, Beat Signal is carried out the A/D sampling, obtain digital signal sequences x (n);
Step 2: digital signal sequences is carried out Fast Fourier Transform (FFT) and calculates the amplitude of each bar spectral line;
Obtain the Fourier transform of digital signal sequences according to the discrete Fourier transformation formula, the transformation relation formula is following
X ( k ) = Σ n = 0 N - 1 x ( n ) W N kn , k = 0,1,2 , . . . , N - 1 x ( n ) = 1 N Σ k = 0 N - 1 X ( k ) W N - kn , n = 0,1,2 , . . . , N - 1 - - - ( 1 )
Wherein: x (n) is a digital signal sequences, and this sequence length is N, and X (k) is the plural expression-form of the Fourier transform of digital signal sequences x (n), W NBe the Fourier transform factor,
Figure FDA00001682909100012
N and k represent the position of the middle each point of Fourier transform sequence X (k) of digital signal sequences x (n) and x (n) sequence respectively; The Fourier transform sequence X (k) that obtains digital signal sequences x (n) afterwards; Through Fourier transform sequence X (k) delivery is obtained amplitude spectrum; And the amplitude of each bar spectral line in the calculating amplitude spectrum, the amplitude threshold value is set simultaneously;
Step 3: adjudicate and measure the difference frequency of digital signal sequences x (n) according to amplitude spectrum spectral line amplitude, if difference frequency is at [f P-Δ f P, f P+ Δ f P] in the scope, write down spectral line positions in Fourier transform sequence X (k) that all surpass the amplitude threshold value, get into step 4; Otherwise the slope of the sawtooth voltage signal that the adjustment saw-toothed wave generator produces is specially: in amplitude spectrum, surpass the corresponding frequency of the spectral line of amplitude threshold value less than f P-Δ f P, then increase the slope of sawtooth wave, in amplitude spectrum, surpass the corresponding frequency of the spectral line of amplitude threshold value greater than f P+ Δ f P, then reduce the slope of sawtooth wave, return step 1 and search for again, until the spectral line difference frequency that in amplitude spectrum, surpasses the amplitude threshold value can search at least one spectral line in this specified hunting zone till;
Step 4: be in [f through the Chirp-Z transfer pair P-Δ f P, f P+ Δ f P] spectral line frequency band in the scope carries out the frequency spectrum refinement, calculates sawtooth wave slope correction, the slope of the sawtooth voltage signal that current saw-toothed wave generator is produced is revised, and makes difference frequency stable to difference frequency ratings f P
Difference frequency is at specified hunting zone [f P-Δ f P, f P+ Δ f P] in, do the frequency spectrum refinement through this frequency band of Chirp-Z transfer pair and analyze, obtain the exact value of difference frequency; Chirp-Z transformation relation formula is:
X ( k ) = Σ n = 0 N - 1 x ( n ) z k - n - - - ( 2 )
z k=AW -k,0≤k≤M-1(3)
Wherein: x (n) is a digital signal sequences, and sequence length is N, and X (k) is the plural expression-form of the transform of digital signal sequences x (n), and M is that the frequency sampling in the frequency band of being analyzed is counted, and A is the starting point of Chirp-Z conversion, and W is the ratio between X (k) each point, z kBe the sampled point of Chirp-Z conversion, k representes the position of X (k) sequence each point, and n representes the position of each point among the digital signal sequences x (n), and A and W are expressed as with polar form:
Figure FDA00001682909100021
A wherein 0Vector radius for the sampling starting point; W 0Be the inside still bandy parameter of control sampling wheel profile; The transform path of Chirp-Z conversion is one section circular arc on the unit circle, θ 0Phase angle for starting point;
Figure FDA00001682909100022
Be the angular frequency rate variance of adjacent spots, following relation arranged with difference frequency:
Figure FDA00001682909100023
F wherein 1Be the initial frequency of refinement analysis frequency band, f 2For stopping frequency, f 1=f P-Δ f P, f 2=f P+ Δ f P, by (5) formula, in this refinement frequency band, frequency resolution does
Figure FDA00001682909100024
Obtain finding the maximum spectral line of amplitude in its amplitude spectrum behind the X (k), the k value that spectral line is corresponding, the exact value that obtains difference frequency is Δ f, Ratings f with Δ f and difference frequency PCompare, if the exact value Δ f of difference frequency equals f P, or the exact value Δ f of difference frequency and ratings f PThe absolute value of difference less than frequency resolution, the slope of record sawtooth voltage signal this moment; If the exact value Δ f of difference frequency is not equal to f P, the exact value Δ f and the ratings f of while difference frequency PThe absolute value of difference more than or equal to frequency resolution, the correction that then calculates the sawtooth voltage signal slope does
Figure FDA00001682909100026
The slope of the sawtooth voltage signal that the adjustment saw-toothed wave generator produces carries out step 4 again and proofreaies and correct, until meeting the demands record sawtooth voltage signal slope γ 0
Step 5: the height value that calculates current height indicator according to the sawtooth voltage signal slope:
The sawtooth voltage signal γ that utilizes step 4 to obtain 0, go out the height value h of height indicator through computes 0:
h 0 = f p γ 0 · c 2 - - - ( 6 )
Wherein, c is a radio wave propagation speed.
2. the digitalized processing method of a kind of linear frequency modulation continuous wave radar altimeter according to claim 1; It is characterized in that: the concrete procurement process of the Beat Signal in the step 1 is: the sawtooth voltage signal controlling radio frequency source that is produced by saw-toothed wave generator produces the linear frequency modulation continuous wave signal; This linear frequency modulation continuous wave signal a part of signal after power amplifier A power amplification is fed to antenna through circulator; By antenna signal radiation is gone out; By antenna receiving signal, signal inputs to receiving front-end through circulator and offsets processing module again, offsets processing with another part signal from power amplifier A; Handle after power amplifier B carries out power amplification through offseting; Signal after amplifying gets into orthogonal mixer and carries out Frequency mixing processing, and the signal of this Frequency mixing processing is that signal and the radio frequency source after power amplifier B power amplification produces the linear frequency modulation continuous wave signal, exports Beat Signal after the orthogonal mixer Frequency mixing processing.
3. the digitalized processing method of a kind of linear frequency modulation continuous wave radar altimeter according to claim 1 is characterized in that: described difference frequency float value Δ f pBe 20%f p
CN 201210166779 2012-05-25 2012-05-25 Digital processing method of altimeter of linear frequency modulation continuous wave radar Expired - Fee Related CN102707275B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210166779 CN102707275B (en) 2012-05-25 2012-05-25 Digital processing method of altimeter of linear frequency modulation continuous wave radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210166779 CN102707275B (en) 2012-05-25 2012-05-25 Digital processing method of altimeter of linear frequency modulation continuous wave radar

Publications (2)

Publication Number Publication Date
CN102707275A true CN102707275A (en) 2012-10-03
CN102707275B CN102707275B (en) 2013-12-18

Family

ID=46900205

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210166779 Expired - Fee Related CN102707275B (en) 2012-05-25 2012-05-25 Digital processing method of altimeter of linear frequency modulation continuous wave radar

Country Status (1)

Country Link
CN (1) CN102707275B (en)

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103616677A (en) * 2013-12-05 2014-03-05 中国航空无线电电子研究所 Radio altimeter and air vehicle altitude measuring method
CN103884317A (en) * 2014-04-14 2014-06-25 西安阿尔特测控技术有限公司 Real-time monitoring system for settlement of roadbed of high-speed rail
CN105005028A (en) * 2015-06-29 2015-10-28 河北工业大学 LFMCW radar range finding early warning method based on genetic algorithm
CN105319544A (en) * 2014-06-26 2016-02-10 霍尼韦尔国际公司 Systems and methods for calibrating and optimizing frequency modulated continuous wave radar altimeters using adjustable self-interference cancellation
CN105487058A (en) * 2015-12-29 2016-04-13 大连楼兰科技股份有限公司 Method for accurate and rapid determination of radar system parameters
CN105675084A (en) * 2014-11-20 2016-06-15 中国航空工业集团公司雷华电子技术研究所 High-accuracy liquid-level measurement method with iterative interpolation
CN105842681A (en) * 2016-03-24 2016-08-10 北京工业大学 One-way FM continuous electromagnetic wave high-precision range finder
CN105866746A (en) * 2016-04-01 2016-08-17 芜湖航飞科技股份有限公司 Application of FMCW system T/R unit in digital phased array
CN106840310A (en) * 2017-02-17 2017-06-13 水利部南京水利水文自动化研究所 A kind of CW/FM continuous wave/frequency modulated modulated indicator water gage measuring method
CN107153189A (en) * 2017-04-18 2017-09-12 上海交通大学 The signal processing method of Modulation Continuous Wave Radar ranging
CN107389154A (en) * 2017-08-09 2017-11-24 武汉隽龙科技股份有限公司 The continuous liquid level sensing device of hollow optic fibre and measuring method based on OFDR
CN107390210A (en) * 2017-06-22 2017-11-24 电子科技大学 The digital processing method of Beat Signal in a kind of level gauging
CN107632292A (en) * 2017-09-21 2018-01-26 北京工业大学 A kind of method that frequency modulation Fourier transform is carried out to radar signal
CN108334682A (en) * 2018-01-24 2018-07-27 北京机械设备研究所 A kind of batch processing method that frequency spectrum refines automatically
CN108802724A (en) * 2018-06-28 2018-11-13 北京豪恩科技有限责任公司 A method of test wave heights
CN109581369A (en) * 2019-01-25 2019-04-05 西安电子科技大学 The radar altimeter of non-homogeneous multichannel perseverance difference frequency system
CN110045341A (en) * 2019-02-28 2019-07-23 西南电子技术研究所(中国电子科技集团公司第十研究所) Radar altimeter low interception performance test method
CN110488272A (en) * 2019-08-30 2019-11-22 成都纳雷科技有限公司 A kind of test the speed distance measuring method and device based on frequency modulated continuous wave radar
CN111025278A (en) * 2019-12-31 2020-04-17 北京中科飞鸿科技股份有限公司 Radio altimeter
CN111766597A (en) * 2020-07-08 2020-10-13 哈尔滨理工大学 Frequency measurement method and device for frequency modulation continuous wave laser ranging beat signal
CN113534145A (en) * 2021-09-16 2021-10-22 北京安达维尔航空设备有限公司 Altimeter height measurement method and system based on linear frequency modulation continuous wave system
CN116106886A (en) * 2023-04-11 2023-05-12 烟台初心航空科技有限公司 Airborne radio altimeter device based on frequency modulation continuous wave

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4028700A (en) * 1972-11-16 1977-06-07 Raytheon Company Pulse compression radar and method of operating such a radar
US5351053A (en) * 1993-07-30 1994-09-27 The United States Of America As Represented By The Secretary Of The Air Force Ultra wideband radar signal processor for electronically scanned arrays
CN101089653A (en) * 2007-07-20 2007-12-19 西安理工大学 Short-range frequency-modulation continuous wave FMCW radar anti-interference method
CN101210965A (en) * 2007-12-21 2008-07-02 当代天启技术(北京)有限公司 Wireless distance finding method, wireless distance finding and position finding method, apparatus and system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4028700A (en) * 1972-11-16 1977-06-07 Raytheon Company Pulse compression radar and method of operating such a radar
US5351053A (en) * 1993-07-30 1994-09-27 The United States Of America As Represented By The Secretary Of The Air Force Ultra wideband radar signal processor for electronically scanned arrays
CN101089653A (en) * 2007-07-20 2007-12-19 西安理工大学 Short-range frequency-modulation continuous wave FMCW radar anti-interference method
CN101210965A (en) * 2007-12-21 2008-07-02 当代天启技术(北京)有限公司 Wireless distance finding method, wireless distance finding and position finding method, apparatus and system

Cited By (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103616677B (en) * 2013-12-05 2015-07-01 中国航空无线电电子研究所 Radio altimeter and air vehicle altitude measuring method
CN103616677A (en) * 2013-12-05 2014-03-05 中国航空无线电电子研究所 Radio altimeter and air vehicle altitude measuring method
CN103884317A (en) * 2014-04-14 2014-06-25 西安阿尔特测控技术有限公司 Real-time monitoring system for settlement of roadbed of high-speed rail
CN105319544A (en) * 2014-06-26 2016-02-10 霍尼韦尔国际公司 Systems and methods for calibrating and optimizing frequency modulated continuous wave radar altimeters using adjustable self-interference cancellation
CN105675084A (en) * 2014-11-20 2016-06-15 中国航空工业集团公司雷华电子技术研究所 High-accuracy liquid-level measurement method with iterative interpolation
CN105675084B (en) * 2014-11-20 2019-03-29 中国航空工业集团公司雷华电子技术研究所 A kind of high-precision liquid level measurement method with iterated interpolation
CN105005028A (en) * 2015-06-29 2015-10-28 河北工业大学 LFMCW radar range finding early warning method based on genetic algorithm
CN105487058B (en) * 2015-12-29 2018-07-17 大连楼兰科技股份有限公司 The method that can accurately, quickly determine radar system parameters
CN105487058A (en) * 2015-12-29 2016-04-13 大连楼兰科技股份有限公司 Method for accurate and rapid determination of radar system parameters
CN105842681A (en) * 2016-03-24 2016-08-10 北京工业大学 One-way FM continuous electromagnetic wave high-precision range finder
CN105866746A (en) * 2016-04-01 2016-08-17 芜湖航飞科技股份有限公司 Application of FMCW system T/R unit in digital phased array
CN106840310A (en) * 2017-02-17 2017-06-13 水利部南京水利水文自动化研究所 A kind of CW/FM continuous wave/frequency modulated modulated indicator water gage measuring method
CN106840310B (en) * 2017-02-17 2020-06-26 水利部南京水利水文自动化研究所 Continuous frequency modulation wave radar level gauge measuring method
CN107153189A (en) * 2017-04-18 2017-09-12 上海交通大学 The signal processing method of Modulation Continuous Wave Radar ranging
CN107153189B (en) * 2017-04-18 2021-08-03 上海交通大学 Signal processing method for linear frequency modulation continuous wave radar ranging
CN107390210A (en) * 2017-06-22 2017-11-24 电子科技大学 The digital processing method of Beat Signal in a kind of level gauging
CN107389154A (en) * 2017-08-09 2017-11-24 武汉隽龙科技股份有限公司 The continuous liquid level sensing device of hollow optic fibre and measuring method based on OFDR
CN107389154B (en) * 2017-08-09 2023-04-25 武汉昊衡科技有限公司 Hollow fiber continuous liquid level sensing device and measuring method based on OFDR
CN107632292B (en) * 2017-09-21 2021-07-30 北京工业大学 Method for carrying out frequency modulation Fourier transform on radar signal
CN107632292A (en) * 2017-09-21 2018-01-26 北京工业大学 A kind of method that frequency modulation Fourier transform is carried out to radar signal
CN108334682A (en) * 2018-01-24 2018-07-27 北京机械设备研究所 A kind of batch processing method that frequency spectrum refines automatically
CN108802724A (en) * 2018-06-28 2018-11-13 北京豪恩科技有限责任公司 A method of test wave heights
CN109581369A (en) * 2019-01-25 2019-04-05 西安电子科技大学 The radar altimeter of non-homogeneous multichannel perseverance difference frequency system
CN109581369B (en) * 2019-01-25 2023-03-24 西安电子科技大学 Radar altimeter with non-uniform multi-channel constant difference beat frequency system
CN110045341A (en) * 2019-02-28 2019-07-23 西南电子技术研究所(中国电子科技集团公司第十研究所) Radar altimeter low interception performance test method
CN110045341B (en) * 2019-02-28 2022-06-14 西南电子技术研究所(中国电子科技集团公司第十研究所) Method for testing low interception performance of radar altimeter
CN110488272A (en) * 2019-08-30 2019-11-22 成都纳雷科技有限公司 A kind of test the speed distance measuring method and device based on frequency modulated continuous wave radar
CN111025278A (en) * 2019-12-31 2020-04-17 北京中科飞鸿科技股份有限公司 Radio altimeter
CN111025278B (en) * 2019-12-31 2023-09-01 北京中科飞鸿科技股份有限公司 Radio altimeter
CN111766597A (en) * 2020-07-08 2020-10-13 哈尔滨理工大学 Frequency measurement method and device for frequency modulation continuous wave laser ranging beat signal
CN111766597B (en) * 2020-07-08 2022-05-13 哈尔滨理工大学 Frequency measuring method and device for frequency-modulated continuous wave laser ranging beat signal
CN113534145A (en) * 2021-09-16 2021-10-22 北京安达维尔航空设备有限公司 Altimeter height measurement method and system based on linear frequency modulation continuous wave system
CN116106886A (en) * 2023-04-11 2023-05-12 烟台初心航空科技有限公司 Airborne radio altimeter device based on frequency modulation continuous wave

Also Published As

Publication number Publication date
CN102707275B (en) 2013-12-18

Similar Documents

Publication Publication Date Title
CN102707275B (en) Digital processing method of altimeter of linear frequency modulation continuous wave radar
US10557933B2 (en) Radar device and position-determination method
CN103616679B (en) Based on difference beam modulation and the PD radar range finding angle-measuring method of wave form analysis
CN102176016B (en) Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method
CN103744068B (en) The moving-target detection formation method of dual pathways Continuous Wave with frequency modulation SAR system
CN106443671A (en) SAR radar moving target detecting and imaging method based on FM continuous wave
CN100501425C (en) High-frequency chirp radar directional diagram measuring method
CN106840310B (en) Continuous frequency modulation wave radar level gauge measuring method
CN102749620A (en) Monopulse foresight imaging processing method of missile-borne/airborne radar
CN103235305A (en) Spaceborne ultrahigh-resolution sliding bunching SAR (synthetic aperture radar) imaging method
CN111220954A (en) Radar angle error correction method based on self-correcting amplitude normalization
CN105572668B (en) A kind of moving target radial velocity method of estimation converted based on Radon
US20190004175A1 (en) Object detecting device and sensor device
Xu et al. Comparison of speed over ground estimation using acoustic and radar Doppler sensors
EP3208633B1 (en) Method and system for fmcw radar altimeter system height measurement resolution improvement
CN109001671B (en) Target detection and parameter estimation method and device for frequency hopping signal
CN109507654B (en) LS-based phase information calculation method in complex environment
JP2010175457A (en) Radar apparatus
CN110988873A (en) Single-channel SAR ship speed estimation method and system based on energy center extraction
CN109901155B (en) Hypersonic target correction PRPFT detection method based on Doppler information feedback regulation
CN105548987B (en) A kind of continuous wave radar aimed acceleration blind estimating method
CN108020835A (en) A kind of strong clutter suppression method and device of spaceborne SAR ground moving object instruction GMTI
RU2608551C1 (en) Pulse-doppler airborne radar station operating method during detecting of aerial target, radio reconnaissance station carrier
CN110456317A (en) A kind of phased array radar system calibrating method based on meteor trail echo
CN113126051B (en) Airborne multichannel SAR interference effective baseline estimation method and device

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20131218

Termination date: 20140525