CN103913739A - 水声通信中用于载体相对运动稳健测速的方法 - Google Patents
水声通信中用于载体相对运动稳健测速的方法 Download PDFInfo
- Publication number
- CN103913739A CN103913739A CN201410087900.4A CN201410087900A CN103913739A CN 103913739 A CN103913739 A CN 103913739A CN 201410087900 A CN201410087900 A CN 201410087900A CN 103913739 A CN103913739 A CN 103913739A
- Authority
- CN
- China
- Prior art keywords
- signal
- infin
- tau
- doppler
- signals
- 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.)
- Pending
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/586—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
-
- 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/523—Details of pulse systems
- G01S7/526—Receivers
- G01S7/527—Extracting wanted echo signals
- G01S7/5273—Extracting wanted echo signals using digital techniques
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
本发明涉及一种水声通信中用于载体相对运动稳健测速的方法。本发明包括:将可控增益滤波电路的通频带设置成增益;信号处理电路通过A/D变换器采集信号,将采集后的信号进行数字滤波,用匹配滤波器对接收的信号进行滤波,对滤波后的信号进行取包络;对接收的信号进行相关处理。本方法解决了传统测速方法所采用的单频信号在衰落信道中信号易于衰落,采用LFM脉冲对信号,利用线性调频信号宽带的特点,即使频带内的部分频点受到影响,仍可对接收信号进行连续速度估计。
Description
技术领域
本发明涉及一种水声通信中用于载体相对运动稳健测速的方法。
背景技术
在水声通信中,由于通信节点的运动导致的多普勒效应使得通信信号被压缩和扩展,一般采用对接收信号的功率谱估计来得到多普勒频偏值,如FFT法,复相关法等,且一般采用CW脉冲信号作为发射信号。但这些方法使用的窄带信号在水平移动的通信场合,由于信道变化较快且信道频谱呈现梳妆结构,使得接收信号可能在某段时刻信号幅度很小,因而不利于多普勒速度的连续估计。有关多普勒测速的本章较多,但方法大多采用CW脉冲信号,未见有关LFM脉冲对的测速的资料。可参见耿志辉、王玉林“实正弦信号频率估计的高精度综合算法”信号与信息处理,2008年28卷第3期;徐世军、许伟杰“复相关法测频的硬件实现”声学技术,2009年第28卷第2期。
发明内容
本发明的目的是为了克服现有的测速方法不利于水声通信衰落信道中连续测速,提出了一种水声通信中用于载体相对运动稳健测速的方法。
本发明的目的是这样实现的:
(1)将可控增益滤波电路的通频带设置成增益;
(2)信号处理电路通过A/D变换器采集信号,将采集后的信号进行数字滤波,用匹配滤波器对接收的信号进行滤波,对滤波后的信号进行取包络,相关峰的包络会有一个主峰和两个侧峰,且在没有多普勒时,主峰和侧峰的间距为Tc=τ,当接收到信号存在多普勒时,主峰和侧峰的间距变为Td,根据多普勒原理可以得到脉冲压缩系数a为
fc为发射信号的中心频率,得到没有模糊的粗的多普勒频移fd=(1-a)fc;
(3)对接收的信号进行相关处理,接收到的信号为W(t),由待测量信号x(t)和加性高斯白噪声n(t)组成:
W(t)=x(t)+n(t),
待测信号x(t)与噪声n(t)之间不相关,则有
E[x*(t)n(t+τ)]=0,
E[n*(t)x(t+τ)]=0,
接收信号的自相关函数为
R(τ)=E[W*(t)W(t+τ)]=E[x*(t)x(t+τ)]+E[n*(t)n(t+τ)]+E[x*(t)n(t+τ)]+E[n*(t)x(t+τ)]=Rx(τ)+Rn(τ)
Rx(τ)为待测信号的自相关函数,Rn(τ)为加性高斯白噪声的自相关函数,
R(τ)=Rx(τ)
信号的自相关函数与功率谱互为傅里叶变换,
Sx(w)为信号x(t)的功率谱,把信号的自相关函数用极坐标表示
Ax(τ)是偶函数,是奇函数,
Ax(τ)是偶函数,则有A'x(0)=0,
令τ=0则有
令τ=0则有
多普勒频移可以表示为功率谱Sx(w)的一阶矩:
由w=2πf,得
测得的频偏为fp,如果fd和fp都在内,没有相位模糊,则真实的多普勒频率为fr=fp;若fd超出出现相位模糊,脉冲对自相关函数求得得相位为,真实的相位为k=±1,±2,...则有
k=(fd-fp)τ,得到无模糊地精确值估计出多普勒速度c为声速,数字带通滤波、匹配滤波处理是线性运算,得到的相关峰峰值点的幅值和接受信号的能量成正比,利用幅值来估计接收信号的声级。
本发明的有益效果在于:采用模拟电路和数字电路混合构成闭环电路,利用线性调频信号时间带宽积远大于1的特点,完成接收信号的检测、频率粗估计、频率细估计、接收增益控制。本发明的优点是能在衰落水声信道中完成多普勒速度的连续跟踪估计。提出了一种水声通信中用于载体相对运动稳健测速的方法。本方法解决了传统测速方法所采用的单频信号在衰落信道中信号易于衰落,采用LFM脉冲对信号,利用线性调频信号宽带的特点,即使频带内的部分频点受到影响,仍可对接收信号进行连续速度估计。
附图说明
图1为水声通信中用于载体相对运动稳健测速的组成框图。
具体实施方式
下面结合附图对本发明做进一步描述:
即使频带内的部分频点受到影响,仍可对接收信号进行连续频偏估计。
本发明采用模拟电路和数字电路混合构成闭环电路,利用线性调频信号时间带宽积远大于1的特点,完成接收信号的检测、频率粗估计、频率细估计、接收增益控制,将采集后的信号进行数字滤波,然后用匹配滤波器对接收的信号进行滤波,对滤波后的信号进行取包络,根据检测门限判断信号的有无,然后利用包络相关峰估计粗的频偏,对脉冲对信号用相互正交的两路正余弦信号做乘法运算,得到基带复信号,接着进行相关运算,根据相关运算得到复结果的相位计算得到带有模糊地精确地频偏值,并用之前粗的频偏来修正,得到精确地测速结果。最后根据相关峰峰值来估计接收信号的声级,调节可控增益滤波电路[1]。
本方法解决了传统测速方法所采用的单频信号在衰落信道中信号易于衰落,采用LFM脉冲对信号,利用线性调频信号宽带的特点,即使频带内的部分频点受到影响,仍可对接收信号进行连续速度估计。本发明的构成主要由可控增益滤波电路[1],A/D变换器[2]和信号处理电路[3]组成,如图1所示。其中可控增益滤波电路[1]、A/D变换器[2]和信号处理电路[3]依次电信号连接,信号处理电路[3]又与可控增益滤波电路[1]构成闭环电路。
本发明的测速方法是:首先将可控增益滤波电路[1]的通频带设置成合适的增益;然后信号处理电路[3]通过A/D变换器[2]采集信号,并将采集后的信号进行数字滤波,用匹配滤波器对接收的信号进行滤波,对滤波后的信号进行取包络,相关峰的包络会有一个主峰和两个侧峰,且在没有多普勒时,主峰和侧峰的间距为Tc=τ,当接收到信号存在多普勒时,由于信号发生压缩或扩展,主峰和侧峰的间距变为Td,根据多普勒原理可以得到脉冲压缩系数a为
fc为发射信号的中心频率,根据(1)式,只要测得时延差Td,可以得到一个没有模糊的粗的多普勒频移fd=(1-a)fc。再对接收的LFM脉冲对信号进行相关处理,假设接收到得信号为W(t),它由待测量信号x(t)和加性高斯白噪声n(t)组成,可表示为:
W(t)=x(t)+n(t) (2)
待测信号x(t)与噪声n(t)之间不相关,则有
E[x*(t)n(t+τ)]=0 (3)
E[n*(t)x(t+τ)]=0 (4)
接收信号的自相关函数可以用下式表示
R(τ)=E[W*(t)W(t+τ)]=E[x*(t)x(t+τ)]+E[n*(t)n(t+τ)]+E[x*(t)n(t+τ)]+E[n*(t)x(t+τ)]=Rx(τ)+Rn(τ) (5)
Rx(τ)为待测信号的自相关函数,Rn(τ)为加性高斯白噪声的自相关函数。因白噪声自相关函数的性质
将式(6)带入式(5)可以得到
R(τ)=Rx(τ) (7)
式(7)说明可以由接收信号的自相关函数来估计得到待测信号的自相关函数,又由维纳-辛钦定理,信号的自相关函数与功率谱互为傅里叶变换,即
Sx(w)为信号x(t)的功率谱。若把信号的自相关函数用极坐标表示,则有
根据自相关函数的性质可知,Ax(τ)是偶函数,是奇函数。对Rx(τ)求导可得
又由于Ax(τ)是偶函数,则其导数为奇函数,则有A'x(0)=0,因此可求得
对式(8)求导可得
令τ=0则有
式(8)中令τ=0则有
多普勒频移可以表示为功率谱Sx(w)的一阶矩:
如果τ≠0且充分小时,根据导数的定义有
代入式(17)可得又由w=2πf,可得
假设采用脉冲对法测得的频偏为fp,如果fd和fp都在内,说明没有相位模糊,则真实的多普勒频率为fr=fp。若fd超出这时出现相位模糊,设脉冲对自相关函数求得得相位为真实的相位为k=±1,±2,...则有
由式(20)、(21)可以求得k=(fd-fp)τ,这样就可以得到无模糊地精确值由测得的频偏即可估计出多普勒速度c为声速。由于数字带通滤波、匹配滤波等处理是线性运算,得到的相关峰峰值点的幅值和接受信号的能量成正比,可利用幅值来估计接受信号的声级,从而动态调整可控增益滤波电路[1]使之进入A/D变换器[2]的信号幅值保持在合适的范围。
本发明水声通信中用于载体相对运动稳健测速的方法所采用的电路器件均为市场上可以购买到得常用器件;其中可控增益滤波电路[1]采用VC810;A/D变换器选用双极性A/D变换器,如AD7663;信号处理电路[3]是由集成数字信号处理芯片(DSP芯片)电路构成,如OMAPL138。其中数字滤波算法采用有限冲击响应滤波算法,匹配滤波器可采用快速傅里叶算法计算。
本发明的测速方法是:将采集后的信号进行数字滤波,然后用匹配滤波器对接收的信号进行滤波,对滤波后的信号进行取包络,根据检测门限判断信号的有无,然后利用包络相关峰估计粗的频偏,然后利用对脉冲对信号用相互正交的两路信号做乘法运算,得到基带复信号,接着进行相关运算,根据相关运算得到的复结果的相位计算得到带有模糊地精确地频偏值,并用之前粗的频偏来修正,得到精确地测频结果。最后根据相关峰峰值来估计接收信号的声级,调节可控增益滤波电路[1]。本发明可以应用于水声通信中载体相对运动稳健测速。
Claims (1)
1.水声通信中用于载体相对运动稳健测速的方法,其特征在于:
(1)将可控增益滤波电路的通频带设置成增益;
(2)信号处理电路通过A/D变换器采集信号,将采集后的信号进行数字滤波,用匹配滤波器对接收的信号进行滤波,对滤波后的信号进行取包络,相关峰的包络会有一个主峰和两个侧峰,且在没有多普勒时,主峰和侧峰的间距为Tc=τ,当接收到信号存在多普勒时,主峰和侧峰的间距变为Td,根据多普勒原理可以得到脉冲压缩系数a为
fc为发射信号的中心频率,得到没有模糊的粗的多普勒频移fd=(1-a)fc;
(3)对接收的信号进行相关处理,接收到的信号为W(t),由待测量信号x(t)和加性高斯白噪声n(t)组成:
W(t)=x(t)+n(t),
待测信号x(t)与噪声n(t)之间不相关,则有
E[x*(t)n(t+τ)]=0,
E[nw(t)x(t+τ)]=0,
接收信号的自相关函数为
R(τ)=E[W*(t)W(t+τ)]=E[x*(t)x(t+τ)]+E[n*(t)n(t+τ)]+E[x*(t)n(t+τ)]+E[n*(t)x(t+τ)]=Rx(τ)+Rn(τ)
Rx(τ)为待测信号的自相关函数,Rn(τ)为加性高斯白噪声的自相关函数,
R(τ)=Rx(τ)
信号的自相关函数与功率谱互为傅里叶变换,
Sx(w)为信号x(t)的功率谱,把信号的自相关函数用极坐标表示
Ax(τ)是偶函数,是奇函数,
Ax(τ)是偶函数,则有A'x(0)=0,
令τ=0则有
令τ=0则有
多普勒频移可以表示为功率谱Sx(w)的一阶矩:
由w=2πf,得
测得的频偏为fp,如果fd和fp都在内,没有相位模糊,则真实的多普勒频率为fr=fp;若fd超出出现相位模糊,脉冲对自相关函数求得得相位为,真实的相位为k=±1,±2,...则有
k=(fd-fp)τ,得到无模糊地精确值估计出多普勒速度c为声速,数字带通滤波、匹配滤波处理是线性运算,得到的相关峰峰值点的幅值和接受信号的能量成正比,利用幅值来估计接收信号的声级。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410087900.4A CN103913739A (zh) | 2014-03-12 | 2014-03-12 | 水声通信中用于载体相对运动稳健测速的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410087900.4A CN103913739A (zh) | 2014-03-12 | 2014-03-12 | 水声通信中用于载体相对运动稳健测速的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103913739A true CN103913739A (zh) | 2014-07-09 |
Family
ID=51039565
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410087900.4A Pending CN103913739A (zh) | 2014-03-12 | 2014-03-12 | 水声通信中用于载体相对运动稳健测速的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103913739A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113627398A (zh) * | 2021-10-11 | 2021-11-09 | 西安瀚博电子科技有限公司 | 一种基于自适应重构滤波的信号特征检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005208068A (ja) * | 2005-02-21 | 2005-08-04 | Keio Gijuku | 超音波流速分布計及び流量計、超音波流速分布及び流量測定方法並びに超音波流速分布及び流量測定処理プログラム |
CN101046488A (zh) * | 2007-03-29 | 2007-10-03 | 哈尔滨工程大学 | Cw脉冲的窄带跟踪滤波器及其跟踪滤波方法 |
-
2014
- 2014-03-12 CN CN201410087900.4A patent/CN103913739A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005208068A (ja) * | 2005-02-21 | 2005-08-04 | Keio Gijuku | 超音波流速分布計及び流量計、超音波流速分布及び流量測定方法並びに超音波流速分布及び流量測定処理プログラム |
CN101046488A (zh) * | 2007-03-29 | 2007-10-03 | 哈尔滨工程大学 | Cw脉冲的窄带跟踪滤波器及其跟踪滤波方法 |
Non-Patent Citations (3)
Title |
---|
张殿伦等: "声相关测速技术仿真与试验研究", 《哈尔滨工程大学学报》 * |
李想: "水下高速运动目标轨迹测量技术研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
盛元平等: "多普勒计程仪检测装置的设计与实现", 《通信电源技术》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113627398A (zh) * | 2021-10-11 | 2021-11-09 | 西安瀚博电子科技有限公司 | 一种基于自适应重构滤波的信号特征检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101813673B (zh) | 检测微量二元混合气体浓度的声学信号处理装置及方法 | |
US11255708B2 (en) | Extended range ADC flow meter | |
CN101957443B (zh) | 声源定位方法 | |
CN109188017B (zh) | 一种基于超声波传感器阵列的风向风速检测装置及检测方法 | |
CN103845079B (zh) | 一种基于盲分离的多普勒胎心音瞬时心率的检测方法 | |
CN102901839B (zh) | 一种低复杂度的声学多普勒测流系统速度估算方法 | |
CN104360251B (zh) | 一种变压器局部放电的超声波信号时延估计方法 | |
CN104168232B (zh) | 一种水声信道中多径时延与多普勒频移的测定方法 | |
CN104777328A (zh) | 一种可听声测量气体流速的系统及测量方法 | |
CN104780127B (zh) | 一种时延‑多普勒r‑l解卷多路径信道估计方法 | |
CN103618575B (zh) | 一种调频水声通信系统的自动实时帧同步方法 | |
CN107800659A (zh) | Alpha稳定分布噪声下LFM信号调制参数估计方法 | |
CN106330342B (zh) | 一种低计算复杂度的水声通信多普勒因子估计方法 | |
CN105745556B (zh) | 用于借助声传感器测量相对速度的方法和设备 | |
CN101576618A (zh) | 基于小波变换的水声定位测量方法和测量系统 | |
CN104901718A (zh) | 基于直接序列扩频信号载波频率测量的多普勒估计方法 | |
CN103913739A (zh) | 水声通信中用于载体相对运动稳健测速的方法 | |
CN104049247B (zh) | 一种基于频率方差加权的时延差估计方法 | |
CN102087313A (zh) | 一种卫星搜救信号的频率估计方法 | |
CN105380681A (zh) | 基于可变频解调的超声多普勒音频信号获取系统及方法 | |
KR101839236B1 (ko) | 천부지층 탐지 장치 및 방법 | |
CN111899749B (zh) | 一种监听混凝土泵送管道作业声音的降噪方法 | |
Duarte et al. | Ultrasonic time of flight estimation for wind speed measurement based on time-frequency domain using STFT | |
CN104199013A (zh) | 一种能在有限水域内降低测试频率的方法 | |
CN111342949B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140709 |
|
WD01 | Invention patent application deemed withdrawn after publication |