CN111505601A - Linear motion demodulation implementation method based on improved differential cross multiplication - Google Patents

Linear motion demodulation implementation method based on improved differential cross multiplication Download PDF

Info

Publication number
CN111505601A
CN111505601A CN202010433815.4A CN202010433815A CN111505601A CN 111505601 A CN111505601 A CN 111505601A CN 202010433815 A CN202010433815 A CN 202010433815A CN 111505601 A CN111505601 A CN 111505601A
Authority
CN
China
Prior art keywords
signal
algorithm
motion
demodulation
cross multiplication
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010433815.4A
Other languages
Chinese (zh)
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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN202010433815.4A priority Critical patent/CN111505601A/en
Publication of CN111505601A publication Critical patent/CN111505601A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415Identification of targets based on measurements of movement associated with the target

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

A linear motion demodulation method based on improved difference cross multiplication is based on trigonometric function geometric theorem, directly performs differential processing on an I/Q signal after direct current calibration without inverse orthotropic transformation, and then performs cross multiplication and logical addition operation on the I/Q signal to extract a motion signal to be detected.

Description

Linear motion demodulation implementation method based on improved differential cross multiplication
Technical Field
The invention relates to a technology in the field of radar detection, in particular to a linear motion demodulation implementation method based on improved differential cross multiplication.
Background
In the radar detection of the target motion track, a signal processing technology based on phase information is needed. The target movement causes the nonlinear phase modulation of the radar electromagnetic wave signal, so a high-linearity demodulation method is required. The most commonly used linear demodulation algorithms at present are both the arctan algorithm and the differential cross multiplication (DACM) algorithm. First, the arctangent algorithm is simple and highly linear, once the arctangent algorithm has attracted a lot of attention for linear demodulation. However, the arctangent algorithm has strict value-taking domain constraint on the motion signal to be measured, and the universal application of the algorithm is influenced. When the phase locus of the modulation signal is not in the unique quadrant or the range of the first quadrant and the fourth quadrant due to the motion information to be detected, the demodulation signal generates phase ambiguity, and the waveform is subjected to jumping truncation.
In order to break the limitation of the value range and make the demodulation algorithm applicable to the whole range, scientists further propose the DACM algorithm. The DACM algorithm is realized by performing differentiation and integration operations again on the basis of the arctan algorithm, which causes that the DACM algorithm needs to perform nested operation and a large amount of approximation, increases the accumulated noise and approximation error of demodulation results, and reduces the accuracy and stability of demodulation signals. In addition, both the arctangent algorithm and the DACM algorithm depend on the accuracy of the pre-dc offset calibration step excessively, and the robustness and the anti-noise capability of the algorithm are affected.
Disclosure of Invention
The invention provides a linear motion demodulation implementation method based on improved differential cross multiplication aiming at the phase ambiguity problem of the demodulation signals in the prior art, which utilizes the geometric theorem of difference and trigonometric function, does not need inverse tangential transformation, only relates to the first-order operation of I/Q signals, greatly simplifies the operation process, reduces operation approximation and accumulated noise, and improves the linearity, stability and noise resistance of the algorithm.
The invention is realized by the following technical scheme:
the invention relates to a linear motion demodulation implementation method based on improved differential cross multiplication, which is based on trigonometric function geometric theorem, directly performs differential processing on an I/Q signal after direct current calibration without inverse tangent transformation, and performs cross multiplication and logical addition operation on the I/Q signal and an original I/Q signal to extract a motion signal to be detected.
The original I/Q signals refer to: and amplifying, filtering and carrying out down-conversion treatment on the radar echo signal received by the radio frequency receiving end to obtain the radar echo signal.
Technical effects
The invention integrally solves the technical problems that in the prior art, the value taking domain of a measurable motion signal is strictly limited, the phase of a demodulated signal is fuzzy, the accuracy of a preposed direct current offset calibration algorithm is excessively depended on, the linearity and the stability of the algorithm are influenced, the conventional DACM algorithm is realized by a large number of nested operations, the operation error and the accumulated noise are introduced, the anti-noise capability of the algorithm is reduced, and the complex algorithm operation occupies too many computer resources and influences the execution efficiency of the algorithm.
Compared with the prior art, the method removes the limit of the value taking domain of the measurable motion signal, does not need tangent and anti-tangent transformation, and directly carries out differential cross multiplication demodulation on the original I/Q signal to obtain the motion signal, thereby greatly reducing the iteration times, the accumulated noise and the approximate error of the algorithm, ensuring the high linearity of the algorithm, improving the operation efficiency of the algorithm and simultaneously showing robustness to the noise.
Drawings
FIG. 1 is a block diagram of the system of the present invention;
FIG. 2 is a flow chart of the method of the present invention;
FIG. 3 is a diagram illustrating simulation results of small-amplitude motion;
in the figure: (a) the I/Q phase track, and a demodulation signal obtained by (b) an arc tangent algorithm and (c) the method; simulation result of large amplitude motion: (d) the I/Q phase track, and a demodulation signal obtained by (e) an arc tangent algorithm and (f) the method;
FIG. 4 is a diagram illustrating the simulation results of motion signals at 25 dB;
in the figure: (a) an I/Q phase locus, and a demodulation signal obtained by processing (b) an arc tangent algorithm, (c) an existing DACM algorithm and (d) the method;
FIG. 5 is a diagram illustrating the results of demodulating a motion signal having an amplitude of about 10mm measured using a 2-4 GHz radar module;
in the figure: (a) I/Q phase locus, (b), (c) is the inverse tangent algorithm and demodulation of this method restores the signal oscillogram separately;
FIG. 6 is a schematic diagram showing the experimental results of measuring about 4mm linear motion using a 120GHz radar module;
in the figure: (a) I/Q phase trajectory, and via (b) an arctangent algorithm, (c) an existing DACM algorithm, and (d) the present method, respectively.
Detailed Description
The invention relates to a linear motion demodulation implementation method based on improved differential cross multiplication of the system, which specifically implements scenes including small-amplitude motion detection, large-amplitude motion detection and motion information demodulation simulation under large noise.
Example 1
Small amplitude motion detection
This embodiment utilizes 2.4 GHz's radar sensor system to detect the mechanical motion of target object within 10mm, based on trigonometric function geometric theorem, does not pass through anti-tangential transformation, directly carries out the cross multiplication to calibration signal and handles in order to realize motion information demodulation, specifically does:
step ① differentiates the calibration signal directly in the time domain to yield:
Figure BDA0002501480860000021
Figure BDA0002501480860000022
step ② utilizes the geometric theorem cos of trigonometric functions2Φ+sin2And (3) directly extracting a motion signal to be measured from the differential result of the step ①, wherein phi is 1:
Figure BDA0002501480860000023
because the differential signal is easily interfered by noise and is not beneficial to the processing of a subsequent circuit, the motion signal to be detected is restored to the motion signal under the digital domain in the next step, and a final digitized algorithm expression is obtained:
Figure BDA0002501480860000031
the results of the experiment are recorded in figure 5. As shown in fig. 5 (a), when the motion signal to be measured is not within the range of the value-taking region of the arctangent algorithm, the waveform of the demodulated signal is subjected to a skip truncation, and a phase ambiguity occurs. In contrast, the motion signal reconstructed by the modified DACM algorithm proposed in fig. 3 (b) is free from the limit of the value range, and has higher algorithm universality and robustness.
Example 2
Large amplitude motion detection
In the embodiment, the measurement of the large-amplitude motion signal is realized by using the 120GHz radar detection module with the wavelength of only 2.5mm, and except for small-amplitude motion which is not in the range of the value taking domain, the large-displacement motion can also cause phase aliasing and cause phase ambiguity. In fig. 3 (a) - (c), the I/Q signal phase trace spans the first and second quadrants and the trace extent is less than pi/2. In this case, the motion signal reconstructed by the arctan algorithm is phase blurred, whereas the demodulated signal obtained by the present method is free from phase distortion. Further, when the I/Q phase trajectory is far beyond 2 π, the arctan demodulation result again exhibits phase discontinuity, as shown in (d) - (f) of FIG. 3. On the contrary, the method can still maintain high linearity, authenticity and integrity of the target motion information.
This embodiment is based on the trigonometric function geometric theorem, and does not pass through inverse tangential transformation, and directly performs cross multiplication processing on the calibration signal to realize motion information demodulation, specifically:
step ① differentiates the calibration signal directly in the time domain to yield:
Figure BDA0002501480860000032
Figure BDA0002501480860000033
step ② utilizes the geometric theorem cos of trigonometric functions2Φ+sin2And (3) directly extracting a motion signal to be measured from the differential result of the step ①, wherein phi is 1:
Figure BDA0002501480860000034
because the differential signal is easily interfered by noise and is not beneficial to the processing of a subsequent circuit, the motion signal to be detected is restored to the motion signal under the digital domain in the next step, and a final digitized algorithm expression is obtained:
Figure BDA0002501480860000035
due to the extremely short wavelength of the large-amplitude motion, the I/Q phase locus can easily overflow by 2 pi by weak displacement, so that a large amount of phase aliasing is caused, and the large-amplitude motion characteristic is achieved. In this experiment, the object to be measured is moved linearly by 4mm in a single direction, and fig. 6 records the I/Q phase locus and the demodulation results obtained by the three algorithms, respectively. As shown in fig. 6 (a), due to the defect of the pre-dc offset calibration algorithm, the I/Q signal after calibration has a mismatch phenomenon, so that part of dc offset remains in the I/Q signal, and the amplitude cannot be completely normalized. This has a significant impact on the current DACM algorithm, which relies heavily on calibration accuracy.
Comparing the demodulation results of the algorithms in fig. 6, it is known that large amplitude motion causes phase ambiguity and waveform distortion of the arctan demodulation result. However, although the existing DACM algorithm does not generate phase ambiguity, due to the imperfection of the pre-calibration algorithm, the demodulation signal still changes with time, and the linearity of the signal is reduced. Compared with the traditional DACM algorithm, the improved DACM algorithm provided by the method not only well solves the problem of phase ambiguity, but also can perfectly restore the motion signal to be detected under large displacement, the Root Mean Square Error (RMSE) of the improved DACM algorithm is improved by 32dB compared with the traditional DACM algorithm, and the linearity of the demodulation algorithm after calibration mismatch is obviously improved.
Example 3
Motion information demodulation simulation under large noise
In the traditional DACM algorithm, Q/I is subjected to inverse tangential transformation as a whole, differentiation and cross multiplication are performed on the basis of the inverse tangential transformation, a large amount of approximation errors and accumulated noise are introduced in the operation process, the anti-noise capability of the algorithm is seriously influenced, and higher requirements are provided for the precision of a pre-calibration step. Compared with the DACM algorithm in the method, the DACM algorithm has a more simplified expression, the operation times and the approximation error are greatly reduced, and the stability and the anti-noise capability of the demodulation process are optimized. Fig. 4 shows the demodulated signal obtained by each algorithm at a signal-to-noise ratio of 25 dB. It can be seen that compared with the existing DACM algorithm, the variance of the arctangent algorithm is minimum, and the algorithm proposed by the method also has higher stability and robustness, and the Root Mean Square Error (RMSE) of the algorithm is close to the arctangent result and is improved by about 9dB compared with the existing DACM algorithm.
In summary, compared with the existing differential cross multiplication algorithm, the present invention firstly has no phase ambiguity, secondly has high linearity under large displacement measurement (the root mean square error is improved by 32dB compared with the traditional DACM algorithm), and simultaneously has high stability and robustness to noise (the root mean square error is improved by 9dB compared with the traditional DACM algorithm under 25dB signal-to-noise ratio). In addition, compared with the existing DACM algorithm, the method does not need tangent and inverse tangent transformation, reduces the iteration times of the algorithm, saves the computing resources and improves the operation efficiency of the algorithm.
The foregoing embodiments may be modified in many different ways by those skilled in the art without departing from the spirit and scope of the invention, which is defined by the appended claims and all changes that come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.

