[go: up one dir, main page]

CN114966656A - Positioning method and device based on millimeter wave equipment - Google Patents

Positioning method and device based on millimeter wave equipment Download PDF

Info

Publication number
CN114966656A
CN114966656A CN202210589994.XA CN202210589994A CN114966656A CN 114966656 A CN114966656 A CN 114966656A CN 202210589994 A CN202210589994 A CN 202210589994A CN 114966656 A CN114966656 A CN 114966656A
Authority
CN
China
Prior art keywords
target
signal
millimeter wave
discrete
frequency band
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
CN202210589994.XA
Other languages
Chinese (zh)
Other versions
CN114966656B (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN202210589994.XA priority Critical patent/CN114966656B/en
Publication of CN114966656A publication Critical patent/CN114966656A/en
Application granted granted Critical
Publication of CN114966656B publication Critical patent/CN114966656B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/42Simultaneous measurement of distance and other co-ordinates
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12Simultaneous equations, e.g. systems of linear equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141Discrete Fourier transforms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Theoretical Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Operations Research (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Probability & Statistics with Applications (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Discrete Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供一种基于毫米波设备的定位方法及装置,方法包括以下步骤:获取待测目标的个数和回波信号,并基于待测目标的个数和回波信号构建频谱信号模型;基于频谱信号模型计算待测目标的采样误差值和相邻天线的相位差;基于采样误差值,计算待测目标到毫米波设备的距离值,并基于相位差,计算待测目标的波达角;根据距离值和波达角,生成待测目标的定位结果。本发明以毫米波设备为基础,利用毫米波波长短、分辨率高的优点,针对不同目标个数的情况,通过建立频谱采样点与信号频率的定量分析模型,对信号频率进行精准计算,从而完成高精度的距离测量;在信号频率的计算过程中提高了信号波达角的测量精度,进而实现高精度的待测目标定位。

Figure 202210589994

The invention provides a positioning method and device based on millimeter wave equipment. The method includes the following steps: acquiring the number of targets to be measured and echo signals, and constructing a spectrum signal model based on the number of targets to be measured and echo signals; The spectral signal model calculates the sampling error value of the target to be measured and the phase difference between adjacent antennas; based on the sampling error value, the distance value from the target to be measured to the millimeter wave device is calculated, and the angle of arrival of the target to be measured is calculated based on the phase difference; According to the distance value and the angle of arrival, the positioning result of the target to be measured is generated. The invention is based on millimeter wave equipment, utilizes the advantages of short wavelength and high resolution of millimeter wave, and according to the situation of different target numbers, by establishing a quantitative analysis model of spectrum sampling points and signal frequency, the signal frequency is accurately calculated, thereby Complete high-precision distance measurement; improve the measurement accuracy of the signal wave arrival angle during the calculation of the signal frequency, thereby achieving high-precision target positioning.

Figure 202210589994

Description

一种基于毫米波设备的定位方法及装置A positioning method and device based on millimeter wave equipment

技术领域technical field

本发明涉及无线感知技术领域,尤其涉及一种基于毫米波设备的定位方法及装置。The present invention relates to the field of wireless sensing technologies, and in particular, to a positioning method and device based on a millimeter wave device.

背景技术Background technique

准确的目标定位技术在日常应用中显得尤为重要,包括安全监控、虚拟现实以及智能家居。同时,一些工业应用也需要高分辨率的定位,比如机械臂、传送带、地铁列车控制系统。结合这些应用,目前诞生了很多基于不同设备的目标定位技术,包括基于可穿戴传感器的定位以及基于视觉信号的定位。然而,这些定位技术都存在着一定的缺陷,基于传感器的定位需要目标携带特定的设备,在安防监控等场景下使用不够方便。基于视觉信号的定位容易对用户的隐私产生威胁,对于家庭场景以及重要的工作单位并不适用;此外,摄像头受光照条件影响较大,在光照条件极端的情况下,将会对测量结果产生很大的干扰。Accurate object location technology is especially important in everyday applications, including security surveillance, virtual reality, and smart homes. At the same time, some industrial applications also require high-resolution positioning, such as robotic arms, conveyor belts, and subway train control systems. Combining these applications, many target positioning technologies based on different devices have been born, including positioning based on wearable sensors and positioning based on visual signals. However, these positioning technologies all have certain defects. Sensor-based positioning requires the target to carry specific equipment, which is not convenient to use in scenarios such as security monitoring. Positioning based on visual signals is easy to threaten the privacy of users, and is not suitable for family scenes and important work units; in addition, the camera is greatly affected by lighting conditions, and in the case of extreme lighting conditions, it will have a great impact on the measurement results. big distraction.

近一段时间以来,基于无线信号的定位技术逐渐受到大家的关注,无线信号主要包括Wi-Fi信号,RFID以及毫米波信号等。近些年来,除了其常见的通信功能,研究者们开始探索利用无线信号的反射、散射等现象来获取其包含的活动信息,并实现目标运动状态的感知。毫米波信号的波长通常在毫米级别。由于毫米波信号的穿透力强、带宽大,具有较高的信号分辨率以及穿透烟雾的能力,受环境影响较小。利用电磁波信号对人员活动进行感知,能够有效的保护人员的隐私;此外,电磁波能够利用人员的反射信号实现活动感知,不需要目标携带特定设备就可实现定位与追踪,使用范围更加广泛。Recently, positioning technology based on wireless signals has gradually attracted everyone's attention. Wireless signals mainly include Wi-Fi signals, RFID and millimeter wave signals. In recent years, in addition to its common communication functions, researchers have begun to explore the use of reflection, scattering and other phenomena of wireless signals to obtain the activity information contained in them, and to realize the perception of the movement state of the target. The wavelength of mmWave signals is usually in the millimeter level. Due to the strong penetrating power, large bandwidth, high signal resolution and ability to penetrate smoke, the millimeter wave signal is less affected by the environment. The use of electromagnetic wave signals to sense personnel activities can effectively protect the privacy of personnel; in addition, electromagnetic waves can use the reflected signals of personnel to realize activity perception, which can achieve positioning and tracking without the need for the target to carry specific equipment, and has a wider range of applications.

然而,由于受到设备的带宽和采样率的影响,目前已有的基于无线信号的定位工作,通常只能达到厘米精度,难以实现针对普适目标的更高精度的定位。在硬件条件受限的情况下,提高信号的分辨率,还面对着如下几个挑战:However, due to the influence of the bandwidth and sampling rate of the device, the existing positioning work based on wireless signals can only achieve centimeter accuracy, and it is difficult to achieve higher-precision positioning for ubiquitous targets. In the case of limited hardware conditions, improving the resolution of the signal also faces the following challenges:

(1)测距的分辨率受限。基于信号处理的分析理论,信号的测距分辨率与带宽成反比关系。目前常见的商用毫米波雷达,带宽通常在4GHz左右,计算得到的距离分辨率通常在4厘米左右。由于硬件的限制,通过增加信号的带宽来提高距离的分辨率并不现实。目前的分辨率还难以达到毫米级别。(1) The resolution of ranging is limited. Based on the analysis theory of signal processing, the ranging resolution of the signal is inversely proportional to the bandwidth. At present, common commercial millimeter-wave radars usually have a bandwidth of about 4 GHz, and the calculated distance resolution is usually about 4 cm. Due to hardware limitations, it is not practical to increase the resolution of the distance by increasing the bandwidth of the signal. The current resolution is still difficult to reach the millimeter level.

(2)角度分辨率受限。已有的信号到达角的计算方法是根据天线阵列中各天线之间的相位差进行计算,容易受到相位噪声的影响。(2) The angular resolution is limited. The existing calculation method of the angle of arrival of the signal is based on the phase difference between the antennas in the antenna array, which is easily affected by phase noise.

(3)多目标的定位。多个目标可能会同时出现在监控区域中,反射信号可能相互混叠,影响目标定位的精度。(3) Multi-target positioning. Multiple targets may appear in the monitoring area at the same time, and the reflected signals may overlap with each other, which affects the accuracy of target positioning.

发明内容SUMMARY OF THE INVENTION

本发明提供一种基于毫米波设备的定位方法及装置,用以解决现有技术中毫米波设备距离分辨率低下、定位精度差的缺陷,实现对待测目标的高精度定位。The invention provides a positioning method and device based on millimeter wave equipment, which are used to solve the defects of low distance resolution and poor positioning accuracy of the millimeter wave equipment in the prior art, and realize high-precision positioning of the target to be measured.

本发明提供一种基于毫米波设备的定位方法,所述方法包括:The present invention provides a positioning method based on a millimeter wave device, the method comprising:

获取待测目标的个数和回波信号,并基于所述待测目标的个数和所述回波信号构建频谱信号模型;Acquire the number of targets to be measured and echo signals, and build a spectrum signal model based on the number of targets to be tested and the echo signals;

基于所述频谱信号模型计算所述待测目标的采样误差值和相邻天线的相位差;Calculate the sampling error value of the target to be measured and the phase difference between adjacent antennas based on the spectral signal model;

基于所述采样误差值,计算所述待测目标到所述毫米波设备的距离值,并基于所述相位差,计算所述待测目标的波达角;Based on the sampling error value, calculate the distance value from the target to be measured to the millimeter wave device, and calculate the angle of arrival of the target to be measured based on the phase difference;

根据所述距离值和所述波达角,生成所述待测目标的定位结果。According to the distance value and the angle of arrival, a positioning result of the target to be detected is generated.

根据本发明提供的一种基于毫米波设备的定位方法,所述基于所述待测目标的个数和所述回波信号构建频谱信号模型,具体包括:According to a positioning method based on a millimeter wave device provided by the present invention, the construction of a spectrum signal model based on the number of the targets to be measured and the echo signals specifically includes:

提取所述回波信号中的目标频段信号,对所述目标频段信号进行离散采样,得到离散目标频段信号;extracting the target frequency band signal in the echo signal, and discretely sampling the target frequency band signal to obtain the discrete target frequency band signal;

对所述离散目标频段信号进行离散傅里叶变换得到处理结果,基于所述处理结果获取所述离散目标频段信号对应的峰值采样点;Performing discrete Fourier transform on the discrete target frequency band signal to obtain a processing result, and obtaining a peak sampling point corresponding to the discrete target frequency band signal based on the processing result;

选取所述峰值采样点预设距离内的若干采样点数据,结合所述处理结果,建立频谱信号模型。Several sampling point data within a preset distance of the peak sampling point are selected, and a spectrum signal model is established in combination with the processing results.

根据本发明提供的一种基于毫米波设备的定位方法,所述待测目标的个数为单个时,所述离散目标频段信号表示为:According to a positioning method based on a millimeter wave device provided by the present invention, when the number of the target to be measured is single, the discrete target frequency band signal is expressed as:

Figure BDA0003664731460000031
Figure BDA0003664731460000031

其中,n为离散目标频段信号的离散时间点,N为信号的采样点数,n=0,1,…,N-1;e为自然对数的底数,j为虚数单位;Among them, n is the discrete time point of the discrete target frequency band signal, N is the number of sampling points of the signal, n=0,1,...,N-1; e is the base of the natural logarithm, and j is the imaginary unit;

Figure BDA0003664731460000032
是信号的幅度,A0和θ0均为实数;
Figure BDA0003664731460000033
为信号的归一化频率,范围为[0,1],δ是采样误差,δ为
Figure BDA0003664731460000034
范围内的实数,π为圆周率,w[n]为高斯白噪声。
Figure BDA0003664731460000032
is the amplitude of the signal, A 0 and θ 0 are real numbers;
Figure BDA0003664731460000033
is the normalized frequency of the signal, in the range [0,1], δ is the sampling error, and δ is
Figure BDA0003664731460000034
A real number in the range, π is pi, and w[n] is white Gaussian noise.

根据本发明提供的一种基于毫米波设备的定位方法,所述对所述离散目标频段信号进行离散傅里叶变换得到处理结果,具体为:According to a positioning method based on a millimeter wave device provided by the present invention, the discrete Fourier transform of the discrete target frequency band signal is performed to obtain a processing result, specifically:

对所述离散目标频段信号的表示r[n]进行离散傅里叶变换,所述处理结果如下:The discrete Fourier transform is performed on the representation r[n] of the discrete target frequency band signal, and the processing result is as follows:

Figure BDA0003664731460000035
Figure BDA0003664731460000035

其中,m=kp-N+1,kp-N+2,…,kp;W[kp-m]是w[n]的离散傅里叶变换结果;kp为峰值采样点。Wherein, m= kp -N+1, kp -N+2,..., kp ; W[ kp -m] is the discrete Fourier transform result of w[n]; kp is the peak sampling point.

根据本发明提供的一种基于毫米波设备的定位方法,所述选取所述峰值采样点预设距离内的若干采样点数据,结合所述处理结果,建立频谱信号模型,具体为:According to a positioning method based on a millimeter wave device provided by the present invention, the selection of several sampling point data within a preset distance of the peak sampling point, and combining the processing results, establish a spectrum signal model, specifically:

选取所述峰值采样点附近的三个点R[kp-1]、R[kp]、R[kp+1],建立频谱信号模型,对模型求解得到所述离散中频信号的复数幅度

Figure BDA0003664731460000048
和归一化频率
Figure BDA0003664731460000041
Select three points R[k p -1], R[k p ], R[k p +1] near the peak sampling point, establish a spectrum signal model, and solve the model to obtain the complex amplitude of the discrete intermediate frequency signal
Figure BDA0003664731460000048
and normalized frequency
Figure BDA0003664731460000041

其中,基于所述归一化频率计算所述待测目标的距离:Wherein, the distance of the target to be measured is calculated based on the normalized frequency:

Figure BDA0003664731460000042
Fs是所述毫米波设备的采样率;
Figure BDA0003664731460000042
F s is the sampling rate of the millimeter wave device;

当所述待测目标的个数为两个时,所述离散目标频段信号的表示为:

Figure BDA0003664731460000043
When the number of the target to be measured is two, the representation of the discrete target frequency band signal is:
Figure BDA0003664731460000043

其中,n为离散目标频段信号的离散时间点,N为信号的采样点数,n=0,1,…,N-1;

Figure BDA0003664731460000044
Figure BDA0003664731460000045
是未知的复数信号幅度,
Figure BDA0003664731460000046
Figure BDA0003664731460000047
是信号的归一化频率分量,w[n]为高斯白噪声;Among them, n is the discrete time point of the discrete target frequency band signal, N is the number of sampling points of the signal, n=0,1,...,N-1;
Figure BDA0003664731460000044
and
Figure BDA0003664731460000045
is the unknown complex signal amplitude,
Figure BDA0003664731460000046
and
Figure BDA0003664731460000047
is the normalized frequency component of the signal, w[n] is Gaussian white noise;

对所述离散目标频段信号t[n]进行离散傅里叶变换,得到运算结果T[k];根据运算结果得到两个峰值采样点kp1、kp2Discrete Fourier transform is performed on the discrete target frequency band signal t[n] to obtain an operation result T[k]; two peak sampling points k p1 and k p2 are obtained according to the operation result;

在所述两个峰值采样点kp1、kp2间选取多个采样点数据,结合所述运算结果T[k]建立所述频谱信号模型,通过模型对未知参数A1,A21212求解;Select a plurality of sampling point data between the two peak sampling points k p1 and k p2 , establish the spectral signal model in combination with the operation result T[k], and use the model to determine the unknown parameters A 1 , A 2 , θ 1 , θ 2 , δ 1 , δ 2 to solve;

基于两个归一化频率f1、f2分别计算所述待测目标的距离。Based on the two normalized frequencies f 1 and f 2 , the distances of the objects to be measured are calculated respectively.

根据本发明提供的一种基于毫米波设备的定位方法,基于所述频谱信号模型计算相邻天线的相位差,具体包括:According to a positioning method based on a millimeter wave device provided by the present invention, the phase difference between adjacent antennas is calculated based on the spectrum signal model, which specifically includes:

获取所述离散目标频段信号模型中的复数信号幅度;obtaining the complex signal amplitude in the discrete target frequency band signal model;

根据所述复数信号幅度,提取所述毫米波设备中相邻的接收天线的相位差。According to the complex signal amplitude, the phase difference between adjacent receiving antennas in the millimeter wave device is extracted.

本发明还提供一种基于毫米波设备的定位装置,所述装置包括:The present invention also provides a positioning device based on a millimeter wave device, the device comprising:

模型构建模块,用于获取待测目标的个数和回波信号,并基于所述待测目标的个数和所述回波信号构建频谱信号模型;a model building module for acquiring the number of targets to be tested and echo signals, and building a spectrum signal model based on the number of targets to be tested and the echo signals;

第一计算模块,用于基于所述频谱信号模型计算所述待测目标的采样误差值和相邻天线的相位差;a first calculation module, configured to calculate the sampling error value of the target to be measured and the phase difference of adjacent antennas based on the spectrum signal model;

第二计算模块,用于基于所述采样误差值,计算所述待测目标到所述毫米波设备的距离值,并基于所述相位差,计算所述待测目标的波达角;a second calculation module, configured to calculate the distance value from the target to be measured to the millimeter wave device based on the sampling error value, and calculate the angle of arrival of the target to be measured based on the phase difference;

结果生成模块,用于根据所述距离值和所述波达角,生成所述待测目标的定位结果。A result generating module, configured to generate a positioning result of the target to be measured according to the distance value and the angle of arrival.

本发明还提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述任一种所述基于毫米波设备的定位方法。The present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and running on the processor, when the processor executes the program, the processor implements the millimeter wave-based method described in any of the above The positioning method of the device.

本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如上述任一种所述基于毫米波设备的定位方法。The present invention also provides a non-transitory computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, implements the positioning method based on any one of the above-mentioned millimeter wave devices.

本发明还提供一种计算机程序产品,包括计算机程序,所述计算机程序被处理器执行时实现如上述任一种所述基于毫米波设备的定位方法。The present invention also provides a computer program product, including a computer program, which, when executed by a processor, implements the positioning method based on any one of the above-mentioned millimeter wave devices.

本发明提供的基于毫米波设备的定位方法及装置,本发明以毫米波设备为基础,利用毫米波波长短、分辨率高的优点,针对不同目标个数的情况,通过建立频谱采样点与信号频率的定量分析模型,对信号频率进行精准计算,从而完成高精度的距离测量;在信号频率的计算过程中提高了信号波达角的测量精度,进而实现高精度的目标定位。The positioning method and device based on millimeter-wave equipment provided by the present invention are based on millimeter-wave equipment and take advantage of the advantages of short wavelength and high resolution of millimeter-wave. The quantitative analysis model of frequency can accurately calculate the signal frequency, so as to complete the high-precision distance measurement; in the process of calculating the signal frequency, the measurement accuracy of the signal arrival angle is improved, thereby achieving high-precision target positioning.

附图说明Description of drawings

为了更清楚地说明本发明或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to explain the present invention or the technical solutions in the prior art more clearly, the following will briefly introduce the accompanying drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are the For some embodiments of the invention, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without any creative effort.

图1是本发明提供的基于毫米波设备的定位方法的流程示意图之一;Fig. 1 is one of the schematic flow charts of the positioning method based on millimeter wave equipment provided by the present invention;

图2是本发明提供的基于毫米波设备的定位方法的流程示意图之二;2 is the second schematic flowchart of the positioning method based on a millimeter wave device provided by the present invention;

图3是本发明提供的基于毫米波设备的定位方法的流程示意图之三;3 is a third schematic flowchart of a positioning method based on a millimeter wave device provided by the present invention;

图4是本发明提供的基于毫米波设备的定位方法的流程示意图之四;4 is a fourth schematic flowchart of a positioning method based on a millimeter wave device provided by the present invention;

图5是本发明提供的频谱采样误差示意图;Fig. 5 is the spectrum sampling error schematic diagram provided by the present invention;

图6是本发明提供的使用不同方法的定位误差的CDF图;6 is a CDF diagram of the positioning error using different methods provided by the present invention;

图7是本发明提供的使用不同方法的AoA计算误差的CDF图;Fig. 7 is the CDF diagram of AoA calculation error using different methods provided by the present invention;

图8是本发明提供的基于毫米波设备的定位装置的结构示意图;8 is a schematic structural diagram of a positioning device based on a millimeter wave device provided by the present invention;

图9是本发明提供的电子设备的结构示意图。FIG. 9 is a schematic structural diagram of an electronic device provided by the present invention.

具体实施方式Detailed ways

为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions in the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are part of the embodiments of the present invention. , not all examples. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

下面结合图1-图4描述本发明的基于毫米波设备的定位方法。The positioning method based on the millimeter wave device of the present invention will be described below with reference to FIG. 1 to FIG. 4 .

图1是本发明实施例提供的基于毫米雷达波设备的定位方法的流程示意图之一。FIG. 1 is one of the schematic flowcharts of a positioning method based on a millimeter radar wave device provided by an embodiment of the present invention.

如图1、4所示,本发明实施例提供的一种基于毫米波设备的定位方法,包括以下步骤:As shown in FIGS. 1 and 4 , a positioning method based on a millimeter wave device provided by an embodiment of the present invention includes the following steps:

步骤110、获取待测目标的个数和回波信号,并基于待测目标的个数和回波信号构建频谱信号模型。具体的,毫米波设备向待测目标发射探测信号,然后接收待测目标反射回来的回波信号。在实际的使用过程中,毫米波雷达通过发射线性调频信号,并接收目标的反射信号来实现高精度的目标定位,毫米波设备从回波信号中提取中频信号,再对中频信号进行离散采样得到中频离散信号。根据理论分析,毫米波设备读取的中频离散信号的频率与目标到设备的距离成正比关系,因此只要得到中频信号的频率,就能够计算出目标到设备的距离,从而实现目标的定位。Step 110: Acquire the number of targets to be tested and echo signals, and build a spectrum signal model based on the number of targets to be tested and echo signals. Specifically, the millimeter wave device transmits a detection signal to the target to be measured, and then receives the echo signal reflected from the target to be measured. In actual use, the millimeter-wave radar achieves high-precision target positioning by transmitting a chirp signal and receiving the reflected signal of the target. The millimeter-wave device extracts the intermediate frequency signal from the echo signal, and then discretely samples the intermediate frequency signal to obtain IF discrete signal. According to theoretical analysis, the frequency of the intermediate frequency discrete signal read by the millimeter wave device is proportional to the distance from the target to the device. Therefore, as long as the frequency of the intermediate frequency signal is obtained, the distance from the target to the device can be calculated, thereby realizing the positioning of the target.

本发明中,针对不同数量的探测目标,毫米波设备对回波信号进行不同的处理,构建用于计算采样误差和相位差的频谱信号模型。In the present invention, for different numbers of detection targets, the millimeter wave device performs different processing on the echo signals, and constructs a spectrum signal model for calculating the sampling error and the phase difference.

步骤120、基于频谱信号模型计算待测目标的采样误差值和相邻天线的相位差。Step 120: Calculate the sampling error value of the target to be measured and the phase difference between adjacent antennas based on the spectral signal model.

常见的提取信号的频率的方法是对目标频段信号进行离散傅里叶变换,从而得到信号的频谱,频谱中的峰值代表信号所包含的频率分量。这样,FMCW信号的频谱分辨率为:A common method for extracting the frequency of a signal is to perform discrete Fourier transform on the target frequency band signal to obtain the frequency spectrum of the signal, and the peaks in the frequency spectrum represent the frequency components contained in the signal. Thus, the spectral resolution of the FMCW signal is:

Δd=c/2B,Δd=c/2B,

其中B为信号的带宽。所以毫米波设备的距离分辨率与带宽成反比。where B is the bandwidth of the signal. So the range resolution of mmWave devices is inversely proportional to the bandwidth.

目前常见的商用毫米波设备的最大带宽为4GHz。考虑到实际应用场景中的带宽损耗问题,距离分辨率通常约为4cm,这样毫米波设备并不能实现较高精度的定位。The maximum bandwidth of common commercial mmWave devices is 4GHz. Considering the problem of bandwidth loss in practical application scenarios, the distance resolution is usually about 4cm, so millimeter-wave devices cannot achieve high-precision positioning.

产生测距误差的一个重要原因是频谱的采样误差。由于频率采样的原因,观察到的信号频率与实际的信号频率之间存在着差别,这样产生了频率误差δ,使得测距结果存在着误差。An important reason for the ranging error is the sampling error of the spectrum. Due to frequency sampling, there is a difference between the observed signal frequency and the actual signal frequency, which results in a frequency error δ, which causes an error in the ranging result.

本发明计算了现有技术中忽略的采样误差。The present invention accounts for sampling errors that are ignored in the prior art.

步骤130、基于采样误差值,计算待测目标到毫米波设备的距离值,并基于相位差,计算待测目标的波达角。Step 130: Calculate the distance value from the target to be measured to the millimeter wave device based on the sampling error value, and calculate the angle of arrival of the target to be measured based on the phase difference.

对中频离散信号的频率估计是实现目标准确测距的关键,本发明通过采样误差值,得到更为精确的频率计算值,从而得到更为精确的距离值,该距离值指的是待测目标到毫米波设备的距离。The frequency estimation of the intermediate frequency discrete signal is the key to realizing the accurate ranging of the target. The present invention obtains a more accurate frequency calculation value by sampling the error value, thereby obtaining a more accurate distance value, and the distance value refers to the target to be measured. Distance to mmWave device.

步骤140、根据距离值和波达角,生成待测目标的定位结果。Step 140: Generate a positioning result of the target to be measured according to the distance value and the angle of arrival.

进一步的,如图2所示,基于待测目标的个数和回波信号构建频谱信号模型,具体包括以下步骤:Further, as shown in Figure 2, building a spectrum signal model based on the number of targets to be measured and echo signals, specifically including the following steps:

步骤210、提取回波信号中的目标频段信号,对目标频段信号进行离散采样,得到离散目标频段信号。本发明中频离散信号即为目标频段信号。Step 210: Extract the target frequency band signal in the echo signal, and perform discrete sampling on the target frequency band signal to obtain the discrete target frequency band signal. The intermediate frequency discrete signal of the present invention is the target frequency band signal.

具体的,当环境中仅存在着单个待测目标时,理论上毫米波设备提取的中频信号中仅存在着一个频率分量,离散目标频段信号可以表示为:Specifically, when there is only a single target to be measured in the environment, theoretically, there is only one frequency component in the intermediate frequency signal extracted by the millimeter wave device, and the discrete target frequency band signal can be expressed as:

Figure BDA0003664731460000081
Figure BDA0003664731460000081

其中,n为离散目标频段信号的离散时间点,N为信号的采样点数,n=0,1,…,N-1;e为自然对数的底数,j为虚数单位;Among them, n is the discrete time point of the discrete target frequency band signal, N is the number of sampling points of the signal, n=0,1,...,N-1; e is the base of the natural logarithm, and j is the imaginary unit;

Figure BDA0003664731460000082
是信号的幅度,A0和θ0均为实数;
Figure BDA0003664731460000083
为信号的归一化频率,范围为[0,1],δ是采样误差,δ为
Figure BDA0003664731460000084
范围内的实数,π为圆周率,w[n]为高斯白噪声。
Figure BDA0003664731460000082
is the amplitude of the signal, A 0 and θ 0 are real numbers;
Figure BDA0003664731460000083
is the normalized frequency of the signal, in the range [0,1], δ is the sampling error, and δ is
Figure BDA0003664731460000084
A real number in the range, π is pi, and w[n] is white Gaussian noise.

如图5、6所示,传统的频率估计方案通常只估计kp而忽视了δ。在本发明中,通过设计频谱峰值重构算法,得到更为准确的频率估计结果。As shown in Figures 5 and 6, traditional frequency estimation schemes usually only estimate k p and ignore δ. In the present invention, a more accurate frequency estimation result is obtained by designing a spectrum peak reconstruction algorithm.

步骤220、对离散目标频段信号进行离散傅里叶变换得到处理结果,基于处理结果获取离散目标频段信号对应的峰值采样点。Step 220: Perform discrete Fourier transform on the discrete target frequency band signal to obtain a processing result, and obtain peak sampling points corresponding to the discrete target frequency band signal based on the processing result.

具体的,通过对r[n]进行DFT,能够得到频谱中的峰值采样点kp。在原信号r[n]中,还有三个实数未知量:A0、θ0和δ。r[n]的DFT运算结果如下:Specifically, by performing DFT on r[n], the peak sampling point k p in the spectrum can be obtained. In the original signal r[n], there are also three real unknowns: A 0 , θ 0 and δ. The DFT operation result of r[n] is as follows:

Figure BDA0003664731460000091
Figure BDA0003664731460000091

其中,m=kp-N+1,kp-N+2,…,kp;W[kp-m]是w[n]的离散傅里叶变换结果;kp为峰值采样点。Wherein, m= kp -N+1, kp -N+2,..., kp ; W[ kp -m] is the discrete Fourier transform result of w[n]; kp is the peak sampling point.

步骤230、选取峰值采样点预设距离内的若干采样点数据,结合处理结果,建立频谱信号模型。Step 230: Select data of several sampling points within a preset distance of the peak sampling point, and combine the processing results to establish a spectrum signal model.

进一步的,选取峰值采样点附近的三个点:Further, select three points near the peak sampling point:

R[kp-1]、R[kp]、R[kp+1],建立频谱信号模型,对模型求解得到离散中频信号的复数幅度

Figure BDA0003664731460000092
和归一化频率
Figure BDA0003664731460000093
R[k p -1], R[k p ], R[k p +1], establish a spectral signal model, and solve the model to obtain the complex amplitude of the discrete intermediate frequency signal
Figure BDA0003664731460000092
and normalized frequency
Figure BDA0003664731460000093

其中,基于归一化频率计算待测目标的距离:Among them, the distance of the target to be measured is calculated based on the normalized frequency:

Figure BDA0003664731460000094
Fs是毫米波设备的采样率。
Figure BDA0003664731460000094
F s is the sampling rate of mmWave devices.

具体的,由于DFT的运算结果是复数,所以上述公式两边的实部和虚部应该分别相等。因此,每个DFT采样点能够提供2个实数方程。Specifically, since the operation result of the DFT is a complex number, the real part and the imaginary part of both sides of the above formula should be equal respectively. Therefore, each DFT sample point can provide 2 real equations.

为了求解3个未知实数,至少需要2个采样点的数据。理论上,使用更多的采样点来进行估计,将得到更为精确的估计结果。然而,幅度较低的采样点容易受到噪声的影响,误差较大。根据频谱计算公式,距离峰值点kp较近的采样点通常有着更高的幅度。因此,使用频谱峰值点附近的三个点R[kp-1],R[kp],R[kp+1]来建立频谱重构非线性方程组,即本发明中的频谱信号模型。To solve for 3 unknown real numbers, at least 2 sample points of data are required. In theory, using more sampling points for estimation will result in more accurate estimation results. However, sampling points with lower amplitudes are susceptible to noise and have larger errors. According to the spectrum calculation formula, the sampling points closer to the peak point k p usually have higher amplitudes. Therefore, three points R[k p -1], R[k p ], R[k p +1] near the spectral peak point are used to establish the spectral reconstruction nonlinear equation system, that is, the spectral signal model in the present invention .

最后,使用常见的非线性方程组求解算法Levenberg-Marquardt对未知参数进行求解,得到信号的复数幅度

Figure BDA0003664731460000101
和归一化频率
Figure BDA0003664731460000102
然后,根据信号的频率计算目标的距离:Finally, use the common nonlinear equation solving algorithm Levenberg-Marquardt to solve the unknown parameters to obtain the complex amplitude of the signal
Figure BDA0003664731460000101
and normalized frequency
Figure BDA0003664731460000102
Then, calculate the distance to the target based on the frequency of the signal:

Figure BDA0003664731460000103
Figure BDA0003664731460000103

其中Fs是毫米波设备的采样率。where F s is the sampling rate of the mmWave device.

与传统的测距方法相比较,获得了更为准确的频率信息,也就得到了更高精度的距离信息。Compared with the traditional ranging method, more accurate frequency information is obtained, and thus higher precision distance information is obtained.

进一步的,本发明对待测目标的个数进行判定,根据个数为单个或多个,分别对目标频段信号进行处理。Further, the present invention determines the number of targets to be measured, and processes the target frequency band signals respectively according to whether the number is single or multiple.

当待测目标的个数为两个时,离散目标频段信号的表示为:When the number of targets to be measured is two, the representation of the discrete target frequency band signal is:

Figure BDA0003664731460000104
Figure BDA0003664731460000104

其中,n为离散目标频段信号的离散时间点,N为信号的采样点数,n=0,1,…,N-1;

Figure BDA0003664731460000105
Figure BDA0003664731460000106
是未知的复数信号幅度,
Figure BDA0003664731460000107
Figure BDA0003664731460000108
是信号的归一化频率分量,w[n]为高斯白噪声;Among them, n is the discrete time point of the discrete target frequency band signal, N is the number of sampling points of the signal, n=0,1,...,N-1;
Figure BDA0003664731460000105
and
Figure BDA0003664731460000106
is the unknown complex signal amplitude,
Figure BDA0003664731460000107
and
Figure BDA0003664731460000108
is the normalized frequency component of the signal, w[n] is Gaussian white noise;

对离散目标频段信号t[n]进行离散傅里叶变换,得到运算结果T[k];根据运算结果得到两个峰值采样点kp1、kp2Discrete Fourier transform is performed on the discrete target frequency band signal t[n] to obtain the operation result T[k]; two peak sampling points k p1 and k p2 are obtained according to the operation result;

在两个峰值采样点kp1、kp2间选取多个采样点数据,结合运算结果T[k]建立所述频谱信号模型,通过模型对未知参数A1,A21212求解;Select a plurality of sampling point data between the two peak sampling points k p1 and k p2 , and combine the operation result T[ k ] to establish the spectrum signal model, and the unknown parameters A1, A2, θ1 , θ2 , δ 12 solve;

基于两个归一化频率f1、f2分别计算待测目标的距离。Based on the two normalized frequencies f 1 and f 2 , the distances of the objects to be measured are calculated respectively.

实际应用场景中,多个目标可能会同时出现,信号中将会同时存在多个频率分量。不同的频率分量之间可能发生混叠,影响各分量参数的估计。为了应对这些复杂情况,提出了多频率分量估计方案。In practical application scenarios, multiple targets may appear at the same time, and there will be multiple frequency components in the signal at the same time. Aliasing may occur between different frequency components, affecting the estimation of the parameters of each component. To deal with these complex situations, a multi-frequency component estimation scheme is proposed.

为了便于理解,首先考虑具有两个频率分量的信号:For ease of understanding, first consider a signal with two frequency components:

Figure BDA0003664731460000111
Figure BDA0003664731460000111

其中,

Figure BDA0003664731460000112
Figure BDA0003664731460000113
是未知的复数信号幅度,
Figure BDA0003664731460000114
Figure BDA0003664731460000115
是信号的归一化频率分量,N为信号的采样点数,w[n]为高斯白噪声。in,
Figure BDA0003664731460000112
and
Figure BDA0003664731460000113
is the unknown complex signal amplitude,
Figure BDA0003664731460000114
and
Figure BDA0003664731460000115
is the normalized frequency component of the signal, N is the number of sampling points of the signal, and w[n] is Gaussian white noise.

根据单目标距离测量的推导内容,计算了t[n]的DFT运算结果T[k],得到信号的频率分量kp1,kp2。在两频率分量信号中,共有6个未知参数:A1,A21212According to the derivation content of single-target distance measurement, the DFT operation result T[k] of t[n] is calculated, and the frequency components k p1 , k p2 of the signal are obtained. In the two-frequency component signal, there are 6 unknown parameters: A 1 , A 2 , θ 1 , θ 2 , δ 1 , δ 2 .

在单频率分量信号中,使用3个DFT采样点的信息来计算未知参数。在两频率分量信号中,不妨设kp1≤kp2,那么使用T[kp1-1]到T[kp2+1]的采样信息对未知参数进行求解。这部分采样点个数不小于3个,足以求解两个分量的6个未知参数。In a single frequency component signal, the information of the 3 DFT sample points is used to calculate the unknown parameters. In the two-frequency component signal, it is advisable to set k p1 ≤k p2 , then use the sampling information from T[k p1 -1] to T[k p2 +1] to solve the unknown parameters. The number of sampling points in this part is not less than 3, which is enough to solve the 6 unknown parameters of the two components.

如果信号中包含有更多的频率分量,使用上述方法,在一个非线性方程组中对其进行求解。为了求解包含有u个频率分量,即3u个未知参数的方程组,需要至少3u个非线性方程,也就是

Figure BDA0003664731460000116
个采样点的信息。If the signal contains more frequency components, solve it in a nonlinear system of equations using the above method. In order to solve a system of equations containing u frequency components, i.e. 3u unknown parameters, at least 3u nonlinear equations are required, i.e.
Figure BDA0003664731460000116
information about sampling points.

如图3所示,在步骤120中,基于频谱信号模型计算相邻天线的相位差,具体包括以下步骤:As shown in Figure 3, in step 120, the phase difference between adjacent antennas is calculated based on the spectrum signal model, which specifically includes the following steps:

步骤310、获取离散目标频段信号模型中的复数信号幅度;Step 310, acquiring the complex signal amplitude in the discrete target frequency band signal model;

步骤320、根据复数信号幅度,提取毫米波设备中相邻的线性接收天线的相位差。Step 320: Extract the phase difference of adjacent linear receiving antennas in the millimeter wave device according to the complex signal amplitude.

具体的,为了获取目标的二维坐标,还需要知道信号的波达角(AoA)。利用毫米波设备的线性接收天线阵列来求解AoA。Specifically, in order to obtain the two-dimensional coordinates of the target, it is also necessary to know the angle of arrival (AoA) of the signal. AoA is solved using linear receive antenna arrays of mmWave devices.

线性阵列中天线之间的间距为d0。考虑到目标到设备的距离d>>d0,反射信号到达各天线时可以近似认为平行。相邻天线的相位差可以表示成:The spacing between antennas in a linear array is do. Considering the distance d>>d 0 from the target to the device, when the reflected signal reaches each antenna, it can be approximately considered to be parallel. The phase difference of adjacent antennas can be expressed as:

Figure BDA0003664731460000121
Figure BDA0003664731460000121

其中λ为信号的波长,θ为信号的AoA。根据上式,可以根据不同天线的相位差来求解AoA。where λ is the wavelength of the signal and θ is the AoA of the signal. According to the above formula, AoA can be solved according to the phase difference of different antennas.

如图3、7所示,实际应用中,由于随机噪声和多目标反射的影响,直接利用相位差求解AoA通常有较大的误差。观察到可以通过DFT来提取各频率分量的初始相位。在单目标的信号建模中,忽略频率偏差δ和高斯白噪声w[n]的影响,其DFT采样点的相位如下:As shown in Figures 3 and 7, in practical applications, due to the influence of random noise and multi-target reflection, there is usually a large error in solving AoA directly by using the phase difference. It is observed that the initial phase of each frequency component can be extracted by DFT. In the single-target signal modeling, ignoring the influence of frequency deviation δ and Gaussian white noise w[n], the phase of the DFT sampling point is as follows:

Figure BDA0003664731460000122
Figure BDA0003664731460000122

其中kp为DFT的峰值点,θ0为对应频率分量的初始相位。根据上式,如果忽略频率偏差δ的影响,那么可以根据DFT的采样点信息提取出各天线的初始相位。where k p is the peak point of the DFT, and θ 0 is the initial phase of the corresponding frequency component. According to the above formula, if the influence of the frequency deviation δ is ignored, the initial phase of each antenna can be extracted according to the sampling point information of the DFT.

但是在大多数情况下,频率偏差δ≠0。这时DFT频谱中峰值点对应的相位如下:But in most cases, the frequency deviation δ≠0. At this time, the phase corresponding to the peak point in the DFT spectrum is as follows:

angle(R[kp])=θ0+angle(F(δ)),angle(R[k p ])=θ 0 +angle(F(δ)),

其中

Figure BDA0003664731460000123
in
Figure BDA0003664731460000123

由于受到F(δ)的影响,DFT峰值点的相位不再对应天线的初始相位,因此使用DFT峰值点的相位计算AoA将会产生较大的系统误差。Due to the influence of F(δ), the phase of the DFT peak point no longer corresponds to the initial phase of the antenna, so using the phase of the DFT peak point to calculate the AoA will generate a large systematic error.

为了减少系统误差的产生,提高AoA计算的精度,需要准确地计算出各频率分量对应的初始相位θ0。在前面的介绍中,在求解出各频率分量的参数δ的同时,也获取了其复数幅度

Figure BDA0003664731460000124
因此,各频率分量的初始相位已经被所解出。因此,可以直接利用频谱重建所获得的参数信息,直接得到各频率分量的初始相位,并利用各接收天线的求解结果,得到更为准确的AoA,实现了AoA的去畸变测量。In order to reduce the generation of systematic errors and improve the accuracy of AoA calculation, it is necessary to accurately calculate the initial phase θ 0 corresponding to each frequency component. In the previous introduction, while solving the parameter δ of each frequency component, its complex amplitude is also obtained.
Figure BDA0003664731460000124
Therefore, the initial phase of each frequency component has been solved. Therefore, the parameter information obtained by spectrum reconstruction can be directly used to directly obtain the initial phase of each frequency component, and the solution result of each receiving antenna can be used to obtain a more accurate AoA, which realizes the AoA de-distortion measurement.

下面对本发明提供的基于毫米波设备的定位装置进行描述,下文描述的基于毫米波设备的定位装置与上文描述的基于毫米波设备的定位方法可相互对应参照。The positioning apparatus based on a millimeter wave device provided by the present invention is described below, and the positioning apparatus based on a millimeter wave device described below and the positioning method based on a millimeter wave device described above may refer to each other correspondingly.

本发明以毫米波设备为基础,利用毫米波波长短、分辨率高的优点,针对不同目标个数的情况,通过建立频谱采样点与信号频率的定量分析模型,设计实现了频谱峰值重构算法,从而完成高精度的距离测量;并通过信号的幅度信息提高信号波达角的测量精度,进而实现高精度的目标定位系统。The invention is based on millimeter wave equipment, utilizes the advantages of short wavelength and high resolution of millimeter wave, and designs and realizes a spectrum peak reconstruction algorithm by establishing a quantitative analysis model of spectrum sampling points and signal frequencies for different target numbers. , so as to complete the high-precision distance measurement; and improve the measurement accuracy of the signal arrival angle through the signal amplitude information, thereby realizing a high-precision target positioning system.

如图8所示,本发明实施例提供的一种基于毫米波设备的定位装置,该装置包括以下模块:模型构建模块810、第一计算模块820、第二计算模块830及结果生成模块840。As shown in FIG. 8 , a positioning apparatus based on a millimeter wave device provided by an embodiment of the present invention includes the following modules: a model building module 810 , a first calculation module 820 , a second calculation module 830 , and a result generation module 840 .

具体的,模型构建模块810用于获取待测目标的个数和回波信号,并基于待测目标的个数和回波信号构建频谱信号模型。第一计算模块820用于基于频谱信号模型计算待测目标的采样误差值和相邻天线的相位差。第二计算模块830用于基于采样误差值,计算待测目标到毫米波设备的距离值,并基于相位差,计算待测目标的波达角。结果生成模块840用于根据距离值和波达角,生成待测目标的定位结果。Specifically, the model building module 810 is configured to acquire the number of the objects to be tested and the echo signals, and build a spectrum signal model based on the number of the objects to be tested and the echo signals. The first calculation module 820 is configured to calculate the sampling error value of the target to be measured and the phase difference between adjacent antennas based on the spectral signal model. The second calculation module 830 is configured to calculate the distance value from the target to be measured to the millimeter wave device based on the sampling error value, and calculate the angle of arrival of the target to be measured based on the phase difference. The result generating module 840 is configured to generate the positioning result of the target to be measured according to the distance value and the angle of arrival.

图9示例了一种电子设备的实体结构示意图,如图9所示,该电子设备可以包括:处理器(processor)910、通信接口(Communications Interface)920、存储器(memory)930和通信总线940,其中,处理器910,通信接口920,存储器930通过通信总线940完成相互间的通信。处理器910可以调用存储器930中的逻辑指令,以执行一种基于毫米波设备的定位方法,该方法包括以下步骤:获取待测目标的个数和回波信号,并基于待测目标的个数和回波信号构建频谱信号模型;基于频谱信号模型计算待测目标的采样误差值和相邻天线的相位差;基于采样误差值,计算待测目标到毫米波设备的距离值,并基于相位差,计算待测目标的波达角;根据距离值和波达角,生成待测目标的定位结果。FIG. 9 illustrates a schematic diagram of the physical structure of an electronic device. As shown in FIG. 9 , the electronic device may include: a processor (processor) 910, a communication interface (Communications Interface) 920, a memory (memory) 930, and a communication bus 940, The processor 910 , the communication interface 920 , and the memory 930 communicate with each other through the communication bus 940 . The processor 910 can call the logic instructions in the memory 930 to execute a positioning method based on a millimeter wave device, the method includes the following steps: obtaining the number of the objects to be measured and the echo signals, and based on the number of the objects to be measured and echo signals to construct a spectrum signal model; based on the spectrum signal model, calculate the sampling error value of the target to be tested and the phase difference between adjacent antennas; , calculate the angle of arrival of the target to be measured; generate the positioning result of the target to be measured according to the distance value and the angle of arrival.

此外,上述的存储器930中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。In addition, the above-mentioned logic instructions in the memory 930 can be implemented in the form of software functional units and can be stored in a computer-readable storage medium when sold or used as an independent product. Based on this understanding, the technical solution of the present invention can be embodied in the form of a software product in essence, or the part that contributes to the prior art or the part of the technical solution. The computer software product is stored in a storage medium, including Several instructions are used to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, Read-Only Memory (ROM, Read-Only Memory), Random Access Memory (RAM, Random Access Memory), magnetic disk or optical disk and other media that can store program codes .

另一方面,本发明还提供一种计算机程序产品,所述计算机程序产品包括计算机程序,计算机程序可存储在非暂态计算机可读存储介质上,所述计算机程序被处理器执行时,计算机能够执行上述各方法所提供的一种基于毫米波设备的定位方法,该方法包括以下步骤:获取待测目标的个数和回波信号,并基于待测目标的个数和回波信号构建频谱信号模型;基于频谱信号模型计算待测目标的采样误差值和相邻天线的相位差;基于采样误差值,计算待测目标到毫米波设备的距离值,并基于相位差,计算待测目标的波达角;根据距离值和波达角,生成待测目标的定位结果。In another aspect, the present invention also provides a computer program product, the computer program product includes a computer program, the computer program can be stored on a non-transitory computer-readable storage medium, and when the computer program is executed by a processor, the computer can Execute a positioning method based on a millimeter wave device provided by the above methods, the method includes the following steps: acquiring the number of targets to be measured and echo signals, and based on the number of targets to be tested and echo signals Constructing a spectrum signal Model; calculate the sampling error value of the target to be tested and the phase difference between adjacent antennas based on the spectrum signal model; Arrival angle: According to the distance value and the arrival angle, the positioning result of the target to be measured is generated.

又一方面,本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各方法提供的一种基于毫米波设备的定位方法,该方法包括以下步骤:获取待测目标的个数和回波信号,并基于待测目标的个数和回波信号构建频谱信号模型;基于频谱信号模型计算待测目标的采样误差值和相邻天线的相位差;基于采样误差值,计算待测目标到毫米波设备的距离值,并基于相位差,计算待测目标的波达角;根据距离值和波达角,生成待测目标的定位结果。In yet another aspect, the present invention also provides a non-transitory computer-readable storage medium on which a computer program is stored, and the computer program is implemented when executed by a processor to perform a millimeter-wave device-based positioning provided by the above methods The method includes the following steps: obtaining the number of the objects to be measured and echo signals, and constructing a spectrum signal model based on the number of objects to be tested and the echo signals; calculating the sampling error value and Phase difference between adjacent antennas; based on the sampling error value, calculate the distance value from the target to be measured to the millimeter wave device, and calculate the angle of arrival of the target to be measured based on the phase difference; generate the target to be measured according to the distance value and the angle of arrival positioning results.

以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。The device embodiments described above are only illustrative, wherein the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in One place, or it can be distributed over multiple network elements. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution in this embodiment. Those of ordinary skill in the art can understand and implement it without creative effort.

通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。From the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and certainly can also be implemented by hardware. Based on this understanding, the above-mentioned technical solutions can be embodied in the form of software products in essence or the parts that make contributions to the prior art, and the computer software products can be stored in computer-readable storage media, such as ROM/RAM, magnetic A disc, an optical disc, etc., includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform the methods described in various embodiments or some parts of the embodiments.

最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, but not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that it can still be The technical solutions described in the foregoing embodiments are modified, or some technical features thereof are equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (10)

1. A positioning method based on millimeter wave equipment is characterized by comprising the following steps:
acquiring the number of targets to be detected and echo signals, and constructing a frequency spectrum signal model based on the number of the targets to be detected and the echo signals;
calculating a sampling error value of the target to be detected and a phase difference of adjacent antennas based on the frequency spectrum signal model;
calculating the distance value from the target to be detected to the millimeter wave equipment based on the sampling error value, and calculating the angle of arrival of the target to be detected based on the phase difference;
and generating a positioning result of the target to be detected according to the distance value and the wave arrival angle.
2. The positioning method based on the millimeter wave device according to claim 1, wherein the constructing a spectrum signal model based on the number of the targets to be measured and the echo signals specifically comprises:
extracting a target frequency band signal in the echo signal, and performing discrete sampling on the target frequency band signal to obtain a discrete target frequency band signal;
performing discrete Fourier transform on the discrete target frequency band signal to obtain a processing result, and acquiring a peak sampling point corresponding to the discrete target frequency band signal based on the processing result;
and selecting a plurality of sampling point data within a preset distance of the peak sampling point, and establishing a frequency spectrum signal model by combining the processing result.
3. The millimeter wave device based positioning method according to claim 2,
when the number of the targets to be detected is single, the discrete target frequency band signal is expressed as:
Figure FDA0003664731450000011
wherein N is a discrete time point of a discrete target frequency band signal, N is a sampling point number of the signal, and N is 0,1, …, N-1; e is the base number of the natural logarithm, and j is the unit of an imaginary number;
Figure FDA0003664731450000012
is the amplitude of the signal, A 0 And theta 0 Are all real numbers;
Figure FDA0003664731450000013
is the normalized frequency of the signal, in the range of [0,1 ]]Delta is the sampling error, delta is
Figure FDA0003664731450000014
Real numbers in the range, pi is the circumferential ratio, w [ n ]]Is gaussian white noise.
4. The millimeter wave device-based positioning method according to claim 3, wherein the discrete target frequency band signal is subjected to discrete Fourier transform to obtain a processing result, specifically:
performing discrete Fourier transform on the representation r [ n ] of the discrete target frequency band signal, wherein the processing result is as follows:
Figure FDA0003664731450000021
wherein m ═ k p -N+1,k p -N+2,…,k p ;W[k p -m]Is w [ n ]]The result of the discrete fourier transform of (a); k is a radical of p Is a peak sampling point;
selecting a plurality of sampling point data within a preset distance of the peak sampling point, and establishing a frequency spectrum signal model by combining the processing result, wherein the specific steps are as follows:
selecting three points R [ k ] p -1]、R[k p ]、R[k p +1]Establishing a spectrum signal model, and solving the model to obtain the complex amplitude of the discrete intermediate frequency signal
Figure FDA0003664731450000022
And normalized frequency
Figure FDA0003664731450000023
Calculating the distance of the target to be measured based on the normalized frequency:
Figure FDA0003664731450000024
F s is the sampling rate of the millimeter wave device.
5. The millimeter wave device based positioning method according to claim 2,
when the number of the targets to be detected is two, the discrete target frequency band signal is represented as:
Figure FDA0003664731450000025
wherein N is a discrete time point of a discrete target frequency band signal, N is a sampling point number of the signal, and N is 0,1, …, N-1;
Figure FDA0003664731450000026
and
Figure FDA0003664731450000027
is the amplitude of the complex signal that is unknown,
Figure FDA0003664731450000028
and
Figure FDA0003664731450000029
is the normalized frequency component of the signal, w n]Is white gaussian noise;
for the discrete target frequency band signal t [ n ]]Performing discrete Fourier transform to obtain operation result T [ k ]](ii) a Obtaining two peak sampling points k according to the operation result p1 、k p2
At the two peak sampling points k p1 、k p2 Selecting a plurality of sampling point data and combining the operation result T [ k ]]Establishing the spectrum signal model, and carrying out model pair on unknown parameters A 1 ,A 21212 Solving;
based on two normalized frequencies f 1 、f 2 And respectively calculating the distance of the target to be measured.
6. The millimeter wave device-based positioning method according to claim 5, wherein calculating the phase difference between adjacent antennas based on the spectrum signal model specifically comprises:
obtaining the complex signal amplitude in the discrete target frequency band signal model;
and extracting the phase difference of adjacent receiving antennas in the linear antenna array of the millimeter wave equipment according to the amplitude of the complex signal.
7. A positioning apparatus based on millimeter wave devices, the apparatus comprising:
the model building module is used for obtaining the number of the targets to be tested and echo signals and building a frequency spectrum signal model based on the number of the targets to be tested and the echo signals;
the first calculation module is used for calculating a sampling error value of the target to be detected and a phase difference of adjacent antennas based on the spectrum signal model;
the second calculation module is used for calculating a distance value from the target to be detected to the millimeter wave equipment based on the sampling error value and calculating a reaching angle of the target to be detected based on the phase difference;
and the result generation module is used for generating a positioning result of the target to be detected according to the distance value and the arrival angle.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements a positioning method based on millimeter wave devices according to any of claims 1 to 6 when executing the program.
9. A non-transitory computer-readable storage medium, on which a computer program is stored, wherein the computer program, when executed by a processor, implements the millimeter wave device-based positioning method according to any one of claims 1 to 6.
10. A computer program product comprising a computer program, wherein the computer program when executed by a processor implements a millimeter wave device based positioning method according to any of claims 1 to 6.
CN202210589994.XA 2022-05-26 2022-05-26 A positioning method and device based on millimeter wave equipment Active CN114966656B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210589994.XA CN114966656B (en) 2022-05-26 2022-05-26 A positioning method and device based on millimeter wave equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210589994.XA CN114966656B (en) 2022-05-26 2022-05-26 A positioning method and device based on millimeter wave equipment

