CN110703217B - 一种基于自相关坐标轴旋转的线性调频信号检测方法及系统 - Google Patents

一种基于自相关坐标轴旋转的线性调频信号检测方法及系统 Download PDF

Info

Publication number
CN110703217B
CN110703217B CN201910804558.8A CN201910804558A CN110703217B CN 110703217 B CN110703217 B CN 110703217B CN 201910804558 A CN201910804558 A CN 201910804558A CN 110703217 B CN110703217 B CN 110703217B
Authority
CN
China
Prior art keywords
signal data
signal
frequency modulation
rotation angle
detected
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
CN201910804558.8A
Other languages
English (en)
Other versions
CN110703217A (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN201910804558.8A priority Critical patent/CN110703217B/zh
Publication of CN110703217A publication Critical patent/CN110703217A/zh
Application granted granted Critical
Publication of CN110703217B publication Critical patent/CN110703217B/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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418—Theoretical aspects
    • 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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/292—Extracting wanted echo-signals
    • 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
    • G06F17/141—Discrete Fourier transforms
    • G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
    • 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Remote Sensing (AREA)
  • Algebra (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Discrete Mathematics (AREA)
  • Computing Systems (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明属于信号处理领域，具体涉及一种基于自相关坐标轴旋转的线性调频信号检测方法及系统，获取待检测线性调频信号；根据待检测线性调频信号得到第一信号数据；根据第一信号数据得到第二信号数据；根据第二信号数据得到第三信号数据；根据待检测线性调频信号得到最优旋转角；根据最优旋转角得到第四信号数据；根据第四信号数据得到补偿第四信号数据；根据补偿第四信号数据得到第五信号数据；根据第五信号数据得到第五信号数据的峰值坐标；根据待检测线性调频信号得到调频率的估计值和中心频率的估计值。本发明可以在低信噪比的情况下可以有效的实现多分量线性调频信号的检测和参数估计。

Description

一种基于自相关坐标轴旋转的线性调频信号检测方法及系统
技术领域
本发明属于信号处理领域，具体涉及一种基于自相关坐标轴旋转的线性调频信号检测方法及系统。
背景技术
线性调频信号因其优良的特性被广泛的应用于雷达、声纳和通信等工程领域。近年来，线性调频信号越发成为研究的热点，各国学者也相继提出了许多算法来解决线性调频信号的检测和参数估计问题。这些算法大体上可以分为以下两类：搜索类算法和非搜索类算法；典型的搜索类算法包括：分数阶傅里叶变换、Radon-Wigner变换和Radon模糊变换等；这类算法通过搜索参数实现在各自域中积累，进而实现信号的检测和参数估计；分数阶傅里叶变换是传统傅里叶变换的一种广义的变形，它通过搜索与调频率对应的阶数来实现信号的积累和检测；Radon-Wigner变换和Radon模糊变换通过二维搜索来对Wigner-Ville平面和模糊函数平面进行积累，但这种积累是非相参积累，积累增益有限，除此之外，多分量线性调频信号带来的交叉项和伪峰问题会大大影响它们的检测性能。
2011年由Lv提出的一种积累算法LVD是一种典型的非搜索类算法，它重构了时频平面，消除了时间和延时变量在原时频平面中的耦合，进而进行新时频平面的二维相参积累，完成检测和参数估计。因为其积累过程是相参的，因此相比于Radon-Wigner变换和Radon模糊变换，LVD有着更好的积累性能。但是，由于该方法重构时频平面时引起了信噪比的损失，其在低信噪比情况下的检测性能有限。
发明内容
为了解决现有技术中存在的上述问题，本发明提供了一种基于自相关坐标轴旋转的线性调频信号检测方法及系统。本发明要解决的技术问题通过以下技术方案实现：
一种基于自相关坐标轴旋转的线性调频信号检测方法，获取待检测线性调频信号；
对所述待检测线性调频信号进行数字采样和预处理，得到第一信号数据；
对所述第一信号数据进行自相关处理，得到第二信号数据；
对所述第二信号数据进行快速傅里叶变换和平移操作，得到第三信号数据；
根据所述待检测线性调频信号得到最优旋转角；
根据所述最优旋转角对所述第三信号数据进行坐标轴旋转，得到第四信号数据；
对所述第四信号数据进行相位补偿，得到补偿第四信号数据；并对所述补偿第四信号数据进行快速傅里叶变换，得到第五信号数据；
根据所述第五信号数据得到所述第五信号数据的峰值坐标；
根据所述最优旋转角和所述第五信号数据的峰值坐标对所述待检测线性调频信号的调频率和中心频率进行估计，得到所述调频率的估计值和所述中心频率的估计值。
在本发明的一个实施例中，根据所述待检测线性调频信号得到最优旋转角，包括：
根据所述待检测线性调频信号得到待搜索旋转角；
根据所述待搜索旋转角对所述第三信号数据进行坐标轴旋转变换得到旋转后的第三信号数据；
对所述旋转后的第三信号数据进行相位补偿快速傅里叶变换，得到傅里叶变换后的第三信号数据；
根据所述傅里叶变换后的第三信号数据得到最大的傅里叶变换后的第三信号数据，并根据所述最大的傅里叶变换后的第三信号数据得到最优旋转角。
在本发明的一个实施例中，所述第一信号数据为长度为P的向量，P＝Lfs，L和fs均大于0，其中fs为采样率，L为信号时长。
在本发明的一个实施例中，所述第二信号数据为一个P行P列的矩阵，表达式为：S2(nt,nτ)＝exp[j2π(2fcen(nτ-floor(0.5P)-1)T+2K(nt-1)(nτ-floor(0.5P)-1)T2)]，
其中，nt为离散的时间变量，nτ为离散的延时变量，fcen为待检测线性调频信号的中心频率，floor表示向下取整操作，K为调频率，T为采样间隔。
在本发明的一个实施例中，所述第三信号数据为一个P行P列的矩阵，表达式为：
Figure GDA0003142047000000041
其中，kt为离散时间频率变量，FFTshift(·)表示平移操作，FFT(·)表示快速傅里叶变换。
在本发明的一个实施例中，所述第四信号数据为一个P行P列的矩阵，记为S4(kt',nτ')，第四信号数据和第三信号数据的对应关系为：
Figure GDA0003142047000000042
其中，kt'为坐标轴旋转后的离散时间频率变量，nτ'为坐标轴旋转后的离散延时变量，round(·)为四舍五入操作，αmax为最优旋转角。
在本发明的一个实施例中，所述傅里叶变换后的第三信号数据表达式如下：
Figure GDA0003142047000000043
其中，kτ'为离散延时频率变量，
Figure GDA0003142047000000044
为使用第i个待搜索旋转角旋转后的第三个信号数据，αi为第i个待搜索旋转角。
在本发明的一个实施例中，所述第五信号数据的表达式为：
Figure GDA0003142047000000045
其中，A表示第五信号数据的幅度，fL＝2fcen+KTP。
在本发明还提供了一种基于自相关坐标轴旋转的线性调频信号检测系统，包括：
信号获取模块，用于获取待检测线性调频信号；
信号处理模块，用于对所述待检测线性调频信号进行数字采样和预处理，得到第一信号数据；用于对所述第一信号数据进行自相关处理，得到第二信号数据；用于对所述第二信号数据进行快速傅里叶变换和平移操作，得到第三信号数据；用于根据所述待检测线性调频信号得到最优旋转角；用于根据所述最优旋转角对所述第三信号数据进行坐标轴旋转，得到第四信号数据；用于对所述第四信号数据进行相位补偿，得到补偿第四信号数据；还用于对所述补偿第四信号数据进行快速傅里叶变换，得到第五信号数据；
信号检测模块，用于根据所述第五信号数据得到所述第五信号数据的峰值坐标；还用于根据所述最优旋转角和所述第五信号数据的峰值坐标对所述待检测线性调频信号的调频率和中心频率进行估计，得到所述调频率的估计值和所述中心频率的估计值。
本发明的有益效果：
本发明通过对第四信号数据的特定行进行相位补偿和快速傅里叶变换可以有效的实现待检测目标的相参积累，在低信噪比的情况下可以有效的实现多分量线性调频信号的检测和参数估计。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测方法的流程框图；
图2是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测方法中仿真实验分量1的检测结果图；
图3是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测方法中仿真实验分量2的检测结果图；
图4是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测方法实测雷达回波数据示意图；
图5是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测方法对实测数据的检测结果；
图6是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测系统的结构框图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述，但本发明的实施方式不限于此。
请参见图1，图1是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测方法的流程框图，获取待检测线性调频信号；
对所述待检测线性调频信号进行数字采样和预处理，得到第一信号数据；
对所述第一信号数据进行自相关处理，得到第二信号数据；
对所述第二信号数据进行快速傅里叶变换和平移操作，得到第三信号数据；
根据所述待检测线性调频信号得到最优旋转角；
根据所述最优旋转角对所述第三信号数据进行坐标轴旋转，得到第四信号数据；
对所述第四信号数据进行相位补偿，得到补偿第四信号数据；并对所述补偿第四信号数据进行快速傅里叶变换，得到第五信号数据；
根据所述第五信号数据得到所述第五信号数据的峰值坐标；
根据所述最优旋转角和所述第五信号数据的峰值坐标对所述待检测线性调频信号的调频率和中心频率进行估计，得到所述调频率的估计值和所述中心频率的估计值。
本发明通过对第四信号数据的特定行进行相位补偿和快速傅里叶变换可以有效的实现待检测目标的相参积累，在低信噪比的情况下可以有效的实现多分量线性调频信号的检测和参数估计。
在本发明的一个实施例中，根据所述待检测线性调频信号得到最优旋转角，包括：
根据所述待检测线性调频信号得到待搜索旋转角；
根据所述待搜索旋转角对所述第三信号数据进行坐标轴旋转变换得到旋转后的第三信号数据；
对所述旋转后的第三信号数据进行相位补偿快速傅里叶变换，得到傅里叶变换后的第三信号数据；
根据所述傅里叶变换后的第三信号数据得到最大的傅里叶变换后的第三信号数据，并根据所述最大的傅里叶变换后的第三信号数据得到最优旋转角。
进一步地，获取最优旋转角的方法具体如下：
(1)根据待检测线性调频信号的调频率范围[Kmin,Kmax]，确定N个待搜索旋转角{α1,α2,...,αN}，其中α1＝arctan(2T2P·Kmin)为第一个待搜索旋转角，αN＝arctan(2T2P·Kmax)为最后一个待搜索旋转角，αw＝α1+(w-1)(αN-α1)/N为第w个待搜索旋转角；
令旋转次数i＝1；
(2)使用第i个待搜索旋转角αi对所述第三信号数据进行坐标轴旋转变换，得到旋转后的第三信号数据；
(3)取出旋转后的第三信号数据的第floor(0.5P)+1行数据进行相位补偿快速傅里叶变换，得到傅里叶变换后的第三信号数据，并从傅里叶变换后的第三信号数据中找出最大数据值记为A(i)；
(4)令循环次数i＝i+1，重复执行上述步骤(2)(3)直到i＝N，得到N个最大数据值{A(1),A(2),…,A(N)}，确定{A(1),A(2),…,A(N)}中的最大值Amax(I)，即最大的傅里叶变换后的第三信号数据，进而确定最大的傅里叶变换后的第三信号数据Amax(I)在{A(1),A(2),…,A(N)}中的索引I，根据该索引I确定该最大值Amax(I)对应的待搜索旋转角αI，I∈{1,2,…,N}；
将待搜索旋转角αI确定为最优旋转角αmax。
在本发明的一个实施例中，所述第一信号数据为长度为P的向量，P＝Lfs，L和fs均大于0，其中fs为采样率，L为信号时长。
在本发明的一个实施例中，所述第二信号数据为一个P行P列的矩阵，表达式为：S2(nt,nτ)＝exp[j2π(2fcen(nτ-floor(0.5P)-1)T+2K(nt-1)(nτ-floor(0.5P)-1)T2)]，
其中，nt为离散的时间变量，nτ为离散的延时变量，fcen为待检测线性调频信号的中心频率，floor表示向下取整操作，K为调频率，T为采样间隔。
在本发明的一个实施例中，所述第三信号数据为一个P行P列的矩阵，表达式为：
Figure GDA0003142047000000091
其中，kt为离散时间频率变量，FFTshift(·)表示平移操作，FFT(·)表示快速傅里叶变换。
在本发明的一个实施例中，所述第四信号数据为一个P行P列的矩阵，记为S4(kt',nτ')，第四信号数据和第三信号数据的对应关系为：
Figure GDA0003142047000000092
其中，kt'为坐标轴旋转后的离散时间频率变量，nτ'为坐标轴旋转后的离散延时变量，round(·)为四舍五入操作，αmax为最优旋转角。
在本发明的一个实施例中，所述傅里叶变换后的第三信号数据表达式如下：
Figure GDA0003142047000000093
其中，kτ'为离散延时频率变量，
Figure GDA0003142047000000094
为使用第i个待搜索旋转角旋转后的第三个信号数据，αi为第i个待搜索旋转角。
在本发明的一个实施例中，所述第五信号数据的表达式为：
Figure GDA0003142047000000101
其中，A表示第五信号数据的幅度，fL＝2fcen+KTP。
具体的，调频率的估计值为Kest，中心频率的估计值为fest，进而得到：
Kest＝tanαmax/(2T2P)，
Figure GDA0003142047000000102
其中，k为第五信号数据的峰值坐标。
以下通过仿真实验对本发明效果作进一步验证说明：
(一)仿真实验参数设置；
本实验设置接收信号数据参数如下：信号采样率fs＝256Hz，采样点数256个，分量1：中心频率f1＝20Hz，调频率K1＝10Hz/s，信噪比SNR1＝-5dB；分量2：中心频率f2＝40Hz，调频率K2＝30Hz/s,信噪比SNR1＝-7dB。
(二)仿真实验内容及结果分析；
本发明实施例对2分量线性调频信号进行检测和参数估计。请参见图2和图3，图2是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测方法中仿真实验分量1的检测结果图，图3是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测方法中仿真实验分量2的检测结果图。
可见，通过本发明实施例的处理，分量1和分量2可以被明显的区分和检测。同时，根据分量1和分量2的峰值坐标和最优旋转角，可以对它们的参数实现准确的估计。
(三)实测数据条件；
实测数据为雷达接收到的数据，所述雷达的载频为8.85GHz，带宽为40MHz，脉冲重复频率为1000Hz。
(四)实测数据处理内容及结果；
采用本发明实施例对所述雷达回波数据进行处理，并分析检测结果。请参见图4，图4是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测方法实测雷达回波数据示意图；图5为本发明实施例对线性调频信号的检测结果。
选取所述雷达回波第54个距离单元进行处理，该段数据中包含一个线性调频信号。由图5可见，通过本发明实施例的处理，该线性调频信号得到有效检测，证明了本发明有效性。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
请参见图6，图6是本发明实施例提供的一种基于自相关坐标轴旋转的线性调频信号检测系统的结构框图，包括：
信号获取模块，用于获取待检测线性调频信号；
信号处理模块，用于对所述待检测线性调频信号进行数字采样和预处理，得到第一信号数据；用于对所述第一信号数据进行自相关处理，得到第二信号数据；用于对所述第二信号数据进行快速傅里叶变换和平移操作，得到第三信号数据；用于根据所述待检测线性调频信号得到最优旋转角；用于根据所述最优旋转角对所述第三信号数据进行坐标轴旋转，得到第四信号数据；用于对所述第四信号数据进行相位补偿，得到补偿第四信号数据；还用于对所述补偿第四信号数据进行快速傅里叶变换，得到第五信号数据；
信号检测模块，用于根据所述第五信号数据得到所述第五信号数据的峰值坐标；还用于根据所述最优旋转角和所述第五信号数据的峰值坐标对所述待检测线性调频信号的调频率和中心频率进行估计，得到所述调频率的估计值和所述中心频率的估计值。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (8)

1.一种基于自相关坐标轴旋转的线性调频信号检测方法，其特征在于，包括：
获取待检测线性调频信号；
对所述待检测线性调频信号进行数字采样和预处理，得到第一信号数据；
对所述第一信号数据进行自相关处理，得到第二信号数据；
对所述第二信号数据进行快速傅里叶变换和平移操作，得到第三信号数据；
根据所述待检测线性调频信号得到最优旋转角；
根据所述最优旋转角对所述第三信号数据进行坐标轴旋转，得到第四信号数据；
对所述第四信号数据进行相位补偿，得到补偿第四信号数据；并对所述补偿第四信号数据进行快速傅里叶变换，得到第五信号数据；
根据所述第五信号数据得到所述第五信号数据的峰值坐标；
根据所述最优旋转角和所述第五信号数据的峰值坐标对所述待检测线性调频信号的调频率和中心频率进行估计，得到所述调频率的估计值和所述中心频率的估计值；
所述第三信号数据为一个P行P列的矩阵，表达式为：
Figure FDA0003097382500000021
其中，kt为离散时间频率变量，FFTshift(·)表示平移操作，FFT(·)表示快速傅里叶变换，其中，nt为离散的时间变量，nτ为离散的延时变量。
2.根据权利要求1所述的基于自相关坐标轴旋转的线性调频信号检测方法，其特征在于，根据所述待检测线性调频信号得到最优旋转角，包括：
根据所述待检测线性调频信号得到待搜索旋转角；
根据所述待搜索旋转角对所述第三信号数据进行坐标轴旋转变换得到旋转后的第三信号数据；
对所述旋转后的第三信号数据进行相位补偿和快速傅里叶变换，得到傅里叶变换后的第三信号数据；
根据所述傅里叶变换后的第三信号数据得到最大的傅里叶变换后的第三信号数据，并根据所述最大的傅里叶变换后的第三信号数据得到最优旋转角。
3.根据权利要求1所述的基于自相关坐标轴旋转的线性调频信号检测方法，其特征在于，所述第一信号数据为长度为P的向量，P＝Lfs，L和fs均大于0，其中fs为采样率，L为信号时长。
4.根据权利要求3所述的基于自相关坐标轴旋转的线性调频信号检测方法，其特征在于，所述第二信号数据为一个P行P列的矩阵，表达式为：S2(nt,nτ)＝exp[j2π(2fcen(nτ-floor(0.5P)-1)T+2K(nt-1)(nτ-floor(0.5P)-1)T2)]，
其中，nt为离散的时间变量，nτ为离散的延时变量，fcen为待检测线性调频信号的中心频率，floor表示向下取整操作，K为调频率，T为采样间隔。
5.根据权利要求4所述的基于自相关坐标轴旋转的线性调频信号检测方法，其特征在于，所述第四信号数据为一个P行P列的矩阵，记为S4(kt',nτ')，第四信号数据和第三信号数据的对应关系为：
Figure FDA0003097382500000031
其中，kt'为坐标轴旋转后的离散时间频率变量，nτ'为坐标轴旋转后的离散延时变量，round(·)为四舍五入操作，αmax为最优旋转角。
6.根据权利要求5所述的基于自相关坐标轴旋转的线性调频信号检测方法，其特征在于，所述傅里叶变换后的第三信号数据表达式如下：
Figure FDA0003097382500000032
其中，kτ'为离散延时频率变量，
Figure FDA0003097382500000033
为使用第i个待搜索旋转角旋转后的第三个信号数据，αi为第i个待搜索旋转角。
7.根据权利要求6所述的基于自相关坐标轴旋转的线性调频信号检测方法，其特征在于，所述第五信号数据的表达式为：
Figure FDA0003097382500000034
其中，A表示第五信号数据的幅度，fL＝2fcen+KTP。
8.一种基于自相关坐标轴旋转的线性调频信号检测系统，其特征在于，包括：
信号获取模块，用于获取待检测线性调频信号；
信号处理模块，用于对所述待检测线性调频信号进行数字采样和预处理，得到第一信号数据；用于对所述第一信号数据进行自相关处理，得到第二信号数据；用于对所述第二信号数据进行快速傅里叶变换和平移操作，得到第三信号数据；用于根据所述待检测线性调频信号得到最优旋转角；用于根据所述最优旋转角对所述第三信号数据进行坐标轴旋转，得到第四信号数据；用于对所述第四信号数据进行相位补偿，得到补偿第四信号数据；还用于对所述补偿第四信号数据进行快速傅里叶变换，得到第五信号数据；
信号检测模块，用于根据所述第五信号数据得到所述第五信号数据的峰值坐标；还用于根据所述最优旋转角和所述第五信号数据的峰值坐标对所述待检测线性调频信号的调频率和中心频率进行估计，得到所述调频率的估计值和所述中心频率的估计值；
所述第三信号数据为一个P行P列的矩阵，表达式为：
Figure FDA0003097382500000041
其中，kt为离散时间频率变量，FFTshift(·)表示平移操作，FFT(·)表示快速傅里叶变换，其中，nt为离散的时间变量，nτ为离散的延时变量。
CN201910804558.8A 2019-08-28 2019-08-28 一种基于自相关坐标轴旋转的线性调频信号检测方法及系统 Active CN110703217B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910804558.8A CN110703217B (zh) 2019-08-28 2019-08-28 一种基于自相关坐标轴旋转的线性调频信号检测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910804558.8A CN110703217B (zh) 2019-08-28 2019-08-28 一种基于自相关坐标轴旋转的线性调频信号检测方法及系统

Publications (2)

Publication Number Publication Date
CN110703217A CN110703217A (zh) 2020-01-17
CN110703217B true CN110703217B (zh) 2021-08-24

Family

ID=69194068

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910804558.8A Active CN110703217B (zh) 2019-08-28 2019-08-28 一种基于自相关坐标轴旋转的线性调频信号检测方法及系统

Country Status (1)

Country Link
CN (1) CN110703217B (zh)

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101094217B (zh) * 2007-07-27 2010-04-21 哈尔滨工程大学 水声信道参数测量方法及移动水声通信同步获取方法
CN103237000B (zh) * 2013-04-22 2016-01-20 北京理工大学 Frft-ofdm系统的低复杂度峰均比抑制方法
CN104007318B (zh) * 2014-06-17 2016-11-09 中国科学院电子学研究所 获取信号时频函数的方法
CN105783974B (zh) * 2014-12-25 2018-08-07 中国科学院声学研究所 一种线性调频信号的检测、参数估计方法及系统
US10180495B2 (en) * 2016-07-05 2019-01-15 The Aerospace Corporation Separating weak and strong moving targets using the fractional fourier transform
CN106338721B (zh) * 2016-08-23 2019-03-29 西安电子科技大学 基于多帧回波相参积累的空中匀速弱小目标检测方法
JP2018119858A (ja) * 2017-01-25 2018-08-02 日本信号株式会社 Ｆｍｃｗ方式レーダ
CN106896358A (zh) * 2017-04-27 2017-06-27 电子科技大学 一种基于位置旋转变换的高速目标相参积累检测方法
CN107167777B (zh) * 2017-06-20 2019-10-18 南京理工大学 锯齿波线性调频信号参数提取方法
CN110146886A (zh) * 2019-05-29 2019-08-20 西安电子科技大学 非均匀旋转目标运动参数的快速估计方法

Also Published As

Publication number Publication date
CN110703217A (zh) 2020-01-17

Similar Documents

Publication Publication Date Title
CN107561508B (zh) 一种用于匀加速运动目标的相参积累检测方法
CN107290589B (zh) 基于短时分数阶傅里叶变换的非线性信号时频分析方法
CN107153189B (zh) 线性调频连续波雷达测距的信号处理方法
CN108196241B (zh) 一种基于Hough变换的高速动目标速度估计方法
Catherall et al. High resolution spectrograms using a component optimized short-term fractional Fourier transform
CN107329115B (zh) 基于gcrbf网络的lfm信号参数估计方法
Liu et al. Frequency estimation of chirp signals based on fractional fourier transform combined with Otsu’s method
Khan et al. Robust sparse reconstruction of signals with gapped missing samples from multi-sensor recordings
US6577968B2 (en) Method of estimating signal frequency
CN110703217B (zh) 一种基于自相关坐标轴旋转的线性调频信号检测方法及系统
CN108957389A (zh) 一种实数域多通道信号目标方位估计方法
CN117368875A (zh) 一种宽带雷达多模态运动目标相参积累检测方法及装置
TWI693419B (zh) 訊號處理方法
CN112255609A (zh) 一种恒加速机动目标相参积累检测方法
CN105652256B (zh) 一种基于极化信息的高频地波雷达tbd方法
CN111044996A (zh) 一种基于降维近似消息传递的lfmcw雷达目标检测方法
CN112130138A (zh) 一种超视距雷达的电离层相位污染校正方法
CN107450058B (zh) 基于FrFT和HT的雷达信号时频参数估计方法
Sun et al. A novel multicomponent LFM signal detection based on fast block FRFT
CN112505640B (zh) 基于参数自适应的扩展b分布脉冲信号时频分析方法
CN109597042B (zh) 基于奇异谱分析的目标进动频率估计方法
Huang et al. Low-observable maneuvering target detection based on Radon-advanced discrete chirp Fourier transform
US6683820B1 (en) Method and apparatus for tracking sonar targets
CN110764062A (zh) 基于分数阶傅里叶域滤波的多分量线性调频信号参数估计方法
CN112327260A (zh) 一种sar回波数据中脉冲式干扰信号的抑制方法和装置

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