[发明专利]一种基于谐波能量分布的通信设备个体识别方法

说明书

技术领域

本发明属于现代通信系统安全和通信信号处理技术领域，更具体地，涉及一种基于谐波能量分布的通信设备个体识别方法。

背景技术

无线网络广泛渗透到各个领域，由此带来的网络系统安全问题也随之日益突出。在各种网络攻击手段中，入侵系统攻击的危害最大；这类攻击不仅窃取机密材料，而且随时都可能令受害系统陷入崩溃。有效抵御各种网络入侵的手段是稳定、及时、准确的实现网络访问对象的安全检测和控制。现有的无线网络安全机制主要依赖密钥系统对接入设备进行身份检测，而一旦发生密钥泄漏，将会形成极大的安全隐患。此外，大多数安全机制的研究仍然局限在无线网络的高层协议，复杂的加密系统和脆弱的软件协议系统难以解决入侵攻击的威胁。如果将安全措施引入到无线网络的物理层，通过无线信道监测和认证通信信号，使非法用户难以探测、伪造设备信号特征，将能够有效监测和控制现有无线网络的大多数非法设备。

国外相关研究主要集中在两个方面：从电台非稳定工作状态下的开机信号(暂态信号)提取特征实现电台个体识别的方法，以及通过测量电台稳态的工作参数来实现电台个体识别的方法。由于暂态信号持续时间短，实际应用中捕获暂态信号难度极大，因此难以应用到实际系统中。目前利用稳态信号进行通信设备识别的方法主要有参数测量法和双谱分析法。参数测量法是通过对调制信号的载频、码元速率等参数的测量来对通信设备进行识别，该方法能够区分不同的通信设备，但是对于不同的调制方式需要采用不同的测量方法，通用性不强，而且测量结果还会收到调制信号本身的干扰，难以得到准确的结果，从而影响区分能力。双谱分析法通过对信号的双谱进行计算，通过提取双谱特征来进行通信设备识别，虽然该方法具有通用性，但是双谱计算量大，耗时比较长，难以满足实时性的要求。

发明内容

针对现有技术的以上缺陷或改进需求，本发明提供了一种基于谐波能量分析的通信设备识别方法，通过计算调制信号中基波和二次谐波能量分布比例来对通信设备进行识别，由于谐波能量比例反映了通信设备功率放大器的参数特性，因此能够唯一的用于标识不同的通信设备。由于调制信息同时包含在基波和二次谐波中，它们的能量分布比例中消除调制信息本身的影响，同时本发明提出的方法可以采用快速傅立叶变换来实现，能够极大的提高运算速度，满足实际系统中实时性的要求。

为实现上述目的，本发明提供了一种基于谐波能量分布的通信设备个体识别方法，包括：

(1)对接收的通信设备的数字调制信号进行傅立叶变换得到信号频谱；

(2)根据信号的基波频谱成分和二次谐波频谱成分，计算谐波能量分布参数；

(3)将计算得到的谐波能量分布参数与预先存储的通信设备标准模板库进行比较，确定接收到的数字调制信号所属的通信设备，完成通信设备的个体识别。

具体地，对于连续形式，所述步骤(1)中对通信设备信号进行傅立叶变换得到信号频谱，根据下式计算：

其中，所述s(t)为接收的通信设备的数字调制信号，S(f)为经过傅立叶变换得到的信号频谱。

具体地，所述步骤(2)包括：

计算3个谐波幅度比例参数，其中A1表示基波幅度对称度，A2表示二次谐波幅度对称度，A3表示基波与二次谐波的幅度比例，其计算公式如下：

计算3个谐波能量比例参数，其中B1表示基波能量对称度，B2表示二次谐波能量对称度，B3表示基波与二次谐波的能量比例，其计算公式如下：

其中，S(f)为经过傅立叶变换得到的信号频谱，f_c为接收的通信设备的数字调制信号的载波频率，T为接收的通信设备的数字调制信号的码元周期。

具体地，所述步骤(3)包括：

将步骤(2)中计算得到的谐波能量分布参数与预先存储的通信设备标准模板库进行对比，计算接收信号参数与标准模板库中的模板偏差量，其计算公式如下：

其中，A1表示基波幅度对称度，A2表示二次谐波幅度对称度，A3表示基波与二次谐波的幅度比例，B1表示基波能量对称度，B2表示二次谐波能量对称度，B3表示基波与二次谐波的能量比例，m是模板编号，分布是第m个模板的幅度比例参数和能量比例参数；

选择D_m最小值对应的模板编号作为通信设备的识别结果。

具体地，对于离散形式，所述步骤(1)中对通信设备信号进行傅立叶变换得到信号频谱，具体根据下式计算：