CN109510787B - 水声信道下线性调频信号参数估计方法及装置 - Google Patents

水声信道下线性调频信号参数估计方法及装置 Download PDF

Info

Publication number
CN109510787B
CN109510787B CN201811196283.6A CN201811196283A CN109510787B CN 109510787 B CN109510787 B CN 109510787B CN 201811196283 A CN201811196283 A CN 201811196283A CN 109510787 B CN109510787 B CN 109510787B
Authority
CN
China
Prior art keywords
frequency
signal
chirp signal
component
underwater acoustic
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
CN201811196283.6A
Other languages
English (en)
Other versions
CN109510787A (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.)
Information Engineering University of PLA Strategic Support Force
Original Assignee
Information Engineering University of PLA Strategic Support Force
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 Information Engineering University of PLA Strategic Support Force filed Critical Information Engineering University of PLA Strategic Support Force
Priority to CN201811196283.6A priority Critical patent/CN109510787B/zh
Publication of CN109510787A publication Critical patent/CN109510787A/zh
Application granted granted Critical
Publication of CN109510787B publication Critical patent/CN109510787B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B11/00Transmission systems employing sonic, ultrasonic or infrasonic waves
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B13/00Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
    • H04B13/02Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明属于水声通信技术领域,特别涉及一种水声信道下线性调频信号参数估计方法及装置,该方法包含:构建水声通道信号模型;接收通过水声通道模型处理后的N分量Chirp信号,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频率;对N分量Chirp信号依次进行调频补偿,并通过快速傅里叶变换对起始频率估计并输出。本发明将FRFT域起始频率估计通过调频率分量处理后旋转到FFT域处理,获得低信噪比下均方误差更小的起始频率估计值;并通过对水声信道建模,在水声多径信道的背景下获得多分量Chirp信号参数估计,提升起始频率估计性能,保证信号重构准确率,提高水声通信质量,保证通信的稳健,性能稳定、运行高效,具有较强的实际应用价值和发展前景。

Description

