CN109031260B - 一种基于分数傅里叶调制率分析的lfm信号时延测量方法 - Google Patents

一种基于分数傅里叶调制率分析的lfm信号时延测量方法 Download PDF

Info

Publication number
CN109031260B
CN109031260B CN201810685402.8A CN201810685402A CN109031260B CN 109031260 B CN109031260 B CN 109031260B CN 201810685402 A CN201810685402 A CN 201810685402A CN 109031260 B CN109031260 B CN 109031260B
Authority
CN
China
Prior art keywords
signal
time delay
method based
fractional fourier
modulation
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
CN201810685402.8A
Other languages
English (en)
Other versions
CN109031260A (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201810685402.8A priority Critical patent/CN109031260B/zh
Publication of CN109031260A publication Critical patent/CN109031260A/zh
Application granted granted Critical
Publication of CN109031260B publication Critical patent/CN109031260B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于分数傅里叶调制率分析的LFM信号时延测量方法，本发明可由FrFT及已估计出的接收信号的调制斜率及多普勒系数，直接推导出时延及幅度参数；在不同的时间带宽的条件下，该方法都具有一定的稳健性，且在0dB左右的低信噪比条件下，依旧保持较好的性能。

Description

一种基于分数傅里叶调制率分析的LFM信号时延测量方法
技术领域
本发明涉及水声信号参数测量技术领域，尤其是一种基于分数傅里叶调制率分析的LFM信号时延测量方法。
背景技术
在主动声纳领域中，接收的水声信号为目标反射的回波，通过对回波信号的多普勒系数、时延及幅度的参数测量，可以计算出目标声纳之间的相对速度及距离。在水声通信领域中，通过发射测试信号，分析出接收信号中参数的变化，可以估计传输信道，完成信道均衡的目的。线性调频信号(LFM)是各水声领域中最常用的信号之一，因此，实现该类接收信号参数的快速、高精度测量具有十分重要的意义。
分数阶傅里叶变换(FrFT)是一种重要的时频分析工具，是傅里叶变换的一种广义形式，信号的FrFT可以解释为信号在时频平面内绕原点旋转任意角度后所构成的分数阶傅里叶域上的表示。LFM信号通过分数傅里叶变换，能在分数域上形成一个汇聚点，且汇聚点的坐标与信号的调制斜率具有对应关系，可用于白噪声下LFM信号的增强。Ozaktas提出了处理离散信号的快速FrFT方法。此方法经常被用来计算LFM信号的调制斜率及多普勒系数，而后续一般利用其他方法得出信号的时延及幅度，或者舍弃对信号时延值的测量。额外的计算将占用系统运算空间，降低效率，因此，需要更加有效便捷的方法测量接收信号的时延及幅度。
发明内容
本发明所要解决的技术问题在于，提供一种基于分数傅里叶调制率分析的LFM信号时延测量方法，能够在不同的时间带宽的条件下，依旧保持较好的性能。
为解决上述技术问题，本发明提供一种基于分数傅里叶调制率分析的LFM信号时延测量方法，包括如下步骤：
(1)发射信号为线型调频信号：
s[n]＝cos(2πf0·n/fs+πk·(n/fs)2)，n∈[0,T·fs]
其中，f0为信号起始频，k为信号的调制斜率，T为信号脉宽，fs为采样率；将接收信号r进行变阶数快速离散FrFT处理得到结果：Fpr，其中阶数p的搜索变化范围为[0,2]，离散化后的时间坐标轴n转化为u轴：u＝(n/fs-TWL/2)/S，n∈[0,TWL·fs]，其中
Figure BDA0001711551110000021
为单位化间隔，TWL为接收信号总时长；
(2)搜索出Fpr,p∈[0,2]中最大值所对应的阶数为最优阶数：
Figure BDA0001711551110000022
(3)在离散归一化坐标系中，提取
Figure BDA0001711551110000023
最大值所对应的u轴坐标：
Figure BDA0001711551110000024
则信号估计时延为：
Figure BDA0001711551110000025
(4)构造还原单位幅度的接收信号：
Figure BDA0001711551110000026
其中
Figure BDA0001711551110000027
表示向上取整，则其估计幅度为：
Figure BDA0001711551110000028
(5)得到线性调频信号的多普勒系数
Figure BDA0001711551110000029
时延
Figure BDA00017115511100000210
及幅度
Figure BDA00017115511100000211
三个接收信号参数的测量值。
优选的，步骤(2)中，最优旋转角为φ*＝πp*/2，则调制斜率估计值为：
Figure BDA00017115511100000212
多普勒系数估计值为
Figure BDA00017115511100000213
本发明的有益效果为：本发明可由FrFT及已估计出的接收信号的调制斜率及多普勒系数，直接推导出时延及幅度参数；在不同的时间带宽的条件下，该方法都具有一定的稳健性，且在0dB左右的低信噪比条件下，依旧保持较好的性能。
附图说明
图1为本发明的方法流程示意图。
图2为本发明的实施例中接收LFM信号时域波形图。
图3为本发明的实施例中接收LFM信号在归一化FrFT坐标系中的示意图。
具体实施方式
如图1所示，一种基于分数傅里叶调制率分析的LFM信号时延测量方法，包括如下步骤：
(1)发射信号为线型调频信号：
s[n]＝cos(2πf0·n/fs+πk·(n/fs)2)，n∈[0,T·fs]
其中，f0为信号起始频，k为信号的调制斜率，T为信号脉宽，fs为采样率；将接收信号r进行变阶数快速离散FrFT处理得到结果：Fpr，其中阶数p的搜索变化范围为[0,2]，离散化后的时间坐标轴n转化为u轴：u＝(n/fs-TWL/2)/S，n∈[0,TWL·fs]，其中
Figure BDA0001711551110000031
为单位化间隔，TWL为接收信号总时长；
(2)搜索出Fpr,p∈[0,2]中最大值所对应的阶数为最优阶数：
Figure BDA0001711551110000032
相应的最优旋转角为φ*＝πp*/2，则调制斜率估计值为：
Figure BDA0001711551110000033
多普勒系数估计值为
Figure BDA0001711551110000034
(3)在离散归一化坐标系中，提取
Figure BDA0001711551110000035
最大值所对应的u轴坐标：
Figure BDA0001711551110000036
则信号估计时延为：
Figure BDA0001711551110000037
(4)构造还原单位幅度的接收信号：
Figure BDA0001711551110000038
其中
Figure BDA0001711551110000039
表示向上取整，则其估计幅度为：
Figure BDA0001711551110000041
(5)得到线性调频信号的多普勒系数
Figure BDA0001711551110000042
时延
Figure BDA0001711551110000043
及幅度
Figure BDA0001711551110000044
三个接收信号参数的测量值。
下面对本实施例进行仿真验证，采用MATLAB软件，设定发射LFM信号参数：脉宽T＝2s，频率范围f0＝1200Hz，调制斜率k＝300，声速c＝1500m/s，采样率fs＝6000Hz，信噪比为5dB，处理窗长TWL＝5·fs。仿真中接收信号为：
r[n]＝r0[n]+ns[n]，n∈[0,TWL·fs]
其中，ns[n]为随机噪声，
Figure BDA0001711551110000045
其中，幅度A＝1，时延τ＝1.06674s，多普勒系数η＝0.99468。具体的时域波形图如图2所示。
利用离散信号变分数阶FrFT的方法计算处理该实施例中的接收LFM信号，如图3所示。图3为信号在归一化FrFT坐标系中的示意图，信号在该平面中呈X型，且最大值为交叉点处，则搜索出最优阶数为p*＝1.2490，相对应的最优旋转角为cscφ*＝1.9619，多普勒系数估计值为
Figure BDA0001711551110000046
在离散归一化坐标系中，单位化间隔S＝0.0289，提取
Figure BDA0001711551110000047
最大值所对应的u轴坐标为u*＝45.4663，代入以下公式中：
Figure BDA0001711551110000048
则信号估计时延为
Figure BDA0001711551110000049
再由以下公式：
Figure BDA00017115511100000410
估计出信号幅度为
Figure BDA00017115511100000411
由
Figure BDA00017115511100000412
计算估计误差，得多普勒系数误差εη＝0.0056％，时延误差ετ＝0.0444％，幅度误差εA＝1.09％。
由实验数据处理结果可以得出结论，该基于FrFT的方法对LFM信号的时延及幅度参数都具有较准确的估计。
本发明的方法在对处理时效性上有较高要求时，更加节约计算空间；信号时延的测量精度将随多普勒系数的测量精度的提高而提高，完成信号参数测量精度的整体性提升。

Claims (2)

1.一种基于分数傅里叶调制率分析的LFM信号时延测量方法，其特征在于，包括如下步骤：
(1)发射信号为线型调频信号：
s[n]＝cos(2πf0·n/fs+πk·(n/fs)2)，n∈[0,T·fs]
其中，f0为信号起始频，k为信号的调制斜率，T为信号脉宽，fs为采样率；将接收信号r进行变阶数快速离散FrFT处理得到结果：Fpr，其中阶数p的搜索变化范围为[0,2]，离散化后的时间坐标轴n转化为u轴：u＝(n/fs-TWL/2)/S，n∈[0,TWL·fs]，其中
Figure FDA0003498037840000011
为单位化间隔，TWL为接收信号总时长；
(2)搜索出Fpr,p∈[0,2]中最大值所对应的阶数为最优阶数：
Figure FDA0003498037840000012
(3)在离散归一化坐标系中，提取
Figure FDA0003498037840000018
最大值所对应的u轴坐标：
Figure FDA0003498037840000013
p∈[0,2]
则信号估计时延为：
Figure FDA0003498037840000014
其中，φ*为最优旋转角；
(4)构造还原单位幅度的接收信号：
Figure FDA0003498037840000015
其中
Figure FDA0003498037840000016
表示向上取整，则其估计幅度为：
Figure FDA0003498037840000017
(5)得到线性调频信号的多普勒系数
Figure FDA0003498037840000021
时延
Figure FDA0003498037840000022
及幅度
Figure FDA0003498037840000023
三个接收信号参数的测量值。
2.如权利要求1所述的基于分数傅里叶调制率分析的LFM信号时延测量方法，其特征在于，步骤(3)中，最优旋转角为φ*＝πp*/2，则调制斜率估计值为：
Figure FDA0003498037840000024
多普勒系数估计值为
Figure FDA0003498037840000025
CN201810685402.8A 2018-06-28 2018-06-28 一种基于分数傅里叶调制率分析的lfm信号时延测量方法 Active CN109031260B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810685402.8A CN109031260B (zh) 2018-06-28 2018-06-28 一种基于分数傅里叶调制率分析的lfm信号时延测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810685402.8A CN109031260B (zh) 2018-06-28 2018-06-28 一种基于分数傅里叶调制率分析的lfm信号时延测量方法

Publications (2)

Publication Number Publication Date
CN109031260A CN109031260A (zh) 2018-12-18
CN109031260B true CN109031260B (zh) 2022-04-26

Family

ID=65520682

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810685402.8A Active CN109031260B (zh) 2018-06-28 2018-06-28 一种基于分数傅里叶调制率分析的lfm信号时延测量方法

Country Status (1)

Country Link
CN (1) CN109031260B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133598B (zh) * 2019-05-09 2023-06-23 西安电子科技大学 基于FrFT的线性调频信号参数快速估计方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102866391A (zh) * 2012-09-05 2013-01-09 中北大学 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法
CN102999473A (zh) * 2012-10-18 2013-03-27 中国人民解放军电子工程学院 一种线性调频信号的检测与参数估计方法
CN103532656A (zh) * 2013-08-08 2014-01-22 北京理工大学 基于分数阶傅里叶域信道化的宽带线性调频信号多假目标干扰方法
CN103901470A (zh) * 2014-03-21 2014-07-02 哈尔滨工程大学 一种海底沉积层反射强度及时延估计方法
CN105527614A (zh) * 2016-02-29 2016-04-27 南京理工大学 基于动目标检测的lfm脉冲雷达抗drfm干扰方法
CN106249208A (zh) * 2016-07-11 2016-12-21 西安电子科技大学 基于分数阶傅里叶变换的噪声调幅干扰下信号检测方法
CN106443588A (zh) * 2016-05-23 2017-02-22 中国人民解放军63892部队 一种lfmcw信号快速检测和估计方法
KR20170029899A (ko) * 2015-09-08 2017-03-16 국방과학연구소 능동 소나 시스템의 lfm신호 파라미터 추정 장치 및 방법
CN106646436A (zh) * 2016-12-09 2017-05-10 东南大学 一种基于信号宽窄带模糊度的侦察信号参数估计方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102866391A (zh) * 2012-09-05 2013-01-09 中北大学 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法
CN102999473A (zh) * 2012-10-18 2013-03-27 中国人民解放军电子工程学院 一种线性调频信号的检测与参数估计方法
CN103532656A (zh) * 2013-08-08 2014-01-22 北京理工大学 基于分数阶傅里叶域信道化的宽带线性调频信号多假目标干扰方法
CN103901470A (zh) * 2014-03-21 2014-07-02 哈尔滨工程大学 一种海底沉积层反射强度及时延估计方法
KR20170029899A (ko) * 2015-09-08 2017-03-16 국방과학연구소 능동 소나 시스템의 lfm신호 파라미터 추정 장치 및 방법
CN105527614A (zh) * 2016-02-29 2016-04-27 南京理工大学 基于动目标检测的lfm脉冲雷达抗drfm干扰方法
CN106443588A (zh) * 2016-05-23 2017-02-22 中国人民解放军63892部队 一种lfmcw信号快速检测和估计方法
CN106249208A (zh) * 2016-07-11 2016-12-21 西安电子科技大学 基于分数阶傅里叶变换的噪声调幅干扰下信号检测方法
CN106646436A (zh) * 2016-12-09 2017-05-10 东南大学 一种基于信号宽窄带模糊度的侦察信号参数估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"一种单频水声信号多径时延估计算法";童峰等;《声学学报》;20080131;第33卷(第1期);第62-68页 *

Also Published As

Publication number Publication date
CN109031260A (zh) 2018-12-18

Similar Documents

Publication Publication Date Title
CN104237871B (zh) 一种基于相位补偿的时延差估计方法
CN107179535A (zh) 一种基于畸变拖曳阵的保真增强波束形成的方法
CN104820218B (zh) 一种基于频域自相关的浅海海底单参数反演方法
CN106019263B (zh) 基于多亮点回波模型的目标径向速度测量方法
Zedel et al. Resolving velocity ambiguity in multifrequency, pulse-to-pulse coherent Doppler sonar
CN107085216A (zh) 一种基于单水听器的深海水声被动测距测深方法
CN103630706A (zh) 一种在声学多普勒流速剖面仪中获取径向水流速度的方法
Wei et al. Doppler estimation based on dual-HFM signal and speed spectrum scanning
CN106546947A (zh) 一种联合波导不变量和线谱的单水听器被动定位方法
CN104748704A (zh) 薄壁结构超声共振测厚频谱分析内插校正方法
CN110488219A (zh) 一种浅海单水听器运动目标距离估计方法
RU2550576C1 (ru) Способ измерения дистанции до шумящего объекта
CN109031260B (zh) 一种基于分数傅里叶调制率分析的lfm信号时延测量方法
CN108646248A (zh) 一种针对低速运动声源的被动声学测速测距方法
KR20170029899A (ko) 능동 소나 시스템의 lfm신호 파라미터 추정 장치 및 방법
CN110048741A (zh) 一种基于短时分数阶傅里叶变换的跳频信号的参数估计方法
CN111830481B (zh) 雷达回波单分量幅度分布模型参数估计方法及装置
CN109471113B (zh) 基于相位法的多波束声呐海底地形测量质量实时评估方法
CN108957416A (zh) 脉冲噪声环境下基于分数阶功率谱密度的线性调频信号参数估计方法
CN110108797B (zh) 利用声阻抗变化信息的介质分界面超声检测方法
CN110470253A (zh) 基于超声波的厚度测量方法、装置、电子设备及存储介质
KR101480170B1 (ko) 음향 도플러를 이용한 수중 이동체속도 추정방법 및 이를 이용한 수중 이동체속도 측정장치
CN112083466B (zh) 一种顾及时间偏差的海底应答器定位方法及系统
CN113126029B (zh) 适用于深海可靠声路径环境的多传感器脉冲声源定位方法
CN106646436B (zh) 一种基于信号宽窄带模糊度的侦察信号参数估计方法

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