CN114095102A - 一种基于时频分析的无人机遥控信号识别方法 - Google Patents
一种基于时频分析的无人机遥控信号识别方法 Download PDFInfo
- Publication number
- CN114095102A CN114095102A CN202111372078.2A CN202111372078A CN114095102A CN 114095102 A CN114095102 A CN 114095102A CN 202111372078 A CN202111372078 A CN 202111372078A CN 114095102 A CN114095102 A CN 114095102A
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- remote control
- aerial vehicle
- unmanned aerial
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 46
- 238000001228 spectrum Methods 0.000 claims abstract description 25
- 238000011084 recovery Methods 0.000 claims abstract description 10
- 230000009466 transformation Effects 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 6
- 230000005540 biological transmission Effects 0.000 claims description 19
- 238000013507 mapping Methods 0.000 claims description 15
- 238000005457 optimization Methods 0.000 claims description 15
- 238000000605 extraction Methods 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 6
- 238000005065 mining Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000007596 consolidation process Methods 0.000 claims description 3
- 238000013480 data collection Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000012806 monitoring device Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 12
- 238000001514 detection method Methods 0.000 description 4
- 241000592183 Eidolon Species 0.000 description 3
- 230000009471 action Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000021615 conjugation Effects 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Selective Calling Equipment (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明属于无人机遥控信号识别技术领域,具体公开了一种基于时频分析的无人机遥控信号识别方法,该方法针对无人机飞行状态时信息传输过程,其中无人机向地面遥控站传输图传信号,遥控站向无人机发射遥控信号,监测装置收集遥控信号和图传信号。包括数据采集及合并、对合并后的时域信号做乔‑威廉姆斯变换得到CWD时频分布、采用基于时频谱能量统计的方法,剔除随机噪声、利用低秩矩阵恢复进行无人机遥控信号提取、采用基于模式距离度量的视频序列相似性探索方法,进行无人机遥控信号识别。采用本发明的方案,摆脱了传统的傅里叶变换单一频域描述信号特征的局限性。
Description
技术领域
本发明属于无人机遥控信号识别技术领域,具体涉及了一种基于时频分析的无人机遥控信号识别方法。
背景技术
近年来伴随社会经济与无线网络的发展,移动无线设备的数量在城市中不断增加,导致无人机所处的环境变得日益复杂,为了提高对低空飞行目标的探测、识别、跟踪和处置能力,各界都在积极开展针对“低慢小”飞行器预警探测、处置拦截技术和设备的研究。对“黑飞”无人机进行捕获方法的研究,如在无人机起飞行时对遥控或图传信号侦收,并实现对无人机的定位,便可以对无人机实现干扰,阻止其在不安全区域飞行,保护低空空域的安全。显然,研究设计出在复杂密集的通信以及无人机涉及的频段信号环境中识别和发现无人机遥控信号,对无人机管控设备进行引导和干扰具有现实意义。
申请号为201910930192.9的专利采用基于频域特征来识别无人机遥控信号,由于无人机遥控信号受到传播损失、周围建筑物反射、噪声和无人机图传信号等的干扰,通常接收到的信号信噪比较低,单一的时域或者频域分析难以体现信号的全部信息。
通常无人机的遥控信号是跳频信号,属于非平稳信号,信号的频率随着时间呈现非线性变化,传统的傅里叶变换不再适用于跳频信号的分析。因此,有必要从时域和频域多角度联合分析,祈望全面展现遥控信号频率随时间变化情况。
发明内容
针对现有技术中所存在的不足,提供一种基于时频分析的无人机遥控信号识别方法,基于时频分析利用时间频率域来表征并且进行处理和分析信号的方法,摆脱了传统的傅里叶变换单一频域描述信号特征的局限性。
为实现上述目的,本发明采用了如下的技术方案:
一种基于时频分析的无人机遥控信号识别方法,该方法针对无人机飞行状态时信息传输过程,其中无人机向地面遥控站传输图传信号,遥控站向无人机发射遥控信号,监测装置收集遥控信号和图传信号。
该无人机遥控信号识别方法包含:
S100:数据采集及合并
采用基于零中频接收机架构的数据采集系统,对无人机信号进行采集。接收机内部的模拟/数据器件对接收的模拟信号进行模数转换,并通过正交解调、数字下变频将接收信号转换为同相分量I和正交分量Q两路数字基带信号,然后将I、Q两路信号进行合并,即时域信号Sig=I+j·Q,j表示虚数单位,
S200:对合并后的时域信号做乔-威廉姆斯(Choi-Williams)变换得到CWD时频分布
采用乔-威廉姆斯变换对合并后的时域信号进行以时间为横轴,频率为纵轴进行分析,通过对采集的Num采数据点进行处理,每组数据为Num组点,共分析Num采/Num组次,每次代表一个时间点。Num组点对应Bw带宽内的信号频率,每一点频率为Bw/Num组。
CWD时频分布为:
式中,μ为时间积分参数,τ为时延积分参数;β为衰减系数,β越小,对交叉项的抑制能力越大,但β的取值要在分辨率和交叉项抑制之间取折中,一般0.1<β<10;t为工作时间,f为频率,Sig为时域信号上角标#表示求共轭。
以上共经过Num采/Num组次乔-威廉姆斯变换处理,每一次输出皆可得到CWD时频谱;再计算每组数据沿着时间轴的CWD时频谱最大值,将各组数据的最大值拼接起来,得到全部样本数据的CWD时频谱最大值图。
S300:剔除随机噪声
在实际采集的无人机信号中会存在许多随机噪声,进行时频分析时,根据随机噪声的幅度通常弱于无人机信号的特点,通过设定能量门限,来消除随机噪声对遥控信号检测的影响。
假设阈值为CWD阈值(t,f),则
阈值设定采用基于时频谱能量统计的方法,即对时频谱进行能量累计,设置合适的阈值,可实现阈值的自适应效果;
式中,Nt为时域样点数,Nf为频域样点数。
S400:无人机遥控信号提取
低秩矩阵恢复先将时频谱矩阵CWD(t,f)表示为低秩矩阵A与稀疏(噪声)矩阵E之和,再通过求解范数优化问题来恢复低秩矩阵。
即求解下列最优化问题:
利用迭代阈值算法(iterative thresholding,IT)求解(1)式得到稀疏矩阵E。
进一步,求解方法如下:
a)将最优化问题式(1)正则化,便得到优化问题:
b)建立优化问题式(2)的拉格朗日(Lagrange)函数:
式中,μ为平衡参数,这里取μ=0.05,Y为拉格朗日乘子矩阵,< >表示求内积。
c)迭代交替更新矩阵A,E和Y
当E=Ek,Y=Yk时,
当A=Ak+1,Y=Yk时,
当A=Ak+1,E=Ek+1时,
Yk+1=Yk+δk(CWD(t,f)-Ak+1-Ek+1)
其中δk为迭代步长,且0<δk<1,k为迭代次数。通过控制迭代次数,使其达到预期结果。或者当Ek+1-Ek≈0时,结束迭代。得到稀疏矩阵E,或称为采样数据时频序列。
S500:无人机遥控信号识别
经过低秩矩阵恢复方法可以进一步抑制图传信号干扰,在稀疏矩阵E中得到清晰的遥控信号时频表示,从而提取出遥控信号的时频特征,进一步与无人机遥控(跳频)信号(图样)特征库中的图样进行比对。
采用基于模式距离度量的时频序列相似性搜索方法,对提取的时频谱特征进行识别,得到该无人机的型号。
在步骤S500中,无人机遥控信号识别,包含:
S501:以时频序列分段表示的每段直线与时间轴的夹角θ作为度量:
将时频序列的模式定义为五元集,记为模式集合M={-2,-1,0,1,2}。
于是,某型无人机遥控信号标准时频序列为:
S1={s1,1,s1,2,…s1,i…s1,N}
式中,s1,i=(m1,i,t1,i),i(i=1,2,…,N)为时频序列序号,N表示时频序列模式的划分数量,s1,i表示时频序列S1中的第i个模式。m1,i为模式集合中的一个元素,这里m1,i∈{-2,-1,0,1,2};
S502:沿时间轴滑动抽取时频序列
从采集数据时频序列E中取与S1同样长度的时频序列S2,假设S2={s2,1,s2,2,…s2,i…s2,N}
式中,s2,i=(m2,i,t2,i),s2,i表示时频序列S2中的第i个模式,且假设t1,i=t2,i,i=1,2,…,N。
采集数据时频序列抽取时,沿时间轴从开始位置取与S1同样长度的时频序列,依次滑动一个时间段再取与S1同样长度的时频序列,直至取完所有的数据。
S503:估计时频序列的模式距离,并判别相似程度
时频序列的模式距离表示具有相同长度两个序列趋势的差异程度,则标准时频序列S1与采集数据时频序列S2的模式距离为:
其中,模式距离DM(s1,i,s2,i)=|m1,i-m2,i|,i=1,2,…,N,| |表示求绝对值;tN表示序列长度,tih表示第i个模式的保持时间,δ为相似度门限值,这里取δ=0.3。
此时采集数据时频序列S2与标准时频序列S1相似,可以判断该机型。
相比于现有技术,本发明具有如下有益效果:
1、本发明的方法基于时频分析利用时间-频率域来表征并且进行处理和分析无人机遥控信号的方法,摆脱了传统的傅里叶变换单一频域描述信号特征的局限性。时频域的分析比单一时域的信号冲击响应或单一频域的信号功率谱密度更为有效。通常把时域和频域结合起来描述和观察信号的时频特征,构成信号的时频谱图,通过时频分析算法分析信号的频谱含量怎样随时间变化。
2、本发明基于低秩矩阵恢复时频滤波方法,采用乔-威廉姆斯变换将无人机信号变换到时频域上,然后依据遥控信号与图传信号在时频域上能量分布的相关性的差异,通过低秩矩阵恢复方法将遥控信号与图传信号分到稀疏矩阵和低秩矩阵中,从而分离遥控信号与图传信号,降低图传信号对遥控信号的干扰。
3、本发明时频序列的模式距离具有模式保持时间加权因子,只有在模式切换点的时间噪声才可能改变时频序列的模式距离,且时间噪声引起的模式保持时间改变量远小于整个序列的时间长度,它对整个距离的影响很小,使得它对时间噪声具有较好的免疫力。
附图说明
图1为本发明的基于时频分析的无人机遥控信号识别方法的流程图。
图2为本发明基于模式距离度量的时频序列相似性搜索方法的流程图。
图3为本发明无人机飞行状态时信息传输的示意图。
图4为本发明无人机采集信号的乔-威廉姆斯时频谱图。
图5为本发明无人机采集信号的乔-威廉姆斯时频谱在时频平面的等高线图。
图6为本发明采样数据时频序列的时频图。
图7为本发明无人机遥控信号标准库中大疆精灵4型的示意图。
图8为本发明时频序列抽取及滑动搜索示意图。
具体实施方式
实施例基本如图1所示,一种基于时频分析的无人机遥控信号识别方法,下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种基于时频分析的无人机遥控信号识别方法,如图3示,为本发明无人机飞行状态时信息传输的示意图,其中无人机向地面遥控站传输图传信号,遥控站向无人机发射遥控信号,监测装置收集遥控信号和图传信号。
如图1所示,为本发明的基于时频分析的无人机遥控信号识别方法的流程图,该识别方法包含以下步骤:
S100:数据采集及合并
采用基于AD9361零中频接收机架构的数据采集系统,对大疆某型号无人机信号进行采集。接收机内部的模拟/数据器件对接收的模拟信号进行模数转换,并通过正交解调、数字下变频将接收信号转换为同相分量I和正交分量Q两路数字基带信号,然后将I、Q两路信号进行合并,得到时域信号Sig=I+j·Q,j表示虚数单位,
S200:对合并后的时域信号做乔-威廉姆斯(Choi-Williams)变换得到CWD时频分布
由于无人机遥控信号属于非平稳信号,因此采用乔-威廉姆斯变换对合并后的时域信号进行以时间为横轴,频率为纵轴进行分析,通过对采集的8000×2048数据点进行处理,每组数据为2048点,共分析8000次,每次代表一个时间点。2048点对应40MHz带宽内的信号频率,每一点频率为0.019532625MHz。
CWD时频分布为:
式中,μ为时间积分参数,τ为时延积分参数;β为衰减系数,β越小,对交叉项的抑制能力越大,但β的取值要在分辨率和交叉项抑制之间取折中,一般0.1<β<10;t为工作时间,f为频率,Sig为时域信号上角标#表示求共轭。
以上共经过8000次乔-威廉姆斯变换处理,每一次输出皆可得到CWD时频谱;再计算每组数据沿着时间轴的CWD时频谱最大值,将各组数据的最大值拼接起来,得到全部样本数据的CWD时频谱最大值图,如图4、图5所示,t为时间坐标轴,f为频率坐标轴,tfr为时频谱坐标轴。图4为无人机采集信号的乔-威廉姆斯时频谱图、图5为无人机采集信号的乔-威廉姆斯时频谱图在时频平面的等高线图。
S300:剔除随机噪声
在实际采集的无人机信号中会存在许多随机噪声,进行时频分析时,根据随机噪声的幅度通常弱于无人机信号的特点,通过设定能量门限,来消除随机噪声对遥控信号检测的影响。
假设阈值为CWD阈值(t,f),则
阈值设定采用基于时频谱能量统计的方法,即对时频谱进行能量累计,设置合适的阈值,可实现阈值的自适应效果;
式中,Nt=8000为时域样点数,Nf=2048为频域样点数。
S400:无人机遥控信号提取
根据在时频域上遥控信号与图传信号相关性存在的差异,利用低秩矩阵恢复的方法将相关性较强部分分配到低秩矩阵中,而稀疏矩阵中存放相关性弱的部分。图传信号在带宽范围内是连续存在的,而遥控信号在时频域上是一条频率随时间变化的阶梯状斜(或折)线,因此在分离两者时,大部分的图传信号分在低秩矩阵中,遥控信号被分在稀疏矩阵中,从而达到分离遥控信号与图传信号的目的。
低秩矩阵恢复先将时频谱矩阵CWD(t,f)表示为低秩矩阵A与稀疏(噪声)矩阵E之和,再通过求解范数优化问题来恢复低秩矩阵。
即求解下列最优化问题:
利用迭代阈值算法求解(1)式得到稀疏矩阵E。
进一步,求解方法如下:
a)将最优化问题式(1)正则化,便得到优化问题:
b)建立优化问题式(2)的拉格朗日函数:
式中,μ为平衡参数,这里取μ=0.05,Y为拉格朗日乘子矩阵,< >表示求内积。
c)迭代交替更新矩阵A,E和Y
当E=Ek,Y=Yk时,
当A=Ak+1,Y=Yk时,
当A=Ak+1,E=Ek+1时,
Yk+1=Yk+δk(CWD(t,f))-Ak+1-Ek+1)
其中δk为迭代步长,且0<δk<1,k为迭代次数。通过控制迭代次数,使其达到预期结果。或者说当Ek+1-Ek≈0时,结束迭代。
这里取k=5,则得到稀疏矩阵E,或称为采样数据时频序列。图6即为采样数据时频序列E的时频图。
S500:无人机遥控信号识别经过低秩矩阵恢复方法可以进一步抑制图传信号干扰,在稀疏矩阵E中得到清晰的遥控信号时频表示,从而提取出遥控信号的时频特征,这里如图6采样数据时频序列E的时频图中的折线,进一步与无人机遥控(跳频)信号(图样)标准库中的图样进行比对。
采用基于模式距离度量的时频序列相似性搜索方法,对提取的时频谱特征进行识别,如图2所示为基于模式距离度量的时频序列相似性搜索方法流程图,得到该无人机为大疆精灵4型。
具体实施方法如下:
S501:确定时频序列的模式集
以时频序列分段表示的每段直线与时间轴的夹角θ作为度量:
将时频序列的模式定义为五元集,记为模式集合M={-2,-1,0,1,2}。
于是,如图7无人机遥控信号标准库中大疆精灵4型的示意图,其时频序列分为10段,则大疆精灵4型无人机遥控信号标准时频序列为:
S1={(1,t1,1),(1,t1,2),(1,t1,3),(1,t1,4),(1,t1,5),(1,t1,6),(1,t1,7),(-2,t1,8)(-2,t1,9),(-2,t1,10)}
={s1,1,s1,2,…s1,i…s1,N}
式中,s1,i=(m1,i,t1,i),i(i=1,2,…N)为时频序列序号,N表示序列模式的划分数量,这里N=10;s1,i表示时频序列S1中的第i个模式。m为模式集合中的一个元素,这里m∈{-2,-1,0,1,2)
S502:沿时间轴滑动抽取时频序列从采集数据时频序列E中取与S1同样长度的时频序列S2,假设
S2={(1,t2,1),(1,t2,2),(1,t2,3),(1,t2,4),(1,t2,5),(1,t2,6),(1,t2,7),(=-2,t2,8)(-2,t2,9),(-2,t2,10)}
={s2,1,s2,2,…s2,i…s2,N}
式中,s2,i=(m2,i,t2,i),s2,i表示时频序列S2中的第i个模式,且假设t1,i=t2,i,i=1,2,…,N。
采集数据时频序列抽取方式如图8所示:沿时间轴从开始位置取与S1同样长度的时频序列,依次滑动一个时间段再取与S1同样长度的时频序列,直至取完所有的数据。
S503:估计时频序列的模式距离,并判别相似程度
时频序列的模式距离表示具有相同长度两个序列趋势的差异程度,则标准时频序列S1与采集数据时频序列S2的模式距离为:
其中,模式距离DM(s1,i,s2,i)=|m1,i-m2,i|,i=1,2,…,N,| |表示求绝对值;tN表示序列长度,tih表示第i个模式的保持时间,δ为相似度门限值,这里取δ=0.3。
此时采集数据时频序列S2与标准时频序列S1相似,可以判断该机型为大疆精灵4。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种基于时频分析的无人机遥控信号识别方法,其特征在于:针对无人机飞行状态时信息传输过程,其中无人机向地面遥控站传输图传信号,遥控站向无人机发射遥控信号,监测装置收集遥控信号和图传信号;该无人机遥控信号识别方法包含:
S100:数据采集及合并
采用基于零中频接收机架构的数据采集系统,对无人机信号进行采集;接收机内部的模拟/数据器件对接收的模拟信号进行模数转换,并通过正交解调、数字下变频将接收信号转换为同相分量I和正交分量Q两路数字基带信号,然后将I、Q两路信号进行合并,即时域信号Sin=I+j·Q,j表示虚数单位,
S200:对合并后的时域信号做乔-威廉姆斯变换得到CWD时频分布
采用乔-威廉姆斯变换对合并后的时域信号进行以时间为横轴,频率为纵轴进行分析,通过对采集的Num采数据点进行处理,每组数据为Num组点,共分析Num采/Num组次,每次代表一个时间点,Num组点对应Bw带宽内的信号频率,每一点频率为Bw/Num组;
CWD时频分布为:
式中,μ为时间积分参数,τ为时延积分参数;t为工作时间,f为频率,Sig为时域信号上角标#表示求共轭;
以上共经过Num采/Num组次乔-威廉姆斯变换处理,每一次输出得到CWD时频谱;再计算每组数据沿着时间轴的CWD时频谱最大值,将各组数据的最大值拼接起来,得到全部样本数据的CWD时频谱最大值图;
S300:剔除随机噪声
假设阈值为CWD阈值(t,f),则
阈值设定采用基于时频谱能量统计的方法,即对时频谱进行能量累计,设置合适的阈值,可实现阈值的自适应效果;
式中,Nt为时域样点数,Nf为频域样点数;
S400:无人机遥控信号提取
低秩矩阵恢复先将时频谱矩阵CWD(t,f)表示为低秩矩阵A与稀疏矩阵E之和,再通过求解范数优化问题来恢复低秩矩阵;
即求解下列最优化问题:
利用迭代阈值算法求解(1)式得到稀疏矩阵E;
S500:无人机遥控信号识别
经过低秩矩阵恢复方法可以进一步抑制图传信号干扰,在稀疏矩阵E中得到清晰的遥控信号时频表示,从而提取出遥控信号的时频特征,进一步与无人机遥控信号特征库中的图样进行比对;
采用基于模式距离度量的时频序列相似性搜索方法,对提取的时频谱特征进行识别,得到该无人机的型号。
2.根据权利要求1所述的一种基于时频分析的无人机遥控信号识别方法,其特征在于:所述步骤S400中,求解方法如下:
a)将最优化问题式(1)正则化,便得到优化问题:
b)建立优化问题式(2)的拉格朗日(Lagrange)函数:
式中,μ为平衡参数,这里取μ=0.05,Y为拉格朗日乘子矩阵,<>表示求内积;
c)迭代交替更新矩阵A,E和Y
当E=Ek,Y=Yk时,
当A=Ak+1,Y=Yk时,
当A=Ak+1,E=Ek+1时,
Yk+1=Yk+δk(CWD(t,f)-Ak+1-Ek+1)
其中δk为迭代步长,且0<δk<1,k为迭代次数;或者当Ek+1-Ek≈0时,结束迭代,得到稀疏矩阵E。
3.根据权利要求1所述的一种基于时频分析的无人机遥控信号识别方法,其特征在于:所述步骤S500中,无人机遥控信号识别,包含:
S501:以时频序列分段表示的每段直线与时间轴的夹角θ作为度量:
将时频序列的模式定义为五元集,记为模式集合M={-2,-1,0,1,2};
于是,某型无人机遥控信号标准时频序列为:
S1={s1,1,s1,2,…s1,i,…s1,N}
式中,s1,i=(m1,i,t1,i),i(i=1,2,…,N)为时频序列序号,N表示时频序列模式的划分数量,s1,i表示时频序列S1中的第i个模式。m1,i为模式集合中的一个元素,这里m1,i∈{-2,-1,0,1,2};
S502:沿时间轴滑动抽取时频序列
从采集数据时频序列E中取与S1同样长度的时频序列S2,假设
S2={s2,1,s2,2,…s2,i…s2,N}
式中,s2,i=(m2,i,t2,i),s2,i表示时频序列S2中的第i个模式,且假设t1,i=t2,i,i=1,2,…,N;
采集数据时频序列抽取时,沿时间轴从开始位置取与S1同样长度的时频序列,依次滑动一个时间段再取与S1同样长度的时频序列,直至取完所有的数据;
S503:估计时频序列的模式距离,并判别相似程度
时频序列的模式距离表示具有相同长度两个序列趋势的差异程度,则标准时频序列S1与采集数据时频序列S2的模式距离为:
其中,模式距离DM(s1,i,s2,i)=|m1,i-m2,i|,i=1,2,…,N,| |表示求绝对值;tN表示序列长度,tih表示第i个模式的保持时间,δ为相似度门限值,这里取δ=0.3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111372078.2A CN114095102B (zh) | 2021-11-18 | 2021-11-18 | 一种基于时频分析的无人机遥控信号识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111372078.2A CN114095102B (zh) | 2021-11-18 | 2021-11-18 | 一种基于时频分析的无人机遥控信号识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114095102A true CN114095102A (zh) | 2022-02-25 |
CN114095102B CN114095102B (zh) | 2023-12-29 |
Family
ID=80301971
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111372078.2A Active CN114095102B (zh) | 2021-11-18 | 2021-11-18 | 一种基于时频分析的无人机遥控信号识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114095102B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115085831A (zh) * | 2022-06-15 | 2022-09-20 | 西京学院 | 一种基于混合时频分析的无人机遥控信号识别系统和方法 |
CN115586798A (zh) * | 2022-12-12 | 2023-01-10 | 广东电网有限责任公司湛江供电局 | 一种无人机防坠毁方法和系统 |
CN117420527A (zh) * | 2023-12-19 | 2024-01-19 | 南京航天工业科技有限公司 | 一种基于宽带实时频谱分析的雷达监测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107301432A (zh) * | 2017-07-11 | 2017-10-27 | 哈尔滨工程大学 | 基于时频分析的自适应辐射源调制识别方法 |
CN110133632A (zh) * | 2019-05-20 | 2019-08-16 | 哈尔滨工程大学 | 一种基于cwd时频分析的复合调制信号识别方法 |
CN112560803A (zh) * | 2021-01-22 | 2021-03-26 | 南京航空航天大学 | 基于时频分析与机器学习的雷达信号调制识别方法 |
-
2021
- 2021-11-18 CN CN202111372078.2A patent/CN114095102B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107301432A (zh) * | 2017-07-11 | 2017-10-27 | 哈尔滨工程大学 | 基于时频分析的自适应辐射源调制识别方法 |
CN110133632A (zh) * | 2019-05-20 | 2019-08-16 | 哈尔滨工程大学 | 一种基于cwd时频分析的复合调制信号识别方法 |
CN112560803A (zh) * | 2021-01-22 | 2021-03-26 | 南京航空航天大学 | 基于时频分析与机器学习的雷达信号调制识别方法 |
Non-Patent Citations (2)
Title |
---|
甘泉;胡来招;: "一种基于短时信号改进CWD的识别方法", 四川大学学报(工程科学版), no. 03 * |
甘泉;胡来招;: "基于改进Choi-Williams分布的多分量信号频率识别方法", 信号处理, no. 06 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115085831A (zh) * | 2022-06-15 | 2022-09-20 | 西京学院 | 一种基于混合时频分析的无人机遥控信号识别系统和方法 |
CN115586798A (zh) * | 2022-12-12 | 2023-01-10 | 广东电网有限责任公司湛江供电局 | 一种无人机防坠毁方法和系统 |
CN115586798B (zh) * | 2022-12-12 | 2023-03-24 | 广东电网有限责任公司湛江供电局 | 一种无人机防坠毁方法和系统 |
CN117420527A (zh) * | 2023-12-19 | 2024-01-19 | 南京航天工业科技有限公司 | 一种基于宽带实时频谱分析的雷达监测方法 |
CN117420527B (zh) * | 2023-12-19 | 2024-03-12 | 南京航天工业科技有限公司 | 一种基于宽带实时频谱分析的雷达监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114095102B (zh) | 2023-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114095102A (zh) | 一种基于时频分析的无人机遥控信号识别方法 | |
CN108280395B (zh) | 一种对低小慢无人机飞控信号的高效识别方法 | |
CN104360332B (zh) | 基于地基合成孔径雷达干涉的大气相位屏提取方法 | |
CN107436427B (zh) | 空间目标运动轨迹与辐射信号关联方法 | |
CN106597404A (zh) | 太赫兹云雷达信号处理方法及系统 | |
CN109947238B (zh) | 一种基于wifi的非合作式手势识别的方法 | |
CN107194936B (zh) | 基于超像素联合稀疏表示的高光谱图像目标检测方法 | |
CN101794437B (zh) | 一种高光谱遥感影像异常探测方法 | |
CN107909018B (zh) | 一种稳健的多模态遥感影像匹配方法和系统 | |
CN111401168B (zh) | 一种无人机的多层雷达特征提取与选择方法 | |
CN109117747A (zh) | 基于循环平稳与深度卷积神经网络的雷达信号分类方法 | |
CN114154545A (zh) | 强互干扰条件下无人机测控信号智能识别方法 | |
CN113095145B (zh) | 像素对匹配与双窗口判别的高光谱异常检测深度学习方法 | |
CN116243273B (zh) | 一种针对植被冠层提取的光子计数激光雷达数据滤波方法 | |
CN114792107A (zh) | 一种基于时频分布的无人机遥控信号识别系统和方法 | |
CN112269182B (zh) | 目标雷达信号的确定方法和装置、存储介质、电子装置 | |
CN114118131A (zh) | 基于注意力机制的多特征融合无线设备射频指纹提取方法 | |
CN116359851A (zh) | 一种基于融合网络的雷达有源干扰检测识别方法及装置 | |
CN114581958B (zh) | 一种基于csi信号到达角估计的静态人体姿势估计方法 | |
CN110324098A (zh) | 一种宽带频谱序列的干扰信号的识别与统计方法 | |
CN107526064A (zh) | 基于二维特征的自适应lfm信号参数估计方法 | |
CN115085831A (zh) | 一种基于混合时频分析的无人机遥控信号识别系统和方法 | |
Schwartzkopf et al. | Radar Generalized Image Quality Equation Applied to Capella Open Dataset | |
CN109446910A (zh) | 一种通信辐射源信号分类识别方法 | |
CN112929141A (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 |