CN107064942B - 一种高精度多普勒信号模拟方法 - Google Patents
一种高精度多普勒信号模拟方法 Download PDFInfo
- Publication number
- CN107064942B CN107064942B CN201611158289.5A CN201611158289A CN107064942B CN 107064942 B CN107064942 B CN 107064942B CN 201611158289 A CN201611158289 A CN 201611158289A CN 107064942 B CN107064942 B CN 107064942B
- Authority
- CN
- China
- Prior art keywords
- signal
- sample
- value
- signal sequence
- simulation
- 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
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
-
- 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
Abstract
本发明属于水声传播及水声信号处理技术领域,具体涉及水声信道中声波在(浮标水听器)接收端的高精度多普勒信号模拟方法。本发明包括:基于浮标阵的运动目标模拟,设计一组发射信号序列,利用时延计算,得到相应的接收信号序列;利用周期性样本逼近重构方法对接收信号序列进行预处理,由接收信号序列得到等间隔采样的信号序列;利用反馈及时延判断,并行多通道信号序列进行拼接整合,按周期依次存储数据,实现高精度多普勒信号的模拟。本发明多通道并行处理,提高了信号生成效率。多周期拼接信号,使长信号大数据量存储简单易行。以上两点使得算法用于实际应用设计时,电脑运行内存不再是主要的影响因子,从而扩展了算法的实用性。
Description
技术领域
本发明属于水声传播及水声信号处理技术领域,具体涉及水声信道中声波在(浮标水听器)接收端的高精度多普勒信号模拟方法。
背景技术
随着海洋科技的不断发展,水下定位技术成为如今发展研究的一大课题,对精度要求的不断提高使得目标模拟器应运而生。高精度多普勒信号是模拟器设计的重要内容,以往多普勒信号的仿真局限在依据多普勒效应在频域上进行频移,以实现时域上变宽窄的方法,计算量大,运算速度慢,使得高速目标信号精度受到很大限制。因此,采用有效措施实现高精度多普勒信号的模拟是模拟器精度提高的关键技术之一。
比较近年来水下目标模拟器信号设计主要存在两大问题。首先是硬件设备负责信号生成工作,并且不同系统之间目标种类和信号类型都有特殊要求,调试阶段要根据需求随之更新,这种情况下,要满足一定精度就需要较高的采样率,但这又使得硬件设备难以负荷,其可扩展性也受到了严重的制约。其次是生成多普勒信号的算法往往局限在利用多普勒效应的方法上,使得信号只能单点添加频移合成,这样大大降低了信号的生成效率和存储速率。李刚提出可以采用多路拼接方法将窄带信号拼接为宽带信号从而提高生成与存储数据的效率([3]李刚.基于多路拼接的宽带雷达信号源的研究.电子科技大学硕士学位论文.2007.),此方法中各路信号是独立的,其幅值、频率、相位均不连续,即使使用同一参考时钟源,信号质量也会受到一定的影响。王丽提出了可将插值方法用于解决少样本低精度的问题([1]王丽.具有抗差特点嫦娥一号异源数据地形滤波与插值方法研究.中国地质大学硕士学位论文.2013.),但此方法主要应用于构造三维空间形态或图像放大处理中,且插值选用的是空间插值思想。范展物理模型角度,提出了可将插值方法用于生成多普勒信号([4]范展.目标高速运动对矢量信号处理的影响及匹配技术研究.哈尔滨工程大学硕士学位论文.2009.),但没有对插值方法做具体的分析介绍以及实践应用。
发明内容
本发明的目的在于提出一种效果良好、计算简单、适用于在目标模拟中对高速运动目标高精度多普勒信号模拟方法。
本发明的目的是这样实现的:
(1)基于浮标阵的运动目标模拟,设计一组发射信号序列,利用时延计算,得到相应的接收信号序列;
(2)利用周期性样本逼近重构方法对接收信号序列进行预处理,由接收信号序列得到等间隔采样的信号序列;
(3)利用反馈及时延判断,并行多通道信号序列进行拼接整合,按周期依次存储数据,实现高精度多普勒信号的模拟。
所述的周期性样本逼近重构方法为:以一个周期为样本重构单位,比较所需采样点处时间轴附近的样本值,确定最接近的样本值从而确定采样点处的值,以此作为逼近样本的重构值,加之接收信号序列的时延零值部分,重构为每周期的样本值。
所述的利用周期性样本逼近重构方法实现由接收信号序列得到等间隔采样的信号序列的方法为:
如果tn为接收点(tr1,y(ts1))和(tr2,y(ts2))之间的样本逼近重构时刻,通过比较tn与tr1和tr2的差值的大小,来确定y(tn)的值,即:
|tr1-tn|≤|tr2-tn|→y(tn)=y(ts1) (a)
或者
|tr1-tn|>|tr2-tn|→y(tn)=y(ts2) (b)
通过式(a)和(b)可以得到均匀化的采样时刻对应的幅值,利用信道模型对该信号进一步处理可得到浮标接收端所需的接收信号序列。
所述的利用反馈信息进行并行多通道周期性的信号拼接方法为:信号进行周期性样本逼近重构步骤之后,将信号首脉冲到达接收端的初始时延t0与信号脉宽t1相加的和t与信号周期T进行比较;若t≤T,则将t时间的信号数据进行存储;若t>T,则截取信号前T时间数据进行存储,后(t-T)数据作为下一周期的信号判断数据;此时判断程序循环次数是否达到最大值,若未达到最大值,则重复以上步骤;若达到最大值,则将数据进行拼接并完成存储,并且拼接存储过程中均是信号截断处再次连接,保证了信号的连续性。
所述的周期性样本逼近重构方法可以对接收信号序列进行预处理。
所述的拼接整合方法可以对多通道信号序列进行预处理。
本发明的有益效果在于:(1)周期性样本逼近重构方法与目标信号的信号形式无关,也不受目标速度及运动方向的影响。其算法精度主要取决于对发射信号的采样率。采样率越高精度也就越高,但同时也增大了运算量。因此在实际的应用中要结合精度和运算量两方面考虑问题。(2)多通道并行处理,提高了信号生成效率。多周期拼接信号,使长信号大数据量存储简单易行。以上两点使得算法用于实际应用设计时,电脑运行内存不再是主要的影响因子,从而扩展了算法的实用性。
附图说明
图1为浮标阵定位系统工作态势图;
图2为周期性样本逼近重构方法过程图;
图3为高精度多普勒信号周期性整合流程图;
图4a为目标发射信号及其放大图;
图4b为含多普勒的接收信号及其放大图。
具体实施方式
下面结合附图对发明作进一步描述:
本发明包括:(1)基于浮标阵的运动目标模拟,设计一组发射信号序列,利用时延计算,得到对应的接收信号序列。
(2)利用周期性样本逼近重构方法,由接收信号序列得到等间隔采样的信号序列。
(3)利用反馈及时延判断,并行多通道信号序列进行拼接整合,按重构周期依次存储,实现高精度多普勒信号的模拟。
本发明的主要技术特征是:
1.多普勒信号合成方法采用周期性样本逼近重构的方法、信号拼接采用多通道并行处理方式。
2.本发明所述的周期性样本逼近重构方法是以一个周期为样本重构单位,比较所需采样点处时间轴附近的样本值,确定最接近的样本值从而确定采样点处的值,以此作为逼近样本的重构值,加之接收信号序列的时延零值部分,重构为每周期的样本值。
3.本发明所述的拼接整合方法是以一个周期为拼接整合单位,每周期数据包括时延零数据部分和信号脉冲部分,即每周期的样本值。首先存储等于一个周期的数据,其次进行判断,以大于周期的数据作为下一周期进行叠加的反馈数据,以此反复叠加,最终实现对高精度多普勒信号多通道数据的存储。
实施条件:浮标阵为800m*800m的方阵,浮标分为位于阵的四个角位置,坐标依次为(0,0)、(800,0)、(800,800)、(0,800),如图1所示。信号为单频余弦信号,频率f0=10kHz,幅值A=1V,脉宽为50ms,采样频率fs=200kHz。目标数目为一个,速度v=10m/s,位于海深200m的位置,且在该平面做直线运动,起点为N(0,300),终点为M(200,700)。
针对上述条件,给出该具体实施例高精度多普勒信号生成步骤,详见图2、图3。
1.基于浮标阵的接收信号序列模拟原理
浮标阵如图1所示,设定一个800m*800m的方阵,四个浮标位于方阵的四个角,目标轨迹设置其起点为N(0,300),终点为M(200,700),且在海深为200m的平面做速度为v=10m/s的匀速直线运动,以浮标1所在位置(0,0)为原点,相邻两浮标方向分别为x方向和y方向建立坐标系,则可依据该数学物理模型利用公式
求出对应目标发射信号的采样时间的浮标接收端接收信号的时刻,式中,ts表示目标发射信号的采样时间,以T为采样间隔,称为信号周期,tr表示浮标接收端接收信号的时刻,(x2(ts),y2(ts))表示目标在ts时刻的位置坐标,c为声速。
2.周期性样本逼近重构方法的基本原理
利用公式(1)求得的浮标端接收信号接收时刻为非均匀分布的时间点,而目标模拟器的目的是模拟浮标等信号接收端模拟信号的仪器,是将目标发射的模拟信号经处理后,得到的数字接收序列经DA转换为模拟信号输出。因此需要对以上信号进行间隔均匀化处理,方法为周期性样本逼近重构方法,如图2所示。
如果tn为接收点(tr1,y(ts1))和(tr2,y(ts2))之间的样本逼近重构时刻,通过比较tn与tr1和tr2的差值的大小,来确定y(tn)的值,即:
|tr1-tn|≤|tr2-tn|→y(tn)=y(ts1) (2)
或者
|tr1-tn|>|tr2-tn|→y(tn)=y(ts2) (3)
通过式(2)和(3)可以得到均匀化的采样时刻对应的幅值,利用信道模型对该信号进一步处理可得到浮标接收端所需的接收信号序列。
3.并行多通道信号拼接原理
以往目标模拟器对信号的仿真受数据量大的影响,会出现占用内存过大存储失败的后果,且单点生成速度慢,影响调试效率。
本发明设计了一种利用反馈信息进行并行多通道周期性的信号拼接方法,过程如图3所示:信号进行周期性样本逼近重构步骤之后,将信号首脉冲到达接收端的初始时延t0与信号脉宽t1相加的和t与信号周期T进行比较。若t≤T,则将t时间的信号数据进行存储;若t>T,则截取信号前T时间数据进行存储,后(t-T)数据作为下一周期的信号判断数据并按下一周期数据长度补零与之相加,得到下一周期的数据。此时判断程序循环次数是否达到最大值,若未达到最大值,则重复以上步骤;若达到最大值,则将数据进行拼接并完成存储,并且拼接存储过程中均是信号截断处再次连接,保证了信号的连续性。
此时已完成对高精度多普勒信号的模拟,结果如图4(b)所示,与图4(a)目标发射信号进行对比,周期显示为1.5s,不含多普勒信息。图4(b)数据生成以浮标1为接收端为例,周期为30个,通过MATLAB仿真验证所得初始时延值0.2s,以及含多普勒信号均与与理论值一致。
以上结果表明多通道并行处理方法,提高了信号生成效率。多周期拼接信号,使长信号大数据量存储简单易行。以上两点使得算法用于实际应用设计时,电脑运行内存不再是主要的影响因子,从而扩展了算法的实用性。
Claims (6)
1.一种高精度多普勒信号模拟方法,其特征在于,包括如下步骤:
(1)基于浮标阵的运动目标模拟,设计一组发射信号序列,利用时延计算,得到相应的接收信号序列;
(2)利用周期性样本逼近重构方法对接收信号序列进行预处理,由接收信号序列得到等间隔采样的信号序列;
(3)利用反馈及时延判断,并行多通道信号序列进行拼接整合,按周期依次存储数据,实现高精度多普勒信号的模拟。
2.根据权利要求1所述的一种高精度多普勒信号模拟方法,其特征在于:所述的周期性样本逼近重构方法为:以一个周期为样本重构单位,比较所需采样点处时间轴附近的样本值,确定最接近的样本值从而确定采样点处的值,以此作为逼近样本的重构值,加之接收信号序列的时延零值部分,重构为每周期的样本值。
3.根据权利要求1所述的一种高精度多普勒信号模拟方法,其特征在于:所述的利用周期性样本逼近重构方法实现由接收信号序列得到等间隔采样的信号序列的方法为:
如果tn为接收点(tr1,y(ts1))和(tr2,y(ts2))之间的样本逼近重构时刻,通过比较tn与tr1和tr2的差值的大小,来确定y(tn)的值,即:
|tr1-tn|≤|tr2-tn|→y(tn)=y(ts1) (a)
或者
|tr1-tn|>|tr2-tn|→y(tn)=y(ts2) (b)
通过式(a)和(b)可以得到均匀化的采样时刻对应的幅值,利用信道模型对该信号进一步处理可得到浮标接收端所需的接收信号序列。
4.根据权利要求1所述的一种高精度多普勒信号模拟方法,其特征在于:所述的并行多通道信号序列进行拼接整合的方法为:信号进行周期性样本逼近重构步骤之后,将信号首脉冲到达接收端的初始时延t0与信号脉宽t1相加的和t与信号周期T进行比较;若t≤T,则将t时间的信号数据进行存储;若t>T,则截取信号前T时间数据进行存储,后(t-T)数据作为下一周期的信号判断数据;此时判断程序循环次数是否达到最大值,若未达到最大值,则重复以上步骤;若达到最大值,则将数据进行拼接并完成存储,并且拼接存储过程中均是信号截断处再次连接,保证了信号的连续性。
5.根据权利要求1所述的一种高精度多普勒信号模拟方法,其特征在于:所述的周期性样本逼近重构方法可以对接收信号序列进行预处理。
6.根据权利要求1所述的一种高精度多普勒信号模拟方法,其特征在于:所述的拼接整合方法可以对多通道信号序列进行预处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611158289.5A CN107064942B (zh) | 2016-12-15 | 2016-12-15 | 一种高精度多普勒信号模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611158289.5A CN107064942B (zh) | 2016-12-15 | 2016-12-15 | 一种高精度多普勒信号模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064942A CN107064942A (zh) | 2017-08-18 |
CN107064942B true CN107064942B (zh) | 2019-12-06 |
Family
ID=59619329
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611158289.5A Active CN107064942B (zh) | 2016-12-15 | 2016-12-15 | 一种高精度多普勒信号模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107064942B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732553B (zh) * | 2018-06-01 | 2022-02-01 | 北京航空航天大学 | 一种激光雷达波形时刻鉴别方法与在线测距系统 |
CN114218764B (zh) * | 2021-11-25 | 2022-10-21 | 哈尔滨工程大学 | 水下运动声源动态声场模拟计算方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101634708A (zh) * | 2009-07-13 | 2010-01-27 | 中国船舶重工集团公司第七一五研究所 | 一种主动声纳目标模拟器信号处理方法 |
CN102183741A (zh) * | 2011-02-28 | 2011-09-14 | 哈尔滨工程大学 | 在长基线异步水声定位系统实现高帧率无模糊定位的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5423713B2 (ja) * | 2011-03-30 | 2014-02-19 | 日本電気株式会社 | アクティブソーナー装置及びその信号処理方法 |
-
2016
- 2016-12-15 CN CN201611158289.5A patent/CN107064942B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101634708A (zh) * | 2009-07-13 | 2010-01-27 | 中国船舶重工集团公司第七一五研究所 | 一种主动声纳目标模拟器信号处理方法 |
CN102183741A (zh) * | 2011-02-28 | 2011-09-14 | 哈尔滨工程大学 | 在长基线异步水声定位系统实现高帧率无模糊定位的方法 |
Non-Patent Citations (1)
Title |
---|
对抗条件下浮标基水声定位技术研究;付进等;《兵工学报》;20120229;第33卷(第2期);第192-196页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107064942A (zh) | 2017-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104007418B (zh) | 一种基于时间同步的大基阵水下宽带扩频信标导航定位系统及方法 | |
Kessinger | Extended split-step Fourier migration | |
CN108089155B (zh) | 一种深海环境下单水听器声源被动定位方法 | |
GB2342165A (en) | 3D prestack seismic data migration method | |
CN101893698B (zh) | 噪声源测试分析方法及其装置 | |
CN103364770B (zh) | 基于矩阵填充的雷达目标检测系统及其检测方法 | |
CN108459087B (zh) | 一种应用于板结构损伤检测的多模态Lamb波模态分离方法 | |
CN107219513A (zh) | 水下运动目标回波模拟方法及收发一体化装置 | |
CN107064942B (zh) | 一种高精度多普勒信号模拟方法 | |
CN103809170A (zh) | 一种多普勒频偏修正方法及系统 | |
CN116068493A (zh) | 一种深海大深度垂直分布式水听器被动声源定位方法 | |
CN103776907A (zh) | 基于sinc插值的超声相控阵接收信号精延时方法 | |
CN101576618A (zh) | 基于小波变换的水声定位测量方法和测量系统 | |
CN110780340B (zh) | 一种浅海弹性软沉积层横波声速或厚度的反演方法 | |
CN109738944B (zh) | 基于广角反射的地震采集参数确定方法及装置 | |
CN113126029B (zh) | 适用于深海可靠声路径环境的多传感器脉冲声源定位方法 | |
CN106918809B (zh) | 快速干涉合成孔径声纳原始回波时域仿真方法 | |
CN108494513A (zh) | 浅海水声信道模型建立及其计算方法 | |
Michalopoulou et al. | Sediment sound speed inversion with time-frequency analysis and modal arrival time probability density functions | |
CN103308900B (zh) | 一种用于目标检测的快速keystone变换方法 | |
CN105631149A (zh) | 一种基于相关性求时差的方法 | |
CN114859404A (zh) | 超采样地震波形匹配方法及装置 | |
CN101762811B (zh) | 一种基于面元散射的合成孔径声纳面目标快速仿真方法 | |
An et al. | Low resolution fourier synthesis modelling for underwater acoustic channel impulse response | |
Bonito | Acoustic system for ground truth underwater positioning in DEEC's test tank |
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 |