水声信道下线性调频信号参数估计方法及装置
技术领域
本发明属于水声通信技术领域,特别涉及一种水声信道下线性调频信号参数估计方法及装置。
背景技术
线性调频信号(Chirp或称LFM)作为一种成熟的非平稳信号,被广泛的应用于雷达、声呐和水声通信中。起始频率和调频率是Chirp信号的关键参数,对于两者的估计取得了广泛的研究。
在传统无线信道的背景下,针对Chirp信号的参数估计方法主要有基于最大似然(Maximum Likelihood,ML)、短时Fourier变换、小波变换、Wigner-Hough变换、Wigner-Ville分布、分数阶傅里叶变换(Fractional Fourier transform,FRFT)等方法,其中最大似然估计方法估计性能好,但其运算复杂度过高,基于Wigner-Ville分布的参数估计方法在估计多分量线性调频信号时存在交叉项问题,小波变换法对多分量线性调频信号估计性能好,但也存在计算量较大的问题,FRFT是以一组正交的Chirp函数作为变换基底,非常适合于对线性调频信号的参数估计。以上方法在无线信道的背景下对参数估计的取得了较好效果,但不适用于水声信道下线性调频信号的参数估计。水声信道是时延和频移双扩散的复杂信道,其多途效应十分严重,水声通信中的时延可达几十毫秒甚至几百毫秒,水声信道时延对Chirp信号的调频率影响较小,但是对起始频率会造成严重影响,导致在无线信道背景下经典的参数估计算法无法直接应用于水声信道,为了获得精确地参数估计性能,迫切需要寻求一种适应于水声信道下针对Chirp信号的参数估计方法。
发明内容
为此,本发明提供一种水声信道下线性调频信号参数估计方法及装置,实现低吸噪比、低占空比、欠采样等恶劣影响下以及多分量Chirp信号的参数估计,便于信号重构,提高系统估计精度和处理效率。
按照本发明所提供的设计方案,一种水声信道下线性调频信号参数估计方法,包含如下内容:
构建用于输入并处理N分量线性调频Chirp信号的水声通道信号模型;
接收通过水声通道模型处理后的N分量Chirp信号,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频率;
对N分量Chirp信号依次进行调频补偿,并通过快速傅里叶变换对起始频率估计并输出。
上述的,构建水声通道模型,包含如下内容:依据单分量Chirp信号,及消除多普勒效应水声多途信道的时变冲激响应,获取N分量Chirp信号经过多途时延的混合叠加信号模型。
优选的,通过信号幅度、起始频率和调频率来表示单分量Chirp信号,通过信道路径的幅度和时延来表示时变冲激响应。
上述的,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频,包含如下内容:设定最佳阶次范围,通过峰值搜索得到不同阶次的峰值,获取N分量Chirp信号的调频率。
优选的,针对接收到的Chirp信号,获取其分数阶傅里叶变换,并通过单分量Chirp信号,获取Chirp信号在分数阶傅里叶变换域内变成一个脉冲的旋转角度,通过该旋转角度得到最佳阶次,依据最佳阶次获取对应分量的调频率。
优选的,设定滑动步长进行不同阶次的峰值搜索得到最佳阶次,依据得到的最佳阶次获取多分量Chirp信号的调频率。
上述的,获取N分量Chirp信号调频率后,对N分量Chirp信号依次进行调频补偿,包含如下内容:首先,将接收到的Chirp信号转换为单频信号;然后,针对每个分量Chirp信号的单频信号进行快速傅里叶变换,得到该分量Chirp信号对应的频率估计值序列,并联立信道时延和起始频率,得到起始频率估计值。
优选的,转换单频信号过程中,通过如下公式实现:
Figure BDA0001828831930000031
其中,N分量Chirp信号对应的调频率为[k1,k2,……kN],R(t)为接收信号。
一种水声信道下线性调频信号参数估计装置,包含模型构建模块、数据获取模块和数据输出模块,其中,
模型构建模块,用于构建用于输入并处理N分量线性调频Chirp信号的水声通道信号模型;
数据获取模块,用于接收通过水声通道模型处理后的N分量Chirp信号,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频率;
数据输出模块,用于对N分量Chirp信号依次进行调频补偿,并通过快速傅里叶变换对起始频率估计并输出。
上述的装置中,数据输出模块包含转换子模块、估值子模块和输出子模块,其中,
转换子模块,用于将接收到的Chirp信号转换为单频信号;
估值子模块,用于针对每个分量Chirp信号的单频信号进行快速傅里叶变换,得到该分量Chirp信号对应的频率估计值序列,并联立信道时延和起始频率,得到起始频率估计值;
输出子模块,用于将得到的起始频率估计值作为线性调频信号参数进行输出。
本发明的有益效果:
本发明针对水声信道下Chirp信号参数估计问题,利用经过水声信道后Chirp信号的调频率不变的特性,通过FRFT变换进行阶次搜索得到多分量Chirp信号的调频率;对接收到的信号依次进行调频率补偿后,借助FFT变换对受不同信道时延影响的起始频率进行区分,得到起始频率估计值;发挥FRFT对多分量Chirp信号处理的优势,获得精确的调频率估计,将FRFT域的起始频率估计通过调频率分量处理后旋转到FFT域处理,从而获得低信噪比下均方误差更小的的起始频率的估计值;并通过对水声信道建模,对Chirp信号经过水声信道后的参数变化进行讨论,在水声多径信道的背景下获得多分量Chirp信号参数估计,提升起始频率的估计性能,保证通信的稳健,性能稳定、运行高效,具有较强的实际应用价值和发展前景。
附图说明:
图1为实施例中线性调频信号参数估计方法流程示意图;
图2为实施例中信号参数估计框架;
图3为实施例中线性调频信号参数估计装置示意图;
图4为实施例中数据输出模块示意图;
图5为实施例中归一化信道冲激响应示意图;
图6为实施例中不同阶数FRFT最大幅值示意图之一;
图7为实施例中通过估计调频率补偿后频域示意图之一;
图8为实施例中不同阶数FRFT最大幅值示意图之二;
图9为实施例中通过估计调频率补偿后频域示意图之二;
图10为实施例中起始频率的均方误差示意图。
具体实施方式:
下面结合附图和技术方案对本发明作进一步详细的说明,并通过优选的实施例详细说明本发明的实施方式,但本发明的实施方式并不限于此。
目前,针对Chirp信号的参数估计方法在无线通信背景下已取得丰硕的成果,针对低信噪比下、低占空比下、欠采样等恶劣影响下以及多分量Chirp信号的参数估计问题进行了丰富的研究。但是,在对Chirp信号参数估计时现有估计方法中大多将Chirp信号默认为雷达信号,假设的背景为无线信道,在实际应用中,Chirp信号作为具有大多普勒容限的信号,在水声通信中广泛的被应用于同步信号以及载波信号中,研究水声信道下Chirp信号的参数估计问题有着重要的意义。为此,本申请实施例中,参见图1所示,提供一种水声信道下线性调频信号参数估计方法,包含如下内容:
构建用于输入并处理N分量线性调频Chirp信号的水声通道信号模型;
接收通过水声通道模型处理后的N分量Chirp信号,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频率;
对N分量Chirp信号依次进行调频补偿,并通过快速傅里叶变换对起始频率估计并输出。
针对水声信道下Chirp信号参数估计问题,首先利用经过水声信道后Chirp信号的调频率不变的特性,通过FRFT变换进行阶次搜索得到多分量Chirp信号的调频率;对接收到的信号依次进行调频率补偿后,借助FFT变换对受不同信道时延影响的起始频率进行区分,得到起始频率估计值,提升对起始频率的估计性能,可靠性高,具有较强的实际应用价值。
构建水声通道模型中,本发明再一个实施例中,依据单分量Chirp信号,及消除多普勒效应水声多途信道的时变冲激响应,获取N分量Chirp信号经过多途时延的混合叠加信号模型。并进一步地,通过信号幅度、起始频率和调频率来表示单分量Chirp信号,通过信道路径的幅度和时延来表示时变冲激响应。
考虑单分量Chirp信号s(t)和消除多普勒效应的水声多途信道的时变冲激响应h(t,τ)表达式为:
Figure BDA0001828831930000051
Figure BDA0001828831930000052
式中,A,f0和k分别为Chirp信号的幅度,起始频率和调频率,Al和τl分别为信道第l条路径的幅度和时延。参数估计的目标就是对f0和k进行估计,假设发送信号为s(t),则Chirp信号通过水声信道后信号表示为:
Figure BDA0001828831930000061
其中n(t)为高斯白噪声,将式(1)代入式(3),得到:
Figure BDA0001828831930000062
其中:
Figure BDA0001828831930000063
显然,r(t)仍是线性调频信号,只是其参数发生了变化,k'和f0'为经过信道后的信号起始频率和调频率,
Figure BDA0001828831930000064
为相位。对于N分量Chirp信号,经过水声信道后接收信号为:
Figure BDA0001828831930000065
可见接收到的信号为N分量Chirp信号经过L条多途时延后的混合叠加,Chirp信号经过多途信道后参数估计问题。借助FRFT实现对多分量信号的调频率估计,在得到调频率的估计后,对接收信号进行调频率补偿将对应调频率的线性调频信号转换为单频信号,借助搜索FFT峰值对应最大频率估计出起始频率。本发明另一个实施例中,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频,设定最佳阶次范围,通过峰值搜索得到不同阶次的峰值,获取N分量Chirp信号的调频率。进一步地,针对接收到的Chirp信号,获取其分数阶傅里叶变换,并通过单分量Chirp信号,获取Chirp信号在分数阶傅里叶变换域内变成一个脉冲的旋转角度,通过该旋转角度得到最佳阶次,依据最佳阶次获取对应分量的调频率。或,设定滑动步长进行不同阶次的峰值搜索得到最佳阶次,依据得到的最佳阶次获取多分量Chirp信号的调频率。
Chirp信号参数联合估计框架如图2所示,发射端发送N分量的线性调频信号,接收端分两步进行处理。上支路对接收到的信号做FRFT变换,最佳阶次的搜索范围为p∈[0,2],通过峰值搜索得到不同阶次的峰值,通过对应的阶次得到多分量Chirp信号的调频率的估计值。下支路接收到的信号与调频率分量做乘积,将对应调频率的信号转换为单频信号,搜索FFT峰值,通过峰值位置可得起始频率的估计值。
FRFT变换是一种时频分析方法,同时融合信号的时域和频域信息,对于给定的Chirp信号存在一个最佳分数阶次使其能量聚集于以最大值而表现为一个冲激函数,因此FRFT特别适用于处理Chirp类信号。时间信号x(t)的分数阶傅里叶变换定义如下:
Figure BDA0001828831930000071
其中,
Figure BDA0001828831930000072
Figure BDA0001828831930000073
称为分数阶傅里叶变换的核函数,其中:
Figure BDA0001828831930000074
是整数。通过变量代换
Figure BDA0001828831930000075
Figure BDA0001828831930000076
式(7)进一步地表示为:
Figure BDA0001828831930000077
式中,
Figure BDA0001828831930000078
称为幅度因子。将Chirp信号式(1)代入式(8)得:
Figure BDA0001828831930000081
当cota=-k时,Fp{s(t)}(u)=ABaexp(jπu2cotα)·δ[2π(ucscα-f0)],Chirp信号在旋转角度α=arccot(-k)的FRFT域内变成一个脉冲,相对应的阶次为p=2α/π称为最佳阶次,通过选取步长进行不同阶次的峰值搜索得到最佳阶次,此时
Figure BDA0001828831930000082
水声通信中选取的Chirp信号的调频率通常比较大,因此对应的最佳旋转角度趋于π,此时FRFT对于接收信号的调频率的微小变化非常敏感,这是借助FRFT精确估计调频率的一个前提。另外给出分数阶傅里叶变换的一个线性变换特性:
{Fp[ax(t)+by(t)]}(u)=aXp(u)+bYp(u) (10)
,借助FRFT在混合的线性调频信号中搜索出不同调频率。本发明再一个实施例中,获取N分量Chirp信号调频率后,对N分量Chirp信号依次进行调频补偿,包含如下内容:首先,将接收到的Chirp信号转换为单频信号;然后,针对每个分量Chirp信号的单频信号进行快速傅里叶变换,得到该分量Chirp信号对应的频率估计值序列,并联立信道时延和起始频率,得到起始频率估计值。
通过FRFT的峰值搜索得到信号的调频率,为了估计起始频率,本发明实施例中单频信号的转换过程可设计如下,在得到N分量Chirp信号的调频率[k1,k2,……kN]后,对接收信号R(t),借助(11)式将线性调频信号转换为单频信号进行处理。
Figure BDA0001828831930000083
以R1(t)为例进行推导说明,其余同理。
Figure BDA0001828831930000091
可见,(12)式中第一项对应于第一个Chirp信号经过转换后的分量,从FFT频域来看,该分量是起始频率为f01的信号经过L条多途时延后在频域形成一系列的峰值,峰值的位置包含初始频率和信道时延的信息。第二项对应于其余N-1个Chirp信号经过L条多途的输出信号,其本质上是多分量线性调频信号的线性叠加,不同之处在于调频率补偿后其调频率发生了变化,由于默认发送的线性调频信号调频率不同,因此ki-k1恒不等于零,Fourier变换后不同调频率的多个信号进行叠加后,其能量较为分散,ni'(t)为高斯白噪声与调频率分量的乘积,分布在整个频域内,而第一项为一系列的单频信号的叠加,其能量集中性较好,通过提取FFT峰值,进而得到起始频率。
对R1(t)进行FFT变换后得到频率估计值由大到小的排序为:
Figure BDA0001828831930000092
列出方程组:
Figure BDA0001828831930000093
取τ1=0,k1通过FRFT二维峰值搜索得到,
Figure BDA0001828831930000094
由FFT频域估计得到,通过联立信道时延和起始频率组成的方程组,得到f01,完成对起始频率的估计。其他分量Chirp信号的起始频率可通过同样估计过程获取。
基于上述的线性调频信号参数估计方法,本发明实施例还提供一种水声信道下线性调频信号参数估计装置,参见图3所示,包含模型构建模块101、数据获取模块102和数据输出模块103,其中,
模型构建模块101,用于构建用于输入并处理N分量线性调频Chirp信号的水声通道信号模型;
数据获取模块102,用于接收通过水声通道模型处理后的N分量Chirp信号,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频率;
数据输出模块103,用于对N分量Chirp信号依次进行调频补偿,并通过快速傅里叶变换对起始频率估计并输出。
上述的装置中,参见图4所示,数据输出模块103包含转换子模块201、估值子模块202和输出子模块203,其中,
转换子模块201,用于将接收到的Chirp信号转换为单频信号;
估值子模块202,用于针对每个分量Chirp信号的单频信号进行快速傅里叶变换,得到该分量Chirp信号对应的频率估计值序列,并联立信道时延和起始频率,得到起始频率估计值;
输出子模块203,用于将得到的起始频率估计值作为线性调频信号参数进行输出。
为验证本发明技术方案的有效性,下面通过仿真实验做进一步解释说明:
仿真实验一:在Matlab R2015b环境下进行了计算机仿真实验,首先验证本发明实施例中技术方案对单分量Chirp信号参数估计的有效性,水声通信中信号频率通常在1-30KHz,选取调频率和起始频率分别为20KHz/s和2000Hz作为发送信号,观测信号时长为T=0.2s,采样频率为8KHz。信道采用Bellhop射线模型模拟浅海水声信道,参数为接收机和发射机水平距离1Km,整个水体平均水深100m,发射机深度10m,接收机深度10m,声波频率50KHz,采用3条声线模拟海洋声信道中的直达声线,海面反射声线和海底反射声线,归一化后的冲激响应如图5所示,信噪比为-5dB,通过峰值搜索后得到对应峰值。
表1单分量Chirp信号参数估计结果
Figure BDA0001828831930000111
图6为不同阶数下FRFT的峰值搜索图,由图6可见Chirp信号在最佳阶次的脉冲压缩特性较好,在较低的信噪比(-5dB)下,借助FRFT可以准确的估计出信号的调频率。图7为用估计调频率补偿后FFT的频域图像,三个峰值的横坐标由右至左分别对应起始频率值、起始频率经过较小时延后的频率值和起始频率经过较大时延后的频率值。表1位单分量Chirp信号参数估计结果,对比了FRFT方法与本发明实施例中技术方案。可以看出,在仿真的水声Bellhop信道和高斯白噪声的背景下,本发明实施例中技术方案可以准确的估计出单分量Chirp信号的参数。与FRFT方法相比,本发明实施例中技术方案对起始频率的估计更加精确,就复杂度而言,FRFT在实际应用时是通过快速FFT方法实现的,因此两种方法的复杂度均为O(NlogN),仿真中两种方法的计算时间基本相同,验证了理论分析。
仿真实验二:对多分量Chirp信号参数估计进行仿真,考虑发送三个Chirp信号,幅度均为1,调频率和起始频率分别为:10000Hz/s、20000Hz/s、30000Hz/s和1500Hz、2000Hz、2500Hz,信道采用仿真实验一中的Bellhop水声信道,在-5dB的信噪比下进行参数估计。
表2多分量Chirp信号参数估计结果
Figure BDA0001828831930000112
由图8可见,在-5dB的信噪比下三个峰值比较突出,对应三个调频率的信号。图9
是表2中信号2的调频率补偿后的频谱图,与图7相比,在同样的信噪比下,其干扰量变多,这部分干扰主要来自于其余未被均衡成单频信号的线性调频信号的频域叠加。由表2可以看出,借助FRFT方法对多分量Chirp信号调频率的估计值较精确,原因在于分数傅里叶变换的对调频率变化的敏感特性和线性变换特性,微弱的调频率的变化会使最佳变换峰值的剧烈变化,使其特别适合处理Chirp类信号。另一方面在-5dB的信噪比下,估计准确是因为经过FRFT变换Chirp信号在最佳阶次的脉冲压缩特性,即形成一个冲激,而白噪声在任何阶次的FRFT都不会出现明显的能量聚集,抗噪能力强是由Chirp信号和FRFT变换的特性决定的,这是本发明实施例中技术方案采用FRFT来估计调频率的两个重要原因。在计算时间比较接近的情况下,本发明实施例中技术方案对起始频率的估计更加准确的原因在于:FRFT方法在对信号进行时频轴旋转到最佳角度时,估计到频域轴在最佳阶次的投影,受水声信道时延影响,频域信息损失较多,本发明实施例中技术方案在调频率补偿后,在频域轴得到起始频率,提升参数估计性能。
仿真实验三:对不同信噪比下参数估计性能的仿真。信道为实验1中的水声信道,发送实验2中的三个信号,信噪比从-10到5dB变化,步长为1dB,对比FRFT方法和简明分数阶傅里叶(CFRFT)方法,得到不同信噪比下起始频率的估计方差曲线图,如图10所示。由图10可见,所有方法的RMSE都随着信噪比的增大而减小。当信噪比低于-8dB时,RMSE随着信噪比的降低而急剧增大,CFRFT方法相对于FRFT方法约提升了0.5dB估计性能,本发明实施例中技术方案较CFRFT约提升了约0.5dB估计性能,因此,本发明实施例中技术方案能够有效提升不同信噪比下的起始频率的估计性能,尤其是低信噪比情况下。
通过上述实验进一步验证了本发明中,利用经过水声信道后Chirp信号的调频率不变的特性,通过FRFT变换进行阶次搜索得到多分量Chirp信号的调频率,并对接收到的信号依次进行调频率补偿后,借助FFT变换对受不同信道时延影响的起始频率进行区分,得到起始频率估计值;能够有效解决在低信噪比、低占空比下、欠采样等恶劣影响下以及多分量Chirp信号参数估计精确度问题,保证信号重构准确率,提高水声通信质量。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
结合本文中所公开的实施例描述的各实例的单元及方法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已按照功能一般性地描述了各示例的组成及步骤。这些功能是以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。本领域普通技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不认为超出本发明的范围。
本领域普通技术人员可以理解上述方法中的全部或部分步骤可通过程序来指令相关硬件完成,所述程序可以存储于计算机可读存储介质中,如:只读存储器、磁盘或光盘等。可选地,上述实施例的全部或部分步骤也可以使用一个或多个集成电路来实现,相应地,上述实施例中的各模块/单元可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。本发明不限制于任何特定形式的硬件和软件的结合。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (6)

