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

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

Detection method of LFM signal based on fractional Fourier transform

【技术实现步骤摘要】
基于分数阶傅里叶变换的线性调频信号的检测方法
本专利技术属于水声通信数字信号处理
，涉及一种同步信号起始位置检测技术。
技术介绍
通信系统中同步方式主要有自同步和外同步两种方式，外同步是水声通信最常用的方式。其做法是在每帧信号的起始位置加入接收双方均已知一段信号。这段信号通常是具有良好自相关特性的线性调频(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|＜η和其中η为最佳阶次判决阈值，则记录这些位置作为经分数阶傅里叶变换检测所得的同步信号起始位置，若不满足，则将不满足的位置记录为疑似信号起始位置。以此循环直至(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]，以每段的中心频率作为判决频率，创建判决频率数组，对步骤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，对相关信号所有点完成此操做后得到的th本文档来自技高网...

【技术保护点】
1.一种基于分数阶傅里叶变换的线性调频信号的检测方法，包括下列的步骤：/n(1)通信过程中以收发两端均已知的一段线性调频LFM信号作为同步信号，设此信号频率范围为f

【技术特征摘要】
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倍，若大于则按照公式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|＜η和其中η为最佳阶次判决阈值，则记录这些位置作为经分数阶傅里叶变换检测所得的同步信号起始位置，若不满足，则将不满足的位...

【专利技术属性】
技术研发人员：苏毅珊，周兆嘉，董立杰，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津;12