CN106019289B - 宽带相控声学多普勒测速系统信号处理方法 - Google Patents

宽带相控声学多普勒测速系统信号处理方法 Download PDF

Info

Publication number
CN106019289B
CN106019289B CN201610338300.XA CN201610338300A CN106019289B CN 106019289 B CN106019289 B CN 106019289B CN 201610338300 A CN201610338300 A CN 201610338300A CN 106019289 B CN106019289 B CN 106019289B
Authority
CN
China
Prior art keywords
signal
sequence
frequency spectrum
calculates
complex
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
Application number
CN201610338300.XA
Other languages
English (en)
Other versions
CN106019289A (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.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi 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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201610338300.XA priority Critical patent/CN106019289B/zh
Publication of CN106019289A publication Critical patent/CN106019289A/zh
Application granted granted Critical
Publication of CN106019289B publication Critical patent/CN106019289B/zh
Active 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
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/50—Systems of measurement, based on relative movement of the target
    • G01S15/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S15/588—Velocity or trajectory determination systems; Sense-of-movement determination systems measuring the velocity vector
    • 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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种宽带相控声学多普勒测速系统的信号处理方法。现有方法实现结构较为复杂，且须对各路信号并行处理，硬件资源开销极大。本发明首先采用基于频域的频谱搬移与正交化处理方式，完成回波信号的正交解调，并构造出一种复包络信号波束合成方法，并行完成信号实部与虚部的波束合成，最后利用复相关技术计算出信号相角，从而推导出所测速度。本发明在保证性能的前提下简化了信号处理流程，实现架构较为简单，较大程度地降低了系统的处理复杂度和硬件资源的规模。

Description

宽带相控声学多普勒测速系统信号处理方法
技术领域
本发明属声学多普勒测速信号处理技术领域，具体涉及一种宽带相控声学多普勒测速系统的信号处理方法，用于声学多普勒水流剖面仪和计程仪等相对于水层或水底运动速度信息的提取。
背景技术
声学多普勒测速系统利用水声换能器向水介质发射一个定向的声脉冲信号，然后接收从海水中的散射体上反射的回波信号，利用声传播路径的时间差异转换成不同的深度单元，测量各层回波信号的频率信息，分析多普勒频移变化，换算得到相对运动速度。目前，采用伪随机序列进行相位编码调制技术的宽带相控声学多普勒测速系统，因具有测量精度高、范围广、抗干扰能力强等优点，受到了日益广泛的应用。
传统的宽带相控声学多普勒测速系统对水声换能器接收到的多路回波进行采样后的信号处理方法为：首先通过频谱搬移及低通滤波器进行基带解调，然后对解调后的多路基带信号进行波束合成，对合成后的信号进行希尔伯特变换以获取复信号，再采用复相关技术计算出回波信号中的多普勒频偏，最后根据多普勒频移值计算出所需测量的相对运算速度。该信号处理方法中的频谱搬移、低通滤波、希尔伯特变换等环节实现结构较为复杂，硬件资源开销极大，且须对各路信号并行处理，进一步增大了硬件资源的规模。
发明内容
本发明的目的在于解决上述方法的缺点和不足，提出一套基于频域处理以及复包络波束合成的宽带相控声学多普勒测速系统信号处理方法，在保证测量性能的前提下简化了信号处理流程，实现架构较为简单，较大程度地降低了处理复杂度并减少了硬件资源的开销。
设系统宽带相控声学多普勒测速系统中进行相位编码调制的伪随机序列周期为Tm，相控换能器发射的声脉冲信号的波束张角为θ，波长为λ，声脉冲信号的载波频率为fc，系统处理的最大多普勒频移绝对值为△f，所述的处理方法处理步骤如下：
步骤(1)对相控换能器阵接收到的四路回波信号进行采样，形成四路数字信号。采样频率fs，采样后的信号为xi(n)，其中i＝1,2,3,4表示信号对应的路数，n＝1,2,3,…,N表示采样点序号，N表示采样后数字序列的长度。
步骤(2)对每路信号各自按样本重叠50％的方式做FFT(快速傅里叶变换)处理，计算出K点数据的频谱Xi(k,m)，这里K的取值根据实际硬件资源情况可选择4096/8192/16384三档。
其中表示FFT运算的批次数，j表示
步骤(3)对Xi(k,m)频谱信号进行频谱搬移。将Xi(k,m)中[fc-△f,fc+△f]频谱搬移到基带两侧，其余部分的频谱数值置0，从而完成基带解调和复包络化处理。搬移后的频谱信号序列为：
步骤(4)将搬移后的基带频谱信号Yi(k,m)先进行IFFT(快速傅里叶逆变换)处理，还原为时域复信号序列：
然后消除50％的重叠部分，得到最终的时域复包络序列，即
步骤(5)对四路时域复包络序列按下式进行复信号波束合成，得到两路合成后的信号：
上式中，Imag()/Real()分别表示取虚部/取实部运算。
步骤(6)利用波束合成信号进行水底跟踪判断，以确定底回波前沿数据序号。判定方法为：计算合成信号中实部序列zR1(n)/zR2(n)中M点的实时平均功率，M取值由水底大致深度决定，范围是120～9600。
计算从回波接收开始至当前时刻，合成信号中实部序列zR1(n)和zR2(n)的实时平均功率：
如果PM1(L)>α×P1(L)且PM2(L)>α×P2(L)，则记录此时的L值作为底回波前沿数据序号，水底跟踪成功，进入步骤(7)；否则重复步骤(6)。这里α为功率检测系数，取值范围为1～20。
步骤(7)分别计算两路合成信号的复相关值：
步骤(8)根据复相关结果计算相角：
φ1,φ2∈[-π,π]，分别为两个复相关值的相角。
步骤(9)根据相角计算出所需测量的速度值：
v数值的正负代表所测速度的方向与参考方向一致或相反。
本发明采用基于频域的频谱搬移与正交化处理方式，结合复包络信号波束合成方法，在保证性能的前提下简化了信号处理流程，大幅度降低了处理复杂度和硬件资源的规模。且可采用提高工作频率的方式，对多路信号的FFT/IFFT运算处理模块加以复用，只需一路信号处理资源按时分方式即可完成多路信号的处理，这将进一步减少硬件资源的开销。
具体实施方式
以下结合实施例对本发明作进一步说明。
本发明方法包括如下步骤：
步骤(1)对相控换能器阵接收到的四路回波信号进行采样，形成四路数字信号。采样频率fs，采样后的信号为xi(n)，其中i＝1,2,3,4表示信号对应的路数，n＝1,2,3,…,N表示采样点序号，N表示采样后数字序列的长度。
步骤(2)对每路信号各自按样本重叠50％的方式做FFT(快速傅里叶变换)处理，计算出K点数据的频谱Xi(k,m)，这里K的取值根据实际硬件资源情况可选择4096/8192/16384三档。
其中表示FFT运算的批次数，j表示
步骤(3)对Xi(k,m)频谱信号进行频谱搬移。将Xi(k,m)中[fc-△f,fc+△f]频谱搬移到基带两侧，其余部分的频谱数值置0，从而完成基带解调和复包络化处理。搬移后的频谱信号序列为：
步骤(4)将搬移后的基带频谱信号Yi(k,m)先进行IFFT(快速傅里叶逆变换)处理，还原为时域复信号序列：
然后消除50％的重叠部分，得到最终的时域复包络序列，即
步骤(5)对四路时域复包络序列按下式进行复信号波束合成，得到两路合成后的信号：
上式中，Imag()/Real()分别表示取虚部/取实部运算。
步骤(6)利用波束合成信号进行水底跟踪判断，以确定底回波前沿数据序号。判定方法为：计算合成信号中实部序列zR1(n)/zR2(n)中M点的实时平均功率，M取值由水底大致深度决定，范围是120～9600。
计算从回波接收开始至当前时刻，合成信号中实部序列zR1(n)和zR2(n)的实时平均功率：
如果PM1(L)>α×P1(L)且PM2(L)>α×P2(L)，则记录此时的L值作为底回波前沿数据序号，水底跟踪成功，进入步骤(7)；否则重复步骤(6)。这里α为功率检测系数，取值范围为1～20。
步骤(7)分别计算两路合成信号的复相关值：
步骤(8)根据复相关结果计算相角：
φ1,φ2∈[-π,π]，分别为两个复相关值的相角。
步骤(9)根据相角计算出所需测量的速度值：
v数值的正负代表所测速度的方向与参考方向一致或相反。

Claims (1)

1.宽带相控声学多普勒测速系统的信号处理方法，该方法中设系统宽带相控声学多普勒测速系统中进行相位编码调制的伪随机序列周期为Tm，相控换能器发射的声脉冲信号的波束张角为θ，波长为λ，声脉冲信号的载波频率为fc，系统处理的最大多普勒频移绝对值为△f，其特征在于该方法的具体步骤是：
步骤(1)对相控换能器阵接收到的四路回波信号进行采样，形成四路数字信号；采样频率fs，采样后的信号为xi(n)，其中i＝1,2,3,4表示信号对应的路数，n＝1,2,3,…,N表示采样点序号，N表示采样后数字序列的长度；
步骤(2)对每路信号各自按样本重叠50％的方式做快速傅里叶变换处理，计算出K点数据的频谱Xi(k,m)；
其中表示快速傅里叶变换运算的批次数，j表示
步骤(3)对Xi(k,m)频谱信号进行频谱搬移；将Xi(k,m)中[fc-△f,fc+△f]频谱搬移到基带两侧，其余部分的频谱数值置0；搬移后的频谱信号序列为：
步骤(4)将搬移后的基带频谱信号Yi(k,m)先进行快速傅里叶逆变换处理，还原为时域复信号序列：
然后消除50％的重叠部分，得到最终的时域复包络序列，即
步骤(5)对四路时域复包络序列按下式进行复信号波束合成，得到两路合成后的信号：
上式中，Imag()/Real()分别表示取虚部/取实部运算；
步骤(6)进行水底跟踪判断，确定底回波前沿数据序号；具体判定为：计算合成信号中实部序列zR1(n)/zR2(n)中M点的实时平均功率，M取值范围：120～9600；
计算从回波接收开始至当前时刻，合成信号中实部序列zR1(n)和zR2(n)的实时平均功率：
如果PM1(L)>α×P1(L)且PM2(L)>α×P2(L)，则记录此时的L值作为底回波前沿数据序号，进入步骤(7)；否则重复步骤(6)；这里α为功率检测系数，取值范围为1～20；
步骤(7)分别计算两路合成信号的复相关值：
步骤(8)根据复相关结果计算相角：
φ1,φ2∈[-π,π]，分别为两个复相关值的相角；
步骤(9)根据相角计算出所需测量的速度值：
v数值的正负代表所测速度的方向与参考方向一致或相反。
CN201610338300.XA 2016-05-19 2016-05-19 宽带相控声学多普勒测速系统信号处理方法 Active CN106019289B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610338300.XA CN106019289B (zh) 2016-05-19 2016-05-19 宽带相控声学多普勒测速系统信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610338300.XA CN106019289B (zh) 2016-05-19 2016-05-19 宽带相控声学多普勒测速系统信号处理方法

Publications (2)

Publication Number Publication Date
CN106019289A CN106019289A (zh) 2016-10-12
CN106019289B true CN106019289B (zh) 2018-06-26

Family

ID=57095596

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610338300.XA Active CN106019289B (zh) 2016-05-19 2016-05-19 宽带相控声学多普勒测速系统信号处理方法

Country Status (1)

Country Link
CN (1) CN106019289B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101296014A (zh) * 2008-06-12 2008-10-29 浙江大学 一种宽带数字波束形成方法
CN102901839A (zh) * 2012-11-05 2013-01-30 杭州电子科技大学 一种低复杂度的声学多普勒测流系统速度估算方法
CN102928619A (zh) * 2012-11-05 2013-02-13 杭州电子科技大学 一种宽带束控声学多普勒测流系统信号处理方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101296014A (zh) * 2008-06-12 2008-10-29 浙江大学 一种宽带数字波束形成方法
CN102901839A (zh) * 2012-11-05 2013-01-30 杭州电子科技大学 一种低复杂度的声学多普勒测流系统速度估算方法
CN102928619A (zh) * 2012-11-05 2013-02-13 杭州电子科技大学 一种宽带束控声学多普勒测流系统信号处理方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Asymptotic Performance of a Pulse-to-Pulse Incoherent Doppler Sonar in an Oceanic Environment;D.SCOTT HANSEN;《IEEE JOURNAL OF OCEANIC ENGINEERING》;19850430;第OE-10卷(第2期);第144-158页 *
一种基于FFT的超宽带数字正交变换算法及实现;王旭东 等;《仪器仪表学报》;20080831;第29卷(第8期);第1714-1718页 *
基于FPGA的多普勒计程仪硬件平台设计;马修准 等;《杭州电子科技大学学报》;20140131;第34卷(第1期);第75-78页 *
相控阵声学多普勒计程仪的信号处理模块设计;曹永红;《声学与电子工程》;20091231(第93期);第44-46页 *

Also Published As

Publication number Publication date
CN106019289A (zh) 2016-10-12

Similar Documents

Publication Publication Date Title
CN102901839B (zh) 一种低复杂度的声学多普勒测流系统速度估算方法
CN112083404B (zh) 一种基于多途特征匹配的单矢量水听器声源深度估计方法
CN102981146B (zh) 一种单矢量水听器被动定位方法
CN103048642B (zh) 基于频域最小二乘法的水声脉冲信号匹配场定位方法
CN107179535A (zh) 一种基于畸变拖曳阵的保真增强波束形成的方法
CN109188387B (zh) 基于插值补偿的分布式相参雷达目标参数估计方法
CN107966688A (zh) 基于相位干涉技术的宽带雷达目标速度解模糊方法
CN104777453A (zh) 舰船线谱噪声源定位的波束域时频分析方法
CN103604944B (zh) 一种基于单站船载高频地波雷达的表面流测量方法
CN111580048B (zh) 一种利用单矢量水听器的宽带声源深度估计方法
CN113109817B (zh) 一种矢量水听器布放深度估计方法
CN104678384B (zh) 一种波束域的声压差互相关谱分析水下目标速度估计方法
CN103323815A (zh) 一种基于等效声速的水下声学定位方法
CN101813772A (zh) 一种快速宽带频域扩展拖曳阵波束形成方法
CN110231616B (zh) 基于北斗卫星辐射源的海面运动目标检测与定位方法
CN109814094A (zh) 一种多目标水声定位时延估计算法
CN108562891A (zh) 一种深海低信噪比条件下声源深度自主实时跟踪方法
CN108107436A (zh) 一种基于可靠声路径的水下目标主动分类与定位方法
CN104808207A (zh) 一种混沌水声定位方法
CN103809170A (zh) 一种多普勒频偏修正方法及系统
CN103845080A (zh) 一种基于线性调频编码的超声脐带血流测量系统及方法
CN102928619B (zh) 一种宽带束控声学多普勒测流系统信号处理方法
CN106019289B (zh) 宽带相控声学多普勒测速系统信号处理方法
CN109541579A (zh) 基于Bezier模型的霍夫变换的多普勒穿墙雷达定位方法
CN103995262B (zh) 用于起伏界面的mimo稀疏阵列的超声测量方法及系统

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
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20161012

Assignee: ZHEJIANG REALLIN ELECTRON Co.,Ltd.

Assignor: HANGZHOU DIANZI University

Contract record no.: X2021330000745

Denomination of invention: Signal processing method of broadband phased acoustic Doppler velocimetry system

Granted publication date: 20180626

License type: Common License

Record date: 20211111