CN111935038A - 基于分数阶傅里叶变换的线性调频干扰消除方法 - Google Patents

基于分数阶傅里叶变换的线性调频干扰消除方法 Download PDF

Info

Publication number
CN111935038A
CN111935038A CN202010768705.3A CN202010768705A CN111935038A CN 111935038 A CN111935038 A CN 111935038A CN 202010768705 A CN202010768705 A CN 202010768705A CN 111935038 A CN111935038 A CN 111935038A
Authority
CN
China
Prior art keywords
interference
fourier transform
received signal
signal
frequency
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.)
Granted
Application number
CN202010768705.3A
Other languages
English (en)
Other versions
CN111935038B (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202010768705.3A priority Critical patent/CN111935038B/zh
Publication of CN111935038A publication Critical patent/CN111935038A/zh
Application granted granted Critical
Publication of CN111935038B publication Critical patent/CN111935038B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03159—Arrangements for removing intersymbol interference operating in the frequency domain
    • 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
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0024—Carrier regulation at the receiver end
    • H04L2027/0026—Correction of carrier offset

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Noise Elimination (AREA)

Abstract

本发明公开了一种基于分数阶傅里叶变换的线性调频干扰消除方法。包括以下步骤：步骤1：对所接收的信号进行量纲归一化；步骤2：估计LFM干扰的调频斜率；步骤3：估计各LFM干扰的初始频率；步骤4：根据所估计得到的各LFM干扰的调频斜率和初始频率，运用子空间投影消除干扰。本发明通过发现在实际采样中，由于信号离散化，为使待处理的信号始终在运算窗口之内，导致在t0时刻，信号x(t)的瞬时频率已不再是初始频率，造成所估计的信号初始频率有偏差，本发明在发现偏差量存在的情况后，给出了一个初始频率的估算方法，使得初始频率的估计更加精确，从而使得在多干扰的环境下，运用子空间投影可以更精确地消除干扰。

Description

基于分数阶傅里叶变换的线性调频干扰消除方法
技术领域
本发明主要涉及到信号处理领域，尤其涉及一种基于分数阶傅里叶变换的线性调频干扰消除方法。
背景技术
战场环境下，接收机可能接收到多个平稳以及非平稳干扰，其中线性调频干扰LFM(Linear Frequency Modulation)是一种普遍的非平稳干扰形式，对于LFM干扰，运用分数阶傅里叶变换进行处理可以得到较好的效果。文献1“王鹏,邱天爽,李景春等.基于高斯加权分数阶傅里叶变换的LFM信号参数估计[J].通信学报,2016,37(4):107-115.”提出了一种高斯加权分数阶傅里叶变换方法，应用于低占空比时的LFM信号参数估计，可以有效提高此种情况下的参数估计精度。文献2“张玉恒,吴启晖,王金龙.基于时频加窗短时傅里叶变换的LFM干扰抑制[J].电子与信息学报,2007,29(6):1361-1364.”提出了一种基于时频加窗短时傅里叶变换(TFW-STFT)的LFM干扰抑制算法，所提出的时频窗对LFM干扰具有较好的频域能量聚集性能，因此加窗对信号的影响要小于无聚集性能的短时傅里叶变换，该算法在信噪比损失和系统误比特率上明显优于基于短时傅里叶变换的算法。但是在实际应用过程中发现，在估计LFM信号的初始频率时，往往会出现偏差。
运用分数阶傅里叶变换可以较准确地估计LFM干扰的参数且运算复杂度比Wigner变换低，在众多文献中得到应用，但是很少有研究频率估计偏差的文献。偏差与误差的性质不同，是由LFM干扰的调频斜率决定的偏差量，不能忽略不计。
分数阶傅里叶变换的定义式如下。设接收信号为x(t)，信号x(t)阶数为p的分数阶傅里叶变换FRFT(Fractional Fourier Transform)定义为
Figure BDA0002615659590000011
式中Kp(t,u)为分数阶傅里叶变换的核函数，其中p为FRFT的的阶次。Kp(t,u)的定义如下：
Figure BDA0002615659590000021
式中n为整数，α表示旋转角度，α与阶次P的关系为α＝Pπ/2。上式的意义是：1、对于一般情况，也即α非π的整数倍时，核函数采取(2)式第一行的形式；2、对于特殊情况，也即α为π的整数倍时，核函数采取冲激函数的形式，即如(2)式第二、三行所示。
以上为分数阶傅里叶变换的定义，其运算复杂难以实现，采用文献3即Ozaktas提出的采样型离散分数阶傅里叶变换(DFRFT:Discrete Fractional Fourier Transform)算法“Haldun M.Ozaktas.Digital computation of the fractional Fourier transform,IEEE Trans on Signal Processing,44(9):2141～2150”，可以化繁为简。
发明内容
本发明要解决的技术问题是提高LFM干扰频率参数估算精度，从而使在此基础上构建的干扰子空间提高消除干扰的精度，提出了一种基于分数阶傅里叶变换的线性调频干扰消除方法。
为解决该问题，所采用的技术方案是：
步骤1：对所接收的信号进行量纲归一化；
步骤2：估计LFM干扰的调频斜率；
步骤2.1：将所述量纲归一化后的接收信号序列通过分数阶傅里叶变换映射到以变换域变量u为横轴，阶数p为纵轴的平面上；
所述分数阶傅里叶变换为：
步骤2.1.1：取接收信号x(t)，以1/Δy为采样间隔进行采样，设长度为L，得到接收信号数据序列
Figure BDA0002615659590000022
步骤2.1.2：对接收信号数据序列
Figure BDA0002615659590000023
的前后补零，补零的长度为数据序列
Figure BDA0002615659590000024
的前后各补
Figure BDA0002615659590000025
长度，其中
Figure BDA0002615659590000026
运算表示向上取整；
步骤2.1.3：取切普信号exp[-jπt2cot(α)]，以1/(2Δy)为采样间隔进行采样，其长度与
Figure BDA0002615659590000031
补零后数据序列相等。
步骤2.1.4：因对于x(t)与exp[-jπt2cot(α)]的采样间隔不同，后者采样间隔是前者一半，所以需要对x(t)采样之后的数据序列进行内插，使用香农内插公式对x(t)采样之后的数据序列
Figure BDA0002615659590000032
进行内插，然后与exp[-jπt2cot(α)]采样之后的结果逐点相乘，不妨设其结果为g(n)。
Figure BDA0002615659590000033
步骤2.1.5：对信号g(n)进行卷积运算；
将信号g(n)与因子
Figure BDA0002615659590000034
相乘后做卷积，将卷积运算求和上下限的中点定到信号g(n)的中点，卷积结果记为
Figure BDA0002615659590000035
步骤2.1.6：将
Figure BDA0002615659590000036
与Aαexp(-jπu2cot(α))的离散表示形式
Figure BDA0002615659590000037
相乘就得到了FRFT表达式；
其中Aα为常数，表示如下：
Figure BDA0002615659590000038
步骤2.2：使分数阶p在区间(0,2)之间递增，每递增一步进行一次分数阶傅里叶变换，得到接收信号在分数阶变换后的变阶分数阶傅里叶变换图谱；
步骤2.3：在变换图谱上搜索谱峰，谱峰位置的纵坐标对应LFM干扰的调频斜率，设某个谱峰所在点纵坐标为p1，则对应的旋转角度为α1＝p1·π/2，调频斜率即为c＝-cotα1；
步骤2.4：通过步骤2.3找到所有谱峰，得到所有LFM干扰的调频斜率；
步骤3：估计各LFM干扰的初始频率；
设某个LFM干扰对应谱峰的横坐标为u1，则其初始频率为
Figure BDA0002615659590000041
上式等号右边没有未知量，其中T取归一化之后的值1；
步骤4：根据所估计得到的各LFM干扰的调频斜率和初始频率，运用子空间投影消除干扰。
进一步地，步骤2.2中，0.5≤|p|≤1.5，当p∈(0,0.5)时，首先对接收信号进行傅里叶反变换，然后对结果运用适用于0.5≤|p|≤1.5时的FRFT处理；当p∈(1.5,2)时，首先对接收信号进行傅里叶变换，然后对变换后的结果运用适用于0.5≤|p|≤1.5时的FRFT处理。
进一步地，运用子空间投影消除干扰的方法是：
步骤5.1：根据估计的LFM干扰参数构建干扰信号序列；
用修正后的频率参数构建的信号为xi(n)＝expjπ(cin2+2fin)，i＝1,2,…,K，其中K表示干扰个数，ci、fi分别表示第i个LFM干扰的调频斜率和初始频率，如ci＝0，则该干扰为单频干扰；
取xi(n)序列的L长度取一段数据作为子空间基矢量：
Figure BDA0002615659590000042
步骤5.2：构建干扰子空间及其正交子空间。
用时域扩展后的干扰信号矢量构建成一个矩阵，
B＝[X1 X2 … XP]
该组矢量张成了干扰的子空间，其正交子空间表示如下:
Figure BDA0002615659590000043
I为单位矩阵；
步骤5.3：进行子空间投影消除干扰；
将接收信号矢量投影到干扰子空间的正交子空间，如下式
Y＝PX＝PXs+PV
式中X为接收信号矢量，通过投影过程消除了接收信号中的干扰分量，剩余有用信号和白噪声分量Xs和V。
本发明还提供了一种基于分数阶傅里叶变换的线性调频干扰消除系统，包括处理器，以及与所述处理器连接的存储器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时实现上述基于分数阶傅里叶变换的线性调频干扰消除方法的步骤。
本发明还提供了一种计算机可读介质，其上存储计算机程序，所述计算机程序可被处理器执行以实现上述基于分数阶傅里叶变换的线性调频干扰消除方法的步骤。
与现有技术相比，本发明所取得的有益效果是：
本发明基于分数阶傅里叶变换的线性调频干扰消除方法，在多干扰环境下，通过对信号进行分数阶傅立叶变换，在分数阶傅立叶变换的阶数对应的旋转角度递增时，得到变阶分数阶傅里叶变换图谱，通过扫描图谱中的谱峰，根据谱峰位置估计得到各干扰的调频斜率和初始频率。然后根据估计得到的干扰参数，消除干扰。本发明通过理论及实际采样的运算过程中分析发现，由于实际采样时，信号离散化，为使待处理的信号始终在运算窗口之内，需要在信号的前后补零，将导致卷积结果的前面一段接近于零。假设卷积结果在t0之前接近于零，而实际的信号流在零时刻已经开始了。而在t0时刻，信号x(t)的瞬时频率已不再是初始频率，在分数阶傅里叶变换的运算中，x(t)在此时的频率将被当作初始频率，因此所估计的信号的初始频率有偏差，不是真实的初始频率，从而带来运用子空间投影消除干扰的不准确，本发明在发现偏差量存在的情况后，给出了一个初始频率的估算方法，通过实验数据表明，在考虑了偏移量后，初始频率的估计更加精确，从而使得在多干扰的环境下，运用子空间投影可以更精确地消除干扰。
附图说明
图1为本发明系统流程图；
图2为变阶分数阶傅里叶变换，2(a)为阶数值变动时的分数阶傅里叶变换平面图，2(b)阶数值变动时的分数阶傅里叶变换立体图；
图3为基于分数阶傅里叶变换频率参数估计的LFM、单频干扰抑制效果。
具体实施方式
图1至图3示出了本发明基于分数阶傅里叶变换的线性调频干扰消除方法的一种具体实施例，如图1所示，包括以下步骤：
步骤1：对所接收的信号进行量纲归一化；
采用文献3中，运用Ozaktas的方法进行处理的过程包含了对时间和频率的采样，将两者量纲归一化便于分析与处理。即对于信号的时间t，有t∈[0,T]，因非对称区间难于计算，需将信号左移半个作用区间，即t∈[-T/2,T/2]。设频率范围为f∈[-F/2,F/2]。根据文献4“杨清红.线性调频信号参数的快速估计技术研究[D].华南理工大学,2015”的方法，通过引入一个尺度因子S来实现量纲的归一化。设尺度变换之后的新坐标为：
y＝t/S,v＝F*S
尺度变化之后，时域和频域表达式将被限定在新时宽T/S和新带宽F*S，现令尺度因子
Figure BDA0002615659590000061
则新的时宽和带宽大小相等，且均等于Δy。故此，原来的时域和频域区间均变成无量纲区间[-Δy/2,Δy/2]，在此基础上可以对归一化的信号以采样间隔1/Δy进行采样。
步骤2：估计LFM干扰的调频斜率；
步骤2.1：将所述量纲归一化后的接收信号离散化，对离散后的接收信号序列通过分数阶傅里叶变换映射到以变换域变量u为横轴，阶数P为纵轴的平面上；
所述分数阶傅里叶变换为：
由背景技术中的式2可知，当α为π的整数倍时，分数阶傅里叶变换的核函数转化成为冲激函数，由此可以将定义式简化，无需讨论，所以讨论一般情况，即α取值范围为(0,π)∪(π,2π)，也即p的取值范围为(0,2)∪(2,4)时。在此情况下，(2)式表示的核函数只需要取第一项，在此情况下信号x(t)的p阶分数阶傅里叶变换为
Figure BDA0002615659590000071
式中Aα为常数，表示如下
Figure BDA0002615659590000072
本实施例中，对于实际运算过程中因采取离散化采样计算，具体实施过程为：
步骤2.1.1：取接收信号x(t)，以1/Δy为采样间隔进行采样，设长度为L，得到离散化后的接收信号数据序列
Figure BDA0002615659590000073
步骤2.1.2：对接收信号数据序列
Figure BDA0002615659590000074
的前后补零，补零的长度为数据序列
Figure BDA0002615659590000075
的前后各补
Figure BDA0002615659590000076
长度，其中
Figure BDA0002615659590000077
运算表示向上取整；
由于实际采样时，信号离散化，为使待处理的信号始终在运算窗口之内，需要在信号的前后补零，将导致卷积结果的前面一段接近于零。假设卷积结果在t0之前接近于零，而实际的信号流在零时刻已经开始了。而在t0时刻，信号x(t)的瞬时频率已不再是初始频率，在分数阶傅里叶变换的运算中，x(t)在此时的频率将被当作初始频率，因此所估计的信号的初始频率有偏差，不是真实的初始频率，从而带来运用子空间投影消除干扰的不准确。
步骤2.1.3：取切普信号exp[-jπt2cot(α)]，以1/(2Δy)为采样间隔进行采样，其长度与
Figure BDA0002615659590000078
补零后数据序列相等；
步骤2.1.4：因对于x(t)与exp[-jπt2cot(α)]的采样间隔不同，后者采样间隔是前者一半，对x(t)采样之后的数据序列进行内插，使用香农内插公式对x(t)采样之后的数据序列
Figure BDA0002615659590000081
进行内插，然后与exp[-jπt2cot(α)]采样之后的结果逐点相乘得到
Figure BDA0002615659590000087
对于接收信号x(t)，将接收信号x(t)与一个初始频率为0，调频斜率为t2cot(α)的线性调频信号exp(jπt2cot(α))相乘，如下式
g(t)＝exp(jπt2cot(α))x(t) (6)
为利于离散化计算g(t)，需在频域上将exp(jπt2cot(α))x(t)限制在无量纲带宽F内，旋转角度需限制在区间0.5≤|P|≤1.5内，此时(6)式可以用香农内插公式表示为：
Figure BDA0002615659590000082
根据文献3的分析，因子exp(jπt2cot(α))的带宽时频积为因子x(t)带宽时频积的两倍，为使之匹配，需要做相应设置。在转换到新的坐标之后，时间和频率共同的归一化采样区间宽度为Δy。如果x(t)的采样间隔为1/Δy，则g(t)的采样间隔应为
Figure BDA0002615659590000083
但如果按照该方式采样，则时间点不匹配。可按如下方式设置：以
Figure BDA0002615659590000084
为采样间隔对g(t)进行采样，仍以1/Δy为间隔对x(t)进行采样，但香农内插公式对其进行两倍插值，使两者时间点匹配，正如(5)式所示。
步骤2.1.5：对信号序列g(n)进行卷积运算，将信号g(n)与因子
Figure BDA0002615659590000085
相乘后做卷积，将卷积运算求和上下限的中点定到信号g(n)的中点，卷积结果记为
Figure BDA0002615659590000086
m为整数；
步骤2.1.6：将
Figure BDA0002615659590000091
与Aαexp(-jπu2cot(α))的离散表示形式
Figure BDA0002615659590000092
相乘就得到分数阶傅里叶变换表达式；
Figure BDA0002615659590000093
其中Aα为常数，表示如下：
Figure BDA0002615659590000094
α表示旋转角度，α与阶次p的关系为α＝pπ/2；
步骤2.1.5和步骤2.1.6也可以通过直接对连续信号g(t)进行计算，然后再进行离散化，具体如下：
对连续信号g(t)进行积分计算，与exp(-j2πutcscα)相乘之后求积分，积分的结果为自变量u的函数，不妨设该函数为g′(u)，
Figure BDA0002615659590000095
按照文献4“杨清红.线性调频信号参数的快速估计技术研究[D].华南理工大学,2015”的运算思路，这一积分类似于傅里叶变换，可推导运算如下：
Figure BDA0002615659590000101
式中rect(·)表示矩形函数，不妨设s为该函数自变量，当矩形函数自变量s满足-1≤s≤1时，有rect(s)≡1。而当0.5≤|P|≤1.5时，因变换域变量u满足|u|≤Δy/2，此时矩形函数的自变量在±1之间，恒等于1。则(9)式可简化为
Figure BDA0002615659590000102
将g′(u)与切普信号Aαexp(-jπu2cot(α))相乘就得到了FRFT表达式；
其中Aα为常数，表示如下：
Figure BDA0002615659590000103
α表示旋转角度，α与阶次p的关系为α＝pπ/2；
将g′(u)与切普信号Aαexp(-jπu2cot(α))相乘，如下式
Figure BDA0002615659590000104
在上式基础上，在u域对其进行采样，使之离散化，可以得到
Figure BDA0002615659590000111
运用半角公式之中的tan(α/2)＝cscα-cotα，可将上式转化为如下形式
Figure BDA0002615659590000112
令
Figure BDA0002615659590000113
则上式求和号之内部分可以看作h(n)与f(n)的卷积，它可以用傅里叶变换来实现。
步骤2.2：使分数阶p在区间(0,2)之间递增，每递增一步进行一次分数阶傅里叶变换，得到接收信号在分数阶变换后的变阶分数阶傅里叶变换图谱；
本实施例中，线性扫频信号的时频脊线在时频分布图上是关于其中点对称的，绕其中点翻转180度形状不变，所以α取值为(0,π)时已经遍历了所有情况，对应的，分数阶p取值区间即(0,2)。
由步骤2.1看出，在做分数阶傅里叶变换时，阶次p限制在区间0.5≤|p|≤1.5内，对于阶次未涵盖到的开区间(0,0.5)和开区间(1.5,2)，不妨设阶次p1∈(0,0.5)，阶次p2∈(1.5,2)，则有
Figure BDA0002615659590000114
式中Fp表示阶次为p的分数阶傅里叶变换，而根据分数阶傅里叶变换的性质，
Figure BDA0002615659590000115
和
Figure BDA0002615659590000116
分别是求取傅里叶变换和傅里叶反变换运算。即当p∈(0,0.5)时，首先对接收信号进行傅里叶反变换，然后对结果运用适用于0.5≤|p|≤1.5时的FRFT处理；当p∈(1.5,2)时，首先对接收信号进行傅里叶变换，然后同样对结果运用适用于0.5≤|p|≤1.5时的FRFT处理。
步骤2.3：在变换图谱上搜索谱峰，谱峰位置的纵坐标对应LFM干扰的调频斜率，设某个谱峰所在点纵坐标为p1，则对应的旋转角度为α1＝p1·π/2，调频斜率即为c＝-cotα1；
一般情况下，p是未知的，可利用分数阶傅里叶变换的性质查找出p来。将p在区间(0,2)之间递增，在每一步进行一次分数阶傅里叶变换，当步进值递增到准确p值时，分数阶傅里叶变换会出现尖峰。
设所需要检测的LFM干扰信号表示式为
Figure BDA0002615659590000121
根据上文所设，p逐渐由0增加到2，在此过程中，LFM干扰对应的p值将会被找到。如图2(a)和2(b)所示阶数值变动时的分数阶傅里叶变换平面图和立体图。在此情况下对LFM干扰做分数阶傅里叶变换，有下式
Figure BDA0002615659590000122
由上式可知，当且仅当c＝-cotα时，积分式(10)内t的二次项可以被消除。因为根据(6)式的统计特性，只有当二次项被消除时，接收信号x(t)的p阶傅里叶变换才会出现谱峰。也即，当分数阶傅里叶变换的阶数对应的旋转角度递增，递增到与LFM干扰在时频面上的脊线垂直时，二次项将会被消除，此时LFM干扰的调频斜率c＝-cotα。
判断p值是否为1，如p值为1，分数阶傅里叶变换转换为傅里叶变换。将接收信号中的单频干扰映射为变阶分数阶傅里叶变换图谱中的谱峰，谱峰对应的横坐标即单频干扰的归一化频率值。
步骤2.4：通过步骤2.3遍历图谱中所有谱峰，得到所有LFM干扰的调频斜率；
步骤3：估计各LFM干扰的初始频率；
设某个LFM干扰对应谱峰的横坐标为u1，谱峰的纵坐标所对应阶数转换的旋转角为α1时，则其初始频率为
Figure BDA0002615659590000131
上式等号右边没有未知量，其中T取归一化之后的值1；
本实施例中，在估测了LFM干扰的调频斜率之后，将c＝-cotα代入(16)式，可以得到
Figure BDA0002615659590000132
但以上分析是在理论情况下，此时积分上下限为无穷，实际中积分区间是有线长度，同前文，设其区间为[0,T]，则有
Figure BDA0002615659590000133
如前文所述，非对称区间难于计算，需将待处理的信号左移半个作用区间，即t∈[-T/2,T/2]。信号左移半个区间之后，以LFM干扰信号而言，相当于表示式发生了变化，应表述如下
Figure BDA0002615659590000134
即调频斜率未发生变化，仍然是
Figure BDA0002615659590000135
但初始频率变为
Figure BDA0002615659590000136
那么对此时的信号x′(t)作阶数为p的分数阶傅里叶变换，为简便起见，将系数Aα省去不影响分析结果，有
Figure BDA0002615659590000137
由上式可以推出，当
Figure BDA0002615659590000138
时，LFM干扰的分数阶傅里叶变换的模可以取到峰值。因此，所估计的初始频率，
Figure BDA0002615659590000141
需要指出的是：单频干扰可以看作特殊的LFM干扰，如接收信号中含有单频干扰，则在FRFT的阶数p变化到”1”时，即可出现谱峰，因为此时的分数阶傅里叶变换已经转化为一般的傅里叶变换。
步骤4：根据所估计得到的各LFM干扰的调频斜率和初始频率，运用子空间投影消除干扰。
运用子空间投影消除干扰的基本思想是：设想接收信号中含有的干扰构成一个子空间，找出该子空间最直接的方法就是找出子空间的基矢量。由于各干扰互不相关，所以构建的干扰序列可作为干扰子空间的基矢量，每一个干扰信号对应一个基矢量，以此可以构建干扰子空间。在构建了干扰子空间的基础上，也就可以运用相应的子空间运算达到消除干扰的目的。
步骤4.1：根据估计的LFM干扰参数构建干扰信号序列；
根据步骤2和步骤3所估计的LFM干扰的频率参数之后，构建干扰信号序列
xi(n)＝exp(j2π(fin+ci/2n2)),i＝1,2,…,K (21)式中K表示干扰个数，fi、ci分别表示第i个干扰的初始频率和调频斜率。取xi(n)序列的L′长度数据作为子空间基矢量，以此构建干扰子空间，在此基础上运用子空间投影消除干扰。
取xi(n)序列的L长度段数据作为子空间基矢量：
Figure BDA0002615659590000142
LFM干扰初始频率的估计结果存在误差，误差尽管微小，但也对子空间投影结果有影响。根据文献5”Zhai Yong-Ping,Zhou Zhu.Linear Frequency ModulationInterference Suppression Based on Error Correction.2019 IEEE9th InternationalConference on Electronics Information and Emergency Communication(ICEIEC),Beijing,China,2019,pp.1-5.”的分析：同样的误差对于子空间投影结果的影响与构建子空间基矢量的长度有关，在同等误差情况下，用长序列作为干扰子空间基矢量引起的干扰抑制效果下降明显，构成干扰子空间基矢量的序列过短也存在频谱泄露的问题。
据上可知，选择合适长度的序列构成干扰子空间基矢量，可以达成良好的抗干扰效果。因为所有的设计均需要经过大量测试实验，大量实验也可以大致确定参数估计的误差范围，也就容易找出最佳的构建子空间的序列长度。这一长度值不需要确切数字。
步骤4.2：构建干扰子空间及其正交子空间。
用时域扩展后的干扰信号矢量构建成一个矩阵，
B＝[X1 X2 … XP] (23)
该组矢量张成了干扰的子空间，其正交子空间表示如下
Figure BDA0002615659590000151
I为单位矩阵；
步骤4.3：进行子空间投影消除干扰；
将接收信号矢量投影到干扰子空间的正交子空间，如下式
Y＝PX＝PXs+PV (25)
式中X为接收信号矢量，通过投影过程消除了接收信号中的干扰分量，剩余有用信号和白噪声分量Xs和V。
抗干扰处理效果分析：
采用场景一，测试本发明的抗干扰效果。将本发明所估计的LFM干扰频率参数用于构建干扰子空间，所估计的LFM干扰的频率参数为：f1＝0.0999、c1＝0.2；f2＝0.2001、c2＝0.2。而单频干扰的频率值估计很精确，两个单频干扰的归一化频率值为：-0.1、0.4。
现对比两种情况：一、用本发明方法估计的频率参数构建干扰子空间，在此基础上通过子空间投影消除干扰；二、用精确的频率参数构建干扰子空间，在此基础上运用子空间投影消除干扰。其中第一种情况即是对本研究最终效果的测试，第二种情况提供对比。因子空间构建与基矢量长度有直接关系，故以子空间长度为横坐标，输出信干噪比为纵坐标，结果如图3所示。
如图3，图中点划线、星划线分别对应于频率参数取精确值、估计值两种情况。由图3可知：1、当用于构建干扰子空间的频率参数取精确值时，序列长度增加对于输出SINR没有影响，始终保持很高的输出SINR；2、用本发明方法估计干扰的频率参数，并据此构建干扰子空间，在此情况下，当用于构建干扰子空间序列较短时，能够获得较高的输出SINR，随着序列长度增加，SINR将会下降，但是序列长度取到32，SINR仍能保持较高，有用信号很容易从噪声中恢复。
由此可知：用本发明方法可以很好的消除单天线接收的数个LFM干扰以及单频干扰。
一种基于分数阶傅里叶变换的线性调频干扰消除系统，包括处理器，以及与所述处理器连接的存储器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时实现上述基于分数阶傅里叶变换的线性调频干扰消除方法的步骤。
一种计算机可读介质，其上存储计算机程序，所述计算机程序可被处理器执行以实现上述基于分数阶傅里叶变换的线性调频干扰消除方法的步骤。
以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。

Claims (5)

1.一种基于分数阶傅里叶变换的线性调频干扰消除方法，其特征在于，包括以下步骤：
步骤1：对所接收的信号进行量纲归一化；
步骤2：估计线性调频干扰的调频斜率；
步骤2.1：将所述量纲归一化后的接收信号离散化，对离散后的接收信号序列通过分数阶傅里叶变换映射到以变换域变量u为横轴，阶数p为纵轴的平面上；
所述分数阶傅里叶变换为：
步骤2.1.1：取接收信号x(t)，以1/Δy为采样间隔进行采样，设长度为L，得到离散化后的接收信号数据序列
Figure FDA0002615659580000011
步骤2.1.2：对接收信号数据序列
Figure FDA0002615659580000012
的前后补零，补零的长度为数据序列
Figure FDA0002615659580000013
的前后各补
Figure FDA0002615659580000014
长度，其中
Figure FDA0002615659580000015
运算表示向上取整；
步骤2.1.3：取切普信号exp[-jπt2cot(α)]，以1/(2Δy)为采样间隔进行采样，其长度与
Figure FDA0002615659580000016
补零后数据序列相等；
步骤2.1.4：因对于x(t)与exp[-jπt2cot(α)]的采样间隔不同，后者采样间隔是前者一半，对x(t)采样之后的数据序列进行内插，使用香农内插公式对x(t)采样之后的数据序列
Figure FDA0002615659580000017
进行内插，然后与exp[-jπt2cot(α)]采样之后的结果逐点相乘，得到
Figure FDA0002615659580000018
步骤2.1.5：对信号g(n)进行卷积运算；将信号g(n)与因子
Figure FDA0002615659580000019
相乘后做卷积，将卷积运算求和上下限的中点定到信号g(n)的中点，卷积结果记为
Figure FDA00026156595800000110
m为整数；
步骤2.1.6：将
Figure FDA00026156595800000111
与Aαexp(-jπu2cot(α))的离散表示形式
Figure FDA0002615659580000021
相乘就得到分数阶傅里叶变换表达式；
Figure FDA0002615659580000022
其中Aα为常数，表示如下：
Figure FDA0002615659580000023
α表示旋转角度，α与阶次p的关系为α＝pπ/2；
步骤2.2：使分数阶p在区间(0,2)之间递增，每递增一步进行一次分数阶傅里叶变换，得到接收信号在分数阶变换后的变阶分数阶傅里叶变换图谱；
步骤2.3：在变换图谱上搜索谱峰，谱峰位置的纵坐标对应LFM干扰的调频斜率，设某个谱峰所在点纵坐标为p1，则对应的旋转角度为α1＝p1·π/2，调频斜率即为c＝-cotα1；
步骤2.4：通过步骤2.3找到所有谱峰，得到所有LFM干扰的调频斜率；
步骤3：估计各LFM干扰的初始频率；
设某个LFM干扰对应谱峰的横坐标为u1，谱峰的纵坐标所对应阶数转换的旋转角为α1时，则其初始频率为
Figure FDA0002615659580000024
其中T取归一化之后的值1；
步骤4：根据所估计得到的各LFM干扰的调频斜率和初始频率，运用子空间投影消除干扰。
2.根据权利要求1所述的方法，其特征在于：步骤2.1中的分数阶傅里叶变换，其阶数p限制在0.5≤|p|≤1.5，当p∈(0,0.5)时，首先对接收信号进行傅里叶反变换，然后对结果运用适用于0.5≤|p|≤1.5时的FRFT处理；当p∈(1.5,2)时，首先对接收信号进行傅里叶变换，然后对变换后的结果运用适用于0.5≤|p|≤1.5时的FRFT处理。
3.根据权利要求1所述的方法，其特征在于：运用子空间投影消除干扰的方法是：
步骤4.1：根据估计的LFM干扰参数构建干扰信号序列；
用修正后的频率参数构建的信号为xi(n)＝expjπ(cin2+2fin)，i＝1,2,…,K，其中K表示干扰个数，ci、fi分别表示第i个LFM干扰的调频斜率和初始频率，如ci＝0，则该干扰为单频干扰；
取xi(n)序列的L长度段数据作为子空间基矢量：
Figure FDA0002615659580000031
步骤4.2：构建干扰子空间及其正交子空间。
用时域扩展后的干扰信号矢量构建成一个矩阵，
B＝[X1 X2 … XP]
该组矢量张成了干扰的子空间，其正交子空间表示如下
Figure FDA0002615659580000032
I为单位矩阵；
步骤4.3：进行子空间投影消除干扰；
将接收信号矢量投影到干扰子空间的正交子空间，如下式
Y＝PX＝PXs+PV
式中X为接收信号矢量，通过投影过程消除了接收信号中的干扰分量，剩余有用信号Xs和白噪声分量V。
4.一种基于分数阶傅里叶变换的线性调频干扰消除系统，其特征在于：包括处理器，以及与所述处理器连接的存储器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时实现权利要求1至3任一项所述方法的步骤。
5.一种计算机可读介质，其特征在于：其上存储计算机程序，所述计算机程序可被处理器执行以实现权利要求1至3任一项所述方法的步骤。
CN202010768705.3A 2020-08-03 2020-08-03 基于分数阶傅里叶变换的线性调频干扰消除方法 Active CN111935038B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010768705.3A CN111935038B (zh) 2020-08-03 2020-08-03 基于分数阶傅里叶变换的线性调频干扰消除方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010768705.3A CN111935038B (zh) 2020-08-03 2020-08-03 基于分数阶傅里叶变换的线性调频干扰消除方法

Publications (2)

Publication Number Publication Date
CN111935038A true CN111935038A (zh) 2020-11-13
CN111935038B CN111935038B (zh) 2022-08-19

Family

ID=73307250

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010768705.3A Active CN111935038B (zh) 2020-08-03 2020-08-03 基于分数阶傅里叶变换的线性调频干扰消除方法

Country Status (1)

Country Link
CN (1) CN111935038B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113406672A (zh) * 2021-06-17 2021-09-17 哈尔滨工程大学 一种基于自适应加窗的分段线性调频干扰消除方法
CN114866159A (zh) * 2022-04-01 2022-08-05 华南理工大学 一种多分量线性调频信号时频分析方法
CN115037329A (zh) * 2022-05-31 2022-09-09 江苏屹信航天科技有限公司 抗多普勒的接收方法、接收装置、终端

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8407020B1 (en) * 2010-09-28 2013-03-26 The United States Of America As Represented By The Secretary Of The Navy Fast method to search for linear frequency-modulated signals
CN102999473A (zh) * 2012-10-18 2013-03-27 中国人民解放军电子工程学院 一种线性调频信号的检测与参数估计方法
CN103063909A (zh) * 2012-12-18 2013-04-24 天津理工大学 一种基于功率谱的线性调频信号参数估值方法
WO2014048193A1 (zh) * 2012-09-28 2014-04-03 北京理工大学 一种用于舰艇编队情况下同型雷达同频干扰抑制方法
CN109510787A (zh) * 2018-10-15 2019-03-22 中国人民解放军战略支援部队信息工程大学 水声信道下线性调频信号参数估计方法及装置
CN110133598A (zh) * 2019-05-09 2019-08-16 西安电子科技大学 基于FrFT的线性调频信号参数快速估计方法
CN110764067A (zh) * 2019-10-15 2020-02-07 哈尔滨工程大学 傅里叶变换的lfm-bpsk复合调制雷达信号参数估计方法
CN110808791A (zh) * 2019-10-15 2020-02-18 天津大学 基于分数阶傅里叶变换的线性调频信号的检测方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8407020B1 (en) * 2010-09-28 2013-03-26 The United States Of America As Represented By The Secretary Of The Navy Fast method to search for linear frequency-modulated signals
WO2014048193A1 (zh) * 2012-09-28 2014-04-03 北京理工大学 一种用于舰艇编队情况下同型雷达同频干扰抑制方法
CN102999473A (zh) * 2012-10-18 2013-03-27 中国人民解放军电子工程学院 一种线性调频信号的检测与参数估计方法
CN103063909A (zh) * 2012-12-18 2013-04-24 天津理工大学 一种基于功率谱的线性调频信号参数估值方法
CN109510787A (zh) * 2018-10-15 2019-03-22 中国人民解放军战略支援部队信息工程大学 水声信道下线性调频信号参数估计方法及装置
CN110133598A (zh) * 2019-05-09 2019-08-16 西安电子科技大学 基于FrFT的线性调频信号参数快速估计方法
CN110764067A (zh) * 2019-10-15 2020-02-07 哈尔滨工程大学 傅里叶变换的lfm-bpsk复合调制雷达信号参数估计方法
CN110808791A (zh) * 2019-10-15 2020-02-18 天津大学 基于分数阶傅里叶变换的线性调频信号的检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ZHAI YANG-PING等: "Linear Frequency Modulation Interference Suppression Based on Subspace Error Correction", 《2019 IEEE 9TH INTERNATIONAL CONFERENCE ON ELECTRONICS INFORMATION AND EMERGENCY COMMUNICATION 》 *
ZHIYUAN REN等: "LFM signal sampling and reconstruction based on compressed sensing and swept-frequency theory", 《6TH INTERNATIONAL CONFERENCE ON WIRELESS, MOBILE AND MULTI-MEDIA》 *
孙富礼等: "低信噪比下LFM-BPSK复合调制信号参数估计", 《制导与引信》 *
渠莹: "基于分数阶傅里叶变换的LFM信号参数估计", 《物联网技术》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113406672A (zh) * 2021-06-17 2021-09-17 哈尔滨工程大学 一种基于自适应加窗的分段线性调频干扰消除方法
CN114866159A (zh) * 2022-04-01 2022-08-05 华南理工大学 一种多分量线性调频信号时频分析方法
CN115037329A (zh) * 2022-05-31 2022-09-09 江苏屹信航天科技有限公司 抗多普勒的接收方法、接收装置、终端
CN115037329B (zh) * 2022-05-31 2024-06-11 江苏屹信航天科技有限公司 抗多普勒的接收方法、接收装置、终端

Also Published As

Publication number Publication date
CN111935038B (zh) 2022-08-19

Similar Documents

Publication Publication Date Title
CN111935038B (zh) 基于分数阶傅里叶变换的线性调频干扰消除方法
CN103630878B (zh) 一种基于迭代自适应方法的频带融合方法及装置
CN107085140B (zh) 基于改进的SmartDFT算法的非平衡系统频率估计方法
Ran et al. An adaptive fast factorized back-projection algorithm with integrated target detection technique for high-resolution and high-squint spotlight SAR imagery
CN110113279B (zh) 一种移动跳频水声通信多普勒因子估计方法
CN107315714B (zh) 一种去卷积功率谱估计方法
Xu et al. An accurate sparse SAR imaging method for enhancing region-based features via nonconvex and TV regularization
JP2014153161A (ja) 画像レーダ装置
Xu et al. Ultrawide-band radar imagery from multiple incoherent frequency subband measurements
CN108205080A (zh) 相干平均法谐波信号功率谱估计方法
Zeng et al. Deconvolution of sparse underwater acoustic multipath channel with a large time-delay spread
CN113640790A (zh) 基于二维自适应谱估计的宽带宽脉冲高速目标检测方法
CN113204005B (zh) 一种提高调频连续波激光雷达距离解算精度的方法及装置
JP2013152112A (ja) 時間差方位探知装置
CN114066749B (zh) 相位相关抗噪位移估计方法、设备及存储介质
Qi et al. Chirp Rates Estimation for Multiple LFM Signals by DPT–SVD
Wang et al. Fast and adaptive method for SAR superresolution imaging based on point scattering model and optimal basis selection
CN112578359A (zh) 一种通过双谱变换处理提取雷达信号脉内特征参数的方法
CN111913162A (zh) 射频干扰检测方法及装置、射频干扰抑制方法及装置
CN111880154A (zh) 基于对称波数谱对消的复图像域运动目标检测方法
CN111060883A (zh) 一种无人机干涉sar时变基线误差估计和补偿方法
Wang et al. Instantaneous frequency estimation of nonlinear FM radar signal based on multi-scale Chirplet Path
Wu et al. Motion compensation for OFDM inverse synthetic aperture sonar imaging based on compressed sensing
CN115356729B (zh) 一种近场非均匀采样的直接稀疏成像方法
US20220034994A1 (en) Time difference of arrival estimator based on a joint-optimization formulation

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