1.一种水声信道下线性调频信号参数估计方法,其特征在于,包含如下内容:
构建用于输入并处理N分量线性调频Chirp信号的水声通道信号模型;
接收通过水声通道模型处理后的N分量Chirp信号,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频率;
对N分量Chirp信号依次进行调频补偿,并通过快速傅里叶变换对起始频率估计并输出;构建水声通道模型,包含如下内容:依据单分量Chirp信号,及消除多普勒效应水声多途信道的时变冲激响应,获取N分量Chirp信号经过多途时延的混合叠加信号模型;
通过信号幅度、起始频率和调频率来表示单分量Chirp信号,通过信道路径的幅度和时延来表示时变冲激响应;
单分量Chirp信号表示为:
Figure FDA0003126772520000011
A,f0和k分别为Chirp信号的幅度,起始频率和调频率;N分量Chirp信号经过水声信道后得到的接收信号表示为
Figure FDA0003126772520000012
k'和f0'为经过信道后的信号起始频率和调频率,
Figure FDA0003126772520000013
为相位,L为多途信道个数,n(t)为高斯白噪声;
为了估计起始频率,单频信号的转换过程设计如下:在得到N分量Chirp信号的调频率[k1,k2,……kN]后,对接收信号R(t),借助如下公式将线性调频信号转换为单频信号进行处理:
Figure FDA0003126772520000014
对R1(t)进行FFT变换后得到频率估计值由大到小的排序为:
Figure FDA0003126772520000015
方程组表示如下:
Figure FDA0003126772520000021
取τ1=0,k1通过FRFT二维峰值搜索得到,
Figure FDA0003126772520000022
由FFT频域估计得到,通过联立信道时延和起始频率组成的方程组,得到f01,完成对起始频率的估计,其他分量Chirp信号的起始频率可通过同样估计过程获取。
2.根据权利要求1所述的水声信道下线性调频信号参数估计方法,其特征在于,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频,包含如下内容:设定最佳阶次范围,通过峰值搜索得到不同阶次的峰值,获取N分量Chirp信号的调频率。
3.根据权利要求2所述的水声信道下线性调频信号参数估计方法,其特征在于,针对接收到的Chirp信号,获取其分数阶傅里叶变换,并通过单分量Chirp信号,获取Chirp信号在分数阶傅里叶变换域内变成一个脉冲的旋转角度,通过该旋转角度得到最佳阶次,依据最佳阶次获取对应分量的调频率。
4.根据权利要求2所述的水声信道下线性调频信号参数估计方法,其特征在于,设定滑动步长进行不同阶次的峰值搜索得到最佳阶次,依据得到的最佳阶次获取多分量Chirp信号的调频率。
5.一种水声信道下线性调频信号参数估计装置,其特征在于,基于权利要求1所述的方法实现,包含模型构建模块、数据获取模块和数据输出模块,其中,
模型构建模块,用于构建用于输入并处理N分量线性调频Chirp信号的水声通道信号模型;
数据获取模块,用于接收通过水声通道模型处理后的N分量Chirp信号,通过分数阶傅里叶变换进行阶次搜索,获取N分量Chirp信号调频率;
数据输出模块,用于对N分量Chirp信号依次进行调频补偿,并通过快速傅里叶变换对起始频率估计并输出。
6.根据权利要求5所述的水声信道下线性调频信号参数估计装置,其特征在于,数据输出模块包含转换子模块、估值子模块和输出子模块,其中,
转换子模块,用于将接收到的Chirp信号转换为单频信号;
估值子模块,用于针对每个分量Chirp信号的单频信号进行快速傅里叶变换,得到该分量Chirp信号对应的频率估计值序列,并联立信道时延和起始频率,得到起始频率估计值;
输出子模块,用于将得到的起始频率估计值作为线性调频信号参数进行输出。
CN201811196283.6A 2018-10-15 2018-10-15 水声信道下线性调频信号参数估计方法及装置 Active CN109510787B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811196283.6A CN109510787B (zh) 2018-10-15 2018-10-15 水声信道下线性调频信号参数估计方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811196283.6A CN109510787B (zh) 2018-10-15 2018-10-15 水声信道下线性调频信号参数估计方法及装置