Claims (2)

1. A linear motion demodulation implementation method based on improved differential cross multiplication is characterized in that a motion signal to be detected can be extracted by directly carrying out differential processing on an I/Q signal subjected to direct current calibration based on trigonometric function geometric theorem without inverse tangent transformation, and then carrying out cross multiplication and logical addition operation on the I/Q signal and an original I/Q signal.
2. The linear motion demodulation implementation method based on the improved differential cross multiplication as claimed in claim 1, characterized by specifically:
step ① differentiates the calibration signal directly in the time domain to yield:
Figure FDA0002501480850000011
Figure FDA0002501480850000012
step ② utilizes the geometric theorem cos of trigonometric functions2Φ+sin2And (3) directly extracting a motion signal to be measured from the differential result of the step ①, wherein phi is 1:
Figure FDA0002501480850000013
because the differential signal is easily interfered by noise and is not beneficial to the processing of a subsequent circuit, the motion signal to be detected is restored to the motion signal under the digital domain in the next step, and a final digitized algorithm expression is obtained:
Figure FDA0002501480850000014
CN202010433815.4A 2020-05-21 2020-05-21 Linear motion demodulation implementation method based on improved differential cross multiplication Pending CN111505601A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010433815.4A CN111505601A (en) 2020-05-21 2020-05-21 Linear motion demodulation implementation method based on improved differential cross multiplication

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010433815.4A CN111505601A (en) 2020-05-21 2020-05-21 Linear motion demodulation implementation method based on improved differential cross multiplication

