CN102236092A - 激光多普勒测速仪的多普勒信号周期分布解调方法 - Google Patents
激光多普勒测速仪的多普勒信号周期分布解调方法 Download PDFInfo
- Publication number
- CN102236092A CN102236092A CN2010101601892A CN201010160189A CN102236092A CN 102236092 A CN102236092 A CN 102236092A CN 2010101601892 A CN2010101601892 A CN 2010101601892A CN 201010160189 A CN201010160189 A CN 201010160189A CN 102236092 A CN102236092 A CN 102236092A
- Authority
- CN
- China
- Prior art keywords
- length
- signal
- period
- doppler
- cycle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种激光多普勒测速仪的多普勒信号周期分布解调方法。该方法采用基于周期计数法的时域处理原理,但是它不计算每个给定时间段内全部周期的平均频率,而是分析各时间段内每个周期的时间值分布,取它们中最频繁的时间值作为该时间段内的周期长度,换算为该时间段的信号频率。本发明的主要优点是,可以提高多普勒频移信号处理的计算精度,有效减小散射光不稳定带来的多普勒频移信号周期变宽或变窄、信号局部丢失和噪声对目标运动速度测量结果的影响。
Description
技术领域
本发明属于激光多普勒测速技术领域,主要涉及一种激光多普勒测速仪的多普勒信号时域解调方法,尤其涉及一种对固定运动目标速度进行测量的外差式干涉计的多普勒信号周期分布解调方法。
背景技术
光源与物体相对运动时会产生多普勒效应,激光多普勒测速仪是依靠运动微粒散射光与照射光之间的频移获得速度信息的,常用的参考光模式的光学外差干涉计原理见附图1。单一频率f0的激光束照射到一个运动物体上时,接收到的散射光频率为f0+ΔfD,其中ΔfD为散射光与照射光之间的频移。另一束参考光通过分光器取自同一发射光束,经布喇格声光单元(Bragg-cell)引入40MHz的附加频移,这一参考光束与接收的散射光束相调制,光检测器(例如光电倍增管)接收调制光的光强变化,外差干涉计产生的多普勒电信号频率为:
fdet=40MHz+ΔfD (1)
对该多普勒信号进行频率解调处理,即可以获得散射光频率ΔfD,从而得到运动物体的速度V:
式中,λ为发射的激光光源的光波波长,β为投射光束与物体运动方向的夹角。
多普勒信号处理可以采用多种方法,时域处理的周期计数法直接计算每个给定时间段tg内的信号周期数N,则
其原理简单,处理简便,是目前激光多普勒测速仪常用的硬件信号处理方法,但是,这种采用周期计数方式的多普勒时域解调方法的抗干扰能力较差。
在实际的激光多普勒测速应用中,运动物体表面微粒分布不均匀或微粒重叠、激光束的会聚点因物体振动而偏离散射面、激光发射和接收器本身的振动、光束路径空间的微尘等外部因素都会影响接收到的散射光波的稳定性,使采集的多普勒频移信号中形成周期变宽或变窄、信号局部丢失和噪声,从而增大频移量处理的动态误差,影响物体运动速度的测量精度。实验证实:采用侧向投射和接收光束测速方法时散射光波的稳定性更差。某气炮激励高冲击校准系统采用双路侧向激光多普勒测速方法测量固体目标的运动速度,示意图见附图2,图中P为运动目标,S1和S2为位于运动目标后侧向的两个激光束投射及散射光接收点。该系统采用两路互相独立的参考光模式外差多普勒测速仪获取目标速度V0在侧后向的速度分量V1和V2,两路测量光束与目标运动方向的夹角分别为β1和β2。附图3为通过外差多普勒测速仪硬件对多普勒信号解调处理后获取的两条速度曲线,图中为了便于观察,将速度分量V2倒置显示,可见所获取的速度曲线存在许多突变的毛刺。观察外差干涉计输出的多普勒信号,其中增速时段内的1μs时间间隔的信号波形在图4中显示,上部为40MHz附加频移信号曲线,中部为S1方向的多普勒信号曲线,下部为S2方向的多普勒信号曲线,可见这一时间段内S2方向的多普勒信号存在周期长度突变。由于采用多普勒信号的周期计数解调方法求出的是一定时间段tg内全部周期N的平均频率,这种散射光波不稳定引起的信号周期变化会影响到频率解调结果,这就是周期计数解调方法抗干扰能力较差并会形成速度测量曲线毛刺的本质原因。
发明内容
本发明要解决的技术问题是,针对现有多普勒信号周期计数解调方法存在的问题,为激光多普勒测速仪提供一种抗干扰能力强的多普勒信号周期分布解调方法,该方法能够在脱机方式下对外差式干涉计输出的多普勒信号进行频率解调,提高对固体目标运动速度的测量精度。
为解决上述技术问题,本发明周期分布解调方法包括以下步骤:
(1)设置多普勒信号解调处理的周期计数时间段长度tg,以及在一个tg内可能出现的最大周期数Nmax和最少周期数Nmin,tg根据测速目标的最大加速度而定,一个时间段长度tg内可能的最多周期数Nmax和最少周期数Nmin通过估算目标运动的最小速度和最大速度确定;
(2)调用存盘的测速多普勒信号,并对该信号进行限幅处理;
(3)求出测速多普勒信号总时间长度ta内所含时间段长度tg的数量M;
(4)设k为调用频率解调处理子程序的循环次数变量,即1≤k≤M,对k从1递增到M,逐次调用频率解调处理子程序进行如下计算和判断:
(4.1)计算当前时间段长度tg内每个信号周期的长度tj,并计数当前时间段长度tg内包括的总周期数N,j=1~N;
(4.2)比较该时间段内的总周期数N与可能的最多周期数Nmax和最少周期数Nmin,在Nmin≤N≤Nmax的情况下,统计当前时间段长度tg内全部信号周期长度tj的分布,并从中求出最频繁的信号周期长度tm,该周期长度对应的周期数为Nm,1≤m≤N;
(4.3)若Nm>1,以周期长度接近tm的所有信号周期的统计平均值为当前时间段长度tg的周期长度td,即设定允差δt,统计周期长度在tm-δt至tm+δt间的所有周期的平均值;若Nm=1,则将N个信号周期长度tj的平均值作为当前时间段长度tg的周期长度td;
(4.4)在N<Nmin或N>Nmax的情况下,若k=1,则令周期长度td为测速目标在静止状态下多普勒信号的周期长度;若k>1,则将第k-1次频率解调计算子程序运行中解算出的周期长度td作为当前时间段长度tg的周期长度td;
(5)将当前时间段tg内的周期长度td转换为当前频率值;
(6)将信号的当前频率值转换为测速目标的当前速度值;
(7)重复第(4)~(6)步,直到k=M。
根据本发明,所述的多普勒信号周期的计数时间段长度tg为0.5μs。
根据本发明,δt的设定值为tm的5%。
本发明的多普勒信号周期分布解调方法采用基于周期计数法的时域处理原理,但是它不计算每个给定时间段tg内全部周期N的平均频率,而是分析各时间段tg内每个周期的时间值分布,取它们中最频繁的时间值td作为该时间段tg内的周期长度,换算为该时间段的信号频率。采用本发明提出的多普勒信号频率解调,可以有效消除信号中个别周期长度异常引起的多普勒频移计算误差,抑制因散射光波不稳定形成的多普勒频移信号周期变宽或变窄、信号局部丢失和噪声,改善周期计数法频率解调的信号处理精度。
附图说明
图1是参考光模式的外差干涉计原理图。
图2是双路侧向多普勒测量固体目标的运动速度示意图。
图3是经测速仪硬件对多普勒信号解调后获取的速度曲线。
图4是外差干涉计输出的实际多普勒信号图。
图5是本发明周期分布解调方法的主程序的工作流程图。
图6是图5中所示频率解调计算子程序的工作流程图。
图7是应用本发明对多普勒信号解调后获取的速度曲线。
具体实施方式
下面结合附图及优选实施例对本发明作进一步的详述。
本发明的优选实施例是通过计算机数据处理系统实现的。当用外差干涉计对测速目标进行激光多普勒测速时,外差干涉计输出的测速多普勒信号被数字示波器实时采集和存储。测量结束后,计算机数据处理系统读取数字示波器保存的测速多普勒信号并转存在计算机硬盘中。随后,可按照图5所示的工作流程进行多普勒信号的频率解调。
第一步,设置解调处理所用到的相关参数。这些参数包括多普勒信号解调处理的周期计数时间段长度tg,以及在一个tg内可能出现的最多周期数Nmax和最少周期数Nmin。如背景技术部分所述,外差干涉计输出的多普勒信号是频率为几十MHz的高频信号,而测速目标的速度变化是机械运动的低频信号。根据应用系统的技术指标估算测速目标运动的最大加速度,根据式(2)、式(3)确定时间段长度tg,使该时间段内的多普勒信号频率改变量在测量精度许可范围内。对于本优选实施例来说,测速目标运动的最大加速度为2,000,000m/s2,因此,设定时间段长度tg为0.5μs。一个时间段长度tg内可能的最多周期数Nmax和最少周期数Nmin通过估算目标运动的最小速度和最大速度确定。在本实施例中,测量目标从静止状态增速,速度上限为130m/s,测量光束与目标运动方向间的夹角β为78°,则在测量光束方向上的速度分量最大值为30m/s,最小值为0m/s。由于散射光束接收点S1和S2位于运动目标侧后向,即目标背离观察点运动,根据多普勒效应原理和式(2)、式(3)可知:当目标从原始的静止状态加速到可能的最大速度,频移信号的频率从40MHz减小到18.6MHz,即在时间段长度tg(0.5μs)内的信号周期数从20减小到9。考虑可能的允差,适当加宽估算的信号频率范围,Nmin取值为小于(0.8×18.6×106×tg)的第一个整数7、Nmax取值为大于(1.2×40×106×tg)的第一个整数24。
第二步,调用存盘的测速多普勒信号,并对该测速多普勒信号进行限幅处理。对一个变频信号的频率解调,只需要对信号每个周期的时间长度进行处理,信号的幅值不影响处理结果,通过限幅处理处理还可以消除信号顶部和底部的高频干扰。对于实施例中的外差干涉计输出的测速多普勒信号,直接取信号半波进行处理,即将信号上限限幅到零电平。
第三步,求出测速多普勒信号总时间长度ta内所含时间段长度tg的数量M,M为小于ta/tg的第一个整数。
第四步,设k为调用频率解调处理子程序的循环次数变量,即1≤k≤M,对k从1递增到M,逐次调用频率解调处理子程序,按照图6所示的工作流程进行如下计算和判断:
(4.1)计算当前时间段长度tg内每个信号周期的长度tj,并对当前时间段长度tg内包含的信号周期计数,j=1~N,这里j对应于时间段tg内第j个信号周期,N为该时间段长度tg内包含的总周期数。信号周期长度的计算根据被处理信号限幅后的上限进行,本实施例是根据零线上每个周期的持续时间计算周期长度。
(4.2)比较该时间段内的总周期数N与可能的最多周期数Nmax和最少周期数Nmin。在Nmin≤N≤Nmax的情况下,统计当前时间段长度tg内全部信号周期长度tj的分布,并从中求出最频繁的信号周期长度tm,该周期长度对应的周期数为Nm,1≤m≤N。所谓信号周期长度分布是指:设当前时间段长度tg内共有L个不同的周期长度ti(i=1~L),统计出对应每个周期长度为ti的周期个数Ni(i=1~L)。所谓最频繁的信号周期长度是指:求出信号周期个数集合Ni(i=1~L)中的最大值Nm,该最大值对应的周期长度即为tm。
(4.3)如果找到tm,即Nm>1时,逐个比较全部信号周期的分布,以周期长度接近tm的所有信号周期的统计平均值为当前时间段长度tg的周期长度td,即:
δt为所述的周期长度接近tm的允差,在本实施例中,δt的值设定为tm的5%。
如果找不到tm,即Nm=1时,说明本时间段tg内每一个信号周期的长度都不相同,则将N个信号周期长度tj的平均值作为当前时间段长度tg的周期长度td,这也就是传统的周期计数法的处理方法:
(4.4)在N<Nmin或N>Nmax的情况下,如果是第一次调用频率解调子程序,即k=1,则令周期长度td为测速目标在静止状态下多普勒信号的周期长度,即取td为1/40μs。若k>1,则将第k-1次频率解调计算子程序运行中解算出的周期长度td作为当前时间段长度tg的周期长度td。
第五步,按照下式将当前时间段tg内的周期长度td转换为当前频率值:
第六步,根据式(1)和式(2)将信号的当前频率值转换为测速目标的当前速度值,并将当前速度值送入计算机的显示器进行显示;
第七步,重复第四~六步,直到k=M,此时,计算机显示屏所显示的是一条速度随时间变化的曲线,保存计算结果。
本发明基于给定时间段长度tg内的周期分布统计,简称为“周期分布解调法”。这一解调方法是从分析应用系统所获取的多普勒信号特性得出的,适用于固体目标从静止状态开始运动的速度测量。
下面结合本实施例对本发明产生的技术效果进行分析。
对于附图4所示的在给定时间段内多普勒信号的大多数周期的长度基本不变,也存在个别周期的长度变宽或变窄的情况,本发明的频率解调处理子程序分析给定时间段tg内每个信号周期的时间长度分布,再取它们中最频繁的时间值td作为该时间段tg内的周期长度,可以有效消除因个别周期长度异常引起的多普勒频率计算误差。在更一般的情况下,散射光波不稳定性也会引起信号局部丢失、噪声等现象,按照本方法所述的解调处理步骤,通过信号限幅处理可以消除信号顶部或底部噪声(在本实施例中不必考虑信号底部噪声)对解调结果的影响;步骤(4.4)对应当前时间段长度tg的信号存在十分显著的周期变宽或变窄、信号局部丢失或噪声等现象,在这些情况下,用前一个信号正常的时间段求出的周期长度作为本时间段的处理结果,由于测速目标的机械运动速度变化频率远低于多普勒信号的频率,这种方法可以排除个别时间段的信号严重失真对多普勒频移处理的影响。因此,本发明提出的多普勒信号周期分布解调方法对目标速度的处理精度比传统的周期计数法高。
将本发明的多普勒信号时域解调方法应用于背景技术部分所述的气炮激励高冲击校准系统,解调处理所获得的两条速度分量曲线V1和V2图示于附图7,图中为了便于观察,将速度分量V2曲线倒置显示。与附图3中同一应用系统通过测速仪硬件解调后获取的速度曲线相比较,可以看出速度曲线上突变的毛刺得到有效抑制。对所获得的速度曲线进行低通滤波,即可获得平滑的目标速度的变化曲线。
Claims (3)
1.一种激光多普勒测速仪的多普勒信号周期分布解调方法,其特征在于,该方法包括以下步骤:
(1)设置多普勒信号解调处理的周期计数时间段长度tg,以及在一个tg内可能出现的最大周期数Nmax和最少周期数Nmin,tg根据测速目标的最大加速度而定,一个时间段长度tg内可能的最多周期数Nmax和最少周期数Nmin通过估算目标运动的最小速度和最大速度确定;
(2)调用存盘的测速多普勒信号,并对该信号进行限幅处理;
(3)求出测速多普勒信号总时间长度ta内所含时间段长度tg的数量M;
(4)设k为调用频率解调处理子程序的循环次数变量,即1≤k≤M,对k从1递增到M,逐次调用频率解调处理子程序进行如下计算和判断:
(4.1)计算当前时间段长度tg内每个信号周期的长度tj,并计数当前时间段长度tg内包括的总周期数N,j=1~N;
(4.2)比较该时间段内的总周期数N与可能的最多周期数Nmax和最少周期数Nmin,在Nmin≤N≤Nmax的情况下,统计当前时间段长度tg内全部信号周期长度tj的分布,并从中求出最频繁的信号周期长度tm,该周期长度对应的周期数记为Nm,1≤m≤N;
(4.3)若Nm>1,以周期长度接近tm的所有信号周期的统计平均值为当前时间段长度tg的周期长度td,即设定允差δt,统计周期长度在tm-δt至tm+δt间的所有周期的平均值;若Nm=1,则将N个信号周期长度tj的平均值作为当前时间段长度tg的周期长度td;
(4.4)在N<Nmin或N>Nmax的情况下,若k=1,则令周期长度td为测速目标在静止状态下多普勒信号的周期长度;若k>1,则将第k-1次频率解调计算子程序运行中解算出的周期长度td作为当前时间段长度tg的周期长度td;
(5)将当前时间段tg内的周期长度td转换为当前频率值;
(6)将信号的当前频率值转换为测速目标的当前速度值;
(7)重复第(4)~(6)步,直到k=M。
2.根据权利要求1所述的激光多普勒测速仪的多普勒信号周期分布解调方法,其特征在于:所述的多普勒信号周期的计数时间段长度tg为0.5μs。
3.根据权利要求1所述的激光多普勒测速仪的多普勒信号周期分布解调方法,其特征在于:δt的设定值为tm的5%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010160189 CN102236092B (zh) | 2010-04-29 | 2010-04-29 | 激光多普勒测速仪的多普勒信号周期分布解调方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010160189 CN102236092B (zh) | 2010-04-29 | 2010-04-29 | 激光多普勒测速仪的多普勒信号周期分布解调方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102236092A true CN102236092A (zh) | 2011-11-09 |
CN102236092B CN102236092B (zh) | 2013-03-06 |
Family
ID=44886928
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010160189 Expired - Fee Related CN102236092B (zh) | 2010-04-29 | 2010-04-29 | 激光多普勒测速仪的多普勒信号周期分布解调方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102236092B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105547632A (zh) * | 2015-12-09 | 2016-05-04 | 安徽省计量科学研究院 | 一种使用冲击速度测量装置检定或校准冲击速度的新方法 |
CN107796510A (zh) * | 2017-09-22 | 2018-03-13 | 中国第汽车股份有限公司 | 一种汽车加速行驶车外噪声测量方法 |
CN110361752A (zh) * | 2019-07-03 | 2019-10-22 | 浙江大学 | 一种基于外差干涉的运动目标运动特征测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2163384Y (zh) * | 1992-11-10 | 1994-04-27 | 天津大学 | 跟踪-计数型激光多普勒测速装置 |
CN1149135A (zh) * | 1995-09-06 | 1997-05-07 | 索尼磁石股份有限公司 | 激光多普勒测速仪 |
WO2002006841A2 (en) * | 2000-07-13 | 2002-01-24 | Honeywell International Inc. | Laser doppler velocimeter with high immunity to phase noise |
-
2010
- 2010-04-29 CN CN 201010160189 patent/CN102236092B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2163384Y (zh) * | 1992-11-10 | 1994-04-27 | 天津大学 | 跟踪-计数型激光多普勒测速装置 |
CN1149135A (zh) * | 1995-09-06 | 1997-05-07 | 索尼磁石股份有限公司 | 激光多普勒测速仪 |
WO2002006841A2 (en) * | 2000-07-13 | 2002-01-24 | Honeywell International Inc. | Laser doppler velocimeter with high immunity to phase noise |
Non-Patent Citations (1)
Title |
---|
刘昌文等: "计数型激光测速仪信号处理器及其应用研究", 《光学精密工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105547632A (zh) * | 2015-12-09 | 2016-05-04 | 安徽省计量科学研究院 | 一种使用冲击速度测量装置检定或校准冲击速度的新方法 |
CN107796510A (zh) * | 2017-09-22 | 2018-03-13 | 中国第汽车股份有限公司 | 一种汽车加速行驶车外噪声测量方法 |
CN110361752A (zh) * | 2019-07-03 | 2019-10-22 | 浙江大学 | 一种基于外差干涉的运动目标运动特征测量方法 |
CN110361752B (zh) * | 2019-07-03 | 2021-05-25 | 浙江大学 | 一种基于外差干涉的运动目标运动特征测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102236092B (zh) | 2013-03-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10429234B2 (en) | Distributed fiber optic acoustic detection device | |
JP7169642B2 (ja) | 光学的測定装置及び測定方法 | |
US8775114B2 (en) | Method and apparatus for determining phase sensitivity of an accelerometer based on an analysis of the harmonic components of the interference signal | |
US20230073827A1 (en) | Signal processing methods for an optical detection system | |
CN106443066B (zh) | 一种激光多普勒测振法冲击加速度测量装置及方法 | |
CN109298457B (zh) | 适用于原子干涉重力仪的振动噪声修正补偿方法 | |
CN109031340B (zh) | 一种测量物体运动速度的连续调频激光雷达装置 | |
Valiere et al. | Acoustic velocity measurements in the air by means of laser Doppler velocimetry: dynamic and frequency range limitations and signal processing improvements | |
Kim et al. | Effect of linewidth enhancement factor on fringe in a self-mixing signal and improved estimation of feedback factor in laser diode | |
CN102236092B (zh) | 激光多普勒测速仪的多普勒信号周期分布解调方法 | |
CN112505719A (zh) | 多普勒频率校正二次补偿的激光测风雷达测风方法及系统 | |
Shang et al. | Experimental study on minimum resolvable velocity for heterodyne laser Doppler vibrometry | |
CN106382979A (zh) | 一种激光测振仪阶跃响应的测量装置及方法 | |
CN108828618A (zh) | 基于等光频间隔重采样的远距离高精度测量装置及方法 | |
US5781283A (en) | Process and arrangement for the measurement of physical values of light scattering moving particles by means of a laser doppler anemometer | |
CN109031341B (zh) | 一种使用连续调频激光雷达装置的物体运动速度测量方法 | |
Le Duff et al. | Analog sensor design proposal for laser Doppler velocimetry | |
Sun et al. | Improving the accuracy of laser self-mixing interferometry for velocity measurement | |
RU107370U1 (ru) | Устройство для определения параметров движения цели | |
CN103323094B (zh) | 一种外差式激光干涉角振动测量方法 | |
CN110726995B (zh) | 激光雷达高精度测距方法及系统 | |
Wang et al. | Signal diversity for the reduction of signal dropouts and speckle noise in a laser-Doppler extensometer | |
Wang et al. | A Fiber Laser Doppler Vibrometer Based on Fringe Counting and Multi-Period Synchronous Frequency Measurement Method | |
RU2334953C2 (ru) | Способ радиоволнового измерения уровня | |
RU2527665C1 (ru) | Способ измерения сдвига фаз |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130306 Termination date: 20170429 |