Publications (2)

Publication Number Publication Date
CN109510787A CN109510787A (zh) 2019-03-22
CN109510787B true CN109510787B (zh) 2021-08-17

Family

ID=65746575

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811196283.6A Active CN109510787B (zh) 2018-10-15 2018-10-15 水声信道下线性调频信号参数估计方法及装置

Country Status (1)

Country Link
CN (1) CN109510787B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110808791B (zh) * 2019-10-15 2021-05-28 天津大学 基于分数阶傅里叶变换的线性调频信号的检测方法
CN110764062B (zh) * 2019-11-08 2020-10-13 中国人民解放军国防科技大学 基于分数阶傅里叶域滤波的多分量线性调频信号参数估计方法
CN110703260B (zh) * 2019-11-12 2023-01-17 南通赛洋电子有限公司 一种基于分数阶傅里叶变换的变频声呐测深方法
CN111935038B (zh) * 2020-08-03 2022-08-19 中国人民解放军国防科技大学 基于分数阶傅里叶变换的线性调频干扰消除方法
CN112866151B (zh) * 2021-01-18 2022-09-13 中国人民解放军战略支援部队信息工程大学 基于信道盲估计的水声MPSK信号盲Turbo均衡方法
CN113093115B (zh) * 2021-04-12 2022-04-15 中国人民解放军陆军工程大学 基于frft的低信噪比lfm信号参数快速估计方法
CN113132281B (zh) * 2021-06-11 2021-09-07 北京国科天迅科技有限公司 线性调频信号跟踪方法、装置、无线通信设备和存储介质
CN113721201B (zh) * 2021-09-08 2023-10-13 哈尔滨工程大学 一种线性调频信号调频率的估计方法
CN114338296B (zh) * 2021-11-24 2023-08-18 北京智芯微电子科技有限公司 信道估计和解调方法、装置及介质、接收设备、发送设备
CN114337982B (zh) * 2022-03-16 2022-07-15 深圳市华普微电子有限公司 一种啁啾导频信号的产生及时频估计方法
CN114866159B (zh) * 2022-04-01 2023-04-21 华南理工大学 一种多分量线性调频信号时频分析方法
CN115037329A (zh) * 2022-05-31 2022-09-09 江苏屹信航天科技有限公司 抗多普勒的接收方法、接收装置、终端
CN117118536B (zh) * 2023-10-25 2023-12-19 南京派格测控科技有限公司 调频稳定性的确定方法、装置、设备及存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101094217A (zh) * 2007-07-27 2007-12-26 哈尔滨工程大学 水声信道参数测量方法及移动水声通信同步获取方法
CN102680948A (zh) * 2012-05-15 2012-09-19 东南大学 一种线性调频信号调频率和起始频率估计方法
CN102832942A (zh) * 2011-06-16 2012-12-19 中国人民解放军海军航空工程学院 基于分数阶Fourier变换的三角线性调频连续波特征提取方法
CN103412295A (zh) * 2013-08-30 2013-11-27 西安电子科技大学 基于回波精确模型的高速机动弱目标检测方法
CN103701492A (zh) * 2014-01-14 2014-04-02 厦门大学 线性调频信号调制解调的水声跳频通信方法
CN103915102A (zh) * 2014-03-12 2014-07-09 哈尔滨工程大学 一种lfm水声多途信号的噪声抑制方法
CN104168232A (zh) * 2014-08-11 2014-11-26 厦门大学 一种水声信道中多径时延与多普勒频移的测定方法
CN105783974A (zh) * 2014-12-25 2016-07-20 中国科学院声学研究所 一种线性调频信号的检测、参数估计方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101094217A (zh) * 2007-07-27 2007-12-26 哈尔滨工程大学 水声信道参数测量方法及移动水声通信同步获取方法
CN102832942A (zh) * 2011-06-16 2012-12-19 中国人民解放军海军航空工程学院 基于分数阶Fourier变换的三角线性调频连续波特征提取方法
CN102680948A (zh) * 2012-05-15 2012-09-19 东南大学 一种线性调频信号调频率和起始频率估计方法
CN103412295A (zh) * 2013-08-30 2013-11-27 西安电子科技大学 基于回波精确模型的高速机动弱目标检测方法
CN103701492A (zh) * 2014-01-14 2014-04-02 厦门大学 线性调频信号调制解调的水声跳频通信方法
CN103915102A (zh) * 2014-03-12 2014-07-09 哈尔滨工程大学 一种lfm水声多途信号的噪声抑制方法
CN104168232A (zh) * 2014-08-11 2014-11-26 厦门大学 一种水声信道中多径时延与多普勒频移的测定方法
CN105783974A (zh) * 2014-12-25 2016-07-20 中国科学院声学研究所 一种线性调频信号的检测、参数估计方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于Radon-Ambiguity变换的机动目标检测方法;何延峰,刘峥,汤建新;《火控雷达技术》;20080630;第37卷(第2期);第21-24页 *