Publications (1)

Publication Number Publication Date
CN111505601A true CN111505601A (en) 2020-08-07

Family

ID=71872098

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010433815.4A Pending CN111505601A (en) 2020-05-21 2020-05-21 Linear motion demodulation implementation method based on improved differential cross multiplication

Country Status (1)

Country Link
CN (1) CN111505601A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112946589A (en) * 2021-02-01 2021-06-11 上海交通大学 Phase self-calibration method for motion measurement of asynchronous FMCW radar system
CN112965035A (en) * 2021-02-01 2021-06-15 上海交通大学 High-linearity phase demodulation implementation method for FMCW radar coherent phase tracking

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102749618A (en) * 2012-07-19 2012-10-24 浙江大学 Method and system for measuring motion trajectory of object by difference and cross-multiplication module
US20150071390A1 (en) * 2013-09-09 2015-03-12 Mstar Semiconductor, Inc. Mixer biasing for intermodulation distortion compensation
CN109965858A (en) * 2019-03-28 2019-07-05 北京邮电大学 Based on ULTRA-WIDEBAND RADAR human body vital sign detection method and device
CN110411486A (en) * 2019-07-26 2019-11-05 浙江理工大学 The PGC-DCDM demodulation method insensitive to phase delay and modulation depth
CN110584631A (en) * 2019-10-10 2019-12-20 重庆邮电大学 Static human heartbeat and respiration signal extraction method based on FMCW radar
CN111091831A (en) * 2020-01-08 2020-05-01 上海交通大学 Silent lip language recognition method and system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102749618A (en) * 2012-07-19 2012-10-24 浙江大学 Method and system for measuring motion trajectory of object by difference and cross-multiplication module
US20150071390A1 (en) * 2013-09-09 2015-03-12 Mstar Semiconductor, Inc. Mixer biasing for intermodulation distortion compensation
CN109965858A (en) * 2019-03-28 2019-07-05 北京邮电大学 Based on ULTRA-WIDEBAND RADAR human body vital sign detection method and device
CN110411486A (en) * 2019-07-26 2019-11-05 浙江理工大学 The PGC-DCDM demodulation method insensitive to phase delay and modulation depth
CN110584631A (en) * 2019-10-10 2019-12-20 重庆邮电大学 Static human heartbeat and respiration signal extraction method based on FMCW radar
CN111091831A (en) * 2020-01-08 2020-05-01 上海交通大学 Silent lip language recognition method and system

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112946589A (en) * 2021-02-01 2021-06-11 上海交通大学 Phase self-calibration method for motion measurement of asynchronous FMCW radar system
CN112965035A (en) * 2021-02-01 2021-06-15 上海交通大学 High-linearity phase demodulation implementation method for FMCW radar coherent phase tracking
CN112946589B (en) * 2021-02-01 2022-09-06 上海交通大学 Phase self-calibration method for motion measurement of asynchronous FMCW radar system

Similar Documents

Publication Publication Date Title
CN111505601A (en) Linear motion demodulation implementation method based on improved differential cross multiplication
CN114910112B (en) Signal error correction method, magnetic encoder, and optical encoder
CN111182705B (en) Time-varying plasma diagnosis method and diagnosis system based on automatic encoder
CN111562438B (en) Sinusoidal signal frequency estimation method and device based on FFT and phase difference
CN110632588B (en) Zero intermediate frequency secondary radar direct current offset compensation algorithm based on FPGA
CN108120888A (en) New cyclotron range stability measuring system
CN105680858B (en) A method of estimation TIADC parallel acquisition system time offset errors
CN111245443A (en) DSADC-based rotary soft decoding processing method and device
US5239273A (en) Digital demodualtor using signal processor to evaluate period measurements
CN113116320A (en) FMCW radar life signal detection method based on VMD
US6794857B2 (en) Apparatus and method for measuring a phase delay characteristic
CN110632563B (en) Intra-pulse frequency coding signal parameter measuring method based on short-time Fourier transform
Fritsch et al. A digital envelope detection filter for real-time operation
US7460979B2 (en) Method and system for enhanced resolution, automatically-calibrated position sensor
US8238506B2 (en) Phase-discriminating device and method
CN113820016B (en) Phase modulation thermal wave signal total variation denoising method
Lee et al. Enhanced pulse amplitude estimation for electronic warfare systems
Huang et al. Frequency estimator based on spectrum correction and remainder sifting for undersampled real-valued waveforms
JP2005140737A (en) Magnetic encoder device
CN107241099B (en) Angle sensor signal processing circuit and processing method
CN117249846B (en) Encoder pre-decoding processing method, system and storage medium
CN114370814A (en) Angle extraction circuit, method and chip
CN114608627B (en) High-precision wide-range phase measurement system based on over-quadrant detection
CN113466859B (en) Spin space debris target ISAR two-dimensional imaging method based on rapid phase interpolation
CN115001345A (en) Sine and cosine encoder subdivision output method and system based on angle interpolation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200807