CN110808791B - 基于分数阶傅里叶变换的线性调频信号的检测方法 - Google Patents

基于分数阶傅里叶变换的线性调频信号的检测方法 Download PDF

Info

Publication number
CN110808791B
CN110808791B CN201910978041.0A CN201910978041A CN110808791B CN 110808791 B CN110808791 B CN 110808791B CN 201910978041 A CN201910978041 A CN 201910978041A CN 110808791 B CN110808791 B CN 110808791B
Authority
CN
China
Prior art keywords
signal
fourier transform
data
positions
lfm
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
CN201910978041.0A
Other languages
English (en)
Other versions
CN110808791A (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN201910978041.0A priority Critical patent/CN110808791B/zh
Publication of CN110808791A publication Critical patent/CN110808791A/zh
Application granted granted Critical
Publication of CN110808791B publication Critical patent/CN110808791B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
    • H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L7/00—Arrangements for synchronising receiver with transmitter
    • H04L7/04—Speed or phase control by synchronisation signals
    • H04L7/10—Arrangements for initial synchronisation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0044—Control loops for carrier regulation
    • H04L2027/0071—Control of loops
    • H04L2027/0073—Detection of synchronisation state
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0083—Signalling arrangements
    • H04L2027/0097—Adaptive synchronisation signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • General Physics & Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种基于分数阶傅里叶变换的线性调频信号的检测方法，包括下列的步骤：通信过程中以收发两端均已知的一段线性调频LFM信号作为同步信号，将接收到的信号与本地信号进行相关检测，将所得相关信号保存；进行p阶次的分数阶傅里叶变换，搜索本地信号对应的最佳阶次；接收端对相关信号中所有数据均进行此操作，确定适合背景噪声的判决门限；遍历相关信号，找出大于判决门限的数据并记录此数据在接收信号中对应的位置；对所记录的每个位置之后的一段信号进行分数阶傅里叶变换，寻找作为经分数阶傅里叶变换检测所得的同步信号起始位置；将记录的疑似信号起始位置之后一段信号再采用分段傅里叶变换进行二次筛选。

Description

基于分数阶傅里叶变换的线性调频信号的检测方法
技术领域
本发明属于水声通信数字信号处理技术领域，涉及一种同步信号起始位置检测技术。
背景技术
通信系统中同步方式主要有自同步和外同步两种方式，外同步是水声通信最常用的方式。其做法是在每帧信号的起始位置加入接收双方均已知一段信号。这段信号通常是具有良好自相关特性的线性调频(LFM)信号，接收端根据已知的线性调频信号相关参数对接收的信号进行相关检测得到相关信号，找到相关信号中的幅值最大值点在信号中对应的位置作为每帧信号的起始位置然后开始解码。但是由于水声信道十分复杂，通信环境也不可预测，水声通信系统的背景噪声存在非平稳特性，并且非平稳特性噪声强度往往高于背景噪声。通过相关之后的输出的相关包络变得起伏较大，影响了信号起始位置的判定的准确性，从而导致接收端按照错误的起始位置开始解调，出现一连串解调错误的现象。
针对上述问题，本发明提出了一种基于分数阶傅里叶变换和分段阶傅里叶变换相结合的线性调频信号检测的方法，对经过相关后的相关信号进行两次起始位置的判决，首先利用设定的自适应门限算法找到门限阈值，然后利用线性调频信号经过分数阶傅里叶变换后出现的能量聚集现象和经分段傅里叶变换后频率的递增特性，对过门限位置之后的一段信号做两次判决，提高了信号起始位置判别的准确度。
发明内容
为了改善上述指出的信号起始位置判决不准确的情况，本发明提出了分数阶傅里叶变换和分段阶傅里叶变换相结合的同步信号起始位置检测方法。利用已知线性调频信号的经过分数阶傅里叶变换出现的最佳阶次和调频斜率，同时将所得的调频斜率与已知的正确斜率进行比对，随后利用LFM信号经过分段傅里叶变换之后出现的递增特性，对信号起始位置做出判断，从而检测到真正的线性调频信号确定信号起始位置。技术方案如下：
一种基于分数阶傅里叶变换的线性调频信号的检测方法，包括下列的步骤：
(1)通信过程中以收发两端均已知的一段线性调频LFM信号作为同步信号，设此信号频率范围为fst～fe,调频斜率为kt，信号时长为tlfm，接收端首先根据此参数生成LFM信号作为本地信号，然后将接收到的信号与本地信号进行相关检测，将所得相关信号c的数据保存。
(2)对(1)中生成的本地信号进行p阶次的分数阶傅里叶变换，搜索本地信号对应的最佳阶次p0即经分数阶傅里叶变换后最大模值对应的阶次；
(3)接收端生成长度为N的数组a[N]初始化全为0，并设置判决门限的初始数值threshold，将相关信号当前数据c(i)赋给数组中最后一个数据a[N-1]，判断此数据是否大于数组a中所有数据中位数的K倍，若大于则按照公式threshold＝0.8*threshold+0.2*c(i)对threshold进行更新；反之，按照公式threshold＝0.3*(threshold+c(i))更新threshold，同时将数组进行左移操作，丢弃已经移出的数据。对相关信号中所有数据均进行此操作，确定适合背景噪声的判决门限最终值；
(4)遍历相关信号c中的所有数据，找出大于经过(3)所得判决门限的数据并记录此数据在接收信号中对应的位置，此记录结果中包含信号真正的起始位置和疑似信号起始位置；
(5)对(4)中所记录的每个位置i之后的一段时长为tlfm的信号进行分数阶傅里叶变换，找到对应的最佳阶次pi,同时计算此位置之后一段时长为tlfm的信号的调频斜率ki，若同时满足|p0-pi|＜η和
Figure BDA0002234287270000021
其中η为最佳阶次判决阈值，则记录这些位置作为经分数阶傅里叶变换检测所得的同步信号起始位置，若不满足，则将不满足的位置记录为疑似信号起始位置。以此循环直至(4)中所有位置均完成此操作，得到所有满足条件的位置；
(6)将(5)中记录的疑似信号起始位置之后一段时长为tlfm的信号再采用分段傅里叶变换进行二次筛选，再次记录作为经分数阶傅里叶变换检测所得的同步信号起始位置；
(7)综合经过(5)和(6)所得的位置作为信号的起始位置。
优选地，步骤(2)中，设定p的步长为pstep＝0.01，并且以此步长从0开始递增到pstop＝2。
步骤(3)中，K取4。
步骤(5)中，η值范围为小于0.03。
步骤(5)中，分数阶傅里叶变换时，阶次从0开始以0.05为步长逐步搜索到2。
步骤(6)的步骤如下：
a.将时长为tlfm的信号分为10小段，取出每小段做FFT得对应的频率，记录所有位置的所有10段FFT的结果；
b.把本地LFM信号的频率范围fst～fe分成10个小段[fst,fs+1],[fst+1,fst+2],......,[fst+10,fe]，以每段的中心频率
Figure BDA0002234287270000022
作为判决频率，创建判决频率数组，对步骤a中所有位置的10个经FFT所得的频率与判决数组中对应的频率做判断，若10个频率中有超过6个频点落在判决频率数组中对应频率的可接受频率范围内，则记录此位置。
本发明提出的一种基于分段傅里叶变换和分数阶傅里叶变换相结合的检测线性调频信号起始位置的方法，利用线性调频信号在不同阶次的分数阶傅里叶变换出现的能量聚集的现象与虚警信号进行区分，同是在此基础上加入阈值判决与分段傅里叶变换，很好的解决了水声通信中由于虚警信号引起的帧起始判决失误导致的码元判决错误的现象，提高了判决的准确性，提升了系统的通信性能。
附图说明
图1为本发明的具体流程图。
具体实施方式
下面对本发明进行说明。
1.生成相关信号：首先根据公式lfm(t)＝sin(2×π×(fst×t+0.5×kt×t2))生成长度为tlfm的本地信号，其中fst为LFM信号的起始频率，kt为LFM信号的调频斜率。利用Hilbert变换将本地的线性调频样本信号构造为解析信号lfmH(t)，利用公式a(t)＝IFFT(FFT(x(t))×FFT(lfmH(t)))将此解析信号与经过滤波和线性放大接收的数据做相关检测生成相关信号c，其中x(t)为接收端接收到的信号，由于线性调频信号具有良好的自相关特性，即与本身做相关检测后会出现很大的相关峰，而其互相关特性很差，即若与其他信号进行相关操作产生的相关信号幅值明显低于自相关信号幅值；
2.搜索本地信号分数阶傅里变换的最佳阶次p0：对本地信号作M(LFM信号长度*采样频率)点P阶的分数阶傅里叶变换，其中p以pstep＝0.01为步长从0开始逐步搜索到pstop＝2，找出经过各阶次分数阶傅里叶变换之后模值最大对应的阶次p0即为样本信号的最佳阶次。
3.计算自适应门限：初始化用作门限计算的数组a[N]全为0，其中N＝300，相关峰门限threshold＝0，对于经过相关检测后所得的相关信号的每一点判断是否大于数组a中所有数据中位数的4倍，若大于则根据公式threshold＝0.8*threshold+0.2*c(i)更新threshold，其中c(i)为相关后相关信号每点的数值；反之，根据公式threshold＝0.3*(threshold+c(i))更新threshold，对相关信号所有点完成此操做后得到的threshold即为相关信号判决门限阈值最终值，对相关信号的每点进行遍历，找到大于此threshold的值在接收信号中对应的位置并记录。
4.搜索最佳阶次和计算调频斜率：对3中记录的所有位置之后的长度为tlfm的信号作M点的分数阶傅里叶变换，同样设定搜索步长pstep＝0.05，终止阶次pstop＝2确定每个记录位置之后时长为tlfm信号的最佳阶次pi(i＝1,2,3,……)，根据公式
Figure BDA0002234287270000031
计算各个最佳阶次对应的旋转角度，同时按照下式计算每个最佳阶次对应的调频斜率：
Figure BDA0002234287270000032
其中fs为采样频率，M为分数阶傅里叶变换的点数。若同时满足|p0-pi|＜0.05和
Figure BDA0002234287270000033
这里将最佳阶次判决阈值η取为0.05，则将此位置记录为经分数阶傅里叶变换搜索所得的信号起始位置，或不满足则将此位置记录为疑似信号起始位置；
6.结合分段傅里叶变换做二次判决：
将本地LFM信号分成10小段每段频率范围为[fst,fst+1],[fst+1,fst+2],......,[fst+10,fe]，以每段的中心频率
Figure BDA0002234287270000034
作为判决频率，创建判决频率数组T[10]其中的10个数值分别为LFM分段后每段的中心频率，设置计数变量CTC＝0。
(1)提所有经过4之后所记录的疑似信号起始位置之后时长为tlfm的信号，将此信号分为10段，每段长度为
Figure BDA0002234287270000041
(其中[]为取整操作)，M为本地LFM信号的点数。对每段分别进行
Figure BDA0002234287270000042
点的傅里叶变换，记录每个位置的10段傅里叶变换结果。
(2)将(1)中的每个位置的分段傅里叶变换后的结果与T[10]中对应的值做对比，若两者之差落在可接受范围内(可接受范围是指T[10]中对应值上下100Hz)则CTC加1，当CTC大于6时，则认为此位置为信号起始位置记录并保存，一个位置判断完成后，将CTC清零，继续执行上述两步操作直至所有自适应门限值的位置均判决完毕。
7.综合6和5筛选所得的位置作为最后的信号起始位置。

Claims (5)

1.一种基于分数阶傅里叶变换的线性调频信号的检测方法，包括下列的步骤：
(1)通信过程中以收发两端均已知的一段线性调频LFM信号作为同步信号，设此信号频率范围为fst～fe,调频斜率为kt，信号时长为tlfm，接收端首先根据此参数生成LFM信号作为本地信号，然后将接收到的信号与本地信号进行相关检测，将所得相关信号c的数据保存；
(2)对(1)中生成的本地信号进行p阶次的分数阶傅里叶变换，搜索本地信号对应的最佳阶次p0即经分数阶傅里叶变换后最大模值对应的阶次；
(3)接收端生成长度为N的数组a[N]初始化全为0，并设置判决门限的初始数值threshold，将相关信号当前数据c(i)赋给数组中最后一个数据a[N-1]，判断此数据是否大于数组a中所有数据中位数的K＝4倍，若大于则按照公式threshold＝0.8*threshold+0.2*c(i)对threshold进行更新；反之，按照公式threshold＝0.3*(threshold+c(i))更新threshold，同时将数组进行左移操作，丢弃已经移出的数据；对相关信号中所有数据均进行此操作，确定适合背景噪声的判决门限最终值；
(4)遍历相关信号c中的所有数据，找出大于经过(3)所得判决门限的数据并记录此数据在接收信号中对应的位置，此记录结果中包含信号真正的起始位置和疑似信号起始位置；
(5)对(4)中所记录的每个位置i之后的一段时长为tlfm的信号进行分数阶傅里叶变换，找到对应的最佳阶次pi,同时计算此位置之后一段时长为tlfm的信号的调频斜率ki，若同时满足|p0-pi|＜η和
Figure FDA0002992114220000011
其中η为最佳阶次判决阈值，则记录这些位置作为经分数阶傅里叶变换检测所得的同步信号起始位置，若不满足，则将不满足的位置记录为疑似信号起始位置；以此循环直至(4)中所有位置均完成此操作，得到所有满足条件的位置；
(6)将(5)中记录的疑似信号起始位置之后一段时长为tlfm的信号再采用分段傅里叶变换进行二次筛选，再次记录作为经分数阶傅里叶变换检测所得的同步信号起始位置，方法如下：
a.将时长为tlfm的信号分为10小段，取出每小段做FFT得对应的频率，记录所有位置的所有10段FFT的结果；
b.把本地LFM信号的频率范围fst～fe分成10个小段[fst,fs+1],[fst+1,fst+2],......,[fst+10,fe]，以每段的中心频率
Figure FDA0002992114220000012
作为判决频率，创建判决频率数组，对步骤a中所有位置的10个经FFT所得的频率与判决数组中对应的频率做判断，若10个频率中有超过6个频点落在判决频率数组中对应频率的可接受频率范围内，则记录此位置；
(7)综合经过(5)和(6)所得的位置作为信号的起始位置。
2.根据权利要求1所述的方法，其特征在于，步骤(2)中，设定p的步长为pstep＝0.01，并且以此步长从0开始递增到pstop＝2。
3.根据权利要求1所述的方法，其特征在于，步骤(5)中，η值范围为小于0.03。
4.根据权利要求1所述的方法，其特征在于，步骤(5)中，分数阶傅里叶变换时，阶次从0开始以0.05为步长逐步搜索到2。
5.根据权利要求1所述的方法，其特征在于，可接收频率范围是，以判决频率数组中对应频率为中心上下浮动100Hz。
CN201910978041.0A 2019-10-15 2019-10-15 基于分数阶傅里叶变换的线性调频信号的检测方法 Active CN110808791B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910978041.0A CN110808791B (zh) 2019-10-15 2019-10-15 基于分数阶傅里叶变换的线性调频信号的检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910978041.0A CN110808791B (zh) 2019-10-15 2019-10-15 基于分数阶傅里叶变换的线性调频信号的检测方法

Publications (2)

Publication Number Publication Date
CN110808791A CN110808791A (zh) 2020-02-18
CN110808791B true CN110808791B (zh) 2021-05-28

Family

ID=69488463

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910978041.0A Active CN110808791B (zh) 2019-10-15 2019-10-15 基于分数阶傅里叶变换的线性调频信号的检测方法

Country Status (1)

Country Link
CN (1) CN110808791B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111935038B (zh) * 2020-08-03 2022-08-19 中国人民解放军国防科技大学 基于分数阶傅里叶变换的线性调频干扰消除方法
CN112559973B (zh) * 2021-02-25 2021-05-14 中国人民解放军空军预警学院 自适应多分量线性调频信号参数估计方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109510787A (zh) * 2018-10-15 2019-03-22 中国人民解放军战略支援部队信息工程大学 水声信道下线性调频信号参数估计方法及装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19734713A1 (de) * 1997-08-11 1999-02-18 Mikrowellen Technologie Und Se Radar-Entfernungsmeßeinrichtung
CN101094217B (zh) * 2007-07-27 2010-04-21 哈尔滨工程大学 水声信道参数测量方法及移动水声通信同步获取方法
CN101237253B (zh) * 2008-02-29 2011-02-09 哈尔滨工业大学 基于分数阶傅立叶变换的自适应干扰分离的信号收发装置
CN104007318B (zh) * 2014-06-17 2016-11-09 中国科学院电子学研究所 获取信号时频函数的方法
CN106443588B (zh) * 2016-05-23 2019-03-12 中国人民解放军63892部队 一种lfmcw信号快速检测和估计方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109510787A (zh) * 2018-10-15 2019-03-22 中国人民解放军战略支援部队信息工程大学 水声信道下线性调频信号参数估计方法及装置

Also Published As

Publication number Publication date
CN110808791A (zh) 2020-02-18

Similar Documents

Publication Publication Date Title
CN110808791B (zh) 基于分数阶傅里叶变换的线性调频信号的检测方法
CN106098079B (zh) 音频信号的信号提取方法与装置
US20040167775A1 (en) Computational effectiveness enhancement of frequency domain pitch estimators
De Seixas et al. Preprocessing passive sonar signals for neural classification
CN101477801A (zh) 一种检测和消除数字音频信号中脉冲噪声的方法
JP4383949B2 (ja) 同期パルスを生成する方法、装置およびその装置を備える受信機
CN109061591B (zh) 一种基于序贯聚类的时频线谱检测方法
CN112751629A (zh) 基于时频图像处理的宽带特定信号检测方法
CN109991520A (zh) 一种电缆振荡波局放检测系统用波速计算新方法
CN114785379A (zh) 一种水声janus信号参数估计方法及系统
CN112462321B (zh) 基于脉冲匹配累加的声信标信号截获处理方法
CN111413702B (zh) 用于宽带探鱼仪的高效目标分割方法
CN110632563B (zh) 一种基于短时傅里叶变换的脉内频率编码信号参数测量方法
CN110233687B (zh) 一种多路数据信号联合判决检测方法
CN115932808A (zh) 一种基于多特征融合的被动声纳智能探测方法
WO2018026329A1 (en) Pitch period and voiced/unvoiced speech marking method and apparatus
CN112333123B (zh) 非合作psk水声通信信号多小波基自动优选基带解调方法
JP5276427B2 (ja) 受信装置及びシンボルタイミング検出方法
CN110808938B (zh) 基于最佳阶次分数阶傅里叶变换的水声信号同步方法
CN100424692C (zh) 音频快速搜索方法
CN114325671B (zh) 改进的基于水下声信标信号多脉冲累积的抗多普勒方法
CN113054945A (zh) 一种基于熵值分析的声表面波谐振器有效激励检测方法
CN116647251B (zh) 适合mc-cdma系统的信号捕获方法
CN111257859B (zh) 一种波束域自相关水下目标识别方法
CN116055262B (zh) 基于同步挤压小波变换的通信信号载频盲估计方法、系统及介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant