CN113608221B - 一种加速度辅助机车雷达测速滤波方法 - Google Patents
一种加速度辅助机车雷达测速滤波方法 Download PDFInfo
- Publication number
- CN113608221B CN113608221B CN202110868586.3A CN202110868586A CN113608221B CN 113608221 B CN113608221 B CN 113608221B CN 202110868586 A CN202110868586 A CN 202110868586A CN 113608221 B CN113608221 B CN 113608221B
- Authority
- CN
- China
- Prior art keywords
- speed
- radar
- measurement
- acceleration
- speed measurement
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 75
- 230000001133 acceleration Effects 0.000 title claims abstract description 41
- 230000003137 locomotive effect Effects 0.000 title claims abstract description 30
- 238000001914 filtration Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 4
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect 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
- G01S13/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/91—Radar or analogous systems specially adapted for specific applications for traffic control
- G01S13/92—Radar or analogous systems specially adapted for specific applications for traffic control for velocity measurement
-
- 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
- G01S13/00—Systems 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/86—Combinations of radar systems with non-radar systems, e.g. sonar, direction finder
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种加速度辅助机车雷达测速滤波方法,其包括如下步骤:步骤S1、根据测速精度和测速实时性要求,选定滤波长度N;步骤S2、利用第N时刻以及之前的N‑1个时刻,共计N个时刻的雷达初始测量速度和加速度测量值计算得到N个速度估计值;步骤S3、对步骤S2得到的N个速度估计值进行均值滤波处理。本发明利用雷达初始测量速度和加速度计测量值得到多个估计速度,并对估计速度进行均值滤波得到最终的速度输出值,达到了减小雷达测速误差、提高测速精度的目的,本发明适用于减小非理想因素引起的雷达测速误差和提高雷达测速精度。
Description
技术领域
本发明属于机车测速系统技术领域,特别涉及一种加速度辅助机车测速雷达提高测速精度的方法。
背景技术
基于多普勒频移效应的雷达测速系统被广泛应用于地铁、火车等机车测速,其测量速度主要受测距精度、多普勒测量精度和雷达垂直高度变化的影响而存在误差。测距精度和多普勒测量精度主要由固有量化误差和系统温度漂移、噪声误差、频率源不稳定等引起,雷达垂直高度变化主要由机车颠簸引起。
为提高雷达测速精度,现有的大多数方法仅能减少单一误差因素的影响,例如双天线雷达测速仪,其通过安装两部载频不同的雷达,根据两部雷达固定的波束指向角度差,抵消因机车颠簸引起的测速误差。但两部雷达增加了测速系统成本,且该方法仅能减小机车颠簸引起的测速误差;或利用同步校频和互备双锁相环提高频率源稳定度,以达到减小测速误差的目的,该方法同样仅能减小因频率源不稳定而产生的误差;或通过对回波频谱加汉宁窗处理,以减弱频谱泄露,减小雷达测距误差,但窗函数参数影响测速精度。
发明内容
本发明借助加速度计辅助雷达测量机车速度,提供了一种加速度辅助机车雷达测速滤波方法,利用加速度信息和雷达初始测量速度实现测量速度的均值滤波,以减少雷达测距精度、多普勒测量精度和雷达垂直高度变化等的影响,从而达到减小测速误差的目的。
为实现上述目的,本发明提供一种加速度辅助机车雷达测速滤波方法,所述机车底部设置有雷达和加速度计,其包含如下步骤:
步骤S1、根据测速精度和测速实时性要求,选定滤波长度N;
步骤S2、利用第N时刻以及之前的N-1个时刻,共计N个时刻的雷达初始测量速度和加速度测量值计算得到N个速度估计值;
步骤S3、对步骤S2得到的N个速度估计值进行均值滤波处理。
优选地,所述步骤S2包括:记第N时刻以及之前的N-1个时刻,共N个时刻的雷达初始速度测量值依次为V1,V2…VN-1,VN,加速度计测量的加速度值依次为a1,a2…aN-1,aN。对加速度关于时间累加可以得到N个时刻的速度估计值
其中,△t表示相邻两个时刻的差值。
优选地,所述步骤S3中,均值滤波处理得到滤波后的速度V′N为:
综上所述,与现有技术相比,本发明提供的一种加速度辅助机车雷达测速滤波方法,具有如下有益效果:利用加速度计辅助机车测速雷达提高测速精度,通过对辅助加速度计的测量值进行积分并结合雷达初始测量速度,得到多个速度估计值,最后对速度估计值进行均值滤波,从而得到更为准确的测量速度。本发明能够改善测距精度、多普勒测量精度和雷达垂直高度变化等因素对雷达测速精度的影响,适用于减小雷达测速误差、提高雷达测速精度。
附图说明
图1为本发明的加速度辅助机车雷达测速滤波方法的实施流程图;
图2为本发明中的雷达测速示意图;
图3为本发明中的雷达测速几何图;
图4为本发明的加速度辅助机车雷达测速滤波方法的加速度辅助雷达测速原理图;
图5为本发明实施例为机车匀速运行,雷达初始测量速度和利用加速度信息滤波后速度的结果及其对比图;
图6为本发明实施例的雷达初始测量速度与滤波后测量速度误差分布图;
图7为本发明实施例在机车匀速运行时,雷达测速结果与激光测速结果对比图。
具体实施方式
以下将结合本发明实施例中的附图1~附图7,对本发明实施例中的技术方案、构造特征、所达成目的及功效予以详细说明。
需要说明的是,附图采用非常简化的形式且均使用非精准的比例,仅用以方便、明晰地辅助说明本发明实施方式的目的,并非用以限定本发明实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容能涵盖的范围内。
本发明提供了一种加速度辅助机车雷达测速滤波方法,在机车底部设置雷达和加速度计,用于测量机车的速度和加速度数据;通过引入加速度数据后,利用加速度能确保短时高精度的特点,结合雷达初始测量速度和加速度能够获取多个速度估计值,对速度估计值进行均值滤波处理得到最终的速度输出值,从而达到提高雷达测速精度的目的,如图1所示,该加速度辅助机车雷达测速滤波方法具体包括如下步骤:
步骤S1、根据测速精度和测速实时性要求,选定滤波长度N;
步骤S2、根据雷达初始测量速度和加速度值,通过计算得到N个速度估计值;
其中,N个速度估计值的具体计算步骤为:
记当前时刻(即第N个时刻)以及之前的N-1个时刻(共计N个时刻)依次为t1,t2,t3…tN-1,tN,两两相邻的时刻差值记为△t,且所有时间差值△t均相等;N个时刻对应的雷达初始测量速度依次为V1,V2,V3…VN,N个时刻对应的加速度计测量得到的加速度值依次为a1,a2,a3…aN。根据加速度值和雷达初始测量速度,通过对加速度关于时间累加可得到N个时刻的速度估计值
其中,上式中△t表示相邻两个时刻的差值。
步骤S3、对步骤S2得到的N个速度估计值进行均值滤波处理,经滤波处理后的速度V′N为:
其中,雷达初始测量速度的测速原理是基于多普勒效应实现的,如图2所示,雷达安装在待测速的机车底部,雷达发射探测波的方向与雷达垂直于地面向下的方向夹角为θ,则当机车以速度V相对于地面运动时,由机车与地面的相对运动而产生的多普勒频率为:
其中,λ为雷达发射的探测波波长,Vr为雷达相对于地面的径向速度。如图3所示,sinθ的计算公式如下:
其中,R为斜距,H为雷达高度。
雷达测速的流程为:雷达发射出探测波,在遇到地面之后反射回传的回波信号(回波),雷达对回波进行二维FFT(傅里叶变换)获取二维图像;提取最强散射点,获取其所在距离门,并计算斜距R;利用回波数据估计其所对应的多普勒频率fd;联合式(1)和式(2)计算机车速度V。
以下通过一个具体实施例详细说明本发明方法。本实施例中,加速度辅助高精度机车雷达测速滤波的参数为:机车运动速度V约为60km/h,雷达初始测量速度的时间间隔和加速度测量的时间间隔均为35.8ms,滤波长度N为100。
如图5所示,图5(a)为雷达初始测量速度图,图5(b)和图5(c)分别为利用加速度信息滤波后测量速度的结果图及对比图;通过对图5(c)的分析发现,由于受到误差因素的影响,雷达初始测量速度波动较大,经加速度信息滤波后速度波动减小。
如图6所示,图6为雷达初始测量速度误差与利用加速度信息滤波后速度误差分布图,雷达初始测量速度误差均值为-0.0128m/s、标准差为-0.1249m/s,滤波后速度测量误差均值为-0.0037m/s、标准差为-0.0531m/s。由此可见,利用加速度信息对雷达初始测量速度滤波,可以减小测距误差、提高测速精度。
如图7所示,以高精度激光测速结果为参考,图7(a)为加速度辅助雷达测量速度结果图,图7(b)为高精度激光测量速度结果图,图7(c)为加速度辅助雷达测量速度与高精度激光测量速度结果的对比图,从图中可见加速度辅助雷达测速结果与激光测速结果整体较吻合。因此,加速度辅助雷达测速能够提高雷达测速精度,且测速精度与激光测速精度相当。
综上所述,与现有提高雷达测速精度的方法相比,本发明所提供的加速度辅助机车雷达测速滤波方法能改善测距精度、多普勒测量精度和雷达垂直高度变化等因素对雷达测速精度的影响,适用于减小雷达测速误差、提高雷达测速精度。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (1)
1.一种加速度辅助机车雷达测速滤波方法,所述机车底部设置有雷达和加速度计,其特征在于,包括如下步骤:
步骤S1、根据测速精度和测速实时性要求,选定滤波长度N;
步骤S2、利用第N时刻及之前共N个时刻的雷达初始测量速度和加速度测量值计算得到N个速度估计值;
步骤S3、对步骤S2得到的N个速度估计值进行均值滤波处理;
所述步骤S2包括:记第N时刻以及之前的N-1个时刻,共计N个时刻的雷达初始速度测量值依次为V1,V2…VN-1,VN,加速度计测量的加速度值依次为a1,a2…aN-1,aN;对加速度关于时间累加可以得到N个时刻的速度估计值
其中,Δt表示相邻两个时刻的差值;
所述步骤S3中,均值滤波处理得到滤波后的速度VN′为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110868586.3A CN113608221B (zh) | 2021-07-29 | 2021-07-29 | 一种加速度辅助机车雷达测速滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110868586.3A CN113608221B (zh) | 2021-07-29 | 2021-07-29 | 一种加速度辅助机车雷达测速滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113608221A CN113608221A (zh) | 2021-11-05 |
CN113608221B true CN113608221B (zh) | 2024-04-26 |
Family
ID=78306092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110868586.3A Active CN113608221B (zh) | 2021-07-29 | 2021-07-29 | 一种加速度辅助机车雷达测速滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113608221B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116626328A (zh) * | 2023-07-26 | 2023-08-22 | 济南嘉宏科技有限责任公司 | 一种井下车辆全程测速设备及测速方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5043737A (en) * | 1990-06-05 | 1991-08-27 | Hughes Aircraft Company | Precision satellite tracking system |
DE102012222715A1 (de) * | 2012-12-11 | 2014-06-12 | Hyundai Motor Company | Vorrichtung und Verfahren zum vorausschauenden Erkennen einer Fahrzeugkollision |
CN104635231A (zh) * | 2015-02-12 | 2015-05-20 | 西安电子科技大学 | 基于自适应卡尔曼滤波的机车雷达测速方法 |
CN108415009A (zh) * | 2018-02-02 | 2018-08-17 | 合肥工大高科信息科技股份有限公司 | 一种列车速度计算方法及系统 |
CN108709521A (zh) * | 2018-06-04 | 2018-10-26 | 中核(天津)科技发展有限公司 | 一种高精度位移测量装置及测量方法 |
CN109633642A (zh) * | 2018-12-11 | 2019-04-16 | 上海无线电设备研究所 | 一种太赫兹高速目标雷达成像方法 |
CN112346045A (zh) * | 2020-10-29 | 2021-02-09 | 上海无线电设备研究所 | 一种多频率点雷达测速方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10641885B2 (en) * | 2017-07-21 | 2020-05-05 | Honeywell International Inc. | Systems and methods for measuring velocity and acceleration with a radar altimeter |
US11436927B2 (en) * | 2018-11-21 | 2022-09-06 | Brunswick Corporation | Proximity sensing system and method for a marine vessel with automated proximity sensor location estimation |
FR3090893B1 (fr) * | 2018-12-20 | 2021-10-08 | Commissariat Energie Atomique | Procede et systeme de mesure de la vitesse d’un porteur par rapport au sol |
US11754672B2 (en) * | 2019-05-13 | 2023-09-12 | Gm Cruise Holdings Llc | RADAR cross section compensation for calibration of vehicle RADAR |
US11699306B2 (en) * | 2019-08-16 | 2023-07-11 | Cambridge Mobile Telematics, Inc. | Method and system for vehicle speed estimation |
-
2021
- 2021-07-29 CN CN202110868586.3A patent/CN113608221B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5043737A (en) * | 1990-06-05 | 1991-08-27 | Hughes Aircraft Company | Precision satellite tracking system |
DE102012222715A1 (de) * | 2012-12-11 | 2014-06-12 | Hyundai Motor Company | Vorrichtung und Verfahren zum vorausschauenden Erkennen einer Fahrzeugkollision |
CN104635231A (zh) * | 2015-02-12 | 2015-05-20 | 西安电子科技大学 | 基于自适应卡尔曼滤波的机车雷达测速方法 |
CN108415009A (zh) * | 2018-02-02 | 2018-08-17 | 合肥工大高科信息科技股份有限公司 | 一种列车速度计算方法及系统 |
CN108709521A (zh) * | 2018-06-04 | 2018-10-26 | 中核(天津)科技发展有限公司 | 一种高精度位移测量装置及测量方法 |
CN109633642A (zh) * | 2018-12-11 | 2019-04-16 | 上海无线电设备研究所 | 一种太赫兹高速目标雷达成像方法 |
CN112346045A (zh) * | 2020-10-29 | 2021-02-09 | 上海无线电设备研究所 | 一种多频率点雷达测速方法 |
Non-Patent Citations (3)
Title |
---|
一种基于卡尔曼滤波的机车测速雷达算法;姬冰冰 等;火控雷达技术;第38卷(第1期);全文 * |
基于双雷达的轨道交通测速方法及误差分析;周媛 等;控制与信息技术(第1期);全文 * |
干涉SAR运动舰船检测及速度估计方法研究;欧阳克威;中国优秀硕士学位论文全文数据库;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113608221A (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107678003B (zh) | 一种地波雷达海杂波背景下的目标检测方法 | |
CN108037498B (zh) | 基于三角波线性调频连续波雷达的高速目标测速测距方法 | |
CN109375202B (zh) | 一种基于车载毫米波雷达的车辆测距测速方法 | |
CN111077524B (zh) | 一种sar-gmti动目标重定位改进方法 | |
CN112051568B (zh) | 一种两坐标雷达的俯仰测角方法 | |
CN104635231B (zh) | 基于自适应卡尔曼滤波的机车雷达测速方法 | |
CN104280728A (zh) | 雷达传感器与用于运行雷达传感器的方法 | |
CN113608221B (zh) | 一种加速度辅助机车雷达测速滤波方法 | |
CN109901162B (zh) | 一种适用于分布式地球同步轨道sar的长基线成像stap方法 | |
JP2014044109A (ja) | 合成開口レーダ装置 | |
CN111505618A (zh) | 一种适用于车载毫米波雷达领域的基于频率估计的去耦合校正测距方法 | |
JP3518363B2 (ja) | Fmcwレーダ装置及び記録媒体並びに車両制御装置 | |
CN110208788B (zh) | 一种基于滑窗的连续帧联合速度解模糊方法 | |
JP3565862B2 (ja) | 車道での電磁波の反射を利用するレーダを用いた車両の対地速度の測定方法 | |
JP5278128B2 (ja) | 変換係数導出装置および変換係数導出方法 | |
CN116719023A (zh) | Gpr能量聚焦和分辨率联合约束下的介质介电常数估计方法 | |
CN110967677A (zh) | 时频域分级多目标分辨与测量方法 | |
CN116699572A (zh) | 基于二次补偿的相干激光雷达振动误差补偿方法 | |
CN115616629A (zh) | 一种基于天基外辐射源信号的运动目标探测补偿方法 | |
CN112068122B (zh) | 基于跟踪信息的相位编码信号的多普勒补偿方法 | |
KR20150100051A (ko) | 레이더 고도계 장치 및 방법 | |
CN111198358B (zh) | 一种基于相关检测的多旋翼无人机回波积累方法 | |
CN110736988B (zh) | 双基地pfa运动目标参数估计和成像方法 | |
CN113253279A (zh) | 一种正反线性调频信号测距测速方法 | |
CN106324590A (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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20230506 Address after: 1555 Zhongchun Road, Minhang District, Shanghai, 201109 Applicant after: SHANGHAI RADIO EQUIPMENT Research Institute Applicant after: SHANGHAI SHENTIAN INDUSTRIAL CO.,LTD. Address before: No. 1555, Zhongchun Road, Minhang District, Shanghai, 200233 Applicant before: SHANGHAI RADIO EQUIPMENT Research Institute |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |