买专利,只认龙图腾
首页 专利交易 科技果 科技人才 科技服务 商标交易 会员权益 IP管家助手 需求市场 关于龙图腾
 /  免费注册
到顶部 到底部
清空 搜索

【发明授权】一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法_北京邮电大学_201710396747.7 

申请/专利权人:北京邮电大学

申请日:2017-05-31

公开(公告)日:2020-09-04

公开(公告)号:CN107092898B

主分类号:G06K9/00(20060101)

分类号:G06K9/00(20060101)

优先权:

专利状态码:有效-授权

法律状态:2020.09.04#授权;2017.09.19#实质审查的生效;2017.08.25#公开

摘要:本发明公开了一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法,涉及无线通信领域,具体是发送端的某个比特流信号,采用QPSK映射得到信号sn;经过上变频后,得到调频信号pn,并输入到功率放大器中输出信号Φn,经过数模转换处理后得到模拟信号,将模拟信号发送出来并在发送过程中加入AWGN,接收端经过模数转换处理后得到数字信号rn,经过下变频后得到基带信号,从基带信号中提取射频指纹特征:双谱能量熵,一阶矩和二阶矩;然后,通过SVM分类器对射频指纹特征进行分类训练和测试,得到测试类别结果;通过将测试类别结果与其实际类别结果做对比,得到分类准确率Pc。本发明通过SVM分类器对射频信号进行有效分类,在低信噪下识别准确率相较于传统方法提升了近20%。

主权项:1.一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法,其特征在于,具体步骤如下:步骤一、针对OFDM系统发送端的某个比特流信号,采用QPSK映射得到QPSK信号sn;n=0,1,2…N-1,N为QPSK信号sn的长度;步骤二、将QPSK信号sn经过上变频后,得到调频信号pn,并输入到功率放大器中,输出信号Φn;调频信号pn计算如下:pn=snej2πnfT其中f为发送端的载波频率,为QPSK信号sn的采样间隔;功率放大器输出信号为: Ls为泰勒多项式的阶数,{al}为多项式系数;步骤三、将输出信号Φn经过数模转换处理后得到模拟信号,将模拟信号发送出来并在发送过程中加入高斯白噪声,接收端经过模数转换处理后得到数字信号rn;rn=Φn+υnυn为高斯白噪声;步骤四、将数字信号rn经过下变频后得到基带信号,从基带信号中提取射频指纹特征;射频的指纹特征包括:双谱能量熵,一阶矩和二阶矩;具体为:步骤401、将数字信号rn经过下变频后得到基带信号yn;yn=rne-j2πnfT步骤402、计算基带信号yn的三阶累积量c3yτ1,τ2,并利用三阶累积量求出信号双谱Bω1,ω2;三阶累积量计算如下:c3yτ1,τ2=E{y*nyn+τ1yn+τ2}*代表复共轭,τ1,τ2≥0代表时间的差值;对三阶累积量c3yτ1,τ2进行二维离散傅里叶变换得到基带信号双谱Bω1,ω2: ω1,ω2∈-π,π]代表数字角频率,分辨率为其中Nfft为傅里叶变换点数;根据傅里叶变换点数构造Nfft*Nfft的矩阵为基带信号双谱矩阵B;双谱矩阵B中第行列的数据为Bω1,ω2的值;步骤403、利用基带信号双谱Bω1,ω2求出射频指纹特征——双谱能量熵;首先,利用基带信号双谱Bω1,ω2计算双谱能量矩阵E中点i,j的双谱能量值Eij;计算如下: 其中i,j=1,2…Nfft;然后,计算双谱能量矩阵E中各点能量之和,并计算点i,j的能量在总能量中的占比pij;双谱能量矩阵E中各点能量之和为:占比pij表示为:最后,利用每个点的能量在总能量中的占比pij计算双谱能量熵En;如下: 步骤404、将双谱矩阵B转化为ζ-bit灰度图像矩阵G; 代表向下取整;Bm',n'代表双谱矩阵B中点m',n'的数据;Gm',n'代表灰度图像矩阵G中的点m',n'的数据,取值范围为0~2ζ-1;m',n'=1,2…Nfft;步骤405、分别求出灰度图像矩阵G的一阶矩μ和二阶矩ξ; 其中NB=Nfft×Nfft代表灰度图像矩阵G中点的总数;步骤五、通过SVM分类器对射频指纹特征进行分类训练和测试,得到测试类别结果;步骤六、将测试类别结果与其实际类别结果做对比,得到分类准确率Pc。

全文数据:一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法技术领域[0001]本发明涉及无线通信领域,具体是一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法。背景技术[0002]在无线通信系统中,由无线网络的开放性带来的安全问题不容忽视。[0003]传统的方法主要基于密码机制的安全协议来实现对数据完整性和机密性的保护,以及提供通信双方身份的认证。但是这样的认证信息很容易被恶意用户通过软件仿造,存在着潜在的威胁。考虑到即使是同一厂家生产的同一型号的不同设备,在制造过程中,由于氧化层厚度、掺杂浓度等不同会形成个体差异,而这些个体差异也会在通信信号中有所体现。[0004]无线通信的射频指纹提取和识别技术就是通过分析无线设备的通信信号来提取设备的“射频指纹”,从而进行设备识别。设备识别准确率很大程度上取决于指纹特征的选择,因此需要研究包含在通信信号中的设备指纹特征,以提高设备识别准确率。[0005]现有的射频指纹识别技术中存在低信噪比情况下,设备的识别准确率不高的问题;比如围线积分双谱中矩形积分双谱、圆周积分双谱和轴线积分双谱等,虽然具有高阶谱抑制高斯噪声、保留信号幅度和相位信息的特征,同时处理的方法比较简便而得到广泛应用。但是,围线积分双谱通过选择不同的围线积分路线将二维双谱矩阵转化为一维矩阵的过程中,忽略了很多双谱矩阵固有的特征,导致在低信噪比下识别准确率并不理想。发明内容[0006]本发明针对现有的射频指纹识别技术中,存在的低信噪比条件下识别准确率不高的问题,提出了一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法;[0007]具体步骤如下:[0008]步骤一、针对OFDM系统发送端的某个比特流信号,采用QPSK映射得到QPSK信号Sη;[0009]η=0,1,2·,·Ν-1,N为QPSK信号sη的长度;[0010]步骤二、将QPSK信号sη经过上变频后,得到调频信号ρη,并输入到功率放大器中,输出信号φη;[0011]调频信号ρη计算如下:[0012]ρn=snej2ltnfT[0013]其中f为发送端的载波频率,为QPSK信号sn的采样间隔。[0014]功率放大器采用泰勒多项式模型,功率放大器输出信号为:L·[0015]Φη=Ζ"[;"]’:1[0016]匕为泰勒多项式的阶数,{m}为多项式系数。[0017]步骤三、将输出信号φη经过数模转换处理后得到模拟信号,将模拟信号发送出来并在发送过程中加入高斯白噪声,接收端经过模数转换处理后得到数字信号rη;[0019]υη为高斯白噪声;[0020]步骤四、将数字信号rη经过下变频后得到基带信号,从基带信号中提取射频指纹特征;[0021]射频的指纹特征包括:双谱能量熵,一阶矩和二阶矩。[0022]具体包括:[0023]步骤401、将数字信号rη经过下变频后得到基带信号yη;[0025]步骤402、计算基带信号yη的三阶累积量C3yG1,τ2,并利用三阶累积量求出信号双谱Βωι,ω2;[0026]三阶累积量计算如下:[0028]*代表复共辄,T1,τ2彡〇代表时间的差值。[0029]对三阶累积量C3yT1,τ2进行二维离散傅里叶变换得到基带信号双谱Bω:,ω2:[0031]代表数字角频率,分辨率为,其中Nfft为傅里叶变换点数;根据傅里叶变换点数构造Nfft*Nfft的矩阵为基带信号双谱矩阵Β;双谱矩阵B中第列的数据为Βωι,ω2的值。[0032]步骤403、利用基带信号双谱Bωι,ω2求出射频指纹特征一一双谱能量熵;[0033]首先,利用基带信号双谱Bωι,ω2计算双谱能量矩阵E中点(i,j的双谱能量值Eij;[0034]计算如下:[0037]然后,计算双谱能量矩阵E中各点能量之和,并计算点(i,j的能量在总能量中的AKpij;[0038]双谱能量矩阵E中各点能量之和为:[0039]占比Pij表示为[0040]最后,利用每个点的能量在总能量中的占比叫计算双谱能量熵En;[0041]如下:[0043]步骤404、将双谱矩阵B转化为ζ-bit灰度图像矩阵G;[0045]L·」代表向下取整、,代表双谱矩阵B中点m’,n’)的数据;Gm,,n,代表灰度图像矩阵G中的点m’,n’)的数据,取值范围为0〜,n’=l,2…Nfft;[0046]步骤405、分别求出灰度图像矩阵G的一阶矩μ和二阶矩ξ;[0049]其中Nb=NfftXNfft代表灰度图像矩阵G中点的总数。[0050]步骤五、通过SVM分类器对射频指纹特征进行分类训练和测试,得到测试类别结果。[0051]具体为:[0052]步骤501、将发送端提取所有比特流的射频指纹特征组成特征向量[Εη,μ,ξ],划分成样本特征向量集合D和测试特征向量集合D'。[0053]步骤502、用样本特征向量集合D中的特征向量训练SVM分类器;[0054]步骤503、将测试特征向量集合D中的特征向量输入到训练好的SVM分类器中,得到测试类别结果。[0055]步骤六、将测试类别结果与其实际类别结果做对比,得到分类准确率Ρ。。[0056]本发明的优点在于:[0057]1、一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法,通过SVM分类器能够对射频信号进行有效分类。[0058]2、一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法,在低信噪情况下的识别准确率相较于传统的围线积分双谱特征提升了近20%。附图说明[0059]图1是本发明基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法信号发送和接收示意图;[0060]图2是本发明基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法的原理图;[0061]图3是本发明基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法的流程图;[0062]图4是本发明实施例中选取IOdB的SNR时,3种设备的射频指纹特征在特征空间的投影图;[0063]图5是本发明实施例中分类准确度与SNR的关系比较图。具体实施例[0064]下面结合附图对本发明的具体实施方法进行详细说明。[0065]一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别技术(ARadioFrequencyFingerprintingidentificationmethodbasedonbispectrumenergyentropyandcolormomentsofQPSKmodulationsignal,结合双谱能量熵,以及双谱矩阵转化为二维灰度数字图像后得到图像一阶矩和二阶矩,组成3维指纹特征进行设备识另IJ,应用于QPSK调制信号的射频指纹识别技术,有效的提高了低信噪比下射频设备的识别准确率,从而保证通信安全。[0066]如图1所示,信号从发送端到接收端的处理流程为:[0067]发送端的某个比特流信号,采用QPSK映射得到QPSK信号sη;经过上变频后,得到调频信号PH,并输入到功率放大器中,输出信号Φη;功率放大器是发送端最末端的一个元件,其非线性特征也是发送设备的主要标识特征;本发明从不同功率放大器的非线性特征的差异在接收端接收的信号中的不同体现入手,选择和提取射频指纹特征,进行设备识别。通过将输出信号Φη经过数模转换处理后得到模拟信号,将模拟信号发送出来并在发送过程中加入AWGN,接收端经过模数转换处理后得到数字信号rη,经过下变频后得到基带信号,从基带信号中提取射频指纹特征;[0068]然后,通过SVM分类器对射频指纹特征进行分类训练和测试,得到测试类别结果。具体过程如图2所示,通过将测试类别结果与其实际类别结果做对比,得到分类准确率Ρ。。[0069]如图3所,具体步骤如下:[0070]步骤一、针对OFDM系统发送端的某个比特流信号,采用QPSK映射得到QPSK信号sη;[0071]η=0,1,2…N-I,N为QPSK信号sη的长度;[0072]步骤二、将QPSK信号sη经过上变频后,得到调频信号ρη,并输入到功率放大器中,输出信号φη;[0073]调频信号ρη计算如下:[0075]其中f为发送端的载波频率,为QPSK信号sη的采样间隔。[0076]功率放大器采用泰勒多项式模型,功率放大器输出信号为:[0078]匕为泰勒多项式的阶数,{m}为多项式系数。[0079]步骤三、将输出信号Φη经过数模转换处理后得到模拟信号,将模拟信号发送出来并在发送过程中加入高斯白噪声,接收端经过模数转换处理后得到数字信号rη;[0081]υη为高斯白噪声;[0082]步骤四、将数字信号rη经过下变频后得到基带信号,从基带信号中提取射频指纹特征;[0083]射频的指纹特征包括:双谱能量熵,一阶矩和二阶矩。[0084]具体包括:[0085]步骤401、将数字信号rη经过下变频后得到基带信号yη;[0087]步骤402、计算基带信号yη的三阶累积量C3yG1,τ2,并利用三阶累积量求出信号双谱Βωι,ω2;[0088]三阶累积量计算如下:[0090]*代表复共辄,T1,τ2彡0代表时间的差值。[0091]对三阶累积量C3yT1,τ2进行二维离散傅里叶变换得到基带信号双谱Bω:,ω2:[0093]._代表数字角频率,分辨率为*其中Nfft为傅里叶变换点数;相掘傅里叶变换点数构造Nfft*Nfft的矩阵为基带信号双谱矩阵B;双谱矩阵B中第行-列的数据为Bω:,ω2的值。[0094]步骤403、利用基带信号双谱Bωι,ω2求出射频指纹特征一一双谱能量熵;[0095]首先,双谱Bωι,ω2的能量矩阵用E表示,利用基带信号双谱Bωι,ω2计算双谱能量矩阵E中点(i,j的双谱能量值Eij;[0096]双谱的;[0097]计算如下:[0100]然后,计算双谱能量矩阵E中各点能量之和,并计算点(i,j的能量在总能量中的占比Pij;[0101]矩阵E中各点能量之和E为:[0102]占比pij表不为:[0103]最后,利用每个点的能量在总能量中的占比Pij计算双谱能量熵En;[0104]为了衡量双谱的分布均匀程度,本发明采用能量熵来表示双谱能量的分布情况。如果能量在二维平面均匀分布,那么其能量熵值最大。反之,若能量主要集中在某些区域内,那么其能量熵值较小。[0105]如下:[0107]步骤404、将双谱矩阵B转化为ζ-bit灰度图像矩阵G;[0108]双谱可以看作一个包含了灰度强度信息的二维灰度数字图像,以ζ-bit灰度图像为例,通过下式将双谱矩阵转化为灰度图像矩阵:[0110]卜」代表向下取整、,代表双谱矩阵B中点m’,n’)的数据;Gm、,代表灰度图像矩阵G中的点m’,n’)的数据,取值范围为0〜,n’=l,2…Nfft;[0111]步骤405、分别求出灰度图像矩阵G的一阶矩μ和二阶矩ξ;[0112]本发明通过提取灰度图像的一阶和二阶颜色矩来描述两个图像灰度强度信息分布情况,一阶矩μ衡量灰度图像的平均强度,二阶矩ξ代表图像灰度强度的标准偏差。[0115]其中Nb=NfftXNfft代表灰度图像矩阵G中点的总数。[0116]步骤五、通过SVM分类器对射频指纹特征进行分类训练和测试,得到测试类别结果。[0117]具体为:[0118]步骤501、将发送端提取所有比特流的射频指纹特征组成特征向量[Εη,μ,ξ],划分成样本特征向量集合D和测试特征向量集合D'。[0119]步骤502、用样本特征向量集合D中的特征向量训练SVM分类器;[0120]步骤503、将测试特征向量集合D中的特征向量输入到训练好的SVM分类器中,得到测试类别结果。[0121]步骤六、将测试类别结果与其实际类别结果做对比,得到分类准确率Ρ。。[0122]本发明将双谱的能量熵和颜色矩结合作为设备的射频指纹特征[Εη,μ,ξ]。功率放大器泰勒多项式模型中泰勒多项式阶数Ls=3,并取三组{ai}作为3种设备功率放大器模型多项式系数。如图4所示,在信噪比SNR=IOdB时,3种设备的射频指纹特征在特征空间的投影图中可知,该射频指纹特征能够较好地对不同设备进行区分。[0123]在提取信号双谱能量熵和颜色矩的基础上,应用SVM分类器实现设备识别,在不同SNR下通过训练和检测设备,其识别情况如图5所示,由图可以发现在低信噪比下,设备正确识别率P。可以达到超过80%,并随着SNR的增加,P。不断上升,在20dB时达到超过95%的正确识别率。

权利要求:1.一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法,其特征在于,具体步骤如下:步骤一、针对OFDM系统发送端的某个比特流信号,采用QPSK映射得到QPSK信号sη;11=0,1,2〜114为0?31信号811的长度;步骤二、将QPSK信号sη经过上变频后,得到调频信号ρη,并输入到功率放大器中,输出信号Φη;调频信号Pη计算如下:ρn=snej2lInfT其中f为发送端的载波频率,为QPSK信号sη的采样间隔;功率放大器输出信号为:I-为泰勒多项式的阶数,{1}为多项式系数;步骤三、将输出信号Φη经过数模转换处理后得到模拟信号,将模拟信号发送出来并在发送过程中加入高斯白噪声,接收端经过模数转换处理后得到数字信号rη;rη=Φη+υηυη为高斯白噪声;步骤四、将数字信号rη经过下变频后得到基带信号,从基带信号中提取射频指纹特征;射频的指纹特征包括:双谱能量熵,一阶矩和二阶矩;步骤五、通过SVM分类器对射频指纹特征进行分类训练和测试,得到测试类别结果;步骤六、将测试类别结果与其实际类别结果做对比,得到分类准确率Ρ。。2.如权利要求1所述的一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法,其特征在于:所述的步骤二中功率放大器采用泰勒多项式模型。3.如权利要求1所述的一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法,其特征在于:所述的步骤四具体为:步骤401、将数字信号rη经过下变频后得到基带信号yη;yn=rne-j2imfT步骤402、计算基带信号yη的三阶累积量c3yτι,τ2,并利用三阶累积量求出信号双谱Bω1,ω2;三阶累积量计算如下:C3yτι,τ2=E{y*nyn+nyη+τ2}*代表复共辄,τι,τ2^0代表时间的差值;对三阶累积量c3yT1,τ2进行二维离散傅里叶变换得到基带信号双谱Bω:,ω2:ωι,ω2ε-π,π]代表数字角频率,分辨率为,其中Nfft为傅里叶变换点数;根据傅里叶变换点数构造Nfft*Nfft的矩阵为基带信号双谱矩阵B;双谱矩阵B中第的数据为Βωι,ω2的值;步骤403、利用基带信号双谱Bω:,ω2求出射频指纹特征一一双谱能量熵;首先,利用基带信号双谱Bωι,ω2计算双谱能量矩阵E中点(i,j的双谱能量值Eij;计算如下:其中i,j=l,2…Nfft;然后,计算双谱能量矩阵E中各点能量之和,并计算点(i,j的能量在总能量中的占比Pij;双谱能量矩阵E中各点能量之和为:占比Pij表不为:最后,利用每个点的能量在总能量中的占比PU计算双谱能量熵En;如下:步骤404、将双谱矩阵B转化为ζ-bit灰度图像矩阵G;L·」代表向下取整、,代表双谱矩阵B中点m’,n’)的数据;Gm,,n,代表灰度图像矩阵G中的点m’,n’)的数据,取值范围为0〜2ζ-1;πι’,n’=l,2-"Nfft;步骤405、分别求出灰度图像矩阵G的一阶矩μ和二阶矩ξ;其中Nb=NfftXNfft代表灰度图像矩阵G中点的总数。4.如权利要求1所述的一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法,其特征在于:所述的步骤五具体为:步骤501、将发送端提取所有比特流的射频指纹特征组成特征向量[Εη,μ,ξ],划分成样本特征向量集合D和测试特征向量集合D';步骤502、用样本特征向量集合D中的特征向量训练SVM分类器;步骤503、将测试特征向量集合D中的特征向量输入到训练好的SVM分类器中,得到测试类别结果。

百度查询: 北京邮电大学 一种基于QPSK信号双谱能量熵和颜色矩的射频指纹识别方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。