Publications (2)

Publication Number Publication Date
CN114966656A true CN114966656A (en) 2022-08-30
CN114966656B CN114966656B (en) 2025-03-11

Family

ID=82956521

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210589994.XA Active CN114966656B (en) 2022-05-26 2022-05-26 A positioning method and device based on millimeter wave equipment

Country Status (1)

Country Link
CN (1) CN114966656B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024060790A1 (en) * 2022-09-21 2024-03-28 加特兰微电子科技(上海)有限公司 Method and apparatus for improving target detection precision, and electronic device
CN118152763A (en) * 2024-05-11 2024-06-07 北京智芯微电子科技有限公司 Distribution network data sampling method, device and electronic equipment

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080309546A1 (en) * 2007-06-13 2008-12-18 Mitsubishi Electric Corporation Radar device
CN103529444A (en) * 2013-09-27 2014-01-22 安徽师范大学 A vehicle-mounted millimeter-wave radar moving target recognizer and recognition method
CN109559525A (en) * 2018-11-26 2019-04-02 厦门精益远达智能科技有限公司 A kind of method for monitoring overspeed based on millimetre-wave radar, device and equipment
CN111308437A (en) * 2020-02-27 2020-06-19 南京慧尔视智能科技有限公司 Entropy-solving and speed-ambiguity-solving method for millimeter wave MIMO traffic radar
CN111736131A (en) * 2020-07-13 2020-10-02 深圳大学 A method for eliminating false targets of one-bit signal harmonics and related components
CN112526474A (en) * 2020-11-23 2021-03-19 哈尔滨工程大学 FMCW radar range-velocity joint estimation method based on full-phase Fourier transform
CN114236525A (en) * 2021-12-27 2022-03-25 浙江大学 Multi-target tracking and respiration detection method and device based on millimeter wave radar
CN114325680A (en) * 2021-12-27 2022-04-12 北京源清慧虹信息科技有限公司 Active power range measurement method in millimeter wave radar

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080309546A1 (en) * 2007-06-13 2008-12-18 Mitsubishi Electric Corporation Radar device
CN103529444A (en) * 2013-09-27 2014-01-22 安徽师范大学 A vehicle-mounted millimeter-wave radar moving target recognizer and recognition method
CN109559525A (en) * 2018-11-26 2019-04-02 厦门精益远达智能科技有限公司 A kind of method for monitoring overspeed based on millimetre-wave radar, device and equipment
CN111308437A (en) * 2020-02-27 2020-06-19 南京慧尔视智能科技有限公司 Entropy-solving and speed-ambiguity-solving method for millimeter wave MIMO traffic radar
CN111736131A (en) * 2020-07-13 2020-10-02 深圳大学 A method for eliminating false targets of one-bit signal harmonics and related components
CN112526474A (en) * 2020-11-23 2021-03-19 哈尔滨工程大学 FMCW radar range-velocity joint estimation method based on full-phase Fourier transform
CN114236525A (en) * 2021-12-27 2022-03-25 浙江大学 Multi-target tracking and respiration detection method and device based on millimeter wave radar
CN114325680A (en) * 2021-12-27 2022-04-12 北京源清慧虹信息科技有限公司 Active power range measurement method in millimeter wave radar

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
苏宏艳 等: "毫米波宽带雷达相位误差校正方法研究", 系统工程与电子技术, vol. 30, no. 02, 15 February 2008 (2008-02-15), pages 253 - 256 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024060790A1 (en) * 2022-09-21 2024-03-28 加特兰微电子科技(上海)有限公司 Method and apparatus for improving target detection precision, and electronic device
CN118152763A (en) * 2024-05-11 2024-06-07 北京智芯微电子科技有限公司 Distribution network data sampling method, device and electronic equipment

Also Published As

Publication number Publication date
CN114966656B (en) 2025-03-11

Similar Documents

Publication Publication Date Title
CN114966656A (en) Positioning method and device based on millimeter wave equipment
CN112630771B (en) Multi-target tracking method and system based on millimeter wave equipment
Long et al. High accuracy unambiguous angle estimation using multi‐scale combination in distributed coherent aperture radar
Dezfuli et al. Reduced complexity and near optimum detector for linear‐frequency‐modulated and phase‐modulated LPI radar signals
Zhang et al. Push the limit of millimeter-wave radar localization
Wei et al. Beam‐space reduced‐dimension space‐time adaptive processing for airborne radar in sample starved heterogeneous environments
RU2380719C2 (en) Method for location finding of radiation sources at one frequency
Fereidoony et al. Model‐based super‐resolution time‐delay estimation with sample rate consideration
Zhao et al. Multiple-target localization by millimeter-wave radars with trapezoid virtual antenna arrays
CN113109797B (en) Method and device for detecting target of frequency modulation continuous wave staring radar and computer equipment
CN113534077A (en) Radar radiation source power inversion method and device and electronic equipment
CN114518562A (en) A target identification method, device, electronic device and storage medium
Maaref et al. Array‐Based Ultrawideband through‐Wall Radar: Prediction and Assessment of Real Radar Abilities
US10371813B2 (en) Systems and methods for using time of flight measurements for imaging target objects
Vermesan et al. Ghost image cancellation algorithm through numeric beamforming for multi‐antenna radar imaging
CN116990777A (en) Terahertz RCS high-precision measurement method, system, device and equipment
Gu et al. Resolution threshold analysis of MUSIC algorithm in radar range imaging
Kim et al. Direction of arrival estimation of humans with a small sensor array using an artificial neural network
Meschino et al. A practical guide on using SPOT-GPR, a freeware tool implementing a SAP-DoA technique
Bai et al. Multi‐target pair‐matching method based on angle information in transmit/receive‐receive synergetic High Frequency Surface Wave Radar
Tian et al. SAPIL: single access point based indoor localisation using Wi‐Fi L‐shaped antenna array
Wen et al. Contactless Micron-Level Vibration Measurement with Millimeter Wave Radar
CN114814746A (en) Detection method, detection radar, equipment and storage medium for radar
CN114355324A (en) A track generation method
CN114280568A (en) Method and device for contour recognition of wall

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
GR01 Patent grant
GR01 Patent grant