CN112379335B - 一种基于幅度差分信息的到达时间测量方法 - Google Patents
一种基于幅度差分信息的到达时间测量方法 Download PDFInfo
- Publication number
- CN112379335B CN112379335B CN202011266345.3A CN202011266345A CN112379335B CN 112379335 B CN112379335 B CN 112379335B CN 202011266345 A CN202011266345 A CN 202011266345A CN 112379335 B CN112379335 B CN 112379335B
- Authority
- CN
- China
- Prior art keywords
- pulse
- arrival time
- envelope
- time
- differential
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000001514 detection method Methods 0.000 claims description 8
- 230000000630 rising effect Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 abstract description 33
- 238000012545 processing Methods 0.000 abstract description 3
- 238000000691 measurement method Methods 0.000 description 7
- 238000009499 grossing Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
- G01S7/2927—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods by deriving and controlling a threshold value
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于幅度差分信息的到达时间测量方法，其特征在于，包括以下步骤：1)对接收机接收到的相邻时刻的脉冲信号包络值进行差分；2)对步骤1)中得到的包络值差分值进行比较大小，并基于得到的差分最大值所在位置处确定脉冲到达时间。本发明方法对脉冲到达时间测量时，通过对脉冲包络进行差分，使得到达时间的测量不受幅度的影响，进而在较高的动态范围情形下，幅度差分的到达时间测量具有更高的测量精度。因此，本发明可以广泛应用于信号处理技术领域。
Description
技术领域
本发明涉及信号处理技术领域，特别是关于一种基于幅度差分信息的到达时间测量方法。
背景技术
随着无线电电子技术和电子对抗技术的迅速发展，各种辐射源的功能不断加强，电磁环境变得越来越复杂，如何高效获取雷达信息，并对多个同时到达、复杂未知调 制的雷达信号进行分选和信号识别处理，已经成为雷达探测领域面临的一个重大难题。 其中，雷达信号的脉冲到达时间是重要的雷达脉冲描述参数，提高脉冲到达时间的测量精度具有重要的意义。
目前，常用的脉冲到达时间测量方法主要基于固定门限，然而在感知雷达信号时，由于雷达波束的扫描特性以及距离的变化等原因，接收机依次感知到雷达主瓣不同位 置的脉冲信号，即由雷达主瓣发射功率最低，接收的脉冲幅度最小处动态变化到发射 功率最高，接收的脉冲幅度最大处，接收雷达脉冲功率动态范围很大，严重影响雷达信号到达时间测量的精度。与此同时，脉冲到达时间的测量精度还会受到随机噪声的 干扰，因此需要一种硬件上实现简单的算法消除上述问题带来的影响。
发明内容
针对上述问题，本发明的目的是提供一种基于幅度差分信息的到达时间测量方法， 通过对脉冲包络进行差分，使得对脉冲到达时间的测量更加精确。
为实现上述目的，本发明采取以下技术方案：一种基于幅度差分信息的到达时间测量方法，其包括以下步骤：
1)对接收机接收到的相邻时刻的脉冲信号包络值进行差分；
2)对步骤1)中得到的包络值差分值进行比较大小，并基于得到的差分最大值所在位置处确定脉冲到达时间。
进一步地，所述步骤1)中，对接收机接收到的相邻时刻的脉冲信号包络值进行 差分时，采用常规差分方法或多点差分方法。
进一步地，所述步骤1)中，采用常规差分方法对接收机接收到的脉冲信号进行 平滑处理时，计算公式为：
式中，k(n)为脉冲到达时间，p(n)和p(n-1)分别为n时刻和n-1时刻的脉冲包络值，n为脉冲包络检测时间序号。
进一步地，所述步骤1)中，采用多点差分方法对接收机接收到的脉冲信号进行 平滑处理时，计算公式为：
式中，kM(n)为M点差分信号；M为差分点数；p(n)为脉冲包络，n为脉冲包络检 测时间序号。
进一步地，采用多点差分值方法对脉冲信号进行处理时，差分点数M取24。
进一步地，所述步骤2)中，脉冲到达时间的计算公式为：
ta＝tr/2
式中，ta为脉冲到达时间，tr为脉冲上升沿时间。
本发明由于采取以上技术方案，其具有以下优点：1、本发明对脉冲到达时间测量时，通过对脉冲包络进行差分，使得到达时间的测量不受幅度的影响，在较高的动态 范围情形下，幅度差分的到达时间测量具有更高的测量精度。2、本发明在对脉冲包络 进行差分之前，首先采用多点差分值平滑的方法，对接收到的脉冲信号进行平滑处理， 解决了差分带来的随机误差增大的问题，使得对到达时间的测量更加精确。本发明可 以广泛应用于信号处理技术领域。
附图说明
图1是现有的固定门限脉冲到达时间测量示意图；
图2是本发明实施方式中脉冲到达时间测量误差和差分点数的关系曲线；
图3是本发明实施方式中动态范围两种到达时间测量方法的测量误差曲线；
图4是本发明实施方式中动态范围、随机噪声同时存在两种到达时间测量方法的测量误差曲线。
具体实施方式
以下结合附图来对本发明进行详细的描绘。然而应当理解，附图的提供仅为了更好地理解本发明，它们不应该理解成对本发明的限制。
首先对常规的固定门限到达时间测量方法进行介绍。常规测量方法中，通常在低仰角接收条件下，采用固定的检测门限，因而脉冲到达时间的测量值tn会引入两种误 差：脉冲幅度变化带来的误差Δt1和随机误差Δt2。随机误差是指环境中噪声等随机因 素给到达时间测量带来的误差。
设脉冲幅度达到固定检测门限的时刻为ta，则脉冲幅度的变化给到达时间的测量带来的误差如图1所示。图1中，tr为脉冲上升沿时间，tf为脉冲下降时间，td表示 幅度保持时间。p1(t)和p2(t)具有相同的包络形状和不同的幅度A1、A2，因此采用固 定门限ATH，将得到两个不同的到达时间测量值ta1和ta2。令脉冲幅度为A，可将脉冲 包络p(t)表示为：
由式(1)可知，当0≤t＜tr时，p(t)的切线斜率为：
考虑到脉冲幅度在实际接收过程中动态范围很大，因此固定门限到达时间测量方法中，脉冲幅度带来的误差Δt1对辐射源定位精度的影响十分显著。
脉冲到达门限时刻ta的表达式为：
由以上公式可知，脉冲到达门限时刻ta随脉冲包络幅度A变化而变化，即脉冲幅度变化会带来到达时间的测量误差Δt1。
实施例一
为了减小这方面的影响，本发明提供一种基于幅度差分信息的到达时间测量方法， 基于幅度差分的方法来确定脉冲到达时间，具体的，包括以下步骤：
1)采用幅度差分方法对接收机接收到的相邻时刻脉冲信号的包络值进行差分。
考虑到脉冲幅度在实际接收过程中动态范围很大，因此对到达时间测量精度的影响十分显著。为了减小这方面的影响，可以考虑采用差分的方法确定脉冲到达时间。 对相邻时刻脉冲的包络值进行差分并比较大小，将最大值所在位置确定为脉冲的到达 时间。
式中，k(n)为脉冲到达时间，p(n)和p(n-1)分别为n时刻和n-1时刻的脉冲包络值，n为脉冲包络检测时间序号。
2)对得到的各差分值进行比较，并基于得到的差分最大值所在位置得到脉冲到达时间，其中，本发明中将脉冲到达时间定义为脉冲上升沿斜率最大处。
对公式(4)进行求导，得到其切线斜率方程，也就是脉冲到达时间为：
ta＝tr/2 (5)
可以看出，差分方法下的脉冲到达时间出现在上升沿的斜率最大处，根据公式(5)可知，差分方法下的脉冲到达时间与脉冲幅度无关，可以大幅减小脉冲幅度变化对到 达时间测量带来的影响。
实施例二
由于差分方法会放大随机噪声对到达时间测量的影响，因此，本实施例提供的一种基于幅度差分信息的脉冲到达时间测量方法，其包括以下步骤：
1)采用多点差分值平滑的方法，对接收机接收到的相邻时刻脉冲信号的包络值进行差分。
采用多点差分值平滑的方法对脉冲信号进行平滑处理时，表达式如下所示：
式中，kM(n)为M点差分信号；M为差分点数；p(n)为脉冲包络，n为脉冲包络检 测时间序号。
2)对得到的各多点差分值进行比较，并基于得到的多点差分最大值所在位置得到脉冲到达时间。
上述步骤1)中，类似的，如果M的值取得很小时，噪声平滑效果不明显；如果M 的值取得很大，脉冲上升沿的信息将被平滑处理，反而导致测量误差增大。因此平滑 长度的上限应当由实际应用中最窄脉冲的上升沿宽度确定。需要结合具体问题进行分 析，确定最优的差分点数。
实施例三
本实施例针对128MHz采样率下脉宽为1μs的窄脉冲进行仿真，得到到达时间测 量误差和差分点数的关系曲线如图3所示。可以看出随着差分点数的增加，到达时间 测量误差先减小后增大，与前面的理论分析相符。为了更准确的获得最优差分点数， 针对16～32点差分进行了步长更小的仿真。可以看出当差分点数取24时，到达时间 测量误差最小，因此本发明采用24点差分实现到达时间测量。
实施例四
本实施例中，设定脉冲宽度为6μs，通过数字信道化后采样频率为128MHz，脉冲 最大动态范围40dB。在没有随机噪声影响的条件下，分别采用固定门限方法和本发明 基于幅度差分信息的方法对不同功率的脉冲到达时间进行测量，结果如图4所示。可 以看出，本发明方法在高动态范围下测量精度几乎不发生变化，而固定门限的测量误 差显著增大，说明基于固定门限的到达时间测量方法受脉冲幅度的影响远大于本发明 提出的基于幅度差分信息的到达时间测量方法。
在脉冲最大动态范围40dB的基础上添加随机噪声，即最大信噪比变化范围为15dB～55dB的条件下，分别采用固定门限方法和包络差分方法对不同功率的脉冲进行 到达时间测量，结果如图4所示。可以看出，在动态范围较低的条件下，固定门限的 到达时间测量具有更高的精度。而随着动态范围增大，包络差分的到达时间测量精度提高，而固定门限精度下降，因此本发明基于包络差分的方法具有更高的测量精度。
上述各实施例仅用于说明本发明，其中各部件的结构、连接方式和制作工艺等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应 排除在本发明的保护范围之外。
Claims (3)
1.一种基于幅度差分信息的到达时间测量方法，其特征在于，包括以下步骤：
1）对接收机接收到的相邻时刻的脉冲信号包络值进行差分；
其中，对接收机接收到的相邻时刻的脉冲信号包络值进行差分时，采用常规差分方法或多点差分方法；
采用常规差分方法对接收机接收到的脉冲信号进行平滑处理时，计算公式为：
式中，为脉冲到达时间，和分别为时刻和时刻的脉冲包络值，n为脉冲包络检测时间序号；
采用多点差分方法对接收机接收到的脉冲信号进行平滑处理时，计算公式为：
式中，为M点差分信号；M为差分点数；为脉冲包络，为脉冲包络检测时间序号；
2）对步骤1）中得到的包络值差分值进行比较大小，并基于得到的差分最大值所在位置处确定脉冲到达时间。
2.如权利要求1所述的一种基于幅度差分信息的到达时间测量方法，其特征在于，采用多点差分值方法对脉冲信号进行处理时，差分点数M取24。
3.如权利要求1所述的一种基于幅度差分信息的到达时间测量方法，其特征在于，所述步骤2）中，脉冲到达时间的计算公式为：
式中，为脉冲到达时间，为脉冲上升沿时间。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011266345.3A CN112379335B (zh) | 2020-11-13 | 2020-11-13 | 一种基于幅度差分信息的到达时间测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011266345.3A CN112379335B (zh) | 2020-11-13 | 2020-11-13 | 一种基于幅度差分信息的到达时间测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112379335A CN112379335A (zh) | 2021-02-19 |
| CN112379335B true CN112379335B (zh) | 2024-07-26 |
Family
ID=74583702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011266345.3A Active CN112379335B (zh) | 2020-11-13 | 2020-11-13 | 一种基于幅度差分信息的到达时间测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112379335B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101790888B1 (ko) * | 2016-08-09 | 2017-10-26 | 서울대학교산학협력단 | 초음파를 이용한 거리 측정 방법 및 이를 이용한 거리 측정 장치 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102495392B (zh) * | 2011-12-05 | 2013-06-12 | 中国电子科技集团公司第三十八研究所 | 高精度测量脉冲到达时间的方法和装置 |
| CN108509377B (zh) * | 2018-03-19 | 2021-07-27 | 东南大学 | 一种基于沿特征提取的脉冲信号到达时间与脉宽估计方法 |
| CN109347549B (zh) * | 2018-08-17 | 2020-11-06 | 中国科学院上海应用物理研究所 | 一种高带宽射频信号到达时间测量装置和方法 |
| CN111901743A (zh) * | 2020-04-01 | 2020-11-06 | 中兴通讯股份有限公司 | 到达时间的确定方法、装置、信号接收设备及存储介质 |
-
2020
- 2020-11-13 CN CN202011266345.3A patent/CN112379335B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101790888B1 (ko) * | 2016-08-09 | 2017-10-26 | 서울대학교산학협력단 | 초음파를 이용한 거리 측정 방법 및 이를 이용한 거리 측정 장치 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112379335A (zh) | 2021-02-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7683825B2 (en) | Method and apparatus for measuring distance between a target and a receiver in a ranging system | |
| CN110515037B (zh) | 非视距环境下能时频多域联合的被动定位方法 | |
| CN112965037B (zh) | 雷达信号波形不确定度测试系统 | |
| CN105223481A (zh) | 基于差值能量函数的局部放电特高频信号起始时刻确定方法 | |
| CN107861112A (zh) | 一种全数字调节的激光测距系统 | |
| CN112462356B (zh) | 一种基于高斯插值提高雷达测距精度的方法 | |
| CN105891795B (zh) | 一种双站rcs测量定标方法 | |
| CN115236465A (zh) | 一种局部放电定位特高频信号时间差提取方法与系统 | |
| CN112379335B (zh) | 一种基于幅度差分信息的到达时间测量方法 | |
| CN112448774B (zh) | 一种基于外辐射信号的宽带射频接收处理系统自检方法 | |
| CN207586422U (zh) | 一种全数字调节的激光测距系统 | |
| CN108718223B (zh) | 一种非合作信号的盲频谱感知方法 | |
| CN114035149B (zh) | 一种提升干涉仪测向系统灵敏度的方法 | |
| CN111693982B (zh) | 一种雷达测速ai算法 | |
| CN212905422U (zh) | 宽带数字测频模块 | |
| CN112763023A (zh) | 基于优化数据模型的雷达物位计高精度测量输出处理方法 | |
| CN103336252B (zh) | 迟滞时间差型磁通门传感器信号检测方法 | |
| Mbouombouo Mboungam et al. | Clutter Map Constant False Alarm Rate Mixed with the Gabor Transform for Target Detection via Monte Carlo Simulation | |
| CN105738878A (zh) | 一种基于频偏分类的精密测量雷达测频方法 | |
| CN105737728A (zh) | 一种金属层脉冲涡流测厚方法 | |
| CN111929671A (zh) | 基于正态匹配的无线电侦测的接收幅度的估计方法及系统 | |
| Wang et al. | Research on a New Model Structure for Unknown Radar Emitter Signal Sorting | |
| CN111404518A (zh) | 一种灵活可调的延时比较方法及系统 | |
| CN113777412B (zh) | 一种提高天线方向图零深位置测试精度的方法 | |
| CN111142076B (zh) | 提高雷达低截获性能的功率控制方法 |
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 |