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 PDFInfo
- 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
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
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:
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,
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:
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:
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;
Be the angular frequency rate variance of adjacent spots, following relation arranged with difference frequency:
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
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
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:
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
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,
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:
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;
Be the angular frequency rate variance of adjacent spots, following relation arranged with difference frequency,
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,
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, can calculate the height value h of height indicator through following formula
0
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
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,
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:
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:
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;
Be the angular frequency rate variance of adjacent spots, following relation arranged with difference frequency:
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
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
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:
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
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)
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)
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 |
-
2012
- 2012-05-25 CN CN 201210166779 patent/CN102707275B/en not_active Expired - Fee Related
Patent Citations (4)
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)
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 |