Also Published As

Publication number Publication date
CN109510787A (zh) 2019-03-22

Similar Documents

Publication Publication Date Title
CN109510787B (zh) 水声信道下线性调频信号参数估计方法及装置
CN204495996U (zh) 宽带声纳接收器
Asl et al. Contrast enhancement and robustness improvement of adaptive ultrasound imaging using forward-backward minimum variance beamforming
Worthmann et al. High frequency source localization in a shallow ocean sound channel using frequency difference matched field processing
Abadi et al. Blind deconvolution for robust signal estimation and approximate source localization
CN104678372B (zh) 正交频分复用雷达超分辨距离与角度值联合估计方法
CN109088835A (zh) 基于时间多重稀疏贝叶斯学习的水声时变信道估计方法
CN102879785B (zh) 基于频谱认知和分段跳频调频的水下物体探测方法及系统
CN107271955B (zh) 一种宽带线性调频信号的时差和尺度差估计方法
CN106936514A (zh) 稀疏信道下基于能量集中的水声前导信号检测方法
Josso et al. Source motion detection, estimation, and compensation for underwater acoustics inversion by wideband ambiguity lag-Doppler filtering
CN103618575B (zh) 一种调频水声通信系统的自动实时帧同步方法
CN116840821B (zh) 基于数据分析的双频测深控制系统
CN202995039U (zh) 基于频谱认知和分段跳频调频的水下物体探测系统
CN108957389A (zh) 一种实数域多通道信号目标方位估计方法
CN112435692A (zh) 一种基于多脉冲相干累加的信标信号高增益检测方法
Jiang et al. Discrete time-frequency characterizations of dispersive linear time-varying systems
CN104793197A (zh) 基于ifft频谱相除法和梯度自适应格型滤波的直达波抑制方法
CN103584886B (zh) 一种基于相位相干信息的自适应变迹方法
RU2677835C1 (ru) Способ разрешения тональных гидроакустических сигналов
CN108318883B (zh) 一种基于全频谱换能器自然加权的距离向计算方法
CN111722178A (zh) 一种基于指向性模型数值求解的远场窄带信号来波方向估计方法
Kochanska Testing the wide-sense stationarity of bandpass signals for underwater acoustic communications
Kebkal et al. Ultrasonic link for improved incoherent data transmission in horizontal shallow water channels
Hamschin et al. Waveform design with time and frequency constraints for optimal detection of elastic objects

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
GR01 Patent grant
GR01 Patent grant