CN105277947B - 一种基于fpga激光测速仪的信号处理系统及方法 - Google Patents

一种基于fpga激光测速仪的信号处理系统及方法 Download PDF

Info

Publication number
CN105277947B
CN105277947B CN201510695808.0A CN201510695808A CN105277947B CN 105277947 B CN105277947 B CN 105277947B CN 201510695808 A CN201510695808 A CN 201510695808A CN 105277947 B CN105277947 B CN 105277947B
Authority
CN
China
Prior art keywords
module
frequency
data
state
fft
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.)
Expired - Fee Related
Application number
CN201510695808.0A
Other languages
English (en)
Other versions
CN105277947A (zh
Inventor
梁文强
周健
聂晓明
樊振方
何鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201510695808.0A priority Critical patent/CN105277947B/zh
Publication of CN105277947A publication Critical patent/CN105277947A/zh
Application granted granted Critical
Publication of CN105277947B publication Critical patent/CN105277947B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/50—Systems of measurement based on relative movement of target
    • G01S17/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/491—Details of non-pulse systems
    • G01S7/4912—Receivers
    • 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/491—Details of non-pulse systems
    • G01S7/4912—Receivers
    • G01S7/4913—Circuits for detection, sampling, integration or read-out

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于FPGA激光测速仪的信号处理系统及方法，上电后由FPGA发出采样时钟给AD采样模块，FPGA接收到AD采样的数据后，对数据进行实时处理，处理结果通过USB传输给PC机存储、显示。本发明省去了数据在两个核心芯片之间的流动，能有效减少数据处理时间；减小数据暴露在外界的机会，增加数据的抗干扰能力；省去了大量储存数据的时间，最大限度的减少了信号处理时间，减小了数据处理延时，提高了实时性；自适应不同频段信号的效果，兼顾了测量精度与测量范围；采用USB的传输方式，数据传输更加稳定。

Description

一种基于FPGA激光测速仪的信号处理系统及方法
技术领域
本发明属于实时多普勒信号处理技术领域，尤其涉及一种基于FPGA激光测速仪的信号处理系统及方法。
背景技术
激光测速仪是利用激光源与目标的相对位移产生的多普勒频移来探测目标的速度信息。激光测速仪相比传统上声波测速、微波测速，具有工作波长短，发散角小，精度高，线性度高，动态响应快等优点。当激光照在运动物体上，被运动物体散射回来的光线相对于入射光线在频率上将会发生一个偏移量即多普勒频移。此频移携带有运动物体的速度信息，能否对此频移进行准确，快速的测量，直接影响着激光测速仪的核心性能。数据处理电路是激光测速仪的核心部分，用于拾取信号，并分析多普勒频移，从而解算出需要的速度信息。此电路是制约测速仪测量精度，测量范围，测量稳定性、测量速度的重要因素，是限制激光测速仪技术发展的重要一环。激光测速仪在高动态环境下测量，测量数据往往具有突发性。先将数据存储然后进行处理，能很好的解决这个问题。但是这会花大量的时间用在数据的存储与提取上，信号处理的整体时间就会拉长。数据突发性与测量实时性就存在矛盾关系，这是激光测速仪信号处理电路存在的一个问题。测量精度是激光测速仪的一个重要指标。传统上采用FPGA与DSP相结合的方式采集与处理数据，能够很好地解决数据处理的精度问题。但是数据在FPGA与DSP之间的流动会耗费大量的时间，且增加数据暴露在外界的机会，数据可靠性下降。所以测量精度与测量时间、数据可靠性之间存在矛盾，这是激光测速仪信号处理电路存在的又一问题。速度测量范围是激光测速仪的又一个重要指标。大的速度测量范围，意味着信号处理器大的频率测量范围。大的频率测量范围意味着采样频率的提高，然而测量精度在一定意义下 与采样频率是反比关系，于是测量范围与测量精度又存在着一定的矛盾，这是激光测速仪信号处理电路需要解决的第三个问题。
传统采用FPGA与DSP相结合的方式采集与处理数据，数据流动会耗费大量时间，数据可靠性下降，测量范围较小。
发明内容
本发明的目的在于提供一种基于FPGA激光测速仪的信号处理系统及方法，旨在解决传统采用FPGA与DSP相结合的方式采集与处理数据，数据流动会耗费大量时间，数据可靠性下降，测量范围较小的问题。
本发明是这样实现的，一种基于FPGA激光测速仪的信号处理方法，所述基于FPGA激光测速仪的信号处理方法包括以下步骤：
步骤一，根据反馈采样频率信息采样频率状态量，产生相应频率的时钟，经过模数转换生成数字信号；
步骤二，对数字信号进行加Hanning窗处理，并将12位的数字信号推高到14位，同时记录此时的采样频率产生相应的采样频率状态量；
步骤三，对数字信号进行实时的FFT变换，同时输出运算的状态信息，状态信息包含每次FFT转换的开始与结束信号，此开始与结束信息协调整个程序的运行，采样频率状态量要延时；
步骤四，接收频谱信息，根据FFT转换的开始与结束状态，连续产生每组频谱相应的功率谱，该步骤结束后转到步骤五；
步骤五，根据开始与结束信息，对每组功率谱进行运算，求得功率谱的中心序号峰值序号，将此信息传递；
步骤六，根据运算状态信息开始与结束信息与采样频率状态量，对每组功率谱中心序号峰值序号乘以相应的系数，得到多普勒频率；
步骤七，控制USB通信电路，将多普勒频率传递出去；
步骤八，通过LABVIEW程序接收、校验、存储、显示多普勒信息；
步骤九，频率反馈函数根据得到的FFT转换的开始与结束信息，由反馈的功率谱中心序号决定下一组的数据采样频率，由提供的多普勒频率决定下一运算周期的采样频率，输出下一运算周期的采样频率状态量；
步骤十，将反馈采样频率信息采样频率状态量进行缓存，根据读取的开始与结束信息，当一次组数据开始采集时，释放反馈采样频率信息的采样频率状态量，结束后转到步骤一至此一次循环运算周期结束。
本发明的另一目的在于提供一种所述的基于FPGA激光测速仪的信号处理方法的信号处理系统，所述信号处理系统包括：
降压偏置模块，由两片运算放大器构成，用于对输入的多普勒模拟信号进行降压与偏执偏置，使信号电压满足AD转换芯片的输入范围1-3v；采用模拟电子学的乘法器与加法器电路；
AD采样模块，与降压偏置模块连接，由一片数模转换芯片构成，对经过降压后的模拟信号进行模数转换，得到相应的数字信号给FPGA处理；由电容电阻简单配置后，形成单端输入模式，输出12位的数字信号，最高采样频率为65M，电压输入范围为1V-3V；
USB通信模块，由一块USB控制芯片组成，屏蔽复杂的USB通信协议，用户只需要对芯片进行简单的配置即可完成USB通信任务；
PC机接收模块，与USB通信模块连接，接收USB传输的数据，进行校验、储存，显示；
FPGA控制模块，与AD采样模块和USB通信模块连接，用于实现组合逻辑功能又可实现时序逻辑。
进一步，所述降压偏置模块采用模拟电子学的乘法器与加法器电路，具体连接如下：
电阻R1连接运算放大器U1的正极，电阻R2与电阻R1并联，与运算放大器U1的输出端连接，电阻R3连接放大器U1的负极，运算放大器U1的正极和负极之间并联有两个二极管，电阻R4与其中一个二极管串联；电阻R3连接运 算放大器U2的输出端，电阻R5与电阻R3并联，并连接运算放大器U2的负极，运算放大器U2的的正极连接GND端；电阻R6与电阻R5并联；电阻R6与运算放大器U3的输出端和负极连接，运算放大器U3的正极加+2V电压；
运算放大器U1、运算放大器U2和运算放大器U3都加载+5V的双电源；电阻R1＝2K，电阻R2＝2K，电阻R3＝18K，电阻R4＝2K，电阻R5＝10K，电阻R6＝10K。
进一步，所述FPGA控制模块进一步包括：
采样频率产生模块，输入为频率缓存模块的反馈采样频率信号采样频率状态量，输出时钟信号给AD采样模块；
窗函数模块，通过14位的数字信号和采样频率状态量与FFT模块相连，输入为AD采样模块生成的12位的数字信号，输出为14位的经过截取的数字信号和反应此刻采样频率的状态量，此状态量将跟随此运算周期的数据经过各模块；
FFT模块，输入为窗函数截取的14位数字信号，输出为信号的频谱数据与运算、表明运算开始与结束时刻的运算状态信息；
功率谱模块，通过频谱数据、FFT运算状态信息和采样频率状态量与FFT模块相连，通过功率谱和采样频率状态量与峰值序号模块相连，用于输入为FFT模块的频谱数据与运算、运算状态信息，输出为信号的功率谱数据；
中心序号峰值序号模块，读取FFT运算状态信息与FFT模块相连，输出功率谱峰值序号和采样频率状态量与系数模块连接；输入为功率谱模块的功率谱数据与FFT模块运算状态信息，输出为功率谱幅值最高处的中心峰值序号；
系数模块，通过读取FFT运算状态信息与FFT模块连接，输出多普勒频率与USB控制输出模块和频率反馈模块连接；输入为功率谱中心序号峰值序号与FFT模块的运算状态信息，输出为多普勒频率；
USB控制模块，输入为系数模块的多普勒频率数据，输出为USB控制信号；
频率反馈模块，通过功率谱峰值序号连接峰值序号模块，通过状态信息与 FFT模块连接，输出反馈采样频率信息采样频率状态量与频率缓存模块连接；输入为中心序号模块的功率谱中心序号系数模块的多普勒频率与FFT模块运算状态信息，输出为将要采取下一周期的采样频率信息状态量；
频率缓存模块，通过状态信息与FFT模块连接，输出反馈采样频率与频率产生模块连接；输入为FFT模块运算状态信息与频率反馈模块的反馈采样频率数据采样频率状态量，输出为缓存后的反馈采样频率状态量。
进一步，所述系数模块中系数即为频谱分辨率Δf＝fs/N，其中N是确定的，根据与峰值序号同时到达的采样频率状态量，知道此周期峰值序号对应采样频率fs，根据此采样频率，系数模块为功率谱峰值序号乘以相应的系数，得到多普勒频率。
本发明的另一目的在于提供一种使用所述基于FPGA激光测速仪的信号处理方法的板材在线切割监控控制系统。
本发明的另一目的在于提供一种使用所述基于FPGA激光测速仪的信号处理方法的管材在线切割监控控制系统。
本发明的另一目的在于提供一种使用所述基于FPGA激光测速仪的信号处理方法的电缆长度速度测量装置。
本发明的另一目的在于提供一种使用所述基于FPGA激光测速仪的信号处理方法的砂纸长度速度测量装置。
本发明提供的基于FPGA激光测速仪的信号处理系统及方法，
与现有技术相比，具有以下优势：
相对于现有技术，本发明的有益效果为：
1、本发明取代传统上FPGA与DSP相结合的方式，完全通过FPGA实现数据采集和数据处理的功能，省去了数据在两个核心芯片之间的流动，能有效减少数据处理时间；减小数据暴露在外界的机会，增加数据的抗干扰能力。
2、本发明的程序采用流水线方式设计，从采集到最终完成传输，数据即刻产生即刻处理，几乎没有任何的停顿，省去了大量储存数据的时间，最大限度 的减少了信号处理时间，减小了数据处理延时，提高了实时性。
3、本发明由输出结果反馈来调整系统采样与处理的频率，已到达自适应不同频段信号的效果，兼顾了测量精度与测量范围。
4、本发明采用USB的传输方式，数据传输更加稳定。
附图说明
图1是本发明实施例提供的基于FPGA激光测速仪的信号处理系统结构示意图；
图中：1、降压偏置模块；2、AD采样模块；3、USB通信模块；4、PC机接收模块；5、FPGA控制模块。
图2是本发明实施例提供的基于FPGA激光测速仪的信号处理方法流程图。
图3是本发明实施例提供的降压偏置电路图。
图4是本发明实施例提供的多普勒功率谱示意图。
图5是本发明实施例提供的数据延迟与数据更新时间示意图。
图6是本发明实施例提供的过滤畸变数据示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明适用的激光测速仪的速度测量范围为0.16m/s～6.1m/s，量程比达到38:1(量程比即最大测量值与最小测量值之比)，该速度范围与532nm波长的激光配合对应的多普勒频率范围为0.6MHz～23MHz，比例系数为2.66×10-3。频率的测量精度测量准确度小于优于8×10-4，测量稳定度优于2.5×10-7，数据处理延迟小于1230us，数据处理结果更新周期小于410us。激光测速仪可以在板材、管材在线切割监控，电缆或砂纸长度速度测量等；也适用于绒布、毛皮等纺织品、涂层或粘胶表面、泡沫橡胶表面物体的测速测长得到广泛的应用。
下面结合附图对本发明的应用原理作详细的描述。
如图1所示，本发明实施例的基于FPGA激光测速仪的信号处理系统主要包括：降压偏置模块1、AD采样模块2、USB通信模块3、PC机接收模块4、FPGA控制模块5。
降压偏置模块1，由两片运算放大器构成，用于对输入的多普勒模拟信号进行降压与偏执偏置，使信号电压满足AD转换芯片的输入范围1-3v；采用模拟电子学的乘法器与加法器电路，降压偏置模块的电连接示意图，如图3所示；
电阻R1连接运算放大器U1的正极，电阻R2与电阻R1并联，与运算放大器U1的输出端连接，电阻R3连接放大器U1的负极，运算放大器U1的正极和负极之间并联有两个二极管，电阻R4与其中一个二极管串联；电阻R3连接运算放大器U2的输出端，电阻R5与电阻R3并联，并连接运算放大器U2的负极，运算放大器U2的的正极连接GND端；电阻R6与电阻R5并联；电阻R6与运算放大器U3的输出端和负极连接，运算放大器U3的正极加+2V电压。
运算放大器U1、运算放大器U2和运算放大器U3都加载+5V的双电源；电阻R1＝2K，电阻R2＝2K，电阻R3＝18K，电阻R4＝2K，电阻R5＝10K，电阻R6＝10K。
AD采样模块2，与降压偏置模块1连接，由一片数模转换芯片构成，对经过降压后的模拟信号进行模数转换，得到相应的数字信号给FPGA处理；由电容电阻简单配置后，形成单端输入模式，输出12位的数字信号，最高采样频率为65M，电压输入范围为1V-3V。
USB通信模块3，由一块USB控制芯片组成，屏蔽了复杂的USB通信协议，用户只需要对芯片进行简单的配置即可完成USB通信任务。
PC机接收模块4，与USB通信模块3连接，接收USB传输的数据，对其进行校验、储存，显示。
FPGA控制模块5，与AD采样模块2和USB通信模块3连接，用于实现 组合逻辑功能又可实现时序逻辑。
FPGA控制模块5进一步包括：
采样频率产生模块，输入为频率缓存模块的反馈采样频率信号采样频率状态量，输出为合适的时钟信号给AD采样模块。
窗函数模块，输入为AD采样模块生成的12位的数字信号，输出为14位的经过截取的数字信号和反应此刻采样频率的状态量，此状态量将跟随此运算周期的数据经过各模块。
FFT模块，输入为窗函数截取的14位数字信号，输出为信号的频谱数据与运算、表明运算开始与结束时刻的运算状态信息。
功率谱模块，输入为FFT模块的频谱数据与运算、运算状态信息，输出为信号的功率谱数据。
中心序号峰值序号模块，输入为功率谱模块的功率谱数据与FFT模块运算状态信息，输出为功率谱幅值最高处的中心峰值序号。
系数模块，输入为功率谱中心序号峰值序号与FFT模块的运算状态信息，输出为多普勒频率。
USB控制模块，输入为系数模块的多普勒频率数据，输出为USB控制信号。
频率反馈模块，输入为中心序号模块的功率谱中心序号系数模块的多普勒频率与FFT模块运算状态信息，输出为将要采取下一周期的采样频率信息状态量。
频率缓存模块，输入为FFT模块运算状态信息与频率反馈模块的反馈采样频率数据采样频率状态量，输出为缓存后的反馈采样频率状态量。
FPGA内部各模块的连接关系：窗函数模块通过14位的数字信号和采样频率状态量与FFT模块相连；功率谱模块通过频谱数据、FFT运算状态信息和采样频率状态量与FFT模块相连，通过功率谱和采样频率状态量与峰值序号模块相连；峰值序号模块读取FFT运算状态信息与FFT模块相连，输出功率谱峰值序号和采样频率状态量与系数模块连接；系数模块通过读取FFT运算状态信息 与FFT模块连接，输出多普勒频率与USB控制输出模块和频率反馈模块连接；频率反馈模块通过功率谱峰值序号连接峰值序号模块，通过状态信息与FFT模块连接，输出反馈采样频率信息采样频率状态量与频率缓存模块连接；频率缓存模块通过状态信息与FFT模块连接，输出反馈采样频率与频率产生模块连接。
如图2所示，本发明实施例的基于FPGA激光测速仪的信号处理方法包括以下步骤：
步骤一，采样频率产生模块根据反馈采样频率信息采样频率状态量，产生相应频率的时钟，控制AD采样电路，AD采样电路经过模数转换生成数字信号，该步骤结束后转到步骤二；
步骤二，窗函数模块对数字信号进行加Hanning窗处理，并将12位的数字信号推高到14位，同时记录此时的采样频率产生相应的采样频率状态量，此状态量在以下的程序中只做相应的延时处理，伴随本周期的数据通过各函数，无特殊情况以下不做累述，该步骤结束后转到步骤三；
步骤三，FFT模块对数字信号进行实时的FFT变换，产生数字信号的频谱给功率谱模块，同时输出运算的状态信息，该状态信息包含每次FFT转换的开始与结束信号，此开始与结束信息协调整个程序的运行，由于FFT模块从采集数据到结果输出有两个运算周期的延时，所以采样频率状态量也要做相应的延时，该步骤结束后转到步骤四；
步骤四，功率谱模块接收频谱信息，根据FFT模块的开始与结束状态，连续产生每组频谱相应的功率谱，该步骤结束后转到步骤五；
步骤五，中心序号峰值序号模块根据开始与结束信息，对每组功率谱进行运算，求得功率谱的中心序号峰值序号，将此信息传递给系数模块与频率反馈模块，该步骤结束后转到步骤六和步骤七；
步骤六，根据运算状态信息开始与结束信息与采样频率状态量，系数模块对每组功率谱中心序号峰值序号乘以相应的系数，得到多普勒频率，该步骤结束后转到步骤步骤七和步骤八；
步骤七，USB控制模块控制USB通信电路，将多普勒频率传递出去。该步骤结束后转到步骤步骤八；
步骤八，PC机通过LABVIEW程序接收、校验、存储、显示多普勒信息；
步骤九，频率反馈函数根据得到的FFT模块的开始与结束信息，由中心序号模块反馈的功率谱中心序号决定下一组的数据采样频率，由系数模块提供的多普勒频率决定下一运算周期的采样频率，输出下一运算周期的采样频率状态量，该步骤结束后转到步骤十；
步骤十，频率缓存模块将反馈采样频率信息采样频率状态量进行缓存，根据读取的开始与结束信息，当一次组数据开始采集时，释放反馈采样频率信息采样频率状态量给采样频率产生采样频率产生模块，该步骤结束后转到步骤一至此一次循环运算周期结束。
下面附图对本发明的方法做进一步的说明。
(1)窗函数模块实现四个功能，将数据位推高；对数据进行加窗截取；记录采样频率；除去因为采样频率变化引起信号畸变。由于FPGA无法进行浮点数运算，所以将数据位数推高，低位看作是小数位，高位看作是整数位，有利于降低“有限字长效应”带来的影响，有利于提高运算准确度；由于硬件数据处理能力有限，只能对有限的数据进行FFT变换，于是要对采样得到的数据进行点数为N(2048)的截取，既对数据进行加窗。然而时域数据的加窗操作会带来频域内的频谱泄露，泄漏后的频谱如图4所示。本发明采用Hanning窗对数据进行截取。Hanning窗具有良好的抑制旁瓣的作用，将能量集中在主瓣中，抑制频谱泄露带来的影响；记录此周期数据的采样频率，输出采样频率状态量，此状态量只经过一些必要的延时处理，将始终跟随此周期数据经过各函数，如无特殊情况以下不再累述。
由于本发明的采样频率是随着信号频率自适应变化。采样频率在一个信号处理周期内是固定的。然而在相邻两周期的交界处，采样频率有可能发生改变。在采样频率变化时，由于AD采样电路有一定的延时，不能随着程序采样频率 的变化而立即变化。在频率变化的衔接处会出现采样数据的畸变。本窗函数将此畸变数据进行了截取，保证参与运算的数据都是正常的。过滤畸变数据如图6所示
(2)FFT模块，采用现有的IP核函数，采用基-4算法，双FFT引擎，数据位宽为14位，旋转因子位宽为14位，完全使用数据流结构。FFT运算长度为N(2048)个点，FFT核函数采样数据需要N个时钟周期，数据运算需要N个时钟周期，运算结果输出需要N个时钟周期。每N个时钟周期为一个运算周期，即每一部分需要一个运算周期。每次解算都包含这三个部分，不同的解算进程并行进行。所以FFT从采样到结果完全输出需要3个运算周期(数据延迟时间)，由于采样与运算可以并行运行，运算结果的更新时间只要一个运算周期。数据更新与延迟时间如图5所示。
结果输出部分是N个点的14位复数结果，每个时钟周期输出一个点。前N/2个点已经包含了所关心的全部信息，所以本发明利用后N/2个周期完成此组数据的后续所有运算。后期处理并没有花费额外的周期，数据延迟时间为3N个时钟周期，数据更新时间为N个周期。本发明时钟周期变化范围为20ns—200ns,故数据处理时间为123—1230us,数据更新时间为41us—410us。由于数据连续输入，结果连续输出，为了表明每个周期的开始与结束时刻，FFT模块输出运算状态信息量。
(3)功率谱模块。由FFT函数得到数据的频谱是复数形式的，分为实部与虚部两部分。将实部谱的各点数值的平方与虚部谱相应各点的平方相加就得到数据的功率谱。多普勒信号由于加窗的影响，功率谱发生了泄露，功率谱不再是两根竖线，而是在原有功率谱处有一定的展宽。多普勒功率谱如图4所示。
(4)峰值序号模块。得到功率谱后，需要检测功率谱峰值最大处的序号。功率谱输出为N(2048)点的结果，根据FFT模块给出的运算状态信息，可以识别功率谱的开始与结束时刻，前N/2是正频率点，包含所有关心的信息。利用后N/2的负频率点的输出时间，便可在新的一轮结果输出之前计算得到所关心的功 率谱最大峰值处的序号。
由于“栅栏效应”的影响，直接求取功率谱峰值序号会出现偏差，需要对此偏差进行频谱校正。在本发明中运用了能量重心法进行频谱校正。考虑到速度性和硬件资源有限，采用谱峰值附近五点参与能量中心校正。功率谱幅值为yk，k为序号，H所求经过校正的序号。能量重心法频率校正公式如式(1)所示。
H＝k(yk/E)+(k+1)(yk+1/E)+(k+2)(yk+2/E)
(1)
＝(kyk+(k+1)yk+1+(k+2)yk+2)/E
(5)系数模块。由峰值序号模块得到的序号乘以相应的系数就可以得到多普勒频率。其中系数即为频谱分辨率Δf＝fs/N。其中N是确定的，由于不同的采样频率fs，系数是不同的。根据与峰值序号同时到达的采样频率状态量，可以知道此周期峰值序号对应采样频率fs。根据此采样频率，系数模块为功率谱峰值序号乘以相应的系数，就得到多普勒频率。
(6)频率反馈模块。根据当前解算得到的的多普勒频率，决定下一组数据的采样频率。由数字信号处理的相关知识，频谱分辨率Δf＝fs/N。N是采样点数。要提高分辨率Δf，就要尽量减小fs，尽量增大采样点数N。由奈奎斯特定理，AD采样频率要大于多普勒频率fd的两倍，即fs≥2fd，为减小频谱混叠的影响，经过实验fs＝3fd～5fd比较合适。由于激光测速仪所用的场合不同，本发明将测速仪的测量范围适当的放大，速度测量范围为0.16m/s～6.1m/s，量程比达到38:1，对应的频率测量范围为0.6MHZ-23MHZ。为保证频率测量准确度优于8×10-4，测量稳定度优于2.5×10-7，采用频率分段测量。根据得到的多普勒频率，调整下一运算周期数据的采样频率，输出下一运算周期的采样频率状态量。
(7)频率缓存模块。由于FFT采样数据需要一个完整的运算周期，一个运算周期的数据采样频率要恒定，所以本模块对反馈的采样频率状态量进行缓存。根据FFT模块的运算状态量，可以识别每一个运算周期的始末，在一个运算周期的开始之际(数据采样的开始之际)释放采样频率状态量。该模块取代传统上的FIFO函数，运用更加灵活，速度更加快。
(8)采样频率产生模块。根据采样频率状态量，从PLL已经产生的多个时钟中挑选出一个合适的时钟频率作为采样时钟。此时钟控制AD采样电路完成采样。
(9)USB控制模块。由于USB控制芯片CH376屏蔽了底层的USB通信协议，用户只需要对芯片进行一些简单的读写操作即可完成USB传输。此模块将得到的多普勒频率加上必要的校验信息，然后发送出去。
(10)PC机接收模块。此模块由LabVIEW软件编写，对接收到的数据进行校验、提取、储存、显示。
本发明的工作原理：
外部多普勒信号经过降压模拟降压偏置模块后输入到AD采样模块，FPGA的I/O引脚输出时钟信号到AD采样模块的时钟输入端，AD采样模块12位数字信号输出引脚接到FPGA的I/O引脚。FPGA的输入输出数据接口与控制信号接口分别连接到CH376芯片的数据端与控制端。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种基于FPGA激光测速仪的信号处理方法，其特征在于，所述基于FPGA激光测速仪的信号处理方法包括以下步骤：
步骤一，根据反馈采样频率信息采样频率状态量，产生相应频率的时钟，经过模数转换生成数字信号；
步骤二，对数字信号进行加Hanning窗处理，并将12位的数字信号推高到14位，同时记录此时的采样频率产生相应的采样频率状态量；
步骤三，对数字信号进行实时的FFT变换，同时输出运算的状态信息，状态信息包含每次FFT转换的开始与结束信号，此开始与结束信息协调整个程序的运行，采样频率状态量要延时；
步骤四，接收频谱信息，根据FFT转换的开始与结束状态，连续产生每组频谱相应的功率谱，该步骤结束后转到步骤五；
步骤五，根据开始与结束信息，对每组功率谱进行运算，求得功率谱的中心序号峰值序号，将此信息传递；
步骤六，根据运算状态信息开始与结束信息与采样频率状态量，对每组功率谱中心序号峰值序号乘以相应的系数，得到多普勒频率；
步骤七，控制USB通信电路，将多普勒频率传递出去；
步骤八，通过LABVIEW程序接收、校验、存储、显示多普勒信息；
步骤九，频率反馈函数根据得到的FFT转换的开始与结束信息，由反馈的功率谱中心序号决定下一组的数据采样频率，由提供的多普勒频率决定下一运算周期的采样频率，输出下一运算周期的采样频率状态量；
步骤十，将反馈采样频率信息采样频率状态量进行缓存，根据读取的开始与结束信息，当一次组数据开始采集时，释放反馈采样频率信息的采样频率状态量，结束后转到步骤一至此一次循环运算周期结束。
2.一种如权利要求1所述的基于FPGA激光测速仪的信号处理方法的信号处理系统，其特征在于，所述信号处理系统包括：
降压偏置模块，由两片运算放大器构成，用于对输入的多普勒模拟信号进行降压与偏执偏置，使信号电压满足AD转换芯片的输入范围1-3v；采用模拟电子学的乘法器与加法器电路；
AD采样模块，与降压偏置模块连接，由一片数模转换芯片构成，对经过降压后的模拟信号进行模数转换，得到相应的数字信号给FPGA处理；由电容电阻简单配置后，形成单端输入模式，输出12位的数字信号，最高采样频率为65M，电压输入范围为1V-3V；
USB通信模块，由一块USB控制芯片组成，屏蔽复杂的USB通信协议，用户只需要对芯片进行简单的配置即可完成USB通信任务；
PC机接收模块，与USB通信模块连接，接收USB传输的数据，进行校验、储存，显示；
FPGA控制模块，与AD采样模块和USB通信模块连接，用于实现组合逻辑功能又可实现时序逻辑。
3.如权利要求2所述的信号处理系统，其特征在于，所述降压偏置模块采用模拟电子学的乘法器与加法器电路，具体连接如下：
电阻R1连接运算放大器U1的正极，电阻R2与电阻R1并联，与运算放大器U1的输出端连接，电阻R3连接放大器U1的负极，运算放大器U1的正极和负极之间并联有两个二极管，电阻R4与其中一个二极管串联；电阻R3连接运算放大器U2的输出端，电阻R5与电阻R3并联，并连接运算放大器U2的负极，运算放大器U2的的正极连接GND端；电阻R6与电阻R5并联；电阻R6与运算放大器U3的输出端和负极连接，运算放大器U3的正极加+2V电压；
运算放大器U1、运算放大器U2和运算放大器U3都加载±5V的双电源；电阻R1＝2K，电阻R2＝2K，电阻R3＝18K，电阻R4＝2K，电阻R5＝10K，电阻R6＝10K。
4.如权利要求2所述的信号处理系统，其特征在于，所述FPGA控制模块进一步包括：
采样频率产生模块，输入为频率缓存模块的反馈采样频率信号采样频率状态量，输出时钟信号给AD采样模块；
窗函数模块，通过14位的数字信号和采样频率状态量与FFT模块相连，输入为AD采样模块生成的12位的数字信号，输出为14位的经过截取的数字信号和反应此刻采样频率的状态量，此状态量将跟随此运算周期的数据经过各模块；
FFT模块，输入为窗函数截取的14位数字信号，输出为信号的频谱数据与运算、表明运算开始与结束时刻的运算状态信息；
功率谱模块，通过频谱数据、FFT运算状态信息和采样频率状态量与FFT模块相连，通过功率谱和采样频率状态量与峰值序号模块相连，用于输入为FFT模块的频谱数据与运算、运算状态信息，输出为信号的功率谱数据；
中心序号峰值序号模块，读取FFT运算状态信息与FFT模块相连，输出功率谱峰值序号和采样频率状态量与系数模块连接；输入为功率谱模块的功率谱数据与FFT模块运算状态信息，输出为功率谱幅值最高处的中心峰值序号；
系数模块，通过读取FFT运算状态信息与FFT模块连接，输出多普勒频率与USB控制输出模块和频率反馈模块连接；输入为功率谱中心序号峰值序号与FFT模块的运算状态信息，输出为多普勒频率；
USB控制模块，输入为系数模块的多普勒频率数据，输出为USB控制信号；
频率反馈模块，通过功率谱峰值序号连接峰值序号模块，通过状态信息与FFT模块连接，输出反馈采样频率信息采样频率状态量与频率缓存模块连接；输入为中心序号模块的功率谱中心序号系数模块的多普勒频率与FFT模块运算状态信息，输出为将要采取下一周期的采样频率信息状态量；
频率缓存模块，通过状态信息与FFT模块连接，输出反馈采样频率与频率产生模块连接；输入为FFT模块运算状态信息与频率反馈模块的反馈采样频率数据采样频率状态量，输出为缓存后的反馈采样频率状态量。
5.如权利要求4所述的信号处理系统，其特征在于，所述系数模块中系数即为频谱分辨率Δf＝fs/N，其中N为采样点数2048，根据与峰值序号同时到达的采样频率状态量，知道此周期峰值序号对应采样频率fs，根据此采样频率，系数模块为功率谱峰值序号乘以相应的系数，得到多普勒频率。
6.一种使用权利要求1所述基于FPGA激光测速仪的信号处理方法的板材在线切割监控控制系统。
7.一种使用权利要求1所述基于FPGA激光测速仪的信号处理方法的管材在线切割监控控制系统。
8.一种使用权利要求1所述基于FPGA激光测速仪的信号处理方法的电缆长度速度测量装置。
9.一种使用权利要求1所述基于FPGA激光测速仪的信号处理方法的砂纸长度速度测量装置。
CN201510695808.0A 2015-10-22 2015-10-22 一种基于fpga激光测速仪的信号处理系统及方法 Expired - Fee Related CN105277947B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510695808.0A CN105277947B (zh) 2015-10-22 2015-10-22 一种基于fpga激光测速仪的信号处理系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510695808.0A CN105277947B (zh) 2015-10-22 2015-10-22 一种基于fpga激光测速仪的信号处理系统及方法

Publications (2)

Publication Number Publication Date
CN105277947A CN105277947A (zh) 2016-01-27
CN105277947B true CN105277947B (zh) 2017-12-19

Family

ID=55147289

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510695808.0A Expired - Fee Related CN105277947B (zh) 2015-10-22 2015-10-22 一种基于fpga激光测速仪的信号处理系统及方法

Country Status (1)

Country Link
CN (1) CN105277947B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106093961A (zh) * 2016-05-30 2016-11-09 天水师范学院 一种基于fpga激光测速仪的信号处理系统及方法
CN107715298A (zh) * 2017-11-16 2018-02-23 陈敏 一种多功能妇产科镇痛电子治疗仪
CN107874337A (zh) * 2017-11-25 2018-04-06 蒋鹏宇 一种防雾霾口罩
CN108317737A (zh) * 2018-02-21 2018-07-24 付志昊 一种家用智能热水器系统
CN109460906A (zh) * 2018-10-26 2019-03-12 赵先 一种新能源汽车充换电管理平台的设计方法
CN110488676A (zh) * 2019-07-19 2019-11-22 刘训祥 一种中医推拿手法力度数据化控制系统及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101320086A (zh) * 2008-06-27 2008-12-10 北京航空航天大学 一种多普勒测速激光雷达的回波信号处理装置和方法
CN102736074A (zh) * 2012-06-25 2012-10-17 中国科学院上海技术物理研究所 基于fpga的激光多普勒雷达信号处理器及处理方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008127335A1 (en) * 2007-04-13 2008-10-23 Raytheon Company Spotlight synthetic aperture radar system using fpga and corresponding method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101320086A (zh) * 2008-06-27 2008-12-10 北京航空航天大学 一种多普勒测速激光雷达的回波信号处理装置和方法
CN102736074A (zh) * 2012-06-25 2012-10-17 中国科学院上海技术物理研究所 基于fpga的激光多普勒雷达信号处理器及处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Design and Implementation of Signal Acquisition System for Laser Doppler;Bing Yang et al.;《2013 Third International Conference on Instrumentation, Measurement, Computer, Communication and Control》;20140623;第1237-1241页 *
用于车载自主惯导系统的激光多普勒测速仪;周健 等;《中国激光》;20100510;第37卷(第5期);第1316-1321页 *

Also Published As

Publication number Publication date
CN105277947A (zh) 2016-01-27

Similar Documents

Publication Publication Date Title
CN105277947B (zh) 一种基于fpga激光测速仪的信号处理系统及方法
CN104181577B (zh) 一种基于全数字化技术的束流位置和相位测量系统及方法
CN106093961A (zh) 一种基于fpga激光测速仪的信号处理系统及方法
CN104075802B (zh) 一种高动态范围的光子计数型微弱光信号测量装置及方法
CN106502308B (zh) 一种超声脉冲激励信号的波形产生系统及产生方法
CN107145468A (zh) 一种信号发生与读出装置及控制方法
CN206193081U (zh) 基于fpga的信号源、示波器、万用表一体机
CN102435847A (zh) 高精度电参量分析仪
CN101672692B (zh) 一种基于虚拟仪器的音叉谐振频率快速测量方法
CN104713575A (zh) 一种闭环光纤陀螺频率特性的测试方法
CN212060395U (zh) 一种支持实时采样和等效采样的数字示波器
CN103075981B (zh) 一种超声波测厚方法
CN108445509A (zh) 基于gpu的相干激光雷达信号处理方法
CN103499386B (zh) 一种精密测量频率的方法和装置
CN205091464U (zh) 一种基于fpga激光测速仪的信号处理系统
CN107664962A (zh) 基于fpga的光纤传感系统的集成型控制方法及装置
CN102829940A (zh) 针对卫星飞轮的扰动仿真的实现方法
CN201622352U (zh) α辐射核信号模拟装置
CN107036704A (zh) 一种基于周期非均匀采样技术的分布式光纤传感系统
CN204594436U (zh) 一种超声波厚度测量的装置
CN201622351U (zh) 核信号随机特性模拟器
CN204241102U (zh) 基于dds的变压器振动监测校验平台
CN102944255B (zh) 数字化光纤光栅解调仪表及方法
CN102426865A (zh) 裂变室输出信号的数字化处理系统及方法
CN202421332U (zh) 高精度电参量分析仪

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
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: 20171219

Termination date: 20181022